MARRVEL: Integration of Human and Model Organism Genetic Resources to Facilitate Functional Annotation of the Human Genome

Functional analysis of conserved C. elegans bHLH family members uncovers lifespan control by a peptidergic hub neuron

Throughout the animal kingdom, several members of the basic helix-loop-helix (bHLH) family act as proneural genes during early steps of nervous system development. Roles of bHLH genes in specifying terminal differentiation of postmitotic neurons have been less extensively studied. We analyze here the function of five C. elegans bHLH genes, falling into three phylogenetically conserved subfamilies, which are continuously expressed in a very small number of postmitotic neurons in the central nervous system. We show (a) that two orthologs of the vertebrate bHLHb4/b5 genes, called hlh-17 and hlh-32, function redundantly to specify the identity of a single head interneuron (AUA), as well as an individual motor neuron (VB2), (b) that the PTF1a ortholog hlh-13 acts as a terminal selector to control terminal differentiation and function of the sole octopaminergic neuron class in C. elegans , RIC, and (c) that the NHLH1/2 ortholog hlh-15 controls terminal differentiation and function of the peptidergic AVK head interneuron class, a known neuropeptidergic signaling hub in the animal. Strikingly, through null mutant analysis and cell-specific rescue experiments, we find that loss of hlh-15/NHLH in the peptidergic AVK neurons and the resulting abrogation of neuropeptide secretion causes a substantially expanded lifespan of the animal, revealing an unanticipated impact of a central, peptidergic hub neuron in regulating lifespan, which we propose to be akin to hypothalamic control of lifespan in vertebrates. Taken together, our functional analysis reveals themes of bHLH gene function during terminal differentiation that are complementary to the earlier lineage specification roles of other bHLH family members. However, such late functions are much more sparsely employed by members of the bHLH transcription factor family, compared to the function of the much more broadly employed homeodomain transcription factor family.

Homozygous missense variants in YKT6 result in loss of function and are associated with developmental delay, with or without severe infantile liver disease and risk for hepatocellular carcinoma

PURPOSE

YKT6 plays important roles in multiple intracellular vesicle trafficking events but has not been associated with Mendelian diseases.

METHODS

We report 3 unrelated individuals with rare homozygous missense variants in YKT6 who exhibited neurological disease with or without a progressive infantile liver disease. We modeled the variants in Drosophila. We generated wild-type and variant genomic rescue constructs of the fly ortholog dYkt6 and compared their ability in rescuing the loss-of-function phenotypes in mutant flies. We also generated a dYkt6KozakGAL4 allele to assess the expression pattern of dYkt6.

RESULTS

Two individuals are homozygous for YKT6 [NM_006555.3:c.554A>G p.(Tyr185Cys)] and exhibited normal prenatal course followed by failure to thrive, developmental delay, and progressive liver disease. Haplotype analysis identified a shared homozygous region flanking the variant, suggesting a common ancestry. The third individual is homozygous for YKT6 [NM_006555.3:c.191A>G p.(Tyr64Cys)] and exhibited neurodevelopmental disorders and optic atrophy. Fly dYkt6 is essential and is expressed in the fat body (analogous to liver) and central nervous system. Wild-type genomic rescue constructs can rescue the lethality and autophagic flux defects, whereas the variants are less efficient in rescuing the phenotypes.

CONCLUSION

The YKT6 variants are partial loss-of-function alleles, and the p.(Tyr185Cys) is more severe than p.(Tyr64Cys).

Dominant missense variants in SREBF2 are associated with complex dermatological, neurological, and skeletal abnormalities

PURPOSE

We identified 2 individuals with de novo variants in SREBF2 that disrupt a conserved site 1 protease (S1P) cleavage motif required for processing SREBP2 into its mature transcription factor. These individuals exhibit complex phenotypic manifestations that partially overlap with sterol regulatory element binding proteins (SREBP) pathway-related disease phenotypes, but SREBF2-related disease has not been previously reported. Thus, we set out to assess the effects of SREBF2 variants on SREBP pathway activation.

METHODS

We undertook ultrastructure and gene expression analyses using fibroblasts from an affected individual and utilized a fly model of lipid droplet (LD) formation to investigate the consequences of SREBF2 variants on SREBP pathway function.

RESULTS

We observed reduced LD formation, endoplasmic reticulum expansion, accumulation of aberrant lysosomes, and deficits in SREBP2 target gene expression in fibroblasts from an affected individual, indicating that the SREBF2 variant inhibits SREBP pathway activation. Using our fly model, we discovered that SREBF2 variants fail to induce LD production and act in a dominant-negative manner, which can be rescued by overexpression of S1P.

CONCLUSION

Taken together, these data reveal a mechanism by which SREBF2 pathogenic variants that disrupt the S1P cleavage motif cause disease via dominant-negative antagonism of S1P, limiting the cleavage of S1P targets, including SREBP1 and SREBP2.

Exome sequencing implicates ancestry-related Mendelian variation at SYNE1 in childhood-onset essential hypertension

Childhood-onset essential hypertension (COEH) is an uncommon form of hypertension that manifests in childhood or adolescence and, in the United States, disproportionately affects children of African ancestry. The etiology of COEH is unknown, but its childhood onset, low prevalence, high heritability, and skewed ancestral demography suggest the potential to identify rare genetic variation segregating in a Mendelian manner among affected individuals and thereby implicate genes important to disease pathogenesis. However, no COEH genes have been reported to date. Here, we identify recessive segregation of rare and putatively damaging missense variation in the spectrin domain of spectrin repeat containing nuclear envelope protein 1 (SYNE1), a cardiovascular candidate gene, in 3 of 16 families with early-onset COEH without an antecedent family history. By leveraging exome sequence data from an additional 48 COEH families, 1,700 in-house trios, and publicly available data sets, we demonstrate that compound heterozygous SYNE1 variation in these COEH individuals occurred more often than expected by chance and that this class of biallelic rare variation was significantly enriched among individuals of African genetic ancestry. Using in vitro shRNA knockdown of SYNE1, we show that reduced SYNE1 expression resulted in a substantial decrease in the elasticity of smooth muscle vascular cells that could be rescued by pharmacological inhibition of the downstream RhoA/Rho-associated protein kinase pathway. These results provide insights into the molecular genetics and underlying pathophysiology of COEH and suggest a role for precision therapeutics in the future.

AI-MARRVEL - A Knowledge-Driven AI System for Diagnosing Mendelian Disorders

BACKGROUND

Diagnosing genetic disorders requires extensive manual curation and interpretation of candidate variants, a labor-intensive task even for trained geneticists. Although artificial intelligence (AI) shows promise in aiding these diagnoses, existing AI tools have only achieved moderate success for primary diagnosis.

METHODS

AI-MARRVEL (AIM) uses a random-forest machine-learning classifier trained on over 3.5 million variants from thousands of diagnosed cases. AIM additionally incorporates expert-engineered features into training to recapitulate the intricate decision-making processes in molecular diagnosis. The online version of AIM is available at https://ai.marrvel.org. To evaluate AIM, we benchmarked it with diagnosed patients from three independent cohorts.

RESULTS

AIM improved the rate of accurate genetic diagnosis, doubling the number of solved cases as compared with benchmarked methods, across three distinct real-world cohorts. To better identify diagnosable cases from the unsolved pools accumulated over time, we designed a confidence metric on which AIM achieved a precision rate of 98% and identified 57% of diagnosable cases out of a collection of 871 cases. Furthermore, AIM's performance improved after being fine-tuned for targeted settings including recessive disorders and trio analysis. Finally, AIM demonstrated potential for novel disease gene discovery by correctly predicting two newly reported disease genes from the Undiagnosed Diseases Network.

