The State of Systems Genetics in 2017

Systems genetics approaches for understanding complex traits with relevance for human disease

Quantitative traits are often complex because of the contribution of many loci, with further complexity added by environmental factors. In medical research, systems genetics is a powerful approach for the study of complex traits, as it integrates intermediate phenotypes, such as RNA, protein, and metabolite levels, to understand molecular and physiological phenotypes linking discrete DNA sequence variation to complex clinical and physiological traits. The primary purpose of this review is to describe some of the resources and tools of systems genetics in humans and rodent models, so that researchers in many areas of biology and medicine can make use of the data.

Leveraging genetic diversity to identify small molecules that reverse mouse skeletal muscle insulin resistance

Systems genetics has begun to tackle the complexity of insulin resistance by capitalising on computational advances to study high-diversity populations. 'Diversity Outbred in Australia (DOz)' is a population of genetically unique mice with profound metabolic heterogeneity. We leveraged this variance to explore skeletal muscle's contribution to whole-body insulin action through metabolic phenotyping and skeletal muscle proteomics of 215 DOz mice. Linear modelling identified 553 proteins that associated with whole-body insulin sensitivity (Matsuda Index) including regulators of endocytosis and muscle proteostasis. To enrich for causality, we refined this network by focusing on negatively associated, genetically regulated proteins, resulting in a 76-protein fingerprint of insulin resistance. We sought to perturb this network and restore insulin action with small molecules by integrating the Broad Institute Connectivity Map platform and in vitro assays of insulin action using the Prestwick chemical library. These complementary approaches identified the antibiotic thiostrepton as an insulin resistance reversal agent. Subsequent validation in ex vivo insulin-resistant mouse muscle and palmitate-induced insulin-resistant myotubes demonstrated potent insulin action restoration, potentially via upregulation of glycolysis. This work demonstrates the value of a drug-centric framework to validate systems-level analysis by identifying potential therapeutics for insulin resistance.

The potential of integrating human and mouse discovery platforms to advance our understanding of cardiometabolic diseases

Cardiometabolic diseases encompass a range of interrelated conditions that arise from underlying metabolic perturbations precipitated by genetic, environmental, and lifestyle factors. While obesity, dyslipidaemia, smoking, and insulin resistance are major risk factors for cardiometabolic diseases, individuals still present in the absence of such traditional risk factors, making it difficult to determine those at greatest risk of disease. Thus, it is crucial to elucidate the genetic, environmental, and molecular underpinnings to better understand, diagnose, and treat cardiometabolic diseases. Much of this information can be garnered using systems genetics, which takes population-based approaches to investigate how genetic variance contributes to complex traits. Despite the important advances made by human genome-wide association studies (GWAS) in this space, corroboration of these findings has been hampered by limitations including the inability to control environmental influence, limited access to pertinent metabolic tissues, and often, poor classification of diseases or phenotypes. A complementary approach to human GWAS is the utilisation of model systems such as genetically diverse mouse panels to study natural genetic and phenotypic variation in a controlled environment. Here, we review mouse genetic reference panels and the opportunities they provide for the study of cardiometabolic diseases and related traits. We discuss how the post-GWAS era has prompted a shift in focus from discovery of novel genetic variants to understanding gene function. Finally, we highlight key advantages and challenges of integrating complementary genetic and multi-omics data from human and mouse populations to advance biological discovery.

Does Therapeutic Repurposing in Cancer Meet the Expectations of Having Drugs at a Lower Price?

Therapeutic repurposing emerged as an alternative to the traditional drug discovery and development model (DDD) of new molecular entities (NMEs). It was anticipated that by being faster, safer, and cheaper, the development would result in lower-cost drugs. As defined in this work, a repurposed cancer drug is one approved by a health regulatory authority against a non-cancer indication that then gains new approval for cancer. With this definition, only three drugs are repurposed for cancer: Bacillus Calmette-Guerin (BCG) vaccine (superficial bladder cancer, thalidomide [multiple myeloma], and propranolol [infantile hemangioma]). Each of these has a different history regarding price and affordability, and it is not yet possible to generalize the impact of drug repurposing on the final price to the patient. However, the development, including the price, does not differ significantly from an NME. For the end consumer, the product's price is unrelated to whether it followed the classical development or repurposing. Economic constraints for clinical development, and drug prescription biases for repurposing drugs, are barriers yet to be overcome. The affordability of cancer drugs is a complex issue that varies from country to country. Many alternatives for having affordable drugs have been put forward, however these measures have thus far failed and are, at best, palliative. There are no immediate solutions to the problem of access to cancer drugs. It is necessary to critically analyze the impact of the current drug development model and be creative in implementing new models that genuinely benefit society.

Proteome-wide systems genetics identifies UFMylation as a regulator of skeletal muscle function

