Lipidomic QTL in Diversity Outbred mice identifies a novel function for α/ß hydrolase domain 2 (Abhd2) as an enzyme that metabolizes phosphatidylcholine and cardiolipin.

We and others have previously shown that genetic association can be used to make causal connections between gene loci and small molecules measured by mass spectrometry in the bloodstream and in tissues. We identified a locus on mouse chromosome 7 where several phospholipids in liver showed strong genetic association to distinct gene loci. In this study, we integrated gene expression data with genetic association data to identify a single gene at the chromosome 7 locus as the driver of the phospholipid phenotypes. The gene encodes α/ß-hydrolase domain 2 (Abhd2), one of 23 members of the ABHD gene family. We validated this observation by measuring lipids in a mouse with a whole-body deletion of Abhd2. The Abhd2KO mice had a significant increase in liver levels of phosphatidylcholine and phosphatidylethanolamine. Unexpectedly, we also found a decrease in two key mitochondrial lipids, cardiolipin and phosphatidylglycerol, in male Abhd2KO mice. These data suggest that Abhd2 plays a role in the synthesis, turnover, or remodeling of liver phospholipids.

Cistrome Explorer: an interactive visual analysis tool for large-scale epigenomic data.

SUMMARY: The regulation of genes by cis-regulatory elements (CREs) is complex and differs between cell types. Visual analysis of large collections of chromatin profiles across diverse cell types, integrated with computational methods, can reveal meaningful biological insights. We developed Cistrome Explorer, a web-based interactive visual analytics tool for exploring thousands of chromatin profiles in diverse cell types. Integrated with the Cistrome Data Browser database which contains thousands of ChIP-seq, DNase-seq and ATAC-seq samples, Cistrome Explorer enables the discovery of patterns of CREs across cell types and the identification of transcription factor binding underlying these patterns. AVAILABILITY AND IMPLEMENTATION: Cistrome Explorer and its source code are available at http://cisvis.gehlenborglab.org/ and released under the MIT License. Documentation can be accessed via http://cisvis.gehlenborglab.org/docs/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Diet-Microbiota Interactions Mediate Global Epigenetic Programming in Multiple Host Tissues.

Histone-modifying enzymes regulate transcription and are sensitive to availability of endogenous small-molecule metabolites, allowing chromatin to respond to changes in environment. The gut microbiota produces a myriad of metabolites that affect host physiology and susceptibility to disease; however, the underlying molecular events remain largely unknown. Here we demonstrate that microbial colonization regulates global histone acetylation and methylation in multiple host tissues in a diet-dependent manner: consumption of a "Western-type" diet prevents many of the microbiota-dependent chromatin changes that occur in a polysaccharide-rich diet. Finally, we demonstrate that supplementation of germ-free mice with short-chain fatty acids, major products of gut bacterial fermentation, is sufficient to recapitulate chromatin modification states and transcriptional responses associated with colonization. These findings have profound implications for understanding the complex functional interactions between diet, gut microbiota, and host health.

Identification of genetic drivers of plasma lipoprotein size in the Diversity Outbred mouse population.

Despite great progress in understanding lipoprotein physiology, there is still much to be learned about the genetic drivers of lipoprotein abundance, composition, and function. We used ion mobility spectrometry to survey 16 plasma lipoprotein subfractions in 500 Diversity Outbred mice maintained on a Western-style diet. We identified 21 quantitative trait loci (QTL) affecting lipoprotein abundance. To refine the QTL and link them to disease risk in humans, we asked if the human homologs of genes located at each QTL were associated with lipid traits in human genome-wide association studies. Integration of mouse QTL with human genome-wide association studies yielded candidate gene drivers for 18 of the 21 QTL. This approach enabled us to nominate the gene encoding the neutral ceramidase, Asah2, as a novel candidate driver at a QTL on chromosome 19 for large HDL particles (HDL-2b). To experimentally validate Asah2, we surveyed lipoproteins in Asah2-/- mice. Compared to wild-type mice, female Asah2-/- mice showed an increase in several lipoproteins, including HDL. Our results provide insights into the genetic regulation of circulating lipoproteins, as well as mechanisms by which lipoprotein subfractions may affect cardiovascular disease risk in humans.

