Multilayered genetic and omics dissection of mitochondrial activity in a mouse reference population

The lysine degradation pathway: Subcellular compartmentalization and enzyme deficiencies

Lysine degradation via formation of saccharopine is a pathway confined to the mitochondria. The second pathway for lysine degradation, the pipecolic acid pathway, is not yet fully elucidated and known enzymes are localized in the mitochondria, cytosol and peroxisome. The tissue-specific roles of these two pathways are still under investigation. The lysine degradation pathway is clinically relevant due to the occurrence of two severe neurometabolic disorders, pyridoxine-dependent epilepsy (PDE) and glutaric aciduria type 1 (GA1). The existence of three other disorders affecting lysine degradation without apparent clinical consequences opens up the possibility to find alternative therapeutic strategies for PDE and GA1 through pathway modulation. A better understanding of the mechanisms, compartmentalization and interplay between the different enzymes and metabolites involved in lysine degradation is of utmost importance.

Inhibition and Crystal Structure of the Human DHTKD1-Thiamin Diphosphate Complex

DHTKD1 is the E1 component of the 2-oxoadipate dehydrogenase complex, which is an enzyme involved in the catabolism of (hydroxy-)lysine and tryptophan. Mutations in DHTKD1 have been associated with 2-aminoadipic and 2-oxoadipic aciduria, Charcot-Marie-Tooth disease type 2Q and eosinophilic esophagitis, but the pathophysiology of these clinically distinct disorders remains elusive. Here, we report the identification of adipoylphosphonic acid and tenatoprazole as DHTKD1 inhibitors using targeted and high throughput screening, respectively. We furthermore elucidate the DHTKD1 crystal structure with thiamin diphosphate bound at 2.25 Å. We also report the impact of 10 disease-associated missense mutations on DHTKD1. Whereas the majority of the DHTKD1 variants displayed impaired folding or reduced thermal stability in combination with absent or reduced enzyme activity, three variants showed no abnormalities. Our work provides chemical and structural tools for further understanding of the function of DHTKD1 and its role in several human pathologies.

A biochemical network modeling of a whole-cell

All cellular processes can be ultimately understood in terms of respective fundamental biochemical interactions between molecules, which can be modeled as networks. Very often, these molecules are shared by more than one process, therefore interconnecting them. Despite this effect, cellular processes are usually described by separate networks with heterogeneous levels of detail, such as metabolic, protein-protein interaction, and transcription regulation networks. Aiming at obtaining a unified representation of cellular processes, we describe in this work an integrative framework that draws concepts from rule-based modeling. In order to probe the capabilities of the framework, we used an organism-specific database and genomic information to model the whole-cell biochemical network of the Mycoplasma genitalium organism. This modeling accounted for 15 cellular processes and resulted in a single component network, indicating that all processes are somehow interconnected. The topological analysis of the network showed structural consistency with biological networks in the literature. In order to validate the network, we estimated gene essentiality by simulating gene deletions and compared the results with experimental data available in the literature. We could classify 212 genes as essential, being 95% of them consistent with experimental results. Although we adopted a relatively simple organism as a case study, we suggest that the presented framework has the potential for paving the way to more integrated studies of whole organisms leading to a systemic analysis of cells on a broader scale. The modeling of other organisms using this framework could provide useful large-scale models for different fields of research such as bioengineering, network biology, and synthetic biology, and also provide novel tools for medical and industrial applications.

Skeletal muscle enhancer interactions identify genes controlling whole-body metabolism

Obesity and type 2 diabetes (T2D) are metabolic disorders influenced by lifestyle and genetic factors that are characterized by insulin resistance in skeletal muscle, a prominent site of glucose disposal. Numerous genetic variants have been associated with obesity and T2D, of which the majority are located in non-coding DNA regions. This suggests that most variants mediate their effect by altering the activity of gene-regulatory elements, including enhancers. Here, we map skeletal muscle genomic enhancer elements that are dynamically regulated after exposure to the free fatty acid palmitate or the inflammatory cytokine TNFα. By overlapping enhancer positions with the location of disease-associated genetic variants, and resolving long-range chromatin interactions between enhancers and gene promoters, we identify target genes involved in metabolic dysfunction in skeletal muscle. The majority of these genes also associate with altered whole-body metabolic phenotypes in the murine BXD genetic reference population. Thus, our combined genomic investigations identified genes that are involved in skeletal muscle metabolism.

Obesity and type 2 diabetes (T2D) are metabolic disorders characterized by insulin resistance in skeletal muscle. Here, the authors map skeletal muscle enhancer elements dynamically regulated after exposure to free fatty acid palmitate or inflammatory cytokine TNFα and identify target genes involved in metabolic dysfunction in skeletal muscle.

Systemic effects of mitochondrial stress

Multicellular organisms are complex biological systems, composed of specialized tissues that require coordination of the metabolic and fitness state of each component. In the cells composing the tissues, one central organelle is the mitochondrion, a compartment essential for many energetic and fundamental biological processes. Beyond serving these functions, mitochondria have emerged as signaling hubs in biological systems, capable of inducing changes to the cell they are in, to cells in distal tissues through secreted factors, and to overall animal physiology. Here, we describe our current understanding of these communication mechanisms in the context of mitochondrial stress. We focus on cellular mechanisms that deal with perturbations to the mitochondrial proteome and outline recent advances in understanding how local perturbations can affect distal tissues and animal physiology in model organisms. Finally, we discuss recent findings of these responses associated with metabolic and age-associated diseases in mammalian systems, and how they may be employed as diagnostic and therapeutic tools.

Omics Integration for Mitochondria Systems Biology

Significance: Elucidation of the central importance of mitophagy in homeostasis of cells and organisms emphasizes that mitochondrial functions extend far beyond short-term needs for energy production. In mitochondria systems biology, the mitochondrial genome, proteome, and metabolome operate as a functional network in coordination of cell activities. Organization occurs through subnetworks that are interconnected by membrane potential, transport activities, allosteric and cooperative interactions, redox signaling mechanisms, rheostatic control by post-translational modifications, and metal ion homeostasis. These subnetworks enable use of varied energy precursors, defense against environmental stressors, and macromolecular rewiring to titrate energy production, biosynthesis, and detoxification according to cell-specific needs. Rewiring mechanisms, termed mitochondrial reprogramming, enhance fitness to respond to metabolic resources and challenges from the environment. Maladaptive responses can cause cell death. Maladaptive rewiring can cause disease. In cancer, adaptive rewiring can interfere with effective treatment. Recent Advances: Many recent advances have been facilitated by the development of new omics tools, which create opportunities to use data-driven analysis of omics data to address these complex adaptive and maladaptive mechanisms of mitochondrial reprogramming in human disease. Critical Issues: Application of omics integration to model systems reveals a critical role for metal ion homeostasis broadly impacting mitochondrial reprogramming. Importantly, data show that trans-omics associations are more robust and biologically relevant than single omics associations. Future Directions: Application of omics integration to mitophagy research creates new opportunities to link the complex, interactive functions of mitochondrial form and function in mitochondria systems biology.

A genetic screen identifies Crat as a regulator of pancreatic beta-cell insulin secretion

Objectives

Glucose-stimulated insulin secretion is a critical function in the regulation of glucose homeostasis, and its deregulation is associated with the development of type 2 diabetes. Here, we performed a genetic screen using islets isolated from the BXD panel of advanced recombinant inbred (RI) lines of mice to search for novel regulators of insulin production and secretion.

Methods

