Modulator-Dependent RBPs Changes Alternative Splicing Outcomes in Kidney Cancer

Alternative splicing alterations can contribute to human disease. The ability of an RNA-binding protein to regulate alternative splicing outcomes can be modulated by a variety of genetic and epigenetic mechanisms. In this study, we use a computational framework to investigate the roles of certain genes, termed modulators, on changing RBPs’ effect on splicing regulation. A total of 1,040,254 modulator-mediated RBP-splicing interactions were identified, including 137 RBPs, 4,309 splicing events and 2,905 modulator candidates from TCGA-KIRC RNA sequencing data. Modulators function categories were defined according to the correlation changes between RBPs expression and their targets splicing outcomes. QKI, as one of the RBPs influencing the most splicing events, attracted our attention in this study: 2,014 changing triplets were identified, including 1,101 modulators and 187 splicing events. Pathway enrichment analysis showed that QKI splicing targets were enriched in tight junction pathway, endocytosis and MAPK signaling pathways, all of which are highly associated with cancer development and progression. This is the first instance of a comprehensive study on how alternative splicing outcomes changes are associated with different expression level of certain proteins, even though they were regulated by the same RBP. Our work may provide a novel view on understanding alternative splicing mechanisms in kidney cancer

Highly robust model of transcription regulator activity predicts breast cancer overall survival

Background: While several multigene signatures are available for predicting breast cancer prognosis, particularly in early stage disease, effective molecular indicators are needed, especially for triple-negative carcinomas, to improve treatments and predict diagnostic outcomes. The objective of this study was to identify transcriptional regulatory networks to better understand mechanisms giving rise to breast cancer development and to incorporate this information into a model for predicting clinical outcomes. Methods: Gene expression profiles from 1097 breast cancer patients were retrieved from The Cancer Genome Atlas (TCGA). Breast cancer-specific transcription regulatory information was identified by considering the binding site information from ENCODE and the top co-expressed targets in TCGA using a nonlinear approach. We then used this information to predict breast cancer patient survival outcome. Result: We built a multiple regulator-based prediction model for breast cancer. This model was validated in more than 5000 breast cancer patients from the Gene Expression Omnibus (GEO) databases. We demonstrated our regulator model was significantly associated with clinical stage and that cell cycle and DNA replication related pathways were significantly enriched in high regulator risk patients. Conclusion: Our findings demonstrate that transcriptional regulator activities can predict patient survival. This finding provides additional biological insights into the mechanisms of breast cancer progression

Sesn3 protects against diet‐induced nonalcoholic steatohepatitis in mice via suppression of the TGFβ signal transduction

Sesn3 belongs to the three‐member sestrin protein family. Sestrins have been implicated in anti‐oxidative stress, AMPK and mTOR signal transduction, and metabolic homeostasis. However, the role of Sesn3 in the development of nonalcoholic steatohepatitis (NASH) has not been previously studied. In this work, we generated Sesn3 whole‐body knockout and liver‐specific transgenic mice to investigate the hepatic function of Sesn3 in diet‐induced NASH. With only 4 weeks of dietary treatment, Sesn3 knockout mice developed severe NASH phenotype as characterized by hepatic steatosis, inflammation, and fibrosis. Strikingly, after 8‐week feeding with a NASH‐inducing diet, Sesn3 transgenic mice were largely protected against NASH development. Transcriptomic analysis revealed that multiple extracellular matrix related processes were upregulated including TGFβ signaling and collagen production. Further biochemical and cell biological analyses have illustrated a critical control of the TGFβ‐Smad pathway by Sesn3 at the TGFβ receptor and Smad3 levels. First, Sesn3 inhibits the TGFβ receptor through an interaction with Smad7; second, Sesn3 directly inhibits the Smad3 function through protein‐protein interaction and cytosolic retention

Aberrant gene expression induced by a high fat diet is linked to H3K9 acetylation in the promoter-proximal region

