Distinct Classes of Complex Structural Variation Uncovered across Thousands of Cancer Genome Graphs.

Cancer genomes often harbor hundreds of somatic DNA rearrangement junctions, many of which cannot be easily classified into simple (e.g., deletion) or complex (e.g., chromothripsis) structural variant classes. Applying a novel genome graph computational paradigm to analyze the topology of junction copy number (JCN) across 2,778 tumor whole-genome sequences, we uncovered three novel complex rearrangement phenomena: pyrgo, rigma, and tyfonas. Pyrgo are "towers" of low-JCN duplications associated with early-replicating regions, superenhancers, and breast or ovarian cancers. Rigma comprise "chasms" of low-JCN deletions enriched in late-replicating fragile sites and gastrointestinal carcinomas. Tyfonas are "typhoons" of high-JCN junctions and fold-back inversions associated with expressed protein-coding fusions, breakend hypermutation, and acral, but not cutaneous, melanomas. Clustering of tumors according to genome graph-derived features identified subgroups associated with DNA repair defects and poor prognosis.

Excess Translation of Epigenetic Regulators Contributes to Fragile X Syndrome and Is Alleviated by Brd4 Inhibition.

Fragile X syndrome (FXS) is a leading genetic cause of intellectual disability and autism. FXS results from the loss of function of fragile X mental retardation protein (FMRP), which represses translation of target transcripts. Most of the well-characterized target transcripts of FMRP are synaptic proteins, yet targeting these proteins has not provided effective treatments. We examined a group of FMRP targets that encode transcriptional regulators, particularly chromatin-associated proteins. Loss of FMRP in mice results in widespread changes in chromatin regulation and aberrant gene expression. To determine if targeting epigenetic factors could reverse phenotypes associated with the disorder, we focused on Brd4, a BET protein and chromatin reader targeted by FMRP. Inhibition of Brd4 function alleviated many of the phenotypes associated with FXS. We conclude that loss of FMRP results in significant epigenetic misregulation and that targeting transcription via epigenetic regulators like Brd4 may provide new treatments for FXS.

Genomic Patterns of De Novo Mutation in Simplex Autism.

To further our understanding of the genetic etiology of autism, we generated and analyzed genome sequence data from 516 idiopathic autism families (2,064 individuals). This resource includes >59 million single-nucleotide variants (SNVs) and 9,212 private copy number variants (CNVs), of which 133,992 and 88 are de novo mutations (DNMs), respectively. We estimate a mutation rate of ∼1.5 × 10-8 SNVs per site per generation with a significantly higher mutation rate in repetitive DNA. Comparing probands and unaffected siblings, we observe several DNM trends. Probands carry more gene-disruptive CNVs and SNVs, resulting in severe missense mutations and mapping to predicted fetal brain promoters and embryonic stem cell enhancers. These differences become more pronounced for autism genes (p = 1.8 × 10-3, OR = 2.2). Patients are more likely to carry multiple coding and noncoding DNMs in different genes, which are enriched for expression in striatal neurons (p = 3 × 10-3), suggesting a path forward for genetically characterizing more complex cases of autism.

Hepatitis C virus RNA functionally sequesters miR-122.

Hepatitis C virus (HCV) uniquely requires the liver-specific microRNA-122 for replication, yet global effects on endogenous miRNA targets during infection are unexplored. Here, high-throughput sequencing and crosslinking immunoprecipitation (HITS-CLIP) experiments of human Argonaute (AGO) during HCV infection showed robust AGO binding on the HCV 5'UTR at known and predicted miR-122 sites. On the human transcriptome, we observed reduced AGO binding and functional mRNA de-repression of miR-122 targets during virus infection. This miR-122 "sponge" effect was relieved and redirected to miR-15 targets by swapping the miRNA tropism of the virus. Single-cell expression data from reporters containing miR-122 sites showed significant de-repression during HCV infection depending on expression level and site number. We describe a quantitative mathematical model of HCV-induced miR-122 sequestration and propose that such miR-122 inhibition by HCV RNA may result in global de-repression of host miR-122 targets, providing an environment fertile for the long-term oncogenic potential of HCV.

YTHDC2 control of gametogenesis requires helicase activity but not m6A binding.

