Search

RNA provides the framework for the assembly of some of the most intricate macromolecular complexes within the cell, including the spliceosome and the mature ribosome. The assembly of these complexes relies on the coordinated association of RNA with hundreds of trans-acting protein factors. While some of these trans-acting factors are RNA-binding proteins (RBPs), others are adaptor proteins, and others still, function as both. Defects in the assembly of these complexes results in a number of human pathologies including neurodegeneration and cancer. Here, we demonstrate that Silencing Defective 2 (SDE2) is both an RNA binding protein and also a trans-acting adaptor protein that functions to regulate RNA splicing and ribosome biogenesis. SDE2 depletion leads to widespread changes in alternative splicing, defects in ribosome biogenesis and ultimately complete loss of cell viability. Our data highlight SDE2 as a previously uncharacterized essential gene required for the assembly and maturation of the complexes that carry out two of the most fundamental processes in mammalian cells., Graphical Abstract Graphical AbstractSDE2 directly interacts with C/D snoRNAs and the methyltransferase FBL to promote ribosome biogenesis. SDE2 also interacts with the U2 snRNP complex and is required for efficient pre-mRNA splicing.

Our understanding of the molecular underpinnings of chronic traumatic encephalopathy (CTE) and its associated pathology in post-mortem brain is incomplete. Factors including years of play and genetic risk variants influence the extent of tau pathology associated with disease expression, but how these factors affect gene expression, and whether those effects are consistent across the development of disease, is unknown. To address these questions, we conducted an analysis of the largest mRNASeq whole-transcriptome dataset available to date. We examined the genes and biological processes associated with disease by comparing individuals with CTE with control individuals with a history of repetitive head impacts that lack CTE pathology. We then identified genes and biological processes associated with total years of play as a measure of exposure, amount of tau pathology present at time of death, and the presence of APOE and TMEM106B risk variants. Samples were stratified into low and high pathology groups based on extent of tau pathology and years of play to model early vs late changes in response to exposure, and the relative effects associated with these factors were compared between these groups. Substantial gene expression changes were associated with severe disease for most of these factors, primarily implicating diverse, highly increased neuroinflammatory and neuroimmune processes. In contrast, low exposure groups had many fewer genes and processes implicated and show striking differences for some factors when compared with severe disease. Specifically, gene expression associated with amount of tau pathology showed a nearly perfect inverse relationship when compared between these two groups. Together, these results suggest the early disease process may differ substantially from that observed in late stages, that total years of play and tau pathology influence disease expression differently, and that related pathology-modifying risk variants may do so via distinct biological pathways.

Single nucleotide variants (SNVs) located in transcriptional regulatory regions can result in gene expression changes that lead to adaptive or detrimental phenotypic outcomes. Here, we predict gain or loss of binding sites for 741 transcription factors (TFs) across the human genome. We calculated ‘gainability’ and ‘disruptability’ scores for each TF that represent the likelihood of binding sites being created or disrupted, respectively. We found that functional cis-eQTL SNVs are more likely to alter TF binding sites than rare SNVs in the human population. In addition, we show that cancer somatic mutations have different effects on TF binding sites from different TF families on a cancer-type basis. Finally, we discuss the relationship between these results and cancer mutational signatures. Altogether, we provide a blueprint to study the impact of SNVs derived from genetic variation or disease association on TF binding to gene regulatory regions.

