Transcriptome content and dynamics at single-nucleotide resolution

Massively parallel short-tag 'RNA sequencing' is giving new insights into the complexity of the eukaryotic transcriptome

The uniqueome: a mappability resource for short-tag sequencing

Summary: Quantification applications of short-tag sequencing data (such as CNVseq and RNAseq) depend on knowing the uniqueness of specific genomic regions at a given threshold of error. Here, we present the ‘uniqueome’, a genomic resource for understanding the uniquely mappable proportion of genomic sequences. Pre-computed data are available for human, mouse, fly and worm genomes in both color-space and nucletotide-space, and we demonstrate the utility of this resource as applied to the quantification of RNAseq data

Profiling Gene Expression Induced by Protease-Activated Receptor 2 (PAR2) Activation in Human Kidney Cells

Protease-Activated Receptor-2 (PAR2) has been implicated through genetic knockout mice with cytokine regulation and arthritis development. Many studies have associated PAR2 with inflammatory conditions (arthritis, airways inflammation, IBD) and key events in tumor progression (angiogenesis, metastasis), but they have relied heavily on the use of single agonists to identify physiological roles for PAR2. However such probes are now known not to be highly selective for PAR2, and thus precisely what PAR2 does and what mechanisms of downstream regulation are truly affected remain obscure. Effects of PAR2 activation on gene expression in Human Embryonic Kidney cells (HEK293), a commonly studied cell line in PAR2 research, were investigated here by comparing 19,000 human genes for intersecting up- or down-regulation by both trypsin (an endogenous protease that activates PAR2) and a PAR2 activating hexapeptide (2f-LIGRLO-NH2). Among 2,500 human genes regulated similarly by both agonists, there were clear associations between PAR2 activation and cellular metabolism (1,000 genes), the cell cycle, the MAPK pathway, HDAC and sirtuin enzymes, inflammatory cytokines, and anti-complement function. PAR-2 activation up-regulated four genes more than 5 fold (DUSP6, WWOX, AREG, SERPINB2) and down-regulated another six genes more than 3 fold (TXNIP, RARG, ITGB4, CTSD, MSC and TM4SF15). Both PAR2 and PAR1 activation resulted in up-regulated expression of several genes (CD44, FOSL1, TNFRSF12A, RAB3A, COPEB, CORO1C, THBS1, SDC4) known to be important in cancer. This is the first widespread profiling of specific activation of PAR2 and provides a valuable platform for better understanding key mechanistic roles of PAR2 in human physiology. Results clearly support the development of both antagonists and agonists of human PAR2 as potential disease modifying therapeutic agents

Integrated Genome Analysis Suggests that Most Conserved Non-Coding Sequences are Regulatory Factor Binding Sites

More than 98% of a typical vertebrate genome does not code for proteins. Although non-coding regions are sprinkled with short (<200 bp) islands of evolutionarily conserved sequences, the function of most of these unannotated conserved islands remains unknown. One possibility is that unannotated conserved islands could encode non-coding RNAs (ncRNAs); alternatively, unannotated conserved islands could serve as promoter-distal regulatory factor binding sites (RFBSs) like enhancers. Here we assess these possibilities by comparing unannotated conserved islands in the human and mouse genomes to transcribed regions and to RFBSs, relying on a detailed case study of one human and one mouse cell type. We define transcribed regions by applying a novel transcript-calling algorithm to RNA-Seq data obtained from total cellular RNA, and we define RFBSs using ChIP-Seq and DNAse-hypersensitivity assays. We find that unannotated conserved islands are four times more likely to coincide with RFBSs than with unannotated ncRNAs. Thousands of conserved RFBSs can be categorized as insulators based on the presence of CTCF or as enhancers based on the presence of p300/CBP and H3K4me1. While many unannotated conserved RFBSs are transcriptionally active to some extent, the transcripts produced tend to be unspliced, non-polyadenylated and expressed at levels 10 to 100-fold lower than annotated coding or ncRNAs. Extending these findings across multiple cell types and tissues, we propose that most conserved non-coding genomic DNA in vertebrate genomes corresponds to promoter-distal regulatory elements

Dynamic transcription programs during ES cell differentiation towards mesoderm in serum versus serum-freeBMP4 culture

