Thousands of loci in the human and mouse genomes give rise to circular RNA transcripts; at many of these loci, the predominant RNA isoform is a circle. Using an improved computational approach for circular RNA identification, we found widespread circular RNA expression in Drosophila melanogaster and estimate that in humans, circular RNA may account for 1% as many molecules as poly(A) RNA. Analysis of data from the ENCODE consortium revealed that the repertoire of genes expressing circular RNA, the ratio of circular to linear transcripts for each gene, and even the pattern of splice isoforms of circular RNAs from each gene were cell-type specific. These results suggest that biogenesis of circular RNA is an integral, conserved, and regulated feature of the gene expression program.

An unexpectedly large fraction of genes in metazoans (human, mouse, zebrafish, worm, fruit fly) express high levels of circularized RNAs containing canonical exons. Here we report that circular RNA isoforms are found in diverse species whose most recent common ancestor existed more than one billion years ago: fungi (Schizosaccharomyces pombe and Saccharomyces cerevisiae), a plant (Arabidopsis thaliana), and protists (Plasmodium falciparum and Dictyostelium discoideum). For all species studied to date, including those in this report, only a small fraction of the theoretically possible circular RNA isoforms from a given gene are actually observed. Unlike metazoans, Arabidopsis, D. discoideum, P. falciparum, S. cerevisiae, and S. pombe have very short introns (∼100 nucleotides or shorter), yet they still produce circular RNAs. A minority of genes in S. pombe and P. falciparum have documented examples of canonical alternative splicing, making it unlikely that all circular RNAs are by-products of alternative splicing or ‘piggyback’ on signals used in alternative RNA processing. In S. pombe, the relative abundance of circular to linear transcript isoforms changed in a gene-specific pattern during nitrogen starvation. Circular RNA may be an ancient, conserved feature of eukaryotic gene expression programs.

Pancreatic cancer is the fifth leading cause of cancer death in the United States. We used cDNA microarrays to analyze global gene expression patterns in 14 pancreatic cancer cell lines, 17 resected infiltrating pancreatic cancer tissues, and 5 samples of normal pancreas to identify genes that are differentially expressed in pancreatic cancer. We found more than 400 cDNAs corresponding to genes that were differentially expressed in the pancreatic cancer tissues and cell lines as compared to normal pancreas. These genes that tended to be expressed at higher levels in pancreatic cancers were associated with a variety of processes, including cell-cell and cell-matrix interactions, cytoskeletal remodeling, proteolytic activity, and Ca(++) homeostasis. Two prominent clusters of genes were related to the high rates of cellular proliferation in pancreatic cancer cell lines and the host desmoplastic response in the resected pancreatic cancer tissues. Of 149 genes identified as more highly expressed in the pancreatic cancers compared with normal pancreas, 103 genes have not been previously reported in association with pancreatic cancer. The expression patterns of 14 of these highly expressed genes were validated by either immunohistochemistry or reverse transcriptase-polymerase chain reaction as being expressed in pancreatic cancer. The overexpression of one gene in particular, 14-3-3 sigma, was found to be associated with aberrant hypomethylation in the majority of pancreatic cancers analyzed. The genes and expressed sequence tags presented in this study provide clues to the pathobiology of pancreatic cancer and implicate a large number of potentially new molecular markers for the detection and treatment of pancreatic cancer.