Obesity, the accumulation of body fat to excess, may cause serious negative health effects, including increased risk of heart disease, type 2 diabetes, stroke and certain cancers. The biology of obesity is complex and not well understood, involving both environmental and genetic factors and affecting metabolic and endocrine mechanisms in tissues of the gut, adipose, and brain. Previous RNA sequencing studies have identified transcripts associated with obesity and body mass index in blood and fat, often using animal models, but RNA sequencing studies in human brain tissue related to obesity have not been previously undertaken. We conducted both large and small RNA sequencing of hypothalamus (207 samples) and nucleus accumbens (276 samples) from individuals defined as consistently obese (124 samples), consistently normal weight as controls (148 samples) or selected without respect to BMI and falling within neither case nor control definition (211 samples), based on longitudinal BMI measures. The samples were provided by three cohort studies with brain donation programs; the Framingham Heart Study (FHS), the Religious Orders Study (ROS) and the Rush Memory and Aging Project (MAP). For each brain region and large/small RNA sequencing set, differential expression of obesity, BMI, brain region and sex was performed. Analyses were done transcriptome-wide as well as with a priori defined sets of obesity or BMI-associated mRNAs and microRNAs (miRNAs). There are sixteen mRNAs and five microRNAs that are differentially expressed (adjusted p < 0.05) by obesity or BMI in these tissues, several of which were validated with qPCR data. The results include many that are BMI-associated, such as APOBR and CES1, as well as many associated with the immune system and some with addiction, such as the gene sets “cytokine signaling in immune system” and “opioid signaling”. In spite of the relatively large number of samples, our study was likely under-powered to detect other transcripts or miRNA with relevant but smaller effects.

Virtually all patients with BRAF-mutant melanoma develop resistance to MAPK inhibitors largely through non-mutational events1,2. Although the epigenetic landscape has been shown to be altered in therapy-resistant melanomas and other cancers3,4, a specific targetable epigenetic mechanism regulating treatment resistance has not been validated to date. Here we evaluate the CoREST repressor complex and the novel inhibitor, corin5, within the context of melanoma phenotype plasticity and therapeutic resistance in order to define epigenetic mechanisms underlying these processes. We find that CoREST is a critical mediator of the major distinct melanoma phenotypes and that corin treatment of melanoma cells leads to phenotype reprogramming. We further demonstrate that treatment of BRAF inhibitor (BRAFi)-resistant melanomas with corin leads to resensitization of tumor cells to BRAFi. Among the transcriptional targets of CoREST in melanoma are the dual-specificity phosphatases (DUSPs). DUSP1 is shown to be consistently downregulated in BRAFi-resistant melanomas which can be reversed by corin treatment, thereby leading to downstream inhibition of p38 MAPK activity and resensitization of resistant cells to targeted BRAFi therapies. These findings identify the CoREST repressor complex as a central mediator of melanoma phenotype plasticity and resistance to targeted therapy and suggest that CoREST inhibitors may prove beneficial to patients with BRAF-mutant melanomas who have acquired BRAFi-resistance.

Previous studies on Parkinson’s disease mechanisms have shown dysregulated extracellular transport of α-synuclein and growth factors in the extracellular space. In the human brain these consist of perineuronal nets, interstitial matrices, and basement membranes, each composed of a set of collagens, non-collagenous glycoproteins, proteoglycans, and hyaluronan. The manner by which amyloidogenic proteins spread extracellularly, become seeded, oligomerize, and are taken up by cells, depends on intricate interactions with extracellular matrix molecules. We sought to assess the alterations to structure of glycosaminoglycans and proteins that occur in PD brain relative to controls of similar age. We found that PD differs markedly from normal brain in upregulation of extracellular matrix structural components including collagens, proteoglycans and glycosaminoglycan binding molecules. We also observed that levels of hemoglobin chains, possibly related to defects in iron metabolism, were enriched in PD brains. These findings shed important new light on disease processes that occur in association with PD.

