Search

Tauopathies are neurodegenerative diseases that affect millions of people worldwide including those with Alzheimer’s disease. While many efforts have focused on understanding the role of tau protein in neurodegeneration, there has been little done to systematically analyze and study the structures within tau’s encoding RNA and their connection to disease pathology. Knowledge of RNA structure can provide insights into disease mechanisms and how to affect protein production for therapeutic benefit. Using computational methods based on thermodynamic stability and evolutionary conservation, we identified structures throughout the tau pre-mRNA, especially at exon-intron junctions and within the 5′ and 3′ untranslated regions (UTRs). In particular, structures were identified at twenty exon-intron junctions. The 5′ UTR contains one structured region, which lies within a known internal ribosome entry site. The 3′ UTR contains eight structured regions, including one that contains a polyadenylation signal. A series of functional experiments were carried out to assess the effects of mutations associated with mis-regulation of alternative splicing of exon 10 and to identify regions of the 3′ UTR that contain cis-regulatory elements. These studies defined novel structural regions within the mRNA that affect stability and pre-mRNA splicing and may lead to new therapeutic targets for treating tau-associated diseases. [ABSTRACT FROM AUTHOR]

The secreted protein, insulin-like growth factor 2 (IGF2), plays a central role in fetal and prenatal growth and development, and is regulated at the genetic level by parental imprinting, being expressed predominantly from the paternally derived chromosome in mice and humans. Here, IGF2/Igf2 and its locus has been examined in 19 mammals from 13 orders spanning ~166 million years of evolutionary development. By using human or mouse DNA segments as queries in genome analyses, and by assessing gene expression using RNA-sequencing libraries, more complexity was identified within IGF2/Igf2 than was annotated previously. Multiple potential 5’ non-coding exons were mapped in most mammals and are presumably linked to distinct IGF2/Igf2 promoters, as shown for several species by interrogating RNA-sequencing libraries. DNA similarity was highest in IGF2/Igf2 coding exons; yet, even though the mature IGF2 protein was conserved, versions of 67 or 70 residues are produced secondary to species-specific maintenance of alternative RNA splicing at a variable intron-exon junction. Adjacent H19 was more divergent than IGF2/Igf2, as expected in a gene for a noncoding RNA, and was identified in only 10/19 species. These results show that common features, including those defining IGF2/Igf2 coding and several non-coding exons, were likely present at the onset of the mammalian radiation, but that others, such as a putative imprinting control region 5’ to H19 and potential enhancer elements 3’ to H19, diversified with speciation. This study also demonstrates that careful analysis of genomic and gene expression repositories can provide new insights into gene structure and regulation. [ABSTRACT FROM AUTHOR]

Gestational alcohol exposure causes fetal alcohol spectrum disorder (FASD) and is a prominent cause of neurodevelopmental disability. Whole transcriptome sequencing (RNA-Seq) offer insights into mechanisms underlying FASD, but gene-level analysis provides limited information regarding complex transcriptional processes such as alternative splicing and non-coding RNAs. Moreover, traditional analytical approaches that use multiple hypothesis testing with a false discovery rate adjustment prioritize genes based on an adjusted p-value, which is not always biologically relevant. We address these limitations with a novel approach and implemented an unsupervised machine learning model, which we applied to an exon-level analysis to reduce data complexity to the most likely functionally relevant exons, without loss of novel information. This was performed on an RNA-Seq paired-end dataset derived from alcohol-exposed neural fold-stage chick crania, wherein alcohol causes facial deficits recapitulating those of FASD. A principal component analysis along with k-means clustering was utilized to extract exons that deviated from baseline expression. This identified 6857 differentially expressed exons representing 1251 geneIDs; 391 of these genes were identified in a prior gene-level analysis of this dataset. It also identified exons encoding 23 microRNAs (miRNAs) having significantly differential expression profiles in response to alcohol. We developed an RDAVID pipeline to identify KEGG pathways represented by these exons, and separately identified predicted KEGG pathways targeted by these miRNAs. Several of these (ribosome biogenesis, oxidative phosphorylation) were identified in our prior gene-level analysis. Other pathways are crucial to facial morphogenesis and represent both novel (focal adhesion, FoxO signaling, insulin signaling) and known (Wnt signaling) alcohol targets. Importantly, there was substantial overlap between the exomes themselves and the predicted miRNA targets, suggesting these miRNAs contribute to the gene-level expression changes. Our novel application of unsupervised machine learning in conjunction with statistical analyses facilitated the discovery of signaling pathways and miRNAs that inform mechanisms underlying FASD. [ABSTRACT FROM AUTHOR]