Pancreatic islets were isolated from 36 RI BXD lines and insulin secretion was measured following exposure to 2.8 or 16.7 mM glucose with or without exendin-4. Islets from the same RI lines were used for RNA extraction and transcript profiling. Quantitative trait loci (QTL) mapping was performed for each secretion condition and combined with transcriptome data to prioritize candidate regulatory genes within the identified QTL regions. Functional studies were performed by mRNA silencing or overexpression in MIN6B1 cells and by studying mice and islets with beta-cell-specific gene inactivation.

Results

Insulin secretion under the 16.7 mM glucose plus exendin-4 condition was mapped significantly to a chromosome 2 QTL. Within this QTL, RNA-Seq data prioritized Crat (carnitine O-acetyl transferase) as a strong candidate regulator of the insulin secretion trait. Silencing Crat expression in MIN6B1 cells reduced insulin content and insulin secretion by ∼30%. Conversely, Crat overexpression enhanced insulin content and secretion by ∼30%. When islets from mice with beta-cell-specific Crat inactivation were exposed to high glucose, they displayed a 30% reduction of insulin content as compared to control islets. We further showed that decreased Crat expression in both MIN6B1 cells and pancreatic islets reduced the oxygen consumption rate in a glucose concentration-dependent manner.

Conclusions

We identified Crat as a regulator of insulin secretion whose action is mediated by an effect on total cellular insulin content; this effect also depends on the genetic background of the RI mouse lines. These data also show that in the presence of the stimulatory conditions used the insulin secretion rate is directly related to the insulin content.

•

A QTL analysis in BXD mice identifies Crat as a regulator of insulin secretion.

•

Crat regulates insulin content in MIN6B1 cells and pancreatic islets.

•

Crat regulates glucose oxidation in MIN6B1 cells and pancreatic islets.

•

Crat links glucose metabolism to the control of beta-cell insulin content.

•

Insulin content limits insulin secretion in response to high glucose and exendin-4 level.

An Integrated Global Analysis of Compartmentalized HRAS Signaling

Modern omics technologies allow us to obtain global information on different types of biological networks. However, integrating these different types of analyses into a coherent framework for a comprehensive biological interpretation remains challenging. Here, we present a conceptual framework that integrates protein interaction, phosphoproteomics, and transcriptomics data. Applying this method to analyze HRAS signaling from different subcellular compartments shows that spatially defined networks contribute specific functions to HRAS signaling. Changes in HRAS protein interactions at different sites lead to different kinase activation patterns that differentially regulate gene transcription. HRAS-mediated signaling is the strongest from the cell membrane, but it regulates the largest number of genes from the endoplasmic reticulum. The integrated networks provide a topologically and functionally resolved view of HRAS signaling. They reveal distinct HRAS functions including the control of cell migration from the endoplasmic reticulum and TP53-dependent cell survival when signaling from the Golgi apparatus.

An adipocyte-specific defect in oxidative phosphorylation increases systemic energy expenditure and protects against diet-induced obesity in mouse models

AIMS/HYPOTHESIS

Mitochondrial oxidative phosphorylation (OxPhos) is essential for energy production and survival. However, the tissue-specific and systemic metabolic effects of OxPhos function in adipocytes remain incompletely understood.

METHODS

We used adipocyte-specific Crif1 (also known as Gadd45gip1) knockout (AdKO) mice with decreased adipocyte OxPhos function. AdKO mice fed a normal chow or high-fat diet were evaluated for glucose homeostasis, weight gain and energy expenditure (EE). RNA sequencing of adipose tissues was used to identify the key mitokines affected in AdKO mice, which included fibroblast growth factor 21 (FGF21) and growth differentiation factor 15 (GDF15). For in vitro analysis, doxycycline was used to pharmacologically decrease OxPhos in 3T3L1 adipocytes. To identify the effects of GDF15 and FGF21 on the metabolic phenotype of AdKO mice, we generated AdKO mice with global Gdf15 knockout (AdGKO) or global Fgf21 knockout (AdFKO).

RESULTS

Under high-fat diet conditions, AdKO mice were resistant to weight gain and exhibited higher EE and improved glucose tolerance. In vitro pharmacological and in vivo genetic inhibition of OxPhos in adipocytes significantly upregulated mitochondrial unfolded protein response-related genes and secretion of mitokines such as GDF15 and FGF21. We evaluated the metabolic phenotypes of AdGKO and AdFKO mice, revealing that GDF15 and FGF21 differentially regulated energy homeostasis in AdKO mice. Both mitokines had beneficial effects on obesity and insulin resistance in the context of decreased adipocyte OxPhos, but only GDF15 regulated EE in AdKO mice.

CONCLUSIONS/INTERPRETATION

The present study demonstrated that the adipose tissue adaptive mitochondrial stress response affected systemic energy homeostasis via cell-autonomous and non-cell-autonomous pathways. We identified novel roles for adipose OxPhos and adipo-mitokines in the regulation of systemic glucose homeostasis and EE, which facilitated adaptation of an organism to local mitochondrial stress.

Proteomic and interactomic insights into the molecular basis of cell functional diversity

The ability of living systems to adapt to changing conditions originates from their capacity to change their molecular constitution. This is achieved by multiple mechanisms that modulate the quantitative composition and the diversity of the molecular inventory. Molecular diversification is particularly pronounced on the proteome level, at which multiple proteoforms derived from the same gene can in turn combinatorially form different protein complexes, thus expanding the repertoire of functional modules in the cell. The study of molecular and modular diversity and their involvement in responses to changing conditions has only recently become possible through the development of new 'omics'-based screening technologies. This Review explores our current knowledge of the mechanisms regulating functional diversification along the axis of gene expression, with a focus on the proteome and interactome. We explore the interdependence between different molecular levels and how this contributes to functional diversity. Finally, we highlight several recent techniques for studying molecular diversity, with specific focus on mass spectrometry-based analysis of the proteome and its organization into functional modules, and examine future directions for this rapidly growing field.

Treatment With Lipopolysaccharide Induces Distinct Changes in Metabolite Profile and Body Weight in 129Sv and Bl6 Mouse Strains

Mouse strains differ significantly in their behaviors and responses to pathogenic and pharmacological agents. This study seeks to characterize behavioral and metabolomic profiles of two widely used mouse lines, 129S6/SvEvTac (129Sv) and C57BL/6NTac (Bl6), to acute administration of lipopolysaccharide (LPS). LPS caused a significant suppression of locomotor activity and a decline in body weight (BW) in both strains within 24 h. However, the BW loss was more pronounced in Bl6 than in 129Sv. Comparison of strains revealed clear differences between their metabolomic profiles. According to the general linear model analysis (GLM), the 1.5 h LPS challenge in Bl6 caused a decrease of propionylcarnitine (C3), glucogenic amino acids, and acetylornithine (Ac-Orn), whereas the response of 129Sv included decreased concentrations of short-chain acylcarnitines (SCACs), citrulline, and elevation of glycerophospholipid (PCaa C42:0) and sphingolipid [SM(OH)C16:1]. 24 h after LPS administration, robust alterations in lipid profile were observed in both strains. LPS treatment caused elevation of sphingolipids, phosphatidylcholine diacyls (PCaa) as well as a decrease in lysophosphatidylcholines (LysoPC). However, the number of elevated PCaa and sphingolipids was considerably higher in 129Sv. In addition to lipids, 24 h LPS challenge in Bl6 mice induced increased levels of kynurenine (KYN), putrescine and decreased levels of citrulline, hexoses, Ac-Orn, and PC acyl-alkyl (PCae 38:2) as well as severe BW loss. In contrast, the 24 h LPS challenge in 129Sv mice induced increased levels of KYN, long-chain acylcarnitines (LCACs) and decreased levels of citrulline as well as moderate BW loss. Altogether, our study revealed both similarities and differences in response to LPS in Bl6 and 129Sv strains. For major differences, Bl6 mice showed stronger reduction of BW 24 h after LPS treatment, accompanied by significantly reduced levels of hexoses, the ratio between LysoPC16:1/LysoPC16:0, and elevated levels of neuroprotective putrescine. In 129Sv mice, the BW loss was milder, accompanied by increased levels of hydroxylated LCACs, probably reflecting shifts in oxidative metabolism of fatty acids. One may suggest that LPS caused stronger hypometabolic state in the Bl6 mice than in the 129Sv strain. Altogether, this study confirms that Bl6 and 129Sv mice display vastly distinct adaptation capacities independent from the nature of stressful challenge.

