An integrative approach to ortholog prediction for disease-focused and other functional studies

Memory consolidation in honey bees is enhanced by down-regulation of Down syndrome cell adhesion molecule and changes its alternative splicing

Down syndrome cell adhesion molecule (Dscam) gene encodes a cell adhesion molecule required for neuronal wiring. A remarkable feature of arthropod Dscam is massive alternative splicing generating thousands of different isoforms from three variable clusters of alternative exons. Dscam expression and diversity arising from alternative splicing have been studied during development, but whether they exert functions in adult brains has not been determined. Here, using honey bees, we find that Dscam expression is critically linked to memory retention as reducing expression by RNAi enhances memory after reward learning in adult worker honey bees. Moreover, alternative splicing of Dscam is altered in all three variable clusters after learning. Since identical Dscam isoforms engage in homophilic interactions, these results suggest a mechanism to alter inclusion of variable exons during memory consolidation to modify neuronal connections for memory retention.

Lipophorin receptors genetically modulate neurodegeneration caused by reduction of Psn expression in the aging Drosophila brain

Mutations in the Presenilin (PSEN) genes are the most common cause of early-onset familial Alzheimer's disease (FAD). Studies in cell culture, in vitro biochemical systems, and knockin mice showed that PSEN mutations are loss-of-function mutations, impairing γ-secretase activity. Mouse genetic analysis highlighted the importance of Presenilin (PS) in learning and memory, synaptic plasticity and neurotransmitter release, and neuronal survival, and Drosophila studies further demonstrated an evolutionarily conserved role of PS in neuronal survival during aging. However, molecular pathways that interact with PS in neuronal survival remain unclear. To identify genetic modifiers that modulate PS-dependent neuronal survival, we developed a new DrosophilaPsn model that exhibits age-dependent neurodegeneration and increases of apoptosis. Following a bioinformatic analysis, we tested top ranked candidate genes by selective knockdown (KD) of each gene in neurons using two independent RNAi lines in Psn KD models. Interestingly, 4 of the 9 genes enhancing neurodegeneration in Psn KD flies are involved in lipid transport and metabolism. Specifically, neuron-specific KD of lipophorin receptors, lpr1 and lpr2, dramatically worsens neurodegeneration in Psn KD flies, and overexpression of lpr1 or lpr2 does not alleviate Psn KD-induced neurodegeneration. Furthermore, lpr1 or lpr2 KD alone also leads to neurodegeneration, increased apoptosis, climbing defects, and shortened lifespan. Lastly, heterozygotic deletions of lpr1 and lpr2 or homozygotic deletions of lpr1 or lpr2 similarly lead to age-dependent neurodegeneration and further exacerbate neurodegeneration in Psn KD flies. These findings show that LpRs modulate Psn-dependent neuronal survival and are critically important for neuronal integrity in the aging brain.

Genome-wide association in Drosophila identifies a role for Piezo and Proc-R in sleep latency

Sleep latency, the amount of time that it takes an individual to fall asleep, is a key indicator of sleep need. Sleep latency varies considerably both among and within species and is heritable, but lacks a comprehensive description of its underlying genetic network. Here we conduct a genome-wide association study of sleep latency. Using previously collected sleep and activity data on a wild-derived population of flies, we calculate sleep latency, confirming significant, heritable genetic variation for this complex trait. We identify 520 polymorphisms in 248 genes contributing to variability in sleep latency. Tests of mutations in 23 candidate genes and additional putative pan-neuronal knockdown of 9 of them implicated CG44153, Piezo, Proc-R and Rbp6 in sleep latency. Two large-effect mutations in the genes Proc-R and Piezo were further confirmed via genetic rescue. This work greatly enhances our understanding of the genetic factors that influence variation in sleep latency.

Ribogenesis boosts controlled by HEATR1-MYC interplay promote transition into brain tumour growth

Cell commitment to tumourigenesis and the onset of uncontrolled growth are critical determinants in cancer development but the early events directing tumour initiating cell (TIC) fate remain unclear. We reveal a single-cell transcriptome profile of brain TICs transitioning into tumour growth using the brain tumour (brat) neural stem cell-based Drosophila model. Prominent changes in metabolic and proteostasis-associated processes including ribogenesis are identified. Increased ribogenesis is a known cell adaptation in established tumours. Here we propose that brain TICs boost ribogenesis prior to tumour growth. In brat-deficient TICs, we show that this dramatic change is mediated by upregulated HEAT-Repeat Containing 1 (HEATR1) to promote ribosomal RNA generation, TIC enlargement and onset of overgrowth. High HEATR1 expression correlates with poor glioma patient survival and patient-derived glioblastoma stem cells rely on HEATR1 for enhanced ribogenesis and tumourigenic potential. Finally, we show that HEATR1 binds the master growth regulator MYC, promotes its nucleolar localisation and appears required for MYC-driven ribogenesis, suggesting a mechanism co-opted in ribogenesis reprogramming during early brain TIC development.

Inferring Interaction Networks from Transcriptomic Data: Methods and Applications

Transcriptomic data is a treasure trove in modern molecular biology, as it offers a comprehensive viewpoint into the intricate nuances of gene expression dynamics underlying biological systems. This genetic information must be utilized to infer biomolecular interaction networks that can provide insights into the complex regulatory mechanisms underpinning the dynamic cellular processes. Gene regulatory networks and protein-protein interaction networks are two major classes of such networks. This chapter thoroughly investigates the wide range of methodologies used for distilling insightful revelations from transcriptomic data that include association-based methods (based on correlation among expression vectors), probabilistic models (using Bayesian and Gaussian models), and interologous methods. We reviewed different approaches for evaluating the significance of interactions based on the network topology and biological functions of the interacting molecules and discuss various strategies for the identification of functional modules. The chapter concludes with highlighting network-based techniques of prioritizing key genes, outlining the centrality-based, diffusion- based, and subgraph-based methods. The chapter provides a meticulous framework for investigating transcriptomic data to uncover assembly of complex molecular networks for their adaptable analyses across a broad spectrum of biological domains.

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.

The effects of methylphenidate and atomoxetine on Drosophila brain at single-cell resolution and potential drug repurposing for ADHD treatment