Missense variants in SORT1 are associated with LDL-C in an Amish population.

Common noncoding variants at the human 1p13.3 locus associated with SORT1 expression are among those most strongly associated with low-density lipoprotein cholesterol (LDL-C) in human genome-wide association studies. However, validation studies in mice and cell lines have produced variable results regarding the directionality of the effect of SORT1 on LDL-C. This, together with the fact that the 1p13.3 variants are associated with expression of several genes, has raised the question of whether SORT1 is the causal gene at this locus. Using whole exome sequencing in members of an Amish population, we identified coding variants in SORT1 that are associated with increased (rs141749679, K302E) and decreased (rs149456022, Q225H) LDL-C. Further, analysis of plasma lipoprotein particle subclasses by ion mobility in a subset of rs141749679 (K302E) carriers revealed higher levels of large LDL particles compared to noncarriers. In contrast to the effect of these variants in the Amish, the sortilin K302E mutation introduced into a C57BL/6J mouse via CRISPR/Cas9 resulted in decreased non-high-density lipoprotein cholesterol, and the sortilin Q225H mutation did not alter cholesterol levels in mice. This is indicative of different effects of these mutations on cholesterol metabolism in the two species. To our knowledge, this is the first evidence that naturally occurring coding variants in SORT1 are associated with LDL-C, thus supporting SORT1 as the gene responsible for the association of the 1p13.3 locus with LDL-C.

OME-Zarr: a cloud-optimized bioimaging file format with international community support.

A growing community is constructing a next-generation file format (NGFF) for bioimaging to overcome problems of scalability and heterogeneity. Organized by the Open Microscopy Environment (OME), individuals and institutes across diverse modalities facing these problems have designed a format specification process (OME-NGFF) to address these needs. This paper brings together a wide range of those community members to describe the cloud-optimized format itself-OME-Zarr-along with tools and data resources available today to increase FAIR access and remove barriers in the scientific process. The current momentum offers an opportunity to unify a key component of the bioimaging domain-the file format that underlies so many personal, institutional, and global data management and analysis tasks.

Extreme Domain Wall Speeds under Ultrafast Optical Excitation.

Time-resolved ultrafast EUV magnetic scattering was used to test a recent prediction of >10 km/s domain wall speeds by optically exciting a magnetic sample with a nanoscale labyrinthine domain pattern. Ultrafast distortion of the diffraction pattern was observed at markedly different timescales compared to the magnetization quenching. The diffraction pattern distortion shows a threshold dependence with laser fluence, not seen for magnetization quenching, consistent with a picture of domain wall motion with pinning sites. Supported by simulations, we show that a speed of ≈66 km/s for highly curved domain walls can explain the experimental data. While our data agree with the prediction of extreme, nonequilibrium wall speeds locally, it differs from the details of the theory, suggesting that additional mechanisms are required to fully understand these effects.

Sorting through the extensive and confusing roles of sortilin in metabolic disease.

Sortilin is a post-Golgi trafficking receptor homologous to the yeast vacuolar protein sorting receptor 10 (VPS10). The VPS10 motif on sortilin is a 10-bladed ß-propeller structure capable of binding more than 50 proteins, covering a wide range of biological functions including lipid and lipoprotein metabolism, neuronal growth and death, inflammation, and lysosomal degradation. Sortilin has a complex cellular trafficking itinerary, where it functions as a receptor in the trans-Golgi network, endosomes, secretory vesicles, multivesicular bodies, and at the cell surface. In addition, sortilin is associated with hypercholesterolemia, Alzheimer's disease, prion diseases, Parkinson's disease, and inflammation syndromes. The 1p13.3 locus containing SORT1, the gene encoding sortilin, carries the strongest association with LDL-C of all loci in human genome-wide association studies. However, the mechanism by which sortilin influences LDL-C is unclear. Here, we review the role sortilin plays in cardiovascular and metabolic diseases and describe in detail the large and often contradictory literature on the role of sortilin in the regulation of LDL-C levels.

