Prospective functional classification of all possible missense variants in PPARG

Sequencing methods, functional characterization, prevalence, and penetrance of rare coding mutations in panels of monogenic obesity genes from the leptin-melanocortin pathway: A systematic review

The recent development of next-generation sequencing (NGS) technologies has led to an increase of mutation screening reports of monogenic obesity genes in diverse experimental designs. However, no study to date has summarized their findings. Two reviewers independently conducted a systematic review of MEDLINE, Embase, and Web of Science Core Collection databases from inception to September 2022 to identify monogenic non-syndromic obesity gene screening studies. Of 1051 identified references, 31 were eligible after title and abstract screening and 28 after full-text reading and risk of bias and quality assessment. Most studies (82%) used NGS methods. The number of genes screened varied from 2 to 12 genes from the leptin-melanocortin pathway. While all the included studies used in silico tools to assess the functional status of mutations, only 2 performed in vitro tests. The prevalence of carriers of pathogenic/likely pathogenic monogenic mutations is 13.24% on average (heterozygous: 12.31%; homozygous/heterozygous composite: 0.93%). As no study reported the penetrance of pathogenic mutations on obesity, we estimated that homozygous carriers exhibited a complete penetrance (100%) and heterozygous carriers a variable penetrance (3-100%). The review provides an exhaustive description of sequencing methods, functional characterization, prevalence, and penetrance of rare coding mutations in monogenic non-syndromic obesity genes.

Integrated multiplexed assays of variant effect reveal determinants of catechol-O-methyltransferase gene expression

Multiplexed assays of variant effect are powerful methods to profile the consequences of rare variants on gene expression and organismal fitness. Yet, few studies have integrated several multiplexed assays to map variant effects on gene expression in coding sequences. Here, we pioneered a multiplexed assay based on polysome profiling to measure variant effects on translation at scale, uncovering single-nucleotide variants that increase or decrease ribosome load. By combining high-throughput ribosome load data with multiplexed mRNA and protein abundance readouts, we mapped the cis-regulatory landscape of thousands of catechol-O-methyltransferase (COMT) variants from RNA to protein and found numerous coding variants that alter COMT expression. Finally, we trained machine learning models to map signatures of variant effects on COMT gene expression and uncovered both directional and divergent impacts across expression layers. Our analyses reveal expression phenotypes for thousands of variants in COMT and highlight variant effects on both single and multiple layers of expression. Our findings prompt future studies that integrate several multiplexed assays for the readout of gene expression.

Rare variant associations with birth weight identify genes involved in adipose tissue regulation, placental function and insulin-like growth factor signalling

Investigating the genetic factors influencing human birth weight may lead to biological insights into fetal growth and long-term health. Genome-wide association studies of birth weight have highlighted associated variants in more than 200 regions of the genome, but the causal genes are mostly unknown. Rare genetic variants with robust evidence of association are more likely to point to causal genes, but to date, only a few rare variants are known to influence birth weight. We aimed to identify genes that harbour rare variants that impact birth weight when carried by either the fetus or the mother, by analysing whole exome sequence data in UK Biobank participants. We annotated rare (minor allele frequency <0.1%) protein-truncating or high impact missense variants on whole exome sequence data in up to 234,675 participants with data on their own birth weight (fetal variants), and up to 181,883 mothers who reported the birth weight of their first child (maternal variants). Variants within each gene were collapsed to perform gene burden tests and for each associated gene, we compared the observed fetal and maternal effects. We identified 8 genes with evidence of rare fetal variant effects on birth weight, of which 2 also showed maternal effects. One additional gene showed evidence of maternal effects only. We observed 10/11 directionally concordant associations in an independent sample of up to 45,622 individuals (sign test P=0.01). Of the genes identified, IGF1R and PAPPA2 (fetal and maternal-acting) have known roles in insulin-like growth factor bioavailability and signalling. PPARG, INHBE and ACVR1C (all fetal-acting) have known roles in adipose tissue regulation and rare variants in the latter two also showed associations with favourable adiposity patterns in adults. We highlight the dual role of PPARG in both adipocyte differentiation and placental angiogenesis. NOS3, NRK, and ADAMTS8 (fetal and maternal-acting) have been implicated in both placental function and hypertension. Analysis of rare coding variants has identified regulators of fetal adipose tissue and fetoplacental angiogenesis as determinants of birth weight, as well as further evidence for the role of insulin-like growth factors.

Rheostats, toggles, and neutrals, Oh my! A new framework for understanding how amino acid changes modulate protein function

Advances in personalized medicine and protein engineering require accurately predicting outcomes of amino acid substitutions. Many algorithms correctly predict that evolutionarily-conserved positions show "toggle" substitution phenotypes, which is defined when a few substitutions at that position retain function. In contrast, predictions often fail for substitutions at the less-studied "rheostat" positions, which are defined when different amino acid substitutions at a position sample at least half of the possible functional range. This review describes efforts to understand the impact and significance of rheostat positions: (1) They have been observed in globular soluble, integral membrane, and intrinsically disordered proteins; within single proteins, their prevalence can be up to 40%. (2) Substitutions at rheostat positions can have biological consequences and ∼10% of substitutions gain function. (3) Although both rheostat and "neutral" (defined when all substitutions exhibit wild-type function) positions are nonconserved, the two classes have different evolutionary signatures. (4) Some rheostat positions have pleiotropic effects on function, simultaneously modulating multiple parameters (e.g., altering both affinity and allosteric coupling). (5) In structural studies, substitutions at rheostat positions appear to cause only local perturbations; the overall conformations appear unchanged. (6) Measured functional changes show promising correlations with predicted changes in protein dynamics; the emergent properties of predicted, dynamically coupled amino acid networks might explain some of the complex functional outcomes observed when substituting rheostat positions. Overall, rheostat positions provide unique opportunities for using single substitutions to tune protein function. Future studies of these positions will yield important insights into the protein sequence/function relationship.

The intrinsically disordered transcriptional activation domain of CIITA is functionally tuneable by single substitutions: An exception or a new paradigm?

During protein evolution, some amino acid substitutions modulate protein function ("tuneability"). In most proteins, the tuneable range is wide and can be sampled by a set of protein variants that each contains multiple amino acid substitutions. In other proteins, the full tuneable range can be accessed by a set of variants that each contains a single substitution. Indeed, in some globular proteins, the full tuneable range can be accessed by the set of site-saturating substitutions at an individual "rheostat" position. However, in proteins with intrinsically disordered regions (IDRs), most functional studies-which would also detect tuneability-used multiple substitutions or small deletions. In disordered transcriptional activation domains (ADs), studies with multiple substitutions led to the "acidic exposure" model, which does not anticipate the existence of rheostat positions. In the few studies that did assess effects of single substitutions on AD function, results were mixed: the ADs of two full-length transcription factors did not show tuneability, whereas a fragment of a third AD was tuneable by single substitutions. In this study, we tested tuneability in the AD of full-length human class II transactivator (CIITA). Sequence analyses and experiments showed that CIITA's AD is an IDR. Functional assays of singly-substituted AD variants showed that CIITA's function was highly tuneable, with outcomes not predicted by the acidic exposure model. Four tested positions showed rheostat behavior for transcriptional activation. Thus, tuneability of different IDRs can vary widely. Future studies are needed to illuminate the biophysical features that govern whether an IDR is tuneable by single substitutions.

SIGMA leverages protein structural information to predict the pathogenicity of missense variants

