A quantitative atlas of polyadenylation in five mammals.

We developed PolyA-seq, a strand-specific and quantitative method for high-throughput sequencing of 3' ends of polyadenylated transcripts, and used it to globally map polyadenylation (polyA) sites in 24 matched tissues in human, rhesus, dog, mouse, and rat. We show that PolyA-seq is as accurate as existing RNA sequencing (RNA-seq) approaches for digital gene expression (DGE), enabling simultaneous mapping of polyA sites and quantitative measurement of their usage. In human, we confirmed 158,533 known sites and discovered 280,857 novel sites (FDR < 2.5%). On average 10% of novel human sites were also detected in matched tissues in other species. Most novel sites represent uncharacterized alternative polyA events and extensions of known transcripts in human and mouse, but primarily delineate novel transcripts in the other three species. A total of 69.1% of known human genes that we detected have multiple polyA sites in their 3'UTRs, with 49.3% having three or more. We also detected polyadenylation of noncoding and antisense transcripts, including constitutive and tissue-specific primary microRNAs. The canonical polyA signal was strongly enriched and positionally conserved in all species. In general, usage of polyA sites is more similar within the same tissues across different species than within a species. These quantitative maps of polyA usage in evolutionarily and functionally related samples constitute a resource for understanding the regulatory mechanisms underlying alternative polyadenylation.

Microarray analysis shows that some microRNAs downregulate large numbers of target mRNAs.

MicroRNAs (miRNAs) are a class of noncoding RNAs that post-transcriptionally regulate gene expression in plants and animals. To investigate the influence of miRNAs on transcript levels, we transfected miRNAs into human cells and used microarrays to examine changes in the messenger RNA profile. Here we show that delivering miR-124 causes the expression profile to shift towards that of brain, the organ in which miR-124 is preferentially expressed, whereas delivering miR-1 shifts the profile towards that of muscle, where miR-1 is preferentially expressed. In each case, about 100 messages were downregulated after 12 h. The 3' untranslated regions of these messages had a significant propensity to pair to the 5' region of the miRNA, as expected if many of these messages are the direct targets of the miRNAs. Our results suggest that metazoan miRNAs can reduce the levels of many of their target transcripts, not just the amount of protein deriving from these transcripts. Moreover, miR-1 and miR-124, and presumably other tissue-specific miRNAs, seem to downregulate a far greater number of targets than previously appreciated, thereby helping to define tissue-specific gene expression in humans.

Genetic validation of whole-transcriptome sequencing for mapping expression affected by cis-regulatory variation.

BACKGROUND: Identifying associations between genotypes and gene expression levels using microarrays has enabled systematic interrogation of regulatory variation underlying complex phenotypes. This approach has vast potential for functional characterization of disease states, but its prohibitive cost, given hundreds to thousands of individual samples from populations have to be genotyped and expression profiled, has limited its widespread application. RESULTS: Here we demonstrate that genomic regions with allele-specific expression (ASE) detected by sequencing cDNA are highly enriched for cis-acting expression quantitative trait loci (cis-eQTL) identified by profiling of 500 animals in parallel, with up to 90% agreement on the allele that is preferentially expressed. We also observed widespread noncoding and antisense ASE and identified several allele-specific alternative splicing variants. CONCLUSION: Monitoring ASE by sequencing cDNA from as little as one sample is a practical alternative to expression genetics for mapping cis-acting variation that regulates RNA transcription and processing.

Plasma MicroRNAs as sensitive and specific biomarkers of tissue injury.

BACKGROUND: MicroRNAs (miRNAs) are endogenous, small noncoding RNAs. Because of their size, abundance, tissue specificity, and relative stability in plasma, miRNAs hold promise as unique accessible biomarkers to monitor tissue injury. METHODS: We investigated the use of liver-, muscle- and brain-specific miRNAs as circulating biomarkers of tissue injury. We used a highly sensitive quantitative PCR assay to measure specific miRNAs (miR-122, miR-133a, and miR-124) in plasma samples from rats treated with liver or muscle toxicants and from a rat surgical model of stroke. RESULTS: We observed increases in plasma concentrations of miR-122, miR-133a, and miR-124 corresponding to injuries in liver, muscle, and brain, respectively. miR-122 and miR-133a illustrated specificity for liver and muscle toxicity, respectively, because they were not detectable in the plasma of animals with toxicity to the other organ. This result contrasted with the results for alanine aminotransferase (ALT) and aspartate aminotransferase, which were both increased with either organ toxicity. Furthermore, miR-122 exhibited a diagnostic sensitivity superior to that of ALT when the results were correlated to the liver histopathologic results. The miR-124 concentration increased in the plasma of rats 8 h after surgery to produce brain injury and peaked at 24 h, while the miR-122 and miR-133a concentrations remained at baseline values. CONCLUSIONS: These results demonstrate that tissue-specific miRNAs may serve as diagnostically sensitive plasma biomarkers of tissue injury.