CONCLUSIONS

AIM achieved superior accuracy compared with existing methods for genetic diagnosis. We anticipate that this tool may aid in primary diagnosis, reanalysis of unsolved cases, and the discovery of novel disease genes. (Funded by the NIH Common Fund and others.).

Cdk8/CDK19 promotes mitochondrial fission through Drp1 phosphorylation and can phenotypically suppress pink1 deficiency in Drosophila

Cdk8 in Drosophila is the orthologue of vertebrate CDK8 and CDK19. These proteins have been shown to modulate transcriptional control by RNA polymerase II. We found that neuronal loss of Cdk8 severely reduces fly lifespan and causes bang sensitivity. Remarkably, these defects can be rescued by expression of human CDK19, found in the cytoplasm of neurons, suggesting a non-nuclear function of CDK19/Cdk8. Here we show that Cdk8 plays a critical role in the cytoplasm, with its loss causing elongated mitochondria in both muscles and neurons. We find that endogenous GFP-tagged Cdk8 can be found in both the cytoplasm and nucleus. We show that Cdk8 promotes the phosphorylation of Drp1 at S616, a protein required for mitochondrial fission. Interestingly, Pink1, a mitochondrial kinase implicated in Parkinson's disease, also phosphorylates Drp1 at the same residue. Indeed, overexpression of Cdk8 significantly suppresses the phenotypes observed in flies with low levels of Pink1, including elevated levels of ROS, mitochondrial dysmorphology, and behavioral defects. In summary, we propose that Pink1 and Cdk8 perform similar functions to promote Drp1-mediated fission.

The relationship between intraflagellar transport and upstream protein trafficking pathways and macrocyclic lactone resistance in Caenorhabditis elegans

Parasitic nematodes are globally important and place a heavy disease burden on infected humans, crops, and livestock, while commonly administered anthelmintics used for treatment are being rendered ineffective by increasing levels of resistance. It has recently been shown in the model nematode Caenorhabditis elegans that the sensory cilia of the amphid neurons play an important role in resistance toward macrocyclic lactones such as ivermectin (an avermectin) and moxidectin (a milbemycin) either through reduced uptake or intertissue signaling pathways. This study interrogated the extent to which ciliary defects relate to macrocyclic lactone resistance and dye-filling defects using a combination of forward genetics and targeted resistance screening approaches and confirmed the importance of intraflagellar transport in this process. This approach also identified the protein trafficking pathways used by the downstream effectors and the components of the ciliary basal body that are required for effector entry into these nonmotile structures. In total, 24 novel C. elegans anthelmintic survival-associated genes were identified in this study. When combined with previously known resistance genes, there are now 46 resistance-associated genes that are directly involved in amphid, cilia, and intraflagellar transport function.

FAM81A is a postsynaptic protein that regulates the condensation of postsynaptic proteins via liquid-liquid phase separation

Proteome analyses of the postsynaptic density (PSD), a proteinaceous specialization beneath the postsynaptic membrane of excitatory synapses, have identified several thousands of proteins. While proteins with predictable functions have been well studied, functionally uncharacterized proteins are mostly overlooked. In this study, we conducted a comprehensive meta-analysis of 35 PSD proteome datasets, encompassing a total of 5,869 proteins. Employing a ranking methodology, we identified 97 proteins that remain inadequately characterized. From this selection, we focused our detailed analysis on the highest-ranked protein, FAM81A. FAM81A interacts with PSD proteins, including PSD-95, SynGAP, and NMDA receptors, and promotes liquid-liquid phase separation of those proteins in cultured cells or in vitro. Down-regulation of FAM81A in cultured neurons causes a decrease in the size of PSD-95 puncta and the frequency of neuronal firing. Our findings suggest that FAM81A plays a crucial role in facilitating the interaction and assembly of proteins within the PSD, and its presence is important for maintaining normal synaptic function. Additionally, our methodology underscores the necessity for further characterization of numerous synaptic proteins that still lack comprehensive understanding.

Loss of the endoplasmic reticulum protein Tmem208 affects cell polarity, development, and viability

Nascent proteins destined for the cell membrane and the secretory pathway are targeted to the endoplasmic reticulum (ER) either posttranslationally or cotranslationally. The signal-independent pathway, containing the protein TMEM208, is one of three pathways that facilitates the translocation of nascent proteins into the ER. The in vivo function of this protein is ill characterized in multicellular organisms. Here, we generated a CRISPR-induced null allele of the fruit fly ortholog CG8320/Tmem208 by replacing the gene with the Kozak-GAL4 sequence. We show that Tmem208 is broadly expressed in flies and that its loss causes lethality, although a few short-lived flies eclose. These animals exhibit wing and eye developmental defects consistent with impaired cell polarity and display mild ER stress. Tmem208 physically interacts with Frizzled (Fz), a planar cell polarity (PCP) receptor, and is required to maintain proper levels of Fz. Moreover, we identified a child with compound heterozygous variants in TMEM208 who presents with developmental delay, skeletal abnormalities, multiple hair whorls, cardiac, and neurological issues, symptoms that are associated with PCP defects in mice and humans. Additionally, fibroblasts of the proband display mild ER stress. Expression of the reference human TMEM208 in flies fully rescues the loss of Tmem208, and the two proband-specific variants fail to rescue, suggesting that they are loss-of-function alleles. In summary, our study uncovers a role of TMEM208 in development, shedding light on its significance in ER homeostasis and cell polarity.

De novo variants in PLCG1 are associated with hearing impairment, ocular pathology, and cardiac defects

Phospholipase C isozymes (PLCs) hydrolyze phosphatidylinositol 4,5-bisphosphate into inositol 1,4,5-trisphosphate and diacylglycerol, important signaling molecules involved in many cellular processes. PLCG1 encodes the PLCγ1 isozyme that is broadly expressed. Hyperactive somatic mutations of PLCG1 are observed in multiple cancers, but only one germline variant has been reported. Here we describe three unrelated individuals with de novo heterozygous missense variants in PLCG1 (p.Asp1019Gly, p.His380Arg, and p.Asp1165Gly) who exhibit variable phenotypes including hearing loss, ocular pathology and cardiac septal defects. To model these variants in vivo, we generated the analogous variants in the Drosophila ortholog, small wing (sl). We created a null allele slT2A and assessed the expression pattern. sl is broadly expressed, including in wing discs, eye discs, and a subset of neurons and glia. Loss of sl causes wing size reductions, ectopic wing veins and supernumerary photoreceptors. We document that mutant flies exhibit a reduced lifespan and age-dependent locomotor defects. Expressing wild-type sl in slT2A mutant rescues the loss-of-function phenotypes whereas expressing the variants causes lethality. Ubiquitous overexpression of the variants also reduces viability, suggesting that the variants are toxic. Ectopic expression of an established hyperactive PLCG1 variant (p.Asp1165His) in the wing pouch causes severe wing phenotypes, resembling those observed with overexpression of the p.Asp1019Gly or p.Asp1165Gly variants, further arguing that these two are gain-of-function variants. However, the wing phenotypes associated with p.His380Arg overexpression are mild. Our data suggest that the PLCG1 de novo heterozygous missense variants are pathogenic and contribute to the features observed in the probands.

Integrating non-mammalian model organisms in the diagnosis of rare genetic diseases in humans