The Cancer Genome Atlas (TCGA) provides a genetic characterization of more than ten thousand tumors, enabling the discovery of novel driver mutations, molecular subtypes, and enticing drug targets across many histologies. Here we investigated why some mutations are common in particular cancer types but absent in others. As an example, we observed that the gene CCDC168 has no mutations in the stomach adenocarcinoma (STAD) cohort despite its common presence in other tumor types. Surprisingly, we found that the lack of called mutations was due to a systematic insufficiency in the number of sequencing reads in the STAD and other cohorts, as opposed to differential driver biology. Using strict filtering criteria, we found similar behavior in four other genes across TCGA cohorts, with each gene exhibiting systematic sequencing depth issues affecting the ability to call mutations. We identified the culprit as the choice of exome capture kit, as kit choice was highly associated with the set of genes that have insufficient reads to call a mutation. Overall, we found that thousands of samples across all cohorts are subject to some capture kit problems. For example, for the 6353 samples using the Broad Institute’s Custom capture kit there are undercalling biases for at least 4833 genes. False negative mutation calls at these genes may obscure biological similarities between tumor types and other important cancer driver effects in TCGA datasets. [ABSTRACT FROM AUTHOR]

Recursive splicing (RS) is an evolutionarily conserved process of removing long introns via multiple steps of splicing. It was first discovered in Drosophila and recently proven to occur also in humans. The detailed mechanism of recursive splicing is not well understood, in particular, whether it is kinetically coupled with transcription. To investigate the dynamic process that underlies recursive splicing, we systematically characterized 342 RS sites in three human cell types using published time-series data that monitored synchronized Pol II elongation and nascent RNA production with 4-thiouridine labeling. We found that half of the RS events occurred post-transcriptionally with long delays. For at least 18–47% RS introns, we detected RS junction reads only after detecting canonical splicing junction reads, supporting the notion that these introns were removed by both recursive splicing and canonical splicing. Furthermore, the choice of which splicing mechanism was used showed cell type specificity. Our results suggest that recursive splicing supplements, rather than replaces, canonical splicing for removing long introns. [ABSTRACT FROM AUTHOR]

Single-cell RNA sequencing (scRNA-seq) technology allows researchers to profile the transcriptomes of thousands of cells simultaneously. Protocols that incorporate both designed and random barcodes have greatly increased the throughput of scRNA-seq, but give rise to a more complex data structure. There is a need for new tools that can handle the various barcoding strategies used by different protocols and exploit this information for quality assessment at the sample-level and provide effective visualization of these results in preparation for higher-level analyses. To this end, we developed scPipe, an R/Bioconductor package that integrates barcode demultiplexing, read alignment, UMI-aware gene-level quantification and quality control of raw sequencing data generated by multiple protocols that include CEL-seq, MARS-seq, Chromium 10X, Drop-seq and Smart-seq. scPipe produces a count matrix that is essential for downstream analysis along with an HTML report that summarises data quality. These results can be used as input for downstream analyses including normalization, visualization and statistical testing. scPipe performs this processing in a few simple R commands, promoting reproducible analysis of single-cell data that is compatible with the emerging suite of open-source scRNA-seq analysis tools available in R/Bioconductor and beyond. The scPipe R package is available for download from . [ABSTRACT FROM AUTHOR]

The exocyst is a molecular tether that retains secretory vesicles at the plasma membrane prior to SNARE-mediated docking and fusion. However, individual exocyst complex components (EXOCs) may also function independently of exocyst assembly. Alternative splice variants of EXOC mRNA and paralogs of EXOC genes have been described and several have been attributed functions that may be independent of the exocyst complex. Here we describe a novel splice variant of murine Exoc3l2, which we term Exoc3l2a. We discuss possible functional implications of the resulting domain excision from this isoform of EXOC3L2 based on structural similarities with its paralog M-Sec (EXOC3L3), which is implicated in tunneling nanotube formation. The identification of this Exoc3l2 splice variant expands the potential for subunit diversity within the exocyst and for alternative functionality of this component independently of the exocyst. [ABSTRACT FROM AUTHOR]

The microtubule associated protein tau plays a critical role in the pathogenesis of neurodegenerative disease. Recent studies suggest that tau also plays a role in disorders of neuronal connectivity, including epilepsy and post-traumatic stress disorder. Animal studies have shown that the MAPT gene, which codes for the tau protein, undergoes complex pre-mRNA alternative splicing to produce multiple isoforms during brain development. Human data, particularly on temporal and regional variation in tau splicing during development are however lacking. In this study, we present the first detailed examination of the temporal and regional sequence of MAPT alternative splicing in the developing human brain. We used a novel computational analysis of large transcriptomic datasets (total n = 502 patients), quantitative polymerase chain reaction (qPCR) and western blotting to examine tau expression and splicing in post-mortem human fetal, pediatric and adult brains. We found that MAPT exons 2 and 10 undergo abrupt shifts in expression during the perinatal period that are unique in the canonical human microtubule-associated protein family, while exon 3 showed small but significant temporal variation. Tau isoform expression may be a marker of neuronal maturation, temporally correlated with the onset of axonal growth. Immature brain regions such as the ganglionic eminence and rhombic lip had very low tau expression, but within more mature regions, there was little variation in tau expression or splicing. We thus demonstrate an abrupt, evolutionarily conserved shift in tau isoform expression during the human perinatal period that may be due to tau expression in maturing neurons. Alternative splicing of the MAPT pre-mRNA may play a vital role in normal brain development across multiple species and provides a basis for future investigations into the developmental and pathological functions of the tau protein. [ABSTRACT FROM AUTHOR]