Evolutionarily conserved and diverged alternative splicing events show different expression and functional profiles.

To better decipher the functional impact of alternative splicing, we classified alternative splicing events in 10,818 pairs of human and mouse genes based on conservation at genome and transcript levels. Expression levels of conserved alternative splices in human and mouse expressed sequence tag databases show strong correlation, indicating that alternative splicing is similarly regulated in both species. A total of 43% (8921) of mouse alternative splices could be found in the human genome but not in human transcripts. Five of eleven tested mouse predictions were observed in human tissues, demonstrating that mouse transcripts provide a valuable resource for identifying alternative splicing events in human genes. Combining gene-specific measures of conserved and diverged alternative splicing with both gene classification based on Gene Ontology (GO) and microarray-determined gene expression in 52 diverse human tissues and cell lines, we found conserved alternative splicing most enriched in brain-expressed signaling pathways. Diverged alternative splicing is more prevalent in testis and cancerous cell line up-regulated processes, including protein biosynthesis, responses to stress and responses to endogenous stimuli. Using conservation as a surrogate for functional significance, these results suggest that alternative splicing plays an important role in enhancing the functional capacity of central nervous systems, while non-functional splicing more frequently occurs in testis and cell lines, possibly as a result of cellular stress and rapid proliferation.

Discovery of biomarkers predictive of GSI response in triple-negative breast cancer and adenoid cystic carcinoma.

UNLABELLED: Next-generation sequencing was used to identify Notch mutations in a large collection of diverse solid tumors. NOTCH1 and NOTCH2 rearrangements leading to constitutive receptor activation were confined to triple-negative breast cancers (TNBC; 6 of 66 tumors). TNBC cell lines with NOTCH1 rearrangements associated with high levels of activated NOTCH1 (N1-ICD) were sensitive to the gamma-secretase inhibitor (GSI) MRK-003, both alone and in combination with paclitaxel, in vitro and in vivo, whereas cell lines with NOTCH2 rearrangements were resistant to GSI. Immunohistochemical staining of N1-ICD in TNBC xenografts correlated with responsiveness, and expression levels of the direct Notch target gene HES4 correlated with outcome in patients with TNBC. Activating NOTCH1 point mutations were also identified in other solid tumors, including adenoid cystic carcinoma (ACC). Notably, ACC primary tumor xenografts with activating NOTCH1 mutations and high N1-ICD levels were sensitive to GSI, whereas N1-ICD-low tumors without NOTCH1 mutations were resistant. SIGNIFICANCE: NOTCH1 mutations, immunohistochemical staining for activated NOTCH1, and HES4 expression are biomarkers that can be used to identify solid tumors that are likely to respond to GSI-based therapies.

Circulating miR-122 as a potential biomarker of liver disease.

AIM: Expression profiles indicate that miR-122 is specifically and abundantly expressed in liver. This study sought to determine miR-122 plasma concentrations in 15 apparently healthy subjects and 30 patients with liver disease, and clarify whether plasma miR-122 correlates with ALT. MATERIALS & METHODS: miR-122 was measured by quantitative PCR in healthy volunteers and patients with liver disease. RESULTS: ALT was increased in two out of 15 (13%) apparently healthy subjects and 17 out of 30 (57%) liver disease patients. In healthy subjects, median miR-122 plasma concentration was 51.7 copies/20 pg RNA (range 16.0-312.0). In liver disease patients, median miR-122 was significantly elevated to 202.3 copies/20 pg RNA (range 20.9-1160.0; Mann-Whitney test between median concentrations; p = 0.0016). CONCLUSION: This small proof-of-principle study suggests that miR-122 may be a potential plasma biomarker of liver damage.

ACC2 is expressed at high levels in human white adipose and has an isoform with a novel N-terminus [corrected].