The Natural Metabolite 4-Cresol Improves Glucose Homeostasis and Enhances β-Cell Function

Exposure to natural metabolites contributes to the risk of cardiometabolic diseases (CMDs). Through metabolome profiling, we identify the inverse correlation between serum concentrations of 4-cresol and type 2 diabetes. The chronic administration of non-toxic doses of 4-cresol in complementary preclinical models of CMD reduces adiposity, glucose intolerance, and liver triglycerides, enhances insulin secretion in vivo, stimulates islet density and size, and pancreatic β-cell proliferation, and increases vascularization, suggesting activated islet enlargement. In vivo insulin sensitivity is not affected by 4-cresol. The incubation of mouse isolated islets with 4-cresol results in enhanced insulin secretion, insulin content, and β-cell proliferation of a magnitude similar to that induced by GLP-1. In both CMD models and isolated islets, 4-cresol is associated with the downregulated expression of the kinase DYRK1A, which may mediate its biological effects. Our findings identify 4-cresol as an effective regulator of β-cell function, which opens up perspectives for therapeutic applications in syndromes of insulin deficiency.

Mouse Systems Genetics as a Prelude to Precision Medicine

Mouse models have been instrumental in understanding human disease biology and proposing possible new treatments. The precise control of the environment and genetic composition of mice allows more rigorous observations, but limits the generalizability and translatability of the results into human applications. In the era of precision medicine, strategies using mouse models have to be revisited to effectively emulate human populations. Systems genetics is one promising paradigm that may promote the transition to novel precision medicine strategies. Here, we review the state-of-the-art resources and discuss how mouse systems genetics helps to understand human diseases and to advance the development of precision medicine, with an emphasis on the existing resources and strategies.

Synthetic analogues of 2-oxo acids discriminate metabolic contribution of the 2-oxoglutarate and 2-oxoadipate dehydrogenases in mammalian cells and tissues

The biological significance of the DHTKD1-encoded 2-oxoadipate dehydrogenase (OADH) remains obscure due to its catalytic redundancy with the ubiquitous OGDH-encoded 2-oxoglutarate dehydrogenase (OGDH). In this work, metabolic contributions of OADH and OGDH are discriminated by exposure of cells/tissues with different DHTKD1 expression to the synthesized phosphonate analogues of homologous 2-oxodicarboxylates. The saccharopine pathway intermediates and phosphorylated sugars are abundant when cellular expressions of DHTKD1 and OGDH are comparable, while nicotinate and non-phosphorylated sugars are when DHTKD1 expression is order(s) of magnitude lower than that of OGDH. Using succinyl, glutaryl and adipoyl phosphonates on the enzyme preparations from tissues with varied DHTKD1 expression reveals the contributions of OADH and OGDH to oxidation of 2-oxoadipate and 2-oxoglutarate in vitro. In the phosphonates-treated cells with the high and low DHTKD1 expression, adipate or glutarate, correspondingly, are the most affected metabolites. The marker of fatty acid β-oxidation, adipate, is mostly decreased by the shorter, OGDH-preferring, phosphonate, in agreement with the known OGDH dependence of β-oxidation. The longest, OADH-preferring, phosphonate mostly affects the glutarate level. Coupled decreases in sugars and nicotinate upon the OADH inhibition link the perturbation in glucose homeostasis, known in OADH mutants, to the nicotinate-dependent NAD metabolism.

Lipoatrophy-Associated Insulin Resistance and Hepatic Steatosis are Attenuated by Intake of Diet Rich in Omega 3 Fatty Acids

SCOPE

Glucose homeostasis and progression of nonalcoholic fatty liver disease (NAFLD) and hepatomegaly in severe lipoatrophic mice and their modulation by intake of a diet rich in omega 3 (n-3) fatty acids (HFO) are evaluated.

METHODS AND RESULTS

Severe lipoatrophic mice induced by PPAR-γ deletion exclusively in adipocytes (A-PPARγ KO) and littermate controls (A-PPARγ WT) are evaluated for glucose homeostasis and liver mass, proteomics, lipidomics, inflammation, and fibrosis. Lipoatrophic mice are heavier than controls, severely glucose intolerant, and hyperinsulinemic, and develop NAFLD characterized by increased liver glycogen, triacylglycerol, and diacylglycerol contents, mitotic index, apoptosis, inflammation, steatosis score, fibrosis, and fatty acid synthase (FAS) content and activity. Lipoatrophic mice also display liver enrichment with monounsaturated in detriment of polyunsaturated fatty acids including n-3 fatty acids, and increased content of cardiolipin, a tetracyl phospholipid exclusively found at the mitochondria inner membrane. Administration of a high-fat diet rich in n-3 fatty acids (HFO) to lipoatrophic mice enriches liver with n-3 fatty acids, reduces hepatic steatosis, FAS content and activity, apoptosis, inflammation, and improves glucose homeostasis.

CONCLUSION

Diet enrichment with n-3 fatty acids improves glucose homeostasis and reduces liver steatosis and inflammation without affecting hepatomegaly in severe lipoatrophic mice.

Integrative omics approaches provide biological and clinical insights: examples from mitochondrial diseases

High-throughput technologies for genomics, transcriptomics, proteomics, and metabolomics, and integrative analysis of these data, enable new, systems-level insights into disease pathogenesis. Mitochondrial diseases are an excellent target for hypothesis-generating omics approaches, as the disease group is mechanistically exceptionally complex. Although the genetic background in mitochondrial diseases is in either the nuclear or the mitochondrial genome, the typical downstream effect is dysfunction of the mitochondrial respiratory chain. However, the clinical manifestations show unprecedented variability, including either systemic or tissue-specific effects across multiple organ systems, with mild to severe symptoms, and occurring at any age. So far, the omics approaches have provided mechanistic understanding of tissue-specificity and potential treatment options for mitochondrial diseases, such as metabolome remodeling. However, no curative treatments exist, suggesting that novel approaches are needed. In this Review, we discuss omics approaches and discoveries with the potential to elucidate mechanisms of and therapies for mitochondrial diseases.

The Human Immunopeptidome Project: A Roadmap to Predict and Treat Immune Diseases

The science that investigates the ensembles of all peptides associated to human leukocyte antigen (HLA) molecules is termed "immunopeptidomics" and is typically driven by mass spectrometry (MS) technologies. Recent advances in MS technologies, neoantigen discovery and cancer immunotherapy have catalyzed the launch of the Human Immunopeptidome Project (HIPP) with the goal of providing a complete map of the human immunopeptidome and making the technology so robust that it will be available in every clinic. Here, we provide a long-term perspective of the field and we use this framework to explore how we think the completion of the HIPP will truly impact the society in the future. In this context, we introduce the concept of immunopeptidome-wide association studies (IWAS). We highlight the importance of large cohort studies for the future and how applying quantitative immunopeptidomics at population scale may provide a new look at individual predisposition to common immune diseases as well as responsiveness to vaccines and immunotherapies. Through this vision, we aim to provide a fresh view of the field to stimulate new discussions within the community, and present what we see as the key challenges for the future for unlocking the full potential of immunopeptidomics in this era of precision medicine.