Expression profiling of embryonic stem (ES) cell differentiation in the presence of serum has been performed previously. It remains unclear if transcriptional activation is dependent on complext growth factor mixtures in serum or whether this process is intrinsic to ES cells once the stem cell program has been inactivated. The aims of this study were to determine the transcriptional programs associated with the stem cell state and to characterize mesoderm differentiation between serum and serum-free culture

Exome-wide association study of pancreatic cancer risk

We conducted a case-control exome-wide association study to discover germline variants in coding regions that affect risk for pancreatic cancer, combining data from 5 studies. We analyzed exome and genome sequencing data from 437 patients with pancreatic cancer (cases) and 1922 individuals not known to have cancer (controls). In the primary analysis, BRCA2 had the strongest enrichment for rare inactivating variants (17/437 cases vs 3/1922 controls) (P=3.27x10(-6); exome-wide statistical significance threshold P&lt;2.5x10(-6)). Cases had more rare inactivating variants in DNA repair genes than controls, even after excluding 13 genes known to predispose to pancreatic cancer (adjusted odds ratio, 1.35, P=.045). At the suggestive threshold (P&lt;.001), 6 genes were enriched for rare damaging variants (UHMK1, AP1G2, DNTA, CHST6, FGFR3, and EPHA1) and 7 genes had associations with pancreatic cancer risk, based on the sequence-kernel association test. We confirmed variants in BRCA2 as the most common high-penetrant genetic factor associated with pancreatic cancer and we also identified candidate pancreatic cancer genes. Large collaborations and novel approaches are needed to overcome the genetic heterogeneity of pancreatic cancer predisposition

The uniqueome: a mappability resource for short-tag sequencing

Summary: Quantification applications of short-tag sequencing data (such as CNVseq and RNAseq) depend on knowing the uniqueness of specific genomic regions at a given threshold of error. Here, we present the ‘uniqueome’, a genomic resource for understanding the uniquely mappable proportion of genomic sequences. Pre-computed data are available for human, mouse, fly and worm genomes in both color-space and nucletotide-space, and we demonstrate the utility of this resource as applied to the quantification of RNAseq data

X-MATE: a flexible system for mapping short read data

Summary: Accurate and complete mapping of short-read sequencing to a reference genome greatly enhances the discovery of biological results and improves statistical predictions. We recently presented RNA-MATE, a pipeline for the recursive mapping of RNA-Seq datasets. With the rapid increase in genome re-sequencing projects, progression of available mapping software and the evolution of file formats, we now present X-MATE, an updated version of RNA-MATE, capable of mapping both RNA-Seq and DNA datasets and with improved performance, output file formats, configuration files, and flexibility in core mapping software

Ferrets exclusively synthesize Neu5Ac and express naturally humanized influenza A virus receptors

Mammals express the sialic acids ​N-acetylneuraminic acid (​Neu5Ac) and ​N-glycolylneuraminic acid (​Neu5Gc) on cell surfaces, where they act as receptors for pathogens, including influenza A virus (IAV). ​Neu5Gc is synthesized from ​Neu5Ac by the enzyme cytidine monophosphate-N-acetylneuraminic acid hydroxylase (CMAH). In humans, this enzyme is inactive and only ​Neu5Ac is produced. Ferrets are susceptible to human-adapted IAV strains and have been the dominant animal model for IAV studies. Here we show that ferrets, like humans, do not synthesize ​Neu5Gc. Genomic analysis reveals an ancient, nine-exon deletion in the ferret CMAH gene that is shared by the Pinnipedia and Musteloidia members of the Carnivora. Interactions between two human strains of IAV with the sialyllactose receptor (sialic acid—α2,6Gal) confirm that the type of terminal sialic acid contributes significantly to IAV receptor specificity. Our results indicate that exclusive expression of ​Neu5Ac contributes to the susceptibility of ferrets to human-adapted IAV strains

The miR-17-5p microRNA is a key regulator of the G1/S phase cell cycle transition

Novel targets of the oncogenic miR-17-92 cluster have been identified and the mechanism of regulation of proliferation at the G1/S phase cell cycle transition via the miR-17-5p microRNA has been elucidated