Leveraging protein structural information to evaluate pathogenicity has been hindered by the scarcity of experimentally determined 3D protein. With the aid of AlphaFold2 predictions, we developed the structure-informed genetic missense mutation assessor (SIGMA) to predict missense variant pathogenicity. In comparison with existing predictors across labeled variant datasets and experimental datasets, SIGMA demonstrates superior performance in predicting missense variant pathogenicity (AUC = 0.933). We found that the relative solvent accessibility of the mutated residue contributed greatly to the predictive ability of SIGMA. We further explored combining SIGMA with other top-tier predictors to create SIGMA+, proving highly effective for variant pathogenicity prediction (AUC = 0.966). To facilitate the application of SIGMA, we pre-computed SIGMA scores for over 48 million possible missense variants across 3,454 disease-associated genes and developed an interactive online platform (https://www.sigma-pred.org/). Overall, by leveraging protein structure information, SIGMA offers an accurate structure-based approach to evaluating the pathogenicity of missense variants.

ProteinGym: Large-Scale Benchmarks for Protein Design and Fitness Prediction

Predicting the effects of mutations in proteins is critical to many applications, from understanding genetic disease to designing novel proteins that can address our most pressing challenges in climate, agriculture and healthcare. Despite a surge in machine learning-based protein models to tackle these questions, an assessment of their respective benefits is challenging due to the use of distinct, often contrived, experimental datasets, and the variable performance of models across different protein families. Addressing these challenges requires scale. To that end we introduce ProteinGym, a large-scale and holistic set of benchmarks specifically designed for protein fitness prediction and design. It encompasses both a broad collection of over 250 standardized deep mutational scanning assays, spanning millions of mutated sequences, as well as curated clinical datasets providing high-quality expert annotations about mutation effects. We devise a robust evaluation framework that combines metrics for both fitness prediction and design, factors in known limitations of the underlying experimental methods, and covers both zero-shot and supervised settings. We report the performance of a diverse set of over 70 high-performing models from various subfields (eg., alignment-based, inverse folding) into a unified benchmark suite. We open source the corresponding codebase, datasets, MSAs, structures, model predictions and develop a user-friendly website that facilitates data access and analysis.

Deep Mutational Scanning in Disease-related Genes with Saturation Mutagenesis-Reinforced Functional Assays (SMuRF)

Interpretation of disease-causing genetic variants remains a challenge in human genetics. Current costs and complexity of deep mutational scanning methods hamper crowd-sourcing approaches toward genome-wide resolution of variants in disease-related genes. Our framework, Saturation Mutagenesis-Reinforced Functional assays (SMuRF), addresses these issues by offering simple and cost-effective saturation mutagenesis, as well as streamlining functional assays to enhance the interpretation of unresolved variants. Applying SMuRF to neuromuscular disease genes FKRP and LARGE1, we generated functional scores for over 99.8% of all possible coding single nucleotide variants and resolved 310 clinically reported variants of uncertain significance with high confidence, enhancing clinical variant interpretation in dystroglycanopathies. SMuRF also demonstrates utility in predicting disease severity, resolving critical structural regions, and providing training datasets for the development of computational predictors. Our approach opens new directions for enabling variant-to-function insights for disease genes in a manner that is broadly useful for crowd-sourcing implementation across standard research laboratories.

MmisAT and MmisP: an efficient and accurate suite of variant analysis toolkit for primary mitochondrial diseases

Recent advances in next-generation sequencing (NGS) technology have greatly accelerated the need for efficient annotation to accurately interpret clinically relevant genetic variants in human diseases. Therefore, it is crucial to develop appropriate analytical tools to improve the interpretation of disease variants. Given the unique genetic characteristics of mitochondria, including haplogroup, heteroplasmy, and maternal inheritance, we developed a suite of variant analysis toolkits specifically designed for primary mitochondrial diseases: the Mitochondrial Missense Variant Annotation Tool (MmisAT) and the Mitochondrial Missense Variant Pathogenicity Predictor (MmisP). MmisAT can handle protein-coding variants from both nuclear DNA and mtDNA and generate 349 annotation types across six categories. It processes 4.78 million variant data in 76 min, making it a valuable resource for clinical and research applications. Additionally, MmisP provides pathogenicity scores to predict the pathogenicity of genetic variations in mitochondrial disease. It has been validated using cross-validation and external datasets and demonstrated higher overall discriminant accuracy with a receiver operating characteristic (ROC) curve area under the curve (AUC) of 0.94, outperforming existing pathogenicity predictors. In conclusion, the MmisAT is an efficient tool that greatly facilitates the process of variant annotation, expanding the scope of variant annotation information. Furthermore, the development of MmisP provides valuable insights into the creation of disease-specific, phenotype-specific, and even gene-specific predictors of pathogenicity, further advancing our understanding of specific fields.

Integrated multiplexed assays of variant effect reveal cis-regulatory determinants of catechol-O-methyltransferase gene expression

Multiplexed assays of variant effect are powerful methods to profile the consequences of rare variants on gene expression and organismal fitness. Yet, few studies have integrated several multiplexed assays to map variant effects on gene expression in coding sequences. Here, we pioneered a multiplexed assay based on polysome profiling to measure variant effects on translation at scale, uncovering single-nucleotide variants that increase and decrease ribosome load. By combining high-throughput ribosome load data with multiplexed mRNA and protein abundance readouts, we mapped the cis-regulatory landscape of thousands of catechol-O-methyltransferase (COMT) variants from RNA to protein and found numerous coding variants that alter COMT expression. Finally, we trained machine learning models to map signatures of variant effects on COMT gene expression and uncovered both directional and divergent impacts across expression layers. Our analyses reveal expression phenotypes for thousands of variants in COMT and highlight variant effects on both single and multiple layers of expression. Our findings prompt future studies that integrate several multiplexed assays for the readout of gene expression.

Genotype-stratified treatment for monogenic insulin resistance: a systematic review

BACKGROUND

Monogenic insulin resistance (IR) includes lipodystrophy and disorders of insulin signalling. We sought to assess the effects of interventions in monogenic IR, stratified by genetic aetiology.

METHODS

Systematic review using PubMed, MEDLINE and Embase (1 January 1987 to 23 June 2021). Studies reporting individual-level effects of pharmacologic and/or surgical interventions in monogenic IR were eligible. Individual data were extracted and duplicates were removed. Outcomes were analysed for each gene and intervention, and in aggregate for partial, generalised and all lipodystrophy.

RESULTS

10 non-randomised experimental studies, 8 case series, and 23 case reports meet inclusion criteria, all rated as having moderate or serious risk of bias. Metreleptin use is associated with the lowering of triglycerides and haemoglobin A1c (HbA1c) in all lipodystrophy (n = 111), partial (n = 71) and generalised lipodystrophy (n = 41), and in LMNA, PPARG, AGPAT2 or BSCL2 subgroups (n = 72,13,21 and 21 respectively). Body Mass Index (BMI) is lowered in partial and generalised lipodystrophy, and in LMNA or BSCL2, but not PPARG or AGPAT2 subgroups. Thiazolidinediones are associated with improved HbA1c and triglycerides in all lipodystrophy (n = 13), improved HbA1c in PPARG (n = 5), and improved triglycerides in LMNA (n = 7). In INSR-related IR, rhIGF-1, alone or with IGFBP3, is associated with improved HbA1c (n = 17). The small size or absence of other genotype-treatment combinations preclude firm conclusions.

CONCLUSIONS

The evidence guiding genotype-specific treatment of monogenic IR is of low to very low quality. Metreleptin and Thiazolidinediones appear to improve metabolic markers in lipodystrophy, and rhIGF-1 appears to lower HbA1c in INSR-related IR. For other interventions, there is insufficient evidence to assess efficacy and risks in aggregated lipodystrophy or genetic subgroups.

The use of precision diagnostics for monogenic diabetes: a systematic review and expert opinion

BACKGROUND

Monogenic diabetes presents opportunities for precision medicine but is underdiagnosed. This review systematically assessed the evidence for (1) clinical criteria and (2) methods for genetic testing for monogenic diabetes, summarized resources for (3) considering a gene or (4) variant as causal for monogenic diabetes, provided expert recommendations for (5) reporting of results; and reviewed (6) next steps after monogenic diabetes diagnosis and (7) challenges in precision medicine field.

METHODS

Pubmed and Embase databases were searched (1990-2022) using inclusion/exclusion criteria for studies that sequenced one or more monogenic diabetes genes in at least 100 probands (Question 1), evaluated a non-obsolete genetic testing method to diagnose monogenic diabetes (Question 2). The risk of bias was assessed using the revised QUADAS-2 tool. Existing guidelines were summarized for questions 3-5, and review of studies for questions 6-7, supplemented by expert recommendations. Results were summarized in tables and informed recommendations for clinical practice.

RESULTS

There are 100, 32, 36, and 14 studies included for questions 1, 2, 6, and 7 respectively. On this basis, four recommendations for who to test and five on how to test for monogenic diabetes are provided. Existing guidelines for variant curation and gene-disease validity curation are summarized. Reporting by gene names is recommended as an alternative to the term MODY. Key steps after making a genetic diagnosis and major gaps in our current knowledge are highlighted.

CONCLUSIONS

We provide a synthesis of current evidence and expert opinion on how to use precision diagnostics to identify individuals with monogenic diabetes.

Accurate proteome-wide missense variant effect prediction with AlphaMissense

The vast majority of missense variants observed in the human genome are of unknown clinical significance. We present AlphaMissense, an adaptation of AlphaFold fine-tuned on human and primate variant population frequency databases to predict missense variant pathogenicity. By combining structural context and evolutionary conservation, our model achieves state-of-the-art results across a wide range of genetic and experimental benchmarks, all without explicitly training on such data. The average pathogenicity score of genes is also predictive for their cell essentiality, capable of identifying short essential genes that existing statistical approaches are underpowered to detect. As a resource to the community, we provide a database of predictions for all possible human single amino acid substitutions and classify 89% of missense variants as either likely benign or likely pathogenic.

PFAS Exposures and the Human Metabolome: A Systematic Review of Epidemiological Studies

Purpose of Review

There is a growing interest in understanding the health effects of exposure to per- and polyfluoroalkyl substances (PFAS) through the study of the human metabolome. In this systematic review, we aimed to identify consistent findings between PFAS and metabolomic signatures. We conducted a search matching specific keywords that was independently reviewed by two authors on two databases (EMBASE and PubMed) from their inception through July 19, 2022 following PRISMA guidelines.

Recent Findings

We identified a total of 28 eligible observational studies that evaluated the associations between 31 different PFAS exposures and metabolomics in humans. The most common exposure evaluated was legacy long-chain PFAS. Population sample sizes ranged from 40 to 1,105 participants at different stages across the lifespan. A total of 19 studies used a non-targeted metabolomics approach, 7 used targeted approaches, and 2 included both. The majority of studies were cross-sectional (n = 25), including four with prospective analyses of PFAS measured prior to metabolomics.

Summary

Most frequently reported associations across studies were observed between PFAS and amino acids, fatty acids, glycerophospholipids, glycerolipids, phosphosphingolipids, bile acids, ceramides, purines, and acylcarnitines. Corresponding metabolic pathways were also altered, including lipid, amino acid, carbohydrate, nucleotide, energy metabolism, glycan biosynthesis and metabolism, and metabolism of cofactors and vitamins. We found consistent evidence across studies indicating PFAS-induced alterations in lipid and amino acid metabolites, which may be involved in energy and cell membrane disruption.

Analysis of the genetic contribution to thoracic aortic aneurysm or dissection in a prospective cohort of patients with familial and sporadic cases in East China

BACKGROUND

Thoracic aortic aneurysm or dissections (TAADs) represent a group of life-threatening diseases. Genetic aetiology can affect the age of onset, clinical phenotype, and timing of intervention. We conducted a prospective trial to determine the prevalence of pathogenic variants in TAAD patients and to elucidate the traits related to harbouring the pathogenic variants. One hundred and one unrelated TAAD patients underwent genetic sequencing and analysis for 23 TAAD-associated genes using a targeted PCR and next-generation sequencing-based panel.

RESULTS

A total of 47 variants were identified in 52 TAAD patients (51.5%), including 5 pathogenic, 1 likely pathogenic and 41 variants of uncertain significance. The pathogenic or likely pathogenic (P/LP) variants in 4 disease-causing genes were carried by 1 patient with familial and 5 patients with sporadic TAAD (5.9%). In addition to harbouring one variant causing familial TAAD, the FBN1 gene harboured half of the P/LP variants causing sporadic TAAD. Individuals with an age of onset less than 50 years or normotension had a significantly increased genetic risk.

CONCLUSIONS

TAAD patients with a younger age at diagnosis or normotension were more likely to carry a P/LP variant; thus, routine genetic testing will be beneficial to a better prognosis through genetically personalized care prior to acute rupture or dissection.

Human gain-of-function variants in HNF1A confer protection from diabetes but independently increase hepatic secretion of atherogenic lipoproteins

Loss-of-function mutations in hepatocyte nuclear factor 1A (HNF1A) are known to cause rare forms of diabetes and alter hepatic physiology through unclear mechanisms. In the general population, 1:100 individuals carry a rare, protein-coding HNF1A variant, most of unknown functional consequence. To characterize the full allelic series, we performed deep mutational scanning of 11,970 protein-coding HNF1A variants in human hepatocytes and clinical correlation with 553,246 exome-sequenced individuals. Surprisingly, we found that ∼1:5 rare protein-coding HNF1A variants in the general population cause molecular gain of function (GOF), increasing the transcriptional activity of HNF1A by up to 50% and conferring protection from type 2 diabetes (odds ratio [OR] = 0.77, p = 0.007). Increased hepatic expression of HNF1A promoted a pro-atherogenic serum profile mediated in part by enhanced transcription of risk genes including ANGPTL3 and PCSK9. In summary, ∼1:300 individuals carry a GOF variant in HNF1A that protects carriers from diabetes but enhances hepatic secretion of atherogenic lipoproteins.

A systematic review on machine learning approaches in the diagnosis and prognosis of rare genetic diseases

BACKGROUND

The diagnosis of rare genetic diseases is often challenging due to the complexity of the genetic underpinnings of these conditions and the limited availability of diagnostic tools. Machine learning (ML) algorithms have the potential to improve the accuracy and speed of diagnosis by analyzing large amounts of genomic data and identifying complex multiallelic patterns that may be associated with specific diseases. In this systematic review, we aimed to identify the methodological trends and the ML application areas in rare genetic diseases.

METHODS

We performed a systematic review of the literature following the PRISMA guidelines to search studies that used ML approaches to enhance the diagnosis of rare genetic diseases. Studies that used DNA-based sequencing data and a variety of ML algorithms were included, summarized, and analyzed using bibliometric methods, visualization tools, and a feature co-occurrence analysis.

FINDINGS

Our search identified 22 studies that met the inclusion criteria. We found that exome sequencing was the most frequently used sequencing technology (59%), and rare neoplastic diseases were the most prevalent disease scenario (59%). In rare neoplasms, the most frequent applications of ML models were the differential diagnosis or stratification of patients (38.5%) and the identification of somatic mutations (30.8%). In other rare diseases, the most frequent goals were the prioritization of rare variants or genes (55.5%) and the identification of biallelic or digenic inheritance (33.3%). The most employed method was the random forest algorithm (54.5%). In addition, the features of the datasets needed for training these algorithms were distinctive depending on the goal pursued, including the mutational load in each gene for the differential diagnosis of patients, or the combination of genotype features and sequence-derived features (such as GC-content) for the identification of somatic mutations.

CONCLUSIONS

ML algorithms based on sequencing data are mainly used for the diagnosis of rare neoplastic diseases, with random forest being the most common approach. We identified key features in the datasets used for training these ML models according to the objective pursued. These features can support the development of future ML models in the diagnosis of rare genetic diseases.

Quantitative determination of SLC2A1 variant functional effects in GLUT1 deficiency syndrome

OBJECTIVE

The goal of this study is to demonstrate the utility of a growth assay to quantify the functional impact of single nucleotide variants (SNVs) in SLC2A1, the gene responsible for Glut1DS.

METHODS

The functional impact of 40 SNVs in SLC2A1 was quantitatively determined in HAP1 cells in which SLC2A1 is required for growth. Donor libraries were introduced into the endogenous SLC2A1 gene in HAP1-Lig4KO cells using CRISPR/Cas9. Cell populations were harvested and sequenced to quantify the effect of variants on growth and generate a functional score. Quantitative functional scores were compared to 3-OMG uptake, SLC2A1 cell surface expression, CADD score, and clinical data, including CSF/blood glucose ratio.

RESULTS

Nonsense variants (N = 3) were reduced in cell culture over time resulting in negative scores (mean score: -1.15 ± 0.17), whereas synonymous variants (N = 10) were not depleted (mean score: 0.25 ± 0.12) (P < 2e-16). Missense variants (N = 27) yielded a range of functional scores including slightly negative scores, supporting a partial function and intermediate phenotype. Several variants with normal results on either cell surface expression (p.N34S and p.W65R) or 3-OMG uptake (p.W65R) had negative functional scores. There is a moderate but significant correlation between our functional scores and CADD scores.

INTERPRETATION

Cell growth is useful to quantitatively determine the functional effects of SLC2A1 variants. Nonsense variants were reliably distinguished from benign variants in this in vitro functional assay. For facilitating early diagnosis and therapeutic intervention, future work is needed to determine the functional effect of every possible variant in SLC2A1.

Systematic review of genotype-stratified treatment for monogenic insulin resistance

Objective

To assess the effects of pharmacologic and/or surgical interventions in monogenic insulin resistance (IR), stratified by genetic aetiology.

Design

Systematic review.

Data sources

PubMed, MEDLINE and Embase, from 1 January 1987 to 23 June 2021.

Review methods

Studies reporting individual-level effects of pharmacologic and/or surgical interventions in monogenic IR were eligible. Individual subject data were extracted and duplicate data removed. Outcomes were analyzed for each affected gene and intervention, and in aggregate for partial, generalised and all lipodystrophy.

Results

10 non-randomised experimental studies, 8 case series, and 21 single case reports met inclusion criteria, all rated as having moderate or serious risk of bias. Metreleptin was associated with lower triglycerides and hemoglobin A1c in aggregated lipodystrophy (n=111), in partial lipodystrophy (n=71) and generalised lipodystrophy (n=41)), and in LMNA , PPARG , AGPAT2 or BSCL2 subgroups (n=72,13,21 and 21 respectively). Body Mass Index (BMI) was lower after treatment in partial and generalised lipodystrophy overall, and in LMNA or BSCL2 , but not PPARG or AGPAT2 subgroups. Thiazolidinedione use was associated with improved hemoglobin A1c and triglycerides in aggregated lipodystrophy (n=13), improved hemoglobin A1c only in the PPARG subgroup (n=5), and improved triglycerides only in the LMNA subgroup (n=7). In INSR -related IR, use of rhIGF-1, alone or with IGFBP3, was associated with improved hemoglobin A1c (n=15). The small size or absence of all other genotype-treatment combinations precluded firm conclusions.

