Impacts of Alternative Splicing Events on the Differentiation of Adipocytes

Transcriptomic analysis reveals regulation of adipogenesis via long non-coding RNA, alternative splicing, and alternative polyadenylation

Obesity is characterized by dysregulated adipogenesis that leads to increased number and/or size of adipocytes. Understanding the molecular mechanisms governing adipogenesis is therefore key to designing therapeutic interventions against obesity. In our study, we analyzed 3'-end sequencing data that we generated from human preadipocytes and adipocytes, as well as previously published RNA-seq datasets, to elucidate mechanisms of regulation via long non-coding RNA (lncRNA), alternative splicing (AS) and alternative polyadenylation (APA). We discovered lncRNAs that have not been previously characterized but may be key regulators of white adipogenesis. We also detected 100 AS events and, using motif enrichment analysis, identified RNA binding proteins (RBPs) that could mediate exon skipping-the most prevalent AS event. In addition, we show that usage of alternative poly(A) sites in introns or 3'-UTRs of key adipogenesis genes leads to isoform diversity, which can have significant biological consequences on differentiation efficiency. We also identified RBPs that may modulate APA and defined how 3'-UTR APA can regulate gene expression through gain or loss of specific microRNA binding sites. Taken together, our bioinformatics-based analysis reveals potential therapeutic avenues for obesity through manipulation of lncRNA levels and the profile of mRNA isoforms via alternative splicing and polyadenylation.

Comparative Genome-Wide Alternative Splicing Analysis between Preadipocytes and Adipocytes

Alternative splicing is a ubiquitous regulatory mechanism in gene expression that allows a single gene to generate multiple messenger RNAs (mRNAs). Adipocyte development is regulated by many processes, and recent studies have found that splicing factors also play an important role in adipogenic development. In the present study, we further investigated the differences in selective shearing during different periods of adipocyte differentiation. We identified five alternative splicing types including skipped exon, mutually exclusive exon, Alternative 5' splice site, Alternative 3' splice site, and Retained intron, with skipped exons being the most abundant type of selective shearing. In total, 641 differentially expressed selective shearing genes were obtained, enriched in 279 pathways, from which we selected and verified the accuracy of the sequencing results. Overall, RNA-seq revealed changes in the splicing and expression levels of these new candidate genes between precursor adipocytes and adipocytes, suggesting that they may be involved in adipocyte generation and differentiation.

TMT Labeling Reveals the Effects of Exercises on the Proteomic Characteristics of the Subcutaneous Adipose Tissue of Growing High-Fat-Diet-Fed Rats

Aim: Growing period is an important period for fat remodeling. High-fat diet and exercise are reasons for adipose tissue (AT) remodeling, but existing evidence is not enough. Therefore, the effects of moderate-intensity continuous training (MICT) and high-intensity interval training (HIIT) on the proteomic characteristics of the subcutaneous AT of growing rats on normal diet or high-fat diet (HFD) were determined. Methods: Four-week-old male Sprague-Dawley rats (n = 48) were subdivided into six groups: normal diet control group, normal diet-MICT group, normal diet-HIIT group, HFD control group, HFD-MICT group, and HFD-HIIT group. Rats in the training group ran on a treadmill 5 days a week for 8 weeks (MICT: 50 min at 60-70% VO_2max intensity; HIIT: 7 min of warm-up and recovery at 70% VO_2max intensity, 6 sets of 3 min of 30% VO_2max followed by 3 min 90% VO_2max). Following physical assessment, inguinal subcutaneous adipose tissue (sWAT) was collected for proteome analysis using tandem mass tag labeling. Results: MICT and HIIT attenuated body fat mass and lean body mass but did not affect weight gain. Proteomics revealed the impact of exercise on ribosome, spliceosome, and the pentose phosphate pathway. However, the effect was reversed on HFD and normal diet. The differentially expressed proteins (DEPs) affected by MICT were related to oxygen transport, ribosome, and spliceosome. In comparison, the DEPs affected by HIIT were related to oxygen transport, mitochondrial electron transport, and mitochondrion protein. In HFD, HIIT was more likely to cause changes in immune proteins than MICT. However, exercise did not seem to reverse the protein effects of HFD. Conclusion: The exercise stress response in the growing period was stronger but increased the energy metabolism and metabolism. MICT and HIIT can reduce fat, increase muscle percentage, and improve maximum oxygen uptake in rats fed with HFD. However, in rats with normal diet, MICT and HIIT triggered more immune responses of sWAT, especially HIIT. In addition, spliceosomes may be the key factors in AT remodeling triggered by exercise and diet.