Non-alcoholic fatty liver disease (NAFLD) is the most common chronic liver disease, with an estimated global prevalence of 1 in 4 individuals. Aberrant transcriptional control of gene expression is central to the pathophysiology of metabolic diseases. However, the molecular mechanisms leading to gene dysregulation are not well understood. Histone modifications play important roles in the control of transcription. Acetylation of histone 3 at lysine 9 (H3K9ac) is associated with transcriptional activity and is implicated in transcript elongation by controlling RNA polymerase II (RNAPII) pause-release. Hence, changes in this histone modification may shed information on novel pathways linking transcription control and metabolic dysfunction. Here, we carried out genome-wide analysis of H3K9ac in the liver of mice fed a control or a high-fat diet (an animal model of NAFLD), and asked whether this histone mark associates with changes in gene expression. We found that over 70% of RNAPII peaks in promoter-proximal regions overlapped with H3K9ac, consistent with a role of H3K9ac in the regulation of transcription. When comparing high-fat with control diet, approximately 17% of the differentially expressed genes were associated with changes in H3K9ac in their promoters, showing a strong correlation between changes in H3K9ac signal and gene expression. Overall, our data indicate that in response to a high-fat diet, dysregulated gene expression of a subset of genes may be attributable to changes in transcription elongation driven by H3K9ac. Our results point at an added mechanism of gene regulation that may be important in the development of metabolic diseases

Single-cell RNA sequencing of intramedullary canal tissue to improve methods for studying fracture repair biology

The ability to study the bone microenvironment of failed fracture healing may lead to biomarkers for fracture nonunion. Herein the authors describe a technique for isolating individual cells suitable for single-cell RNA sequencing analyses from intramedullary canal tissue collected by reaming during surgery. The purpose was to detail challenges and solutions inherent to the collection and processing of intramedullary canal tissue samples. The authors then examined single-cell RNA sequencing data from fresh and reanimated samples to demonstrate the feasibility of this approach for prospective studies

A transcriptome analysis of mitten crab testes (Eriocheir sinensis)

The identification of expressed genes involved in sexual precocity of the mitten crab (Eriocheir sinensis) is critical for a better understanding of its reproductive development. To this end, we constructed a cDNA library from the rapid developmental stage of testis of E. sinensis and sequenced 3,388 randomly picked clones. After processing, 2,990 high-quality expressed sequence tags (ESTs) were clustered into 2,415 unigenes including 307 contigs and 2,108 singlets, which were then compared to the NCBI non-redundant (nr) protein and nucleotide (nt) database for annotation with Blastx and Blastn, respectively. After further analysis, 922 unigenes were obtained with concrete annotations and 30 unigenes were found to have functions possibly related to the process of reproduction in male crabs – six transcripts relevant to spermatogenesis (especially Cyclin K and RecA homolog DMC1), two transcripts involved in nuclear protein transformation, two heat-shock protein genes, eleven transcription factor genes (a series of zinc-finger proteins), and nine cytoskeleton protein-related genes. Our results, besides providing valuable information related to crustacean reproduction, can also serve as a base for future studies of reproductive and developmental biology

RNA Polymerase II Binding Patterns Reveal Genomic Regions Involved in MicroRNA Gene Regulation