Conclusions

The evidence guiding genotype-specific treatment of monogenic IR is of low to very low quality. Metreleptin and Thiazolidinediones appear to have beneficial metabolic effects in lipodystrophy, and rhIGF-1 appears to lower hemoglobin A1c in INSR-related IR. For other interventions there is insufficient evidence to assess efficacy and risks either in aggregated lipodystrophy or in genetic subgroups. There is a pressing need to improve the evidence base for management of monogenic IR.

satmut_utils: a simulation and variant calling package for multiplexed assays of variant effect

The impact of millions of individual genetic variants on molecular phenotypes in coding sequences remains unknown. Multiplexed assays of variant effect (MAVEs) are scalable methods to annotate relevant variants, but existing software lacks standardization, requires cumbersome configuration, and does not scale to large targets. We present satmut_utils as a flexible solution for simulation and variant quantification. We then benchmark MAVE software using simulated and real MAVE data. We finally determine mRNA abundance for thousands of cystathionine beta-synthase variants using two experimental methods. The satmut_utils package enables high-performance analysis of MAVEs and reveals the capability of variants to alter mRNA abundance.

A Systematic Review of the use of Precision Diagnostics in Monogenic Diabetes

Monogenic forms of diabetes present opportunities for precision medicine as identification of the underlying genetic cause has implications for treatment and prognosis. However, genetic testing remains inconsistent across countries and health providers, often resulting in both missed diagnosis and misclassification of diabetes type. One of the barriers to deploying genetic testing is uncertainty over whom to test as the clinical features for monogenic diabetes overlap with those for both type 1 and type 2 diabetes. In this review, we perform a systematic evaluation of the evidence for the clinical and biochemical criteria used to guide selection of individuals with diabetes for genetic testing and review the evidence for the optimal methods for variant detection in genes involved in monogenic diabetes. In parallel we revisit the current clinical guidelines for genetic testing for monogenic diabetes and provide expert opinion on the interpretation and reporting of genetic tests. We provide a series of recommendations for the field informed by our systematic review, synthesizing evidence, and expert opinion. Finally, we identify major challenges for the field and highlight areas for future research and investment to support wider implementation of precision diagnostics for monogenic diabetes.