Acetyl-CoA carboxylases ACC1 and ACC2 catalyze the carboxylation of acetyl-CoA to malonyl-CoA, regulating fatty-acid synthesis and oxidation, and are potential targets for treatment of metabolic syndrome. Expression of ACC1 in rodent lipogenic tissues and ACC2 in rodent oxidative tissues, coupled with the predicted localization of ACC2 to the mitochondrial membrane, have suggested separate functional roles for ACC1 in lipogenesis and ACC2 in fatty acid oxidation. We find, however, that human adipose tissue, unlike rodent adipose, expresses more ACC2 mRNA relative to the oxidative tissues muscle and heart. Human adipose, along with human liver, expresses more ACC2 than ACC1. Using RT-PCR, real-time PCR, and immunoprecipitation we report a novel isoform of ACC2 (ACC2.v2) that is expressed at significant levels in human adipose. The protein generated by this isoform has enzymatic activity, is endogenously expressed in adipose, and lacks the N-terminal sequence. Both ACC2 isoforms are capable of de novo lipogenesis, suggesting that ACC2, in addition to ACC1, may play a role in lipogenesis. The results demonstrate a significant difference in ACC expression between human and rodents, which may introduce difficulties for the use of rodent models for development of ACC inhibitors.

MicroRNA targeting specificity in mammals: determinants beyond seed pairing.

Mammalian microRNAs (miRNAs) pair to 3'UTRs of mRNAs to direct their posttranscriptional repression. Important for target recognition are approximately 7 nt sites that match the seed region of the miRNA. However, these seed matches are not always sufficient for repression, indicating that other characteristics help specify targeting. By combining computational and experimental approaches, we uncovered five general features of site context that boost site efficacy: AU-rich nucleotide composition near the site, proximity to sites for coexpressed miRNAs (which leads to cooperative action), proximity to residues pairing to miRNA nucleotides 13-16, positioning within the 3'UTR at least 15 nt from the stop codon, and positioning away from the center of long UTRs. A model combining these context determinants quantitatively predicts site performance both for exogenously added miRNAs and for endogenous miRNA-message interactions. Because it predicts site efficacy without recourse to evolutionary conservation, the model also identifies effective nonconserved sites and siRNA off-targets.

PDE4B5, a novel, super-short, brain-specific cAMP phosphodiesterase-4 variant whose isoform-specifying N-terminal region is identical to that of cAMP phosphodiesterase-4D6 (PDE4D6).

The cAMP-specific phosphodiesterase-4 (PDE4) gene family is the target of several potential selective therapeutic inhibitors. The four PDE4 genes generate several distinct protein-coding isoforms through the use of alternative promoters and 5'-coding exons. Using mouse transcripts, we identified a novel, super-short isoform of human PDE4B encoding a novel 5' terminus, which we label PDE4B5. The protein-coding region of the novel 5' exon is conserved across vertebrates, chicken, zebrafish, and fugu. Reverse-transcription-polymerase chain reaction (PCR) and quantitative (PCR) measurements show that this isoform is brain-specific. The novel protein is 58 +/- 2 kDa; it has cAMP hydrolyzing enzymatic activity and is inhibited by PDE4-selective inhibitors rolipram and cilomilast (Ariflo). Confocal and subcellular fractionation analyses show that it is distributed predominantly and unevenly within the cytosol. The 16 novel N-terminal residues of PDE4B5 are identical to the 16 N-terminal residues of the super-short isoform of PDE4D (PDE4D6), which is also brain-specific. PDE4B5 is able to bind the scaffold protein DISC1, whose gene has been linked to schizophrenia. Microarray expression profiling of the PDE4 gene family shows that specific PDE4 genes are enriched in muscle and blood fractions; however, only by monitoring the individual isoforms is the brain specificity of the super-short PDE4D and PDE4B isoforms revealed. Understanding the distinct tissue specificity of PDE4 isoforms will be important for understanding phosphodiesterase biology and opportunities for therapeutic intervention.

Expression of alternatively spliced sodium channel alpha-subunit genes. Unique splicing patterns are observed in dorsal root ganglia.

Molecular medicine requires the precise definition of drug targets, and tools are now in place to provide genome-wide information on the expression and alternative splicing patterns of any known gene. DNA microarrays were used to monitor transcript levels of the nine well-characterized alpha-subunit sodium channel genes across a broad range of tissues from cynomolgus monkey, a non-human primate model. Alternative splicing of human transcripts for a subset of the genes that are expressed in dorsal root ganglia, SCN8A (Na(v)1.6), SCN9A (Na(v)1.7), and SCN11A (Na(v)1.9) was characterized in detail. Genomic sequence analysis among gene family paralogs and between cross-species orthologs suggested specific alternative splicing events within transcripts of these genes, all of which were experimentally confirmed in human tissues. Quantitative PCR revealed that certain alternative splice events are uniquely expressed in dorsal root ganglia. In addition to characterization of human transcripts, alternatively spliced sodium channel transcripts were monitored in a rat model for neuropathic pain. Consistent down-regulation of all transcripts was observed, as well as significant changes in the splicing patterns of SCN8A and SCN9A.