Transcriptomes of Testes at Different Developmental Stages in the Opsariichthys bidens Predict Key Genes for Testis Development and Spermatogenesis

Testis development is a complex process involving multiple genes, and the molecular mechanisms underlying testis development in Opsariichthys bidens remain unclear. We performed transcriptome sequencing analysis on a total of 12 samples of testes from stages II, III, IV, and V of O. bidens and obtained a total of 79.52 Gb clean data, as well as 288,573 transcripts and 116,215 unigenes. Differential expression analysis showed that 22,857 differentially expressed genes (DEGs) were screened in six comparison groups (III vs. II, IV vs. II, V vs. II, IV vs. III, V vs. III, and V vs. IV). Kyoto Encyclopedia of Genes and Genomes enrichment analysis of DEGs showed that six comparison groups were significantly enriched for a total of 20 significantly up- or down-regulated pathways, including six pathways related to signal transduction, three pathways related to energy metabolism, five pathways related to disease, and two pathways related to ribosomes. Furthermore, our investigation revealed that DEGs were enriched in several important functional pathways, such as Huntington's disease signaling pathway, TGF-β signaling pathway, and ribosome signaling pathway. Protein-protein interaction network analysis of DEGs identified 63 up-regulated hub genes, including 9 kinesin genes and 2 cytoplasmic dynein genes, and 39 down-regulated hub genes, including 13 ribosomal protein genes. This result contributes to the knowledge of spermatogenesis and testis development in O. bidens.

The Role of Alternative Splicing in Cancer: Regulatory Mechanism, Therapeutic Strategy, and Bioinformatics Application

[Formula: see text] Alternative splicing (AS) can generate distinct transcripts and subsequent isoforms that play differential functions from the same pre-mRNA. Recently, increasing numbers of studies have emerged, unmasking the association between AS and cancer. In this review, we arranged AS events that are closely related to cancer progression and presented promising treatments based on AS for cancer therapy. Obtaining proliferative capacity, acquiring invasive properties, gaining angiogenic features, shifting metabolic ability, and getting immune escape inclination are all splicing events involved in biological processes. Spliceosome-targeted and antisense oligonucleotide technologies are two novel strategies that are hopeful in tumor therapy. In addition, bioinformatics applications based on AS were summarized for better prediction and elucidation of regulatory routines mingled in. Together, we aimed to provide a better understanding of complicated AS events associated with cancer biology and reveal AS a promising target of cancer treatment in the future.

Identification and characterization of alternative splicing variants of buffalo LXRα expressed in mammary gland

Liver X receptor α (LXRα) is a ligand-dependent transcription factor and plays an important role in the regulation of cholesterol homeostasis, fatty acid biosynthesis and glucose metabolism. In this study, transcripts of LXRα gene were cloned and characterized from buffalo mammary gland, and three alternative splicing transcripts of buffalo LXRα gene were identified, named LXRα1, LXRα2 and LXRα3. The structure of the LXRα transcripts of buffalo and cattle was highly similar. Bioinformatics analysis showed that LXRα1 contains two complete functional domains of LXRα, one is the DNA-binding domain (NR_DBD_LXR) and the other is the ligand-binding domain (NR_LBD_LXR). The reading frame of LXRα2 is altered due to the skipping of exon 9, which truncates its encoding protein prematurely at the 400th amino acid residue, making it contain a complete DNA-binding domain and part of a ligand-binding domain. Due to the deletion of exon 4, the protein encoded by LXRα3 lacks 89 amino acid residues and contains only a complete ligand-binding domain, which makes it lose its transcriptional regulation function. In addition, motifs and conserved domains of three LXRα variants of buffalo were highly consistent with those of corresponding transcripts from other mammal species. Subcellular localization analysis showed that LXRα1 plays a functional role in the nucleus of buffalo mammary epithelial cells, while LXRα2 and LXRα3 are distributed in the nucleus and cytoplasm. Compared with non-lactating period, the mRNA abundance of the three LXRα transcripts in the mammary gland tissue of buffalo increased during lactating period, revealing that they play a key role in the synthesis of buffalo milk fat. Among the three LXRα transcripts, LXRα1 has the highest expression in the mammary gland, indicating that it is the major transcript in the mammary gland and has important regulatory functions, while LXRα2 and LXRα3 may have regulatory effects on the function of LXRα1. This study highlights the key role of LXRα alternative splicing in the post-transcriptional regulation of buffalo lactation.

