Minimum information and guidelines for reporting a multiplexed assay of variant effect.

Multiplexed assays of variant effect (MAVEs) have emerged as a powerful approach for interrogating thousands of genetic variants in a single experiment. The flexibility and widespread adoption of these techniques across diverse disciplines have led to a heterogeneous mix of data formats and descriptions, which complicates the downstream use of the resulting datasets. To address these issues and promote reproducibility and reuse of MAVE data, we define a set of minimum information standards for MAVE data and metadata and outline a controlled vocabulary aligned with established biomedical ontologies for describing these experimental designs.

Multi-ancestry polygenic mechanisms of type 2 diabetes.

Type 2 diabetes (T2D) is a multifactorial disease with substantial genetic risk, for which the underlying biological mechanisms are not fully understood. In this study, we identified multi-ancestry T2D genetic clusters by analyzing genetic data from diverse populations in 37 published T2D genome-wide association studies representing more than 1.4 million individuals. We implemented soft clustering with 650 T2D-associated genetic variants and 110 T2D-related traits, capturing known and novel T2D clusters with distinct cardiometabolic trait associations across two independent biobanks representing diverse genetic ancestral populations (African, n = 21,906; Admixed American, n = 14,410; East Asian, n =2,422; European, n = 90,093; and South Asian, n = 1,262). The 12 genetic clusters were enriched for specific single-cell regulatory regions. Several of the polygenic scores derived from the clusters differed in distribution among ancestry groups, including a significantly higher proportion of lipodystrophy-related polygenic risk in East Asian ancestry. T2D risk was equivalent at a body mass index (BMI) of 30 kg m-2 in the European subpopulation and 24.2 (22.9-25.5) kg m-2 in the East Asian subpopulation; after adjusting for cluster-specific genetic risk, the equivalent BMI threshold increased to 28.5 (27.1-30.0) kg m-2 in the East Asian group. Thus, these multi-ancestry T2D genetic clusters encompass a broader range of biological mechanisms and provide preliminary insights to explain ancestry-associated differences in T2D risk profiles.

Multi-ancestry Polygenic Mechanisms of Type 2 Diabetes Elucidate Disease Processes and Clinical Heterogeneity.

We identified genetic subtypes of type 2 diabetes (T2D) by analyzing genetic data from diverse groups, including non-European populations. We implemented soft clustering with 650 T2D-associated genetic variants, capturing known and novel T2D subtypes with distinct cardiometabolic trait associations. The twelve genetic clusters were distinctively enriched for single-cell regulatory regions. Polygenic scores derived from the clusters differed in distribution between ancestry groups, including a significantly higher proportion of lipodystrophy-related polygenic risk in East Asian ancestry. T2D risk was equivalent at a BMI of 30 kg/m2 in the European subpopulation and 24.2 (22.9-25.5) kg/m2 in the East Asian subpopulation; after adjusting for cluster-specific genetic risk, the equivalent BMI threshold increased to 28.5 (27.1-30.0) kg/m2 in the East Asian group, explaining about 75% of the difference in BMI thresholds. Thus, these multi-ancestry T2D genetic subtypes encompass a broader range of biological mechanisms and help explain ancestry-associated differences in T2D risk profiles.

Multi-ancestry Polygenic Mechanisms of Type 2 Diabetes Elucidate Disease Processes and Clinical Heterogeneity.

We identified genetic subtypes of type 2 diabetes (T2D) by analyzing genetic data from diverse groups, including non-European populations. We implemented soft clustering with 650 T2D-associated genetic variants, capturing known and novel T2D subtypes with distinct cardiometabolic trait associations. The twelve genetic clusters were distinctively enriched for single-cell regulatory regions. Polygenic scores derived from the clusters differed in distribution between ancestry groups, including a significantly higher proportion of lipodystrophy-related polygenic risk in East Asian ancestry. T2D risk was equivalent at a BMI of 30 kg/m2 in the European subpopulation and 24.2 (22.9-25.5) kg/m2 in the East Asian subpopulation; after adjusting for cluster-specific genetic risk, the equivalent BMI threshold increased to 28.5 (27.1-30.0) kg/m2 in the East Asian group, explaining about 75% of the difference in BMI thresholds. Thus, these multi-ancestry T2D genetic subtypes encompass a broader range of biological mechanisms and help explain ancestry-associated differences in T2D risk profiles.