Mechanisms regulating meiotic progression in mammals are poorly understood. The N6-methyladenosine (m6A) reader and 3' → 5' RNA helicase YTHDC2 switches cells from mitotic to meiotic gene expression programs and is essential for meiotic entry, but how this critical cell fate change is accomplished is unknown. Here, we provide insight into its mechanism and implicate YTHDC2 in having a broad role in gene regulation during multiple meiotic stages. Unexpectedly, mutation of the m6A-binding pocket of YTHDC2 had no detectable effect on gametogenesis and mouse fertility, suggesting that YTHDC2 function is m6A-independent. Supporting this conclusion, CLIP data defined YTHDC2-binding sites on mRNA as U-rich and UG-rich motif-containing regions within 3' UTRs and coding sequences, distinct from the sites that contain m6A during spermatogenesis. Complete loss of YTHDC2 during meiotic entry did not substantially alter translation of its mRNA binding targets in whole-testis ribosome profiling assays but did modestly affect their steady-state levels. Mutation of the ATPase motif in the helicase domain of YTHDC2 did not affect meiotic entry, but it blocked meiotic prophase I progression, causing sterility. Our findings inform a model in which YTHDC2 binds transcripts independent of m6A status and regulates gene expression during multiple stages of meiosis by distinct mechanisms.

AGO HITS-CLIP reveals distinct miRNA regulation of white and brown adipose tissue identity.

MicroRNAs (miRNAs) are short, noncoding RNAs that associate with Argonaute (AGO) to influence mRNA stability and translation, thereby regulating cellular determination and phenotype. While several individual miRNAs have been shown to control adipocyte function, including energy storage in white fat and energy dissipation in brown fat, a comprehensive analysis of miRNA activity in these tissues has not been performed. We used high-throughput sequencing of RNA isolated by cross-linking immunoprecipitation (HITS-CLIP) to comprehensively characterize the network of high-confidence, in vivo mRNA:miRNA interactions across white and brown fat, revealing >20,000 unique AGO binding sites. When coupled with miRNA and mRNA sequencing, we found an inverse correlation between depot-enriched miRNAs and their targets. To illustrate the functionality of our HITS-CLIP data set in identifying specific miRNA:mRNA interactions, we show that miR-29 is a novel regulator of leptin, an adipocyte-derived hormone that coordinates food intake and energy homeostasis. Two independent miR-29 binding sites in the leptin 3' UTR were validated using luciferase assays, and miR-29 gain and loss of function modulated leptin mRNA and protein secretion in primary adipocytes. This work represents the only experimentally generated miRNA targetome in adipose tissue and identifies multiple regulatory pathways that may specify the unique identities of white and brown fat.

The Genetic Control of Stoichiometry Underlying Autism.

Autism is a common and complex neurologic disorder whose scientific underpinnings have begun to be established in the past decade. The essence of this breakthrough has been a focus on families, where genetic analyses are strongest, versus large-scale, case-control studies. Autism genetics has progressed in parallel with technology, from analyses of copy number variation to whole-exome sequencing (WES) and whole-genome sequencing (WGS). Gene mutations causing complete loss of function account for perhaps one-third of cases, largely detected through WES. This limitation has increased interest in understanding the regulatory variants of genes that contribute in more subtle ways to the disorder. Strategies combining biochemical analysis of gene regulation, WGS analysis of the noncoding genome, and machine learning have begun to succeed. The emerging picture is that careful control of the amounts of transcription, mRNA, and proteins made by key brain genes-stoichiometry-plays a critical role in defining the clinical features of autism.

NOVA1 acts as an oncogenic RNA-binding protein to regulate cholesterol homeostasis in human glioblastoma cells.

NOVA1 is a neuronal RNA-binding protein identified as the target antigen of a rare autoimmune disorder associated with cancer and neurological symptoms, termed paraneoplastic opsoclonus-myoclonus ataxia. Despite the strong association between NOVA1 and cancer, it has been unclear how NOVA1 function might contribute to cancer biology. In this study, we find that NOVA1 acts as an oncogenic factor in a GBM (glioblastoma multiforme) cell line established from a patient. Interestingly, NOVA1 and Argonaute (AGO) CLIP identified common 3' untranslated region (UTR) targets, which were down-regulated in NOVA1 knockdown GBM cells, indicating a transcriptome-wide intersection of NOVA1 and AGO-microRNA (miRNA) targets regulation. NOVA1 binding to 3'UTR targets stabilized transcripts including those encoding cholesterol homeostasis related proteins. Selective inhibition of NOVA1-RNA interactions with antisense oligonucleotides disrupted GBM cancer cell fitness. The precision of our GBM CLIP studies point to both mechanism and precise RNA sequence sites to selectively inhibit oncogenic NOVA1-RNA interactions. Taken together, we find that NOVA1 is commonly overexpressed in GBM, where it can antagonize AGO2-miRNA actions and consequently up-regulates cholesterol synthesis, promoting cell viability.

