Challenges and potential solutions to health disparities in genomic medicine.

Significant disparities in the clinical usefulness of genomic information across diverse groups are due to underrepresentation in genetic databases and inequitable access to genetic services. Remedying disparities is immediately needed to ensure that genomic medicine is more equitable but will take a long-term commitment and active engagement of diverse communities.

De novo variants in FRYL are associated with developmental delay, intellectual disability, and dysmorphic features.

FRY-like transcription coactivator (FRYL) belongs to a Furry protein family that is evolutionarily conserved from yeast to humans. The functions of FRYL in mammals are largely unknown, and variants in FRYL have not previously been associated with a Mendelian disease. Here, we report fourteen individuals with heterozygous variants in FRYL who present with developmental delay, intellectual disability, dysmorphic features, and other congenital anomalies in multiple systems. The variants are confirmed de novo in all individuals except one. Human genetic data suggest that FRYL is intolerant to loss of function (LoF). We find that the fly FRYL ortholog, furry (fry), is expressed in multiple tissues, including the central nervous system where it is present in neurons but not in glia. Homozygous fry LoF mutation is lethal at various developmental stages, and loss of fry in mutant clones causes defects in wings and compound eyes. We next modeled four out of the five missense variants found in affected individuals using fry knockin alleles. One variant behaves as a severe LoF variant, whereas two others behave as partial LoF variants. One variant does not cause any observable defect in flies, and the corresponding human variant is not confirmed to be de novo, suggesting that this is a variant of uncertain significance. In summary, our findings support that fry is required for proper development in flies and that the LoF variants in FRYL cause a dominant disorder with developmental and neurological symptoms due to haploinsufficiency.

B56δ long-disordered arms form a dynamic PP2A regulation interface coupled with global allostery and Jordan's syndrome mutations.

Intrinsically disordered regions (IDR) and short linear motifs (SLiMs) play pivotal roles in the intricate signaling networks governed by phosphatases and kinases. B56δ (encoded by PPP2R5D) is a regulatory subunit of protein phosphatase 2A (PP2A) with long IDRs that harbor a substrate-mimicking SLiM and multiple phosphorylation sites. De novo missense mutations in PPP2R5D cause intellectual disabilities (ID), macrocephaly, Parkinsonism, and a broad range of neurological symptoms. Our single-particle cryo-EM structures of the PP2A-B56δ holoenzyme reveal that the long, disordered arms at the B56δ termini fold against each other and the holoenzyme core. This architecture suppresses both the phosphatase active site and the substrate-binding protein groove, thereby stabilizing the enzyme in a closed latent form with dual autoinhibition. The resulting interface spans over 190 Šand harbors unfavorable contacts, activation phosphorylation sites, and nearly all residues with ID-associated mutations. Our studies suggest that this dynamic interface is coupled to an allosteric network responsive to phosphorylation and altered globally by mutations. Furthermore, we found that ID mutations increase the holoenzyme activity and perturb the phosphorylation rates, and the severe variants significantly increase the mitotic duration and error rates compared to the normal variant.

Studying the impact of translational genomic research: Lessons from eMERGE.

Two major goals of the Electronic Medical Record and Genomics (eMERGE) Network are to learn how best to return research results to patient/participants and the clinicians who care for them and also to assess the impact of placing these results in clinical care. Yet since its inception, the Network has confronted a host of challenges in achieving these goals, many of which had ethical, legal, or social implications (ELSIs) that required consideration. Here, we share impediments we encountered in recruiting participants, returning results, and assessing their impact, all of which affected our ability to achieve the goals of eMERGE, as well as the steps we took to attempt to address these obstacles. We divide the domains in which we experienced challenges into four broad categories: (1) study design, including recruitment of more diverse groups; (2) consent; (3) returning results to participants and their health care providers (HCPs); and (4) assessment of follow-up care of participants and measuring the impact of research on participants and their families. Since most phases of eMERGE have included children as well as adults, we also address the particular ELSI posed by including pediatric populations in this research. We make specific suggestions for improving translational genomic research to ensure that future projects can effectively return results and assess their impact on patient/participants and providers if the goals of genomic-informed medicine are to be achieved.

