Statistical integration of multi-omics and drug screening data from cell lines.

Data integration methods are used to obtain a unified summary of multiple datasets. For multi-modal data, we propose a computational workflow to jointly analyze datasets from cell lines. The workflow comprises a novel probabilistic data integration method, named POPLS-DA, for multi-omics data. The workflow is motivated by a study on synucleinopathies where transcriptomics, proteomics, and drug screening data are measured in affected LUHMES cell lines and controls. The aim is to highlight potentially druggable pathways and genes involved in synucleinopathies. First, POPLS-DA is used to prioritize genes and proteins that best distinguish cases and controls. For these genes, an integrated interaction network is constructed where the drug screen data is incorporated to highlight druggable genes and pathways in the network. Finally, functional enrichment analyses are performed to identify clusters of synaptic and lysosome-related genes and proteins targeted by the protective drugs. POPLS-DA is compared to other single- and multi-omics approaches. We found that HSPA5, a member of the heat shock protein 70 family, was one of the most targeted genes by the validated drugs, in particular by AT1-blockers. HSPA5 and AT1-blockers have been previously linked to α-synuclein pathology and Parkinson's disease, showing the relevance of our findings. Our computational workflow identified new directions for therapeutic targets for synucleinopathies. POPLS-DA provided a larger interpretable gene set than other single- and multi-omic approaches. An implementation based on R and markdown is freely available online.

Loss of a gluconeogenic muscle enzyme contributed to adaptive metabolic traits in hummingbirds.

Hummingbirds possess distinct metabolic adaptations to fuel their energy-demanding hovering flight, but the underlying genomic changes are largely unknown. Here, we generated a chromosome-level genome assembly of the long-tailed hermit and screened for genes that have been specifically inactivated in the ancestral hummingbird lineage. We discovered that FBP2 (fructose-bisphosphatase 2), which encodes a gluconeogenic muscle enzyme, was lost during a time period when hovering flight evolved. We show that FBP2 knockdown in an avian muscle cell line up-regulates glycolysis and enhances mitochondrial respiration, coincident with an increased mitochondria number. Furthermore, genes involved in mitochondrial respiration and organization have up-regulated expression in hummingbird flight muscle. Together, these results suggest that FBP2 loss was likely a key step in the evolution of metabolic muscle adaptations required for true hovering flight.

Phenotype loss is associated with widespread divergence of the gene regulatory landscape in evolution.

Detecting the genomic changes underlying phenotypic changes between species is a main goal of evolutionary biology and genomics. Evolutionary theory predicts that changes in cis-regulatory elements are important for morphological changes. We combined genome sequencing, functional genomics and genome-wide comparative analyses to investigate regulatory elements in lineages that lost morphological traits. We first show that limb loss in snakes is associated with widespread divergence of limb regulatory elements. We next show that eye degeneration in subterranean mammals is associated with widespread divergence of eye regulatory elements. In both cases, sequence divergence results in an extensive loss of transcription factor binding sites. Importantly, diverged regulatory elements are associated with genes required for normal limb patterning or normal eye development and function, suggesting that regulatory divergence contributed to the loss of these phenotypes. Together, our results show that genome-wide decay of the phenotype-specific cis-regulatory landscape is a hallmark of lost morphological traits.

Isogenic FUS-eGFP iPSC Reporter Lines Enable Quantification of FUS Stress Granule Pathology that Is Rescued by Drugs Inducing Autophagy.

Perturbations in stress granule (SG) dynamics may be at the core of amyotrophic lateral sclerosis (ALS). Since SGs are membraneless compartments, modeling their dynamics in human motor neurons has been challenging, thus hindering the identification of effective therapeutics. Here, we report the generation of isogenic induced pluripotent stem cells carrying wild-type and P525L FUS-eGFP. We demonstrate that FUS-eGFP is recruited into SGs and that P525L profoundly alters their dynamics. With a screening campaign, we demonstrate that PI3K/AKT/mTOR pathway inhibition increases autophagy and ameliorates SG phenotypes linked to P525L FUS by reducing FUS-eGFP recruitment into SGs. Using a Drosophila model of FUS-ALS, we corroborate that induction of autophagy significantly increases survival. Finally, by screening clinically approved drugs for their ability to ameliorate FUS SG phenotypes, we identify a number of brain-penetrant anti-depressants and anti-psychotics that also induce autophagy. These drugs could be repurposed as potential ALS treatments.

Protective efficacy of phosphodiesterase-1 inhibition against alpha-synuclein toxicity revealed by compound screening in LUHMES cells.

α-synuclein-induced neurotoxicity is a core pathogenic event in neurodegenerative synucleinopathies such as Parkinson's disease, dementia with Lewy bodies, or multiple system atrophy. There is currently no disease-modifying therapy available for these diseases. We screened 1,600 FDA-approved drugs for their efficacy to protect LUHMES cells from degeneration induced by wild-type α-synuclein and identified dipyridamole, a non-selective phosphodiesterase inhibitor, as top hit. Systematic analysis of other phosphodiesterase inhibitors identified a specific phosphodiesterase 1 inhibitor as most potent to rescue from α-synuclein toxicity. Protection was mediated by an increase of cGMP and associated with the reduction of a specific α-synuclein oligomeric species. RNA interference experiments confirmed PDE1A and to a smaller extent PDE1C as molecular targets accounting for the protective efficacy. PDE1 inhibition also rescued dopaminergic neurons from wild-type α-synuclein induced degeneration in the substantia nigra of mice. In conclusion, this work identifies inhibition of PDE1A in particular as promising target for neuroprotective treatment of synucleinopathies.

A conserved RhoGAP limits M phase contractility and coordinates with microtubule asters to confine RhoA during cytokinesis.

During animal cell cytokinesis, the spindle directs contractile ring assembly by activating RhoA in a narrow equatorial zone. Rapid GTPase activating protein (GAP)-mediated inactivation (RhoA flux) is proposed to limit RhoA zone dimensions. Testing the significance of RhoA flux has been hampered by the fact that the GAP targeting RhoA is not known. Here, we identify M phase GAP (MP-GAP) as the primary GAP targeting RhoA during mitosis and cytokinesis. MP-GAP inhibition caused excessive RhoA activation in M phase, leading to the uncontrolled formation of large cortical protrusions and late cytokinesis failure. RhoA zone width was broadened by attenuation of the centrosomal asters but was not affected by MP-GAP inhibition alone. Simultaneous aster attenuation and MP-GAP inhibition led to RhoA accumulation around the entire cell periphery. These results identify the major GAP restraining RhoA during cell division and delineate the relative contributions of RhoA flux and centrosomal asters in controlling RhoA zone dimensions.