FMRP stalls ribosomal translocation on mRNAs linked to synaptic function and autism.

FMRP loss of function causes Fragile X syndrome (FXS) and autistic features. FMRP is a polyribosome-associated neuronal RNA-binding protein, suggesting that it plays a key role in regulating neuronal translation, but there has been little consensus regarding either its RNA targets or mechanism of action. Here, we use high-throughput sequencing of RNAs isolated by crosslinking immunoprecipitation (HITS-CLIP) to identify FMRP interactions with mouse brain polyribosomal mRNAs. FMRP interacts with the coding region of transcripts encoding pre- and postsynaptic proteins and transcripts implicated in autism spectrum disorders (ASD). We developed a brain polyribosome-programmed translation system, revealing that FMRP reversibly stalls ribosomes specifically on its target mRNAs. Our results suggest that loss of a translational brake on the synthesis of a subset of synaptic proteins contributes to FXS. In addition, they provide insight into the molecular basis of the cognitive and allied defects in FXS and ASD and suggest multiple targets for clinical intervention.

Argonaute CLIP Defines a Deregulated miR-122-Bound Transcriptome that Correlates with Patient Survival in Human Liver Cancer.

MicroRNA-122, an abundant and conserved liver-specific miRNA, regulates hepatic metabolism and functions as a tumor suppressor, yet systematic and direct biochemical elucidation of the miR-122 target network remains incomplete. To this end, we performed Argonaute crosslinking immunoprecipitation (Argonaute [Ago]-CLIP) sequencing in miR-122 knockout and control mouse livers, as well as in matched human hepatocellular carcinoma (HCC) and benign liver tissue to identify miRNA target sites transcriptome-wide in two species. We observed a majority of miR-122 binding on 3' UTRs and coding exons followed by extensive binding to other genic and non-genic sites. Motif analysis of miR-122-dependent binding revealed a G-bulged motif in addition to canonical motifs. A large number of miR-122 targets were found to be species specific. Upregulation of several common mouse and human targets, most notably BCL9, predicted survival in HCC patients. These results broadly define the molecular consequences of miR-122 downregulation in hepatocellular carcinoma.

m6A mRNA modifications are deposited in nascent pre-mRNA and are not required for splicing but do specify cytoplasmic turnover.

Understanding the biologic role of N6-methyladenosine (m6A) RNA modifications in mRNA requires an understanding of when and where in the life of a pre-mRNA transcript the modifications are made. We found that HeLa cell chromatin-associated nascent pre-mRNA (CA-RNA) contains many unspliced introns and m6A in exons but very rarely in introns. The m6A methylation is essentially completed upon the release of mRNA into the nucleoplasm. Furthermore, the content and location of each m6A modification in steady-state cytoplasmic mRNA are largely indistinguishable from those in the newly synthesized CA-RNA or nucleoplasmic mRNA. This result suggests that quantitatively little methylation or demethylation occurs in cytoplasmic mRNA. In addition, only ∼10% of m6As in CA-RNA are within 50 nucleotides of 5' or 3' splice sites, and the vast majority of exons harboring m6A in wild-type mouse stem cells is spliced the same in cells lacking the major m6A methyltransferase Mettl3. Both HeLa and mouse embryonic stem cell mRNAs harboring m6As have shorter half-lives, and thousands of these mRNAs have increased half-lives (twofold or more) in Mettl3 knockout cells compared with wild type. In summary, m6A is added to exons before or soon after exon definition in nascent pre-mRNA, and while m6A is not required for most splicing, its addition in the nascent transcript is a determinant of cytoplasmic mRNA stability.

Vaccine Breakthrough Infections with SARS-CoV-2 Variants.

Emerging variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are of clinical concern. In a cohort of 417 persons who had received the second dose of BNT162b2 (Pfizer-BioNTech) or mRNA-1273 (Moderna) vaccine at least 2 weeks previously, we identified 2 women with vaccine breakthrough infection. Despite evidence of vaccine efficacy in both women, symptoms of coronavirus disease 2019 developed, and they tested positive for SARS-CoV-2 by polymerase-chain-reaction testing. Viral sequencing revealed variants of likely clinical importance, including E484K in 1 woman and three mutations (T95I, del142-144, and D614G) in both. These observations indicate a potential risk of illness after successful vaccination and subsequent infection with variant virus, and they provide support for continued efforts to prevent and diagnose infection and to characterize variants in vaccinated persons. (Funded by the National Institutes of Health and others.).