MetaNetwork Enhances Biological Insights from Quantitative Proteomics Differences by Combining Clustering and Enrichment Analyses.

Interpreting proteomics data remains challenging due to the large number of proteins that are quantified by modern mass spectrometry methods. Weighted gene correlation network analysis (WGCNA) can identify groups of biologically related proteins using only protein intensity values by constructing protein correlation networks. However, WGCNA is not widespread in proteomic analyses due to challenges in implementing workflows. To facilitate the adoption of WGCNA by the proteomics field, we created MetaNetwork, an open-source, R-based application to perform sophisticated WGCNA workflows with no coding skill requirements for the end user. We demonstrate MetaNetwork's utility by employing it to identify groups of proteins associated with prostate cancer from a proteomic analysis of tumor and adjacent normal tissue samples. We found a decrease in cytoskeleton-related protein expression, a known hallmark of prostate tumors. We further identified changes in module eigenproteins indicative of dysregulation in protein translation and trafficking pathways. These results demonstrate the value of using MetaNetwork to improve the biological interpretation of quantitative proteomics experiments with 15 or more samples.

SurvMaximin: Robust federated approach to transporting survival risk prediction models.

OBJECTIVE: For multi-center heterogeneous Real-World Data (RWD) with time-to-event outcomes and high-dimensional features, we propose the SurvMaximin algorithm to estimate Cox model feature coefficients for a target population by borrowing summary information from a set of health care centers without sharing patient-level information. MATERIALS AND METHODS: For each of the centers from which we want to borrow information to improve the prediction performance for the target population, a penalized Cox model is fitted to estimate feature coefficients for the center. Using estimated feature coefficients and the covariance matrix of the target population, we then obtain a SurvMaximin estimated set of feature coefficients for the target population. The target population can be an entire cohort comprised of all centers, corresponding to federated learning, or a single center, corresponding to transfer learning. RESULTS: Simulation studies and a real-world international electronic health records application study, with 15 participating health care centers across three countries (France, Germany, and the U.S.), show that the proposed SurvMaximin algorithm achieves comparable or higher accuracy compared with the estimator using only the information of the target site and other existing methods. The SurvMaximin estimator is robust to variations in sample sizes and estimated feature coefficients between centers, which amounts to significantly improved estimates for target sites with fewer observations. CONCLUSIONS: The SurvMaximin method is well suited for both federated and transfer learning in the high-dimensional survival analysis setting. SurvMaximin only requires a one-time summary information exchange from participating centers. Estimated regression vectors can be very heterogeneous. SurvMaximin provides robust Cox feature coefficient estimates without outcome information in the target population and is privacy-preserving.

Genetic determinants of gut microbiota composition and bile acid profiles in mice.

The microbial communities that inhabit the distal gut of humans and other mammals exhibit large inter-individual variation. While host genetics is a known factor that influences gut microbiota composition, the mechanisms underlying this variation remain largely unknown. Bile acids (BAs) are hormones that are produced by the host and chemically modified by gut bacteria. BAs serve as environmental cues and nutrients to microbes, but they can also have antibacterial effects. We hypothesized that host genetic variation in BA metabolism and homeostasis influence gut microbiota composition. To address this, we used the Diversity Outbred (DO) stock, a population of genetically distinct mice derived from eight founder strains. We characterized the fecal microbiota composition and plasma and cecal BA profiles from 400 DO mice maintained on a high-fat high-sucrose diet for ~22 weeks. Using quantitative trait locus (QTL) analysis, we identified several genomic regions associated with variations in both bacterial and BA profiles. Notably, we found overlapping QTL for Turicibacter sp. and plasma cholic acid, which mapped to a locus containing the gene for the ileal bile acid transporter, Slc10a2. Mediation analysis and subsequent follow-up validation experiments suggest that differences in Slc10a2 gene expression associated with the different strains influences levels of both traits and revealed novel interactions between Turicibacter and BAs. This work illustrates how systems genetics can be utilized to generate testable hypotheses and provide insight into host-microbe interactions.