Molecular Mechanisms and Clinical Phenotypes of GJB2 Missense Variants

The GJB2 gene is the most common gene responsible for hearing loss (HL) worldwide, and missense variants are the most abundant type. GJB2 pathogenic missense variants cause nonsyndromic HL (autosomal recessive and dominant) and syndromic HL combined with skin diseases. However, the mechanism by which these different missense variants cause the different phenotypes is unknown. Over 2/3 of the GJB2 missense variants have yet to be functionally studied and are currently classified as variants of uncertain significance (VUS). Based on these functionally determined missense variants, we reviewed the clinical phenotypes and investigated the molecular mechanisms that affected hemichannel and gap junction functions, including connexin biosynthesis, trafficking, oligomerization into connexons, permeability, and interactions between other coexpressed connexins. We predict that all possible GJB2 missense variants will be described in the future by deep mutational scanning technology and optimizing computational models. Therefore, the mechanisms by which different missense variants cause different phenotypes will be fully elucidated.

Perinatal, metabolic, and reproductive features in PPARG-related lipodystrophy

OBJECTIVE

The adipogenic PPARG-encoded PPARγ nuclear receptor also displays essential placental functions. We evaluated the metabolic, reproductive, and perinatal features of patients with PPARG-related lipodystrophy.

METHODS

Current and retrospective data were collected in patients referred to a National Rare Diseases Reference Centre.

RESULTS