RNA Identification of PRIME Cells Predicting Rheumatoid Arthritis Flares.

BACKGROUND: Rheumatoid arthritis, like many inflammatory diseases, is characterized by episodes of quiescence and exacerbation (flares). The molecular events leading to flares are unknown. METHODS: We established a clinical and technical protocol for repeated home collection of blood in patients with rheumatoid arthritis to allow for longitudinal RNA sequencing (RNA-seq). Specimens were obtained from 364 time points during eight flares over a period of 4 years in our index patient, as well as from 235 time points during flares in three additional patients. We identified transcripts that were differentially expressed before flares and compared these with data from synovial single-cell RNA-seq. Flow cytometry and sorted-blood-cell RNA-seq in additional patients were used to validate the findings. RESULTS: Consistent changes were observed in blood transcriptional profiles 1 to 2 weeks before a rheumatoid arthritis flare. B-cell activation was followed by expansion of circulating CD45-CD31-PDPN+ preinflammatory mesenchymal, or PRIME, cells in the blood from patients with rheumatoid arthritis; these cells shared features of inflammatory synovial fibroblasts. Levels of circulating PRIME cells decreased during flares in all 4 patients, and flow cytometry and sorted-cell RNA-seq confirmed the presence of PRIME cells in 19 additional patients with rheumatoid arthritis. CONCLUSIONS: Longitudinal genomic analysis of rheumatoid arthritis flares revealed PRIME cells in the blood during the period before a flare and suggested a model in which these cells become activated by B cells in the weeks before a flare and subsequently migrate out of the blood into the synovium. (Funded by the National Institutes of Health and others.).

NOVA2 regulates neural circRNA biogenesis.

Circular RNAs (circRNAs) are highly expressed in the brain and their expression increases during neuronal differentiation. The factors regulating circRNAs in the developing mouse brain are unknown. NOVA1 and NOVA2 are neural-enriched RNA-binding proteins with well-characterized roles in alternative splicing. Profiling of circRNAs from RNA-seq data revealed that global circRNA levels were reduced in embryonic cortex of Nova2 but not Nova1 knockout mice. Analysis of isolated inhibitory and excitatory cortical neurons lacking NOVA2 revealed an even more dramatic reduction of circRNAs and establishes a widespread role for NOVA2 in enhancing circRNA biogenesis. To investigate the cis-elements controlling NOVA2-regulation of circRNA biogenesis, we generated a backsplicing reporter based on the Efnb2 gene. We found that NOVA2-mediated backsplicing of circEfnb2 was impaired when YCAY clusters located in flanking introns were mutagenized. CLIP (cross-linking and immunoprecipitation) and additional reporter analyses demonstrated the importance of NOVA2 binding sites located in both flanking introns of circRNA loci. NOVA2 is the first RNA-binding protein identified to globally promote circRNA biogenesis in the developing brain.

A majority of m6A residues are in the last exons, allowing the potential for 3' UTR regulation.

We adapted UV CLIP (cross-linking immunoprecipitation) to accurately locate tens of thousands of m(6)A residues in mammalian mRNA with single-nucleotide resolution. More than 70% of these residues are present in the 3'-most (last) exons, with a very sharp rise (sixfold) within 150-400 nucleotides of the start of the last exon. Two-thirds of last exon m(6)A and >40% of all m(6)A in mRNA are present in 3' untranslated regions (UTRs); contrary to earlier suggestions, there is no preference for location of m(6)A sites around stop codons. Moreover, m(6)A is significantly higher in noncoding last exons than in next-to-last exons harboring stop codons. We found that m(6)A density peaks early in the 3' UTR and that, among transcripts with alternative polyA (APA) usage in both the brain and the liver, brain transcripts preferentially use distal polyA sites, as reported, and also show higher proximal m(6)A density in the last exons. Furthermore, when we reduced m6A methylation by knocking down components of the methylase complex and then examined 661 transcripts with proximal m6A peaks in last exons, we identified a set of 111 transcripts with altered (approximately two-thirds increased proximal) APA use. Taken together, these observations suggest a role of m(6)A modification in regulating proximal alternative polyA choice.