Proteogenomics methods have identified many non-annotated protein-coding genes in the human genome. Many of the newly discovered protein-coding genes encode peptides and small proteins, referred to collectively as microproteins. Microproteins are produced through ribosome translation of small open reading frames (smORFs). The discovery of many smORFs reveals a blind spot in traditional gene-finding algorithms for these genes. Biological studies have found roles for microproteins in cell biology and physiology, and the potential that there exists additional bioactive microproteins drives the interest in detection and discovery of these molecules. A key step in any proteogenomics workflow is the assembly of RNA-Seq data into likely mRNA transcripts that are then used to create a searchable protein database. Here we demonstrate that specific features of the assembled transcriptome impact microprotein detection by shotgun proteomics. By tailoring transcript assembly for downstream mass spectrometry searching, we show that we can detect more than double the number of high-quality microprotein candidates and introduce a novel open-source mRNA assembler for proteogenomics (MAPS) that incorporates all of these features. By integrating our specialized assembler, MAPS, and a popular generalized assembler into our proteogenomics pipeline, we detect 45 novel human microproteins from a high quality proteogenomics dataset of a human cell line. We then characterize the features of the novel microproteins, identifying two classes of microproteins. Our work highlights the importance of specialized transcriptome assembly upstream of proteomics validation when searching for short and potentially rare and poorly conserved proteins. [ABSTRACT FROM AUTHOR]

Eukaryotic translation initiation factor 2-alpha kinase (EIF2AK) proteins inhibit protein synthesis at translation initiation level, in response to various stress conditions, including oxidative stress, heme deficiency, osmotic shock, and heat shock. Origin and functional diversification of EIF2AK sequences remain ambiguous. Here we determine the origin and molecular evolution of EIF2AK proteins in lower eukaryotes and studied the molecular basis of divergence among sub-family sequences. Present work emphasized primitive origin of EIF2AK4 sub-family gene in lower eukaryotes of protozoan lineage. Phylogenetic analysis supported common origin and sub-family based classification of EIF2AKs. Functional divergence studies across sub-families revealed several putative amino acid sites, which assist in altered protein interactions of kinase domains. The data can facilitate designing site-directed experimental studies aiming at elucidating diverse functional aspects of kinase domains regarding down-regulation of protein synthesis. [ABSTRACT FROM AUTHOR]

Substitutions that disrupt pre-mRNA splicing are a common cause of genetic disease. On average, 13.4% of all hereditary disease alleles are classified as splicing mutations mapping to the canonical 5′ and 3′ splice sites. However, splicing mutations present in exons and deeper intronic positions are vastly underreported. A recent re-analysis of coding mutations in exon 10 of the Lynch Syndrome gene, MLH1, revealed an extremely high rate (77%) of mutations that lead to defective splicing. This finding is confirmed by extending the sampling to five other exons in the MLH1 gene. Further analysis suggests a more general phenomenon of defective splicing driving Lynch Syndrome. Of the 36 mutations tested, 11 disrupted splicing. Furthermore, analyzing past reports suggest that MLH1 mutations in canonical splice sites also occupy a much higher fraction (36%) of total mutations than expected. When performing a comprehensive analysis of splicing mutations in human disease genes, we found that three main causal genes of Lynch Syndrome, MLH1, MSH2, and PMS2, belonged to a class of 86 disease genes which are enriched for splicing mutations. Other cancer genes were also enriched in the 86 susceptible genes. The enrichment of splicing mutations in hereditary cancers strongly argues for additional priority in interpreting clinical sequencing data in relation to cancer and splicing. [ABSTRACT FROM AUTHOR]

Nephronophthisis-related ciliopathies (NPHP-RC) are a group of disorders that present with end-stage renal failure in childhood/adolescence, kidney cysts, retinal degeneration, and cerebellar hypoplasia. One disorder that shares clinical features with NPHP-RC is Bardet-Biedl Syndrome (BBS). Serologically defined colon cancer antigen 8 (SDCCAG8; also known as NPHP10 and BBS16) is an NPHP gene that is also associated with BBS. To better understand the patho-mechanisms of NPHP and BBS caused by loss of SDCCAG8 function, we characterized an SDCCAG8 mouse model (Sdccag8Tn(sb-Tyr)2161B . CA1C2Ove) generated by Sleeping Beauty Transposon (SBT)-mediated insertion mutagenesis. Consistent with the previously reported, independent SDCCAG8 mouse models, our mutant mice display pre-axial polydactyly in their hind limbs. In addition, we report patterning defects in the secondary palate, brain abnormalities, as well as neonatal lethality associated with developmental defects in the lung in our mouse model. The neonatal lethality phenotype is genetic background dependent and rescued by introducing 129S6/SvEvTac background. Genetic modifier(s) responsible for this effect were mapped to a region between SNPs rs3714172 and rs3141832 on chromosome 11. While determining the precise genetic lesion in our mouse model, we found that SBT insertion resulted in a deletion of multiple exons from both Sdccag8 and its neighboring gene Akt3. We ascribe the patterning defects in the limb and the secondary palate as well as lung abnormalities to loss of SDCCAG8, while the developmental defects in the brain are likely due to the loss of AKT3. This mouse model may be useful to study features not observed in other SDCCAG8 models but cautions are needed in interpreting data. [ABSTRACT FROM AUTHOR]

