The genetic regulatory signature of type 2 diabetes in human skeletal muscle

CRISPR Screen Identifies the RNA-Binding Protein Eef1a1 as a Key Regulator of Myogenesis

Skeletal muscle myogenesis hinges on gene regulation, meticulously orchestrated by molecular mechanisms. While the roles of transcription factors and non-coding RNAs in myogenesis are widely known, the contribution of RNA-binding proteins (RBPs) has remained unclear until now. Therefore, to investigate the functions of post-transcriptional regulators in myogenesis and uncover new functional RBPs regulating myogenesis, we employed CRISPR high-throughput RBP-KO (RBP-wide knockout) library screening. Through this approach, we successfully identified Eef1a1 as a novel regulatory factor in myogenesis. Using CRISPR knockout (CRISPRko) and CRISPR interference (CRISPRi) technologies, we successfully established cellular models for both CRISPRko and CRISPRi. Our findings demonstrated that Eef1a1 plays a crucial role in promoting proliferation in C2C12 myoblasts. Through siRNA inhibition and overexpression methods, we further elucidated the involvement of Eef1a1 in promoting proliferation and suppressing differentiation processes. RIP (RNA immunoprecipitation), miRNA pull-down, and Dual-luciferase reporter assays confirmed that miR-133a-3p targets Eef1a1. Co-transfection experiments indicated that miR-133a-3p can rescue the effect of Eef1a1 on C2C12 myoblasts. In summary, our study utilized CRISPR library high-throughput screening to unveil a novel RBP, Eef1a1, involved in regulating myogenesis. Eef1a1 promotes the proliferation of myoblasts while inhibiting the differentiation process. Additionally, it acts as an antagonist to miR-133a-3p, thus modulating the process of myogenesis.

Cardiometabolic characteristics of people with metabolically healthy and unhealthy obesity

There is considerable heterogeneity in the cardiometabolic abnormalities associated with obesity. We evaluated multi-organ system metabolic function in 20 adults with metabolically healthy obesity (MHO; normal fasting glucose and triglycerides, oral glucose tolerance, intrahepatic triglyceride content, and whole-body insulin sensitivity), 20 adults with metabolically unhealthy obesity (MUO; prediabetes, hepatic steatosis, and whole-body insulin resistance), and 15 adults who were metabolically healthy lean. Compared with MUO, people with MHO had (1) altered skeletal muscle biology (decreased ceramide content and increased expression of genes involved in BCAA catabolism and mitochondrial structure/function); (2) altered adipose tissue biology (decreased expression of genes involved in inflammation and extracellular matrix remodeling and increased expression of genes involved in lipogenesis); (3) lower 24-h plasma glucose, insulin, non-esterified fatty acids, and triglycerides; (4) higher plasma adiponectin and lower plasma PAI-1 concentrations; and (5) decreased oxidative stress. These findings provide a framework of potential mechanisms responsible for MHO and the metabolic heterogeneity of obesity. This study was registered at ClinicalTrials.gov (NCT02706262).

Common and ethnic-specific derangements in skeletal muscle transcriptome associated with obesity

BACKGROUND

Obesity is a common disease with a higher prevalence among African Americans. Obesity alters cellular function in many tissues, including skeletal muscle, and is a risk factor for many life-threatening diseases, including cardiovascular disease and diabetes. The similarities and differences in molecular mechanisms that may explain ethnic disparities in obesity between African and European ancestry individuals have not been studied.

METHODS

In this study, data from transcriptome-wide analyses on skeletal muscle tissues from well-powered human cohorts were used to compare genes and biological pathways affected by obesity in European and African ancestry populations. Data on obesity-induced differentially expressed transcripts and GWAS-identified SNPs were integrated to prioritize target genes for obesity-associated genetic variants.

RESULTS

Linear regression analysis in the FUSION (European, N = 301) and AAGMEx (African American, N = 256) cohorts identified a total of 2569 body mass index (BMI)-associated transcripts (q < 0.05), of which 970 genes (at p < 0.05) are associated in both cohorts, and the majority showed the same direction of effect on BMI. Biological pathway analyses, including over-representation and gene-set enrichment analyses, identified enrichment of protein synthesis pathways (e.g., ribosomal function) and the ceramide signaling pathway in both cohorts among BMI-associated down- and up-regulated transcripts, respectively. A comparison using the IPA-tool suggested the activation of inflammation pathways only in Europeans with obesity. Interestingly, these analyses suggested repression of the mitochondrial oxidative phosphorylation pathway in Europeans but showed its activation in African Americans. Integration of SNP-to-Gene analyses-predicted target genes for obesity-associated genetic variants (GWAS-identified SNPs) and BMI-associated transcripts suggested that these SNPs might cause obesity by altering the expression of 316 critical target genes (e.g., GRB14) in the muscle.

CONCLUSIONS

This study provides a replication of obesity-associated transcripts and biological pathways in skeletal muscle across ethnicities, but also identifies obesity-associated processes unique in either African or European ancestry populations.

Mitochondrial Dysfunction, Oxidative Stress, and Inter-Organ Miscommunications in T2D Progression

Type 2 diabetes (T2D) is a heterogenous disease, and conventionally, peripheral insulin resistance (IR) was thought to precede islet β-cell dysfunction, promoting progression from prediabetes to T2D. New evidence suggests that T2D-lean individuals experience early β-cell dysfunction without significant IR. Regardless of the primary event (i.e., IR vs. β-cell dysfunction) that contributes to dysglycemia, significant early-onset oxidative damage and mitochondrial dysfunction in multiple metabolic tissues may be a driver of T2D onset and progression. Oxidative stress, defined as the generation of reactive oxygen species (ROS), is mediated by hyperglycemia alone or in combination with lipids. Physiological oxidative stress promotes inter-tissue communication, while pathological oxidative stress promotes inter-tissue mis-communication, and new evidence suggests that this is mediated via extracellular vesicles (EVs), including mitochondria containing EVs. Under metabolic-related stress conditions, EV-mediated cross-talk between β-cells and skeletal muscle likely trigger mitochondrial anomalies leading to prediabetes and T2D. This article reviews the underlying molecular mechanisms in ROS-related pathogenesis of prediabetes, including mitophagy and mitochondrial dynamics due to oxidative stress. Further, this review will describe the potential of various therapeutic avenues for attenuating oxidative damage, reversing prediabetes and preventing progression to T2D.

Population-scale skeletal muscle single-nucleus multi-omic profiling reveals extensive context specific genetic regulation

Skeletal muscle, the largest human organ by weight, is relevant to several polygenic metabolic traits and diseases including type 2 diabetes (T2D). Identifying genetic mechanisms underlying these traits requires pinpointing the relevant cell types, regulatory elements, target genes, and causal variants. Here, we used genetic multiplexing to generate population-scale single nucleus (sn) chromatin accessibility (snATAC-seq) and transcriptome (snRNA-seq) maps across 287 frozen human skeletal muscle biopsies representing 456,880 nuclei. We identified 13 cell types that collectively represented 983,155 ATAC summits. We integrated genetic variation to discover 6,866 expression quantitative trait loci (eQTL) and 100,928 chromatin accessibility QTL (caQTL) (5% FDR) across the five most abundant cell types, cataloging caQTL peaks that atlas-level snATAC maps often miss. We identified 1,973 eGenes colocalized with caQTL and used mediation analyses to construct causal directional maps for chromatin accessibility and gene expression. 3,378 genome-wide association study (GWAS) signals across 43 relevant traits colocalized with sn-e/caQTL, 52% in a cell-specific manner. 77% of GWAS signals colocalized with caQTL and not eQTL, highlighting the critical importance of population-scale chromatin profiling for GWAS functional studies. GWAS-caQTL colocalization showed distinct cell-specific regulatory paradigms. For example, a C2CD4A/B T2D GWAS signal colocalized with caQTL in muscle fibers and multiple chromatin loop models nominated VPS13C, a glucose uptake gene. Sequence of the caQTL peak overlapping caSNP rs7163757 showed allelic regulatory activity differences in a human myocyte cell line massively parallel reporter assay. These results illuminate the genetic regulatory architecture of human skeletal muscle at high-resolution epigenomic, transcriptomic, and cell state scales and serve as a template for population-scale multi-omic mapping in complex tissues and traits.

Adipose tissue eQTL meta-analysis reveals the contribution of allelic heterogeneity to gene expression regulation and cardiometabolic traits