A genome-wide atlas of human cell morphology.

A key challenge of the modern genomics era is developing data-driven representations of gene function. Here, we present the first unbiased morphology-based genome-wide perturbation atlas in human cells, containing three genome-scale genotype-phenotype maps comprising >20,000 single-gene CRISPR-Cas9-based knockout experiments in >30 million cells. Our optical pooled cell profiling approach (PERISCOPE) combines a de-stainable high-dimensional phenotyping panel (based on Cell Painting1,2) with optical sequencing of molecular barcodes and a scalable open-source analysis pipeline to facilitate massively parallel screening of pooled perturbation libraries. This approach provides high-dimensional phenotypic profiles of individual cells, while simultaneously enabling interrogation of subcellular processes. Our atlas reconstructs known pathways and protein-protein interaction networks, identifies culture media-specific responses to gene knockout, and clusters thousands of human genes by phenotypic similarity. Using this atlas, we identify the poorly-characterized disease-associated transmembrane protein TMEM251/LYSET as a Golgi-resident protein essential for mannose-6-phosphate-dependent trafficking of lysosomal enzymes, showing the power of these representations. In sum, our atlas and screening technology represent a rich and accessible resource for connecting genes to cellular functions at scale.

Minimum information and guidelines for reporting a Multiplexed Assay of Variant Effect.

Multiplexed Assays of Variant Effect (MAVEs) have emerged as a powerful approach for interrogating thousands of genetic variants in a single experiment. The flexibility and widespread adoption of these techniques across diverse disciplines has led to a heterogeneous mix of data formats and descriptions, which complicates the downstream use of the resulting datasets. To address these issues and promote reproducibility and reuse of MAVE data, we define a set of minimum information standards for MAVE data and metadata and outline a controlled vocabulary aligned with established biomedical ontologies for describing these experimental designs.

Discovering cellular programs of intrinsic and extrinsic drivers of metabolic traits using LipocyteProfiler.

A primary obstacle in translating genetic associations with disease into therapeutic strategies is elucidating the cellular programs affected by genetic risk variants and effector genes. Here, we introduce LipocyteProfiler, a cardiometabolic-disease-oriented high-content image-based profiling tool that enables evaluation of thousands of morphological and cellular profiles that can be systematically linked to genes and genetic variants relevant to cardiometabolic disease. We show that LipocyteProfiler allows surveillance of diverse cellular programs by generating rich context- and process-specific cellular profiles across hepatocyte and adipocyte cell-state transitions. We use LipocyteProfiler to identify known and novel cellular mechanisms altered by polygenic risk of metabolic disease, including insulin resistance, fat distribution, and the polygenic contribution to lipodystrophy. LipocyteProfiler paves the way for large-scale forward and reverse deep phenotypic profiling in lipocytes and provides a framework for the unbiased identification of causal relationships between genetic variants and cellular programs relevant to human disease.

A non-coding variant linked to metabolic obesity with normal weight affects actin remodelling in subcutaneous adipocytes.

Recent large-scale genomic association studies found evidence for a genetic link between increased risk of type 2 diabetes and decreased risk for adiposity-related traits, reminiscent of metabolically obese normal weight (MONW) association signatures. However, the target genes and cellular mechanisms driving such MONW associations remain to be identified. Here, we systematically identify the cellular programmes of one of the top-scoring MONW risk loci, the 2q24.3 risk locus, in subcutaneous adipocytes. We identify a causal genetic variant, rs6712203, an intronic single-nucleotide polymorphism in the COBLL1 gene, which changes the conserved transcription factor motif of POU domain, class 2, transcription factor 2, and leads to differential COBLL1 gene expression by altering the enhancer activity at the locus in subcutaneous adipocytes. We then establish the cellular programme under the genetic control of the 2q24.3 MONW risk locus and the effector gene COBLL1, which is characterized by impaired actin cytoskeleton remodelling in differentiating subcutaneous adipocytes and subsequent failure of these cells to accumulate lipids and develop into metabolically active and insulin-sensitive adipocytes. Finally, we show that perturbations of the effector gene Cobll1 in a mouse model result in organismal phenotypes matching the MONW association signature, including decreased subcutaneous body fat mass and body weight along with impaired glucose tolerance. Taken together, our results provide a mechanistic link between the genetic risk for insulin resistance and low adiposity, providing a potential therapeutic hypothesis and a framework for future identification of causal relationships between genome associations and cellular programmes in other disorders.