26 patients from 15 unrelated families were studied (18 women, median age 43 years). They carried monoallelic PPARG variants except a homozygous patient with congenital generalized lipodystrophy. Among heterozygous patients aged 16 or more (n = 24), 92% had diabetes, 96% partial lipodystrophy (median age at diagnosis 24 and 37 years), 78% hypertriglyceridaemia, 71% liver steatosis, and 58% hypertension. The mean BMI was 26 ± 5.0 kg/m2. Women (n = 16) were frequently affected by acute pancreatitis (n = 6) and/or polycystic ovary syndrome (n = 12). Eleven women obtained one or several pregnancies, all complicated by diabetes (n = 8), hypertension (n = 4), and/or hypertriglyceridaemia (n = 10). We analysed perinatal data of patients according to the presence (n = 8) or absence (n = 9) of a maternal dysmetabolic environment. The median gestational age at birth was low in both groups (37 and 36 weeks of amenorrhea, respectively). As expected, the birth weight was higher in patients exposed to a foetal dysmetabolic environment of maternal origin. In contrast, 85.7% of non-exposed patients, in whom the variant is, or is very likely to be, paternally-inherited, were small for gestational age.

CONCLUSIONS

Lipodystrophy-related PPARG variants induce early metabolic complications. Our results suggest that placental expression of PPARG pathogenic variants carried by affected foetuses could impair prenatal growth and parturition. This justifies careful pregnancy monitoring in affected families.

Hacking hematopoiesis - emerging tools for examining variant effects

Hematopoiesis is a continuous process of blood and immune cell production. It is orchestrated by thousands of gene products that respond to extracellular signals by guiding cell fate decisions to meet the needs of the organism. Although much of our knowledge of this process comes from work in model systems, we have learned a great deal from studies on human genetic variation. Considerable insight has emerged from studies on presumed monogenic blood disorders, which continue to provide key insights into the mechanisms critical for hematopoiesis. Furthermore, the emergence of large-scale biobanks and cohorts has uncovered thousands of genomic loci associated with blood cell traits and diseases. Some of these blood cell trait-associated loci act as modifiers of what were once thought to be monogenic blood diseases. However, most of these loci await functional validation. Here, we discuss the validation bottleneck and emerging methods to more effectively connect variant to function. In particular, we highlight recent innovations in genome editing, which have paved the path forward for high-throughput functional assessment of loci. Finally, we discuss existing barriers to progress, including challenges in manipulating the genomes of primary hematopoietic cells.

Deep mutational scanning: A versatile tool in systematically mapping genotypes to phenotypes

Unveiling how genetic variations lead to phenotypic variations is one of the key questions in evolutionary biology, genetics, and biomedical research. Deep mutational scanning (DMS) technology has allowed the mapping of tens of thousands of genetic variations to phenotypic variations efficiently and economically. Since its first systematic introduction about a decade ago, we have witnessed the use of deep mutational scanning in many research areas leading to scientific breakthroughs. Also, the methods in each step of deep mutational scanning have become much more versatile thanks to the oligo-synthesizing technology, high-throughput phenotyping methods and deep sequencing technology. However, each specific possible step of deep mutational scanning has its pros and cons, and some limitations still await further technological development. Here, we discuss recent scientific accomplishments achieved through the deep mutational scanning and describe widely used methods in each step of deep mutational scanning. We also compare these different methods and analyze their advantages and disadvantages, providing insight into how to design a deep mutational scanning study that best suits the aims of the readers' projects.

Disabling Uncompetitive Inhibition of Oncogenic IDH Mutations Drives Acquired Resistance

Mutations in IDH genes occur frequently in acute myeloid leukemia (AML) and other human cancers to generate the oncometabolite R-2HG. Allosteric inhibition of mutant IDH suppresses R-2HG production in a subset of patients with AML; however, acquired resistance emerges as a new challenge, and the underlying mechanisms remain incompletely understood. Here we establish isogenic leukemia cells containing common IDH oncogenic mutations by CRISPR base editing. By mutational scanning of IDH single amino acid variants in base-edited cells, we describe a repertoire of IDH second-site mutations responsible for therapy resistance through disabling uncompetitive enzyme inhibition. Recurrent mutations at NADPH binding sites within IDH heterodimers act in cis or trans to prevent the formation of stable enzyme-inhibitor complexes, restore R-2HG production in the presence of inhibitors, and drive therapy resistance in IDH-mutant AML cells and patients. We therefore uncover a new class of pathogenic mutations and mechanisms for acquired resistance to targeted cancer therapies.

SIGNIFICANCE

Comprehensive scanning of IDH single amino acid variants in base-edited leukemia cells uncovers recurrent mutations conferring resistance to IDH inhibition through disabling NADPH-dependent uncompetitive inhibition. Together with targeted sequencing, structural, and functional studies, we identify a new class of pathogenic mutations and mechanisms for acquired resistance to IDH-targeting cancer therapies. This article is highlighted in the In This Issue feature, p. 1.

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.

Insights on variant analysis in silico tools for pathogenicity prediction

Molecular biology is currently a fast-advancing science. Sequencing techniques are getting cheaper, but the interpretation of genetic variants requires expertise and computational power, therefore is still a challenge. Next-generation sequencing releases thousands of variants and to classify them, researchers propose protocols with several parameters. Here we present a review of several in silico pathogenicity prediction tools involved in the variant prioritization/classification process used by some international protocols for variant analysis and studies evaluating their efficiency.

PPARγ lipodystrophy mutants reveal intermolecular interactions required for enhancer activation

Peroxisome proliferator-activated receptor γ (PPARγ) is the master regulator of adipocyte differentiation, and mutations that interfere with its function cause lipodystrophy. PPARγ is a highly modular protein, and structural studies indicate that PPARγ domains engage in several intra- and inter-molecular interactions. How these interactions modulate PPARγ's ability to activate target genes in a cellular context is currently poorly understood. Here we take advantage of two previously uncharacterized lipodystrophy mutations, R212Q and E379K, that are predicted to interfere with the interaction of the hinge of PPARγ with DNA and with the interaction of PPARγ ligand binding domain (LBD) with the DNA-binding domain (DBD) of the retinoid X receptor, respectively. Using biochemical and genome-wide approaches we show that these mutations impair PPARγ function on an overlapping subset of target enhancers. The hinge region-DNA interaction appears mostly important for binding and remodelling of target enhancers in inaccessible chromatin, whereas the PPARγ-LBD:RXR-DBD interface stabilizes the PPARγ:RXR:DNA ternary complex. Our data demonstrate how in-depth analyses of lipodystrophy mutants can unravel molecular mechanisms of PPARγ function.

Defining the sensitivity landscape of EGFR variants to tyrosine kinase inhibitors

Tyrosine kinase inhibitor (TKI) is a standard treatment for patients with NSCLC harboring constitutively active epidermal growth factor receptor (EGFR) mutations. However, most rare EGFR mutations lack treatment regimens except for the well-studied ones. We constructed two EGFR variant libraries containing substitutions, deletions, or insertions using the saturation mutagenesis method. All the variants were located in the EGFR mutation hotspot (exons 18-21). The sensitivity of these variants to afatinib, erlotinib, gefitinib, icotinib, and osimertinib was systematically studied by determining their enrichment in massively parallel cytotoxicity assays using an endogenous EGFR-depleted cell line. A total of 3914 and 70,475 variants were detected in the constructed EGFR Substitution-Deletion (Sub-Del) and exon 20 Insertion (Ins) libraries. Of the 3914 Sub-Del variants, 221 proliferated fast in the control assay and were sensitive to EGFR-TKIs. For the 70,475 Ins variants, insertions at amino acid positions 770-774 were highly enriched in all 5 TKI cytotoxicity assays. Moreover, the top 5% of the enriched insertion variants included a glycine or serine insertion at high frequency. We present a comprehensive reference for the sensitivity of EGFR variants to five commonly used TKIs. The approach used here should be applicable to other genes and targeted drugs. BACKGROUND: Tyrosine kinase inhibitors (TKIs) therapy is a standard treatment for patients with advanced non-small-cell lung carcinoma (NSCLC) when activating epidermal growth factor receptor (EGFR) mutations are detected. However, except for the well-studied EGFR mutations, most EGFR mutations lack treatment regimens. TRANSLATIONAL SIGNIFICANCE: The results demonstrated that patients with rare EGFR mutations were most likely to benefit from osimertinib therapy compared to afatinib, erlotinib, gefitinib, or icotinib therapy. This study provides a case of deep mutational scanning that simultaneously assayed substitution, deletion, and insertion variants. This approach is applicable for other oncogenes and targeted drugs.