Comprehensive characterization of 11 prognostic alternative splicing events in ovarian cancer interacted with the immune microenvironment

Alternative splicing (AS) events play a crucial role in the tumorigenesis and progression of cancer. Transcriptome data and Percent Spliced In (PSI) values of ovarian cancer patients were downloaded from TCGA database and TCGA SpliceSeq. Totally we identified 1472 AS events that were associated with survival of ovarian serous cystadenocarcinoma (OC) and exon skipping (ES) was the most important type. Univariate and multivariate Cox regression analysis were performed to identify survival-associated AS events and developed the prognostic model based on 11-AS events. The immune cells and different response to cytotoxic T lymphocyte associated antigen 4 (CTLA-4) and programmed cell death protein 1 (PD-1) blockers in low-risk and high-risk group of OC patients were analyzed. Ten kinds of immune cells were found up-regulated in low-risk group. Activated B cell, natural killer T cell, natural killer cell and regulatory T cell were associated with survival of OC. The patients in low-risk group had good response to CTLA-4 and PD-1 blockers treatment. Moreover, a regulatory network was established according to the correlation between AS events and splicing factors (SFs). The present study provided valuable insights into the underlying mechanisms of OC. AS events that were correlated with the immune system might be potential therapeutic targets.

Impaired mRNA splicing and proteostasis in preadipocytes in obesity-related metabolic disease

Preadipocytes are crucial for healthy adipose tissue expansion. Preadipocyte differentiation is altered in obese individuals, which has been proposed to contribute to obesity-associated metabolic disturbances. Here, we aimed at identifying the pathogenic processes underlying impaired adipocyte differentiation in obese individuals with insulin resistance (IR)/type 2 diabetes (T2D). We report that down-regulation of a key member of the major spliceosome, PRFP8/PRP8, as observed in IR/T2D preadipocytes from subcutaneous (SC) fat, prevented adipogenesis by altering both the expression and splicing patterns of adipogenic transcription factors and lipid droplet-related proteins, while adipocyte differentiation was restored upon recovery of PRFP8/PRP8 normal levels. Adipocyte differentiation was also compromised under conditions of endoplasmic reticulum (ER)-associated protein degradation (ERAD) hyperactivation, as occurs in SC and omental (OM) preadipocytes in IR/T2D obesity. Thus, targeting mRNA splicing and ER proteostasis in preadipocytes could improve adipose tissue function and thus contribute to metabolic health in obese individuals.

Differential expression and alternative splicing of transcripts in orbital adipose/connective tissue of thyroid-associated ophthalmopathy

Thyroid-associated ophthalmopathy is a typical autoimmune disease of orbital tissues. Alternative splicing significantly influences many diseases progression, including cancer, age-related macular degeneration, and multiple sclerosis, by modulating the expression of transcripts. However, its role in thyroid-associated ophthalmopathy is still unclear. In this study, differential expression transcripts and differential alternative splicing genes in orbital adipose/connective tissues of thyroid-associated ophthalmopathy patients were detected using RNA sequencing, Cuffdiff, and replicate multivariate analysis of transcript splicing. Three thousand ninety six differential expression transcripts and 2355 differential alternative splicing genes were screened out, while functional enrichment analysis indicated that differential expression transcript and differential alternative splicing genes were associated with immune modulation, extracellular matrix remodeling, and adipogenesis. The expression of the SORBS1, SEPT2, COL12A1, and VCAN gene transcripts was verified by qRT-PCR. In conclusion, prevalent alternative splicing is involved in the disease development in thyroid-associated ophthalmopathy. More attention should be paid to the mechanism of alternative splicing to explore more potential therapeutic targets in thyroid-associated ophthalmopathy.