Complete characterization of the genetic effects on gene expression is needed to elucidate tissue biology and the etiology of complex traits. Here, we analyzed 2,344 subcutaneous adipose tissue samples and identified 34K conditionally distinct expression quantitative trait locus (eQTL) signals in 18K genes. Over half of eQTL genes exhibited at least two eQTL signals. Compared to primary signals, non-primary signals had lower effect sizes, lower minor allele frequencies, and less promoter enrichment; they corresponded to genes with higher heritability and higher tolerance for loss of function. Colocalization of eQTL with conditionally distinct genome-wide association study signals for 28 cardiometabolic traits identified 3,605 eQTL signals for 1,861 genes. Inclusion of non-primary eQTL signals increased colocalized signals by 46%. Among 30 genes with ≥2 pairs of colocalized signals, 21 showed a mediating gene dosage effect on the trait. Thus, expanded eQTL identification reveals more mechanisms underlying complex traits and improves understanding of the complexity of gene expression regulation.

Dysregulation of early gene response to a mixed meal in skeletal muscle in obesity and type 2 diabetes

Obesity- and type 2 diabetes mellitus-induced changes in the expression of protein-coding genes in human skeletal muscle were extensively examined at baseline (after an overnight fast). We aimed to compare the early transcriptomic response to a typical single meal in skeletal muscle of metabolically healthy subjects and obese individuals without and with type 2 diabetes. Transcriptomic response (RNA-seq) to a mixed meal (nutritional drink, ∼25 kJ/kg of body mass) was examined in the vastus lateralis muscle (1 h after a meal) in 7 healthy subjects and 14 obese individuals without or with type 2 diabetes. In all obese individuals, the transcriptome response to a meal was dysregulated (suppressed and altered) and associated with different biological processes compared with healthy control. To search for potential transcription factors regulating transcriptomic response to a meal, the enrichment of transcription factor-binding sites in individual promoters of the human skeletal muscle was examined. In obese individuals, the transcriptomic response is associated with a different set of transcription factors than that in healthy subjects. In conclusion, metabolic disorders are associated with a defect in the regulation of mixed meal/insulin-mediated gene expression-insulin resistance in terms of gene expression. Importantly, this dysregulation occurs in obese individuals without type 2 diabetes, i.e., at the first stage of the development of metabolic disorders.NEW & NOTEWORTHY In skeletal muscle of metabolically healthy subjects, a typical single meal normalized to body mass induces activation of various transcription factors, expression of numerous receptor tyrosine kinases associated with the insulin signaling cascade, and transcription regulators. In skeletal muscle of obese individuals without and with type 2 diabetes, this signaling network is poorly regulated at the transcriptional level, indicating dysregulation of the early gene response to a mixed meal.

Toward Ameliorating Insulin Resistance: Targeting a Novel PAK1 Signaling Pathway Required for Skeletal Muscle Mitochondrial Function

The p21-activated kinase 1 (PAK1) is required for insulin-stimulated glucose uptake in skeletal muscle cells. However, whether PAK1 regulates skeletal muscle mitochondrial function, which is a central determinant of insulin sensitivity, is unknown. Here, the effect of modulating PAK1 levels (knockdown via siRNA, overexpression via adenoviral transduction, and/or inhibition of activation via IPA3) on mitochondrial function was assessed in normal and/or insulin-resistant rat L6.GLUT4myc and human muscle (LHCN-M2) myotubes. Human type 2 diabetes (T2D) and non-diabetic (ND) skeletal muscle samples were also used for validation of the identified signaling elements. PAK1 depletion in myotubes decreased mitochondrial copy number, respiration, altered mitochondrial structure, downregulated PGC1α (a core regulator of mitochondrial biogenesis and oxidative metabolism) and PGC1α activators, p38 mitogen-activated protein kinase (p38MAPK) and activating transcription factor 2 (ATF2). PAK1 enrichment in insulin-resistant myotubes improved mitochondrial function and rescued PGC1α expression levels. Activated PAK1 was localized to the cytoplasm, and PAK1 enrichment concurrent with p38MAPK inhibition did not increase PGC1α levels. PAK1 inhibition and enrichment also modified nuclear phosphorylated-ATF2 levels. T2D human samples showed a deficit for PGC1α, and PAK1 depletion in LHCN-M2 cells led to reduced mitochondrial respiration. Overall, the results suggest that PAK1 regulates muscle mitochondrial function upstream of the p38MAPK/ATF2/PGC1α-axis pathway.

An integrative profiling of metabolome and transcriptome in the plasma and skeletal muscle following an exercise intervention in diet-induced obese mice

Exercise intervention at the early stage of type 2 diabetes mellitus (T2DM) can aid in the maintenance of blood glucose homeostasis and prevent the development of macrovascular and microvascular complications. However, the exercise-regulated pathways that prevent the development of T2DM remain largely unclear. In this study, two forms of exercise intervention, treadmill training and voluntary wheel running, were conducted for high-fat diet (HFD)-induced obese mice. We observed that both forms of exercise intervention alleviated HFD-induced insulin resistance and glucose intolerance. Skeletal muscle is recognized as the primary site for postprandial glucose uptake and for responsive alteration beyond exercise training. Metabolomic profiling of the plasma and skeletal muscle in Chow, HFD, and HFD-exercise groups revealed robust alterations in metabolic pathways by exercise intervention in both cases. Overlapping analysis identified nine metabolites, including beta-alanine, leucine, valine, and tryptophan, which were reversed by exercise treatment in both the plasma and skeletal muscle. Transcriptomic analysis of gene expression profiles in the skeletal muscle revealed several key pathways involved in the beneficial effects of exercise on metabolic homeostasis. In addition, integrative transcriptomic and metabolomic analyses uncovered strong correlations between the concentrations of bioactive metabolites and the expression levels of genes involved in energy metabolism, insulin sensitivity, and immune response in the skeletal muscle. This work established two models of exercise intervention in obese mice and provided mechanistic insights into the beneficial effects of exercise intervention on systemic energy homeostasis.

A holistic view of muscle metabolic reprogramming through personalized metabolic modeling in newly diagnosed diabetic patients

Type 2 diabetes mellitus (T2DM) is a challenging and progressive metabolic disease caused by insulin resistance. Skeletal muscle is the major insulin-sensitive tissue that plays a pivotal role in blood sugar homeostasis. Dysfunction of muscle metabolism is implicated in the disturbance of glucose homeostasis, the development of insulin resistance, and T2DM. Understanding metabolism reprogramming in newly diagnosed patients provides opportunities for early diagnosis and treatment of T2DM as a challenging disease to manage. Here, we applied a system biology approach to investigate metabolic dysregulations associated with the early stage of T2DM. We first reconstructed a human muscle-specific metabolic model. The model was applied for personalized metabolic modeling and analyses in newly diagnosed patients. We found that several pathways and metabolites, mainly implicating in amino acids and lipids metabolisms, were dysregulated. Our results indicated the significance of perturbation of pathways implicated in building membrane and extracellular matrix (ECM). Dysfunctional metabolism in these pathways possibly interrupts the signaling process and develops insulin resistance. We also applied a machine learning method to predict potential metabolite markers of insulin resistance in skeletal muscle. 13 exchange metabolites were predicted as the potential markers. The efficiency of these markers in discriminating insulin-resistant muscle was successfully validated.

Transcriptomics for Clinical and Experimental Biology Research: Hang on a Seq

Sequencing the human genome empowers translational medicine, facilitating transcriptome-wide molecular diagnosis, pathway biology, and drug repositioning. Initially, microarrays are used to study the bulk transcriptome; but now short-read RNA sequencing (RNA-seq) predominates. Positioned as a superior technology, that makes the discovery of novel transcripts routine, most RNA-seq analyses are in fact modeled on the known transcriptome. Limitations of the RNA-seq methodology have emerged, while the design of, and the analysis strategies applied to, arrays have matured. An equitable comparison between these technologies is provided, highlighting advantages that modern arrays hold over RNA-seq. Array protocols more accurately quantify constitutively expressed protein coding genes across tissue replicates, and are more reliable for studying lower expressed genes. Arrays reveal long noncoding RNAs (lncRNA) are neither sparsely nor lower expressed than protein coding genes. Heterogeneous coverage of constitutively expressed genes observed with RNA-seq, undermines the validity and reproducibility of pathway analyses. The factors driving these observations, many of which are relevant to long-read or single-cell sequencing are discussed. As proposed herein, a reappreciation of bulk transcriptomic methods is required, including wider use of the modern high-density array data-to urgently revise existing anatomical RNA reference atlases and assist with more accurate study of lncRNAs.