Next-generation sequencing technology has rapidly accelerated the discovery of genetic variants of interest in individuals with rare diseases. However, showing that these variants are causative of the disease in question is complex and may require functional studies. Use of non-mammalian model organisms - mainly fruitflies (Drosophila melanogaster), nematode worms (Caenorhabditis elegans) and zebrafish (Danio rerio) - enables the rapid and cost-effective assessment of the effects of gene variants, which can then be validated in mammalian model organisms such as mice and in human cells. By probing mechanisms of gene action and identifying interacting genes and proteins in vivo, recent studies in these non-mammalian model organisms have facilitated the diagnosis of numerous genetic diseases and have enabled the screening and identification of therapeutic options for patients. Studies in non-mammalian model organisms have also shown that the biological processes underlying rare diseases can provide insight into more common mechanisms of disease and the biological functions of genes. Here, we discuss the opportunities afforded by non-mammalian model organisms, focusing on flies, worms and fish, and provide examples of their use in the diagnosis of rare genetic diseases.

Matching variants for functional characterization of genetic variants

Rapid and low-cost sequencing, as well as computer analysis, have facilitated the diagnosis of many genetic diseases, resulting in a substantial rise in the number of disease-associated genes. However, genetic diagnosis of many disorders remains problematic due to the lack of interpretation for many genetic variants, especially missenses, the infeasibility of high-throughput experiments on mammals, and the shortcomings of computational prediction technologies. Additionally, the available mutant databases are not well-utilized. Toward this end, we used Caenorhabditis elegans mutant resources to delineate the functions of eight missense variants (V444I, V517D, E610K, L732F, E817K, H873P, R1105K, and G1205E) and two stop codons (W937stop and Q1434stop), including several matching variants (MatchVar) with human in ciliopathy associated IFT-140 (also called CHE-11)//IFT140 (intraflagellar transport protein 140). Moreover, MatchVars carrying C. elegans mutants, including IFT-140(G680S) and IFT-140(P702A) for the human (G704S) (dbSNP: rs150745099) and P726A (dbSNP: rs1057518064 and a conflicting variation) were created using CRISPR/Cas9. IFT140 is a key component of IFT complex A (IFT-A), which is involved in the retrograde transport of IFT along cilia and the entrance of G protein-coupled receptors into cilia. Functional analysis of all 10 variants revealed that P702A and W937stop, but not others phenocopied the ciliary phenotypes (short cilia, IFT accumulations, mislocalization of membrane proteins, and cilia entry of nonciliary proteins) of the IFT-140 null mutant, indicating that both P702A and W937stop are phenotypic in C. elegans. Our functional data offered experimental support for interpreting human variants, by using ready-to-use mutants carrying MatchVars and generating MatchVars with CRISPR/Cas9.

Finding information about uncharacterized Drosophila melanogaster genes

Genes that have been identified in the genome but remain uncharacterized with regards to function offer an opportunity to uncover novel biological information. Novelty is exciting but can also be a barrier. If nothing is known, how does one start planning and executing experiments? Here, we provide a recommended information-mining workflow and a corresponding guide to accessing information about uncharacterized Drosophila melanogaster genes, such as those assigned only a systematic coding gene identifier. The available information can provide insights into where and when the gene is expressed, what the function of the gene might be, whether there are similar genes in other species, whether there are known relationships to other genes, and whether any other features have already been determined. In addition, available information about relevant reagents can inspire and facilitate experimental studies. Altogether, mining available information can help prioritize genes for further study, as well as provide starting points for experimental assays and other analyses.

Molecular and functional characterization of the Drosophila melanogaster conserved smORFome

Short polypeptides encoded by small open reading frames (smORFs) are ubiquitously found in eukaryotic genomes and are important regulators of physiology, development, and mitochondrial processes. Here, we focus on a subset of 298 smORFs that are evolutionarily conserved between Drosophila melanogaster and humans. Many of these smORFs are conserved broadly in the bilaterian lineage, and ∼182 are conserved in plants. We observe remarkably heterogeneous spatial and temporal expression patterns of smORF transcripts-indicating wide-spread tissue-specific and stage-specific mitochondrial architectures. In addition, an analysis of annotated functional domains reveals a predicted enrichment of smORF polypeptides localizing to mitochondria. We conduct an embryonic ribosome profiling experiment and find support for translation of 137 of these smORFs during embryogenesis. We further embark on functional characterization using CRISPR knockout/activation, RNAi knockdown, and cDNA overexpression, revealing diverse phenotypes. This study underscores the importance of identifying smORF function in disease and phenotypic diversity.

Rare de novo gain-of-function missense variants in DOT1L are associated with developmental delay and congenital anomalies

Misregulation of histone lysine methylation is associated with several human cancers and with human developmental disorders. DOT1L is an evolutionarily conserved gene encoding a lysine methyltransferase (KMT) that methylates histone 3 lysine-79 (H3K79) and was not previously associated with a Mendelian disease in OMIM. We have identified nine unrelated individuals with seven different de novo heterozygous missense variants in DOT1L through the Undiagnosed Disease Network (UDN), the SickKids Complex Care genomics project, and GeneMatcher. All probands had some degree of global developmental delay/intellectual disability, and most had one or more major congenital anomalies. To assess the pathogenicity of the DOT1L variants, functional studies were performed in Drosophila and human cells. The fruit fly DOT1L ortholog, grappa, is expressed in most cells including neurons in the central nervous system. The identified DOT1L variants behave as gain-of-function alleles in flies and lead to increased H3K79 methylation levels in flies and human cells. Our results show that human DOT1L and fly grappa are required for proper development and that de novo heterozygous variants in DOT1L are associated with a Mendelian disease.

A defect in mitochondrial fatty acid synthesis impairs iron metabolism and causes elevated ceramide levels

In most eukaryotic cells, fatty acid synthesis (FAS) occurs in the cytoplasm and in mitochondria. However, the relative contribution of mitochondrial FAS (mtFAS) to the cellular lipidome is not well defined. Here we show that loss of function of Drosophila mitochondrial enoyl coenzyme A reductase (Mecr), which is the enzyme required for the last step of mtFAS, causes lethality, while neuronal loss of Mecr leads to progressive neurodegeneration. We observe a defect in Fe-S cluster biogenesis and increased iron levels in flies lacking mecr, leading to elevated ceramide levels. Reducing the levels of either iron or ceramide suppresses the neurodegenerative phenotypes, indicating an interplay between ceramide and iron metabolism. Mutations in human MECR cause pediatric-onset neurodegeneration, and we show that human-derived fibroblasts display similar elevated ceramide levels and impaired iron homeostasis. In summary, this study identifies a role of mecr/MECR in ceramide and iron metabolism, providing a mechanistic link between mtFAS and neurodegeneration.

A de novo missense variant in EZH1 associated with developmental delay exhibits functional deficits in Drosophila melanogaster

EZH1, a polycomb repressive complex-2 component, is involved in a myriad of cellular processes. EZH1 represses transcription of downstream target genes through histone 3 lysine27 (H3K27) trimethylation (H3K27me3). Genetic variants in histone modifiers have been associated with developmental disorders, while EZH1 has not yet been linked to any human disease. However, the paralog EZH2 is associated with Weaver syndrome. Here we report a previously undiagnosed individual with a novel neurodevelopmental phenotype identified to have a de novo missense variant in EZH1 through exome sequencing. The individual presented in infancy with neurodevelopmental delay and hypotonia and was later noted to have proximal muscle weakness. The variant, p.A678G, is in the SET domain, known for its methyltransferase activity, and an analogous somatic or germline mutation in EZH2 has been reported in patients with B-cell lymphoma or Weaver syndrome, respectively. Human EZH1/2 are homologous to fly Enhancer of zeste (E(z)), an essential gene in Drosophila, and the affected residue (p.A678 in humans, p.A691 in flies) is conserved. To further study this variant, we obtained null alleles and generated transgenic flies expressing wildtype [E(z)WT] and the variant [E(z)A691G]. When expressed ubiquitously the variant rescues null-lethality similar to the wildtype. Overexpression of E(z)WT induces homeotic patterning defects but notably the E(z)A691G variant leads to dramatically stronger morphological phenotypes. We also note a dramatic loss of H3K27me2 and a corresponding increase in H3K27me3 in flies expressing E(z)A691G, suggesting this acts as a gain-of-function allele. In conclusion, here we present a novel EZH1 de novo variant associated with a neurodevelopmental disorder. Furthermore, we found that this variant has a functional impact in Drosophila.