LncRNA Ctcflos orchestrates transcription and alternative splicing in thermogenic adipogenesis

The recruitment of thermogenic brite adipocytes within white adipose tissue attenuates obesity and metabolic comorbidities, arousing interest in understanding the underlying regulatory mechanisms. The molecular network of brite adipogenesis, however, remains largely unresolved. In this light, long noncoding RNAs (lncRNAs) emerged as a versatile class of modulators that control many steps within the differentiation machinery. Leveraging the naturally varying propensities of different inbred mouse strains for white adipose tissue browning, we identify the nuclear lncRNA Ctcflos as a pivotal orchestrator of thermogenic gene expression during brite adipocyte differentiation. Mechanistically, Ctcflos acts as a pleiotropic regulator, being essential for the transcriptional recruitment of the early core thermogenic regulatory program and the modulation of alternative splicing to drive brite adipogenesis. This is showcased by Ctcflos‐regulated gene transcription and splicing of the key browning factor Prdm16 toward the isoform that is specific for the thermogenic gene program. Conclusively, our findings emphasize the mechanistic versatility of lncRNAs acting at several independent levels of gene expression for effective regulation of key differentiation factors to direct cell fate and function.

An alternative splicing signature model for predicting hepatocellular carcinoma-specific survival

BACKGROUND

Alternative splicing (AS) is a transcriptional regulation mechanism, which can expand the coding ability of genome and contribute to the occurrence and development of cancer. A systematic analysis of AS in hepatocellular carcinoma (HCC) is lacking and urgently needed.

METHODS

Univariate and multivariate Cox regression analyses were used to distinguish survival-related AS events and to calculate the risk score. Kaplan-Meier analysis and receiver operating characteristic (ROC) curves were used to evaluate the AS events' clinical significance to build a risk model in HCC.

RESULTS

Data of AS events was obtained from the Splice-Seq database. The corresponding clinical information of HCC was downloaded from The Cancer Genome Atlas (TCGA) data portal. We analyzed 78,878 AS events from 13,045 genes in HCC patients. A total of 2,440 and 2,888 AS events were significantly related to HCC patients' disease-free survival (DFS) and overall survival (OS). The two prognostic models (DFS and OS) were constructed based on a total of seven AS types from survival-related AS events above. The area under the curve (AUC) of the ROC curves was 0.769 in the DFS cohort and 0.886 in the OS cohort.

CONCLUSIONS

The prognostic model constructed by AS events can be used to predict the prognosis of HCC patients and provide potential therapeutic targets for further validation.

Alternative splicing events during adipogenesis from hMSCs

Adipogenesis, the developmental process of progenitor-cell differentiating into adipocytes, leads to fat metabolic disorders. Alternative splicing (AS), a ubiquitous regulatory mechanism of gene expression, allows the generation of more than one unique messenger RNA (mRNA) species from a single gene. Till now, alternative splicing events during adipogenesis from human mesenchymal stem cells (hMSCs) are not yet fully elucidated. We performed RNA-Seq coupled with bioinformatics analysis to identify the differentially expressed AS genes and events during adipogenesis from hMSCs. A global survey separately identified 1262, 1181, 1167, and 1227 ASE involved in the most common types of AS including cassette exon, alt3, and alt5, especially with cassette exon the most prevalent, at 7, 14, 21, and 28 days during adipogenesis. Interestingly, 122 differentially expressed ASE referred to 118 genes, and the three genes including ACTN1 (alt3 and cassette), LRP1 (alt3 and alt5), and LTBP4 (cassette, cassette_multi, and unknown), appeared in multiple AS types of ASE during adipogenesis. Except for all the identified ASE of LRP1 occurred in the extracellular topological domain, alt3 (84) in transmembrane domain significantly differentially expressed was the potential key event during adipogenesis. Overall, we have, for the first time, conducted the global transcriptional profiling during adipogenesis of hMSCs to identify differentially expressed ASE and ASE-related genes. This finding would provide extensive ASE as the regulator of adipogenesis and the potential targets for future molecular research into adipogenesis-related metabolic disorders.