Insulin-like growth factor 1 (IGF1), a small, secreted peptide growth factor, is involved in a variety of physiological and patho-physiological processes, including somatic growth, tissue repair, and metabolism of carbohydrates, proteins, and lipids. IGF1 gene expression appears to be controlled by several different signaling cascades in the few species in which it has been evaluated, with growth hormone playing a major role by activating a pathway involving the Stat5b transcription factor. Here, genes encoding IGF1 have been evaluated in 25 different mammalian species representing 15 different orders and ranging over ~180 million years of evolutionary diversification. Parts of the IGF1 gene have been fairly well conserved. Like rat Igf1 and human IGF1, 21 of 23 other genes are composed of 6 exons and 5 introns, and all 23 also contain recognizable tandem promoters, each with a unique leader exon. Exon and intron lengths are similar in most species, and DNA sequence conservation is moderately high in orthologous exons and proximal promoter regions. In contrast, putative growth hormone-activated Stat5b-binding enhancers found in analogous locations in rodent Igf1 and in human IGF1 loci, have undergone substantial variation in other mammals, and a processed retro-transposed IGF1 pseudogene is found in the sloth locus, but not in other mammalian genomes. Taken together, the fairly high level of organizational and nucleotide sequence similarity in the IGF1 gene among these 25 species supports the contention that some common regulatory pathways had existed prior to the beginning of mammalian speciation. [ABSTRACT FROM AUTHOR]

ciRS-7 is an intensely studied, highly expressed and conserved circRNA. Essentially nothing is known about its biogenesis, including the location of its promoter. A prevailing assumption has been that ciRS-7 is an exceptional circRNA because it is transcribed from a locus lacking any mature linear RNA transcripts of the same sense. To study the biogenesis of ciRS-7, we developed an algorithm to define its promoter and predicted that the human ciRS-7 promoter coincides with that of the long non-coding RNA, LINC00632. We validated this prediction using multiple orthogonal experimental assays. We also used computational approaches and experimental validation to establish that ciRS-7 exonic sequence is embedded in linear transcripts that are flanked by cryptic exons in both human and mouse. Together, this experimental and computational evidence generate a new model for regulation of this locus: (a) ciRS-7 is like other circRNAs, as it is spliced into linear transcripts; (b) expression of ciRS-7 is primarily determined by the chromatin state of LINC00632 promoters; (c) transcription and splicing factors sufficient for ciRS-7 biogenesis are expressed in cells that lack detectable ciRS-7 expression. These findings have significant implications for the study of the regulation and function of ciRS-7, and the analytic framework we developed to jointly analyze RNA-seq and ChIP-seq data reveal the potential for genome-wide discovery of important biological regulation missed in current reference annotations. [ABSTRACT FROM AUTHOR]

Troxerutin, a semi-synthetic derivative of the natural bioflavanoid rutin, has been reported to possess many beneficial effects in human bodies, such as vasoprotection, immune support, anti-inflammation and anti-aging. However, the effects of troxerutin on genome-wide transcription in blood cells are still unknown. In order to find out effects of troxerutin on gene transcription, a high-throughput RNA sequencing was employed to analysis differential gene expression in blood cells consisting of leucocytes, erythrocytes and platelets isolated from the mice received subcutaneous injection of troxerutin. Transcriptome analysis demonstrated that the expression of only fifteen genes was significantly changed by the treatment with troxerutin, among which 5 genes were up-regulated and 10 genes were down-regulated. Bioinformatic analysis of the fifteen differentially expressed genes was made by utilizing the Gene Ontology (GO), and the differential expression induced by troxerutin was further evaluated by real-time quantitative PCR (Q-PCR). [ABSTRACT FROM AUTHOR]

Alternative splicing is an ubiquitous phenomenon in most human genes and has important functions. The switch-like exon is the type of exon that has a high level of usage in some tissues, but has a low level of usage in the other tissues. They usually undergo strong tissue-specific regulations. There is still a lack a systematic method to identify switch-like exons from multiple RNA-seq samples. We proposed a novel method called iterative Tertile Absolute Deviation around the mode (iTAD) to profile the distribution of exon relative usages among multiple samples and to identify switch-like exons and other types of exons using a robust statistic estimator. We validated the method with simulation data, and applied it on RNA-seq data of 16 human body tissues and detected 3,100 switch-like exons. We found that switch-like exons tend to be more associated with Alu elements in their flanking intron regions than other types of exons. [ABSTRACT FROM AUTHOR]