Functional consequences of C-terminal mutations in RUNX2

Cleidocranial dysplasia (CCD) is a genetic disorder caused by mutations in the RUNX2 gene, affecting bone and teeth development. Previous studies focused on mutations in the RUNX2 RHD domain, with limited investigation of mutations in the C-terminal domain. This study aimed to investigate the functional consequences of C-terminal mutations in RUNX2. Eight mutations were analyzed, and their effects on transactivation activity, protein expression, subcellular localization, and osteogenic potential were studied. Truncating mutations in the PST region and a missense mutation in the NMTS region resulted in increased transactivation activity, while missense mutations in the PST showed activity comparable to the control. Truncating mutations produced truncated proteins, while missense mutations produced normal-sized proteins. Mutant proteins were mislocalized, with six mutant proteins detected in both the nucleus and cytoplasm. CCD patient bone cells exhibited mislocalization of RUNX2, similar to the generated mutant. Mislocalization of RUNX2 and reduced expression of downstream genes were observed in MSCs from a CCD patient with the p.Ser247Valfs*3 mutation, leading to compromised osteogenic potential. This study provides insight into the functional consequences of C-terminal mutations in RUNX2, including reduced expression, mislocalization, and aberrant transactivation of downstream genes, contributing to the compromised osteogenic potential observed in CCD.

PANGEA: a new gene set enrichment tool for Drosophila and common research organisms