Improving muscle function has great potential to improve the quality of life. To identify novel regulators of skeletal muscle metabolism and function, we performed a proteomic analysis of gastrocnemius muscle from 73 genetically distinct inbred mouse strains, and integrated the data with previously acquired genomics and >300 molecular/phenotypic traits via quantitative trait loci mapping and correlation network analysis. These data identified thousands of associations between protein abundance and phenotypes and can be accessed online (https://muscle.coffeeprot.com/) to identify regulators of muscle function. We used this resource to prioritize targets for a functional genomic screen in human bioengineered skeletal muscle. This identified several negative regulators of muscle function including UFC1, an E2 ligase for protein UFMylation. We show UFMylation is up-regulated in a mouse model of amyotrophic lateral sclerosis, a disease that involves muscle atrophy. Furthermore, in vivo knockdown of UFMylation increased contraction force, implicating its role as a negative regulator of skeletal muscle function.

BAPST. A Combo of Common use drugs as metabolic therapy of cancer-a theoretical proposal

Advances in cancer therapy have yet to impact worldwide cancer mortality. Poor cancer drug affordability is one of the factors limiting mortality burden strikes. Up to now, cancer drug repurposing had no meet expectations concerning drug affordability. The three FDA-approved cancer drugs developed under repurposing -all-trans-retinoic acid, arsenic trioxide, and thalidomide- do not differ in price from other drugs developed under the classical model. Though additional factors affect the whole process from inception to commercialization, the repurposing of widely used, commercially available, and cheap drugs may help. This work reviews the concept of the malignant metabolic phenotype and its exploitation by simultaneously blocking key metabolic processes altered in cancer. We elaborate on a combination called BAPST, which stands for the following drugs and pathways they inhibit: Benserazide (glycolysis), Apomorphine (glutaminolysis), Pantoprazole (Fatty-acid synthesis), Simvastatin (mevalonate pathway), and Trimetazidine (Fatty-acid oxidation). Their respective primary indications are: • Parkinson's disease (benserazide and apomorphine). • Peptic ulcer disease (pantoprazole). • Hypercholesterolemia (simvastatin). • Ischemic heart disease (trimetazidine). When used for their primary indication, the literature review on each of these drugs shows they have a good safety profile and lack predicted pharmacokinetic interaction among them. Most importantly, the inhibitory enzymatic concentrations required for inhibiting their cancer targets enzymes are below the plasma concentrations observed when these drugs are used for their primary indication. Based on that, we propose that the regimen BAPTS merits preclinical testing.

AI models and the future of genomic research and medicine: True sons of knowledge?: Artificial intelligence needs to be integrated with causal conceptions in biomedicine to harness its societal benefits for the field

The increasing availability of large-scale, complex data has made research into how human genomes determine physiology in health and disease, as well as its application to drug development and medicine, an attractive field for artificial intelligence (AI) approaches. Looking at recent developments, we explore how such approaches interconnect and may conflict with needs for and notions of causal knowledge in molecular genetics and genomic medicine. We provide reasons to suggest that-while capable of generating predictive knowledge at unprecedented pace and scale-if and how these approaches will be integrated with prevailing causal concepts will not only determine the future of scientific understanding and self-conceptions in these fields. But these questions will also be key to develop differentiated policies, such as for education and regulation, in order to harness societal benefits of AI for genomic research and medicine.

Genome-environment association methods comparison supports omnigenic adaptation to ecological niche in malaria vector mosquitoes

The concept of a fundamental ecological niche is central to questions of geographic distribution, population demography, species conservation, and evolutionary potential. However, robust inference of genomic regions associated with evolutionary adaptation to particular environmental conditions remains difficult due to the myriad of potential confounding processes that can generate heterogeneous patterns of variation across the genome. Here, we interrogate the potential role of genome environment association (GEA) testing as an initial step in building an understanding of the genetic basis of ecological niche. We leverage publicly available genomic data from the Anopheles gambiae 1000 Genomes (Ag1000g) Consortium to test the ability of multiple analytically unique GEA methods to handle confounding patterns of genetic variation, control false positive rates, and discern associations with broadly relevant climate variables from random allele frequency patterns throughout the genome. We found evidence supporting the ability of commonly implemented GEA methods to account for confounding patterns of spatial and genetic variation, and control false positive rates. However, we fail to find evidence supporting the ability of GEA tests to reject signals of adaptation to randomly simulated environmental variables, indicating that discerning between true signals of genome environment adaptation and genome environment correlations resulting from alternative evolutionary processes, remains challenging. Because signals of environmental adaptation are so diffuse and confounded throughout the genome, we argue that genomic adaptation to ecological niche is likely best understood under an omnigenic model wherein highly interconnected, genome-wide gene regulatory networks shape genomic adaptation to key environmental conditions.

Perspectives on Drug Repurposing

Complex common diseases are a significant burden for our societies and demand not only preventive measures but also more effective, safer, and more affordable treatments. The whole process of the current model of drug discovery and development implies a high investment by the pharmaceutical industry, which ultimately impact in high drug prices. In this sense, drug repurposing would help meet the needs of patients to access useful and novel treatments. Unlike the traditional approach, drug repurposing enters both the preclinical evaluation and clinical trials of the compound of interest faster, budgeting research and development costs, and limiting potential biosafety risks. The participation of government, society, and private investors is needed to secure the funds for experimental design and clinical development of repurposing candidates to have affordable, effective, and safe repurposed drugs. Moreover, extensive advertising of repurposing as a concept in the health community, could reduce prescribing bias when enough clinical evidence exists, which will support the employment of cheaper and more accessible repurposed compounds for common conditions.

Mapping Robust Genetic Variants Associated with Exercise Responses

This review summarised robust and consistent genetic variants associated with aerobic-related and resistance-related phenotypes. In total we highlight 12 SNPs and 7 SNPs that are robustly associated with variance in aerobic-related and resistance-related phenotypes respectively. To date, there is very little literature ascribed to understanding the interplay between genes and environmental factors and the development of physiological traits. We discuss future directions, including large-scale exercise studies to elucidate the functional relevance of the discovered genomic markers. This approach will allow more rigour and reproducible research in the field of exercise genomics.

A trans-eQTL network regulates osteoclast multinucleation and bone mass

Functional characterisation of cell-type-specific regulatory networks is key to establish a causal link between genetic variation and phenotype. The osteoclast offers a unique model for interrogating the contribution of co-regulated genes to in vivo phenotype as its multinucleation and resorption activities determine quantifiable skeletal traits. Here we took advantage of a trans-regulated gene network (MMnet, macrophage multinucleation network) which we found to be significantly enriched for GWAS variants associated with bone-related phenotypes. We found that the network hub gene Bcat1 and seven other co-regulated MMnet genes out of 13, regulate bone function. Specifically, global (Pik3cb^-/-, Atp8b2^+/-, Igsf8^-/-, Eml1^-/-, Appl2^-/-, Deptor^-/-) and myeloid-specific Slc40a1 knockout mice displayed abnormal bone phenotypes. We report opposing effects of MMnet genes on bone mass in mice and osteoclast multinucleation/resorption in humans with strong correlation between the two. These results identify MMnet as a functionally conserved network that regulates osteoclast multinucleation and bone mass.

Systems Genetics for Mechanistic Discovery in Heart Diseases

Cardiovascular diseases are the leading cause of death worldwide. Complex diseases with highly heterogenous disease progression among patient populations, cardiovascular diseases feature multifactorial contributions from both genetic and environmental stressors. Despite significant effort utilizing multiple approaches from molecular biology to genome-wide association studies, the genetic landscape of cardiovascular diseases, particularly for the nonfamilial forms of heart failure, is still poorly understood. In the past decade, systems-level approaches based on omics technologies have become an important approach for the study of complex traits in large populations. These advances create opportunities to integrate genetic variation with other biological layers to identify and prioritize candidate genes, understand pathogenic pathways, and elucidate gene-gene and gene-environment interactions. In this review, we will highlight some of the recent progress made using systems genetics approaches to uncover novel mechanisms and molecular bases of cardiovascular pathophysiological manifestations. The key technology and data analysis platforms necessary to implement systems genetics will be described, and the current major challenges and future directions will also be discussed. For complex cardiovascular diseases, such as heart failure, systems genetics represents a powerful strategy to obtain mechanistic insights and to develop individualized diagnostic and therapeutic regiments, paving the way for precision cardiovascular medicine.

A multi-omics digital research object for the genetics of sleep regulation

With the aim to uncover the molecular pathways underlying the regulation of sleep, we recently assembled an extensive and comprehensive systems genetics dataset interrogating a genetic reference population of mice at the levels of the genome, the brain and liver transcriptomes, the plasma metabolome, and the sleep-wake phenome. To facilitate a meaningful and efficient re-use of this public resource by others we designed, describe in detail, and made available a Digital Research Object (DRO), embedding data, documentation, and analytics. We present and discuss both the advantages and limitations of our multi-modal resource and analytic pipeline. The reproducibility of the results was tested by a bioinformatician not implicated in the original project and the robustness of results was assessed by re-annotating genetic and transcriptome data from the mm9 to the mm10 mouse genome assembly.

Common Disease Is More Complex Than Implied by the Core Gene Omnigenic Model

The evidence that most adult-onset common diseases have a polygenic genetic architecture fully consistent with robust biological systems supported by multiple back-up mechanisms is now overwhelming. In this context, we consider the recent "omnigenic" or "core genes" model. A key assumption of the model is that there is a relatively small number of core genes relevant to any disease. While intuitively appealing, this model may underestimate the biological complexity of common disease, and therefore, the goal to discover core genes should not guide experimental design. We consider other implications of polygenicity, concluding that a focus on patient stratification is needed to achieve the goals of precision medicine.

Systems genetics identifies a macrophage cholesterol network associated with physiological wound healing

Among other cells, macrophages regulate the inflammatory and reparative phases during wound healing but genetic determinants and detailed molecular pathways that modulate these processes are not fully elucidated. Here, we took advantage of normal variation in wound healing in 1,378 genetically outbred mice, and carried out macrophage RNA-sequencing profiling of mice with extreme wound healing phenotypes (i.e., slow and fast healers, n = 146 in total). The resulting macrophage coexpression networks were genetically mapped and led to the identification of a unique module under strong trans-acting genetic control by the Runx2 locus. This macrophage-mediated healing network was specifically enriched for cholesterol and fatty acid biosynthetic processes. Pharmacological blockage of fatty acid synthesis with cerulenin resulted in delayed wound healing in vivo, and increased macrophage infiltration in the wounded skin, suggesting the persistence of an unresolved inflammation. We show how naturally occurring sequence variation controls transcriptional networks in macrophages, which in turn regulate specific metabolic pathways that could be targeted in wound healing.

A systems genetics resource and analysis of sleep regulation in the mouse

Sleep is essential for optimal brain functioning and health, but the biological substrates through which sleep delivers these beneficial effects remain largely unknown. We used a systems genetics approach in the BXD genetic reference population (GRP) of mice and assembled a comprehensive experimental knowledge base comprising a deep "sleep-wake" phenome, central and peripheral transcriptomes, and plasma metabolome data, collected under undisturbed baseline conditions and after sleep deprivation (SD). We present analytical tools to interactively interrogate the database, visualize the molecular networks altered by sleep loss, and prioritize candidate genes. We found that a one-time, short disruption of sleep already extensively reshaped the systems genetics landscape by altering 60%-78% of the transcriptomes and the metabolome, with numerous genetic loci affecting the magnitude and direction of change. Systems genetics integrative analyses drawing on all levels of organization imply α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptor trafficking and fatty acid turnover as substrates of the negative effects of insufficient sleep. Our analyses demonstrate that genetic heterogeneity and the effects of insufficient sleep itself on the transcriptome and metabolome are far more widespread than previously reported.

Interrogation of Mammalian Protein Complex Structure, Function, and Membership Using Genome-Scale Fitness Screens

Protein complexes are assemblies of subunits that have co-evolved to execute one or many coordinated functions in the cellular environment. Functional annotation of mammalian protein complexes is critical to understanding biological processes, as well as disease mechanisms. Here, we used genetic co-essentiality derived from genome-scale RNAi- and CRISPR-Cas9-based fitness screens performed across hundreds of human cancer cell lines to assign measures of functional similarity. From these measures, we systematically built and characterized functional similarity networks that recapitulate known structural and functional features of well-studied protein complexes and resolve novel functional modules within complexes lacking structural resolution, such as the mammalian SWI/SNF complex. Finally, by integrating functional networks with large protein-protein interaction networks, we discovered novel protein complexes involving recently evolved genes of unknown function. Taken together, these findings demonstrate the utility of genetic perturbation screens alone, and in combination with large-scale biophysical data, to enhance our understanding of mammalian protein complexes in normal and disease states.

Protein Moonlighting Revealed by Noncatalytic Phenotypes of Yeast Enzymes

A single gene can partake in several biological processes, and therefore gene deletions can lead to different-sometimes unexpected-phenotypes. However, it is not always clear whether such pleiotropy reflects the loss of a unique molecular activity involved in different processes or the loss of a multifunctional protein. Here, using Saccharomyces cerevisiae metabolism as a model, we systematically test the null hypothesis that enzyme phenotypes depend on a single annotated molecular function, namely their catalysis. We screened a set of carefully selected genes by quantifying the contribution of catalysis to gene deletion phenotypes under different environmental conditions. While most phenotypes were explained by loss of catalysis, slow growth was readily rescued by a catalytically inactive protein in about one-third of the enzymes tested. Such noncatalytic phenotypes were frequent in the Alt1 and Bat2 transaminases and in the isoleucine/valine biosynthetic enzymes Ilv1 and Ilv2, suggesting novel "moonlighting" activities in these proteins. Furthermore, differential genetic interaction profiles of gene deletion and catalytic mutants indicated that ILV1 is functionally associated with regulatory processes, specifically to chromatin modification. Our systematic study shows that gene loss phenotypes and their genetic interactions are frequently not driven by the loss of an annotated catalytic function, underscoring the moonlighting nature of cellular metabolism.