William R. Jeck

Detecting and characterizing circular RNAs

Circular RNA transcripts were first identified in the early 1990s but knowledge of these species has remained limited, as their study through traditional methods of RNA analysis has been difficult. Now, novel bioinformatic approaches coupled with biochemical enrichment strategies and deep sequencing have allowed comprehensive studies of circular RNA species. Recent studies have revealed thousands of endogenous circular RNAs in mammalian cells, some of which are highly abundant and evolutionarily conserved. Evidence is emerging that some circRNAs might regulate microRNA (miRNA) function, and roles in transcriptional control have also been suggested. Therefore, study of this class of noncoding RNAs has potential implications for therapeutic and research applications. We believe the key future challenge for the field will be to understand the regulation and function of these unusual molecules.

Expression of linear and novel circular forms of an INK4/ARF-associated non-coding RNA correlates with atherosclerosis risk.

Human genome-wide association studies have linked single nucleotide polymorphisms (SNPs) on chromosome 9p21.3 near the INK4/ARF (CDKN2a/b) locus with susceptibility to atherosclerotic vascular disease (ASVD). Although this locus encodes three well-characterized tumor suppressors, p16INK4a, p15INK4b, and ARF, the SNPs most strongly associated with ASVD are ∼120 kb from the nearest coding gene within a long non-coding RNA (ncRNA) known as ANRIL (CDKN2BAS). While individuals homozygous for the atherosclerotic risk allele show decreased expression of ANRIL and the coding INK4/ARF transcripts, the mechanism by which such distant genetic variants influence INK4/ARF expression is unknown. Here, using rapid amplification of cDNA ends (RACE) and analysis of next-generation RNA sequencing datasets, we determined the structure and abundance of multiple ANRIL species. Each of these species was present at very low copy numbers in primary and cultured cells; however, only the expression of ANRIL isoforms containing exons proximal to the INK4/ARF locus correlated with the ASVD risk alleles. Surprisingly, RACE also identified transcripts containing non-colinear ANRIL exonic sequences, whose expression also correlated with genotype and INK4/ARF expression. These non-polyadenylated RNAs resisted RNAse R digestion and could be PCR amplified using outward-facing primers, suggesting they represent circular RNA structures that could arise from by-products of mRNA splicing. Next-generation DNA sequencing and splice prediction algorithms identified polymorphisms within the ASVD risk interval that may regulate ANRIL splicing and circular ANRIL (cANRIL) production. These results identify novel circular RNA products emanating from the ANRIL locus and suggest causal variants at 9p21.3 regulate INK4/ARF expression and ASVD risk by modulating ANRIL expression and/or structure.

Extending assembly of short DNA sequences to handle error

UNLABELLED Inexpensive de novo genome sequencing, particularly in organisms with small genomes, is now possible using several new sequencing technologies. Some of these technologies such as that from Illuminas Solexa Sequencing, produce high genomic coverage by generating a very large number of small reads ( approximately 30 bp). While prior work shows that partial assembly can be performed by k-mer extension in error-free reads, this algorithm is unsuccessful with the sequencing error rates found in practice. We present VCAKE (Verified Consensus Assembly by K-mer Extension), a modification of simple k-mer extension that overcomes error by using high depth coverage. Though it is a simple modification of a previous approach, we show significant improvements in assembly results on simulated and experimental datasets that include error. AVAILABILITY http://152.2.15.114/~labweb/VCAKE

Mutations in isocitrate dehydrogenase 1 and 2 occur frequently in intrahepatic cholangiocarcinomas and share hypermethylation targets with glioblastomas

Mutations in the genes encoding isocitrate dehydrogenase, IDH1 and IDH2, have been reported in gliomas, myeloid leukemias, chondrosarcomas and thyroid cancer. We discovered IDH1 and IDH2 mutations in 34 of 326 (10%) intrahepatic cholangiocarcinomas. Tumor with mutations in IDH1 or IDH2 had lower 5-hydroxymethylcytosine and higher 5-methylcytosine levels, as well as increased dimethylation of histone H3 lysine 79 (H3K79). Mutations in IDH1 or IDH2 were associated with longer overall survival (P=0.028) and were independently associated with a longer time to tumor recurrence after intrahepatic cholangiocarcinoma resection in multivariate analysis (P=0.021). IDH1 and IDH2 mutations were significantly associated with increased levels of p53 in intrahepatic cholangiocarcinomas, but no mutations in the p53 gene were found, suggesting that mutations in IDH1 and IDH2 may cause a stress that leads to p53 activation. We identified 2309 genes that were significantly hypermethylated in 19 cholangiocarcinomas with mutations in IDH1 or IDH2, compared with cholangiocarcinomas without these mutations. Hypermethylated CpG sites were significantly enriched in CpG shores and upstream of transcription start sites, suggesting a global regulation of transcriptional potential. Half of the hypermethylated genes overlapped with DNA hypermethylation in IDH1-mutant gliobastomas, suggesting the existence of a common set of genes whose expression may be affected by mutations in IDH1 or IDH2 in different types of tumors.