Engineered allele substitution at PPARGC1A rs8192678 alters human white adipocyte differentiation, lipogenesis, and PGC-1α content and turnover.

AIMS/HYPOTHESIS: PPARGC1A encodes peroxisome proliferator-activated receptor γ coactivator 1-α (PGC-1α), a central regulator of energy metabolism and mitochondrial function. A common polymorphism in PPARGC1A (rs8192678, C/T, Gly482Ser) has been associated with obesity and related metabolic disorders, but no published functional studies have investigated direct allele-specific effects in adipocyte biology. We examined whether rs8192678 is a causal variant and reveal its biological function in human white adipose cells. METHODS: We used CRISPR-Cas9 genome editing to perform an allelic switch (C-to-T or T-to-C) at rs8192678 in an isogenic human pre-adipocyte white adipose tissue (hWAs) cell line. Allele-edited single-cell clones were expanded and screened to obtain homozygous T/T (Ser482Ser), C/C (Gly482Gly) and heterozygous C/T (Gly482Ser) isogenic cell populations, followed by functional studies of the allele-dependent effects on white adipocyte differentiation and mitochondrial function. RESULTS: After differentiation, the C/C adipocytes were visibly less BODIPY-positive than T/T and C/T adipocytes, and had significantly lower triacylglycerol content. The C allele presented a dose-dependent lowering effect on lipogenesis, as well as lower expression of genes critical for adipogenesis, lipid catabolism, lipogenesis and lipolysis. Moreover, C/C adipocytes had decreased oxygen consumption rate (OCR) at basal and maximal respiration, and lower ATP-linked OCR. We determined that these effects were a consequence of a C-allele-driven dysregulation of PGC-1α protein content, turnover rate and transcriptional coactivator activity. CONCLUSIONS/INTERPRETATION: Our data show allele-specific causal effects of the rs8192678 variant on adipogenic differentiation. The C allele confers lower levels of PPARGC1A mRNA and PGC-1α protein, as well as disrupted dynamics of PGC-1α turnover and activity, with downstream effects on cellular differentiation and mitochondrial function. Our study provides the first experimentally deduced insights on the effects of rs8192678 on adipocyte phenotype.

Specific miRNAs are associated with human cancer cachexia in an organ-specific manner.

