Alcohol use disorder causes global changes in splicing in the human brain

Alcohol use disorder (AUD) is a widespread disease leading to the deterioration of cognitive and other functions. Mechanisms by which alcohol affects the brain are not fully elucidated. Splicing constitutes a nuclear process of RNA maturation, which results in the formation of the transcriptome. We tested the hypothesis as to whether AUD impairs splicing in the superior frontal cortex (SFC), nucleus accumbens (NA), basolateral amygdala (BLA), and central nucleus of the amygdala (CNA). To evaluate splicing, bam files from STAR alignments were indexed with samtools for use by rMATS software. Computational analysis of affected pathways was performed using Gene Ontology Consortium, Gene Set Enrichment Analysis, and LncRNA Ontology databases. Surprisingly, AUD was associated with limited changes in the transcriptome: expression of 23 genes was altered in SFC, 14 in NA, 102 in BLA, and 57 in CNA. However, strikingly, mis-splicing in AUD was profound: 1421 mis-splicing events were detected in SFC, 394 in NA, 1317 in BLA, and 469 in CNA. To determine the mechanism of mis-splicing, we analyzed the elements of the spliceosome: small nuclear RNAs (snRNAs) and splicing factors. While snRNAs were not affected by alcohol, expression of splicing factor heat shock protein family A (Hsp70) member 6 (HSPA6) was drastically increased in SFC, BLA, and CNA. Also, AUD was accompanied by aberrant expression of long noncoding RNAs (lncRNAs) related to splicing. In summary, alcohol is associated with genome-wide changes in splicing in multiple human brain regions, likely due to dysregulation of splicing factor(s) and/or altered expression of splicing-related lncRNAs

Oscillatory cAMP signaling rapidly alters H3K4 methylation

receptors (GPCRs) alter H3K4 methylation via oscillatory intracellular cAMP. Activation of Gs-coupled receptors caused a rapid decrease of H3K4me3 by elevating cAMP, whereas stimulation of Gi-coupled receptors increased H3K4me3 by diminishing cAMP. H3K4me3 gradually recovered towards baseline levels after the removal of GPCR ligands, indicating that H3K4me3 oscillates in tandem with GPCR activation. cAMP increased intracellular labile Fe(II), the cofactor for histone demethylases, through a non-canonical cAMP target—Rap guanine nucleotide exchange factor-2 (RapGEF2), which subsequently enhanced endosome acidification and Fe(II) release from the endosome via vacuolar H+-ATPase assembly. Removing Fe(III) from the media blocked intracellular Fe(II) elevation after stimulation of Gs-coupled receptors. Iron chelators and inhibition of KDM5 demethylases abolished cAMP-mediated H3K4me3 demethylation. Taken together, these results suggest a novel function of cAMP signaling in modulating histone demethylation through labile Fe(II)

Targeted massively parallel sequencing of autism spectrum disorder-associated genes in a case control cohort reveals rare loss-of-function risk variants

BACKGROUND: Autism spectrum disorder (ASD) is highly heritable, yet genome-wide association studies (GWAS), copy number variation screens, and candidate gene association studies have found no single factor accounting for a large percentage of genetic risk. ASD trio exome sequencing studies have revealed genes with recurrent de novo loss-of-function variants as strong risk factors, but there are relatively few recurrently affected genes while as many as 1000 genes are predicted to play a role. As such, it is critical to identify the remaining rare and low-frequency variants contributing to ASD. METHODS: We have utilized an approach of prioritization of genes by GWAS and follow-up with massively parallel sequencing in a case-control cohort. Using a previously reported ASD noise reduction GWAS analyses, we prioritized 837 RefSeq genes for custom targeting and sequencing. We sequenced the coding regions of those genes in 2071 ASD cases and 904 controls of European white ancestry. We applied comprehensive annotation to identify single variants which could confer ASD risk and also gene-based association analysis to identify sets of rare variants associated with ASD. RESULTS: We identified a significant over-representation of rare loss-of-function variants in genes previously associated with ASD, including a de novo premature stop variant in the well-established ASD candidate gene RBFOX1. Furthermore, ASD cases were more likely to have two damaging missense variants in candidate genes than controls. Finally, gene-based rare variant association implicates genes functioning in excitatory neurotransmission and neurite outgrowth and guidance pathways including CACNAD2, KCNH7, and NRXN1. CONCLUSIONS: We find suggestive evidence that rare variants in synaptic genes are associated with ASD and that loss-of-function mutations in ASD candidate genes are a major risk factor, and we implicate damaging mutations in glutamate signaling receptors and neuronal adhesion and guidance molecules. Furthermore, the role of de novo mutations in ASD remains to be fully investigated as we identified the first reported protein-truncating variant in RBFOX1 in ASD. Overall, this work, combined with others in the field, suggests a convergence of genes and molecular pathways underlying ASD etiology. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s13229-015-0034-z) contains supplementary material, which is available to authorized users