De novo assembly using low-coverage short read sequence data from the rice pathogen Pseudomonas syringae pv. oryzae

We developed a novel approach for de novo genome assembly using only sequence data from high-throughput short read sequencing technologies. By combining data generated from 454 Life Sciences (Roche) and Illumina (formerly known as Solexa sequencing) sequencing platforms, we reliably assembled genomes into large scaffolds at a fraction of the traditional cost and without use of a reference sequence. We applied this method to two isolates of the phytopathogenic bacteria Pseudomonas syringae. Sequencing and reassembly of the well-studied tomato and Arabidopsis pathogen, Pto(DC3000), facilitated development and testing of our method. Sequencing of a distantly related rice pathogen, Por(1_)(6), demonstrated our methods efficacy for de novo assembly of novel genomes. Our assembly of Por(1_6) yielded an N50 scaffold size of 531,821 bp with >75% of the predicted genome covered by scaffolds over 100,000 bp. One of the critical phenotypic differences between strains of P. syringae is the range of plant hosts they infect. This is largely determined by their complement of type III effector proteins. The genome of Por(1_6) is the first sequenced for a P. syringae isolate that is a pathogen of monocots, and, as might be predicted, its complement of type III effectors differs substantially from the previously sequenced isolates of this species. The genome of Por(1_6) helps to define an expansion of the P. syringae pan-genome, a corresponding contraction of the core genome, and a further diversification of the type III effector complement for this important plant pathogen species.

Identification of Driver Genes in Hepatocellular Carcinoma by Exome Sequencing

Genetic alterations in specific driver genes lead to disruption of cellular pathways and are critical events in the instigation and progression of hepatocellular carcinoma (HCC). As a prerequisite for individualized cancer treatment, we sought to characterize the landscape of recurrent somatic mutations in HCC. We performed whole‐exome sequencing on 87 HCCs and matched normal adjacent tissues to an average coverage of 59×. The overall mutation rate was roughly two mutations per Mb, with a median of 45 nonsynonymous mutations that altered the amino acid sequence (range, 2‐381). We found recurrent mutations in several genes with high transcript levels: TP53 (18%); CTNNB1 (10%); KEAP1 (8%); C16orf62 (8%); MLL4 (7%); and RAC2 (5%). Significantly affected gene families include the nucleotide‐binding domain and leucine‐rich repeat‐containing family, calcium channel subunits, and histone methyltransferases. In particular, the MLL family of methyltransferases for histone H3 lysine 4 were mutated in 20% of tumors. Conclusion: The NFE2L2‐KEAP1 and MLL pathways are recurrently mutated in multiple cohorts of HCC. (Hepatology 2013;58:1693–1702)

Review: a meta-analysis of GWAS and age-associated diseases

Genome‐Wide Association studies (GWAS) offer an unbiased means to understand the genetic basis of traits by identifying single nucleotide polymorphisms (SNPs) linked to causal variants of complex phenotypes. GWAS have identified a host of susceptibility SNPs associated with many important human diseases, including diseases associated with aging. In an effort to understand the genetics of broad resistance to age‐associated diseases (i.e., ‘wellness’), we performed a meta‐analysis of human GWAS. Toward that end, we compiled 372 GWAS that identified 1775 susceptibility SNPs to 105 unique diseases and used these SNPs to create a genomic landscape of disease susceptibility. This map was constructed by partitioning the genome into 200 kb ‘bins’ and mapping the 1775 susceptibility SNPs to bins based on their genomic location. Investigation of these data revealed significant heterogeneity of disease association within the genome, with 92% of bins devoid of disease‐associated SNPs. In contrast, 10 bins (0.06%) were significantly (P < 0.05) enriched for susceptibility to multiple diseases, 5 of which formed two highly significant peaks of disease association (P < 0.0001). These peaks mapped to the Major Histocompatibility (MHC) locus on 6p21 and the INK4/ARF (CDKN2a/b) tumor suppressor locus on 9p21.3. Provocatively, all 10 significantly enriched bins contained genes linked to either inflammation or cellular senescence pathways, and SNPs near regulators of senescence were particularly associated with disease of aging (e.g., cancer, atherosclerosis, type 2 diabetes, glaucoma). This analysis suggests that germline genetic heterogeneity in the regulation of immunity and cellular senescence influences the human healthspan.