Prospective, multi-site study of healthcare utilization after actionable monogenic findings from clinical sequencing.

As large-scale genomic screening becomes increasingly prevalent, understanding the influence of actionable results on healthcare utilization is key to estimating the potential long-term clinical impact. The eMERGE network sequenced individuals for actionable genes in multiple genetic conditions and returned results to individuals, providers, and the electronic health record. Differences in recommended health services (laboratory, imaging, and procedural testing) delivered within 12 months of return were compared among individuals with pathogenic or likely pathogenic (P/LP) findings to matched individuals with negative findings before and after return of results. Of 16,218 adults, 477 unselected individuals were found to have a monogenic risk for arrhythmia (n = 95), breast cancer (n = 96), cardiomyopathy (n = 95), colorectal cancer (n = 105), or familial hypercholesterolemia (n = 86). Individuals with P/LP results more frequently received services after return (43.8%) compared to before return (25.6%) of results and compared to individuals with negative findings (24.9%; p < 0.0001). The annual cost of qualifying healthcare services increased from an average of $162 before return to $343 after return of results among the P/LP group (p < 0.0001); differences in the negative group were non-significant. The mean difference-in-differences was $149 (p < 0.0001), which describes the increased cost within the P/LP group corrected for cost changes in the negative group. When stratified by individual conditions, significant cost differences were observed for arrhythmia, breast cancer, and cardiomyopathy. In conclusion, less than half of individuals received billed health services after monogenic return, which modestly increased healthcare costs for payors in the year following return.

PLS3 missense variants affecting the actin-binding domains cause X-linked congenital diaphragmatic hernia and body-wall defects.

Congenital diaphragmatic hernia (CDH) is a relatively common and genetically heterogeneous structural birth defect associated with high mortality and morbidity. We describe eight unrelated families with an X-linked condition characterized by diaphragm defects, variable anterior body-wall anomalies, and/or facial dysmorphism. Using linkage analysis and exome or genome sequencing, we found that missense variants in plastin 3 (PLS3), a gene encoding an actin bundling protein, co-segregate with disease in all families. Loss-of-function variants in PLS3 have been previously associated with X-linked osteoporosis (MIM: 300910), so we used in silico protein modeling and a mouse model to address these seemingly disparate clinical phenotypes. The missense variants in individuals with CDH are located within the actin-binding domains of the protein but are not predicted to affect protein structure, whereas the variants in individuals with osteoporosis are predicted to result in loss of function. A mouse knockin model of a variant identified in one of the CDH-affected families, c.1497G>C (p.Trp499Cys), shows partial perinatal lethality and recapitulates the key findings of the human phenotype, including diaphragm and abdominal-wall defects. Both the mouse model and one adult human male with a CDH-associated PLS3 variant were observed to have increased rather than decreased bone mineral density. Together, these clinical and functional data in humans and mice reveal that specific missense variants affecting the actin-binding domains of PLS3 might have a gain-of-function effect and cause a Mendelian congenital disorder.

A framework for an evidence-based gene list relevant to autism spectrum disorder.

Autism spectrum disorder (ASD) is often grouped with other brain-related phenotypes into a broader category of neurodevelopmental disorders (NDDs). In clinical practice, providers need to decide which genes to test in individuals with ASD phenotypes, which requires an understanding of the level of evidence for individual NDD genes that supports an association with ASD. Consensus is currently lacking about which NDD genes have sufficient evidence to support a relationship to ASD. Estimates of the number of genes relevant to ASD differ greatly among research groups and clinical sequencing panels, varying from a few to several hundred. This Roadmap discusses important considerations necessary to provide an evidence-based framework for the curation of NDD genes based on the level of information supporting a clinically relevant relationship between a given gene and ASD.

Exome sequencing and characterization of 49,960 individuals in the UK Biobank.