Defective HNF4alpha-dependent gene expression as a driver of hepatocellular failure in alcoholic hepatitis

Alcoholic hepatitis (AH) is a life-threatening condition characterized by profound hepatocellular dysfunction for which targeted treatments are urgently needed. Identification of molecular drivers is hampered by the lack of suitable animal models. By performing RNA sequencing in livers from patients with different phenotypes of alcohol-related liver disease (ALD), we show that development of AH is characterized by defective activity of liver-enriched transcription factors (LETFs). TGFβ1 is a key upstream transcriptome regulator in AH and induces the use of HNF4α P2 promoter in hepatocytes, which results in defective metabolic and synthetic functions. Gene polymorphisms in LETFs including HNF4α are not associated with the development of AH. In contrast, epigenetic studies show that AH livers have profound changes in DNA methylation state and chromatin remodeling, affecting HNF4α-dependent gene expression. We conclude that targeting TGFβ1 and epigenetic drivers that modulate HNF4α-dependent gene expression could be beneficial to improve hepatocellular function in patients with AH

HeatmapGenerator: high performance RNAseq and microarray visualization software suite to examine differential gene expression levels using an R and C++ hybrid computational pipeline

BACKGROUND: The graphical visualization of gene expression data using heatmaps has become an integral component of modern-day medical research. Heatmaps are used extensively to plot quantitative differences in gene expression levels, such as those measured with RNAseq and microarray experiments, to provide qualitative large-scale views of the transcriptonomic landscape. Creating high-quality heatmaps is a computationally intensive task, often requiring considerable programming experience, particularly for customizing features to a specific dataset at hand. METHODS: Software to create publication-quality heatmaps is developed with the R programming language, C++ programming language, and OpenGL application programming interface (API) to create industry-grade high performance graphics. RESULTS: We create a graphical user interface (GUI) software package called HeatmapGenerator for Windows OS and Mac OS X as an intuitive, user-friendly alternative to researchers with minimal prior coding experience to allow them to create publication-quality heatmaps using R graphics without sacrificing their desired level of customization. The simplicity of HeatmapGenerator is that it only requires the user to upload a preformatted input file and download the publicly available R software language, among a few other operating system-specific requirements. Advanced features such as color, text labels, scaling, legend construction, and even database storage can be easily customized with no prior programming knowledge. CONCLUSION: We provide an intuitive and user-friendly software package, HeatmapGenerator, to create high-quality, customizable heatmaps generated using the high-resolution color graphics capabilities of R. The software is available for Microsoft Windows and Apple Mac OS X. HeatmapGenerator is released under the GNU General Public License and publicly available at: http://sourceforge.net/projects/heatmapgenerator/. The Mac OS X direct download is available at: http://sourceforge.net/projects/heatmapgenerator/files/HeatmapGenerator_MAC_OSX.tar.gz/download. The Windows OS direct download is available at: http://sourceforge.net/projects/heatmapgenerator/files/HeatmapGenerator_WINDOWS.zip/download

Role of nitric oxide-based immunotherapy in augmenting prostate cancer progression by targeting androgen receptor heterogeneity

e17537 Background: A significant proportion of men with Prostate Cancer (PCa) develop castration resistant prostate cancer (CRPC) and do not respond to hormonal agents that decrease androgens. In trying to understand the causes of androgen resistance that develop in CRPC, it is considered most relevant to study the role of Androgen receptor (AR) in the development and progression of PCa from androgen dependent to androgen independent state. Recent studies have highlighted the significance of tumor microenvironment (TME) in regulation of PCa progression in addition to AR. A key molecule in the regulation of TME interactions is nitric oxide (NO). We have shown in our recent study, the critical association of NO with the TME in CRPC. However, the effects of NO to modulate the progression of PCa to CRPC with respect to AR still remains largely unexplored. Methods: 22RV1, LNCaP, LNCaP APIPC (cells expressing no AR), and LNCaP shAR/pATK (cells expressing low AR), cells were used for the study. Cell proliferation was first assessed by MTT assay. The castrated SCID mice were grafted with 22RV1 cells and were treated with GSNO at the dosage of 10mg/kg/day IP. After treatment, animals were humanely sacrificing. Tumor RNA and proteins were analysed for markers that are important for PCa progression using qPCR, western blot and cytokine antibody array. Animal experiments were carried out in compliance with the IACUC of University of Miami. GraphPad Prism (GraphPad Software) was used for statistical analysis. Results: In addition to reducing the tumor burden, the expression of anti-inflammatory (M2) macrophages (CD206 and Arginase1) is decreased and that of the pro-inflammatory (M1) macrophage (iNOS) is increased in mice which received increased NO levels. Furthermore, to study the effects of NO on progression of PCa from androgen dependent to androgen independent stage, we characterized the LNCAP cell models with differential extent of AR knockdown (LNCaP, LNCaP shAR/pATK and LNCaP APIPC ) for the effects of increased NO levels. Results showed that NO had significant impact on cell proliferation on androgen dependent PCa cells however the effects were negligible in cells expressing low or no AR, suggesting that effects of NO on PCa cell proliferation are AR dependent. Conclusions: Our results suggest that during PCa progression, NO suppresses TAMs to target the TME in an AR dependent manner. Further studies are undergoing to establish the impacts of NO in PCa progression