The phylum Placozoa holds a key position for our understanding of the evolution of mitochondrial genomes in Metazoa. Placozoans possess large mitochondrial genomes which harbor several remarkable characteristics such as a fragmented cox1 gene and trans-splicing cox1 introns. A previous study also suggested the existence of cox1 mRNA editing in Trichoplax adhaerens, yet the only formally described species in the phylum Placozoa. We have analyzed RNA-seq data of the undescribed sister species, Placozoa sp. H2 (“Panama” clone), with special focus on the mitochondrial mRNA. While we did not find support for a previously postulated cox1 mRNA editing mechanism, we surprisingly found two independent transcripts representing intermediate cox1 mRNA splicing stages. Both transcripts consist of partial cox1 exon as well as overlapping intron fragments. The data suggest that the cox1 gene harbors a single base pair (cytosine) micro exon. Furthermore, conserved group I intron structures flank this unique micro exon also in other placozoans. We discuss the evolutionary origin of this micro exon in the context of a self-splicing intron gain in the cox1 gene of the last common ancestor of extant placozoans. [ABSTRACT FROM AUTHOR]

The human cytomegalovirus major immediate early proteins IE1 and IE2 are critical drivers of virus replication and are considered pivotal in determining the balance between productive and latent infection. IE1 and IE2 are derived from the same primary transcript by alternative splicing and regulation of their expression likely involves a complex interplay between cellular and viral factors. Here we show that knockdown of the host ubiquitin-dependent segregase VCP/p97, results in loss of IE2 expression, subsequent suppression of early and late gene expression and, ultimately, failure in virus replication. RNAseq analysis showed increased levels of IE1 splicing, with a corresponding decrease in IE2 splicing following VCP knockdown. Global analysis of viral transcription showed the expression of a subset of viral genes is not reduced despite the loss of IE2 expression, including UL112/113. Furthermore, Immunofluorescence studies demonstrated that VCP strongly colocalised with the viral replication compartments in the nucleus. Finally, we show that NMS-873, a small molecule inhibitor of VCP, is a potent HCMV antiviral with potential as a novel host targeting therapeutic for HCMV infection. [ABSTRACT FROM AUTHOR]

Mutation screening of the breast cancer genes BRCA1 and BRCA2 identifies a large fraction of variants of uncertain clinical significance (VUS) whose functional and clinical interpretations pose a challenge for genomic medicine. Likewise, an increasing amount of evidence indicates that genetic variants can have deleterious effects on pre-mRNA splicing. Our goal was to investigate the impact on splicing of a set of reported variants of BRCA2 exons 17 and 18 to assess their role in hereditary breast cancer and to identify critical regulatory elements that may constitute hotspots for spliceogenic variants. A splicing reporter minigene with BRCA2 exons 14 to-20 (MGBR2_ex14-20) was constructed in the pSAD vector. Fifty-two candidate variants were selected with splicing prediction programs, introduced in MGBR2_ex14-20 by site-directed mutagenesis and assayed in triplicate in MCF-7 cells. Wild type MGBR2_ex14-20 produced a stable transcript of the expected size (1,802 nucleotides) and structure (V1-[BRCA2_exons_14–20]–V2). Functional mapping by microdeletions revealed essential sequences for exon recognition on the 3’ end of exon 17 (c.7944-7973) and the 5’ end of exon 18 (c.7979-7988, c.7999-8013). Thirty out of the 52 selected variants induced anomalous splicing in minigene assays with >16 different aberrant transcripts, where exon skipping was the most common event. A wide range of splicing motifs were affected including the canonical splice sites (15 variants), novel alternative sites (3 variants), the polypyrimidine tract (3 variants) and enhancers/silencers (9 variants). According to the guidelines of the American College of Medical Genetics and Genomics (ACMG), 20 variants could be classified as pathogenic (c.7806-2A>G, c.7806-1G>A, c.7806-1G>T, c.7806-1_7806-2dup, c.7976+1G>A, c.7977-3_7978del, c.7977-2A>T, c.7977-1G>T, c.7977-1G>C, c.8009C>A, c.8331+1G>T and c.8331+2T>C) or likely pathogenic (c.7806-9T>G, c.7976G>C, c.7976G>A, c.7977-7C>G, c.7985C>G, c.8023A>G, c.8035G>T and c.8331G>A), accounting for 30.8% of all pathogenic/likely pathogenic variants of exons 17–18 at the BRCA Share database. The remaining 8 variants (c.7975A>G, c.7977-6T>G, c.7988A>T, c.7992T>A, c.8007A>G, c.8009C>T, c.8009C>G, and c.8072C>T) induced partial splicing anomalies with important ratios of the full-length transcript (≥70%), so that they remained classified as VUS. Aberrant splicing is therefore especially prevalent in BRCA2 exons 17 and 18 due to the presence of active ESEs involved in exon recognition. Splicing functional assays with minigenes are a valuable strategy for the initial characterization of the splicing outcomes and the subsequent clinical interpretation of variants of any disease-gene, although these results should be checked, whenever possible, against patient RNA. [ABSTRACT FROM AUTHOR]