A Cross-Species Systems Genetics Analysis Links APBB1IP as a Candidate for Schizophrenia and Prepulse Inhibition

Background: Prepulse inhibition (PPI) of the startle response is a highly conserved form of sensorimotor gating, disruption of which is found in schizophrenia patients and their unaffected first-degree relatives. PPI can be measured in many species, and shows considerable phenotypic variation between and within rodent models. This makes PPI a useful endophenotype. Genome-wide association studies (GWAS) have been carried out to identify genetic variants underlying schizophrenia, and these suggest that schizophrenia is highly polygenic. GWAS have been unable to account for the high heritability of schizophrenia seen in family studies, partly because of the low power of GWAS due to multiple comparisons. By contrast, complementary mouse model linkage studies often have high statistical power to detect variants for behavioral traits but lower resolution, producing loci that include tens or hundreds of genes. To capitalize on the advantages of both GWAS and genetic mouse models, our study uses a cross-species approach to identify novel genes associated with PPI regulation, which thus may contribute to the PPI deficits seen in schizophrenia.

Results: Using experimental data from the recombinant inbred (RI) mouse panel BXD, we identified two significant loci affecting PPI. These genomic regions contain genetic variants which influence PPI in mice and are therefore candidates that may be influencing aspects of schizophrenia in humans. We next investigated these regions in whole-genome data from the Psychiatric Genomics Consortium (PGC) schizophrenia GWAS and identify one novel candidate gene (ABPP1IP) that was significantly associated with PPI in mice and risk of schizophrenia in humans. A systems genetics approach demonstrates that APBB1IP coexpresses with several other genes related to schizophrenia in several brain regions. Gene coexpression and enrichment analysis shows clear links between APBB1IP and the immune system.

Conclusion: The combination of human GWAS and mouse quantitative trait loci (QTL) from some of the largest study systems available has enabled us to identify a novel gene, APBB1IP, which influences schizophrenia in humans and PPI in mice.

Hepatic gene expression explains primary drug toxicity in bipolar disorder

In bipolar disorder (BPD), long-term psychotropic drug treatment is often necessary to prevent relapse or recurrence. Nevertheless, adverse drug effects including disturbances in hepatic metabolism are observed and still poorly understood. Here, the association between hepatic gene expression and histopathological changes of the liver was investigated. By the use of microarrays (Affymetrix U133 plus2.0), a genome-wide expression study was performed on BPD patients with psychotropic drug treatment (n = 29) compared to unaffected controls (n = 20) and validated by quantitative real-time PCR. WebGestalt was used to identify over-represented functional pathways of the Reactome database. Association analyses between histopathological changes and differentially expressed genes comprised in the over-represented functional pathways were performed using regression analyses, from which feature-expression heatmaps were drawn. The majority of identified genes were underexpressed and involved in energy supply, metabolism of lipids and proteins, and the innate immune system. Positive associations were found for genes involved in all pathways and degenerative changes. The strongest negative association was observed between genes involved in energy supply and hepatic activity, as well as inflammation. In summary, we found a possible association between gene expression involved in various biological pathways and histopathological changes of the liver in BPD. Further, we found support for the probable primary toxic effect of psychotropic drugs on hepatic injury in BPD. Even if the safety of psychotropic drugs improves, adverse effects especially on hepatic function should not be underestimated.

Identifying gene function and module connections by the integration of multispecies expression compendia

The functions of many eukaryotic genes are still poorly understood. Here, we developed and validated a new method, termed GeneBridge, which is based on two linked approaches to impute gene function and bridge genes with biological processes. First, Gene-Module Association Determination (G-MAD) allows the annotation of gene function. Second, Module-Module Association Determination (M-MAD) allows predicting connectivity among modules. We applied the GeneBridge tools to large-scale multispecies expression compendia—1700 data sets with over 300,000 samples from human, mouse, rat, fly, worm, and yeast—collected in this study. G-MAD identifies novel functions of genes—for example, DDT in mitochondrial respiration and WDFY4 in T cell activation—and also suggests novel components for modules, such as for cholesterol biosynthesis. By applying G-MAD on data sets from respective tissues, tissue-specific functions of genes were identified—for instance, the roles of EHHADH in liver and kidney, as well as SLC6A1 in brain and liver. Using M-MAD, we identified a list of module-module associations, such as those between mitochondria and proteasome, mitochondria and histone demethylation, as well as ribosomes and lipid biosynthesis. The GeneBridge tools together with the expression compendia are available as an open resource, which will facilitate the identification of connections linking genes, modules, phenotypes, and diseases.

Cell stress management by the mitochondrial LonP1 protease - Insights into mitigating developmental, oncogenic and cardiac stress

Mitochondrial LonP1 is an essential stress response protease that mediates mitochondrial proteostasis, metabolism and bioenergetics. Homozygous and compound heterozygous variants in the LONP1 gene encoding the LonP1 protease have recently been shown to cause a diverse spectrum of human pathologies, ranging from classical mitochondrial disease phenotypes, profound neurologic impairment and multi-organ dysfunctions, some of which are uncommon to mitochondrial disorders. In this review, we focus primarily on human LonP1 and discuss findings, which demonstrate its multidimensional roles in maintaining mitochondrial proteostasis and adapting cells to metabolic flux and stress during normal physiology and disease processes. We also discuss emerging roles of LonP1 in responding to developmental, oncogenic and cardiac stress.

Proteomic, metabolic and immunological changes in Biomphalaria glabrata infected with Schistosoma mansoni

Mansonic schistosomiasis is a neglected disease transmitted by Biomphalaria spp. snails. Understanding what happens inside the intermediate host is important to develop more efficient ways of reducing schistosomiasis prevalence. Our purpose was to characterize metabolic and immunological changes in Biomphalaria glabrata 24 h after exposure to Schistosoma mansoni. For this purpose, proteins were extracted from snails' whole tissue with Tris-Urea buffer and digested with tripsin. Mass spectrometry was performed and analyzed with MaxQuant and Perseus software. Also, the hemolymph of five snails 24 h post exposure was collected, and the numbers of hemocytes, levels of urea, uric acid, nitric oxide, calcium, glycogen and alanine and aspartate aminotransferases activities were assessed. Snails were also dissected for measurement of glycogen content in the cephalopodal region and gonoda-digestive gland complex. Globin domain proteins were found to be up-regulated; also the number of circulating hemocytes was significantly higher after 24 h of exposure to the parasite. NO levels were higher 24 h post exposure. Several proteins associated with energy metabolism were found to be up-regulated. Glycogen analysis showed a significant decrease in the gonad-digestive gland complex glycogen content. We found several proteins which seem to be associated with the host immune response, most of which were up-regulated, however some were down-regulated, which may represent an important clue in understanding B. glabrata - S. mansoni compatibility.

Mitocellular communication: Shaping health and disease

Throughout the animal kingdom, mitochondria are the only organelles that retain their own genome and the transcription and translation machineries that are all essential for energy harvesting. Mitochondria have developed a complex communication network, allowing them to stay in tune with cellular needs and nuclear transcriptional programs and to alleviate mitochondrial dysfunction. Here, we review recent findings on the wide array of mechanisms that contribute to these mitocellular communication networks, spanning from well-studied messenger molecules to mitonuclear genetic interactions. Based on these observations and developments, we advocate a broad and inclusive view on mitocellular interactions, which can have profound impacts on physiological, pathological, and evolutionary processes.