Community-acquired pneumonia is a widespread disease with significant morbidity and mortality. Alveolar macrophages are tissue-resident lung cells that play a crucial role in innate immunity against bacteria that cause pneumonia. We hypothesized that alveolar macrophages display adaptive characteristics after resolution of bacterial pneumonia. We studied mice 1 to 6 months after self-limiting lung infections with Streptococcus pneumoniae, the most common cause of bacterial pneumonia. Alveolar macrophages, but not other myeloid cells, recovered from the lung showed long-term modifications of their surface marker phenotype. The remodeling of alveolar macrophages was (a) long-lasting (still observed 6 months after infection), (b) regionally localized (observed only in the affected lobe after lobar pneumonia), and (c) associated with macrophage-dependent enhanced protection against another pneumococcal serotype. Metabolomic and transcriptomic profiling revealed that alveolar macrophages of mice that recovered from pneumonia had new baseline activities and altered responses to infection that better resembled those of adult humans. The enhanced lung protection after mild and self-limiting bacterial respiratory infections includes a profound remodeling of the alveolar macrophage pool that is long-lasting; compartmentalized; and manifest across surface receptors, metabolites, and both resting and stimulated transcriptomes.

The Alternative Lengthening of Telomeres (ALT) pathway stimulates telomere elongation and prevents cellular senescence in approximately 60% of osteosarcoma. While the precise mechanism underlying activation of the ALT pathway is unclear, mutations in the chromatin remodeling protein ATRX, histone chaperone DAXX, and the histone variant H3.3 correlate with ALT status. ATRX and DAXX facilitate deposition of the histone variant H3.3 within heterochromatic regions suggesting that loss of ATRX, DAXX, and/or H3.3 lead to defects in the stability of telomeric heterochromatin. Genetic mutations in ATRX, DAXX, and H3.3 have been detected in ALT positive cancers, however, a subset of ALT samples show loss of ATRX or DAXX protein expression or localization without evidence of genetic alterations suggesting additional uncharacterized defects in ATRX/DAXX/H3.3 function. Here, using Next Generation Sequencing we identified a novel gene fusion event between DAXX and the kinesin motor protein, KIFC3, leading to the translation of a chimeric DAXX-KIFC3 fusion protein. Moreover, we demonstrate that the fusion of KIFC3 to DAXX causes defects in DAXX function likely promoting ALT activity. These data highlight a potentially unrecognized mechanism of DAXX inactivation in ALT positive osteosarcoma and provide rationale for thorough and comprehensive analyses of ATRX/DAXX/H3.3 proteins in ALT positive cancers.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Most cellular processes are regulated by transcriptional modulation of specific gene programs. Such modulation is achieved through the combined actions of a wide range of transcription factors (TFs) and cofactors mediating transcriptional activation or repression via changes in chromatin structure. Chromatin immunoprecipitation (ChIP) is a useful molecular biology approach for mapping histone modifications and profiling transcription factors/cofactors binding to DNA, thus providing a snapshot of the dynamic nuclear changes occurring during different biological processes. To study transcriptional regulation in adipose tissue, samples derived from in vitro cell cultures of immortalized or primary cell lines are often favored in ChIP assays because of the abundance of starting material and reduced biological variability. However, these models represent a limited snapshot of the actual chromatin state in living organisms. Thus, there is a critical need for optimized protocols to perform ChIP on adipose tissue samples derived from animal models. Here we describe a protocol for efficient ChIP-seq of both histone modifications and non-histone proteins in brown adipose tissue (BAT) isolated from a mouse. The protocol is optimized for investigating genome-wide localization of proteins of interest and epigenetic markers in the BAT, which is a morphologically and physiologically distinct tissue amongst fat depots.