Reversal of hypertriglyceridemia in diabetic BTBR ob/ob mice does not prevent nephropathy.

The etiology of diabetic nephropathy in type 2 diabetes is multifactorial. Sustained hyperglycemia is a major contributor, but additional contributions come from the hypertension, obesity, and hyperlipidemia that are also commonly present in patients with type 2 diabetes and nephropathy. The leptin deficient BTBR ob/ob mouse is a model of type 2 diabetic nephropathy in which hyperglycemia, obesity, and hyperlipidemia, but not hypertension, are present. We have shown that reversal of the constellation of these metabolic abnormalities with leptin replacement can reverse the morphologic and functional manifestations of diabetic nephropathy. Here we tested the hypothesis that reversal specifically of the hypertriglyceridemia, using an antisense oligonucleotide directed against ApoC-III, an apolipoprotein that regulates the interactions of VLDL (very low density lipoproteins) with the LDL receptor, is sufficient to ameliorate the nephropathy of Type 2 diabetes. Antisense treatment resulted in reduction of circulating ApoC-III protein levels and resulted in substantial lowering of triglycerides to near-normal levels in diabetic mice versus controls. Antisense treatment did not ameliorate proteinuria or pathologic manifestations of diabetic nephropathy, including podocyte loss. These studies indicate that pathologic manifestations of diabetic nephropathy are unlikely to be reduced by lipid-lowering therapeutics alone, but does not preclude a role for such interventions to be used in conjunction with other therapeutics commonly employed in the treatment of diabetes and its complications.

Distinguishing attosecond electron-electron scattering and screening in transition metals.

Electron-electron interactions are the fastest processes in materials, occurring on femtosecond to attosecond timescales, depending on the electronic band structure of the material and the excitation energy. Such interactions can play a dominant role in light-induced processes such as nano-enhanced plasmonics and catalysis, light harvesting, or phase transitions. However, to date it has not been possible to experimentally distinguish fundamental electron interactions such as scattering and screening. Here, we use sequences of attosecond pulses to directly measure electron-electron interactions in different bands of different materials with both simple and complex Fermi surfaces. By extracting the time delays associated with photoemission we show that the lifetime of photoelectrons from the d band of Cu are longer by ∼100 as compared with those from the same band of Ni. We attribute this to the enhanced electron-electron scattering in the unfilled d band of Ni. Using theoretical modeling, we can extract the contributions of electron-electron scattering and screening in different bands of different materials with both simple and complex Fermi surfaces. Our results also show that screening influences high-energy photoelectrons (≈20 eV) significantly less than low-energy photoelectrons. As a result, high-energy photoelectrons can serve as a direct probe of spin-dependent electron-electron scattering by neglecting screening. This can then be applied to quantifying the contribution of electron interactions and screening to low-energy excitations near the Fermi level. The information derived here provides valuable and unique information for a host of quantum materials.

Islet proteomics reveals genetic variation in dopamine production resulting in altered insulin secretion.

The mouse is a critical model in diabetes research, but most research in mice has been limited to a small number of mouse strains and limited genetic variation. Using the eight founder strains and both sexes of the Collaborative Cross (C57BL/6J (B6), A/J, 129S1/SvImJ (129), NOD/ShiLtJ (NOD), NZO/HILtJ (NZO), PWK/PhJ (PWK), WSB/EiJ (WSB), and CAST/EiJ (CAST)), we investigated the genetic dependence of diabetes-related metabolic phenotypes and insulin secretion. We found that strain background is associated with an extraordinary range in body weight, plasma glucose, insulin, triglycerides, and insulin secretion. Our whole-islet proteomic analysis of the eight mouse strains demonstrates that genetic background exerts a strong influence on the islet proteome that can be linked to the differences in diabetes-related metabolic phenotypes and insulin secretion. We computed protein modules consisting of highly correlated proteins that enrich for biological pathways and provide a searchable database of the islet protein expression profiles. To validate the data resource, we identified tyrosine hydroxylase (Th), a key enzyme in catecholamine synthesis, as a protein that is highly expressed in ß-cells of PWK and CAST islets. We show that CAST islets synthesize elevated levels of dopamine, which suppresses insulin secretion. Prior studies, using only the B6 strain, concluded that adult mouse islets do not synthesize l-3,4-dihydroxyphenylalanine (l-DOPA), the product of Th and precursor of dopamine. Thus, the choice of the CAST strain, guided by our islet proteomic survey, was crucial for these discoveries. In summary, we provide a valuable data resource to the research community, and show that proteomic analysis identified a strain-specific pathway by which dopamine synthesized in ß-cells inhibits insulin secretion.