An oncogenic Ezh2 mutation induces tumors through global redistribution of histone 3 lysine 27 trimethylation.

B cell lymphoma and melanoma harbor recurrent mutations in the gene encoding the EZH2 histone methyltransferase (EZH2), but the carcinogenic role of these mutations is unclear. Here we describe a mouse model in which the most common somatic Ezh2 gain-of-function mutation (EZH2Y646F in human; Ezh2Y641F in mouse) is conditionally expressed. Expression of Ezh2Y641F in mouse B cells or melanocytes caused high-penetrance lymphoma or melanoma, respectively. Overexpression of the anti-apoptotic protein Bcl2, but not the oncoprotein Myc, or loss of the tumor suppressor protein p53 (encoded by Trp53 in mice) further accelerated lymphoma progression. Expression of the mutant Braf but not the mutant Nras oncoprotein further accelerated melanoma progression. Although expression of Ezh2Y641F globally increased the abundance of trimethylated Lys27 of histone H3 (H3K27me3), it also caused a widespread redistribution of this repressive mark, including a loss of H3K27me3 that was associated with increased transcription at many loci. These results suggest that Ezh2Y641F induces lymphoma and melanoma through a vast reorganization of chromatin structure, inducing both repression and activation of polycomb-regulated loci.B-cell lymphoma and melanoma harbor recurrent mutations in the gene encoding the EZH2 histone methyltransferase, but the carcinogenic role of these mutations is unclear. Here we describe a mouse model in which the most common somatic EZH2 gain-of-function mutation (Y646F in human, Y641F in the mouse) can be conditionally expressed. Expression of Ezh2Y641F in mouse B-cells or melanocytes caused high-penetrance lymphoma or melanoma, respectively. Bcl2 overexpression or p53 loss, but not c-Myc overexpression, further accelerated lymphoma progression, and expression of mutant B-Raf but not mutant N-Ras further accelerated melanoma progression. Although expression of Ezh2Y641F increased abundance of global H3K27 trimethylation (H3K27me3), it also caused a widespread redistribution of this repressive mark, including a loss of H3K27me3 associated with increased transcription at many loci. These results suggest that Ezh2Y641F induces lymphoma and melanoma through a vast reorganization of chromatin structure inducing both repression and activation of polycomb-regulated loci.

Targeted Next Generation Sequencing Identifies Clinically Actionable Mutations in Patients with Melanoma

Somatic sequencing of cancers has produced new insight into tumorigenesis, tumor heterogeneity, and disease progression, but the vast majority of genetic events identified are of indeterminate clinical significance. Here, we describe a NextGen sequencing approach to fully analyzing 248 genes, including all those of known clinical significance in melanoma. This strategy features solution capture of DNA followed by multiplexed, high‐throughput sequencing and was evaluated in 31 melanoma cell lines and 18 tumor tissues from patients with metastatic melanoma. Mutations in melanoma cell lines correlated with their sensitivity to corresponding small molecule inhibitors, confirming, for example, lapatinib sensitivity in ERBB4 mutant lines and identifying a novel activating mutation of BRAF. The latter event would not have been identified by clinical sequencing and was associated with responsiveness to a BRAF kinase inhibitor. This approach identified focal copy number changes of PTEN not found by standard methods, such as comparative genomic hybridization (CGH). Actionable mutations were found in 89% of the tumor tissues analyzed, 56% of which would not be identified by standard‐of‐care approaches. This work shows that targeted sequencing is an attractive approach for clinical use in melanoma.

Subtyping sub-Saharan esophageal squamous cell carcinoma by comprehensive molecular analysis

Esophageal squamous cell carcinoma (ESCC) is endemic in regions of sub-Saharan Africa (SSA), where it is the third most common cancer. Here, we describe whole-exome tumor/normal sequencing and RNA transcriptomic analysis of 59 patients with ESCC in Malawi. We observed similar genetic aberrations as reported in Asian and North American cohorts, including mutations of TP53, CDKN2A, NFE2L2, CHEK2, NOTCH1, FAT1, and FBXW7. Analyses for nonhuman sequences did not reveal evidence for infection with HPV or other occult pathogens. Mutational signature analysis revealed common signatures associated with aging, cytidine deaminase activity (APOBEC), and a third signature of unknown origin, but signatures of inhaled tobacco use, aflatoxin and mismatch repair were notably absent. Based on RNA expression analysis, ESCC could be divided into 3 distinct subtypes, which were distinguished by their expression of cell cycle and neural transcripts. This study demonstrates discrete subtypes of ESCC in SSA, and suggests that the endemic nature of this disease reflects exposure to a carcinogen other than tobacco and oncogenic viruses.