Huntington’s and Parkinson’s Diseases (HD and PD) are neurodegenerative disorders that share some pathological features but are disparate in others. For example, while both diseases are marked by aberrant protein aggregation in the brain, the specific proteins that aggregate and types of neurons affected differ. A better understanding of the molecular similarities and differences between these two diseases may lead to a more complete mechanistic picture of both the individual diseases and the neurodegenerative process in general. We sought to characterize the common transcriptional signature of HD and PD as well as genes uniquely implicated in each of these diseases using mRNA-Seq data from post mortem human brains in comparison to neuropathologically normal controls. The enriched biological pathways implicated by HD differentially expressed genes show remarkable consistency with those for PD differentially expressed genes and implicate the common biological processes of neuroinflammation, apoptosis, transcriptional dysregulation, and neuron-associated functions. Comparison of the differentially expressed (DE) genes highlights a set of consistently altered genes that span both diseases. In particular, processes involving nuclear factor kappa-light-chain-enhancer of activated B cells (NFkB) and transcription factor cAMP response element-binding protein (CREB) are the most prominent among the genes common to HD and PD. When the combined HD and PD data are compared to controls, relatively few additional biological processes emerge as significantly enriched suggesting that most pathways are independently seen within each disorder. Despite showing comparable numbers of DE genes, DE genes unique to HD are enriched in far more coherent biological processes than the DE genes unique to PD, suggesting that PD may represent a more heterogeneous disorder.

Background Next generation sequencing provides a count of RNA molecules in the form of short reads, yielding discrete, often highly non-normally distributed gene expression measurements. Although Negative Binomial (NB) regression has been generally accepted in the analysis of RNA sequencing (RNA-Seq) data, its appropriateness has not been exhaustively evaluated. We explore logistic regression as an alternative method for RNA-Seq studies designed to compare cases and controls, where disease status is modeled as a function of RNA-Seq reads using simulated and Huntington disease data. We evaluate the effect of adjusting for covariates that have an unknown relationship with gene expression. Finally, we incorporate the data adaptive method in order to compare false positive rates. Results When the sample size is small or the expression levels of a gene are highly dispersed, the NB regression shows inflated Type-I error rates but the Classical logistic and Bayes logistic (BL) regressions are conservative. Firth’s logistic (FL) regression performs well or is slightly conservative. Large sample size and low dispersion generally make Type-I error rates of all methods close to nominal alpha levels of 0.05 and 0.01. However, Type-I error rates are controlled after applying the data adaptive method. The NB, BL, and FL regressions gain increased power with large sample size, large log2 fold-change, and low dispersion. The FL regression has comparable power to NB regression. Conclusions We conclude that implementing the data adaptive method appropriately controls Type-I error rates in RNA-Seq analysis. Firth’s logistic regression provides a concise statistical inference process and reduces spurious associations from inaccurately estimated dispersion parameters in the negative binomial framework. Electronic supplementary material The online version of this article (doi:10.1186/s12859-017-1498-y) contains supplementary material, which is available to authorized users.

Diet-induced obesity (DIO) predisposes individuals to insulin resistance, and adipose tissue has a major role in the disease. Insulin resistance can be induced in cultured adipocytes by a variety of treatments, but what aspects of the in vivo responses are captured by these models remains unknown. We use global RNA sequencing to investigate changes induced by TNF-α, hypoxia, dexamethasone, high insulin, and a combination of TNF-α and hypoxia, comparing the results to the changes in white adipose tissue from DIO mice. We found that different in vitro models capture distinct features of DIO adipose insulin resistance, and a combined treatment of TNF-α and hypoxia is most able to mimic the in vivo changes. Using genome-wide DNase I hypersensitivity followed by sequencing, we further examined the transcriptional regulation of TNF-α-induced insulin resistance, and we found that C/EPBβ is a potential key regulator of adipose insulin resistance., National Institutes of Health (U.S.) (Grant DK-068348), National Institutes of Health (U.S.) (Grant R24 DK-090963), National Institutes of Health (U.S.) (Grant R01GM-089903)

Genome-wide association studies (GWAS) have identified the GAK/DGKQ/IDUA region on 4p16.3 among the top three risk loci for Parkinson’s disease (PD), but the specific gene and risk mechanism are unclear. Here, we report transcripts containing the 3’ clathrin-binding domain of GAK identified by RNA deep-sequencing in post-mortem human brain tissue as having increased expression in PD. Furthermore, carriers of 4p16.3 PD GWAS risk SNPs show decreased expression of one of these transcripts, GAK25 (Gencode Transcript 009), which correlates with the expression of genes functioning in the synaptic vesicle membrane. Together, these findings provide strong evidence for GAK clathrin-binding- and J-domain transcripts’ influence on PD pathogenicity, and for a role for GAK in regulating synaptic function in PD.