Exploiting Prophage-Mediated Lysis for Biotherapeutic Release by Lactobacillus reuteri.

Lactobacillus reuteri has the potential to be developed as a microbial therapeutic delivery platform because of an established safety profile, health-promoting properties, and available genome editing tools. Here, we show that L. reuteri VPL1014 exhibits a low mutation rate compared to other Gram-positive bacteria, which we expect will contribute to the stability of genetically modified strains. VPL1014 encodes two biologically active prophages, which are induced during gastrointestinal transit. We hypothesized that intracellularly accumulated recombinant protein can be released following bacteriophage-mediated lysis. To test this, we engineered VPL1014 to accumulate leptin, our model protein, inside the cell. In vitro prophage induction of recombinant VPL1014 released leptin into the extracellular milieu, which corresponded to bacteriophage production. We also employed a plasmid system that does not require antibiotic in the growth medium for plasmid maintenance. Collectively, these data provide new avenues to exploit native prophages to deliver therapeutic molecules.IMPORTANCE Lactic acid bacteria (LAB) have been explored as potential biotherapeutic vehicles for the past 20 years. To secrete a therapeutic in the extracellular milieu, one typically relies on the bacterial secretion pathway, i.e., the Sec pathway. Overexpression of a secreted protein can overload the secretory pathway and impact the organism's fitness, and optimization of the signal peptide is also required to maximize the efficiency of the release of mature protein. Here, we describe a previously unexplored approach to release therapeutics from the probiotic Lactobacillus reuteri We demonstrate that an intracellularly accumulated recombinant protein is released following prophage activation. Since we recently demonstrated that prophages are activated during gastrointestinal transit, we propose that this method will provide a straightforward and efficient approach to deliver therapeutics in vivo.

The Transcription Factor Nfatc2 Regulates ß-Cell Proliferation and Genes Associated with Type 2 Diabetes in Mouse and Human Islets.

Human genome-wide association studies (GWAS) have shown that genetic variation at >130 gene loci is associated with type 2 diabetes (T2D). We asked if the expression of the candidate T2D-associated genes within these loci is regulated by a common locus in pancreatic islets. Using an obese F2 mouse intercross segregating for T2D, we show that the expression of ~40% of the T2D-associated genes is linked to a broad region on mouse chromosome (Chr) 2. As all but 9 of these genes are not physically located on Chr 2, linkage to Chr 2 suggests a genomic factor(s) located on Chr 2 regulates their expression in trans. The transcription factor Nfatc2 is physically located on Chr 2 and its expression demonstrates cis linkage; i.e., its expression maps to itself. When conditioned on the expression of Nfatc2, linkage for the T2D-associated genes was greatly diminished, supporting Nfatc2 as a driver of their expression. Plasma insulin also showed linkage to the same broad region on Chr 2. Overexpression of a constitutively active (ca) form of Nfatc2 induced ß-cell proliferation in mouse and human islets, and transcriptionally regulated more than half of the T2D-associated genes. Overexpression of either ca-Nfatc2 or ca-Nfatc1 in mouse islets enhanced insulin secretion, whereas only ca-Nfatc2 was able to promote ß-cell proliferation, suggesting distinct molecular pathways mediating insulin secretion vs. ß-cell proliferation are regulated by NFAT. Our results suggest that many of the T2D-associated genes are downstream transcriptional targets of NFAT, and may act coordinately in a pathway through which NFAT regulates ß-cell proliferation in both mouse and human islets.