Parallel analysis of tagged deletion mutants efficiently identifies genes involved in endoplasmic reticulum biogenesis.

Increased levels of HMG-CoA reductase induce cell type- and isozyme-specific proliferation of the endoplasmic reticulum. In yeast, the ER proliferations induced by Hmg1p consist of nuclear-associated stacks of smooth ER membranes known as karmellae. To identify genes required for karmellae assembly, we compared the composition of populations of homozygous diploid S. cerevisiae deletion mutants following 20 generations of growth with and without karmellae. Using an initial population of 1,557 deletion mutants, 120 potential mutants were identified as a result of three independent experiments. Each experiment produced a largely non-overlapping set of potential mutants, suggesting that differences in specific growth conditions could be used to maximize the comprehensiveness of similar parallel analysis screens. Only two genes, UBC7 and YAL011W, were identified in all three experiments. Subsequent analysis of individual mutant strains confirmed that each experiment was identifying valid mutations, based on the mutant's sensitivity to elevated HMG-CoA reductase and inability to assemble normal karmellae. The largest class of HMG-CoA reductase-sensitive mutations was a subset of genes that are involved in chromatin structure and transcriptional regulation, suggesting that karmellae assembly requires changes in transcription or that the presence of karmellae may interfere with normal transcriptional regulation.

Genome-wide survey of human alternative pre-mRNA splicing with exon junction microarrays.

Alternative pre-messenger RNA (pre-mRNA) splicing plays important roles in development, physiology, and disease, and more than half of human genes are alternatively spliced. To understand the biological roles and regulation of alternative splicing across different tissues and stages of development, systematic methods are needed. Here, we demonstrate the use of microarrays to monitor splicing at every exon-exon junction in more than 10,000 multi-exon human genes in 52 tissues and cell lines. These genome-wide data provide experimental evidence and tissue distributions for thousands of known and novel alternative splicing events. Adding to previous studies, the results indicate that at least 74% of human multi-exon genes are alternatively spliced.

Discovering modes of action for therapeutic compounds using a genome-wide screen of yeast heterozygotes.

Modern medicine faces the challenge of developing safer and more effective therapies to treat human diseases. Many drugs currently in use were discovered without knowledge of their underlying molecular mechanisms. Understanding their biological targets and modes of action will be essential to design improved second-generation compounds. Here, we describe the use of a genome-wide pool of tagged heterozygotes to assess the cellular effects of 78 compounds in Saccharomyces cerevisiae. Specifically, lanosterol synthase in the sterol biosynthetic pathway was identified as a target of the antianginal drug molsidomine, which may explain its cholesterol-lowering effects. Further, the rRNA processing exosome was identified as a potential target of the cell growth inhibitor 5-fluorouracil. This genome-wide screen validated previously characterized targets or helped identify potentially new modes of action for over half of the compounds tested, providing proof of this principle for analyzing the modes of action of clinically relevant compounds.

A comprehensive transcript index of the human genome generated using microarrays and computational approaches.

BACKGROUND: Computational and microarray-based experimental approaches were used to generate a comprehensive transcript index for the human genome. Oligonucleotide probes designed from approximately 50,000 known and predicted transcript sequences from the human genome were used to survey transcription from a diverse set of 60 tissues and cell lines using ink-jet microarrays. Further, expression activity over at least six conditions was more generally assessed using genomic tiling arrays consisting of probes tiled through a repeat-masked version of the genomic sequence making up chromosomes 20 and 22. RESULTS: The combination of microarray data with extensive genome annotations resulted in a set of 28,456 experimentally supported transcripts. This set of high-confidence transcripts represents the first experimentally driven annotation of the human genome. In addition, the results from genomic tiling suggest that a large amount of transcription exists outside of annotated regions of the genome and serves as an example of how this activity could be measured on a genome-wide scale. CONCLUSIONS: These data represent one of the most comprehensive assessments of transcriptional activity in the human genome and provide an atlas of human gene expression over a unique set of gene predictions. Before the annotation of the human genome is considered complete, however, the previously unannotated transcriptional activity throughout the genome must be fully characterized.