Co-regulation map of the human proteome enables identification of protein functions

Assigning functions to the vast array of proteins present in eukaryotic cells remains challenging. To identify relationships between proteins, and thereby enable functional annotation of proteins, we determined changes in abundance of 10,323 human proteins in response to 294 biological perturbations using isotope-labeling mass spectrometry. We applied the machine learning algorithm treeClust to reveal functional associations between co-regulated human proteins from ProteomeHD, a compilation of our own data and datasets from the Proteomics Identifications database. This produced a co-regulation map of the human proteome. Co-regulation was able to capture relationships between proteins that do not physically interact or colocalize. For example, co-regulation of the peroxisomal membrane protein PEX11β with mitochondrial respiration factors led us to discover an organelle interface between peroxisomes and mitochondria in mammalian cells. We also predicted the functions of microproteins that are difficult to study with traditional methods. The co-regulation map can be explored at www.proteomeHD.net .

Sex-dependent and sex-independent regulatory systems of size variation in natural populations

Size of organs/organisms is a polygenic trait. Many of the growth-regulatory genes constitute conserved growth signaling pathways. However, how these multiple genes are orchestrated at the systems level to attain the natural variation in size including sexual size dimorphism is mostly unknown. Here we take a multi-layered systems omics approach to study size variation in the Drosophila wing. We show that expression levels of many critical growth regulators such as Wnt and TGFβ pathway components significantly differ between sexes but not between lines exhibiting size differences within each sex, suggesting a primary role of these regulators in sexual size dimorphism. Only a few growth genes including a receptor of steroid hormone ecdysone exhibit association with between-line size differences. In contrast, we find that between-line size variation is largely regulated by genes with a diverse range of cellular functions, most of which have never been implicated in growth. In addition, we show that expression quantitative trait loci (eQTLs) linked to these novel growth regulators accurately predict population-wide, between-line wing size variation. In summary, our study unveils differential gene regulatory systems that control wing size variation between and within sexes.

Diet modulates cecum bacterial diversity and physiological phenotypes across the BXD mouse genetic reference population

The BXD family has become one of the preeminent genetic reference populations to understand the genetic and environmental control of phenotypic variation. Here we evaluate the responses to different levels of fat in the diet using both chow diet (CD, 13–18% fat) and a high-fat diet (HFD, 45–60% fat). We studied cohorts of BXD strains, both inbred parents C57BL/6J and DBA/2J (commonly known as B6 and D2, respectively), as well as B6D2 and D2B6 reciprocal F1 hybrids. The comparative impact of genetic and dietary factors was analyzed by profiling a range of phenotypes, most prominently their cecum bacterial composition. The parents of the BXDs and F1 hybrids express limited differences in terms of weight and body fat gain on CD. In contrast, the strain differences on HFD are substantial for percent body fat, with DBA/2J accumulating 12.5% more fat than C57BL/6J (P < 0.0001). The F1 hybrids born to DBA/2J dams (D2B6F1) have 10.6% more body fat (P < 0.001) than those born to C57BL/6J dams. Sequence analysis of the cecum microbiota reveals important differences in bacterial composition among BXD family members with a substantial shift in composition caused by HFD. Relative to CD, the HFD induces a decline in diversity at the phylum level with a substantial increase in Firmicutes (+13.8%) and a reduction in Actinobacteria (-7.9%). In the majority of BXD strains, the HFD also increases cecal sIgA (P < 0.0001)—an important component of the adaptive immunity response against microbial pathogens. Host genetics modulates variation in cecum bacterial composition at the genus level in CD, with significant quantitative trait loci (QTLs) for Oscillibacter mapped to Chr 3 (18.7–19.2 Mb, LRS = 21.4) and for Bifidobacterium mapped to Chr 6 (89.21–89.37 Mb, LRS = 19.4). Introduction of HFD served as an environmental suppressor of these QTLs due to a reduction in the contribution of both genera (P < 0.001). Relations among liver metabolites and cecum bacterial composition were predominant in CD cohort, but these correlations do not persist following the shift to HFD. Overall, these findings demonstrate the important impact of environmental/dietary manipulation on the relationships between host genetics, gastrointestinal bacterial composition, immunological parameters, and metabolites—knowledge that will help in the understanding of the causal sources of metabolic disorders.

Broad Views of Non-alcoholic Fatty Liver Disease

Multi-omics multi-tissue data are used to interpret genome-wide association study results from mice to identify key driver genes of non-alcoholic fatty liver disease. Non-alcoholic fatty liver disease (NAFLD) is the accumulation of fat (steatosis) in the liver due to causes other than excessive alcohol consumption. The disease may progress to more severe forms of liver diseases, including non-alcoholic steatohepatitis, cirrhosis, and hepatocellular carcinoma. In this issue of Cell Systems, Krishnan et al. (2018) reveal mechanisms underlying NAFLD by generating multi-omics data using liver and adipose tissues obtained from the Hybrid Mouse Diversity Panel, consisting of 113 mouse strains with various degrees of NAFLD. The study identified key driver genes of NAFLD that can be used in the development of efficient treatment strategies and illustrates the potential utility of systematic analysis of multi-layer biological networks.

Metabolic and Stress Response Changes Precede Disease Onset in the Spinal Cord of Mutant SOD1 ALS Mice

Many Amyotrophic Lateral Sclerosis (ALS) patients experience hypermetabolism, or an increase in measured vs. calculated metabolic rate. The cause of hypermetabolism and the effects on neuronal metabolism in ALS are currently unknown, but the efficacy of dietary interventions shows promise for metabolism as an ALS therapeutic target. The goal of this study is to measure changes in metabolic pathways as a function of disease progression in spinal cords of the SOD1^G93A mouse model of ALS. We conducted a comprehensive assessment of protein expression for metabolic pathways, antioxidants, chaperones, and proteases in lumbar spinal cord from male SOD1^G93A mice at pre-onset, onset, and end-stages of the disease using targeted proteomic analysis. These results reveal that protein content of metabolic proteins including proteins involved in glycolysis, β-oxidation, and mitochondrial metabolism is altered in SOD1^G93A mouse spinal cord well before disease onset. The changes in mitochondrial metabolism proteins are associated with decreased maximal respiration and glycolytic flux in SOD1^G93A dermal fibroblasts and increased hydrogen peroxide and lipid hydroperoxide production in mitochondria from sciatic nerve and gastrocnemius muscle fibers at end stage of disease. Consistent with redox dysregulation, expression of the glutathione antioxidant system is decreased, and peroxiredoxins and catalase expression are increased. In addition, stress response proteases and chaperones, including those involved in the mitochondrial unfolded protein response (UPR^mt), are induced before disease onset. In summary, we report that metabolic and stress response changes occur in SOD1^G93A lumbar spinal cord before motor symptom onset, and are primarily caused by SOD1^G93A expression and do not vary greatly as a function of disease course.

Biology is the root of variability: cautionary tales in Caenorhabditis elegans biology

Reproducibility is critical for the standardization, interpretation, and progression of research. However, many factors increase variability and reduce reproducibility. In Caenorhabditis elegans research, there are many possible causes of variability that may explain why experimental outcomes sometimes differ between laboratories and between experiments. Factors contributing to experimental variability include the genetic background of both C. elegans and its bacterial diet, differences in media composition, intergenerational and transgenerational effects that may be carried over for generations, and the use of chemicals or reagents that may have unexpected consequences. This review summarizes sources of variability in C. elegans research and serves to identify laboratory practices that could influence reproducibility.