The UK Biobank is a prospective study of 502,543 individuals, combining extensive phenotypic and genotypic data with streamlined access for researchers around the world1. Here we describe the release of exome-sequence data for the first 49,960 study participants, revealing approximately 4 million coding variants (of which around 98.6% have a frequency of less than 1%). The data include 198,269 autosomal predicted loss-of-function (LOF) variants, a more than 14-fold increase compared to the imputed sequence. Nearly all genes (more than 97%) had at least one carrier with a LOF variant, and most genes (more than 69%) had at least ten carriers with a LOF variant. We illustrate the power of characterizing LOF variants in this population through association analyses across 1,730 phenotypes. In addition to replicating established associations, we found novel LOF variants with large effects on disease traits, including PIEZO1 on varicose veins, COL6A1 on corneal resistance, MEPE on bone density, and IQGAP2 and GMPR on blood cell traits. We further demonstrate the value of exome sequencing by surveying the prevalence of pathogenic variants of clinical importance, and show that 2% of this population has a medically actionable variant. Furthermore, we characterize the penetrance of cancer in carriers of pathogenic BRCA1 and BRCA2 variants. Exome sequences from the first 49,960 participants highlight the promise of genome sequencing in large population-based studies and are now accessible to the scientific community.

Dominantly acting KIF5B variants with pleiotropic cellular consequences cause variable clinical phenotypes.

Kinesins are motor proteins involved in microtubule (MT)-mediated intracellular transport. They contribute to key cellular processes, including intracellular trafficking, organelle dynamics and cell division. Pathogenic variants in kinesin-encoding genes underlie several human diseases characterized by an extremely variable clinical phenotype, ranging from isolated neurodevelopmental/neurodegenerative disorders to syndromic phenotypes belonging to a family of conditions collectively termed as 'ciliopathies.' Among kinesins, kinesin-1 is the most abundant MT motor for transport of cargoes towards the plus end of MTs. Three kinesin-1 heavy chain isoforms exist in mammals. Different from KIF5A and KIF5C, which are specifically expressed in neurons and established to cause neurological diseases when mutated, KIF5B is an ubiquitous protein. Three de novo missense KIF5B variants were recently described in four subjects with a syndromic skeletal disorder characterized by kyphomelic dysplasia, hypotonia and DD/ID. Here, we report three dominantly acting KIF5B variants (p.Asn255del, p.Leu498Pro and p.Leu537Pro) resulting in a clinically wide phenotypic spectrum, ranging from dilated cardiomyopathy with adult-onset ophthalmoplegia and progressive skeletal myopathy to a neurodevelopmental condition characterized by severe hypotonia with or without seizures. In vitro and in vivo analyses provide evidence that the identified disease-associated KIF5B variants disrupt lysosomal, autophagosome and mitochondrial organization, and impact cilium biogenesis. All variants, and one of the previously reported missense changes, were shown to affect multiple developmental processes in zebrafish. These findings document pleiotropic consequences of aberrant KIF5B function on development and cell homeostasis, and expand the phenotypic spectrum resulting from altered kinesin-mediated processes.

Is there a way to reduce the inequity in variant interpretation on the basis of ancestry?

The underrepresentation of non-European ancestry groups in current genomic databases complicates interpretation of their genetic test results, yielding a much higher prevalence of variants of uncertain significance (VUSs). Such VUS findings can frustrate the goals of genetic testing, create anxiety in patients, and lead to unnecessary medical interventions. Approaches to addressing underrepresentation of people with genetic ancestries other than European are being undertaken by broad-based recruitment efforts. However, some underrepresented groups have concerns that might preclude participation in such efforts. We describe here two initiatives aimed at meeting the needs of underrepresented ancestry groups in genomic datasets. The two communities, the Sephardi Jewish community in New York and First Peoples of Canada, have very different concerns about contributing to genomic research and datasets. Sephardi concerns focus on the possible negative effects of genetic findings on the marriage prospects of family members. Canadian Indigenous populations seek control over the research uses to which their genetic data would be put. Both cases involve targeted efforts to respond to the groups' concerns; these efforts include governance models aimed at ensuring that the data are used primarily to inform clinical test analyses and at achieving successful engagement and participation of community members. We suggest that these initiatives could provide models for other ancestral groups seeking to improve the accuracy and utility of clinical genetic testing while respecting the underlying preferences and values of community members with regard to the use of their genetic data.