Skeletal muscle overexpression of sAnk1.5 in transgenic mice does not predispose to type 2 diabetes

Genome-wide association studies (GWAS) and cis-expression quantitative trait locus (cis-eQTL) analyses indicated an association of the rs508419 single nucleotide polymorphism (SNP) with type 2 diabetes (T2D). rs508419 is localized in the muscle-specific internal promoter (P2) of the ANK1 gene, which drives the expression of the sAnk1.5 isoform. Functional studies showed that the rs508419 C/C variant results in increased transcriptional activity of the P2 promoter, leading to higher levels of sAnk1.5 mRNA and protein in skeletal muscle biopsies of individuals carrying the C/C genotype. To investigate whether sAnk1.5 overexpression in skeletal muscle might predispose to T2D development, we generated transgenic mice (TgsAnk1.5/+) in which the sAnk1.5 coding sequence was selectively overexpressed in skeletal muscle tissue. TgsAnk1.5/+ mice expressed up to 50% as much sAnk1.5 protein as wild-type (WT) muscles, mirroring the difference reported between individuals with the C/C or T/T genotype at rs508419. However, fasting glucose levels, glucose tolerance, insulin levels and insulin response in TgsAnk1.5/+ mice did not differ from those of age-matched WT mice monitored over a 12-month period. Even when fed a high-fat diet, TgsAnk1.5/+ mice only presented increased caloric intake, but glucose disposal, insulin tolerance and weight gain were comparable to those of WT mice fed a similar diet. Altogether, these data indicate that sAnk1.5 overexpression in skeletal muscle does not predispose mice to T2D susceptibility.

Prioritization of genes associated with type 2 diabetes mellitus for functional studies

Existing therapies for type 2 diabetes mellitus (T2DM) show limited efficacy or have adverse effects. Numerous genetic variants associated with T2DM have been identified, but progress in translating these findings into potential drug targets has been limited. Here, we describe the tools and platforms available to identify effector genes from T2DM-associated coding and non-coding variants and prioritize them for functional studies. We discuss QSER1 and SLC12A8 as examples of genes that have been identified as possible T2DM candidate genes using these tools and platforms. We suggest further approaches, including the use of sequencing data with increased sample size and ethnic diversity, single-cell omics data for analyses, glycaemic trait associations to predict gene function and, potentially, human induced pluripotent stem cell 'village' cultures, to strengthen current gene functionalization workflows. Effective prioritization of T2DM-associated genes for experimental validation could expedite our understanding of the genetic mechanisms responsible for T2DM to facilitate the use of precision medicine in its treatment.

Effect of tissue-grouped regulatory variants associated to type 2 diabetes in related secondary outcomes

Genome-wide association studies have identified over five hundred loci that contribute to variation in type 2 diabetes (T2D), an established risk factor for many diseases. However, the mechanisms and extent through which these loci contribute to subsequent outcomes remain elusive. We hypothesized that combinations of T2D-associated variants acting on tissue-specific regulatory elements might account for greater risk for tissue-specific outcomes, leading to diversity in T2D disease progression. We searched for T2D-associated variants acting on regulatory elements and expression quantitative trait loci (eQTLs) in nine tissues. We used T2D tissue-grouped variant sets as genetic instruments to conduct 2-Sample Mendelian Randomization (MR) in ten related outcomes whose risk is increased by T2D using the FinnGen cohort. We performed PheWAS analysis to investigate whether the T2D tissue-grouped variant sets had specific predicted disease signatures. We identified an average of 176 variants acting in nine tissues implicated in T2D, and an average of 30 variants acting on regulatory elements that are unique to the nine tissues of interest. In 2-Sample MR analyses, all subsets of regulatory variants acting in different tissues were associated with increased risk of the ten secondary outcomes studied on similar levels. No tissue-grouped variant set was associated with an outcome significantly more than other tissue-grouped variant sets. We did not identify different disease progression profiles based on tissue-specific regulatory and transcriptome information. Bigger sample sizes and other layers of regulatory information in critical tissues may help identify subsets of T2D variants that are implicated in certain secondary outcomes, uncovering system-specific disease progression.

Roles of super enhancers and enhancer RNAs in skeletal muscle development and disease

Skeletal muscle development is a multistep biological process regulated by a variety of myogenic regulatory factors, including MyoG, MyoD, Myf5, and Myf6 (also known as MRF4), as well as members of the FoxO subfamily. Differentiation and regeneration during skeletal muscle myogenesis contribute to the physiological function of muscles. Super enhancers (SEs) and enhancer RNAs (eRNAs) are involved in the regulation of development and diseases. Few studies have identified the roles of SEs and eRNAs in muscle development and pathophysiology. To develop approaches to enhance skeletal muscle mass and function, a more comprehensive understanding of the key processes underlying muscular diseases is needed. In this review, we summarize the roles of SEs and eRNAs in muscle development and disease through affecting of DNA methylation, FoxO subfamily, RAS-MEK signaling, chromatin modifications and accessibility, MyoD and cis regulating target genes. The summary could inform strategies to increase muscle mass and treat muscle-related diseases.

Differential DNA methylation of steatosis and non-alcoholic fatty liver disease in adolescence

BACKGROUND AND AIMS

Epigenetic modifications are associated with hepatic fat accumulation and non-alcoholic fatty liver disease (NAFLD). However, few epigenetic modifications directly implicated in such processes have been identified during adolescence, a critical developmental window where physiological changes could influence future disease trajectory. To investigate the association between DNA methylation and NAFLD in adolescence, we undertook discovery and validation of novel methylation marks, alongside replication of previously reported marks.

APPROACH AND RESULTS

We performed a DNA methylation epigenome-wide association study (EWAS) on DNA from whole blood from 707 Raine Study adolescents phenotyped for steatosis score and NAFLD by ultrasound at age 17. Next, we performed pyrosequencing validation of loci within the most 100 strongly associated differentially methylated CpG sites (dmCpGs) for which ≥ 2 probes per gene remained significant across four statistical models with a nominal p value < 0.007. EWAS identified dmCpGs related to three genes (ANK1, MIR10a, PTPRN2) that met our criteria for pyrosequencing. Of the dmCpGs and surrounding loci that were pyrosequenced (ANK1 n = 6, MIR10a n = 7, PTPRN2 n = 3), three dmCpGs in ANK1 and two in MIR10a were significantly associated with NAFLD in adolescence. After adjustment for waist circumference only dmCpGs in ANK1 remained significant. These ANK1 CpGs were also associated with γ-glutamyl transferase and alanine aminotransferase concentrations. Three of twenty-two differentially methylated dmCpGs previously associated with adult NAFLD were associated with NAFLD in adolescence (all adjusted p < 2.3 × 10^-3).

CONCLUSIONS

We identified novel DNA methylation loci associated with NAFLD and serum liver biochemistry markers during adolescence, implicating putative dmCpG/gene regulatory pathways and providing insights for future mechanistic studies.

Loci for insulin processing and secretion provide insight into type 2 diabetes risk

Insulin secretion is critical for glucose homeostasis, and increased levels of the precursor proinsulin relative to insulin indicate pancreatic islet beta-cell stress and insufficient insulin secretory capacity in the setting of insulin resistance. We conducted meta-analyses of genome-wide association results for fasting proinsulin from 16 European-ancestry studies in 45,861 individuals. We found 36 independent signals at 30 loci (p value < 5 × 10-8), which validated 12 previously reported loci for proinsulin and ten additional loci previously identified for another glycemic trait. Half of the alleles associated with higher proinsulin showed higher rather than lower effects on glucose levels, corresponding to different mechanisms. Proinsulin loci included genes that affect prohormone convertases, beta-cell dysfunction, vesicle trafficking, beta-cell transcriptional regulation, and lysosomes/autophagy processes. We colocalized 11 proinsulin signals with islet expression quantitative trait locus (eQTL) data, suggesting candidate genes, including ARSG, WIPI1, SLC7A14, and SIX3. The NKX6-3/ANK1 proinsulin signal colocalized with a T2D signal and an adipose ANK1 eQTL signal but not the islet NKX6-3 eQTL. Signals were enriched for islet enhancers, and we showed a plausible islet regulatory mechanism for the lead signal in the MADD locus. These results show how detailed genetic studies of an intermediate phenotype can elucidate mechanisms that may predispose one to disease.

Network reconstruction for trans acting genetic loci using multi-omics data and prior information

BACKGROUND