Distinct Patterns of PPARγ Promoter Usage, Lipid Degradation Activity, and Gene Expression in Subcutaneous Adipose Tissue of Lean and Obese Swine

Subcutaneous adipose tissue is a loose connective tissue specializing in the regulation of energy storage and metabolization. In domesticated pigs (Sus scrofa), the temporal development of subcutaneous adipose tissue is critical for meat production. However, the regulation of adipose tissue development remains unclear. Here, the subcutaneous adipose tissue development was characterized and compared in lean (Danish-Landrace) and obese (Min) pigs at juvenile and the juvenile-to-adult growth stages. Using RNA sequencing, we profiled the transcriptome of subcutaneous adipose tissue isolated from 4- and 16-week-old pigs and identified 24,718 expressed transcription units. Of them, 6327 genes were differentially expressed between the breeds and/or developmental stages. Compared with obese pigs, upregulated genes in lean pigs showed significant function and pathway enrichment in fatty acid degradation and mitochondrial functions. Further analysis uncovered the distinct usage preferences of the three alternative peroxisomeproliferator-activatedreceptorγ (PPARγ) promoters associated with the development of subcutaneous adipose tissue in both breeds. Transcriptome analysis of subcutaneous adipose tissue in lean and obese pigs suggested that marker-assisted selection of fatty acid degradation and PPARγ signaling pathways could be important directions for future pork quality improvement and modern breeding.

Elucidation of the 14-3-3ζ interactome reveals critical roles of RNA-splicing factors during adipogenesis

Adipogenesis involves a complex signaling network requiring strict temporal and spatial organization of effector molecules. Molecular scaffolds, such as 14-3-3 proteins, facilitate such organization, and we have previously identified 14-3-3ζ as an essential scaffold in adipocyte differentiation. The interactome of 14-3-3ζ is large and diverse, and it is possible that novel adipogenic factors may be present within it, but this possibility has not yet been tested. Herein, we generated mouse embryonic fibroblasts from mice overexpressing a tandem affinity purification (TAP) epitope-tagged 14-3-3ζ molecule. After inducing adipogenesis, TAP-14-3-3ζ complexes were purified, followed by MS analysis to determine the 14-3-3ζ interactome. We observed more than 100 proteins that were unique to adipocyte differentiation, 56 of which were novel interacting partners. Among these, we were able to identify previously established regulators of adipogenesis (i.e. Ptrf/Cavin1) within the 14-3-3ζ interactome, confirming the utility of this approach to detect adipogenic factors. We found that proteins related to RNA metabolism, processing, and splicing were enriched in the interactome. Analysis of transcriptomic data revealed that 14-3-3ζ depletion in 3T3-L1 cells affected alternative splicing of mRNA during adipocyte differentiation. siRNA-mediated depletion of RNA-splicing factors within the 14-3-3ζ interactome, that is, of Hnrpf, Hnrpk, Ddx6, and Sfpq, revealed that they have essential roles in adipogenesis and in the alternative splicing of Pparg and the adipogenesis-associated gene Lpin1 In summary, we have identified novel adipogenic factors within the 14-3-3ζ interactome. Further characterization of additional proteins within the 14-3-3ζ interactome may help identify novel targets to block obesity-associated expansion of adipose tissues.

Insulin Receptor Isoforms in Physiology and Disease: An Updated View