Metabolic and Innate Immune Cues Merge into a Specific Inflammatory Response via the UPR

Innate immune responses are intricately linked with intracellular metabolism of myeloid cells. Toll-like receptor (TLR) stimulation shifts intracellular metabolism toward glycolysis, while anti-inflammatory signals depend on enhanced mitochondrial respiration. How exogenous metabolic signals affect the immune response is unknown. We demonstrate that TLR-dependent responses of dendritic cells (DCs) are exacerbated by a high-fatty-acid (FA) metabolic environment. FAs suppress the TLR-induced hexokinase activity and perturb tricarboxylic acid cycle metabolism. These metabolic changes enhance mitochondrial reactive oxygen species (mtROS) production and, in turn, the unfolded protein response (UPR), leading to a distinct transcriptomic signature with IL-23 as hallmark. Interestingly, chemical or genetic suppression of glycolysis was sufficient to induce this specific immune response. Conversely, reducing mtROS production or DC-specific deficiency in XBP1 attenuated IL-23 expression and skin inflammation in an IL-23-dependent model of psoriasis. Thus, fine-tuning of innate immunity depends on optimization of metabolic demands and minimization of mtROS-induced UPR.

Mitochondrial stress-dependent regulation of cellular protein synthesis

The production of newly synthesized proteins is vital for all cellular functions and is a determinant of cell growth and proliferation. The synthesis of polypeptide chains from mRNA molecules requires sophisticated machineries and mechanisms that need to be tightly regulated, and adjustable to current needs of the cell. Failures in the regulation of translation contribute to the loss of protein homeostasis, which can have deleterious effects on cellular function and organismal health. Unsurprisingly, the regulation of translation appears to be a crucial element in stress response mechanisms. This review provides an overview of mechanisms that modulate cytosolic protein synthesis upon cellular stress, with a focus on the attenuation of translation in response to mitochondrial stress. We then highlight links between mitochondrion-derived reactive oxygen species and the attenuation of reversible cytosolic translation through the oxidation of ribosomal proteins at their cysteine residues. We also discuss emerging concepts of how cellular mechanisms to stress are adapted, including the existence of alternative ribosomes and stress granules, and the regulation of co-translational import upon organelle stress.

Disentangling Genetic and Environmental Effects on the Proteotypes of Individuals

Proteotypes, like genotypes, have been found to vary between individuals in several studies, but consistent molecular functional traits across studies remain to be quantified. In a meta-analysis of 11 proteomics datasets from humans and mice, we use co-variation of proteins in known functional modules across datasets and individuals to obtain a consensus landscape of proteotype variation. We find that individuals differ considerably in both protein complex abundances and stoichiometry. We disentangle genetic and environmental factors impacting these metrics, with genetic sex and specific diets together explaining 13.5% and 11.6% of the observed variation of complex abundance and stoichiometry, respectively. Sex-specific differences, for example, include various proteins and complexes, where the respective genes are not located on sex-specific chromosomes. Diet-specific differences, added to the individual genetic backgrounds, might become a starting point for personalized proteotype modulation toward desired features.

•

Benchmarking of datasets on human and mouse proteotypes

•

Consistent co-variation landscape of functional modules across individuals

•

Protein complexes vary in their stoichiometry across individuals

•

Quantifying effects of genetic sex and specific diets on complexes

Mild inborn errors of metabolism in commonly used inbred mouse strains

Inbred mouse strains are a cornerstone of translational research but paradoxically many strains carry mild inborn errors of metabolism. For example, α-aminoadipic acidemia and branched-chain ketoacid dehydrogenase deficiency are known in C57BL/6J mice. Using RNA sequencing, we now reveal the causal variants in Dhtkd1 and Bckdhb, and the molecular mechanism underlying these metabolic defects. C57BL/6J mice have decreased Dhtkd1 mRNA expression due to a solitary long terminal repeat (LTR) in intron 4 of Dhtkd1. This LTR harbors an alternate splice donor site leading to a partial splicing defect and as a consequence decreased total and functional Dhtkd1 mRNA, decreased DHTKD1 protein and α-aminoadipic acidemia. Similarly, C57BL/6J mice have decreased Bckdhb mRNA expression due to an LTR retrotransposon in intron 1 of Bckdhb. This transposable element encodes an alternative exon 1 causing aberrant splicing, decreased total and functional Bckdhb mRNA and decreased BCKDHB protein. Using a targeted metabolomics screen, we also reveal elevated plasma C5-carnitine in 129 substrains. This biochemical phenotype resembles isovaleric acidemia and is caused by an exonic splice mutation in Ivd leading to partial skipping of exon 10 and IVD protein deficiency. In summary, this study identifies three causal variants underlying mild inborn errors of metabolism in commonly used inbred mouse strains.

Evaluation of Sirtuin-3 probe quality and co-expressed genes using literature cohesion

BACKGROUND

Gene co-expression studies can provide important insights into molecular and cellular signaling pathways. The GeneNetwork database is a unique resource for co-expression analysis using data from a variety of tissues across genetically distinct inbred mice. However, extraction of biologically meaningful co-expressed gene sets is challenging due to variability in microarray platforms, probe quality, normalization methods, and confounding biological factors. In this study, we tested whether literature derived functional cohesion could be used as an objective metric in lieu of 'ground truth' to evaluate the quality of probes and microarray datasets.

RESULTS

We examined Sirtuin-3 (Sirt3) co-expressed gene sets extracted from either liver or brain tissues of BXD recombinant inbred mice in the GeneNetwork database. Depending on the microarray platform, there were as many as 26 probes that targeted different regions of Sirt3 primary transcript. Co-expressed gene sets (ranging from 100-1000 genes) associated with each Sirt3 probe were evaluated using the previously developed literature-derived cohesion p-value (LPv) and benchmarked against 'gold standards' derived from proteomic studies or Gene Ontology classifications. We found that the maximal F-measure was obtained at an average window size of 535 genes. Using set size of 500 genes, the Pearson correlations between LPv and F-measure as well as between LPv and mitochondrial gene enrichment p-values were 0.90 and 0.93, respectively. Importantly, we found that the LPv approach can distinguish high quality Sirt3 probes. Analysis of the most functionally cohesive Sirt3 co-expressed gene set revealed core metabolic pathways that were shared between hippocampus and liver as well as distinct pathways which were unique to each tissue. These results are consistent with other studies that suggest Sirt3 is a key metabolic regulator and has distinct functions in energy-producing vs. energy-demanding tissues.

CONCLUSIONS

Our results provide proof-of-concept that literature cohesion analysis is useful for evaluating the quality of probes and microarray datasets, particularly when experimentally derived gold standards are unavailable. Our approach would enable researchers to rapidly identify biologically meaningful co-expressed gene sets and facilitate discovery from high throughput genomic data.

An integrative systems genetic analysis of mammalian lipid metabolism

Dysregulation of lipid homeostasis is a precipitating event in the pathogenesis and progression of hepatosteatosis and metabolic syndrome. These conditions are highly prevalent in developed societies and currently have limited options for diagnostic and therapeutic intervention. Here, using a proteomic and lipidomic-wide systems genetic approach, we interrogated lipid regulatory networks in 107 genetically distinct mouse strains to reveal key insights into the control and network structure of mammalian lipid metabolism. These include the identification of plasma lipid signatures that predict pathological lipid abundance in the liver of mice and humans, defining subcellular localization and functionality of lipid-related proteins, and revealing functional protein and genetic variants that are predicted to modulate lipid abundance. Trans-omic analyses using these datasets facilitated the identification and validation of PSMD9 as a previously unknown lipid regulatory protein. Collectively, our study serves as a rich resource for probing mammalian lipid metabolism and provides opportunities for the discovery of therapeutic agents and biomarkers in the setting of hepatic lipotoxicity.