Molecular measurements of the genome, the transcriptome, and the epigenome, often termed multi-omics data, provide an in-depth view on biological systems and their integration is crucial for gaining insights in complex regulatory processes. These data can be used to explain disease related genetic variants by linking them to intermediate molecular traits (quantitative trait loci, QTL). Molecular networks regulating cellular processes leave footprints in QTL results as so-called trans-QTL hotspots. Reconstructing these networks is a complex endeavor and use of biological prior information can improve network inference. However, previous efforts were limited in the types of priors used or have only been applied to model systems. In this study, we reconstruct the regulatory networks underlying trans-QTL hotspots using human cohort data and data-driven prior information.

METHODS

We devised a new strategy to integrate QTL with human population scale multi-omics data. State-of-the art network inference methods including BDgraph and glasso were applied to these data. Comprehensive prior information to guide network inference was manually curated from large-scale biological databases. The inference approach was extensively benchmarked using simulated data and cross-cohort replication analyses. Best performing methods were subsequently applied to real-world human cohort data.

RESULTS

Our benchmarks showed that prior-based strategies outperform methods without prior information in simulated data and show better replication across datasets. Application of our approach to human cohort data highlighted two novel regulatory networks related to schizophrenia and lean body mass for which we generated novel functional hypotheses.

CONCLUSIONS

We demonstrate that existing biological knowledge can improve the integrative analysis of networks underlying trans associations and generate novel hypotheses about regulatory mechanisms.

Signaling and Gene Expression in Skeletal Muscles in Type 2 Diabetes: Current Results and OMICS Perspectives

Skeletal muscles mainly contribute to the emergence of insulin resistance, impaired glucose tolerance and the development of type 2 diabetes. Molecular mechanisms that regulate glucose uptake are diverse, including the insulin-dependent as most important, and others as also significant. They involve a wide range of proteins that control intracellular traffic and exposure of glucose transporters on the cell surface to create an extensive regulatory network. Here, we highlight advantages of the omics approaches to explore the insulin-regulated proteins and genes in human skeletal muscle with varying degrees of metabolic disorders. We discuss methodological aspects of the assessment of metabolic dysregulation and molecular responses of human skeletal muscle to insulin. The known molecular mechanisms of glucose uptake regulation and the first results of phosphoproteomic and transcriptomic studies are reviewed, which unveiled a large-scale array of insulin targets in muscle cells. They demonstrate that a clear depiction of changes that occur during metabolic dysfunction requires systemic and combined analysis at different levels of regulation, including signaling pathways, transcription factors, and gene expression. Such analysis seems promising to explore yet undescribed regulatory mechanisms of glucose uptake by skeletal muscle and identify the key regulators as potential therapeutic targets.

ACE2 expression in adipose tissue is associated with cardio-metabolic risk factors and cell type composition-implications for COVID-19

BACKGROUND

COVID-19 severity varies widely. Although some demographic and cardio-metabolic factors, including age and obesity, are associated with increasing risk of severe illness, the underlying mechanism(s) are uncertain.

SUBJECTS/METHODS

In a meta-analysis of three independent studies of 1471 participants in total, we investigated phenotypic and genetic factors associated with subcutaneous adipose tissue expression of Angiotensin I Converting Enzyme 2 (ACE2), measured by RNA-Seq, which acts as a receptor for SARS-CoV-2 cellular entry.

RESULTS

Lower adipose tissue ACE2 expression was associated with multiple adverse cardio-metabolic health indices, including type 2 diabetes (T2D) (P = 9.14 × 10-6), obesity status (P = 4.81 × 10-5), higher serum fasting insulin (P = 5.32 × 10-4), BMI (P = 3.94 × 10-4), and lower serum HDL levels (P = 1.92 × 10-7). ACE2 expression was also associated with estimated proportions of cell types in adipose tissue: lower expression was associated with a lower proportion of microvascular endothelial cells (P = 4.25 × 10-4) and higher proportion of macrophages (P = 2.74 × 10-5). Despite an estimated heritability of 32%, we did not identify any proximal or distal expression quantitative trait loci (eQTLs) associated with adipose tissue ACE2 expression.

CONCLUSIONS

Our results demonstrate that individuals with cardio-metabolic features known to increase risk of severe COVID-19 have lower background ACE2 levels in this highly relevant tissue. Reduced adipose tissue ACE2 expression may contribute to the pathophysiology of cardio-metabolic diseases, as well as the associated increased risk of severe COVID-19.

Overview of Transcriptomic Research on Type 2 Diabetes: Challenges and Perspectives

Type 2 diabetes (T2D) is a common chronic disease whose etiology is known to have a strong genetic component. Standard genetic approaches, although allowing for the detection of a number of gene variants associated with the disease as well as differentially expressed genes, cannot fully explain the hereditary factor in T2D. The explosive growth in the genomic sequencing technologies over the last decades provided an exceptional impetus for transcriptomic studies and new approaches to gene expression measurement, such as RNA-sequencing (RNA-seq) and single-cell technologies. The transcriptomic analysis has the potential to find new biomarkers to identify risk groups for developing T2D and its microvascular and macrovascular complications, which will significantly affect the strategies for early diagnosis, treatment, and preventing the development of complications. In this article, we focused on transcriptomic studies conducted using expression arrays, RNA-seq, and single-cell sequencing to highlight recent findings related to T2D and challenges associated with transcriptome experiments.

Grainyhead 1 acts as a drug-inducible conserved transcriptional regulator linked to insulin signaling and lifespan

Aging is impacted by interventions across species, often converging on metabolic pathways. Transcription factors regulate longevity yet approaches for their pharmacological modulation to exert geroprotection remain sparse. We show that increased expression of the transcription factor Grainyhead 1 (GRH-1) promotes lifespan and pathogen resistance in Caenorhabditis elegans. A compound screen identifies FDA-approved drugs able to activate human GRHL1 and promote nematodal GRH-1-dependent longevity. GRHL1 activity is regulated by post-translational lysine methylation and the phosphoinositide (PI) 3-kinase C2A. Consistently, nematodal longevity following impairment of the PI 3-kinase or insulin/IGF-1 receptor requires grh-1. In BXD mice, Grhl1 expression is positively correlated with lifespan and insulin sensitivity. In humans, GRHL1 expression positively correlates with insulin receptor signaling and also with lifespan. Fasting blood glucose levels, including in individuals with type 2 diabetes, are negatively correlated with GRHL1 expression. Thereby, GRH-1/GRHL1 is identified as a pharmacologically malleable transcription factor impacting insulin signaling and lifespan.

Life- and healthspan of organisms can be modulated by dietary, genetic, or pharmacological interventions, which often affect metabolic pathways. Here the authors report that Grainyhead 1 is an evolutionarily conserved, drug-inducible transcription factor that promotes longevity in C. elegans, and thus a potential target for the development of geroprotective drugs.

Cis-regulation of antisense noncoding RNA at the JAZF1 locus in type 2 diabetes

BACKGROUND

Several genomic loci of type 2 diabetes (T2D) nominated in genome-wide association studies (GWASs) have been suggested to regulate metabolism in muscle. However, a large portion of the genetic risk and the underlying regulation remain unexplained. This study aimed to localize the potentially functional regions or genes at juxtaposed with another zinc finger protein 1 (JAZF1) locus and interpret their possible biological mechanisms in the muscle of T2D.

METHODS AND RESULTS

With a cross-population meta-analysis of 7 GWASs, we identified a linkage disequilibrium (LD) block within intron 1 of JAZF1 that was significantly associated with T2D (FDR < 0.05). The colocalization analysis showed a significant association between genetically determined expression of JAZF1 in skeletal muscle and T2D with a strong probability of colocalization (PP4=75.09%). This region also encodes the upstream regulatory region (URR) of the antisense noncoding RNA JAZF1-AS1. Expression-QTL (e-QTL) analysis detected a regulatory SNP within this LD block, rs864745, that is associated with the expression of JAZF1-AS1 and JAZF1. With in vitro cloning, we further reported the role of JAZF1-AS1 in cis-regulating JAZF1 by directly forming RNA double strands. Downregulation of JAZF1, caused by JAZF1-AS1 depletion, inhibited the glucose uptake and lipid oxidation in skeletal muscle.

CONCLUSIONS

This study proposes a strategy to identify a novel T2D gene at the reported locus and generated a model in which polymorphisms at JAZF1 influence T2D risk through antisense-mediated gene regulation.

Study of Water- and Organic-Soluble Extracts from Trichosanthes on Type 1 Diabetes Mellitus