Polycomb repressive complexes (PRCs) play key roles in developmental epigenetic regulation. Yet the mechanisms that target PRCs to specific loci in mammalian cells remain incompletely understood. In this study we show that Bmi1, a core component of Polycomb Repressive Complex 1 (PRC1), binds directly to the Runx1/CBFβ transcription factor complex. Genome-wide studies in megakaryocytic cells demonstrate significant chromatin occupancy overlap between the PRC1 core component Ring1b and Runx1/CBFβ and functional regulation of a considerable fraction of commonly bound genes. Bmi1/Ring1b and Runx1/CBFβ deficiencies generate partial phenocopies of one another in vivo. We also show that Ring1b occupies key Runx1 binding sites in primary murine thymocytes and that this occurs via PRC2-independent mechanisms. Genetic depletion of Runx1 results in reduced Ring1b binding at these sites in vivo. These findings provide evidence for site-specific PRC1 chromatin recruitment by core binding transcription factors in mammalian cells.

Huntington’s Disease (HD) is a devastating neurodegenerative disorder that is caused by an expanded CAG trinucleotide repeat in the Huntingtin (HTT) gene. Transcriptional dysregulation in the human HD brain has been documented but is incompletely understood. Here we present a genome-wide analysis of mRNA expression in human prefrontal cortex from 20 HD and 49 neuropathologically normal controls using next generation high-throughput sequencing. Surprisingly, 19% (5,480) of the 28,087 confidently detected genes are differentially expressed (FDR

Background Parkinson disease (PD) is a neurodegenerative disease characterized by the accumulation of alpha-synuclein (SNCA) and other proteins in aggregates termed “Lewy Bodies” within neurons. PD has both genetic and environmental risk factors, and while processes leading to aberrant protein aggregation are unknown, past work points to abnormal levels of SNCA and other proteins. Although several genome-wide studies have been performed for PD, these have focused on DNA sequence variants by genome-wide association studies (GWAS) and on RNA levels (microarray transcriptomics), while genome-wide proteomics analysis has been lacking. Methods This study employed two state-of-the-art technologies, three-stage Mass Spectrometry Tandem Mass Tag Proteomics (12 PD, 12 controls) and RNA-sequencing transcriptomics (29 PD, 44 controls), evaluated in the context of PD GWAS implicated loci and microarray transcriptomics (19 PD, 24 controls). The technologies applied for this study were performed in a set of overlapping prefrontal cortex (Brodmann area 9) samples obtained from PD patients and sex and age similar neurologically healthy controls. Results After appropriate filters, proteomics robustly identified 3558 unique proteins, with 283 of these (7.9 %) significantly different between PD and controls (q-value

Background MicroRNAs (miRNAs) are small non-coding RNAs that recognize sites of complementarity of target messenger RNAs, resulting in transcriptional regulation and translational repression of target genes. In Huntington’s disease (HD), a neurodegenerative disease caused by a trinucleotide repeat expansion, miRNA dyregulation has been reported, which may impact gene expression and modify the progression and severity of HD. Methods We performed next-generation miRNA sequence analysis in prefrontal cortex (Brodmann Area 9) from 26 HD, 2 HD gene positive, and 36 control brains. Neuropathological information was available for all HD brains, including age at disease onset, CAG-repeat size, Vonsattel grade, and Hadzi-Vonsattel striatal and cortical scores, a continuous measure of the extent of neurodegeneration. Linear models were performed to examine the relationship of miRNA expression to these clinical features, and messenger RNA targets of associated miRNAs were tested for gene ontology term enrichment. Results We identified 75 miRNAs differentially expressed in HD brain (FDR q-value

Despite 20 years since its discovery, the gene responsible for Huntington’s Disease, HTT, has still not had its function or transcriptional profile completely characterized. In response to a recent report by Ruzo et al. of several novel splice forms of HTT in human embryonic stem cell lines, we have analyzed a set of mRNA sequencing datasets from post mortem human brain from Huntington’s disease, Parkinson’s disease, and neurologically normal control subjects to evaluate support for previously observed and to identify novel splice patterns. A custom analysis pipeline produced supporting evidence for some of the results reported by two previous studies of alternative isoforms as well as identifying previously unreported splice patterns. All of the alternative splice patterns were of relatively low abundance compared to the canonical splice form.

Cellular signal transduction generally involves cascades of post-translational protein modifications that rapidly catalyze changes in protein-DNA interactions and gene expression. High-throughput measurements are improving our ability to study each of these stages individually, but do not capture the connections between them. Here we present an approach for building a network of physical links among these data that can be used to prioritize targets for pharmacological intervention. Our method recovers the critical missing links between proteomic and transcriptional data by relating changes in chromatin accessibility to changes in expression and then uses these links to connect proteomic and transcriptome data. We applied our approach to integrate epigenomic, phosphoproteomic and transcriptome changes induced by the variant III mutation of the epidermal growth factor receptor (EGFRvIII) in a cell line model of glioblastoma multiforme (GBM). To test the relevance of the network, we used small molecules to target highly connected nodes implicated by the network model that were not detected by the experimental data in isolation and we found that a large fraction of these agents alter cell viability. Among these are two compounds, ICG-001, targeting CREB binding protein (CREBBP), and PKF118–310, targeting β-catenin (CTNNB1), which have not been tested previously for effectiveness against GBM. At the level of transcriptional regulation, we used chromatin immunoprecipitation sequencing (ChIP-Seq) to experimentally determine the genome-wide binding locations of p300, a transcriptional co-regulator highly connected in the network. Analysis of p300 target genes suggested its role in tumorigenesis. We propose that this general method, in which experimental measurements are used as constraints for building regulatory networks from the interactome while taking into account noise and missing data, should be applicable to a wide range of high-throughput datasets., National Science Foundation (U.S.) (DB1-0821391), National Institutes of Health (U.S.) (Grant U54-CA112967), National Institutes of Health (U.S.) (Grant R01-GM089903), National Institutes of Health (U.S.) (P30-ES002109)

The earliest stages of Huntington disease are marked by changes in gene expression that are caused in an indirect and poorly understood manner by polyglutamine expansions in the huntingtin (HTT) protein. To explore the hypothesis that DNA methylation may be altered in cells expressing mutated HTT, we use reduced representation bisulfite sequencing (RRBS) to map sites of DNA methylation in cells carrying either wild-type or mutant HTT. We find that a large fraction of the genes that change in expression in the presence of mutant huntingtin demonstrate significant changes in DNA methylation. Regions with low CpG content, which have previously been shown to undergo methylation changes in response to neuronal activity, are disproportionately affected. On the basis of the sequence of regions that change in methylation, we identify AP-1 and SOX2 as transcriptional regulators associated with DNA methylation changes, and we confirm these hypotheses using genome-wide chromatin immunoprecipitation sequencing (ChIP-Seq). Our findings suggest new mechanisms for the effects of polyglutamine-expanded HTT. These results also raise important questions about the potential effects of changes in DNA methylation on neurogenesis and cognitive decline in patients with Huntington disease.

Huntington’s disease (HD) is an autosomal-dominant neurodegenerative disorder resulting from expansion of CAG repeats in the Huntingtin (HTT) gene. Previous studies have shown mutant HTT can alter expression of genes associated with dysregulated epigenetic modifications. One of the most widely studied chromatin modifications is trimethylated lysine 4 of histone 3 (H3K4me3). Here, we conducted the first comprehensive study of H3K4me3 ChIP-sequencing in neuronal chromatin from the prefrontal cortex of six HD cases and six non-neurologic controls, and its association with gene expression measured by RNA-sequencing. We detected 2,830 differentially enriched H3K4me3 peaks between HD and controls, with 55% of them down-regulated in HD. Although H3K4me3 signals are expected to be associated with mRNA levels, we found an unexpected discordance between altered H3K4me3 peaks and mRNA levels. Gene ontology (GO) term enrichment analysis of the genes with differential H3K4me3 peaks, revealed statistically significantly enriched GO terms only in the genes with down-regulated signals in HD. The most frequently implicated biological process terms are organ morphogenesis and positive regulation of gene expression. More than 9,000 H3K4me3 peaks were located not near any recognized transcription start sites and approximately 36% of these “distal” peaks co-localized to known enhancer sites. Six transcription factors and chromatin remodelers are differentially enriched in HD H3K4me3 distal peaks, including EZH2 and SUZ12, two core subunits of the polycomb repressive complex 2 (PRC2). Moreover, PRC2 repressive state was significantly depleted in HD-enriched peaks, suggesting the epigenetic role of PRC2 inhibition associated with up-regulated H3K4me3 in Huntington’s disease. In summary, our study provides new insights into transcriptional dysregulation of Huntington’s disease by analyzing the differentiation of H3K4me3 enrichment.

The mitogen-activated protein kinase (MAPK) extracellular signal–regulated kinase 2 (ERK2) is ubiquitously expressed in mammalian tissues and is involved in a wide range of biological processes. Although MAPKs have been intensely studied, identification of their substrates remains challenging. We have optimized a chemical genetic system using analog-sensitive ERK2, a form of ERK2 engineered to use an analog of adenosine 5′-triphosphate (ATP), to tag and isolate ERK2 substrates in vitro. This approach identified 80 proteins phosphorylated by ERK2, 13 of which are known ERK2 substrates. The 80 substrates are associated with diverse cellular processes, including regulation of transcription and translation, mRNA processing, and regulation of the activity of the Rho family guanosine triphosphatases. We found that one of the newly identified substrates, ETV3 (a member of the E twenty-six family of transcriptional regulators), was extensively phosphorylated on sites within canonical and noncanonical ERK motifs. Phosphorylation of ETV3 regulated transcription by preventing its binding to DNA at promoters for several thousand genes, including some involved in negative feedback regulation of itself and of upstream signals., Massachusetts Institute of Technology (Eugene Bell Career Development Chair), David H. Koch Institute for Integrative Cancer Research at MIT (Graduate Fellowship), Massachusetts Institute of Technology (Whitaker Health Science Fellowship), National Science Foundation (U.S.) (Graduate Research Fellowship), Pfizer Inc., National Institutes of Health (U.S.) (grant ES002109), National Institutes of Health (U.S.) (grant R01DK42816), National Institutes of Health (U.S.) (grant R01CA118705), National Institutes of Health (U.S.) (grant U54CA112967)

Background Genome-wide computational analysis of alternative splicing (AS) in several flowering plants has revealed that pre-mRNAs from about 30% of genes undergo AS. Chlamydomonas, a simple unicellular green alga, is part of the lineage that includes land plants. However, it diverged from land plants about one billion years ago. Hence, it serves as a good model system to study alternative splicing in early photosynthetic eukaryotes, to obtain insights into the evolution of this process in plants, and to compare splicing in simple unicellular photosynthetic and non-photosynthetic eukaryotes. We performed a global analysis of alternative splicing in Chlamydomonas reinhardtii using its recently completed genome sequence and all available ESTs and cDNAs. Results Our analysis of AS using BLAT and a modified version of the Sircah tool revealed AS of 498 transcriptional units with 611 events, representing about 3% of the total number of genes. As in land plants, intron retention is the most prevalent form of AS. Retained introns and skipped exons tend to be shorter than their counterparts in constitutively spliced genes. The splice site signals in all types of AS events are weaker than those in constitutively spliced genes. Furthermore, in alternatively spliced genes, the prevalent splice form has a stronger splice site signal than the non-prevalent form. Analysis of constitutively spliced introns revealed an over-abundance of motifs with simple repetitive elements in comparison to introns involved in intron retention. In almost all cases, AS results in a truncated ORF, leading to a coding sequence that is around 50% shorter than the prevalent splice form. Using RT-PCR we verified AS of two genes and show that they produce more isoforms than indicated by EST data. All cDNA/EST alignments and splice graphs are provided in a website at http://combi.cs.colostate.edu/as/chlamy. Conclusions The extent of AS in Chlamydomonas that we observed is much smaller than observed in land plants, but is much higher than in simple unicellular heterotrophic eukaryotes. The percentage of different alternative splicing events is similar to flowering plants. Prevalence of constitutive and alternative splicing in Chlamydomonas, together with its simplicity, many available public resources, and well developed genetic and molecular tools for this organism make it an excellent model system to elucidate the mechanisms involved in regulated splicing in photosynthetic eukaryotes.

Alternative splicing (AS) of pre-mRNA is a fundamental molecular process that generates diversity in the transcriptome and proteome of eukaryotic organisms. SR proteins, a family of splicing regulators with one or two RNA recognition motifs (RRMs) at the N-terminus and an arg/ser-rich domain at the C-terminus, function in both constitutive and alternative splicing. We identified SR proteins in 27 eukaryotic species, which include plants, animals, fungi and "basal" eukaryotes that lie outside of these lineages. Using RNA recognition motifs (RRMs) as a phylogenetic marker, we classified 272 SR genes into robust sub-families. The SR gene family can be split into five major groupings, which can be further separated into 11 distinct sub-families. Most flowering plants have double or nearly double the number of SR genes found in vertebrates. The majority of plant SR genes are under purifying selection. Moreover, in all paralogous SR genes in Arabidopsis, rice, soybean and maize, one of the two paralogs is preferentially expressed throughout plant development. We also assessed the extent of AS in SR genes based on a splice graph approach (http://combi.cs.colostate.edu/as/gmap_SRgenes). AS of SR genes is a widespread phenomenon throughout multiple lineages, with alternative 3' or 5' splicing events being the most prominent type of event. However, plant-enriched sub-families have 57%-88% of their SR genes experiencing some type of AS compared to the 40%-54% seen in other sub-families. The SR gene family is pervasive throughout multiple eukaryotic lineages, conserved in sequence and domain organization, but differs in gene number across lineages with an abundance of SR genes in flowering plants. The higher number of alternatively spliced SR genes in plants emphasizes the importance of AS in generating splice variants in these organisms.

Background MicroRNAs (miRNAs) are short, non-coding RNAs that regulate gene expression mainly through translational repression of target mRNA molecules. More than 2700 human miRNAs have been identified and some are known to be associated with disease phenotypes and to display tissue-specific patterns of expression. Methods We used high-throughput small RNA sequencing to discover novel miRNAs in 93 human post-mortem prefrontal cortex samples from individuals with Huntington’s disease (n = 28) or Parkinson’s disease (n = 29) and controls without neurological impairment (n = 36). A custom miRNA identification analysis pipeline was built, which utilizes miRDeep* miRNA identification and result filtering based on false positive rate estimates. Results Ninety-nine novel miRNA candidates with a false positive rate of less than 5 % were identified. Thirty-four of the candidate miRNAs show sequence similarity with known mature miRNA sequences and may be novel members of known miRNA families, while the remaining 65 may constitute previously undiscovered families of miRNAs. Nineteen of the 99 candidate miRNAs were replicated using independent, publicly-available human brain RNA-sequencing samples, and seven were experimentally validated using qPCR. Conclusions We have used small RNA sequencing to identify 99 putative novel miRNAs that are present in human brain samples. Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-3114-3) contains supplementary material, which is available to authorized users.