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 has been difficult through traditional methods of RNA analysis. 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 (circRNAs) 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 non-coding RNAs has potential implications for therapeutic and research applications. We believe the key future challenge to the field will be to understand the regulation and function of these unusual molecules

Review: A Meta-Analysis of GWAS Studies 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 1,775 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 1,775 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 (pINK4/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 health span

BET Protein Inhibition Regulates Macrophage Chromatin Accessibility and Microbiota-Dependent Colitis

Introduction In colitis, macrophage functionality is altered compared to normal homeostatic conditions. Loss of IL-10 signaling results in an inappropriate chronic inflammatory response to bacterial stimulation. It remains unknown if inhibition of bromodomain and extra-terminal domain (BET) proteins alters usage of DNA regulatory elements responsible for driving inflammatory gene expression. We determined if the BET inhibitor, (+)-JQ1, could suppress inflammatory activation of macrophages in Il10-/- mice. Methods We performed ATAC-seq and RNA-seq on Il10-/- bone marrow-derived macrophages (BMDMs) cultured in the presence and absence of lipopolysaccharide (LPS) with and without treatment with (+)-JQ1 and evaluated changes in chromatin accessibility and gene expression. Germ-free Il10-/- mice were treated with (+)-JQ1, colonized with fecal slurries and underwent histological and molecular evaluation 14-days post colonization. Results Treatment with (+)-JQ1 suppressed LPS-induced changes in chromatin at distal regulatory elements associated with inflammatory genes, particularly in regions that contain motifs for AP-1 and IRF transcription factors. This resulted in attenuation of inflammatory gene expression. Treatment with (+)-JQ1 in vivo resulted in a mild reduction in colitis severity as compared with vehicle-treated mice. Conclusion We identified the mechanism of action associated with a new class of compounds that may mitigate aberrant macrophage responses to bacteria in colitis

An oncogenic Ezh2 mutation cooperates with particular genetic alterations to induce tumors in mice and redistributes H3K27 trimethylation throughout the genome

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

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, 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 analyze 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

Mutation-Specific RAS Oncogenicity Explains NRAS Codon 61 Selection in Melanoma

N-RAS mutation at codon 12, 13 or 61 is associated with transformation; yet, in melanoma, such alterations are nearly exclusive to codon 61. Here, we compared the melanoma susceptibility of an N-RasQ61R knock-in allele to similarly designed K-RasG12D and N-RasG12D alleles. With concomitant p16INK4a inactivation, K-RasG12D or N-RasQ61R expression efficiently promoted melanoma in vivo, whereas N-RasG12D did not. Additionally, N-RasQ61R mutation potently cooperated with Lkb1/Stk11 loss to induce highly metastatic disease. Functional comparisons of N-RasQ61R and N-RasG12D revealed little difference in the ability of these proteins to engage PI3K or RAF. Instead, N-RasQ61R showed enhanced nucleotide binding, decreased intrinsic GTPase activity and increased stability when compared to N-RasG12D. This work identifies a faithful model of human N-RAS mutant melanoma, and suggests that the increased melanomagenecity of N-RasQ61R over N-RasG12D is due to heightened abundance of the active, GTP-bound form rather than differences in the engagement of downstream effector pathways

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

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, PtoDC3000, facilitated development and testing of our method. Sequencing of a distantly related rice pathogen, Por1_6, demonstrated our method's efficacy for de novo assembly of novel genomes. Our assembly of Por1_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 Por1_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 Por1_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