OARD: Open annotations for rare diseases and their phenotypes based on real-world data.

Diagnosis for rare genetic diseases often relies on phenotype-driven methods, which hinge on the accuracy and completeness of the rare disease phenotypes in the underlying annotation knowledgebase. Existing knowledgebases are often manually curated with additional annotations found in published case reports. Despite their potential, real-world data such as electronic health records (EHRs) have not been fully exploited to derive rare disease annotations. Here, we present open annotation for rare diseases (OARD), a real-world-data-derived resource with annotation for rare-disease-related phenotypes. This resource is derived from the EHRs of two academic health institutions containing more than 10 million individuals spanning wide age ranges and different disease subgroups. By leveraging ontology mapping and advanced natural-language-processing (NLP) methods, OARD automatically and efficiently extracts concepts for both rare diseases and their phenotypic traits from billing codes and lab tests as well as over 100 million clinical narratives. The rare disease prevalence derived by OARD is highly correlated with those annotated in the original rare disease knowledgebase. By performing association analysis, we identified more than 1 million novel disease-phenotype association pairs that were previously missed by human annotation, and >60% were confirmed true associations via manual review of a list of sampled pairs. Compared to the manual curated annotation, OARD is 100% data driven and its pipeline can be shared across different institutions. By supporting privacy-preserving sharing of aggregated summary statistics, such as term frequencies and disease-phenotype associations, it fills an important gap to facilitate data-driven research in the rare disease community.

Bi-allelic variants in DOHH, catalyzing the last step of hypusine biosynthesis, are associated with a neurodevelopmental disorder.

Deoxyhypusine hydroxylase (DOHH) is the enzyme catalyzing the second step in the post-translational synthesis of hypusine [Nε-(4-amino-2-hydroxybutyl)lysine] in the eukaryotic initiation factor 5A (eIF5A). Hypusine is formed exclusively in eIF5A by two sequential enzymatic steps catalyzed by deoxyhypusine synthase (DHPS) and deoxyhypusine hydroxylase (DOHH). Hypusinated eIF5A is essential for translation and cell proliferation in eukaryotes, and all three genes encoding eIF5A, DHPS, and DOHH are highly conserved throughout eukaryotes. Pathogenic variants affecting either DHPS or EIF5A have been previously associated with neurodevelopmental disorders. Using trio exome sequencing, we identified rare bi-allelic pathogenic missense and truncating DOHH variants segregating with disease in five affected individuals from four unrelated families. The DOHH variants are associated with a neurodevelopmental phenotype that is similar to phenotypes caused by DHPS or EIF5A variants and includes global developmental delay, intellectual disability, facial dysmorphism, and microcephaly. A two-dimensional gel analyses revealed the accumulation of deoxyhypusine-containing eIF5A [eIF5A(Dhp)] and a reduction in the hypusinated eIF5A in fibroblasts derived from affected individuals, providing biochemical evidence for deficiency of DOHH activity in cells carrying the bi-allelic DOHH variants. Our data suggest that rare bi-allelic variants in DOHH result in reduced enzyme activity, limit the hypusination of eIF5A, and thereby lead to a neurodevelopmental disorder.

Loss-of-function variants in TIAM1 are associated with developmental delay, intellectual disability, and seizures.