Systems Biology Approaches Toward Understanding Primary Mitochondrial Diseases

Primary mitochondrial diseases form one of the most common and severe groups of genetic disease, with a birth prevalence of at least 1 in 5000. These disorders are multi-genic and multi-phenotypic (even within the same gene defect) and span the entire age range from prenatal to late adult onset. Mitochondrial disease typically affects one or multiple high-energy demanding organs, and is frequently fatal in early life. Unfortunately, to date there are no known curative therapies, mostly owing to the rarity and heterogeneity of individual mitochondrial diseases, leading to diagnostic odysseys and difficulties in clinical trial design. This review aims to discuss recent advances and challenges of systems approaches for the study of primary mitochondrial diseases. Although there has been an explosion in the generation of omics data, few studies have progressed toward the integration of multiple levels of omics. It is evident that the integration of different types of data to create a more complete representation of biology remains challenging, perhaps due to the scarcity of available integrative tools and the complexity inherent in their use. In addition, "bottom-up" systems approaches have been adopted for use in the iterative cycle of systems biology: from data generation to model prediction and validation. Primary mitochondrial diseases, owing to their complex nature, will most likely benefit from a multidisciplinary approach encompassing clinical, molecular and computational studies integrated together by systems biology to elucidate underlying pathomechanisms for better diagnostics and therapeutic discovery. Just as next generation sequencing has rapidly increased diagnostic rates from approximately 5% up to 60% over two decades, more recent advancing technologies are encouraging; the generation of multi-omics, the integration of multiple types of data, and the ability to predict perturbations will, ultimately, be translated into improved patient care.

Characterization of Proteome Variation During Modern Maize Breeding

The success of modern maize breeding has been demonstrated by remarkable increases in productivity with tremendous modification of agricultural phenotypes over the last century. Although the underlying genetic changes of the maize adaptation from tropical to temperate regions have been extensively studied, our knowledge is limited regarding the accordance of protein and mRNA expression levels accompanying such adaptation. Here we conducted an integrative analysis of proteomic and transcriptomic changes in a maize association panel. The minimum extent of correlation between protein and RNA levels suggests that variation in mRNA expression is often not indicative of protein expression at a population scale. This is corroborated by the observation that mRNA- and protein-based coexpression networks are relatively independent of each other, and many pQTLs arise without the presence of corresponding eQTLs. Importantly, compared with transcriptome, the subtypes categorized by the proteome show a markedly high accuracy to resemble the genomic subpopulation. These findings suggest that proteome evolved under a greater evolutionary constraint than transcriptome during maize adaptation from tropical to temperate regions. Overall, the integrated multi-omics analysis provides a functional context to interpret gene expression variation during modern maize breeding.

Uncoupling of transcriptomic and cytological differentiation in mouse spermatocytes with impaired meiosis

Cell differentiation is driven by changes in gene expression that manifest as changes in cellular phenotype or function. Altered cellular phenotypes, stemming from genetic mutations or other perturbations, are widely assumed to directly correspond to changes in the transcriptome and vice versa. Here, we exploited the cytologically well-defined Prdm9 mutant mouse as a model of developmental arrest to test whether parallel programs of cellular differentiation and gene expression are tightly coordinated, or can be disassociated. By comparing cytological phenotype markers and transcriptomes in wild-type and mutant spermatocytes, we identified multiple instances of cellular and molecular uncoupling in Prdm9^–/– mutants. Most notably, although Prdm9^–/– germ cells undergo cytological arrest in a late-leptotene/zygotene stage, they nevertheless develop gene expression signatures characteristic of later developmental substages. These findings suggest that transcriptomic changes may not reliably map to cellular phenotypes in developmentally perturbed systems.

Mitochondrial ribosomal protein PTCD3 mutations cause oxidative phosphorylation defects with Leigh syndrome

Pentatricopeptide repeat domain proteins are a large family of RNA-binding proteins involved in mitochondrial RNA editing, stability, and translation. Mitochondrial translation machinery defects are an expanding group of genetic diseases in humans. We describe a patient who presented with low birth weight, mental retardation, and optic atrophy. Brain MRI showed abnormal bilateral signals at the basal ganglia and brainstem, and the patient was diagnosed as Leigh syndrome. Exome sequencing revealed two potentially loss-of-function variants [c.415-2A>G, and c.1747_1748insCT (p.Phe583Serfs*3)] in PTCD3 (also known as MRPS39). PTCD3, a member of the pentatricopeptide repeat domain protein family, is a component of the small mitoribosomal subunit. The patient had marked decreases in mitochondrial complex I and IV levels and activities, oxygen consumption and ATP biosynthesis, and generalized mitochondrial translation defects in fibroblasts. Quantitative proteomic analysis revealed decreased levels of the small mitoribosomal subunits. Complementation experiments rescued oxidative phosphorylation complex I and IV levels and activities, ATP biosynthesis, and MT-RNR1 rRNA transcript level, providing functional validation of the pathogenicity of identified variants. This is the first report of an association of PTCD3 mutations with Leigh syndrome along with combined oxidative phosphorylation deficiencies caused by defects in the mitochondrial translation machinery.

MicroRNA-382 silencing induces a mitonuclear protein imbalance and activates the mitochondrial unfolded protein response in muscle cells

Proper mitochondrial function plays a central role in cellular metabolism. Various diseases as well as aging are associated with diminished mitochondrial function. Previously, we identified 19 miRNAs putatively involved in the regulation of mitochondrial metabolism in skeletal muscle, a highly metabolically active tissue. In the current study, these 19 miRNAs were individually silenced in C2C12 myotubes using antisense oligonucleotides, followed by measurement of the expression of 27 genes known to play a major role in regulating mitochondrial metabolism. Based on the outcomes, we then focused on miR‐382‐5p and identified pathways affected by its silencing using microarrays, investigated protein expression, and studied cellular respiration. Silencing of miRNA‐382‐5p significantly increased the expression of several genes involved in mitochondrial dynamics and biogenesis. Conventional microarray analysis in C2C12 myotubes silenced for miRNA‐382‐5p revealed a collective downregulation of mitochondrial ribosomal proteins and respiratory chain proteins. This effect was accompanied by an imbalance between mitochondrial proteins encoded by the nuclear and mitochondrial DNA (1.35‐fold, p < 0.01) and an induction of HSP60 protein (1.31‐fold, p < 0.05), indicating activation of the mitochondrial unfolded protein response (mtUPR). Furthermore, silencing of miR‐382‐5p reduced basal oxygen consumption rate by 14% ( p < 0.05) without affecting mitochondrial content, pointing towards a more efficient mitochondrial function as a result of improved mitochondrial quality control. Taken together, silencing of miR‐382‐5p induces a mitonuclear protein imbalance and activates the mtUPR in skeletal muscle, a phenomenon that was previously associated with improved longevity.

Module-detection approaches for the integration of multilevel omics data highlight the comprehensive response of Aspergillus fumigatus to caspofungin

Background

Omics data provide deep insights into overall biological processes of organisms. However, integration of data from different molecular levels such as transcriptomics and proteomics, still remains challenging. Analyzing lists of differentially abundant molecules from diverse molecular levels often results in a small overlap mainly due to different regulatory mechanisms, temporal scales, and/or inherent properties of measurement methods. Module-detecting algorithms identifying sets of closely related proteins from protein-protein interaction networks (PPINs) are promising approaches for a better data integration.