Achievements, prospects and challenges in precision care for monogenic insulin-deficient and insulin-resistant diabetes

Integration of genomic and other data has begun to stratify type 2 diabetes in prognostically meaningful ways, but this has yet to impact on mainstream diabetes practice. The subgroup of diabetes caused by single gene defects thus provides the best example to date of the vision of 'precision diabetes'. Monogenic diabetes may be divided into primary pancreatic beta cell failure, and primary insulin resistance. In both groups, clear examples of genotype-selective responses to therapy have been advanced. The benign trajectory of diabetes due to pathogenic GCK mutations, and the sulfonylurea-hyperresponsiveness conferred by activating KCNJ11 or ABCC8 mutations, or loss-of-function HNF1A or HNF4A mutations, often decisively guide clinical management. In monogenic insulin-resistant diabetes, subcutaneous leptin therapy is beneficial in some severe lipodystrophy. Increasing evidence also supports use of 'obesity therapies' in lipodystrophic people even without obesity. In beta cell diabetes the main challenge is now implementation of the precision diabetes vision at scale. In monogenic insulin-resistant diabetes genotype-specific benefits are proven in far fewer patients to date, although further genotype-targeted therapies are being evaluated. The conceptual paradigm established by the insulin-resistant subgroup with 'adipose failure' may have a wider influence on precision therapy for common type 2 diabetes, however. For all forms of monogenic diabetes, population-wide genome sequencing is currently forcing reappraisal of the importance assigned to pathogenic mutations when gene sequencing is uncoupled from prior suspicion of monogenic diabetes.

Intragenic compensation through the lens of deep mutational scanning

A significant fraction of mutations in proteins are deleterious and result in adverse consequences for protein function, stability, or interaction with other molecules. Intragenic compensation is a specific case of positive epistasis when a neutral missense mutation cancels effect of a deleterious mutation in the same protein. Permissive compensatory mutations facilitate protein evolution, since without them all sequences would be extremely conserved. Understanding compensatory mechanisms is an important scientific challenge at the intersection of protein biophysics and evolution. In human genetics, intragenic compensatory interactions are important since they may result in variable penetrance of pathogenic mutations or fixation of pathogenic human alleles in orthologous proteins from related species. The latter phenomenon complicates computational and clinical inference of an allele's pathogenicity. Deep mutational scanning is a relatively new technique that enables experimental studies of functional effects of thousands of mutations in proteins. We review the important aspects of the field and discuss existing limitations of current datasets. We reviewed ten published DMS datasets with quantified functional effects of single and double mutations and described rates and patterns of intragenic compensation in eight of them.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-022-01005-w.

Embeddings from protein language models predict conservation and variant effects

The emergence of SARS-CoV-2 variants stressed the demand for tools allowing to interpret the effect of single amino acid variants (SAVs) on protein function. While Deep Mutational Scanning (DMS) sets continue to expand our understanding of the mutational landscape of single proteins, the results continue to challenge analyses. Protein Language Models (pLMs) use the latest deep learning (DL) algorithms to leverage growing databases of protein sequences. These methods learn to predict missing or masked amino acids from the context of entire sequence regions. Here, we used pLM representations (embeddings) to predict sequence conservation and SAV effects without multiple sequence alignments (MSAs). Embeddings alone predicted residue conservation almost as accurately from single sequences as ConSeq using MSAs (two-state Matthews Correlation Coefficient-MCC-for ProtT5 embeddings of 0.596 ± 0.006 vs. 0.608 ± 0.006 for ConSeq). Inputting the conservation prediction along with BLOSUM62 substitution scores and pLM mask reconstruction probabilities into a simplistic logistic regression (LR) ensemble for Variant Effect Score Prediction without Alignments (VESPA) predicted SAV effect magnitude without any optimization on DMS data. Comparing predictions for a standard set of 39 DMS experiments to other methods (incl. ESM-1v, DeepSequence, and GEMME) revealed our approach as competitive with the state-of-the-art (SOTA) methods using MSA input. No method outperformed all others, neither consistently nor statistically significantly, independently of the performance measure applied (Spearman and Pearson correlation). Finally, we investigated binary effect predictions on DMS experiments for four human proteins. Overall, embedding-based methods have become competitive with methods relying on MSAs for SAV effect prediction at a fraction of the costs in computing/energy. Our method predicted SAV effects for the entire human proteome (~ 20 k proteins) within 40 min on one Nvidia Quadro RTX 8000. All methods and data sets are freely available for local and online execution through bioembeddings.com, https://github.com/Rostlab/VESPA , and PredictProtein.

Multiancestry exome sequencing reveals INHBE mutations associated with favorable fat distribution and protection from diabetes

Body fat distribution is a major, heritable risk factor for cardiometabolic disease, independent of overall adiposity. Using exome-sequencing in 618,375 individuals (including 160,058 non-Europeans) from the UK, Sweden and Mexico, we identify 16 genes associated with fat distribution at exome-wide significance. We show 6-fold larger effect for fat-distribution associated rare coding variants compared with fine-mapped common alleles, enrichment for genes expressed in adipose tissue and causal genes for partial lipodystrophies, and evidence of sex-dimorphism. We describe an association with favorable fat distribution (p = 1.8 × 10^-09), favorable metabolic profile and protection from type 2 diabetes (~28% lower odds; p = 0.004) for heterozygous protein-truncating mutations in INHBE, which encodes a circulating growth factor of the activin family, highly and specifically expressed in hepatocytes. Our results suggest that inhibin βE is a liver-expressed negative regulator of adipose storage whose blockade may be beneficial in fat distribution-associated metabolic disease.

Ion-ing out the genetic variants of Kir2.1

Deep mutational scanning provides new insights into how mutations alter the expression and activity of the potassium ion channel Kir2.1, which is associated with many diseases.

Loss-of-function, gain-of-function and dominant-negative mutations have profoundly different effects on protein structure

Most known pathogenic mutations occur in protein-coding regions of DNA and change the way proteins are made. Taking protein structure into account has therefore provided great insight into the molecular mechanisms underlying human genetic disease. While there has been much focus on how mutations can disrupt protein structure and thus cause a loss of function (LOF), alternative mechanisms, specifically dominant-negative (DN) and gain-of-function (GOF) effects, are less understood. Here, we investigate the protein-level effects of pathogenic missense mutations associated with different molecular mechanisms. We observe striking differences between recessive vs dominant, and LOF vs non-LOF mutations, with dominant, non-LOF disease mutations having much milder effects on protein structure, and DN mutations being highly enriched at protein interfaces. We also find that nearly all computational variant effect predictors, even those based solely on sequence conservation, underperform on non-LOF mutations. However, we do show that non-LOF mutations could potentially be identified by their tendency to cluster in three-dimensional space. Overall, our work suggests that many pathogenic mutations that act via DN and GOF mechanisms are likely being missed by current variant prioritisation strategies, but that there is considerable scope to improve computational predictions through consideration of molecular disease mechanisms.