BACKGROUND: Cancer cachexia (CCx) is a complex and multi-organ wasting syndrome characterized by substantial weight loss and poor prognosis. An improved understanding of the mechanisms involved in the onset and progression of cancer cachexia is essential. How microRNAs contribute to the clinical manifestation and progression of CCx remains elusive. The aim of this study was to identify specific miRNAs related to organ-specific CCx and explore their functional role in humans. METHODS: miRNA patterns in serum and in cachexia target organs (liver, muscle and adipose tissue) from weight stable (N ≤ 12) and cachectic patients (N ≤ 23) with gastrointestinal cancer were analysed. As a first step, a miRNA array (158 miRNAs) was performed in pooled serum samples. Identified miRNAs were validated in serum and corresponding tissue samples. Using in silico prediction, related genes were identified and evaluated. The findings were confirmed in vitro by siRNA knock-down experiments in human visceral preadipocytes and C2C12 myoblast cells and consecutive gene expression analyses. RESULTS: Validating the results of the array, a 2-fold down-regulation of miR-122-5p (P = 0.0396) and a 4.5-fold down-regulation of miR-194-5p (P < 0.0001) in serum of CCx patients in comparison with healthy controls were detected. Only miR-122-5p correlated with weight loss and CCx status (P = 0.0367). Analysing corresponding tissues six muscle and eight visceral adipose tissue (VAT) cachexia-associated miRNAs were identified. miR-27b-3p, miR-375 and miR-424-5p were the most consistently affected miRNAs in tissues of CCx patients correlating negatively with the severity of body weight loss (P = 0.0386, P = 0.0112 and P = 0.0075, respectively). We identified numerous putative target genes of the miRNAs in association with muscle atrophy and lipolysis pathways. Knock-down experiments in C2C12 myoblast cells revealed an association of miR-27b-3p and the in silico predicted atrophy-related target genes IL-15 and TRIM63. Both were up-regulated in miR-27b-3p knock-down cells (P < 0.05). Concordantly, in muscle tissue of CCx individuals, significant higher expression levels of IL-15 (P = 0.0237) and TRIM63 (P = 0.0442) were detected. miR-424-5p was identified to regulate the expression of lipase genes. Knock-down experiments in human visceral preadipocytes revealed an inverse association of miR-424-5p with its predicted target genes LIPE, PNPLA2, MGLL and LPL (P < 0.01). CONCLUSIONS: The identified miRNAs, in particular miR-122-5p, miR-27b-3p, miR-375 and miR-424-5p, represent features of human CCx and may contribute to tissue wasting and skeletal muscle atrophy through the regulation of catabolic signals. Further studies are needed to explore the potential of the identified miRNAs as a screening tool for early detection of cancer cachexia.

FALCON systematically interrogates free fatty acid biology and identifies a novel mediator of lipotoxicity.

Cellular exposure to free fatty acids (FFAs) is implicated in the pathogenesis of obesity-associated diseases. However, there are no scalable approaches to comprehensively assess the diverse FFAs circulating in human plasma. Furthermore, assessing how FFA-mediated processes interact with genetic risk for disease remains elusive. Here, we report the design and implementation of fatty acid library for comprehensive ontologies (FALCON), an unbiased, scalable, and multimodal interrogation of 61 structurally diverse FFAs. We identified a subset of lipotoxic monounsaturated fatty acids associated with decreased membrane fluidity. Furthermore, we prioritized genes that reflect the combined effects of harmful FFA exposure and genetic risk for type 2 diabetes (T2D). We found that c-MAF-inducing protein (CMIP) protects cells from FFA exposure by modulating Akt signaling. In sum, FALCON empowers the study of fundamental FFA biology and offers an integrative approach to identify much needed targets for diverse diseases associated with disordered FFA metabolism.

FALCON systematically interrogates free fatty acid biology and identifies a novel mediator of lipotoxicity.

Cellular exposure to free fatty acids (FFA) is implicated in the pathogenesis of obesity-associated diseases. However, studies to date have assumed that a few select FFAs are representative of broad structural categories, and there are no scalable approaches to comprehensively assess the biological processes induced by exposure to diverse FFAs circulating in human plasma. Furthermore, assessing how these FFA- mediated processes interact with genetic risk for disease remains elusive. Here we report the design and implementation of FALCON (Fatty Acid Library for Comprehensive ONtologies) as an unbiased, scalable and multimodal interrogation of 61 structurally diverse FFAs. We identified a subset of lipotoxic monounsaturated fatty acids (MUFAs) with a distinct lipidomic profile associated with decreased membrane fluidity. Furthermore, we developed a new approach to prioritize genes that reflect the combined effects of exposure to harmful FFAs and genetic risk for type 2 diabetes (T2D). Importantly, we found that c-MAF inducing protein (CMIP) protects cells from exposure to FFAs by modulating Akt signaling and we validated the role of CMIP in human pancreatic beta cells. In sum, FALCON empowers the study of fundamental FFA biology and offers an integrative approach to identify much needed targets for diverse diseases associated with disordered FFA metabolism. Highlights: FALCON (Fatty Acid Library for Comprehensive ONtologies) enables multimodal profiling of 61 free fatty acids (FFAs) to reveal 5 FFA clusters with distinct biological effectsFALCON is applicable to many and diverse cell typesA subset of monounsaturated FAs (MUFAs) equally or more toxic than canonical lipotoxic saturated FAs (SFAs) leads to decreased membrane fluidityNew approach prioritizes genes that represent the combined effects of environmental (FFA) exposure and genetic risk for diseaseC-Maf inducing protein (CMIP) is identified as a suppressor of FFA-induced lipotoxicity via Akt-mediated signaling.