Results

Here, we made use of transcriptome, proteome and secretome data from the human pathogenic fungus Aspergillus fumigatus challenged with the antifungal drug caspofungin. Caspofungin targets the fungal cell wall which leads to a compensatory stress response. We analyzed the omics data using two different approaches: First, we applied a simple, classical approach by comparing lists of differentially expressed genes (DEGs), differentially synthesized proteins (DSyPs) and differentially secreted proteins (DSePs); second, we used a recently published module-detecting approach, ModuleDiscoverer, to identify regulatory modules from PPINs in conjunction with the experimental data. Our results demonstrate that regulatory modules show a notably higher overlap between the different molecular levels and time points than the classical approach. The additional structural information provided by regulatory modules allows for topological analyses. As a result, we detected a significant association of omics data with distinct biological processes such as regulation of kinase activity, transport mechanisms or amino acid metabolism. We also found a previously unreported increased production of the secondary metabolite fumagillin by A. fumigatus upon exposure to caspofungin. Furthermore, a topology-based analysis of potential key factors contributing to drug-caused side effects identified the highly conserved protein polyubiquitin as a central regulator. Interestingly, polyubiquitin UbiD neither belonged to the groups of DEGs, DSyPs nor DSePs but most likely strongly influenced their levels.

Conclusion

Module-detecting approaches support the effective integration of multilevel omics data and provide a deep insight into complex biological relationships connecting these levels. They facilitate the identification of potential key players in the organism’s stress response which cannot be detected by commonly used approaches comparing lists of differentially abundant molecules.

Electronic supplementary material

The online version of this article (10.1186/s12918-018-0620-8) contains supplementary material, which is available to authorized users.

Determining Allele-Specific Protein Expression (ASPE) Using a Novel Quantitative Concatamer Based Proteomics Method

Measuring allele-specific expression (ASE) is a powerful approach for identifying cis-regulatory genetic variants. Here, we developed a novel targeted proteomics method for the quantification of allele-specific protein expression (ASPE) based on scheduled parallel reaction monitoring (PRM) with a heavy stable isotope-labeled quantitative concatamer (QconCAT) internal protein standard. This strategy was applied to the determination of the ASPE of UGT2B15 in human livers using the common UGT2B15 nonsynonymous variant rs1902023 (i.e., Y85D) as the marker to differentiate expressions from the two alleles. The QconCAT standard contains both the wild-type tryptic peptide and the Y85D mutant peptide at a ratio of 1:1 to ensure accurate measurement of the ASPE of UGT2B15. The results from 18 UGT2B15 Y85D heterozygotes revealed that the ratios between the wild-type Y allele and the mutant D allele varied from 0.60 to 1.46, indicating the presence of cis-regulatory variants. In addition, we observed no significant correlations between the ASPE and mRNA ASE of UGT2B15, suggesting the involvement of different cis-acting variants in regulating the transcription and translation processes of the gene. This novel ASPE approach provides a powerful tool for capturing cis-genetic variants involved in post-transcription processes, an important yet understudied area of research.

Skeletal muscle mitochondrial remodeling in exercise and diseases

Skeletal muscle fitness and plasticity is an important determinant of human health and disease. Mitochondria are essential for maintaining skeletal muscle energy homeostasis by adaptive re-programming to meet the demands imposed by a myriad of physiologic or pathophysiological stresses. Skeletal muscle mitochondrial dysfunction has been implicated in the pathogenesis of many diseases, including muscular dystrophy, atrophy, type 2 diabetes, and aging-related sarcopenia. Notably, exercise counteracts the effects of many chronic diseases on skeletal muscle mitochondrial function. Recent studies have revealed a finely tuned regulatory network that orchestrates skeletal muscle mitochondrial biogenesis and function in response to exercise and in disease states. In addition, increasing evidence suggests that mitochondria also serve to "communicate" with the nucleus and mediate adaptive genomic re-programming. Here we review the current state of knowledge relevant to the dynamic remodeling of skeletal muscle mitochondria in response to exercise and in disease states.

Inactivation of the E3 Ubiquitin Ligase IDOL Attenuates Diet-Induced Obesity and Metabolic Dysfunction in Mice

Objective- The E3 ubiquitin ligase IDOL (inducible degrader of the LDLR [LDL (low-density lipoprotein) receptor]) is a post-transcriptional regulator of LDLR abundance. Model systems and human genetics support a role for IDOL in regulating circulating LDL levels. Whether IDOL plays a broader metabolic role and affects development of metabolic syndrome-associated comorbidities is unknown. Approach and Results- We studied WT (wild type) and Idol(-/-) (Idol-KO) mice in 2 models: physiological aging and diet-induced obesity. In both models, deletion of Idol protected mice from metabolic dysfunction. On a Western-type diet, Idol loss resulted in decreased circulating levels of cholesterol, triglycerides, glucose, and insulin. This was accompanied by protection from weight gain in short- and long-term dietary challenges, which could be attributed to reduced hepatosteatosis and fat mass in Idol-KO mice. Although feeding and intestinal fat uptake were unchanged in Idol-KO mice, their brown adipose tissue was protected from lipid accumulation and had elevated expression of UCP1 (uncoupling protein 1) and TH (tyrosine hydroxylase). Indirect calorimetry indicated a marked increase in locomotion and suggested a trend toward increased cumulative energy expenditure and fat oxidation. An increase in in vivo clearance of reconstituted lipoprotein particles in Idol-KO mice may sustain this energetic demand. In the BXD mouse genetic reference population, hepatic Idol expression correlates with multiple metabolic parameters, thus providing support for findings in the Idol-KO mice. Conclusions- Our study uncovers an unrecognized role for Idol in regulation of whole body metabolism in physiological aging and on a Western-type diet. These findings support Idol inhibition as a therapeutic strategy to target multiple metabolic syndrome-associated comorbidities.

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.

GM-CSF driven myeloid cells in adipose tissue link weight gain and insulin resistance via formation of 2-aminoadipate

In a GM-CSF driven myeloid cell deficient mouse model (Csf2^−/−) that has preserved insulin sensitivity despite increased adiposity, we used unbiased three-dimensional integration of proteome profiles, metabolic profiles, and gene regulatory networks to understand adipose tissue proteome-wide changes and their metabolic implications. Multi-dimensional liquid chromatography mass spectrometry and extended multiplex mass labeling was used to analyze proteomes of epididymal adipose tissues isolated from Csf2^+/+ and Csf2^−/− mice that were fed low fat, high fat, or high fat plus cholesterol diets for 8 weeks. The metabolic health (as measured by body weight, adiposity, plasma fasting glucose, insulin, triglycerides, phospholipids, total cholesterol levels, and glucose and insulin tolerance tests) deteriorated with diet for both genotypes, while mice lacking Csf2 were protected from insulin resistance. Regardless of diet, 30 mostly mitochondrial, branch chain amino acids (BCAA), and lysine metabolism proteins were altered between Csf2^−/− and Csf2^+/+ mice (FDR < 0.05). Lack of GM-CSF driven myeloid cells lead to reduced adipose tissue 2-oxoglutarate dehydrogenase complex (DHTKD1) levels and subsequent increase in plasma 2-aminoadipate (2-AA) levels, both of which are reported to correlate with insulin resistance. Tissue DHTKD1 levels were >4-fold upregulated and plasma 2-AA levels were >2 fold reduced in Csf2^−/− mice (p < 0.05). GM-CSF driven myeloid cells link peripheral insulin sensitivity to adiposity via lysine metabolism involving DHTKD1/2-AA axis in a diet independent manner.