TIAM Rac1-associated GEF 1 (TIAM1) regulates RAC1 signaling pathways that affect the control of neuronal morphogenesis and neurite outgrowth by modulating the actin cytoskeletal network. To date, TIAM1 has not been associated with a Mendelian disorder. Here, we describe five individuals with bi-allelic TIAM1 missense variants who have developmental delay, intellectual disability, speech delay, and seizures. Bioinformatic analyses demonstrate that these variants are rare and likely pathogenic. We found that the Drosophila ortholog of TIAM1, still life (sif), is expressed in larval and adult central nervous system (CNS) and is mainly expressed in a subset of neurons, but not in glia. Loss of sif reduces the survival rate, and the surviving adults exhibit climbing defects, are prone to severe seizures, and have a short lifespan. The TIAM1 reference (Ref) cDNA partially rescues the sif loss-of-function (LoF) phenotypes. We also assessed the function associated with three TIAM1 variants carried by two of the probands and compared them to the TIAM1 Ref cDNA function in vivo. TIAM1 p.Arg23Cys has reduced rescue ability when compared to TIAM1 Ref, suggesting that it is a partial LoF variant. In ectopic expression studies, both wild-type sif and TIAM1 Ref are toxic, whereas the three variants (p.Leu862Phe, p.Arg23Cys, and p.Gly328Val) show reduced toxicity, suggesting that they are partial LoF variants. In summary, we provide evidence that sif is important for appropriate neural function and that TIAM1 variants observed in the probands are disruptive, thus implicating loss of TIAM1 in neurological phenotypes in humans.

Neither cardiac mitochondrial DNA variation nor copy number contribute to congenital heart disease risk.

The well-established manifestation of mitochondrial mutations in functional cardiac disease (e.g., mitochondrial cardiomyopathy) prompted the hypothesis that mitochondrial DNA (mtDNA) sequence and/or copy number (mtDNAcn) variation contribute to cardiac defects in congenital heart disease (CHD). MtDNAcns were calculated and rare, non-synonymous mtDNA mutations were identified in 1,837 CHD-affected proband-parent trios, 116 CHD-affected singletons, and 114 paired cardiovascular tissue/blood samples. The variant allele fraction (VAF) of heteroplasmic variants in mitochondrial RNA from 257 CHD cardiovascular tissue samples was also calculated. On average, mtDNA from blood had 0.14 rare variants and 52.9 mtDNA copies per nuclear genome per proband. No variation with parental age at proband birth or CHD-affected proband age was seen. mtDNAcns in valve/vessel tissue (320 ± 70) were lower than in atrial tissue (1,080 ± 320, p = 6.8E-21), which were lower than in ventricle tissue (1,340 ± 280, p = 1.4E-4). The frequency of rare variants in CHD-affected individual DNA was indistinguishable from the frequency in an unaffected cohort, and proband mtDNAcns did not vary from those of CHD cohort parents. In both the CHD and the comparison cohorts, mtDNAcns were significantly correlated between mother-child, father-child, and mother-father. mtDNAcns among people with European (mean = 52.0), African (53.0), and Asian haplogroups (53.5) were calculated and were significantly different for European and Asian haplogroups (p = 2.6E-3). Variant heteroplasmic fraction (HF) in blood correlated well with paired cardiovascular tissue HF (r = 0.975) and RNA VAF (r = 0.953), which suggests blood HF is a reasonable proxy for HF in heart tissue. We conclude that mtDNA mutations and mtDNAcns are unlikely to contribute significantly to CHD risk.

The genetic architecture of pediatric cardiomyopathy.

To understand the genetic contribution to primary pediatric cardiomyopathy, we performed exome sequencing in a large cohort of 528 children with cardiomyopathy. Using clinical interpretation guidelines and targeting genes implicated in cardiomyopathy, we identified a genetic cause in 32% of affected individuals. Cardiomyopathy sub-phenotypes differed by ancestry, age at diagnosis, and family history. Infants < 1 year were less likely to have a molecular diagnosis (p < 0.001). Using a discovery set of 1,703 candidate genes and informatic tools, we identified rare and damaging variants in 56% of affected individuals. We see an excess burden of damaging variants in affected individuals as compared to two independent control sets, 1000 Genomes Project (p < 0.001) and SPARK parental controls (p < 1 × 10-16). Cardiomyopathy variant burden remained enriched when stratified by ancestry, variant type, and sub-phenotype, emphasizing the importance of understanding the contribution of these factors to genetic architecture. Enrichment in this discovery candidate gene set suggests multigenic mechanisms underlie sub-phenotype-specific causes and presentations of cardiomyopathy. These results identify important information about the genetic architecture of pediatric cardiomyopathy and support recommendations for clinical genetic testing in children while illustrating differences in genetic architecture by age, ancestry, and sub-phenotype and providing rationale for larger studies to investigate multigenic contributions.