Optimized functional annotation of ChIP-seq data [version 1; peer review: 3 approved with reservations]

Different ChIP-seq peak callers often produce different output results from the same input. Since different peak callers are known to produce differentially enriched peaks with a large variance in peak length distribution and total peak count, accurately annotating peak lists with their nearest genes can be an arduous process. Functional genomic annotation of histone modification ChIP-seq data can be a particularly challenging task, as chromatin marks that have inherently broad peaks with a diffuse range of signal enrichment (e.g., H3K9me1, H3K27me3) differ significantly from narrow peaks that exhibit a compact and localized enrichment pattern (e.g., H3K4me3, H3K9ac). In addition, varying degrees of tissue-dependent broadness of an epigenetic mark can make it difficult to accurately and reliably link sequencing data to biological function. Thus, there exists an unmet need to develop a software program that can precisely tailor the computational analysis of a ChIP-seq dataset to the specific peak coordinates of the data and its surrounding genomic features. geneXtendeR optimizes the functional annotation of ChIP-seq peaks by exploring relative differences in annotating ChIP-seq peak sets to variable-length gene bodies. In contrast to prior techniques, geneXtendeR considers peak annotations beyond just the closest gene, allowing users to investigate peak summary statistics for the first-closest gene, second-closest gene, ..., n th -closest gene whilst ranking the output according to biologically relevant events and iteratively comparing the fidelity of peak-to-gene overlap across a user-defined range of upstream and downstream extensions on the original boundaries of each gene's coordinates. We tested geneXtendeR on 547 human transcription factor ChIP-seq ENCODE datasets and 198 human histone modification ChIP-seq ENCODE datasets, providing the analysis results as case studies. The geneXtendeR R/Bioconductor package (including detailed introductory vignettes) is available under the GPL-3 Open Source license and is freely available to download from Bioconductor at: https://bioconductor.org/packages/geneXtendeR

Absence of both MGME1 and POLG EXO abolishes mtDNA whereas absence of either creates unique mtDNA duplications

Both POLG and MGME1 are needed for mitochondrial DNA (mtDNA) maintenance in animal cells. POLG, the primary replicative polymerase of the mitochondria, has an exonuclease activity (3’→5’) that corrects for the misincorporation of bases. MGME1 serves as an exonuclease (5’→3’), producing ligatable DNA ends. Although both have a critical role in mtDNA replication and elimination of linear fragments, these mechanisms are still not fully understood. Using digital PCR to evaluate and compare mtDNA integrity, we show that Mgme1 knock out (Mgme1 KK) tissue mtDNA is more fragmented than POLG exonuclease deficient “Mutator” (Polg MM) or WT tissue. In addition, next generation sequencing of mutant hearts showed abundant duplications in/nearby the D-loop region and unique 100bp duplications evenly spaced throughout the genome only in Mgme1 KK hearts. However, despite these unique mtDNA features at steady-state, we observed a similar delay in the degradation of mtDNA after an induced double strand DNA break in both Mgme1 KK and Polg MM models. Lastly, we characterized double mutant (Polg MM/Mgme1 KK) cells and show that mtDNA cannot be maintained without at least one of these enzymatic activities. We propose a model for the generation of these genomic abnormalities which suggests a role for MGME1 outside of nascent mtDNA end ligation. Our results highlight the role of MGME1 in and outside of the D-loop region during replication, support the involvement of MGME1 in dsDNA degradation and demonstrate that POLG EXO and MGME1 can partially compensate for each other in maintaining mtDNA