Antiviral defenses are inappropriately activated in systemic lupus erythematosus (SLE) and association between SLE and the antiviral helicase gene, IFIH1, is well established. We sought to extend the previously reported association of pathogenic soluble mediators and autoantibodies with mouse Mda5 to its human ortholog, IFIH1. To better understand the role this gene plays in human lupus, we assessed association of IFIH1 variants with soluble mediators and autoantibodies in 357 European-American SLE patients, first-degree relatives, and unrelated, unaffected healthy controls. Association between each of 135 genotyped SNPs in IFIH1 and four lupus-associated plasma mediators, IL-6, TNF-α, IFN-β, and IP-10, were investigated via linear regression. No significant associations were found to SNPs orthologous to those identified in exon 13 of the mouse. However, outside of this region there were significant associations between IL-6 and rs76162067 (p = 0.008), as well as IP-10 and rs79711023 (p = 0.003), located in a region of IFIH1 previously shown to directly influence MDA-5 mediated IP-10 and IL-6 secretion. SLE patients and FDRs carrying the minor allele for rs79711023 demonstrated lower levels of IP-10, while only FDRs carrying the minor allele for rs76162067 demonstrated an increased level of IL-6. This would suggest that the change in IP-10 is genotypically driven, while the change in IL-6 may be reflective of SLE transition status. These data suggest that IFIH1 may contribute to SLE pathogenesis via altered inflammatory mechanisms. [ABSTRACT FROM AUTHOR]

RNA-seq reads containing part of the poly(A) tail of transcripts (denoted as poly(A) reads) provide the most direct evidence for the position of poly(A) sites in the genome. However, due to reduced coverage of poly(A) tails by reads, poly(A) reads are not routinely identified during RNA-seq mapping. Nevertheless, recent studies for several herpesviruses successfully employed mapping of poly(A) reads to identify herpesvirus poly(A) sites using different strategies and customized programs. To more easily allow such analyses without requiring additional programs, we integrated poly(A) read mapping and prediction of poly(A) sites into our RNA-seq mapping program ContextMap 2. The implemented approach essentially generalizes previously used poly(A) read mapping approaches and combines them with the context-based approach of ContextMap 2 to take into account information provided by other reads aligned to the same location. Poly(A) read mapping using ContextMap 2 was evaluated on real-life data from the ENCODE project and compared against a competing approach based on transcriptome assembly (KLEAT). This showed high positive predictive value for our approach, evidenced also by the presence of poly(A) signals, and considerably lower runtime than KLEAT. Although sensitivity is low for both methods, we show that this is in part due to a high extent of spurious results in the gold standard set derived from RNA-PET data. Sensitivity improves for poly(A) sites of known transcripts or determined with a more specific poly(A) sequencing protocol and increases with read coverage on transcript ends. Finally, we illustrate the usefulness of the approach in a high read coverage scenario by a re-analysis of published data for herpes simplex virus 1. Thus, with current trends towards increasing sequencing depth and read length, poly(A) read mapping will prove to be increasingly useful and can now be performed automatically during RNA-seq mapping with ContextMap 2. [ABSTRACT FROM AUTHOR]

Alternative splicing is a vital process for regulating gene expression and promoting proteomic diversity. It plays a key role in tissue-specific expressed genes. This specificity is mainly regulated by splicing factors that bind to specific sequences called splicing regulatory elements (SREs). Here, we report a genome-wide analysis to study alternative splicing on multiple tissues, including brain, heart, liver, and muscle. We propose a pipeline to identify differential exons across tissues and hence tissue-specific SREs. In our pipeline, we utilize the DEXSeq package along with our previously reported algorithms. Utilizing the publicly available RNA-Seq data set from the Human BodyMap project, we identified 28,100 differentially used exons across the four tissues. We identified tissue-specific exonic splicing enhancers that overlap with various previously published experimental and computational databases. A complicated exonic enhancer regulatory network was revealed, where multiple exonic enhancers were found across multiple tissues while some were found only in specific tissues. Putative combinatorial exonic enhancers and silencers were discovered as well, which may be responsible for exon inclusion or exclusion across tissues. Some of the exonic enhancers are found to be co-occurring with multiple exonic silencers and vice versa, which demonstrates a complicated relationship between tissue-specific exonic enhancers and silencers. [ABSTRACT FROM AUTHOR]