This study investigates the effects of the water-soluble and organic-soluble Trichosanthes extracts on the hyperglycemic condition in streptozotocin- (STZ-) induced diabetic rats. The blood glucose levels, body weights, water intake, and urine volumes of rats in different experimental groups were monitored throughout the experiment, and the results obtained indicate that the two extracts can effectively reduce blood sugar levels, increase body weights, and improve water intake and urine volumes in diabetic rats. Based on blood biochemical analyses, the two extracts play an important role in regulating the diabetes-induced lipid metabolism disorder, increasing the levels of insulin and C-peptide, and alleviating the symptoms of diabetes. The variation in the liver glycogen contents of the water-soluble fraction and ethanol fraction groups suggests that the mechanisms underlying the hypoglycemic effects of the two extracts are different. Indeed, the water-soluble fraction alleviates diabetes symptoms in rats mainly by antioxidative activity, unlike the ethanol fraction.

Human and rat skeletal muscle single-nuclei multi-omic integrative analyses nominate causal cell types, regulatory elements, and SNPs for complex traits

Skeletal muscle accounts for the largest proportion of human body mass, on average, and is a key tissue in complex diseases and mobility. It is composed of several different cell and muscle fiber types. Here, we optimize single-nucleus ATAC-seq (snATAC-seq) to map skeletal muscle cell-specific chromatin accessibility landscapes in frozen human and rat samples, and single-nucleus RNA-seq (snRNA-seq) to map cell-specific transcriptomes in human. We additionally perform multi-omics profiling (gene expression and chromatin accessibility) on human and rat muscle samples. We capture type I and type II muscle fiber signatures, which are generally missed by existing single-cell RNA-seq methods. We perform cross-modality and cross-species integrative analyses on 33,862 nuclei and identify seven cell types ranging in abundance from 59.6% to 1.0% of all nuclei. We introduce a regression-based approach to infer cell types by comparing transcription start site-distal ATAC-seq peaks to reference enhancer maps and show consistency with RNA-based marker gene cell type assignments. We find heterogeneity in enrichment of genetic variants linked to complex phenotypes from the UK Biobank and diabetes genome-wide association studies in cell-specific ATAC-seq peaks, with the most striking enrichment patterns in muscle mesenchymal stem cells (∼3.5% of nuclei). Finally, we overlay these chromatin accessibility maps on GWAS data to nominate causal cell types, SNPs, transcription factor motifs, and target genes for type 2 diabetes signals. These chromatin accessibility profiles for human and rat skeletal muscle cell types are a useful resource for nominating causal GWAS SNPs and cell types.

Differential analysis of chromatin accessibility and gene expression profiles identifies cis-regulatory elements in rat adipose and muscle

Chromatin accessibility is a key factor influencing gene expression. We optimized the Omni-ATAC-seq protocol and used it together with RNA-seq to investigate cis-regulatory elements in rat white adipose and skeletal muscle, two tissues with contrasting metabolic functions. While promoter accessibility correlated with RNA expression, integration of the two datasets identified tissue-specific differentially accessible regions (DARs) that predominantly localized in intergenic and intron regions. DARs were mapped to differentially expressed (DE) genes enriched in distinct biological processes in each tissue. Randomly selected DE genes were validated by qPCR. Top enriched motifs in DARs predicted binding sites for transcription factors (TFs) showing tissue-specific up-regulation. The correlation between differential chromatin accessibility at a given TF binding motif and differential expression of target genes further supported the functional relevance of that motif. Our study identified cis-regulatory regions that likely play a major role in the regulation of tissue-specific gene expression in adipose and muscle.

Chromatin accessibility and gene expression during adipocyte differentiation identify context-dependent effects at cardiometabolic GWAS loci

Chromatin accessibility and gene expression in relevant cell contexts can guide identification of regulatory elements and mechanisms at genome-wide association study (GWAS) loci. To identify regulatory elements that display differential activity across adipocyte differentiation, we performed ATAC-seq and RNA-seq in a human cell model of preadipocytes and adipocytes at days 4 and 14 of differentiation. For comparison, we created a consensus map of ATAC-seq peaks in 11 human subcutaneous adipose tissue samples. We identified 58,387 context-dependent chromatin accessibility peaks and 3,090 context-dependent genes between all timepoint comparisons (log2 fold change>1, FDR<5%) with 15,919 adipocyte- and 18,244 preadipocyte-dependent peaks. Adipocyte-dependent peaks showed increased overlap (60.1%) with Roadmap Epigenomics adipocyte nuclei enhancers compared to preadipocyte-dependent peaks (11.5%). We linked context-dependent peaks to genes based on adipocyte promoter capture Hi-C data, overlap with adipose eQTL variants, and context-dependent gene expression. Of 16,167 context-dependent peaks linked to a gene, 5,145 were linked by two or more strategies to 1,670 genes. Among GWAS loci for cardiometabolic traits, adipocyte-dependent peaks, but not preadipocyte-dependent peaks, showed significant enrichment (LD score regression P<0.005) for waist-to-hip ratio and modest enrichment (P < 0.05) for HDL-cholesterol. We identified 659 peaks linked to 503 genes by two or more approaches and overlapping a GWAS signal, suggesting a regulatory mechanism at these loci. To identify variants that may alter chromatin accessibility between timepoints, we identified 582 variants in 454 context-dependent peaks that demonstrated allelic imbalance in accessibility (FDR<5%), of which 55 peaks also overlapped GWAS variants. At one GWAS locus for palmitoleic acid, rs603424 was located in an adipocyte-dependent peak linked to SCD and exhibited allelic differences in transcriptional activity in adipocytes (P = 0.003) but not preadipocytes (P = 0.09). These results demonstrate that context-dependent peaks and genes can guide discovery of regulatory variants at GWAS loci and aid identification of regulatory mechanisms.

Cardiovascular and metabolic diseases are widespread, and an increased understanding of genetic mechanisms behind these diseases could improve treatment. Chromatin accessibility and gene expression in relevant cell contexts can guide identification of regulatory elements and genetic mechanisms for disease traits. A relevant context for cardiovascular and metabolic disease traits is adipocyte differentiation. To identify regulatory elements and genes that display differences in activity during adipocyte differentiation, we profiled chromatin accessibility and gene expression in a human cell model of preadipocytes and adipocytes. We identified chromatin regions that change accessibility during differentiation and predicted genes they may affect. We also linked these chromatin regions to genetic variants associated with risk of disease. At one genomic region linked to fatty acids, a chromatin region more accessible in adipocytes linked to a fatty acid synthesis gene and exhibited allelic differences in transcriptional activity in adipocytes but not preadipocytes. These results demonstrate that chromatin regions and genes that change during cell context can guide discovery of regulatory variants and aid identification of disease mechanisms.

Multi-omics profiling: the way towards precision medicine in metabolic diseases

Metabolic diseases including type 2 diabetes mellitus (T2DM), non-alcoholic fatty liver disease (NAFLD), and metabolic syndrome (MetS) are alarming health burdens around the world, while therapies for these diseases are far from satisfying as their etiologies are not completely clear yet. T2DM, NAFLD, and MetS are all complex and multifactorial metabolic disorders based on the interactions between genetics and environment. Omics studies such as genetics, transcriptomics, epigenetics, proteomics, and metabolomics are all promising approaches in accurately characterizing these diseases. And the most effective treatments for individuals can be achieved via omics pathways, which is the theme of precision medicine. In this review, we summarized the multi-omics studies of T2DM, NAFLD, and MetS in recent years, provided a theoretical basis for their pathogenesis and the effective prevention and treatment, and highlighted the biomarkers and future strategies for precision medicine.

RNA sequencing reveals potential interacting networks between the altered transcriptome and ncRNome in the skeletal muscle of diabetic mice