The insulin receptor (IR) gene undergoes differential splicing that generates two IR isoforms, IR-A and IR-B. The physiological roles of IR isoforms are incompletely understood and appear to be determined by their different binding affinities for insulin-like growth factors (IGFs), particularly for IGF-2. Predominant roles of IR-A in prenatal growth and development and of IR-B in metabolic regulation are well established. However, emerging evidence indicates that the differential expression of IR isoforms may also help explain the diversification of insulin and IGF signaling and actions in various organs and tissues by involving not only different ligand-binding affinities but also different membrane partitioning and trafficking and possibly different abilities to interact with a variety of molecular partners. Of note, dysregulation of the IR-A/IR-B ratio is associated with insulin resistance, aging, and increased proliferative activity of normal and neoplastic tissues and appears to sustain detrimental effects. This review discusses novel information that has generated remarkable progress in our understanding of the physiology of IR isoforms and their role in disease. We also focus on novel IR ligands and modulators that should now be considered as an important strategy for better and safer treatment of diabetes and cancer and possibly other IR-related diseases.

Global transcriptome analysis identifies differentially expressed genes related to lipid metabolism in Wagyu and Holstein cattle

Fat deposition of beef cattle varies between breeds. However, the regulation mechanism is still not elucidated completely at molecular level. In the present study, we comparatively analyzed transcriptome of subcutaneous adipose tissue between Wagyu and Holstein cattle with a significant difference in fat deposition to identify key genes associated with fat metabolism and adipogenesis by high-throughput RNA-seq technology. A total of 59,149,852 and 69,947,982 high quality reads were generated, respectively. With further analysis, 662 differentially expressed genes were identified. Gene Ontology and KEGG pathway analysis revealed that many differentially expressed genes were enriched in several biological processes and pathways relevant to adipogenesis and lipid metabolism, in which PPAR and fatty acid metabolism signaling pathways with related genes such as PPARγ, PLIN2 and ELOVL6 et al. play a critical role. Protein-protein interaction network analysis showed EGR1, FOS, SERPINE1, AGT, MMP2 may have great impact on adipocyte differentiation and adipogenesis. Moreover, potential alternative splicing events and single nucleotide polymorphisms (SNPs) were also identified. In summary, we comprehensively analyzed and discussed the transcriptome of subcutaneous adipose tissue of Wagyu and Holstein cattle, which might provide a theoretical basis for better understanding molecular mechanism of fat metabolism and deposition in beef cattle.

Alternative splicing switches: Important players in cell differentiation

Alternative splicing (AS) greatly expands the coding capacities of genomes by allowing the generation of multiple mature mRNAs from a limited number of genes. Although the massive switch in AS profiles that often accompanies variations in gene expression patterns occurring during cell differentiation has been characterized for a variety of models, their causes and mechanisms remain largely unknown. Here, we integrate foundational and recent studies indicating the AS switches that govern the processes of cell fate determination. We include some distinct AS events in pluripotent cells and somatic reprogramming and discuss new progresses on alternative isoform expression in adipogenesis, myogenic differentiation and stimulation of immune cells. Finally, we cover novel insights on AS mechanisms during neuronal differentiation, paying special attention to the role of chromatin structure.

ARN: Analysis and Visualization System for Adipogenic Regulation Network Information

Adipogenesis is the process of cell differentiation through which preadipocytes become adipocytes. Lots of research is currently ongoing to identify genes, including their gene products and microRNAs, that correlate with fat cell development. However, information fragmentation hampers the identification of key regulatory genes and pathways. Here, we present a database of literature-curated adipogenesis-related regulatory interactions, designated the Adipogenesis Regulation Network (ARN, http://210.27.80.93/arn/), which currently contains 3101 nodes (genes and microRNAs), 1863 regulatory interactions, and 33,969 expression records associated with adipogenesis, based on 1619 papers. A sentence-based text-mining approach was employed for efficient manual curation of regulatory interactions from approximately 37,000 PubMed abstracts. Additionally, we further determined 13,103 possible node relationships by searching miRGate, BioGRID, PAZAR and TRRUST. ARN also has several useful features: i) regulatory map information; ii) tests to examine the impact of a query node on adipogenesis; iii) tests for the interactions and modes of a query node; iv) prediction of interactions of a query node; and v) analysis of experimental data or the construction of hypotheses related to adipogenesis. In summary, ARN can store, retrieve and analyze adipogenesis-related information as well as support ongoing adipogenesis research and contribute to the discovery of key regulatory genes and pathways.