‘Microtubule-associated protein tau’ (MAPT), ‘granulin’ (GRN) and ‘chromosome 9 open reading frame72’ (C9ORF72) gene mutations are the major known genetic causes of frontotemporal dementia (FTD). Recent studies suggest that mutations in these genes may also be associated with other forms of dementia. Therefore we investigated whether MAPT, GRN and C9ORF72 gene mutations are major contributors to dementia in a random, unselected Turkish cohort of dementia patients. A combination of whole-exome sequencing, Sanger sequencing and fragment analysis/Southern blot was performed in order to identify pathogenic mutations and novel variants in these genes as well as other FTD-related genes such as the ‘charged multivesicular body protein 2B’ (CHMP2B), the ‘FUS RNA binding protein’ (FUS), the ‘TAR DNA binding protein’ (TARDBP), the ‘sequestosome1’ (SQSTM1), and the ‘valosin containing protein’ (VCP). We determined one pathogenic MAPT mutation (c.1906C>T, p.P636L) and one novel missense variant (c.38A>G, p.D13G). In GRN we identified a probably pathogenic TGAG deletion in the splice donor site of exon 6. Three patients were found to carry the GGGGCC expansions in the non-coding region of the C9ORF72 gene. In summary, a complete screening for mutations in MAPT, GRN and C9ORF72 genes revealed a frequency of 5.4% of pathogenic mutations in a random cohort of 93 Turkish index patients with dementia. [ABSTRACT FROM AUTHOR]

Repetitive elements (REs) comprise 40–60% of the mammalian genome and have been shown to epigenetically influence the expression of genes through the formation of fusion transcript (FTs). We previously showed that an intracisternal A particle forms an FT with the agouti gene in mice, causing obesity/type 2 diabetes. To determine the frequency of FTs genome-wide, we developed a TopHat-Fusion-based analytical pipeline to identify FTs with high specificity. We applied it to an RNA-seq dataset from the nucleus accumbens (NAc) of mice repeatedly exposed to cocaine. Cocaine was previously shown to increase the expression of certain REs in this brain region. Using this pipeline that can be applied to single- or paired-end reads, we identified 438 genes expressing 813 different FTs in the NAc. Although all types of studied repeats were present in FTs, simple sequence repeats were underrepresented. Most importantly, reverse-transcription and quantitative PCR validated the expression of selected FTs in an independent cohort of animals, which also revealed that some FTs are the prominent isoforms expressed in the NAc by some genes. In other RNA-seq datasets, developmental expression as well as tissue specificity of some FTs differed from their corresponding non-fusion counterparts. Finally, in silico analysis predicted changes in the structure of proteins encoded by some FTs, potentially resulting in gain or loss of function. Collectively, these results indicate the robustness of our pipeline in detecting these new isoforms of genes, which we believe provides a valuable tool to aid in better understanding the broad role of REs in mammalian cellular biology. [ABSTRACT FROM AUTHOR]

Mutations in the X-linked Cyclin-Dependent Kinase-Like 5 gene (CDKL5) cause early onset infantile spasms and subsequent severe developmental delay in affected children. Deleterious mutations have been reported to occur throughout the CDKL5 coding region. Several studies point to a complex CDKL5 gene structure in terms of exon usage and transcript expression. Improvements in molecular diagnosis and more extensive research into the neurobiology of CDKL5 and pathophysiology of CDKL5 disorders necessitate an updated analysis of the gene. In this study, we have analysed human and mouse CDKL5 transcript patterns both bioinformatically and experimentally. We have characterised the predominant brain isoform of CDKL5, a 9.7 kb transcript comprised of 18 exons with a large 6.6 kb 3’-untranslated region (UTR), which we name hCDKL5_1. In addition we describe new exonic regions and a range of novel splice and UTR isoforms. This has enabled the description of an updated gene model in both species and a standardised nomenclature system for CDKL5 transcripts. Profiling revealed tissue- and brain development stage-specific differences in expression between transcript isoforms. These findings provide an essential backdrop for the diagnosis of CDKL5-related disorders, for investigations into the basic biology of this gene and its protein products, and for the rational design of gene-based and molecular therapies for these disorders. [ABSTRACT FROM AUTHOR]

Analysis of splice variants from short read RNA-seq data remains a challenging problem. Here we present a novel method for the genome-guided prediction and quantification of splice events from RNA-seq data, which enables the analysis of unannotated and complex splice events. Splice junctions and exons are predicted from reads mapped to a reference genome and are assembled into a genome-wide splice graph. Splice events are identified recursively from the graph and are quantified locally based on reads extending across the start or end of each splice variant. We assess prediction accuracy based on simulated and real RNA-seq data, and illustrate how different read aligners (GSNAP, HISAT2, STAR, TopHat2) affect prediction results. We validate our approach for quantification based on simulated data, and compare local estimates of relative splice variant usage with those from other methods (MISO, Cufflinks) based on simulated and real RNA-seq data. In a proof-of-concept study of splice variants in 16 normal human tissues (Illumina Body Map 2.0) we identify 249 internal exons that belong to known genes but are not related to annotated exons. Using independent RNA samples from 14 matched normal human tissues, we validate 9/9 of these exons by RT-PCR and 216/249 by paired-end RNA-seq (2 x 250 bp). These results indicate that de novo prediction of splice variants remains beneficial even in well-studied systems. An implementation of our method is freely available as an R/Bioconductor package . [ABSTRACT FROM AUTHOR]