Germline variants in tumor suppressor FBXW7 lead to impaired ubiquitination and a neurodevelopmental syndrome.

Neurodevelopmental disorders are highly heterogenous conditions resulting from abnormalities of brain architecture and/or function. FBXW7 (F-box and WD-repeat-domain-containing 7), a recognized developmental regulator and tumor suppressor, has been shown to regulate cell-cycle progression and cell growth and survival by targeting substrates including CYCLIN E1/2 and NOTCH for degradation via the ubiquitin proteasome system. We used a genotype-first approach and global data-sharing platforms to identify 35 individuals harboring de novo and inherited FBXW7 germline monoallelic chromosomal deletions and nonsense, frameshift, splice-site, and missense variants associated with a neurodevelopmental syndrome. The FBXW7 neurodevelopmental syndrome is distinguished by global developmental delay, borderline to severe intellectual disability, hypotonia, and gastrointestinal issues. Brain imaging detailed variable underlying structural abnormalities affecting the cerebellum, corpus collosum, and white matter. A crystal-structure model of FBXW7 predicted that missense variants were clustered at the substrate-binding surface of the WD40 domain and that these might reduce FBXW7 substrate binding affinity. Expression of recombinant FBXW7 missense variants in cultured cells demonstrated impaired CYCLIN E1 and CYCLIN E2 turnover. Pan-neuronal knockdown of the Drosophila ortholog, archipelago, impaired learning and neuronal function. Collectively, the data presented herein provide compelling evidence of an F-Box protein-related, phenotypically variable neurodevelopmental disorder associated with monoallelic variants in FBXW7.

SHINE: protein language model-based pathogenicity prediction for short inframe insertion and deletion variants.

Accurate variant pathogenicity predictions are important in genetic studies of human diseases. Inframe insertion and deletion variants (indels) alter protein sequence and length, but not as deleterious as frameshift indels. Inframe indel Interpretation is challenging due to limitations in the available number of known pathogenic variants for training. Existing prediction methods largely use manually encoded features including conservation, protein structure and function, and allele frequency to infer variant pathogenicity. Recent advances in deep learning modeling of protein sequences and structures provide an opportunity to improve the representation of salient features based on large numbers of protein sequences. We developed a new pathogenicity predictor for SHort Inframe iNsertion and dEletion (SHINE). SHINE uses pretrained protein language models to construct a latent representation of an indel and its protein context from protein sequences and multiple protein sequence alignments, and feeds the latent representation into supervised machine learning models for pathogenicity prediction. We curated training data from ClinVar and gnomAD, and created two test datasets from different sources. SHINE achieved better prediction performance than existing methods for both deletion and insertion variants in these two test datasets. Our work suggests that unsupervised protein language models can provide valuable information about proteins, and new methods based on these models can improve variant interpretation in genetic analyses.

A new type of blood-brain barrier aminoacidopathy underlies metabolic microcephaly associated with SLC1A4 mutations.

Mutations in the SLC1A4 transporter lead to neurodevelopmental impairments, spastic tetraplegia, thin corpus callosum, and microcephaly in children. SLC1A4 catalyzes obligatory amino acid exchange between neutral amino acids, but the physiopathology of SLC1A4 disease mutations and progressive microcephaly remain unclear. Here, we examined the phenotype and metabolic profile of three Slc1a4 mouse models, including a constitutive Slc1a4-KO mouse, a knock-in mouse with the major human Slc1a4 mutation (Slc1a4-K256E), and a selective knockout of Slc1a4 in brain endothelial cells (Slc1a4tie2-cre). We show that Slc1a4 is a bona fide L-serine transporter at the BBB and that acute inhibition or deletion of Slc1a4 leads to a decrease in serine influx into the brain. This results in microcephaly associated with decreased L-serine content in the brain, accumulation of atypical and cytotoxic 1-deoxysphingolipids in the brain, neurodegeneration, synaptic and mitochondrial abnormalities, and behavioral impairments. Prenatal and early postnatal oral administration of L-serine at levels that replenish the serine pool in the brain rescued the observed biochemical and behavioral changes. Administration of L-serine till the second postnatal week also normalized brain weight in Slc1a4-E256 K mice. Our observations suggest that the transport of "non-essential" amino acids from the blood through the BBB is at least as important as that of essential amino acids for brain metabolism and development. We proposed that SLC1A4 mutations cause a BBB aminoacidopathy with deficits in serine import across the BBB required for optimal brain growth and leads to a metabolic microcephaly, which may be amenable to treatment with L-serine.