RNA protein interaction in neurons.

Neurons have their own systems for regulating RNA. Several multigene families encode RNA binding proteins (RNABPs) that are uniquely expressed in neurons, including the well-known neuron-specific markers ELAV and NeuN and the disease antigen NOVA. New technologies have emerged in recent years to assess the function of these proteins in vivo, and the answers are yielding insights into how and why neurons may regulate RNA in special ways-to increase cellular complexity, to localize messenger RNA (mRNA) spatially, and to regulate their expression in response to synaptic stimuli. The functions of such restricted neuronal proteins are likely to be complemented by more widely expressed RNABPs that may themselves have developed specialized functions in neurons, including Argonaute/microRNAs (miRNAs). Here we review what is known about such RNABPs and explore the potential biologic and neurologic significance of neuronal RNA regulatory systems.

Pre-mRNA processing includes N6 methylation of adenosine residues that are retained in mRNA exons and the fallacy of "RNA epigenetics".

By using a cell fraction technique that separates chromatin-associated nascent RNA, newly completed nucleoplasmic mRNA and cytoplasmic mRNA, we have shown in a previous study that residues in exons are methylated (m6A) in nascent pre-mRNA and remain methylated in the same exonic residues in nucleoplasmic and cytoplasmic mRNA. Thus, there is no evidence of a substantial degree of demethylation in mRNA exons that would correspond to so-called "epigenetic" demethylation. The turnover rate of mRNA molecules is faster, depending on m6A content in HeLa cell mRNA, suggesting that specification of mRNA stability may be the major role of m6A exon modification. In mouse embryonic stem cells (mESCs) lacking Mettl3, the major mRNA methylase, the cells continue to grow, making the same mRNAs with unchanged splicing profiles in the absence (>90%) of m6A in mRNA, suggesting no common obligatory role of m6A in splicing. All these data argue strongly against a commonly used "reversible dynamic methylation/demethylation" of mRNA, calling into question the concept of "RNA epigenetics" that parallels the well-established role of dynamic DNA epigenetics.

Transcriptome-wide miR-155 binding map reveals widespread noncanonical microRNA targeting.

MicroRNAs (miRNAs) are essential components of gene regulation, but identification of miRNA targets remains a major challenge. Most target prediction and discovery relies on perfect complementarity of the miRNA seed to the 3' untranslated region (UTR). However, it is unclear to what extent miRNAs target sites without seed matches. Here, we performed a transcriptome-wide identification of the endogenous targets of a single miRNA-miR-155-in a genetically controlled manner. We found that approximately 40% of miR-155-dependent Argonaute binding occurs at sites without perfect seed matches. The majority of these noncanonical sites feature extensive complementarity to the miRNA seed with one mismatch. These noncanonical sites confer regulation of gene expression, albeit less potently than canonical sites. Thus, noncanonical miRNA binding sites are widespread, often contain seed-like motifs, and can regulate gene expression, generating a continuum of targeting and regulation.

Ptbp2 represses adult-specific splicing to regulate the generation of neuronal precursors in the embryonic brain.

Two polypyrimidine tract RNA-binding proteins (PTBs), one near-ubiquitously expressed (Ptbp1) and another highly tissue-restricted (Ptbp2), regulate RNA in interrelated but incompletely understood ways. Ptbp1, a splicing regulator, is replaced in the brain and differentiated neuronal cell lines by Ptbp2. To define the roles of Ptbp2 in the nervous system, we generated two independent Ptbp2-null strains, unexpectedly revealing that Ptbp2 is expressed in neuronal progenitors and is essential for postnatal survival. A HITS-CLIP (high-throughput sequencing cross-linking immunoprecipitation)-generated map of reproducible Ptbp2-RNA interactions in the developing mouse neocortex, combined with results from splicing-sensitive microarrays, demonstrated that the major action of Ptbp2 is to inhibit adult-specific alternative exons by binding pyrimidine-rich sequences upstream of and/or within them. These regulated exons are present in mRNAs encoding proteins associated with control of cell fate, proliferation, and the actin cytoskeleton, suggesting a role for Ptbp2 in neurogenesis. Indeed, neuronal progenitors in the Ptbp2-null brain exhibited an aberrant polarity and were associated with regions of premature neurogenesis and reduced progenitor pools. Thus, Ptbp2 inhibition of a discrete set of adult neuronal exons underlies early brain development prior to neuronal differentiation and is essential for postnatal survival.