Most known pathogenic mutations occur in protein-coding regions of DNA and change the way proteins are made. Here the authors analyse the locations of thousands of human disease mutations and their predicted effects on protein structure and show that,while loss-of-function mutations tend to be highly disruptive, non-loss-of-function mutations are in general much milder at a protein structural level.

High-throughput functional evaluation of human cancer-associated mutations using base editors

Comprehensive phenotypic characterization of the many mutations found in cancer tissues is one of the biggest challenges in cancer genomics. In this study, we evaluated the functional effects of 29,060 cancer-related transition mutations that result in protein variants on the survival and proliferation of non-tumorigenic lung cells using cytosine and adenine base editors and single guide RNA (sgRNA) libraries. By monitoring base editing efficiencies and outcomes using surrogate target sequences paired with sgRNA-encoding sequences on the lentiviral delivery construct, we identified sgRNAs that induced a single primary protein variant per sgRNA, enabling linking those mutations to the cellular phenotypes caused by base editing. The functions of the vast majority of the protein variants (28,458 variants, 98%) were classified as neutral or likely neutral; only 18 (0.06%) and 157 (0.5%) variants caused outgrowing and likely outgrowing phenotypes, respectively. We expect that our approach can be extended to more variants of unknown significance and other tumor types.

The genetic basis for adult onset glaucoma: Recent advances and future directions

Glaucoma, a diverse group of eye disorders that results in the degeneration of retinal ganglion cells, is the world's leading cause of irreversible blindness. Apart from age and ancestry, the major risk factor for glaucoma is increased intraocular pressure (IOP). In primary open-angle glaucoma (POAG), the anterior chamber angle is open but there is resistance to aqueous outflow. In primary angle-closure glaucoma (PACG), crowding of the anterior chamber angle due to anatomical alterations impede aqueous drainage through the angle. In exfoliation syndrome and exfoliation glaucoma, deposition of white flaky material throughout the anterior chamber directly interfere with aqueous outflow. Observational studies have established that there is a strong hereditable component for glaucoma onset and progression. Indeed, a succession of genome wide association studies (GWAS) that were centered upon single nucleotide polymorphisms (SNP) have yielded more than a hundred genetic markers associated with glaucoma risk. However, a shortcoming of GWAS studies is the difficulty in identifying the actual effector genes responsible for disease pathogenesis. Building on the foundation laid by GWAS studies, research groups have recently begun to perform whole exome-sequencing to evaluate the contribution of protein-changing, coding sequence genetic variants to glaucoma risk. The adoption of this technology in both large population-based studies as well as family studies are revealing the presence of novel, protein-changing genetic variants that could enrich our understanding of the pathogenesis of glaucoma. This review will cover recent advances in the genetics of primary open-angle glaucoma, primary angle-closure glaucoma and exfoliation glaucoma, which collectively make up the vast majority of all glaucoma cases in the world today. We will discuss how recent advances in research methodology have uncovered new risk genes, and how follow up biological investigations could be undertaken in order to define how the risk encoded by a genetic sequence variant comes into play in patients. We will also hypothesise how data arising from characterising these genetic variants could be utilized to predict glaucoma risk and the manner in which new therapeutic strategies might be informed.

Missense Variants Reveal Functional Insights Into the Human ARID Family of Gene Regulators

Missense variants are alterations to protein coding sequences that result in amino acid substitutions. They can be deleterious if the amino acid is required for maintaining structure or/and function, but are likely to be tolerated at other sites. Consequently, missense variation within a healthy population can mirror the effects of negative selection on protein structure and function, such that functional sites on proteins are often depleted of missense variants. Advances in high-throughput sequencing have dramatically increased the sample size of available human variation data, allowing for population-wide analysis of selective pressures. In this study, we developed a convenient set of tools, called 1D-to-3D, for visualizing the positions of missense variants on protein sequences and structures. We used these tools to characterize human homologues of the ARID family of gene regulators. ARID family members are implicated in multiple cancer types, developmental disorders, and immunological diseases but current understanding of their mechanistic roles is incomplete. Combined with phylogenetic and structural analyses, our approach allowed us to characterise sites important for protein-protein interactions, histone modification recognition, and DNA binding by the ARID proteins. We find that comparing missense depletion patterns among paralogs can reveal sub-functionalization at the level of domains. We propose that visualizing missense variants and their depletion on structures can serve as a valuable tool for complementing evolutionary and experimental findings.

When Adipose Tissue Lets You Down: Understanding the Functions of Genes Disrupted in Lipodystrophy

Lipodystrophy syndromes are conditions in which the adipose tissue mass of an individual is altered inappropriately. The change in adipose mass can range from a relatively modest and subtle redistribution in some individuals with partial lipodystrophy to a near-complete absence of adipose tissue in the most severe forms of generalized lipodystrophy. The common feature is a disconnection between the need of the individual for a safe, healthy lipid storage capacity and the available adipose mass to perform this critical role. The inability to partition lipids for storage in appropriately functioning adipocytes leads to lipid accumulation in other tissues, which typically results in conditions such as diabetes, dyslipidemia, fatty liver, and cardiovascular disease. Several genes have been identified whose disruption leads to inherited forms of lipodystrophy. There is a link between some of these genes and adipose dysfunction, so the molecular basis of disease pathophysiology appears clear. However, for other lipodystrophy genes, it is not evident why their disruption should affect adipose development or function or, in the case of partial lipodystrophy, why only some adipose depots should be affected. Elucidating the molecular functions of these genes and their cellular and physiological effects has the capacity to uncover fundamental new insights regarding the development and functions of adipose tissue. This information is also likely to inform better management of lipodystrophy and improved treatments for patients. In addition, the findings will often be relevant to other conditions featuring adipose tissue dysfunction, including the more common metabolic disease associated with obesity.

Benchmarking of SpCas9 variants enables deeper base editor screens of BRCA1 and BCL2

Numerous rationally-designed and directed-evolution variants of SpCas9 have been reported to expand the utility of CRISPR technology. Here, we assess the activity and specificity of WT-Cas9 and 10 SpCas9 variants by benchmarking their PAM preferences, on-target activity, and off-target susceptibility in cell culture assays with thousands of guides targeting endogenous genes. To enhance the coverage and thus utility of base editing screens, we demonstrate that the SpCas9-NG and SpG variants are compatible with both A > G and C > T base editors, more than tripling the number of guides and assayable residues. We demonstrate the performance of these technologies by screening for loss-of-function mutations in BRCA1 and Venetoclax-resistant mutations in BCL2, identifying both known and new mutations that alter function. We anticipate that the tools and methodologies described here will facilitate the investigation of genetic variants at a finer and deeper resolution for any locus of interest.

Numerous rationally-designed and directed-evolution variants of SpCas9 have been reported to expand the utility of CRISPR technology. Here the authors make comparisons of numerous Cas9 variants, nominate options for base editing screens with denser coverage with A>G and C>T base editing screens and identify loss-of-function mutations in BRCA1 and Venetoclax-resistant mutations in BCL2.

Isoform-specific functions of PPARγ in gene regulation and metabolism

In this study, Hu et al. investigated the specific functions of the two main PPARγ isoforms by generating mouse lines in which endogenous PPARγ1 and PPARγ2 were epitope-tagged to interrogate isoform-specific genomic binding, and mice deficient in either PPARγ1 or PPARγ2 to assess isoform-specific gene regulation. They show that PPARγ isoforms have specific and separable metabolic functions that may be targeted to improve therapy for insulin resistance and diabetes.