Idiopathic pulmonary fibrosis (IPF) is a progressive, debilitating disease for which two medications, pirfenidone and nintedanib, have only recently been approved for treatment. The cytokine TGF-β has been shown to be a central mediator in the disease process. We investigated the role of a novel kinase, MAP3K19, upregulated in IPF tissue, in TGF-β-induced signal transduction and in bleomycin-induced pulmonary fibrosis. MAP3K19 has a very limited tissue expression, restricted primarily to the lungs and trachea. In pulmonary tissue, expression was predominantly localized to alveolar and interstitial macrophages, bronchial epithelial cells and type II pneumocytes of the epithelium. MAP3K19 was also found to be overexpressed in bronchoalveolar lavage macrophages from IPF patients compared to normal patients. Treatment of A549 or THP-1 cells with either MAP3K19 siRNA or a highly potent and specific inhibitor reduced phospho-Smad2 & 3 nuclear translocation following TGF-β stimulation. TGF-β-induced gene transcription was also strongly inhibited by both the MAP3K19 inhibitor and nintedanib, whereas pirfenidone had a much less pronounced effect. In combination, the MAP3K19 inhibitor appeared to act synergistically with either pirfenidone or nintedanib, at the level of target gene transcription or protein production. Finally, in an animal model of IPF, inhibition of MAP3K19 strongly attenuated bleomycin-induced pulmonary fibrosis when administered either prophylactically ortherapeutically. In summary, these results strongly suggest that inhibition of MAP3K19 may have a beneficial therapeutic effect in the treatment of IPF and represents a novel strategy to target this disease. [ABSTRACT FROM AUTHOR]

The phase III Thai RV144 vaccine trial showed an estimated vaccine efficacy (VE) to prevent HIV-1 infection of 31.2%, which has motivated the search for immune correlates of vaccine protection. In a recent report, several single nucleotide polymorphisms (SNPs) in FCGR2C were identified to associate with the level of VE in the RV144 trial. To investigate the functional significance of these SNPs, we utilized a large scale B cell RNA sequencing database of 462 individuals from the 1000 Genomes Project to examine associations between FCGR2C SNPs and gene expression. We found that the FCGR2C SNPs that associated with vaccine efficacy in RV144 also strongly associated with the expression of FCGR2A/C and one of them also associated with the expression of Fc receptor-like A (FCRLA), another Fc-γ receptor (FcγR) gene that was not examined in the previous report. These results suggest that the expression of FcγR genes is influenced by these SNPs either directly or in linkage with other causal variants. More importantly, these results motivate further investigations into the potential for a causal association of expression and alternative splicing of FCGR2C and other FcγR genes with the HIV-1 vaccine protection in the RV144 trial and other similar studies. [ABSTRACT FROM AUTHOR]

The unintended consequences of gene targeting in mouse models have not been thoroughly studied and a more systematic analysis is needed to understand the frequency and characteristics of off-target effects. Using RNA-seq, we evaluated targeted and neighboring gene expression in tissues from 44 homozygous mutants compared with C57BL/6N control mice. Two allele types were evaluated: 15 targeted trap mutations (TRAP); and 29 deletion alleles (DEL), usually a deletion between the translational start and the 3’ UTR. Both targeting strategies insert a bacterial beta-galactosidase reporter (LacZ) and a neomycin resistance selection cassette. Evaluating transcription of genes in +/- 500 kb of flanking DNA around the targeted gene, we found up-regulated genes more frequently around DEL compared with TRAP alleles, however the frequency of alleles with local down-regulated genes flanking DEL and TRAP targets was similar. Down-regulated genes around both DEL and TRAP targets were found at a higher frequency than expected from a genome-wide survey. However, only around DEL targets were up-regulated genes found with a significantly higher frequency compared with genome-wide sampling. Transcriptome analysis confirms targeting in 97% of DEL alleles, but in only 47% of TRAP alleles probably due to non-functional splice variants, and some splicing around the gene trap. Local effects on gene expression are likely due to a number of factors including compensatory regulation, loss or disruption of intragenic regulatory elements, the exogenous promoter in the neo selection cassette, removal of insulating DNA in the DEL mutants, and local silencing due to disruption of normal chromatin organization or presence of exogenous DNA. An understanding of local position effects is important for understanding and interpreting any phenotype attributed to targeted gene mutations, or to spontaneous indels. [ABSTRACT FROM AUTHOR]