Identification of a weight loss-associated causal eQTL in MTIF3 and the effects of MTIF3 deficiency on human adipocyte function.

Genetic variation at the MTIF3 (Mitochondrial Translational Initiation Factor 3) locus has been robustly associated with obesity in humans, but the functional basis behind this association is not known. Here, we applied luciferase reporter assay to map potential functional variants in the haplotype block tagged by rs1885988 and used CRISPR-Cas9 to edit the potential functional variants to confirm the regulatory effects on MTIF3 expression. We further conducted functional studies on MTIF3-deficient differentiated human white adipocyte cell line (hWAs-iCas9), generated through inducible expression of CRISPR-Cas9 combined with delivery of synthetic MTIF3-targeting guide RNA. We demonstrate that rs67785913-centered DNA fragment (in LD with rs1885988, r2 > 0.8) enhances transcription in a luciferase reporter assay, and CRISPR-Cas9-edited rs67785913 CTCT cells show significantly higher MTIF3 expression than rs67785913 CT cells. Perturbed MTIF3 expression led to reduced mitochondrial respiration and endogenous fatty acid oxidation, as well as altered expression of mitochondrial DNA-encoded genes and proteins, and disturbed mitochondrial OXPHOS complex assembly. Furthermore, after glucose restriction, the MTIF3 knockout cells retained more triglycerides than control cells. This study demonstrates an adipocyte function-specific role of MTIF3, which originates in the maintenance of mitochondrial function, providing potential explanations for why MTIF3 genetic variation at rs67785913 is associated with body corpulence and response to weight loss interventions.

Impaired Adipocyte SLC7A10 Promotes Lipid Storage in Association With Insulin Resistance and Altered BCAA Metabolism.

CONTEXT: The neutral amino acid transporter SLC7A10/ASC-1 is an adipocyte-expressed gene with reduced expression in insulin resistance and obesity. Inhibition of SLC7A10 in adipocytes was shown to increase lipid accumulation despite decreasing insulin-stimulated uptake of glucose, a key substrate for de novo lipogenesis. These data imply that alternative lipogenic substrates to glucose fuel continued lipid accumulation during insulin resistance in obesity. OBJECTIVE: We examined whether increased lipid accumulation during insulin resistance in adipocytes may involve alter flux of lipogenic amino acids dependent on SLC7A10 expression and activity, and whether this is reflected by extracellular and circulating concentrations of marker metabolites. METHODS: In adipocyte cultures with impaired SLC7A10, we performed RNA sequencing and relevant functional assays. By targeted metabolite analyses (GC-MS/MS), flux of all amino acids and selected metabolites were measured in human and mouse adipose cultures. Additionally, SLC7A10 mRNA levels in human subcutaneous adipose tissue (SAT) were correlated to candidate metabolites and adiposity phenotypes in 2 independent cohorts. RESULTS: SLC7A10 impairment altered expression of genes related to metabolic processes, including branched-chain amino acid (BCAA) catabolism, lipogenesis, and glyceroneogenesis. In 3T3-L1 adipocytes, SLC7A10 inhibition increased fatty acid uptake and cellular content of glycerol and cholesterol. SLC7A10 impairment in SAT cultures altered uptake of aspartate and glutamate, and increased net uptake of BCAAs, while increasing the net release of the valine catabolite 3- hydroxyisobutyrate (3-HIB). In human cohorts, SLC7A10 mRNA correlated inversely with total fat mass, circulating triacylglycerols, BCAAs, and 3-HIB. CONCLUSION: Reduced SLC7A10 activity strongly affects flux of BCAAs in adipocytes, which may fuel continued lipogenesis during insulin resistance, and be reflected in increased circulating levels of the valine-derived catabolite 3-HIB.