The stimulant methylphenidate (MPH) and the non-stimulant atomoxetine (ATX) are frequently used for the treatment of attention-deficit/hyperactivity disorder (ADHD); however, the function of these drugs in different types of brain cells and their effects on related genes remain largely unknown. To address these questions, we built a pipeline for the simultaneous examination of the activity behavior and transcriptional responses of Drosophila melanogaster at single-cell resolution following drug treatment. We selected the Drosophila with significantly increased locomotor activities (hyperactivity-like behavior) following the administration of each drug in comparison with the control (same food as the drug-treated groups with 5% sucrose, yeast, and blue food dye solution) using EasyFlyTracker. Subsequently, single cell RNA sequencing (scRNASEQ) was used to capture the transcriptome of 82,917 cells, unsupervised clustering analysis of which yielded 28 primary cell clusters representing the major cell types in adult Drosophila brain. Indeed, both neuronal and glial cells responded to MPH and ATX. Further analysis of differentially expressed genes (DEGs) revealed distinct transcriptional changes associated with these two drugs, such as two well-studied dopamine receptor genes (Dop2R and DopEcR) were responsive to MPH but not to ATX at their optimal doses, in addition to genes involved in dopamine metabolism pathways such as Syt1, Sytalpha, Syt7, and Ih in different cell types. More importantly, MPH also suppressed the expression of genes encoding other neurotransmitter receptors and synaptic signaling molecules in many cell types, especially those for Glu and GABA, while the responsive effects of ATX were much weaker. In addition to monoaminergic neuronal transmitters, other neurotransmitters have also shown a similar pattern with respect to a stronger effect associated with MPH than with ATX. Moreover, we identified four distinct glial cell subtypes responsive to the two drugs and detected a greater number of differentially expressed genes associated with ensheathing and astrocyte-like glia. Furthermore, our study provides a rich resource of candidate target genes, supported by drug set enrichment analysis (P = 2.10E-4; hypergeometric test), for the further exploration of drug repurposing. The whole list of candidates can be found at ADHDrug ( http://adhdrug.cibr.ac.cn/ ). In conclusion, we propose a fast and cost-efficient pipeline to explore the underlying molecular mechanisms of ADHD drug treatment in Drosophila brain at single-cell resolution, which may further facilitate drug repurposing applications.

Drosophila as a Model to Understand Second Heart Field Development

The genetic model system Drosophila has contributed fundamentally to our understanding of mammalian heart specification, development, and congenital heart disease. The relatively simple Drosophila heart is a linear muscular tube that is specified and develops in the embryo and persists throughout the life of the animal. It functions at all stages to circulate hemolymph within the open circulatory system of the body. During Drosophila metamorphosis, the cardiac tube is remodeled, and a new layer of muscle fibers spreads over the ventral surface of the heart to form the ventral longitudinal muscles. The formation of these fibers depends critically upon genes known to be necessary for mammalian second heart field (SHF) formation. Here, we review the prior contributions of the Drosophila system to the understanding of heart development and disease, discuss the importance of the SHF to mammalian heart development and disease, and then discuss how the ventral longitudinal adult cardiac muscles can serve as a novel model for understanding SHF development and disease.

Nicotinic acetylcholine receptor signaling maintains epithelial barrier integrity

Disruption of epithelial barriers is a common disease manifestation in chronic degenerative diseases of the airways, lung, and intestine. Extensive human genetic studies have identified risk loci in such diseases, including in chronic obstructive pulmonary disease (COPD) and inflammatory bowel diseases. The genes associated with these loci have not fully been determined, and functional characterization of such genes requires extensive studies in model organisms. Here, we report the results of a screen in Drosophila melanogaster that allowed for rapid identification, validation, and prioritization of COPD risk genes that were selected based on risk loci identified in human genome-wide association studies (GWAS). Using intestinal barrier dysfunction in flies as a readout, our results validate the impact of candidate gene perturbations on epithelial barrier function in 56% of the cases, resulting in a prioritized target gene list. We further report the functional characterization in flies of one family of these genes, encoding for nicotinic acetylcholine receptor (nAchR) subunits. We find that nAchR signaling in enterocytes of the fly gut promotes epithelial barrier function and epithelial homeostasis by regulating the production of the peritrophic matrix. Our findings identify COPD-associated genes critical for epithelial barrier maintenance, and provide insight into the role of epithelial nAchR signaling for homeostasis.

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.

Conserved autism-associated genes tune social feeding behavior in C. elegans

Animal foraging is an essential and evolutionarily conserved behavior that occurs in social and solitary contexts, but the underlying molecular pathways are not well defined. We discover that conserved autism-associated genes (NRXN1(nrx-1), NLGN3(nlg-1), GRIA1,2,3(glr-1), GRIA2(glr-2), and GLRA2,GABRA3(avr-15)) regulate aggregate feeding in C. elegans, a simple social behavior. NRX-1 functions in chemosensory neurons (ADL and ASH) independently of its postsynaptic partner NLG-1 to regulate social feeding. Glutamate from these neurons is also crucial for aggregate feeding, acting independently of NRX-1 and NLG-1. Compared to solitary counterparts, social animals show faster presynaptic release and more presynaptic release sites in ASH neurons, with only the latter requiring nrx-1. Disruption of these distinct signaling components additively converts behavior from social to solitary. Aggregation induced by circuit activation is also dependent on nrx-1. Collectively, we find that aggregate feeding is tuned by conserved autism-associated genes through complementary synaptic mechanisms, revealing molecular principles driving social feeding.

Molecular traces of Drosophila hemocytes reveal transcriptomic conservation with vertebrate myeloid cells

Drosophila hemocytes serve as the primary defense system against harmful threats, allowing the animals to thrive. Hemocytes are often compared to vertebrate innate immune system cells due to the observed functional similarities between the two. However, the similarities have primarily been established based on a limited number of genes and their functional homologies. Thus, a systematic analysis using transcriptomic data could offer novel insights into Drosophila hemocyte function and provide new perspectives on the evolution of the immune system. Here, we performed cross-species comparative analyses using single-cell RNA sequencing data from Drosophila and vertebrate immune cells. We found several conserved markers for the cluster of differentiation (CD) genes in Drosophila hemocytes and validated the role of CG8501 (CD59) in phagocytosis by plasmatocytes, which function much like macrophages in vertebrates. By comparing whole transcriptome profiles in both supervised and unsupervised analyses, we showed that Drosophila hemocytes are largely homologous to vertebrate myeloid cells, especially plasmatocytes to monocytes/macrophages and prohemocyte 1 (PH1) to hematopoietic stem cells. Furthermore, a small subset of prohemocytes with hematopoietic potential displayed homology with hematopoietic progenitor populations in vertebrates. Overall, our results provide a deeper understanding of molecular conservation in the Drosophila immune system.

Transcriptome analysis reveals temporally regulated genetic networks during Drosophila border cell collective migration

BACKGROUND

Collective cell migration underlies many essential processes, including sculpting organs during embryogenesis, wound healing in the adult, and metastasis of cancer cells. At mid-oogenesis, Drosophila border cells undergo collective migration. Border cells round up into a small group at the pre-migration stage, detach from the epithelium and undergo a dynamic and highly regulated migration at the mid-migration stage, and stop at the oocyte, their final destination, at the post-migration stage. While specific genes that promote cell signaling, polarization of the cluster, formation of protrusions, and cell-cell adhesion are known to regulate border cell migration, there may be additional genes that promote these distinct active phases of border cell migration. Therefore, we sought to identify genes whose expression patterns changed during border cell migration.

RESULTS

We performed RNA-sequencing on border cells isolated at pre-, mid-, and post-migration stages. We report that 1,729 transcripts, in nine co-expression gene clusters, are temporally and differentially expressed across the three migration stages. Gene ontology analyses and constructed protein-protein interaction networks identified genes expected to function in collective migration, such as regulators of the cytoskeleton, adhesion, and tissue morphogenesis, but also uncovered a notable enrichment of genes involved in immune signaling, ribosome biogenesis, and stress responses. Finally, we validated the in vivo expression and function of a subset of identified genes in border cells.

CONCLUSIONS

Overall, our results identified differentially and temporally expressed genetic networks that may facilitate the efficient development and migration of border cells. The genes identified here represent a wealth of new candidates to investigate the molecular nature of dynamic collective cell migrations in developing tissues.

The alliance of genome resources: transforming comparative genomics

Comparing genomic and biological characteristics across multiple species is essential to using model systems to investigate the molecular and cellular mechanisms underlying human biology and disease and to translate mechanistic insights from studies in model organisms for clinical applications. Building a scalable knowledge commons platform that supports cross-species comparison of rich, expertly curated knowledge regarding gene function, phenotype, and disease associations available for model organisms and humans is the primary mission of the Alliance of Genome Resources (the Alliance). The Alliance is a consortium of seven model organism knowledgebases (mouse, rat, yeast, nematode, zebrafish, frog, fruit fly) and the Gene Ontology resource. The Alliance uses a common set of gene ortholog assertions as the basis for comparing biological annotations across the organisms represented in the Alliance. The major types of knowledge associated with genes that are represented in the Alliance database currently include gene function, phenotypic alleles and variants, human disease associations, pathways, gene expression, and both protein-protein and genetic interactions. The Alliance has enhanced the ability of researchers to easily compare biological annotations for common data types across model organisms and human through the implementation of shared programmatic access mechanisms, data-specific web pages with a unified "look and feel", and interactive user interfaces specifically designed to support comparative biology. The modular infrastructure developed by the Alliance allows the resource to serve as an extensible "knowledge commons" capable of expanding to accommodate additional model organisms.

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.

A conserved role for frizzled in sleep architecture

Previous studies of natural variants in Drosophila melanogaster implicated the Wnt signaling receptor frizzled in sleep. Given that the Wnt signaling pathway is highly conserved across species, we hypothesized that frizzled class receptor 1 (Fzd1), the murine homolog of frizzled, would also have a role in sleep. Using a CRISPR transgenic approach, we removed most of the Fzd1 coding region from C57BL/6N mice. We used a video assay to measure sleep characteristics in Fzd1-deficient mice. As Wnt signaling is known to affect visuospatial memory, we also examined the impact of the deletion on learning and memory using the novel object recognition (NOR) paradigm. Fzd1-deficient mice had altered sleep compared to littermate controls. The mice did not respond differently to the NOR paradigm compared to controls but did display anxiety-like behavior. Our strategy demonstrates that the study of natural variation in Drosophila sleep translates into candidate genes for sleep in vertebrate species such as the mouse.

Updates to the Alliance of Genome Resources Central Infrastructure Alliance of Genome Resources Consortium

The Alliance of Genome Resources (Alliance) is an extensible coalition of knowledgebases focused on the genetics and genomics of intensively-studied model organisms. The Alliance is organized as individual knowledge centers with strong connections to their research communities and a centralized software infrastructure, discussed here. Model organisms currently represented in the Alliance are budding yeast, C. elegans, Drosophila, zebrafish, frog, laboratory mouse, laboratory rat, and the Gene Ontology Consortium. The project is in a rapid development phase to harmonize knowledge, store it, analyze it, and present it to the community through a web portal, direct downloads, and APIs. Here we focus on developments over the last two years. Specifically, we added and enhanced tools for browsing the genome (JBrowse), downloading sequences, mining complex data (AllianceMine), visualizing pathways, full-text searching of the literature (Textpresso), and sequence similarity searching (SequenceServer). We enhanced existing interactive data tables and added an interactive table of paralogs to complement our representation of orthology. To support individual model organism communities, we implemented species-specific "landing pages" and will add disease-specific portals soon; in addition, we support a common community forum implemented in Discourse. We describe our progress towards a central persistent database to support curation, the data modeling that underpins harmonization, and progress towards a state-of-the art literature curation system with integrated Artificial Intelligence and Machine Learning (AI/ML).

Cross species systems biology discovers glial DDR2, STOM, and KANK2 as therapeutic targets in progressive supranuclear palsy

Progressive supranuclear palsy (PSP) is a neurodegenerative parkinsonian disorder characterized by cell-type-specific tau lesions in neurons and glia. Prior work uncovered transcriptome changes in human PSP brains, although their cell-specificity is unknown. Further, systematic data integration and experimental validation platforms to prioritize brain transcriptional perturbations as therapeutic targets in PSP are currently lacking. In this study, we combine bulk tissue (n = 408) and single nucleus RNAseq (n = 34) data from PSP and control brains with transcriptome data from a mouse tauopathy and experimental validations in Drosophila tau models for systematic discovery of high-confidence expression changes in PSP with therapeutic potential. We discover, replicate, and annotate thousands of differentially expressed genes in PSP, many of which reside in glia-enriched co-expression modules and cells. We prioritize DDR2, STOM, and KANK2 as promising therapeutic targets in PSP with striking cross-species validations. We share our findings and data via our interactive application tool PSP RNAseq Atlas ( https://rtools.mayo.edu/PSP_RNAseq_Atlas/ ). Our findings reveal robust glial transcriptome changes in PSP, provide a cross-species systems biology approach, and a tool for therapeutic target discoveries in PSP with potential application in other neurodegenerative diseases.

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.

angptl4 gene expression as a marker of adaptive homeostatic response to social isolation across the lifespan in zebrafish

Social isolation has detrimental health effects, but the underlying mechanisms are unclear. Here, we investigated the impact of 2 weeks of isolation on behavior and gene expression in the central nervous system at different life stages of zebrafish. Results showed that socially deprived young adult zebrafish experienced increased anxiety, accompanied by changes in gene expression. Most gene expression patterns returned to normal within 24 hours of reintroduction to a social environment, except angptl4, which was upregulated after reintroduction, suggesting an adaptive mechanism. Similarly, aging zebrafish displayed heightened anxiety and increased central nervous system expression of angptl4 during isolation, but effects were reversed upon reintroduction to a social group. The findings imply that angptl4 plays a homeostatic role in response to social isolation, which varies across the lifespan. The study emphasizes the importance of social interactions for psychological well-being and highlights the negative consequences of isolation, especially in older individuals. Further research may unravel how social isolation affects angptl4 expression and its developmental and aging effects.

Insulin-mediated endothelin signaling is antiviral during West Nile virus infection

IMPORTANCE

Arboviruses, particularly those transmitted by mosquitoes, pose a significant threat to humans and are an increasing concern because of climate change, human activity, and expanding vector-competent populations. West Nile virus is of significant concern as the most frequent mosquito-borne disease transmitted annually within the continental United States. Here, we identify a previously uncharacterized signaling pathway that impacts West Nile virus infection, namely endothelin signaling. Additionally, we demonstrate that we can successfully translate results obtained from D. melanogaster into the more relevant human system. Our results add to the growing field of insulin-mediated antiviral immunity and identify potential biomarkers or intervention targets to better address West Nile virus infection and severe disease.

Genome dilution by cell growth drives starvation-like proteome remodeling in mammalian and yeast cells

Cell size is tightly controlled in healthy tissues and single-celled organisms, but it remains unclear how size influences cell physiology. Increasing cell size was recently shown to remodel the proteomes of cultured human cells, demonstrating that large and small cells of the same type can be biochemically different. Here, we corroborate these results in mouse hepatocytes and extend our analysis using yeast. We find that size-dependent proteome changes are highly conserved and mostly independent of metabolic state. As eukaryotic cells grow larger, the dilution of the genome elicits a starvation-like proteome phenotype, suggesting that growth in large cells is limited by the genome in a manner analogous to a limiting nutrient. We also demonstrate that the proteomes of replicatively-aged yeast are primarily determined by their large size. Overall, our data suggest that genome concentration is a universal determinant of proteome content in growing cells.

Drosophila Mpv17 forms an ion channel and regulates energy metabolism

Mutations in MPV17 are a major contributor to mitochondrial DNA (mtDNA) depletion syndromes, a group of inherited genetic conditions due to mtDNA instability. To investigate the role of MPV17 in mtDNA maintenance, we generated and characterized a Drosophila melanogaster Mpv17 (dMpv17) KO model showing that the absence of dMpv17 caused profound mtDNA depletion in the fat body but not in other tissues, increased glycolytic flux and reduced lifespan in starvation. Accordingly, the expression of key genes of glycogenolysis and glycolysis was upregulated in dMpv17 KO flies. In addition, we demonstrated that dMpv17 formed a channel in planar lipid bilayers at physiological ionic conditions, and its electrophysiological hallmarks were affected by pathological mutations. Importantly, the reconstituted channel translocated uridine but not orotate across the membrane. Our results indicate that dMpv17 forms a channel involved in translocation of key metabolites and highlight the importance of dMpv17 in energy homeostasis and mitochondrial function.

Drosophila CASK regulates brain size and neuronal morphogenesis, providing a genetic model of postnatal microcephaly suitable for drug discovery

BACKGROUND

CASK-related neurodevelopmental disorders are untreatable. Affected children show variable severity, with microcephaly, intellectual disability (ID), and short stature as common features. X-linked human CASK shows dosage sensitivity with haploinsufficiency in females. CASK protein has multiple domains, binding partners, and proposed functions at synapses and in the nucleus. Human and Drosophila CASK show high amino-acid-sequence similarity in all functional domains. Flies homozygous for a hypomorphic CASK mutation (∆18) have motor and cognitive deficits. A Drosophila genetic model of CASK-related disorders could have great scientific and translational value.

METHODS

We assessed the effects of CASK loss of function on morphological phenotypes in Drosophila using established genetic, histological, and primary neuronal culture approaches. NeuronMetrics software was used to quantify neurite-arbor morphology. Standard nonparametric statistics methods were supplemented by linear mixed effects modeling in some cases. Microfluidic devices of varied dimensions were fabricated and numerous fluid-flow parameters were used to induce oscillatory stress fields on CNS tissue. Dissociation into viable neurons and neurite outgrowth in vitro were assessed.

RESULTS

We demonstrated that ∆18 homozygous flies have small brains, small heads, and short bodies. When neurons from developing CASK-mutant CNS were cultured in vitro, they grew small neurite arbors with a distinctive, quantifiable "bushy" morphology that was significantly rescued by transgenic CASK+. As in humans, the bushy phenotype showed dosage-sensitive severity. To overcome the limitations of manual tissue trituration for neuronal culture, we optimized the design and operation of a microfluidic system for standardized, automated dissociation of CNS tissue into individual viable neurons. Neurons from CASK-mutant CNS dissociated in the microfluidic system recapitulate the bushy morphology. Moreover, for any given genotype, device-dissociated neurons grew larger arbors than did manually dissociated neurons. This automated dissociation method is also effective for rodent CNS.

CONCLUSIONS

These biological and engineering advances set the stage for drug discovery using the Drosophila model of CASK-related disorders. The bushy phenotype provides a cell-based assay for compound screening. Nearly a dozen genes encoding CASK-binding proteins or transcriptional targets also have brain-development mutant phenotypes, including ID. Hence, drugs that improve CASK phenotypes might also benefit children with disorders due to mutant CASK partners.

Literature-based predictions of Mendelian disease therapies

In the effort to treat Mendelian disorders, correcting the underlying molecular imbalance may be more effective than symptomatic treatment. Identifying treatments that might accomplish this goal requires extensive and up-to-date knowledge of molecular pathways-including drug-gene and gene-gene relationships. To address this challenge, we present "parsing modifiers via article annotations" (PARMESAN), a computational tool that searches PubMed and PubMed Central for information to assemble these relationships into a central knowledge base. PARMESAN then predicts putatively novel drug-gene relationships, assigning an evidence-based score to each prediction. We compare PARMESAN's drug-gene predictions to all of the drug-gene relationships displayed by the Drug-Gene Interaction Database (DGIdb) and show that higher-scoring relationship predictions are more likely to match the directionality (up- versus down-regulation) indicated by this database. PARMESAN had more than 200,000 drug predictions scoring above 8 (as one example cutoff), for more than 3,700 genes. Among these predicted relationships, 210 were registered in DGIdb and 201 (96%) had matching directionality. This publicly available tool provides an automated way to prioritize drug screens to target the most-promising drugs to test, thereby saving time and resources in the development of therapeutics for genetic disorders.

PINK1, Keap1, and Rtnl1 regulate selective clearance of endoplasmic reticulum during development

Selective clearance of organelles, including endoplasmic reticulum (ER) and mitochondria, by autophagy plays an important role in cell health. Here, we describe a developmentally programmed selective ER clearance by autophagy. We show that Parkinson's disease-associated PINK1, as well as Atl, Rtnl1, and Trp1 receptors, regulate ER clearance by autophagy. The E3 ubiquitin ligase Parkin functions downstream of PINK1 and is required for mitochondrial clearance while having the opposite function in ER clearance. By contrast, Keap1 and the E3 ubiquitin ligase Cullin3 function downstream of PINK1 to regulate ER clearance by influencing Rtnl1 and Atl. PINK1 regulates a change in Keap1 localization and Keap1-dependent ubiquitylation of the ER-phagy receptor Rtnl1 to facilitate ER clearance. Thus, PINK1 regulates the selective clearance of ER and mitochondria by influencing the balance of Keap1- and Parkin-dependent ubiquitylation of substrates that determine which organelle is removed by autophagy.

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.

Using Drosophila to identify naturally occurring genetic modifiers of amyloid beta 42- and tau-induced toxicity

Alzheimer's disease is characterized by 2 pathological proteins, amyloid beta 42 and tau. The majority of Alzheimer's disease cases in the population are sporadic and late-onset Alzheimer's disease, which exhibits high levels of heritability. While several genetic risk factors for late-onset Alzheimer's disease have been identified and replicated in independent studies, including the ApoE ε4 allele, the great majority of the heritability of late-onset Alzheimer's disease remains unexplained, likely due to the aggregate effects of a very large number of genes with small effect size, as well as to biases in sample collection and statistical approaches. Here, we present an unbiased forward genetic screen in Drosophila looking for naturally occurring modifiers of amyloid beta 42- and tau-induced ommatidial degeneration. Our results identify 14 significant SNPs, which map to 12 potential genes in 8 unique genomic regions. Our hits that are significant after genome-wide correction identify genes involved in neuronal development, signal transduction, and organismal development. Looking more broadly at suggestive hits (P < 10-5), we see significant enrichment in genes associated with neurogenesis, development, and growth as well as significant enrichment in genes whose orthologs have been identified as significantly or suggestively associated with Alzheimer's disease in human GWAS studies. These latter genes include ones whose orthologs are in close proximity to regions in the human genome that are associated with Alzheimer's disease, but where a causal gene has not been identified. Together, our results illustrate the potential for complementary and convergent evidence provided through multitrait GWAS in Drosophila to supplement and inform human studies, helping to identify the remaining heritability and novel modifiers of complex diseases.

A comprehensive Drosophila resource to identify key functional interactions between SARS-CoV-2 factors and host proteins

Development of effective therapies against SARS-CoV-2 infections relies on mechanistic knowledge of virus-host interface. Abundant physical interactions between viral and host proteins have been identified, but few have been functionally characterized. Harnessing the power of fly genetics, we develop a comprehensive Drosophila COVID-19 resource (DCR) consisting of publicly available strains for conditional tissue-specific expression of all SARS-CoV-2 encoded proteins, UAS-human cDNA transgenic lines encoding established host-viral interacting factors, and GAL4 insertion lines disrupting fly homologs of SARS-CoV-2 human interacting proteins. We demonstrate the utility of the DCR to functionally assess SARS-CoV-2 genes and candidate human binding partners. We show that NSP8 engages in strong genetic interactions with several human candidates, most prominently with the ATE1 arginyltransferase to induce actin arginylation and cytoskeletal disorganization, and that two ATE1 inhibitors can reverse NSP8 phenotypes. The DCR enables parallel global-scale functional analysis of SARS-CoV-2 components in a prime genetic model system.

REPTOR and CREBRF encode key regulators of muscle energy metabolism

Metabolic flexibility of muscle tissue describes the adaptive capacity to use different energy substrates according to their availability. The disruption of this ability associates with metabolic disease. Here, using a Drosophila model of systemic metabolic dysfunction triggered by yorkie-induced gut tumors, we show that the transcription factor REPTOR is an important regulator of energy metabolism in muscles. We present evidence that REPTOR is activated in muscles of adult flies with gut yorkie-tumors, where it modulates glucose metabolism. Further, in vivo studies indicate that sustained activity of REPTOR is sufficient in wildtype muscles to repress glycolysis and increase tricarboxylic acid (TCA) cycle metabolites. Consistent with the fly studies, higher levels of CREBRF, the mammalian ortholog of REPTOR, reduce glycolysis in mouse myotubes while promoting oxidative metabolism. Altogether, our results define a conserved function for REPTOR and CREBRF as key regulators of muscle energy metabolism.

Single-cell profiling of the developing embryonic heart in Drosophila

Drosophila is an important model for studying heart development and disease. Yet, single-cell transcriptomic data of its developing heart have not been performed. Here, we report single-cell profiling of the entire fly heart using ∼3000 Hand-GFP embryos collected at five consecutive developmental stages, ranging from bilateral migrating rows of cardiac progenitors to a fused heart tube. The data revealed six distinct cardiac cell types in the embryonic fly heart: cardioblasts, both Svp+ and Tin+ subtypes; and five types of pericardial cell (PC) that can be distinguished by four key transcription factors (Eve, Odd, Ct and Tin) and include the newly described end of the line PC. Notably, the embryonic fly heart combines transcriptional signatures of the mammalian first and second heart fields. Using unique markers for each heart cell type, we defined their number and location during heart development to build a comprehensive 3D cell map. These data provide a resource to track the expression of any gene in the developing fly heart, which can serve as a reference to study genetic perturbations and cardiac diseases.

Genetic and genomic analyses of Drosophila melanogaster models of chromatin modification disorders

Switch/sucrose nonfermentable (SWI/SNF)-related intellectual disability disorders (SSRIDDs) and Cornelia de Lange syndrome are rare syndromic neurodevelopmental disorders with overlapping clinical phenotypes. SSRIDDs are associated with the BAF (Brahma-Related Gene-1 associated factor) complex, whereas CdLS is a disorder of chromatin modification associated with the cohesin complex. Here, we used RNA interference in Drosophila melanogaster to reduce the expression of six genes (brm, osa, Snr1, SMC1, SMC3, vtd) orthologous to human genes associated with SSRIDDs and CdLS. These fly models exhibit changes in sleep, activity, startle behavior (a proxy for sensorimotor integration), and brain morphology. Whole genome RNA sequencing identified 9,657 differentially expressed genes (FDR < 0.05), 156 of which are differentially expressed in both sexes in SSRIDD- and CdLS-specific analyses, including Bap60, which is orthologous to SMARCD1, an SSRIDD-associated BAF component. k-means clustering reveals genes co-regulated within and across SSRIDD and CdLS fly models. RNAi-mediated reduction of expression of six genes co-regulated with focal genes brm, osa, and/or Snr1 recapitulated changes in the behavior of the focal genes. Based on the assumption that fundamental biological processes are evolutionarily conserved, Drosophila models can be used to understand underlying molecular effects of variants in chromatin-modification pathways and may aid in the discovery of drugs that ameliorate deleterious phenotypic effects.

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.

A new metabolic model of Drosophila melanogaster and the integrative analysis of Parkinson's disease

High conservation of the disease-associated genes between flies and humans facilitates the common use of Drosophila melanogaster to study metabolic disorders under controlled laboratory conditions. However, metabolic modeling studies are highly limited for this organism. We here report a comprehensively curated genome-scale metabolic network model of Drosophila using an orthology-based approach. The gene coverage and metabolic information of the draft model derived from a reference human model were expanded via Drosophila-specific KEGG and MetaCyc databases, with several curation steps to avoid metabolic redundancy and stoichiometric inconsistency. Furthermore, we performed literature-based curations to improve gene-reaction associations, subcellular metabolite locations, and various metabolic pathways. The performance of the resulting Drosophila model (8,230 reactions, 6,990 metabolites, and 2,388 genes), iDrosophila1 (https://github.com/SysBioGTU/iDrosophila), was assessed using flux balance analysis in comparison with the other currently available fly models leading to superior or comparable results. We also evaluated the transcriptome-based prediction capacity of iDrosophila1, where differential metabolic pathways during Parkinson's disease could be successfully elucidated. Overall, iDrosophila1 is promising to investigate system-level metabolic alterations in response to genetic and environmental perturbations.

Harmonized cross-species cell atlases of trigeminal and dorsal root ganglia

Peripheral sensory neurons in the dorsal root ganglion (DRG) and trigeminal ganglion (TG) are specialized to detect and transduce diverse environmental stimuli including touch, temperature, and pain to the central nervous system. Recent advances in single-cell RNA-sequencing (scRNA-seq) have provided new insights into the diversity of sensory ganglia cell types in rodents, non-human primates, and humans, but it remains difficult to compare transcriptomically defined cell types across studies and species. Here, we built cross-species harmonized atlases of DRG and TG cell types that describe 18 neuronal and 11 non-neuronal cell types across 6 species and 19 studies. We then demonstrate the utility of this harmonized reference atlas by using it to annotate newly profiled DRG nuclei/cells from both human and the highly regenerative axolotl. We observe that the transcriptomic profiles of sensory neuron subtypes are broadly similar across vertebrates, but the expression of functionally important neuropeptides and channels can vary notably. The new resources and data presented here can guide future studies in comparative transcriptomics, simplify cell type nomenclature differences across studies, and help prioritize targets for future pain therapy development.

Dynamically regulated transcription factors are encoded by highly unstable mRNAs in the Drosophila larval brain

The level of each RNA species depends on the balance between its rates of production and decay. Although previous studies have measured RNA decay across the genome in tissue culture and single-celled organisms, few experiments have been performed in intact complex tissues and organs. It is therefore unclear whether the determinants of RNA decay found in cultured cells are preserved in an intact tissue, and whether they differ between neighboring cell types and are regulated during development. To address these questions, we measured RNA synthesis and decay rates genome wide via metabolic labeling of whole cultured Drosophila larval brains using 4-thiouridine. Our analysis revealed that decay rates span a range of more than 100-fold, and that RNA stability is linked to gene function, with mRNAs encoding transcription factors being much less stable than mRNAs involved in core metabolic functions. Surprisingly, among transcription factor mRNAs there was a clear demarcation between more widely used transcription factors and those that are expressed only transiently during development. mRNAs encoding transient transcription factors are among the least stable in the brain. These mRNAs are characterized by epigenetic silencing in most cell types, as shown by their enrichment with the histone modification H3K27me3. Our data suggest the presence of an mRNA destabilizing mechanism targeted to these transiently expressed transcription factors to allow their levels to be regulated rapidly with high precision. Our study also demonstrates a general method for measuring mRNA transcription and decay rates in intact organs or tissues, offering insights into the role of mRNA stability in the regulation of complex developmental programs.

Yeast NDI1 reconfigures neuronal metabolism and prevents the unfolded protein response in mitochondrial complex I deficiency

Mutations in subunits of the mitochondrial NADH dehydrogenase cause mitochondrial complex I deficiency, a group of severe neurological diseases that can result in death in infancy. The pathogenesis of complex I deficiency remain poorly understood, and as a result there are currently no available treatments. To better understand the underlying mechanisms, we modelled complex I deficiency in Drosophila using knockdown of the mitochondrial complex I subunit ND-75 (NDUFS1) specifically in neurons. Neuronal complex I deficiency causes locomotor defects, seizures and reduced lifespan. At the cellular level, complex I deficiency does not affect ATP levels but leads to mitochondrial morphology defects, reduced endoplasmic reticulum-mitochondria contacts and activation of the endoplasmic reticulum unfolded protein response (UPR) in neurons. Multi-omic analysis shows that complex I deficiency dramatically perturbs mitochondrial metabolism in the brain. We find that expression of the yeast non-proton translocating NADH dehydrogenase NDI1, which reinstates mitochondrial NADH oxidation but not ATP production, restores levels of several key metabolites in the brain in complex I deficiency. Remarkably, NDI1 expression also reinstates endoplasmic reticulum-mitochondria contacts, prevents UPR activation and rescues the behavioural and lifespan phenotypes caused by complex I deficiency. Together, these data show that metabolic disruption due to loss of neuronal NADH dehydrogenase activity cause UPR activation and drive pathogenesis in complex I deficiency.

Cross-Species Transcriptomics Analysis Highlights Conserved Molecular Responses to Per- and Polyfluoroalkyl Substances

In recent decades, per- and polyfluoroalkyl substances (PFASs) have garnered widespread public attention due to their persistence in the environment and detrimental effects on the health of living organisms, spurring the generation of several transcriptome-centered investigations to understand the biological basis of their mechanism. In this study, we collected 2144 publicly available samples from seven distinct animal species to examine the molecular responses to PFAS exposure and to determine if there are conserved responses. Our comparative transcriptional analysis revealed that exposure to PFAS is conserved across different tissues, molecules and species. We identified and reported several genes exhibiting consistent and evolutionarily conserved transcriptional response to PFASs, such as ESR1, HADHA and ID1, as well as several pathways including lipid metabolism, immune response and hormone pathways. This study provides the first evidence that distinct PFAS molecules induce comparable transcriptional changes and affect the same metabolic processes across inter-species borders. Our findings have significant implications for understanding the impact of PFAS exposure on living organisms and the environment. We believe that this study offers a novel perspective on the molecular responses to PFAS exposure and provides a foundation for future research into developing strategies for mitigating the detrimental effects of these substances in the ecosystem.

The Drosophila gene encoding JIG protein (CG14850) is critical for CrebA nuclear trafficking during development

Coordination of mitochondrial and nuclear processes is key to the cellular health; however, very little is known about the molecular mechanisms regulating nuclear-mitochondrial crosstalk. Here, we report a novel molecular mechanism controlling the shuttling of CREB (cAMP response element-binding protein) protein complex between mitochondria and nucleoplasm. We show that a previously unknown protein, herein termed as Jig, functions as a tissue-specific and developmental timing-specific coregulator in the CREB pathway. Our results demonstrate that Jig shuttles between mitochondria and nucleoplasm, interacts with CrebA protein and controls its delivery to the nucleus, thus triggering CREB-dependent transcription in nuclear chromatin and mitochondria. Ablating the expression of Jig prevents CrebA from localizing to the nucleoplasm, affecting mitochondrial functioning and morphology and leads to Drosophila developmental arrest at the early third instar larval stage. Together, these results implicate Jig as an essential mediator of nuclear and mitochondrial processes. We also found that Jig belongs to a family of nine similar proteins, each of which has its own tissue- and time-specific expression profile. Thus, our results are the first to describe the molecular mechanism regulating nuclear and mitochondrial processes in a tissue- and time-specific manner.

Expression of retrotransposons contributes to aging in Drosophila

Retrotransposons are a class of transposable elements capable of self-replication and insertion into new genomic locations. Across species, the mobilization of retrotransposons in somatic cells has been suggested to contribute to the cell and tissue functional decline that occurs during aging. Retrotransposons are broadly expressed across cell types, and de novo insertions have been observed to correlate with tumorigenesis. However, the extent to which new retrotransposon insertions occur during normal aging and their effect on cellular and animal function remains understudied. Here, we use a single nucleus whole genome sequencing approach in Drosophila to directly test whether transposon insertions increase with age in somatic cells. Analyses of nuclei from thoraces and indirect flight muscles using a newly developed pipeline, Retrofind, revealed no significant increase in the number of transposon insertions with age. Despite this, reducing the expression of two different retrotransposons, 412 and Roo, extended lifespan, but did not alter indicators of health such as stress resistance. This suggests a key role for transposon expression and not insertion in regulating longevity. Transcriptomic analyses revealed similar changes to gene expression in 412 and Roo knockdown flies and highlighted changes to genes involved in proteolysis and immune function as potential contributors to the observed changes in longevity. Combined, our data show a clear link between retrotransposon expression and aging.

A novel neurodevelopmental syndrome caused by loss-of-function of the Zinc Finger Homeobox 3 (ZFHX3) gene

Neurodevelopmental disorders (NDDs) result from impaired development and functioning of the brain. Here, we identify loss-of-function variation in ZFHX3 as a novel cause for syndromic intellectual disability (ID). ZFHX3, previously known as ATBF1, is a zinc-finger homeodomain transcription factor involved in multiple biological processes including cell differentiation and tumorigenesis. Through international collaboration, we collected clinical and morphometric data (Face2Gene) of 41 individuals with protein truncating variants (PTVs) or (partial) deletions of ZFHX3 . We used data mining, RNA and protein analysis to identify the subcellular localization and spatiotemporal expression of ZFHX3 in multiple in vitro models. We identified the DNA targets of ZFHX3 using ChIP seq. Immunoprecipitation followed by mass spectrometry indicated potential binding partners of endogenous ZFHX3 in neural stem cells that were subsequently confirmed by reversed co-immunoprecipitation and western blot. We evaluated a DNA methylation profile associated with ZFHX3 haploinsufficiency using DNA methylation analysis on whole blood extracted DNA of six individuals with ZFHX3 PTVs and four with a (partial) deletion of ZFHX3 . A reversed genetic approach characterized the ZFHX3 orthologue in Drosophila melanogaster . Loss-of-function variation of ZFHX3 consistently associates with (mild) ID and/or behavioural problems, postnatal growth retardation, feeding difficulties, and recognizable facial characteristics, including the rare occurrence of cleft palate. Nuclear abundance of ZFHX3 increases during human brain development and neuronal differentiation in neural stem cells and SH-SY5Y cells, ZFHX3 interacts with the chromatin remodelling BRG1/Brm-associated factor complex and the cleavage and polyadenylation complex. In line with a role for chromatin remodelling, ZFHX3 haploinsufficiency associates with a specific DNA methylation profile in leukocyte-derived DNA. The target genes of ZFHX3 are implicated in neuron and axon development. In Drosophila melanogaster , z fh2, considered to be the ZFHX3 orthologue, is expressed in the third instar larval brain. Ubiquitous and neuron-specific knockdown of zfh2 results in adult lethality underscoring a key role for zfh2 in development and neurodevelopment. Interestingly, ectopic expression of zfh2 as well as ZFHX3 in the developing wing disc results in a thoracic cleft phenotype. Collectively, our data shows that loss-of-function variants in ZFHX3 are a cause of syndromic ID, that associates with a specific DNA methylation profile. Furthermore, we show that ZFHX3 participates in chromatin remodelling and mRNA processing.

Tau polarizes an aging transcriptional signature to excitatory neurons and glia

Aging is a major risk factor for Alzheimer's disease (AD), and cell-type vulnerability underlies its characteristic clinical manifestations. We have performed longitudinal, single-cell RNA-sequencing in Drosophila with pan-neuronal expression of human tau, which forms AD neurofibrillary tangle pathology. Whereas tau- and aging-induced gene expression strongly overlap (93%), they differ in the affected cell types. In contrast to the broad impact of aging, tau-triggered changes are strongly polarized to excitatory neurons and glia. Further, tau can either activate or suppress innate immune gene expression signatures in a cell-type-specific manner. Integration of cellular abundance and gene expression pinpoints nuclear factor kappa B signaling in neurons as a marker for cellular vulnerability. We also highlight the conservation of cell-type-specific transcriptional patterns between Drosophila and human postmortem brain tissue. Overall, our results create a resource for dissection of dynamic, age-dependent gene expression changes at cellular resolution in a genetically tractable model of tauopathy.

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.

Genes-to-Pathways Species Conservation Analysis: Enabling the Exploration of Conservation of Biological Pathways and Processes Across Species

The last two decades have witnessed a strong momentum toward integration of cell-based and computational approaches in safety assessments. This is fueling a global regulatory paradigm shift toward reduction and replacement of the use of animals in toxicity tests while promoting the use of new approach methodologies. The understanding of conservation of molecular targets and pathways provides an opportunity to extrapolate effects across species and ultimately to determine the taxonomic applicability domain of assays and biological effects. Despite the wealth of genome-linked data available, there is a compelling need for improved accessibility, while ensuring that it reflects the underpinning biology. We present the novel pipeline Genes-to-Pathways Species Conservation Analysis (G2P-SCAN) to further support understanding on cross-species extrapolation of biological processes. This R package extracts, synthetizes, and structures the data available from different databases, that is, gene orthologs, protein families, entities, and reactions, linked to human genes and respective pathways across six relevant model species. The use of G2P-SCAN enables the overall analysis of orthology and functional families to substantiate the identification of conservation and susceptibility at the pathway level. In the present study we discuss five case studies, demonstrating the validity of the developed pipeline and its potential use as species extrapolation support. We foresee this pipeline will provide valuable biological insights and create space for the use of mechanistically based data to inform potential species susceptibility for research and safety decision purposes. Environ Toxicol Chem 2023;42:1152-1166. © 2023 UNILEVER GLOBAL IP LTD. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Characterizing the targets of transcription regulators by aggregating ChIP-seq and perturbation expression data sets

Mapping the gene targets of chromatin-associated transcription regulators (TRs) is a major goal of genomics research. ChIP-seq of TRs and experiments that perturb a TR and measure the differential abundance of gene transcripts are a primary means by which direct relationships are tested on a genomic scale. It has been reported that there is a poor overlap in the evidence across gene regulation strategies, emphasizing the need for integrating results from multiple experiments. Although research consortia interested in gene regulation have produced a valuable trove of high-quality data, there is an even greater volume of TR-specific data throughout the literature. In this study, we show a workflow for the identification, uniform processing, and aggregation of ChIP-seq and TR perturbation experiments for the ultimate purpose of ranking human and mouse TR-target interactions. Focusing on an initial set of eight regulators (ASCL1, HES1, MECP2, MEF2C, NEUROD1, PAX6, RUNX1, and TCF4), we identified 497 experiments suitable for analysis. We used this corpus to examine data concordance, to identify systematic patterns of the two data types, and to identify putative orthologous interactions between human and mouse. We build upon commonly used strategies to forward a procedure for aggregating and combining these two genomic methodologies, assessing these rankings against independent literature-curated evidence. Beyond a framework extensible to other TRs, our work also provides empirically ranked TR-target listings, as well as transparent experiment-level gene summaries for community use.

A phosphate-sensing organelle regulates phosphate and tissue homeostasis

Inorganic phosphate (Pi) is one of the essential molecules for life. However, little is known about intracellular Pi metabolism and signalling in animal tissues1. Following the observation that chronic Pi starvation causes hyperproliferation in the digestive epithelium of Drosophila melanogaster, we determined that Pi starvation triggers the downregulation of the Pi transporter PXo. In line with Pi starvation, PXo deficiency caused midgut hyperproliferation. Interestingly, immunostaining and ultrastructural analyses showed that PXo specifically marks non-canonical multilamellar organelles (PXo bodies). Further, by Pi imaging with a Förster resonance energy transfer (FRET)-based Pi sensor2, we found that PXo restricts cytosolic Pi levels. PXo bodies require PXo for biogenesis and undergo degradation following Pi starvation. Proteomic and lipidomic characterization of PXo bodies unveiled their distinct feature as an intracellular Pi reserve. Therefore, Pi starvation triggers PXo downregulation and PXo body degradation as a compensatory mechanism to increase cytosolic Pi. Finally, we identified connector of kinase to AP-1 (Cka), a component of the STRIPAK complex and JNK signalling3, as the mediator of PXo knockdown- or Pi starvation-induced hyperproliferation. Altogether, our study uncovers PXo bodies as a critical regulator of cytosolic Pi levels and identifies a Pi-dependent PXo-Cka-JNK signalling cascade controlling tissue homeostasis.

PTPA variants and impaired PP2A activity in early-onset parkinsonism with intellectual disability

The protein phosphatase 2A complex (PP2A), the major Ser/Thr phosphatase in the brain, is involved in a number of signalling pathways and functions, including the regulation of crucial proteins for neurodegeneration, such as alpha-synuclein, tau and LRRK2. Here, we report the identification of variants in the PTPA/PPP2R4 gene, encoding a major PP2A activator, in two families with early-onset parkinsonism and intellectual disability. We carried out clinical studies and genetic analyses, including genome-wide linkage analysis, whole-exome sequencing, and Sanger sequencing of candidate variants. We next performed functional studies on the disease-associated variants in cultured cells and knock-down of ptpa in Drosophila melanogaster. We first identified a homozygous PTPA variant, c.893T>G (p.Met298Arg), in patients from a South African family with early-onset parkinsonism and intellectual disability. Screening of a large series of additional families yielded a second homozygous variant, c.512C>A (p.Ala171Asp), in a Libyan family with a similar phenotype. Both variants co-segregate with disease in the respective families. The affected subjects display juvenile-onset parkinsonism and intellectual disability. The motor symptoms were responsive to treatment with levodopa and deep brain stimulation of the subthalamic nucleus. In overexpression studies, both the PTPA p.Ala171Asp and p.Met298Arg variants were associated with decreased PTPA RNA stability and decreased PTPA protein levels; the p.Ala171Asp variant additionally displayed decreased PTPA protein stability. Crucially, expression of both variants was associated with decreased PP2A complex levels and impaired PP2A phosphatase activation. PTPA orthologue knock-down in Drosophila neurons induced a significant impairment of locomotion in the climbing test. This defect was age-dependent and fully reversed by L-DOPA treatment. We conclude that bi-allelic missense PTPA variants associated with impaired activation of the PP2A phosphatase cause autosomal recessive early-onset parkinsonism with intellectual disability. Our findings might also provide new insights for understanding the role of the PP2A complex in the pathogenesis of more common forms of neurodegeneration.

Next-generation large-scale binary protein interaction network for Drosophila melanogaster

Generating reference maps of interactome networks illuminates genetic studies by providing a protein-centric approach to finding new components of existing pathways, complexes, and processes. We apply state-of-the-art methods to identify binary protein-protein interactions (PPIs) for Drosophila melanogaster. Four all-by-all yeast two-hybrid (Y2H) screens of > 10,000 Drosophila proteins result in the 'FlyBi' dataset of 8723 PPIs among 2939 proteins. Testing subsets of data from FlyBi and previous PPI studies using an orthogonal assay allows for normalization of data quality; subsequent integration of FlyBi and previous data results in an expanded binary Drosophila reference interaction network, DroRI, comprising 17,232 interactions among 6511 proteins. We use FlyBi data to generate an autophagy network, then validate in vivo using autophagy-related assays. The deformed wings (dwg) gene encodes a protein that is both a regulator and a target of autophagy. Altogether, these resources provide a foundation for building new hypotheses regarding protein networks and function.