Gene set enrichment analysis (GSEA) plays an important role in large-scale data analysis, helping scientists discover the underlying biological patterns over-represented in a gene list resulting from, for example, an 'omics' study. Gene Ontology (GO) annotation is the most frequently used classification mechanism for gene set definition. Here we present a new GSEA tool, PANGEA (PAthway, Network and Gene-set Enrichment Analysis; https://www.flyrnai.org/tools/pangea/), developed to allow a more flexible and configurable approach to data analysis using a variety of classification sets. PANGEA allows GO analysis to be performed on different sets of GO annotations, for example excluding high-throughput studies. Beyond GO, gene sets for pathway annotation and protein complex data from various resources as well as expression and disease annotation from the Alliance of Genome Resources (Alliance). In addition, visualizations of results are enhanced by providing an option to view network of gene set to gene relationships. The tool also allows comparison of multiple input gene lists and accompanying visualisation tools for quick and easy comparison. This new tool will facilitate GSEA for Drosophila and other major model organisms based on high-quality annotated information available for these species.

Unlocking the Mysteries of Alpha-N-Terminal Methylation and its Diverse Regulatory Functions

Protein post-translation modifications (PTMs) are a critical regulatory mechanism of protein function. Protein α-N-terminal (Nα) methylation is a conserved PTM across prokaryotes and eukaryotes. Studies of the Nα methyltransferases responsible for Να methylation and their substrate proteins have shown that the PTM involves diverse biological processes, including protein synthesis and degradation, cell division, DNA damage response, and transcription regulation. This review provides an overview of the progress toward the regulatory function of Να methyltransferases and their substrate landscape. More than 200 proteins in humans and 45 in yeast are potential substrates for protein Nα methylation based on the canonical recognition motif, XP[KR]. Based on recent evidence for a less stringent motif requirement, the number of substrates might be increased, but further validation is needed to solidify this concept. A comparison of the motif in substrate orthologs in selected eukaryotic species indicates intriguing gain and loss of the motif across the evolutionary landscape. We discuss the state of knowledge in the field that has provided insights into the regulation of protein Να methyltransferases and their role in cellular physiology and disease. We also outline the current research tools that are key to understanding Να methylation. Finally, challenges are identified and discussed that would aid in unlocking a system-level view of the roles of Να methylation in diverse cellular pathways.

Bi-allelic variants in INTS11 are associated with a complex neurological disorder

The Integrator complex is a multi-subunit protein complex that regulates the processing of nascent RNAs transcribed by RNA polymerase II (RNAPII), including small nuclear RNAs, enhancer RNAs, telomeric RNAs, viral RNAs, and protein-coding mRNAs. Integrator subunit 11 (INTS11) is the catalytic subunit that cleaves nascent RNAs, but, to date, mutations in this subunit have not been linked to human disease. Here, we describe 15 individuals from 10 unrelated families with bi-allelic variants in INTS11 who present with global developmental and language delay, intellectual disability, impaired motor development, and brain atrophy. Consistent with human observations, we find that the fly ortholog of INTS11, dIntS11, is essential and expressed in the central nervous systems in a subset of neurons and most glia in larval and adult stages. Using Drosophila as a model, we investigated the effect of seven variants. We found that two (p.Arg17Leu and p.His414Tyr) fail to rescue the lethality of null mutants, indicating that they are strong loss-of-function variants. Furthermore, we found that five variants (p.Gly55Ser, p.Leu138Phe, p.Lys396Glu, p.Val517Met, and p.Ile553Glu) rescue lethality but cause a shortened lifespan and bang sensitivity and affect locomotor activity, indicating that they are partial loss-of-function variants. Altogether, our results provide compelling evidence that integrity of the Integrator RNA endonuclease is critical for brain development.

Very-long-chain fatty acids induce glial-derived sphingosine-1-phosphate synthesis, secretion, and neuroinflammation

VLCFAs (very-long-chain fatty acids) are the most abundant fatty acids in myelin. Hence, during demyelination or aging, glia are exposed to higher levels of VLCFA than normal. We report that glia convert these VLCFA into sphingosine-1-phosphate (S1P) via a glial-specific S1P pathway. Excess S1P causes neuroinflammation, NF-κB activation, and macrophage infiltration into the CNS. Suppressing the function of S1P in fly glia or neurons, or administration of Fingolimod, an S1P receptor antagonist, strongly attenuates the phenotypes caused by excess VLCFAs. In contrast, elevating the VLCFA levels in glia and immune cells exacerbates these phenotypes. Elevated VLCFA and S1P are also toxic in vertebrates based on a mouse model of multiple sclerosis (MS), experimental autoimmune encephalomyelitis (EAE). Indeed, reducing VLCFA with bezafibrate ameliorates the phenotypes. Moreover, simultaneous use of bezafibrate and fingolimod synergizes to improve EAE, suggesting that lowering VLCFA and S1P is a treatment avenue for MS.

Lateral thinking in syndromic congenital cardiovascular disease

Syndromic birth defects are rare diseases that can present with seemingly pleiotropic comorbidities. Prime examples are rare congenital heart and cardiovascular anomalies that can be accompanied by forelimb defects, kidney disorders and more. Whether such multi-organ defects share a developmental link remains a key question with relevance to the diagnosis, therapeutic intervention and long-term care of affected patients. The heart, endothelial and blood lineages develop together from the lateral plate mesoderm (LPM), which also harbors the progenitor cells for limb connective tissue, kidneys, mesothelia and smooth muscle. This developmental plasticity of the LPM, which founds on multi-lineage progenitor cells and shared transcription factor expression across different descendant lineages, has the potential to explain the seemingly disparate syndromic defects in rare congenital diseases. Combining patient genome-sequencing data with model organism studies has already provided a wealth of insights into complex LPM-associated birth defects, such as heart-hand syndromes. Here, we summarize developmental and known disease-causing mechanisms in early LPM patterning, address how defects in these processes drive multi-organ comorbidities, and outline how several cardiovascular and hematopoietic birth defects with complex comorbidities may be LPM-associated diseases. We also discuss strategies to integrate patient sequencing, data-aggregating resources and model organism studies to mechanistically decode congenital defects, including potentially LPM-associated orphan diseases. Eventually, linking complex congenital phenotypes to a common LPM origin provides a framework to discover developmental mechanisms and to anticipate comorbidities in congenital diseases affecting the cardiovascular system and beyond.

PANGEA: A New Gene Set Enrichment Tool for Drosophila and Common Research Organisms

Gene set enrichment analysis (GSEA) plays an important role in large-scale data analysis, helping scientists discover the underlying biological patterns over-represented in a gene list resulting from, for example, an 'omics' study. Gene Ontology (GO) annotation is the most frequently used classification mechanism for gene set definition. Here we present a new GSEA tool, PANGEA (PAthway, Network and Gene-set Enrichment Analysis; https://www.flyrnai.org/tools/pangea/ ), developed to allow a more flexible and configurable approach to data analysis using a variety of classification sets. PANGEA allows GO analysis to be performed on different sets of GO annotations, for example excluding high-throughput studies. Beyond GO, gene sets for pathway annotation and protein complex data from various resources as well as expression and disease annotation from the Alliance of Genome Resources (Alliance). In addition, visualisations of results are enhanced by providing an option to view network of gene set to gene relationships. The tool also allows comparison of multiple input gene lists and accompanying visualisation tools for quick and easy comparison. This new tool will facilitate GSEA for Drosophila and other major model organisms based on high-quality annotated information available for these species.

A retrospective analysis of phosphatase catalytic subunit gene variants in patients with rare disorders identifies novel candidate neurodevelopmental disease genes

Rare genetic disorders represent some of the most severe and life-limiting conditions that constitute a considerable burden on global healthcare systems and societies. Most individuals affected by rare disorders remain undiagnosed, highlighting the unmet need for improved disease gene discovery and novel variant interpretation. Aberrant (de) phosphorylation can have profound pathological consequences underpinning many disease processes. Numerous phosphatases and associated proteins have been identified as disease genes, with many more likely to have gone undiscovered thus far. To begin to address these issues, we have performed a systematic survey of de novo variants amongst 189 genes encoding phosphatase catalytic subunits found in rare disease patients recruited to the 100,000 Genomes Project (100 kGP), the largest national sequencing project of its kind in the United Kingdom. We found that 49% of phosphatases were found to carry de novo mutation(s) in this cohort. Only 25% of these phosphatases have been previously linked to genetic disorders. A gene-to-patient approach matching variants to phenotypic data identified 9 novel candidate rare-disease genes: PTPRD, PTPRG, PTPRT, PTPRU, PTPRZ1, MTMR3, GAK, TPTE2, PTPN18. As the number of patients undergoing whole genome sequencing increases and information sharing improves, we anticipate that reiterative analysis of genomic and phenotypic data will continue to identify candidate phosphatase disease genes for functional validation. This is the first step towards delineating the aetiology of rare genetic disorders associated with altered phosphatase function, leading to new biological insights and improved clinical outcomes for the affected individuals and their families.

The fly homolog of SUPT16H, a gene associated with neurodevelopmental disorders, is required in a cell-autonomous fashion for cell survival

SUPT16H encodes the large subunit of the FAcilitate Chromatin Transcription (FACT) complex, which functions as a nucleosome organizer during transcription. We identified two individuals from unrelated families carrying de novo missense variants in SUPT16H. The probands exhibit global developmental delay, intellectual disability, epilepsy, facial dysmorphism and brain structural abnormalities. We used Drosophila to characterize two variants: p.T171I and p.G808R. Loss of the fly ortholog, dre4, causes lethality at an early developmental stage. RNAi-mediated knockdown of dre4 in either glia or neurons causes severely reduced eclosion and longevity. Tissue-specific knockdown of dre4 in the eye or wing leads to the loss of these tissues, whereas overexpression of SUPT16H has no dominant effect. Moreover, expression of the reference SUPT16H significantly rescues the loss-of-function phenotypes in the nervous system as well as wing and eye. In contrast, expression of SUPT16H p.T171I or p.G808R rescues the phenotypes poorly, indicating that the variants are partial loss-of-function alleles. While previous studies argued that the developmental arrest caused by loss of dre4 is due to impaired ecdysone production in the prothoracic gland, our data show that dre4 is required for proper cell growth and survival in multiple tissues in a cell-autonomous manner. Altogether, our data indicate that the de novo loss-of-function variants in SUPT16H are indeed associated with developmental and neurological defects observed in the probands.

Bi-allelic ATG4D variants are associated with a neurodevelopmental disorder characterized by speech and motor impairment

Autophagy regulates the degradation of damaged organelles and protein aggregates, and is critical for neuronal development, homeostasis, and maintenance, yet few neurodevelopmental disorders have been associated with pathogenic variants in genes encoding autophagy-related proteins. We report three individuals from two unrelated families with a neurodevelopmental disorder characterized by speech and motor impairment, and similar facial characteristics. Rare, conserved, bi-allelic variants were identified in ATG4D, encoding one of four ATG4 cysteine proteases important for autophagosome biogenesis, a hallmark of autophagy. Autophagosome biogenesis and induction of autophagy were intact in cells from affected individuals. However, studies evaluating the predominant substrate of ATG4D, GABARAPL1, demonstrated that three of the four ATG4D patient variants functionally impair ATG4D activity. GABARAPL1 is cleaved or "primed" by ATG4D and an in vitro GABARAPL1 priming assay revealed decreased priming activity for three of the four ATG4D variants. Furthermore, a rescue experiment performed in an ATG4 tetra knockout cell line, in which all four ATG4 isoforms were knocked out by gene editing, showed decreased GABARAPL1 priming activity for the two ATG4D missense variants located in the cysteine protease domain required for priming, suggesting that these variants impair the function of ATG4D. The clinical, bioinformatic, and functional data suggest that bi-allelic loss-of-function variants in ATG4D contribute to the pathogenesis of this syndromic neurodevelopmental disorder.

A de novo missense variant in EZH1 associated with developmental delay exhibits functional deficits in Drosophila melanogaster

EZH1 ( Enhancer of Zeste, homolog 1) , a Polycomb Repressive Complex-2 (PRC2) component, is involved in a myriad of cellular processes through modifying histone 3 lysine27 (H3K27) residues. EZH1 represses transcription of downstream target genes through H3K27 trimethylation (H3K27me3). Genetic mutations in histone modifiers have been associated with developmental disorders, while EZH1 has not yet been linked to any human disease. However, the paralog EZH2 is associated with Weaver syndrome. Here we report a previously undiagnosed individual with a novel neurodevelopmental phenotype identified to have a de novo variant in EZH1 , p.Ala678Gly, through exome sequencing. The individual presented in infancy with neurodevelopmental delay and hypotonia and was later noted to have proximal muscle weakness. The variant, p.A678G, is in the SET domain, known for its methyltransferase activity, and was the best candidate variant found in the exome. Human EZH1 / 2 are homologous to fly Enhancer of zeste E(z) , an essential gene in flies, and the residue (A678 in humans, A691 in Drosophila ) is conserved. To further study this variant, we obtained Drosophila null alleles and generated transgenic flies expressing wild-type (E(z) WT ) and the variant (E(z) A691G ) . The E(z) A691G variant led to hyper H3K27me3 while the E(z) WT did not, suggesting this is as a gain-of-function allele. When expressed under the tubulin promotor in vivo the variant rescued null-lethality similar to wild-type but the E(z) A691G flies exhibit bang sensitivity and shortened lifespan. In conclusion, here we present a novel EZH1 de novo variant associated with a neurodevelopmental disorder. Furthermore, we found that this variant has a functional impact in Drosophila . Biochemically this allele leads to increased H3K27me3 suggesting gain-of-function, but when expressed in adult flies the E(z) A691G has some characteristics of partial loss-of-function which may suggest it is a more complex allele in vivo .

High-throughput sequencing analysis of nuclear-encoded mitochondrial genes reveals a genetic signature of human longevity

Mitochondrial dysfunction is a well-known contributor to aging and age-related diseases. The precise mechanisms through which mitochondria impact human lifespan, however, remain unclear. We hypothesize that humans with exceptional longevity harbor rare variants in nuclear-encoded mitochondrial genes (mitonuclear genes) that confer resistance against age-related mitochondrial dysfunction. Here we report an integrated functional genomics study to identify rare functional variants in ~ 660 mitonuclear candidate genes discovered by target capture sequencing analysis of 496 centenarians and 572 controls of Ashkenazi Jewish descent. We identify and prioritize longevity-associated variants, genes, and mitochondrial pathways that are enriched with rare variants. We provide functional gene variants such as those in MTOR (Y2396Lfs*29), CPS1 (T1406N), and MFN2 (G548*) as well as LRPPRC (S1378G) that is predicted to affect mitochondrial translation. Taken together, our results suggest a functional role for specific mitonuclear genes and pathways in human longevity.

Evolutionary Landscape of SOX Genes to Inform Genotype-to-Phenotype Relationships

The SOX transcription factor family is pivotal in controlling aspects of development. To identify genotype-phenotype relationships of SOX proteins, we performed a non-biased study of SOX using 1890 open-reading frame and 6667 amino acid sequences in combination with structural dynamics to interpret 3999 gnomAD, 485 ClinVar, 1174 Geno2MP, and 4313 COSMIC human variants. We identified, within the HMG (High Mobility Group)- box, twenty-seven amino acids with changes in multiple SOX proteins annotated to clinical pathologies. These sites were screened through Geno2MP medical phenotypes, revealing novel SOX15 R104G associated with musculature abnormality and SOX8 R159G with intellectual disability. Within gnomAD, SOX18 E137K (rs201931544), found within the HMG box of ~0.8% of Latinx individuals, is associated with seizures and neurological complications, potentially through blood-brain barrier alterations. A total of 56 highly conserved variants were found at sites outside the HMG-box, including several within the SOX2 HMG-box-flanking region with neurological associations, several in the SOX9 dimerization region associated with Campomelic Dysplasia, SOX14 K88R (rs199932938) flanking the HMG box associated with cardiovascular complications within European populations, and SOX7 A379V (rs143587868) within an SOXF conserved far C-terminal domain heterozygous in 0.716% of African individuals with associated eye phenotypes. This SOX data compilation builds a robust genotype-to-phenotype association for a gene family through more robust ortholog data integration.

Variant-to-gene mapping followed by cross-species genetic screening identifies GPI-anchor biosynthesis as a regulator of sleep

Genome-wide association studies (GWAS) in humans have identified loci robustly associated with several heritable diseases or traits, yet little is known about the functional roles of the underlying causal variants in regulating sleep duration or quality. We applied an ATAC-seq/promoter focused Capture C strategy in human iPSC-derived neural progenitors to carry out a "variant-to-gene" mapping campaign that identified 88 candidate sleep effector genes connected to relevant GWAS signals. To functionally validate the role of the implicated effector genes in sleep regulation, we performed a neuron-specific RNA interference screen in the fruit fly, Drosophila melanogaster, followed by validation in zebrafish. This approach identified a number of genes that regulate sleep including a critical role for glycosylphosphatidylinositol (GPI)-anchor biosynthesis. These results provide the first physical variant-to-gene mapping of human sleep genes followed by a model organism-based prioritization, revealing a conserved role for GPI-anchor biosynthesis in sleep regulation.

Drosophila melanogaster as a versatile model organism to study genetic epilepsies: An overview

Epilepsy is one of the most prevalent neurological disorders, affecting more than 45 million people worldwide. Recent advances in genetic techniques, such as next-generation sequencing, have driven genetic discovery and increased our understanding of the molecular and cellular mechanisms behind many epilepsy syndromes. These insights prompt the development of personalized therapies tailored to the genetic characteristics of an individual patient. However, the surging number of novel genetic variants renders the interpretation of pathogenetic consequences and of potential therapeutic implications ever more challenging. Model organisms can help explore these aspects in vivo. In the last decades, rodent models have significantly contributed to our understanding of genetic epilepsies but their establishment is laborious, expensive, and time-consuming. Additional model organisms to investigate disease variants on a large scale would be desirable. The fruit fly Drosophila melanogaster has been used as a model organism in epilepsy research since the discovery of "bang-sensitive" mutants more than half a century ago. These flies respond to mechanical stimulation, such as a brief vortex, with stereotypic seizures and paralysis. Furthermore, the identification of seizure-suppressor mutations allows to pinpoint novel therapeutic targets. Gene editing techniques, such as CRISPR/Cas9, are a convenient way to generate flies carrying disease-associated variants. These flies can be screened for phenotypic and behavioral abnormalities, shifting of seizure thresholds, and response to anti-seizure medications and other substances. Moreover, modification of neuronal activity and seizure induction can be achieved using optogenetic tools. In combination with calcium and fluorescent imaging, functional alterations caused by mutations in epilepsy genes can be traced. Here, we review Drosophila as a versatile model organism to study genetic epilepsies, especially as 81% of human epilepsy genes have an orthologous gene in Drosophila. Furthermore, we discuss newly established analysis techniques that might be used to further unravel the pathophysiological aspects of genetic epilepsies.

The impact of variance in carnitine palmitoyltransferase-1 expression on breast cancer prognosis is stratified by clinical and anthropometric factors

CPT1A is a rate-limiting enzyme in fatty acid oxidation and is upregulated in high-risk breast cancer. Obesity and menopausal status' relationship with breast cancer prognosis is well established, but its connection with fatty acid metabolism is not. We utilized RNA sequencing data in the Xena Functional Genomics Explorer, to explore CPT1A's effect on breast cancer patients' survival probability. Using [18F]-fluorothymidine positron emission tomography-computed tomography images from The Cancer Imaging Archive, we segmented these analyses by obesity and menopausal status. In 1214 patients, higher CPT1A expression is associated with lower breast cancer survivability. We confirmed a previously observed protective relationship between obesity and breast cancer in pre-menopausal patients and supported this data using two-sided Pearson correlations. Taken together, these analyses using open-access databases bolster the potential role of CPT1A-dependent fatty acid metabolism as a pathogenic factor in breast cancer.

MetGENE: gene-centric metabolomics information retrieval tool

BACKGROUND

Biomedical research often involves contextual integration of multimodal and multiomic data in search of mechanisms for improved diagnosis, treatment, and monitoring. Researchers need to access information from diverse sources, comprising data in various and sometimes incongruent formats. The downstream processing of the data to decipher mechanisms by reconstructing networks and developing quantitative models warrants considerable effort.

RESULTS

MetGENE is a knowledge-based, gene-centric data aggregator that hierarchically retrieves information about the gene(s), their related pathway(s), reaction(s), metabolite(s), and metabolomic studies from standard data repositories under one dashboard to enable ease of access through centralization of relevant information. We note that MetGENE focuses only on those genes that encode for proteins directly associated with metabolites. All other gene-metabolite associations are beyond the current scope of MetGENE. Further, the information can be contextualized by filtering by species, anatomy (tissue), and condition (disease or phenotype).

CONCLUSIONS

MetGENE is an open-source tool that aggregates metabolite information for a given gene(s) and presents them in different computable formats (e.g., JSON) for further integration with other omics studies. MetGENE is available at https://bdcw.org/MetGENE/index.php.

Pharos 2023: an integrated resource for the understudied human proteome

The Illuminating the Druggable Genome (IDG) project aims to improve our understanding of understudied proteins and our ability to study them in the context of disease biology by perturbing them with small molecules, biologics, or other therapeutic modalities. Two main products from the IDG effort are the Target Central Resource Database (TCRD) (http://juniper.health.unm.edu/tcrd/), which curates and aggregates information, and Pharos (https://pharos.nih.gov/), a web interface for fusers to extract and visualize data from TCRD. Since the 2021 release, TCRD/Pharos has focused on developing visualization and analysis tools that help reveal higher-level patterns in the underlying data. The current iterations of TCRD and Pharos enable users to perform enrichment calculations based on subsets of targets, diseases, or ligands and to create interactive heat maps and UpSet charts of many types of annotations. Using several examples, we show how to address disease biology and drug discovery questions through enrichment calculations and UpSet charts.

RG/RGG repeats in the C. elegans homologs of Nucleolin and GAR1 contribute to sub-nucleolar phase separation

The intrinsically disordered RG/RGG repeat domain is found in several nucleolar and P-granule proteins, but how it influences their phase separation into biomolecular condensates is unclear. We survey all RG/RGG repeats in C. elegans and uncover nucleolar and P-granule-specific RG/RGG motifs. An uncharacterized protein, K07H8.10, contains the longest nucleolar-like RG/RGG domain in C. elegans. Domain and sequence similarity, as well as nucleolar localization, reveals K07H8.10 (NUCL-1) to be the homolog of Nucleolin, a protein conserved across animals, plants, and fungi, but previously thought to be absent in nematodes. Deleting the RG/RGG repeats within endogenous NUCL-1 and a second nucleolar protein, GARR-1 (GAR1), demonstrates these domains are dispensable for nucleolar accumulation. Instead, their RG/RGG repeats contribute to the phase separation of proteins into nucleolar sub-compartments. Despite this common RG/RGG repeat function, only removal of the GARR-1 RG/RGG domain affects worm fertility and development, decoupling precise sub-nucleolar structure from nucleolar function.

ConVarT: Search Engine for Missense Variants Between Humans and Other Organisms

ConVarT (https://convart.org/) is a search engine for searching for conjugate variants between humans and other species. The search engine is based on matching conjugate variants called MatchVars between species. Matching equivalent variants requires correct alignment of orthologous proteins with the use of multiple sequence alignments (MSA). Indeed, the ConVarT pipeline has performed over a million MSAs and integrated variants and variant-specific annotations (pathogenicity, phenotypic variants; etc.) into the corresponding positions on MSAs. When a clinically relevant variant is discovered whose functional relevance is unknown, ConVarT offers clinician scientists the possibility to search for a MatchVar in other species and to look for functional data on that variant. Fortunately, ConVarT enables users to paste a protein sequence in FASTA format to search for human orthologous proteins. A pairwise sequence alignment (PSA) is then performed between the provided protein sequence and the human orthologous protein, allowing users to visualize human variants on the PSA. Here, we describe the step-by-step usage of ConVarT. © 2022 Wiley Periodicals LLC. Basic Protocol 1: Searching matching variants (MatchVar) with gene/protein identifiers. Basic Protocol 2: Searching with a FASTA sequence. Alternate Protocol: Search with gene name in multiple species. Basic Protocol 3: Search genes associated with a disease.

De novo variants in FRMD5 are associated with developmental delay, intellectual disability, ataxia, and abnormalities of eye movement

Proteins containing the FERM (four-point-one, ezrin, radixin, and moesin) domain link the plasma membrane with cytoskeletal structures at specific cellular locations and have been implicated in the localization of cell-membrane-associated proteins and/or phosphoinositides. FERM domain-containing protein 5 (FRMD5) localizes at cell adherens junctions and stabilizes cell-cell contacts. To date, variants in FRMD5 have not been associated with a Mendelian disease in OMIM. Here, we describe eight probands with rare heterozygous missense variants in FRMD5 who present with developmental delay, intellectual disability, ataxia, seizures, and abnormalities of eye movement. The variants are de novo in all for whom parental testing was available (six out of eight probands), and human genetic datasets suggest that FRMD5 is intolerant to loss of function (LoF). We found that the fly ortholog of FRMD5, CG5022 (dFrmd), is expressed in the larval and adult central nervous systems where it is present in neurons but not in glia. dFrmd LoF mutant flies are viable but are extremely sensitive to heat shock, which induces severe seizures. The mutants also exhibit defective responses to light. The human FRMD5 reference (Ref) cDNA rescues the fly dFrmd LoF phenotypes. In contrast, all the FRMD5 variants tested in this study (c.340T>C, c.1051A>G, c.1053C>G, c.1054T>C, c.1045A>C, and c.1637A>G) behave as partial LoF variants. In addition, our results indicate that two variants that were tested have dominant-negative effects. In summary, the evidence supports that the observed variants in FRMD5 cause neurological symptoms in humans.

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.

Case Report of a Juvenile Patient with Autism Spectrum Disorder with a Novel Combination of Copy Number Variants in ADGRL3 (LPHN3) and Two Pseudogenes

We report the finding of two copy number variants (CNVs) in a 12-year-old boy presenting both with autism spectrum disorder (ASD) and attention deficit/hyperactivity disorder (ADHD). Clinical features included aggressive behavior, mood instability, suicidal statements, repetitive and restrictive behavior, sensitivity to noise, learning problems and dyslexia, though no intellectual disability was present. Using array-based comparative genomic hybridization (array-CGH), we identified two CNVs, both triplex duplications of 324 kb on 3p26.3, and 284 kb on 4q13.1, respectively. One of the CNVs is located on chromosome 4q13.1 in the region of the gene encoding for adhesion G protein-coupled receptor L3 (ADGRL3, former name: latrophilin-3, LPHN3), the other on chromosome 3p26.3 in the region of the two pseudogenes AC090043.1 and RPL23AP39. The patient described in the present study showed increased symptoms under methylphenidate treatment but responded positively to 3 mg per day of the atypical neuroleptic drug aripiprazole. To our knowledge, this is the first report of a CNV in the ADGRL3 gene and its first association with ASD in humans.

Promoting validation and cross-phylogenetic integration in model organism research

Model organism (MO) research provides a basic understanding of biology and disease due to the evolutionary conservation of the molecular and cellular language of life. MOs have been used to identify and understand the function of orthologous genes, proteins, cells and tissues involved in biological processes, to develop and evaluate techniques and methods, and to perform whole-organism-based chemical screens to test drug efficacy and toxicity. However, a growing richness of datasets and the rising power of computation raise an important question: How do we maximize the value of MOs? In-depth discussions in over 50 virtual presentations organized by the National Institutes of Health across more than 10 weeks yielded important suggestions for improving the rigor, validation, reproducibility and translatability of MO research. The effort clarified challenges and opportunities for developing and integrating tools and resources. Maintenance of critical existing infrastructure and the implementation of suggested improvements will play important roles in maintaining productivity and facilitating the validation of animal models of human biology and disease.

'Fly-ing' from rare to common neurodegenerative disease mechanisms

Advances in genome sequencing have enabled researchers and clinicians to probe vast numbers of human variants to distinguish pathogenic from benign variants. Model organisms have been crucial in variant assessment and in delineating the molecular mechanisms of some of the diseases caused by these variants. The fruit fly, Drosophila melanogaster, has played a valuable role in this endeavor, taking advantage of its genetic technologies and established biological knowledge. We highlight the utility of the fly in studying the function of genes associated with rare neurological diseases that have led to a better understanding of common disease mechanisms. We emphasize that shared themes emerge among disease mechanisms, including the importance of lipids, in two prominent neurodegenerative diseases: Alzheimer's disease (AD) and Parkinson's disease (PD).

The microRNA processor DROSHA is a candidate gene for a severe progressive neurological disorder

DROSHA encodes a ribonuclease that is a subunit of the Microprocessor complex and is involved in the first step of microRNA (miRNA) biogenesis. To date, DROSHA has not yet been associated with a Mendelian disease. Here, we describe two individuals with profound intellectual disability, epilepsy, white matter atrophy, microcephaly and dysmorphic features, who carry damaging de novo heterozygous variants in DROSHA. DROSHA is constrained for missense variants and moderately intolerant to loss-of-function (o/e = 0.24). The loss of the fruit fly ortholog drosha causes developmental arrest and death in third instar larvae, a severe reduction in brain size and loss of imaginal discs in the larva. Loss of drosha in eye clones causes small and rough eyes in adult flies. One of the identified DROSHA variants (p.Asp1219Gly) behaves as a strong loss-of-function allele in flies, while another variant (p.Arg1342Trp) is less damaging in our assays. In worms, a knock-in that mimics the p.Asp1219Gly variant at a worm equivalent residue causes loss of miRNA expression and heterochronicity, a phenotype characteristic of the loss of miRNA. Together, our data show that the DROSHA variants found in the individuals presented here are damaging based on functional studies in model organisms and likely underlie the severe phenotype involving the nervous system.

Novel dominant and recessive variants in human ROBO1 cause distinct neurodevelopmental defects through different mechanisms

The Roundabout (Robo) receptors, located on growth cones of neurons, induce axon repulsion in response to the extracellular ligand Slit. The Robo family of proteins controls midline crossing of commissural neurons during development in flies. Mono- and bi-allelic variants in human ROBO1 (HGNC: 10249) have been associated with incomplete penetrance and variable expressivity for a breath of phenotypes, including neurodevelopmental defects such as strabismus, pituitary defects, intellectual impairment, as well as defects in heart and kidney. Here, we report two novel ROBO1 variants associated with very distinct phenotypes. A homozygous missense p.S1522L variant in three affected siblings with nystagmus; and a monoallelic de novo p.D422G variant in a proband who presented with early-onset epileptic encephalopathy. We modeled these variants in Drosophila and first generated a null allele by inserting a CRIMIC T2A-GAL4 in an intron. Flies that lack robo1 exhibit reduced viability but have very severe midline crossing defects in the central nervous system. The fly wild-type cDNA driven by T2A-Gal4 partially rescues both defects. Overexpression of the human reference ROBO1 with T2A-GAL4 is toxic and reduces viability, whereas the recessive p.S1522L variant is less toxic, suggesting that it is a partial loss-of-function allele. In contrast, the dominant variant in fly robo1 (p.D413G) affects protein localization, impairs axonal guidance activity and induces mild phototransduction defects, suggesting that it is a neomorphic allele. In summary, our studies expand the phenotypic spectrum associated with ROBO1 variant alleles.

FertilityOnline: A Straightforward Pipeline for Functional Gene Annotation and Disease Mutation Discovery

Exploring the genetic basis of human infertility is currently under intensive investigation. However, only a handful of genes have been validated in animal models as disease-causing genes in infertile men. Thus, to better understand the genetic basis of human spermatogenesis and bridge the knowledge gap between humans and other animal species, we construct the FertilityOnline, a database integrating the literature-curated functional genes during spermatogenesis into an existing spermatogenic database, SpermatogenesisOnline 1.0. Additional features, including the functional annotation and genetic variants of human genes, are also incorporated into FertilityOnline. By searching this database, users can browse the functional genes involved in spermatogenesis and instantly narrow down the number of candidates of genetic mutations underlying male infertility in a user-friendly web interface. Clinical application of this database was exampled by the identification of novel causative mutations in synaptonemal complex central element protein 1 (SYCE1) and stromal antigen 3 (STAG3) in azoospermic men. In conclusion, FertilityOnline is not only an integrated resource for spermatogenic genes but also a useful tool facilitating the exploration of the genetic basis of male infertility. FertilityOnline can be freely accessed at http://mcg.ustc.edu.cn/bsc/spermgenes2.0/index.html.

CTH exacerbates Helicobacter pylori immunopathogenesis by promoting macrophage metabolic remodeling and activation

Macrophages play a crucial role in the inflammatory response to the human stomach pathogen Helicobacter pylori, which infects half of the world's population and causes gastric cancer. Recent studies have highlighted the importance of macrophage immunometabolism in their activation state and function. We have demonstrated that the cysteine-producing enzyme, cystathionine g-lyase (CTH), is upregulated in humans and mice with H. pylori infection. Here we show that induction of CTH in macrophages by H. pylori promotes persistent inflammation. Cth-/- mice have reduced macrophage and T-cell activation in H. pylori-infected tissues, an altered metabolome, and decreased enrichment of immune-associated gene networks, culminating in decreased H. pylori-induced-gastritis. CTH is downstream of the proposed anti-inflammatory molecule, S-adenosylmethionine (SAM). While Cth-/- mice exhibit gastric SAM accumulation, WT mice treated with SAM did not display protection against H. pylori-induced inflammation. Instead, we demonstrate that Cth-deficient macrophages exhibit alterations in the proteome, decreased NF-kB activation, diminished expression of macrophage activation markers, and impaired oxidative phosphorylation and glycolysis. Thus, through altering cellular respiration, CTH is a key enhancer of macrophage activation contributing to a pathogenic inflammatory response that is the universal precursor for the development of H. pylori-induced gastric disease.

Comprehensive microRNA-seq transcriptomic profiling across 11 organs, 4 ages, and 2 sexes of Fischer 344 rats

Rat is one of the most widely-used models in chemical safety evaluation and biomedical research. However, the knowledge about its microRNA (miRNA) expression patterns across multiple organs and various developmental stages is still limited. Here, we constructed a comprehensive rat miRNA expression BodyMap using a diverse collection of 320 RNA samples from 11 organs of both sexes of juvenile, adolescent, adult and aged Fischer 344 rats with four biological replicates per group. Following the Illumina TruSeq Small RNA protocol, an average of 5.1 million 50 bp single-end reads was generated per sample, yielding a total of 1.6 billion reads. The quality of the resulting miRNA-seq data was deemed to be high from raw sequences, mapped sequences, and biological reproducibility. Importantly, aliquots of the same RNA samples have previously been used to construct the mRNA BodyMap. The currently presented miRNA-seq dataset along with the existing mRNA-seq dataset from the same RNA samples provides a unique resource for studying the expression characteristics of existing and novel miRNAs, and for integrative analysis of miRNA-mRNA interactions, thereby facilitating better utilization of rats for biomarker discovery.

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.

Disease Modeling of Rare Neurological Disorders in Zebrafish

Rare diseases are those which affect a small number of people compared to the general population. However, many patients with a rare disease remain undiagnosed, and a large majority of rare diseases still have no form of viable treatment. Approximately 40% of rare diseases include neurologic and neurodevelopmental disorders. In order to understand the characteristics of rare neurological disorders and identify causative genes, various model organisms have been utilized extensively. In this review, the characteristics of model organisms, such as roundworms, fruit flies, and zebrafish, are examined, with an emphasis on zebrafish disease modeling in rare neurological disorders.