For a global epidemic like Type 2 diabetes mellitus (T2DM), while impaired gene regulation is identified as a primary cause of aberrant cellular physiology; in the past few years, non-coding RNAs (ncRNAs) have emerged as important regulators of cellular metabolism. However, there are no reports of comprehensive in-depth cross-talk between these regulatory elements and the potential consequences in the skeletal muscle during diabetes. Here, using RNA sequencing, we identified 465 mRNAs and 12 long non-coding RNAs (lncRNAs), to be differentially regulated in the skeletal muscle of diabetic mice and pathway enrichment analysis of these altered transcripts revealed pathways of insulin, FOXO and AMP-activated protein kinase (AMPK) signaling to be majorly over-represented. Construction of networks showed that these pathways significantly interact with each other that might underlie aberrant skeletal muscle metabolism during diabetes. Gene-gene interaction network depicted strong interactions among several differentially expressed genes (DEGs) namely, Prkab2, Irs1, Pfkfb3, Socs2 etc. Seven altered lncRNAs depicted multiple interactions with the altered transcripts, suggesting possible regulatory roles of these lncRNAs. Inverse patterns of expression were observed between several of the deregulated microRNAs (miRNAs) and the differentially expressed transcripts in the tissues. Towards validation, overexpression of miR-381-3p and miR-539-5p in skeletal muscle C2C12 cells significantly decreased the transcript levels of their targets, Nfkbia, Pik3r1 and Pi3kr1, Cdkn2d, respectively. Collectively, the findings provide a comprehensive understanding of the interactions and cross-talk between the ncRNome and transcriptome in the skeletal muscle during diabetes and put forth potential therapeutic options for improving insulin sensitivity.

Genetic effects on liver chromatin accessibility identify disease regulatory variants

Identifying the molecular mechanisms by which genome-wide association study (GWAS) loci influence traits remains challenging. Chromatin accessibility quantitative trait loci (caQTLs) help identify GWAS loci that may alter GWAS traits by modulating chromatin structure, but caQTLs have been identified in a limited set of human tissues. Here we mapped caQTLs in human liver tissue in 20 liver samples and identified 3,123 caQTLs. The caQTL variants are enriched in liver tissue promoter and enhancer states and frequently disrupt binding motifs of transcription factors expressed in liver. We predicted target genes for 861 caQTL peaks using proximity, chromatin interactions, correlation with promoter accessibility or gene expression, and colocalization with expression QTLs. Using GWAS signals for 19 liver function and/or cardiometabolic traits, we identified 110 colocalized caQTLs and GWAS signals, 56 of which contained a predicted caPeak target gene. At the LITAF LDL-cholesterol GWAS locus, we validated that a caQTL variant showed allelic differences in protein binding and transcriptional activity. These caQTLs contribute to the epigenomic characterization of human liver and help identify molecular mechanisms and genes at GWAS loci.

Skeletal muscle proteomes reveal downregulation of mitochondrial proteins in transition from prediabetes into type 2 diabetes

Skeletal muscle insulin resistance is a central defect in the pathogenesis of type 2 diabetes (T2D). Here, we analyzed skeletal muscle proteome in 148 vastus lateralis muscle biopsies obtained from men covering all glucose tolerance phenotypes: normal, impaired fasting glucose (IFG), impaired glucose tolerance (IGT) and T2D. Skeletal muscle proteome was analyzed by a sequential window acquisition of all theoretical mass spectra (SWATH-MS) proteomics technique. Our data indicate a downregulation in several proteins involved in mitochondrial electron transport or respiratory chain complex assembly already in IFG and IGT muscles, with most profound decreases observed in T2D. Additional phosphoproteomic analysis reveals altered phosphorylation in several signaling pathways in IFG, IGT, and T2D muscles, including those regulating glucose metabolic processes, and the structure of muscle cells. These data reveal several alterations present in skeletal muscle already in prediabetes and highlight impaired mitochondrial energy metabolism in the trajectory from prediabetes into T2D.

•

Skeletal muscle proteome from men with all stages of glucose tolerance was analyzed

•

Phosphoproteomics reveal altered phosphorylation in IFG, IGT, and T2D muscles

•

OXPHOS proteins are decreased in prediabetic muscles, with most decrease in T2D

The trans-ancestral genomic architecture of glycemic traits

Glycemic traits are used to diagnose and monitor type 2 diabetes and cardiometabolic health. To date, most genetic studies of glycemic traits have focused on individuals of European ancestry. Here we aggregated genome-wide association studies comprising up to 281,416 individuals without diabetes (30% non-European ancestry) for whom fasting glucose, 2-h glucose after an oral glucose challenge, glycated hemoglobin and fasting insulin data were available. Trans-ancestry and single-ancestry meta-analyses identified 242 loci (99 novel; P < 5 × 10-8), 80% of which had no significant evidence of between-ancestry heterogeneity. Analyses restricted to individuals of European ancestry with equivalent sample size would have led to 24 fewer new loci. Compared with single-ancestry analyses, equivalent-sized trans-ancestry fine-mapping reduced the number of estimated variants in 99% credible sets by a median of 37.5%. Genomic-feature, gene-expression and gene-set analyses revealed distinct biological signatures for each trait, highlighting different underlying biological pathways. Our results increase our understanding of diabetes pathophysiology by using trans-ancestry studies for improved power and resolution.

Progress in Defining the Genetic Contribution to Type 2 Diabetes in Individuals of East Asian Ancestry

PURPOSE OF REVIEW

Prevalence of type 2 diabetes (T2D) and progression of complications differ between worldwide populations. While obesity is a major contributing risk factor, variations in physiological manifestations, e.g., developing T2D at lower body mass index in some populations, suggest other contributing factors. Early T2D genetic associations were mostly discovered in European ancestry populations. This review describes the progression of genetic discoveries associated with T2D in individuals of East Asian ancestry in the last 10 years and highlights the shared genetic susceptibility between the population groups and additional insights into genetic contributions to T2D.

RECENT FINDINGS

Through increased sample size and power, new genetic associations with T2D were discovered in East Asian ancestry populations, often with higher allele frequencies than European ancestry populations. As we continue to generate maps of T2D-associated variants across diverse populations, there will be a critical need to expand and diversify other omics resources to enable integration for clinical translation.

MicroRNA-129-5p inhibits C2C12 myogenesis and represses slow fiber gene expression in vitro

The miR-129 family is widely reported as tumor repressors, although their roles in skeletal muscle have not been fully investigated. Here, the function and mechanism of miR-129-5p in skeletal muscle, a member of the miR-129 family, were explored using C2C12 cell line. Our study showed that miR-129-5p was widely detected in mouse tissues, with the highest expression in skeletal muscle. Gain- and loss-of-function study showed that miR-129-5p could negatively regulate myogenic differentiation, indicated by reduced ratio of MyHC-positive myofibers and repressed expression of myogenic genes, such as MyoD, MyoG, and MyHC. Furthermore, miR-129-5p was more enriched in fast extensor digitorum longus (EDL) than in slow soleus (SOL). Enhanced miR-129-5p could significantly reduce the expression of mitochondrial cox family, together with that of MyHC I, and knockdown of miR-129-5p conversely increased the expression of cox genes and MyHC I. Mechanistically, miR-129-5p directly targeted the 3'-UTR of Mef2a, which was suppressed by miR-129-5p agomir at both mRNA and protein levels in C2C12 cells. Moreover, overexpression of Mef2a could rescue the inhibitory effects of miR-129-5p on the expression of myogenic factors and MyHC I. Collectively, our data revealed that miR-129-5p is a negative regulator of myogenic differentiation and slow fiber gene expression, thus affecting body metabolic homeostasis.

Molecular pathways behind acquired obesity: Adipose tissue and skeletal muscle multiomics in monozygotic twin pairs discordant for BMI

Tissue-specific mechanisms prompting obesity-related development complications in humans remain unclear. We apply multiomics analyses of subcutaneous adipose tissue and skeletal muscle to examine the effects of acquired obesity among 49 BMI-discordant monozygotic twin pairs. Overall, adipose tissue appears to be more affected by excess body weight than skeletal muscle. In heavier co-twins, we observe a transcriptional pattern of downregulated mitochondrial pathways in both tissues and upregulated inflammatory pathways in adipose tissue. In adipose tissue, heavier co-twins exhibit lower creatine levels; in skeletal muscle, glycolysis- and redox stress-related protein and metabolite levels remain higher. Furthermore, metabolomics analyses in both tissues reveal that several proinflammatory lipids are higher and six of the same lipid derivatives are lower in acquired obesity. Finally, in adipose tissue, but not in skeletal muscle, mitochondrial downregulation and upregulated inflammation are associated with a fatty liver, insulin resistance, and dyslipidemia, suggesting that adipose tissue dominates in acquired obesity.

Multiomics analyses of adipose tissue and skeletal muscle in BMI-discordant twins

Excess body weight downregulates mitochondrial pathways in both tissues

Excess body weight upregulates proinflammatory pathways in both tissues

Adipose tissue alterations are associated with metabolic health in acquired obesity

Chromatin information content landscapes inform transcription factor and DNA interactions