Peroxisome proliferator-activated receptor γ (PPARγ) is a nuclear receptor that is a vital regulator of adipogenesis, insulin sensitivity, and lipid metabolism. Activation of PPARγ by antidiabetic thiazolidinediones (TZD) reverses insulin resistance but also leads to weight gain that limits the use of these drugs. There are two main PPARγ isoforms, but the specific functions of each are not established. Here we generated mouse lines in which endogenous PPARγ1 and PPARγ2 were epitope-tagged to interrogate isoform-specific genomic binding, and mice deficient in either PPARγ1 or PPARγ2 to assess isoform-specific gene regulation. Strikingly, although PPARγ1 and PPARγ2 contain identical DNA binding domains, we uncovered isoform-specific genomic binding sites in addition to shared sites. Moreover, PPARγ1 and PPARγ2 regulated a different set of genes in adipose tissue depots, suggesting distinct roles in adipocyte biology. Indeed, mice with selective deficiency of PPARγ1 maintained body temperature better than wild-type or PPARγ2-deficient mice. Most remarkably, although TZD treatment improved glucose tolerance in mice lacking either PPARγ1 or PPARγ2, the PPARγ1-deficient mice were protected from TZD-induced body weight gain compared with PPARγ2-deficient mice. Thus, PPARγ isoforms have specific and separable metabolic functions that may be targeted to improve therapy for insulin resistance and diabetes.

Current and Future Approaches to Classify VUSs in LGMD-Related Genes

Next-generation sequencing (NGS) has revealed large numbers of genetic variants in LGMD-related genes, with most of them classified as variants of uncertain significance (VUSs). VUSs are genetic changes with unknown pathological impact and present a major challenge in genetic test interpretation and disease diagnosis. Understanding the phenotypic consequences of VUSs can provide clinical guidance regarding LGMD risk and therapy. In this review, we provide a brief overview of the subtypes of LGMD, disease diagnosis, current classification systems for investigating VUSs, and a potential deep mutational scanning approach to classify VUSs in LGMD-related genes.

Contributions of microRNAs to Peripheral Insulin Sensitivity

An extensive literature base combined with advances in sequencing technologies demonstrate microRNA levels correlate with various metabolic diseases. Mechanistic studies also establish microRNAs regulate central metabolic pathways and thus play vital roles in maintaining organismal energy balance and metabolic homeostasis. This review highlights research progress on the roles and regulation of microRNAs in the peripheral tissues that confer insulin sensitivity. We discuss sequencing technologies used to comprehensively define the target spectrum of microRNAs in metabolic disease that complement studies reporting physiologic roles for microRNAs in the regulation of glucose and lipid metabolism in animal models. We also discuss the emerging roles of exosomal microRNAs as endocrine signals to regulate lipid and carbohydrate metabolism.

Directed mutational scanning reveals a balance between acidic and hydrophobic residues in strong human activation domains

Acidic activation domains are intrinsically disordered regions of the transcription factors that bind coactivators. The intrinsic disorder and low evolutionary conservation of activation domains have made it difficult to identify the sequence features that control activity. To address this problem, we designed thousands of variants in seven acidic activation domains and measured their activities with a high-throughput assay in human cell culture. We found that strong activation domain activity requires a balance between the number of acidic residues and aromatic and leucine residues. These findings motivated a predictor of acidic activation domains that scans the human proteome for clusters of aromatic and leucine residues embedded in regions of high acidity. This predictor identifies known activation domains and accurately predicts previously unidentified ones. Our results support a flexible acidic exposure model of activation domains in which the acidic residues solubilize hydrophobic motifs so that they can interact with coactivators. A record of this paper’s transparent peer review process is included in the supplemental information.

Transcriptional activation domains are poorly conserved, intrinsically disordered regions of the transcription factors that remain difficult to predict from protein sequences. A high-throughput method reveals how strong activation domains require a balance between acidic and hydrophobic residues. This balance powers an accurate predictor of activation domains on human transcription factors.

Case Report: A New Peroxisome Proliferator-Activated Receptor Gamma Mutation Causes Familial Partial Lipodystrophy Type 3 in a Chinese Patient

PURPOSE

Familial partial lipodystrophy type 3 (FPLD3) is an autosomal dominant disease. Patients typically present with loss of adipose tissue and metabolic complications. Here, we reported a Chinese FPLD3 patient with a novel PPARG gene mutation.

METHODS

A 16-year-old female patient and her relatives were assessed by detailed clinical and biochemical examinations. Sequencing was performed by using the extracted DNA. Moreover, we identified FPLD3 patients from previous studies, and according to the protein region affected by the gene mutation. We divided the patients into the DNA-binding domain (DBD) group or the ligand-binding domain (LBD) group, and compared the clinical features between the two groups.

RESULTS

We identified a novel gene mutation affecting the LBD of PPARγ c.929T > C (p.F310S). This mutation leads to the substitution of a phenylalanine by a serine. In our case, subcutaneous fat was significantly diminished in her face, hips and limbs. The patient was also presented with insulin resistance, diabetes mellitus, hypertriglyceridemia, fatty liver, liver dysfunction, albuminuria and diabetic peripheral neuropathy. After literature review, a total of 58 FPLD3 patients were identified and we found no difference in clinical features between the DBD group and LBD group (all P > 0.05).

CONCLUSIONS

A Chinese FPLD3 patient with a novel PPARG gene mutation is described. Our case emphasized the importance of physical examination and genetic testing in young patients with severe metabolic syndromes.

Impact of the Mutational Landscape of the Sodium/Iodide Symporter in Congenital Hypothyroidism

Background: Iodide transport defect is an uncommon cause of dyshormonogenic congenital hypothyroidism due to homozygous or compound heterozygous pathogenic variants in the SLC5A5 gene, which encodes the sodium/iodide symporter (NIS), causing deficient iodide accumulation in thyroid follicular cells, thus impairing thyroid hormonogenesis. Methods:SLC5A5 gene variants were compiled from public databases and research articles exploring the molecular bases of congenital hypothyroidism. Using a dataset of 198 missense NIS variants classified as either benign or pathogenic, we developed and validated a machine learning-based NIS-specific variant classifier to predict the impact of missense NIS variants. Results: We generated a manually curated dataset containing 7793 unique SLC5A5 variants. As most databases compiled exome sequencing data, variant mapping revealed an increased density of variants in SLC5A5 coding exons. Based on allele frequency (AF) analysis, we established an AF threshold of 1:10,000 above which a variant should be considered benign. Most pathogenic NIS variants were located in the protein-coding region, as most patients were genetically diagnosed by using a candidate gene strategy limited to this region. Significantly, we evidenced that 94.5% of missense NIS variants were classified as of uncertain significance. Therefore, we developed an NIS-specific variant classifier to improve the prediction of pathogenicity of missense variants. Our classifier predicted the clinical outcome of missense variants with high accuracy (90%), outperforming state-of-the-art pathogenicity predictors, such as REVEL, PolyPhen-2, and SIFT. Based on the excellent performance of our classifier, we predicted the mutational landscape of NIS. The analysis of the mutational landscape revealed that most missense variants located in transmembrane segments are frequently pathogenic. Moreover, we predicted that ∼28% of all single-nucleotide variants that could cause missense NIS variants are pathogenic, thus putatively leading to congenital hypothyroidism if present in homozygous or compound heterozygous state. Conclusions: We reported the first NIS-specific variant classifier aiming at improving the interpretation of missense NIS variants in clinical practice. Deciphering the mutational landscape for every protein involved in thyroid hormonogenesis is a relevant task for a deep understanding of the molecular mechanisms causing dyshormonogenic congenital hypothyroidism.