De novo and bi-allelic variants in AP1G1 cause neurodevelopmental disorder with developmental delay, intellectual disability, and epilepsy.

Adaptor protein (AP) complexes mediate selective intracellular vesicular trafficking and polarized localization of somatodendritic proteins in neurons. Disease-causing alleles of various subunits of AP complexes have been implicated in several heritable human disorders, including intellectual disabilities (IDs). Here, we report two bi-allelic (c.737C>A [p.Pro246His] and c.1105A>G [p.Met369Val]) and eight de novo heterozygous variants (c.44G>A [p.Arg15Gln], c.103C>T [p.Arg35Trp], c.104G>A [p.Arg35Gln], c.229delC [p.Gln77Lys∗11], c.399_400del [p.Glu133Aspfs∗37], c.747G>T [p.Gln249His], c.928-2A>C [p.?], and c.2459C>G [p.Pro820Arg]) in AP1G1, encoding gamma-1 subunit of adaptor-related protein complex 1 (AP1γ1), associated with a neurodevelopmental disorder (NDD) characterized by mild to severe ID, epilepsy, and developmental delay in eleven families from different ethnicities. The AP1γ1-mediated adaptor complex is essential for the formation of clathrin-coated intracellular vesicles. In silico analysis and 3D protein modeling simulation predicted alteration of AP1γ1 protein folding for missense variants, which was consistent with the observed altered AP1γ1 levels in heterologous cells. Functional studies of the recessively inherited missense variants revealed no apparent impact on the interaction of AP1γ1 with other subunits of the AP-1 complex but rather showed to affect the endosome recycling pathway. Knocking out ap1g1 in zebrafish leads to severe morphological defect and lethality, which was significantly rescued by injection of wild-type AP1G1 mRNA and not by transcripts encoding the missense variants. Furthermore, microinjection of mRNAs with de novo missense variants in wild-type zebrafish resulted in severe developmental abnormalities and increased lethality. We conclude that de novo and bi-allelic variants in AP1G1 are associated with neurodevelopmental disorder in diverse populations.

Rare and de novo variants in 827 congenital diaphragmatic hernia probands implicate LONP1 as candidate risk gene.

Congenital diaphragmatic hernia (CDH) is a severe congenital anomaly that is often accompanied by other anomalies. Although the role of genetics in the pathogenesis of CDH has been established, only a small number of disease-associated genes have been identified. To further investigate the genetics of CDH, we analyzed de novo coding variants in 827 proband-parent trios and confirmed an overall significant enrichment of damaging de novo variants, especially in constrained genes. We identified LONP1 (lon peptidase 1, mitochondrial) and ALYREF (Aly/REF export factor) as candidate CDH-associated genes on the basis of de novo variants at a false discovery rate below 0.05. We also performed ultra-rare variant association analyses in 748 affected individuals and 11,220 ancestry-matched population control individuals and identified LONP1 as a risk gene contributing to CDH through both de novo and ultra-rare inherited largely heterozygous variants clustered in the core of the domains and segregating with CDH in affected familial individuals. Approximately 3% of our CDH cohort who are heterozygous with ultra-rare predicted damaging variants in LONP1 have a range of clinical phenotypes, including other anomalies in some individuals and higher mortality and requirement for extracorporeal membrane oxygenation. Mice with lung epithelium-specific deletion of Lonp1 die immediately after birth, most likely because of the observed severe reduction of lung growth, a known contributor to the high mortality in humans. Our findings of both de novo and inherited rare variants in the same gene may have implications in the design and analysis for other genetic studies of congenital anomalies.