Interactions between transcription factors and chromatin are fundamental to genome organization and regulation and, ultimately, cell state. Here, we use information theory to measure signatures of organized chromatin resulting from transcription factor-chromatin interactions encoded in the patterns of the accessible genome, which we term chromatin information enrichment (CIE). We calculate CIE for hundreds of transcription factor motifs across human samples and identify two classes: low and high CIE. The 10-20% of common and tissue-specific high CIE transcription factor motifs, associate with higher protein-DNA residence time, including different binding site subclasses of the same transcription factor, increased nucleosome phasing, specific protein domains, and the genetic control of both chromatin accessibility and gene expression. These results show that variations in the information encoded in chromatin architecture reflect functional biological variation, with implications for cell state dynamics and memory.

Super enhancers-Functional cores under the 3D genome

Complex biochemical reactions take place in the nucleus all the time. Transcription machines must follow the rules. The chromatin state, especially the three-dimensional structure of the genome, plays an important role in gene regulation and expression. The super enhancers are important for defining cell identity in mammalian developmental processes and human diseases. It has been shown that the major components of transcriptional activation complexes are recruited by super enhancer to form phase-separated condensates. We summarize the current knowledge about super enhancer in the 3D genome. Furthermore, a new related transcriptional regulation model from super enhancer is outlined to explain its role in the mammalian cell progress.

A Multi-omic Integrative Scheme Characterizes Tissues of Action at Loci Associated with Type 2 Diabetes

Resolving the molecular processes that mediate genetic risk remains a challenge because most disease-associated variants are non-coding and functional characterization of these signals requires knowledge of the specific tissues and cell-types in which they operate. To address this challenge, we developed a framework for integrating tissue-specific gene expression and epigenomic maps to obtain "tissue-of-action" (TOA) scores for each association signal by systematically partitioning posterior probabilities from Bayesian fine-mapping. We applied this scheme to credible set variants for 380 association signals from a recent GWAS meta-analysis of type 2 diabetes (T2D) in Europeans. The resulting tissue profiles underscored a predominant role for pancreatic islets and, to a lesser extent, adipose and liver, particularly among signals with greater fine-mapping resolution. We incorporated resulting TOA scores into a rule-based classifier and validated the tissue assignments through comparison with data from cis-eQTL enrichment, functional fine-mapping, RNA co-expression, and patterns of physiological association. In addition to implicating signals with a single TOA, we found evidence for signals with shared effects in multiple tissues as well as distinct tissue profiles between independent signals within heterogeneous loci. Lastly, we demonstrated that TOA scores can be directly coupled with eQTL colocalization to further resolve effector transcripts at T2D signals. This framework guides mechanistic inference by directing functional validation studies to the most relevant tissues and can gain power as fine-mapping resolution and cell-specific annotations become richer. This method is generalizable to all complex traits with relevant annotation data and is made available as an R package.

Are Genome-Wide Association Study Identified Single-Nucleotide Polymorphisms Associated With Sprint Athletic Status? A Replication Study With 3 Different Cohorts

PURPOSE

To replicate previous genome-wide association study identified sprint-related polymorphisms in 3 different cohorts of top-level sprinters and to further validate the obtained results in functional studies.

METHODS

A total of 240 Japanese, 290 Russians, and 593 Brazilians were evaluated in a case-control approach. Of these, 267 were top-level sprint/power athletes. In addition, the relationship between selected polymorphisms and muscle fiber composition was evaluated in 203 Japanese and 287 Finnish individuals.

RESULTS

The G allele of the rs3213537 polymorphism was overrepresented in Japanese (odds ratio [OR]: 2.07, P = .024) and Russian (OR: 1.93, P = .027) sprinters compared with endurance athletes and was associated with an increased proportion of fast-twitch muscle fibers in Japanese (P = .02) and Finnish (P = .041) individuals. A meta-analysis of the data from 4 athlete cohorts confirmed that the presence of the G/G genotype rather than the G/A+A/A genotypes increased the OR of being a sprinter compared with controls (OR: 1.49, P = .01), endurance athletes (OR: 1.79, P = .001), or controls + endurance athletes (OR: 1.58, P = .002). Furthermore, male sprinters with the G/G genotype were found to have significantly faster personal times in the 100-m dash than those with G/A+A/A genotypes (10.50 [0.26] vs 10.76 [0.31], P = .014).

CONCLUSION

The rs3213537 polymorphism found in the CPNE5 gene was identified as a highly replicable variant associated with sprinting ability and the increased proportion of fast-twitch muscle fibers, in which the homozygous genotype for the major allele (ie, the G/G genotype) is preferable for performance.

CHD7 promotes neural progenitor differentiation in embryonic stem cells via altered chromatin accessibility and nascent gene expression

CHARGE syndrome, a rare multiple congenital anomaly condition, is caused by haploinsufficiency of the chromatin remodeling protein gene CHD7 (Chromodomain helicase DNA binding protein 7). Brain abnormalities and intellectual disability are commonly observed in individuals with CHARGE, and neuronal differentiation is reduced in CHARGE patient-derived iPSCs and conditional knockout mouse brains. However, the mechanisms of CHD7 function in nervous system development are not well understood. In this study, we asked whether CHD7 promotes gene transcription in neural progenitor cells via changes in chromatin accessibility. We used Chd7 null embryonic stem cells (ESCs) derived from Chd7 mutant mouse blastocysts as a tool to investigate roles of CHD7 in neuronal and glial differentiation. Loss of Chd7 significantly reduced neuronal and glial differentiation. Sholl analysis showed that loss of Chd7 impaired neuronal complexity and neurite length in differentiated neurons. Genome-wide studies demonstrated that loss of Chd7 leads to modified chromatin accessibility (ATAC-seq) and differential nascent expression (Bru-Seq) of neural-specific genes. These results suggest that CHD7 acts preferentially to alter chromatin accessibility of key genes during the transition of NPCs to neurons to promote differentiation. Our results form a basis for understanding the cell stage-specific roles for CHD7-mediated chromatin remodeling during cell lineage acquisition.

Unique Transcriptome Signature Distinguishes Patients With Heart Failure With Myopathy

Background People with chronic heart failure (CHF) experience severe skeletal muscle dysfunction, characterized by mitochondrial abnormalities, which exacerbates the primary symptom of exercise intolerance. However, the molecular triggers and characteristics underlying mitochondrial abnormalities caused by CHF remain poorly understood. Methods and Results We recruited 28 patients with CHF caused by reduced ejection fraction and 9 controls. We simultaneously biopsied skeletal muscle from the pectoralis major in the upper limb and from the vastus lateralis in the lower limb. We phenotyped mitochondrial function in permeabilized myofibers from both sites and followed this by complete RNA sequencing to identify novel molecular abnormalities in CHF skeletal muscle. Patients with CHF presented with upper and lower limb skeletal muscle impairments to mitochondrial function that were of a similar deficit and indicative of a myopathy. Mitochondrial abnormalities were strongly correlated to symptoms. Further RNA sequencing revealed a unique transcriptome signature in CHF skeletal muscle characterized by a novel triad of differentially expressed genes related to deficits in energy metabolism including adenosine monophosphate deaminase 3, pyridine nucleotide-disulphide oxidoreductase domain 2, and lactate dehydrogenase C. Conclusions Our data suggest an upper and lower limb metabolic myopathy that is characterized by a unique transcriptome signature in skeletal muscle of humans with CHF.

Unraveling the molecular heterogeneity in type 2 diabetes: a potential subtype discovery followed by metabolic modeling

Background

Type 2 diabetes mellitus (T2DM) is a complex multifactorial disease with a high prevalence worldwide. Insulin resistance and impaired insulin secretion are the two major abnormalities in the pathogenesis of T2DM. Skeletal muscle is responsible for over 75% of the glucose uptake and plays a critical role in T2DM. Here, we sought to provide a better understanding of the abnormalities in this tissue.

Methods

The muscle gene expression patterns were explored in healthy and newly diagnosed T2DM individuals using supervised and unsupervised classification approaches. Moreover, the potential of subtyping T2DM patients was evaluated based on the gene expression patterns.

Results

A machine-learning technique was applied to identify a set of genes whose expression patterns could discriminate diabetic subjects from healthy ones. A gene set comprising of 26 genes was found that was able to distinguish healthy from diabetic individuals with 94% accuracy. In addition, three distinct clusters of diabetic patients with different dysregulated genes and metabolic pathways were identified.

Conclusions

This study indicates that T2DM is triggered by different cellular/molecular mechanisms, and it can be categorized into different subtypes. Subtyping of T2DM patients in combination with their real clinical profiles will provide a better understanding of the abnormalities in each group and more effective therapeutic approaches in the future.