MicroRNAs are small non-coding RNAs involved in post-transcriptional regulation of gene expression. Due to the poor annotation of primary microRNA (pri-microRNA) transcripts, the precise location of promoter regions driving expression of many microRNA genes is enigmatic. This deficiency hinders our understanding of microRNA-mediated regulatory networks. In this study, we develop a computational approach to identify the promoter region and transcription start site (TSS) of pri-microRNAs actively transcribed using genome-wide RNA Polymerase II (RPol II) binding patterns derived from ChIP-seq data. Based upon the assumption that the distribution of RPol II binding patterns around the TSS of microRNA and protein coding genes are similar, we designed a statistical model to mimic RPol II binding patterns around the TSS of highly expressed, well-annotated promoter regions of protein coding genes. We used this model to systematically scan the regions upstream of all intergenic microRNAs for RPol II binding patterns similar to those of TSS from protein coding genes. We validated our findings by examining the conservation, CpG content, and activating histone marks in the identified promoter regions. We applied our model to assess changes in microRNA transcription in steroid hormone-treated breast cancer cells. The results demonstrate many microRNA genes have lost hormone-dependent regulation in tamoxifen-resistant breast cancer cells. MicroRNA promoter identification based upon RPol II binding patterns provides important temporal and spatial measurements regarding the initiation of transcription, and therefore allows comparison of transcription activities between different conditions, such as normal and disease states

A study of alterations in DNA epigenetic modifications (5mC and 5hmC) and gene expression influenced by simulated microgravity in human lymphoblastoid cells

Cells alter their gene expression in response to exposure to various environmental changes. Epigenetic mechanisms such as DNA methylation are believed to regulate the alterations in gene expression patterns. In vitro and in vivo studies have documented changes in cellular proliferation, cytoskeletal remodeling, signal transduction, bone mineralization and immune deficiency under the influence of microgravity conditions experienced in space. However microgravity induced changes in the epigenome have not been well characterized. In this study we have used Next-generation Sequencing (NGS) to profile ground-based “simulated” microgravity induced changes on DNA methylation (5-methylcytosine or 5mC), hydroxymethylation (5-hydroxymethylcytosine or 5hmC), and simultaneous gene expression in cultured human lymphoblastoid cells. Our results indicate that simulated microgravity induced alterations in the methylome (~60% of the differentially methylated regions or DMRs are hypomethylated and ~92% of the differentially hydroxymethylated regions or DHMRs are hyperhydroxymethylated). Simulated microgravity also induced differential expression in 370 transcripts that were associated with crucial biological processes such as oxidative stress response, carbohydrate metabolism and regulation of transcription. While we were not able to obtain any global trend correlating the changes of methylation/ hydroxylation with gene expression, we have been able to profile the simulated microgravity induced changes of 5mC over some of the differentially expressed genes that includes five genes undergoing differential methylation over their promoters and twenty five genes undergoing differential methylation over their gene-bodies. To the best of our knowledge, this is the first NGS-based study to profile epigenomic patterns induced by short time exposure of simulated microgravity and we believe that our findings can be a valuable resource for future explorations

Integration of Alzheimer’s disease genetics and myeloid genomics identifies disease risk regulatory elements and genes

Genome-wide association studies (GWAS) have identified more than 40 loci associated with Alzheimer’s disease (AD), but the causal variants, regulatory elements, genes and pathways remain largely unknown, impeding a mechanistic understanding of AD pathogenesis. Previously, we showed that AD risk alleles are enriched in myeloid-specific epigenomic annotations. Here, we show that they are specifically enriched in active enhancers of monocytes, macrophages and microglia. We integrated AD GWAS with myeloid epigenomic and transcriptomic datasets using analytical approaches to link myeloid enhancer activity to target gene expression regulation and AD risk modification. We identify AD risk enhancers and nominate candidate causal genes among their likely targets (including AP4E1, AP4M1, APBB3, BIN1, MS4A4A, MS4A6A, PILRA, RABEP1, SPI1, TP53INP1, and ZYX) in twenty loci. Fine-mapping of these enhancers nominates candidate functional variants that likely modify AD risk by regulating gene expression in myeloid cells. In the MS4A locus we identified a single candidate functional variant and validated it in human induced pluripotent stem cell (hiPSC)-derived microglia and brain. Taken together, this study integrates AD GWAS with multiple myeloid genomic datasets to investigate the mechanisms of AD risk alleles and nominates candidate functional variants, regulatory elements and genes that likely modulate disease susceptibility