Genetics of sexually dimorphic adipose distribution in humans.

Obesity-associated morbidity is exacerbated by abdominal obesity, which can be measured as the waist-to-hip ratio adjusted for the body mass index (WHRadjBMI). Here we identify genes associated with obesity and WHRadjBMI and characterize allele-sensitive enhancers that are predicted to regulate WHRadjBMI genes in women. We found that several waist-to-hip ratio-associated variants map within primate-specific Alu retrotransposons harboring a DNA motif associated with adipocyte differentiation. This suggests that a genetic component of adipose distribution in humans may involve co-option of retrotransposons as adipose enhancers. We evaluated the role of the strongest female WHRadjBMI-associated gene, SNX10, in adipose biology. We determined that it is required for human adipocyte differentiation and function and participates in diet-induced adipose expansion in female mice, but not males. Our data identify genes and regulatory mechanisms that underlie female-specific adipose distribution and mediate metabolic dysfunction in women.

BMI-adjusted adipose tissue volumes exhibit depot-specific and divergent associations with cardiometabolic diseases.

For any given body mass index (BMI), individuals vary substantially in fat distribution, and this variation may have important implications for cardiometabolic risk. Here, we study disease associations with BMI-independent variation in visceral (VAT), abdominal subcutaneous (ASAT), and gluteofemoral (GFAT) fat depots in 40,032 individuals of the UK Biobank with body MRI. We apply deep learning models based on two-dimensional body MRI projections to enable near-perfect estimation of fat depot volumes (R2 in heldout dataset = 0.978-0.991 for VAT, ASAT, and GFAT). Next, we derive BMI-adjusted metrics for each fat depot (e.g. VAT adjusted for BMI, VATadjBMI) to quantify local adiposity burden. VATadjBMI is associated with increased risk of type 2 diabetes and coronary artery disease, ASATadjBMI is largely neutral, and GFATadjBMI is associated with reduced risk. These results - describing three metabolically distinct fat depots at scale - clarify the cardiometabolic impact of BMI-independent differences in body fat distribution.

Scalable Functional Assays for the Interpretation of Human Genetic Variation.

Scalable sequence-function studies have enabled the systematic analysis and cataloging of hundreds of thousands of coding and noncoding genetic variants in the human genome. This has improved clinical variant interpretation and provided insights into the molecular, biophysical, and cellular effects of genetic variants at an astonishing scale and resolution across the spectrum of allele frequencies. In this review, we explore current applications and prospects for the field and outline the principles underlying scalable functional assay design, with a focus on the study of single-nucleotide coding and noncoding variants.

miR-375 is cold exposure sensitive and drives thermogenesis in visceral adipose tissue derived stem cells.

Activation of brown adipose tissue may increase energy expenditure by non-shivering thermogenesis. Cold exposure is one of the options to activate brown adipocytes. To link changes in energy metabolism with microRNA expression (miRNAs), we analyzed 158 miRNAs in serum of 169 healthy individuals before and after cold exposure. Validating the results of a miRNA array, a significant down-regulation of miR-375 after cold exposure (P < 0.0001) was detected. These changes went along with a significant negative correlation between miR-375 and visceral adipose tissue (VAT) mass (P < 0.0001), implicating a specific function of miR-375 in this depot. Significantly higher expression levels of miR-375 were found in VAT in comparison to subcutaneous fat (SAT). Using in silico prediction, we identified putative miR-375 target genes involved in the thermogenesis pathway. Cold-stimulation of subcutaneous and visceral pre-adipocytes (PACs) led to significantly higher expression levels of FABP4, FGF21, PPARGC1A and PRDM16 in VC-PACs. Analyzing miR-375 knock down and cold stimulated VC-PACs revealed a significant up-regulation of thermogenesis associated genes PPARGC1A, ELOVL3 and PRDM16. In summary, our findings identified miR-375 as a potential adipogenic and thermogenesis-associated miRNA exclusively acting in visceral adipose tissue.