ACE2 expression in adipose tissue is associated with COVID-19 cardio-metabolic risk factors and cell type composition

COVID-19 severity has varied widely, with demographic and cardio-metabolic factors increasing risk of severe reactions to SARS-CoV-2 infection, but the underlying mechanisms for this remain uncertain. We investigated phenotypic and genetic factors associated with subcutaneous adipose tissue expression of Angiotensin I Converting Enzyme 2 ( ACE2 ), which has been shown to act as a receptor for SARS-CoV-2 cellular entry. In a meta-analysis of three independent studies including up to 1,471 participants, lower adipose tissue ACE2 expression was associated with adverse cardio-metabolic health indices including type 2 diabetes (T2D) and obesity status, higher serum fasting insulin and BMI, and lower serum HDL levels (P<5.32x10 -4 ). ACE2 expression levels were also associated with estimated proportions of cell types in adipose tissue; lower ACE2 expression was associated with a lower proportion of microvascular endothelial cells (P=4.25x10 -4 ) and higher macrophage proportion (P=2.74x10 -5 ), suggesting a link to inflammation. Despite an estimated heritability of 32%, we did not identify any proximal or distal genetic variants (eQTLs) associated with adipose tissue ACE2 expression. Our results demonstrate that at-risk individuals have lower background ACE2 levels in this highly relevant tissue. Further studies will be required to establish how this may contribute to increased COVID-19 severity.

ATAC-pipe: general analysis of genome-wide chromatin accessibility

Assay of Transposase-Accessible Chromatin by deep sequencing (ATAC-seq) has been widely used to profile the chromatin accessibility genome-wide. For the absence of an integrated scheme for deep data mining of specific biological issues, here we present ATAC-pipe, an efficient pipeline for general analysis of chromatin accessibility data obtained from ATAC-seq experiments. ATAC-pipe captures information includes not only the quality of original data and genome-wide chromatin accessibility but also signatures of significant differential peaks, transcription factor (TF) occupancy and nucleosome positions around regulatory sites. In addition, ATAC-pipe automatically converts statistic results into intuitive plots at publication quality, such as the read length distribution, heatmaps of sample clustering and cell-type-specific regulatory elements, enriched TF occupancy with motifs footprints and TF-driven regulatory networks. ATAC-pipe provides convenient workflow for researchers to study chromatin accessibility and gene regulation. Availability https://github.com/QuKunLab/ATAC-pipe.

A Long Non-coding RNA, LOC157273, Is an Effector Transcript at the Chromosome 8p23.1-PPP1R3B Metabolic Traits and Type 2 Diabetes Risk Locus

AIMS

Causal transcripts at genomic loci associated with type 2 diabetes (T2D) are mostly unknown. The chr8p23.1 variant rs4841132, associated with an insulin-resistant diabetes risk phenotype, lies in the second exon of a long non-coding RNA (lncRNA) gene, LOC157273, located 175 kilobases from PPP1R3B, which encodes a key protein regulating insulin-mediated hepatic glycogen storage in humans. We hypothesized that LOC157273 regulates expression of PPP1R3B in human hepatocytes.

METHODS

We tested our hypothesis using Stellaris fluorescent in situ hybridization to assess subcellular localization of LOC157273; small interfering RNA (siRNA) knockdown of LOC157273, followed by RT-PCR to quantify LOC157273 and PPP1R3B expression; RNA-seq to quantify the whole-transcriptome gene expression response to LOC157273 knockdown; and an insulin-stimulated assay to measure hepatocyte glycogen deposition before and after knockdown.

RESULTS

We found that siRNA knockdown decreased LOC157273 transcript levels by approximately 80%, increased PPP1R3B mRNA levels by 1.7-fold, and increased glycogen deposition by >50% in primary human hepatocytes. An A/G heterozygous carrier (vs. three G/G carriers) had reduced LOC157273 abundance due to reduced transcription of the A allele and increased PPP1R3B expression and glycogen deposition.

CONCLUSION

We show that the lncRNA LOC157273 is a negative regulator of PPP1R3B expression and glycogen deposition in human hepatocytes and a causal transcript at an insulin-resistant T2D risk locus.

Identification of type 2 diabetes loci in 433,540 East Asian individuals

Meta-analyses of genome-wide association studies (GWAS) have identified more than 240 loci that are associated with type 2 diabetes (T2D)1,2; however, most of these loci have been identified in analyses of individuals with European ancestry. Here, to examine T2D risk in East Asian individuals, we carried out a meta-analysis of GWAS data from 77,418 individuals with T2D and 356,122 healthy control individuals. In the main analysis, we identified 301 distinct association signals at 183 loci, and across T2D association models with and without consideration of body mass index and sex, we identified 61 loci that are newly implicated in predisposition to T2D. Common variants associated with T2D in both East Asian and European populations exhibited strongly correlated effect sizes. Previously undescribed associations include signals in or near GDAP1, PTF1A, SIX3, ALDH2, a microRNA cluster, and genes that affect the differentiation of muscle and adipose cells3. At another locus, expression quantitative trait loci at two overlapping T2D signals affect two genes-NKX6-3 and ANK1-in different tissues4-6. Association studies in diverse populations identify additional loci and elucidate disease-associated genes, biology, and pathways.

Multiomic Profiling Identifies cis-Regulatory Networks Underlying Human Pancreatic β Cell Identity and Function

EndoC-βH1 is emerging as a critical human β cell model to study the genetic and environmental etiologies of β cell (dys)function and diabetes. Comprehensive knowledge of its molecular landscape is lacking, yet required, for effective use of this model. Here, we report chromosomal (spectral karyotyping), genetic (genotyping), epigenomic (ChIP-seq and ATAC-seq), chromatin interaction (Hi-C and Pol2 ChIA-PET), and transcriptomic (RNA-seq and miRNA-seq) maps of EndoC-βH1. Analyses of these maps define known (e.g., PDX1 and ISL1) and putative (e.g., PCSK1 and mir-375) β cell-specific transcriptional cis-regulatory networks and identify allelic effects on cis-regulatory element use. Importantly, comparison with maps generated in primary human islets and/or β cells indicates preservation of chromatin looping but also highlights chromosomal aberrations and fetal genomic signatures in EndoC-βH1. Together, these maps, and a web application we created for their exploration, provide important tools for the design of experiments to probe and manipulate the genetic programs governing β cell identity and (dys)function in diabetes.

EndoC-βH1 is becoming an important cellular model to study genes and pathways governing human β cell identity and function, but its (epi)genomic similarity to primary human islets is unknown. Lawlor et al. complete and compare extensive EndoC and primary human islet multiomic maps to identify shared and distinct genomic circuitry.

Genome-wide chromatin accessibility and transcriptome profiling show minimal epigenome changes and coordinated transcriptional dysregulation of hedgehog signaling in Danforth's short tail mice

Danforth's short tail (Sd) mice provide an excellent model for investigating the underlying etiology of human caudal birth defects, which affect 1 in 10 000 live births. Sd animals exhibit aberrant axial skeleton, urogenital and gastrointestinal development similar to human caudal malformation syndromes including urorectal septum malformation, caudal regression, vertebral-anal-cardiac-tracheo-esophageal fistula-renal-limb (VACTERL) association and persistent cloaca. Previous studies have shown that the Sd mutation results from an endogenous retroviral (ERV) insertion upstream of the Ptf1a gene resulting in its ectopic expression at E9.5. Though the genetic lesion has been determined, the resulting epigenomic and transcriptomic changes driving the phenotype have not been investigated. Here, we performed ATAC-seq experiments on isolated E9.5 tailbud tissue, which revealed minimal changes in chromatin accessibility in Sd/Sd mutant embryos. Interestingly, chromatin changes were localized to a small interval adjacent to the Sd ERV insertion overlapping a known Ptf1a enhancer region, which is conserved in mice and humans. Furthermore, mRNA-seq experiments revealed increased transcription of Ptf1a target genes and, importantly, downregulation of hedgehog pathway genes. Reduced sonic hedgehog (SHH) signaling was confirmed by in situ hybridization and immunofluorescence suggesting that the Sd phenotype results, in part, from downregulated SHH signaling. Taken together, these data demonstrate substantial transcriptome changes in the Sd mouse, and indicate that the effect of the ERV insertion on Ptf1a expression may be mediated by increased chromatin accessibility at a conserved Ptf1a enhancer. We propose that human caudal dysgenesis disorders may result from dysregulation of hedgehog signaling pathways.