Perilipin 5 and liver fatty acid binding protein function to restore quiescence in mouse hepatic stellate cells.

Hepatic stellate cell (HSC) activation occurs along with decreased Perilipin5 (Plin5) and liver fatty acid-binding protein (L-Fabp) expression and coincident lipid droplet (LD) depletion. Conversely, the activated phenotype is reversible in WT HSCs upon forced expression of Plin5. Here, we asked if L-Fabp expression is required for Plin5-mediated rescue of the quiescent phenotype. Lentiviral Plin5 transduction of passaged L-Fabp-/- HSCs failed to reverse activation markers or restore lipogenic gene expression and LD formation. However, adenoviral L-Fabp infection of lentiviral Plin5 transduced L-Fabp-/- HSCs restored both the quiescent phenotype and LD formation, an effect also mediated by adenoviral intestine-Fabp or adipocyte-Fabp. Expression of exogenous Plin5 in activated WT HSCs induced a transcriptional program of lipogenic gene expression including endogenous L-Fabp, but none of the other FABPs. We further demonstrated that selective, small molecule inhibition of endogenous L-Fabp also eliminated the ability of exogenous Plin5 to rescue LD formation and reverse activation of WT HSCs. This functional coordination of L-Fabp with Plin5 was 5'-AMP-activated protein kinase (AMPK)-dependent and was eliminated by AMPK inhibition. Taken together, our results indicate that L-Fabp is required for Plin5 to activate a transcriptional program that restores LD formation and reverses HSC activation.

Revealing the Nature of the Ultrafast Magnetic Phase Transition in Ni by Correlating Extreme Ultraviolet Magneto-Optic and Photoemission Spectroscopies.

By correlating time- and angle-resolved photoemission and time-resolved transverse magneto-optical Kerr effect measurements, both at extreme ultraviolet wavelengths, we uncover the universal nature of the ultrafast photoinduced magnetic phase transition in Ni. This allows us to explain the ultrafast magnetic response of Ni at all laser fluences-from a small reduction of the magnetization at low laser fluences, to complete quenching at high laser fluences. Both probe methods exhibit the same demagnetization and recovery timescales. The spin system absorbs the energy required to proceed through a magnetic phase transition within 20 fs after the peak of the pump pulse. However, the spectroscopic signatures of demagnetization of the material appear only after ≈200 fs and the subsequent recovery of magnetization on timescales ranging from 500 fs to >70 ps. We also provide evidence of two competing channels with two distinct timescales in the recovery process that suggest the presence of coexisting phases in the material.

Targeted Mass Spectrometry Approach Enabled Discovery of O-Glycosylated Insulin and Related Signaling Peptides in Mouse and Human Pancreatic Islets.

O-Linked glycosylation often involves the covalent attachment of sugar moieties to the hydroxyl group of serine or threonine on proteins/peptides. Despite growing interest in glycoproteins, little attention has been directed to glycosylated signaling peptides, largely due to lack of enabling analytical tools. Here we explore the occurrence of naturally O-linked glycosylation on the signaling peptides extracted from mouse and human pancreatic islets using mass spectrometry (MS). A novel targeted MS-based method is developed to increase the likelihood of capturing these modified signaling peptides and to provide improved sequence coverage and accurate glycosite localization, enabling the first large-scale discovery of O-glycosylation on signaling peptides. Several glycosylated signaling peptides with multiple glycoforms are identified, including the first report of glycosylated insulin-B chain and insulin-C peptide and BigLEN. This discovery may reveal potential novel functions as glycosylation could influence their conformation and biostability. Given the importance of insulin and its related peptide hormones and previous studies of glycosylated insulin analogues, this natural glycosylation may provide important insights into diabetes research and therapeutic treatments.