Multiplexed transcriptome discovery of RNA-binding protein binding sites by antibody-barcode eCLIP.

Ultraviolet crosslinking and immunoprecipitation (CLIP) methodologies enable the identification of RNA binding sites of RNA-binding proteins (RBPs). Despite improvements in the library preparation of RNA fragments, the enhanced CLIP (eCLIP) protocol requires 4 days of hands-on time and lacks the ability to process several RBPs in parallel. We present a new method termed antibody-barcode eCLIP that utilizes DNA-barcoded antibodies and proximity ligation of the DNA oligonucleotides to RBP-protected RNA fragments to interrogate several RBPs simultaneously. We observe performance comparable with that of eCLIP with the advantage of dramatically increased scaling while maintaining the same material requirement of a single eCLIP experiment.

Nuclear receptor-induced chromosomal proximity and DNA breaks underlie specific translocations in cancer.

Chromosomal translocations are a hallmark of leukemia/lymphoma and also appear in solid tumors, but the underlying mechanism remains elusive. By establishing a cellular model that mimics the relative frequency of authentic translocation events without proliferation selection, we report mechanisms of nuclear receptor-dependent tumor translocations. Intronic binding of liganded androgen receptor (AR) first juxtaposes translocation loci by triggering intra- and interchromosomal interactions. AR then promotes site-specific DNA double-stranded breaks (DSBs) at translocation loci by recruiting two types of enzymatic activities induced by genotoxic stress and liganded AR, including activation-induced cytidine deaminase and the LINE-1 repeat-encoded ORF2 endonuclease. These enzymes synergistically generate site-selective DSBs at juxtaposed translocation loci that are ligated by nonhomologous end joining pathway for specific translocations. Our data suggest that the confluence of two parallel pathways initiated by liganded nuclear receptor and genotoxic stress underlies nonrandom tumor translocations, which may function in many types of tumors and pathological processes.

Nuclear receptor-enhanced transcription requires motor- and LSD1-dependent gene networking in interchromatin granules.

While the transcriptional machinery has been extensively dissected at the molecular level, little is known about regulation of chromosomal organization in the three-dimensional space of the nucleus to achieve integrated transcriptional responses to diverse signaling events. Here, we report that ligand induces rapid interchromosomal interactions among subsets of estrogen receptor alpha-bound transcription units, with a dramatic reorganization of nuclear territories requiring nuclear actin/myosin-I transport machinery, dynein light chain 1 (DLC1), and a specific subset of transcriptional coactivators and chromatin remodeling complexes. We establish a requirement for the histone lysine demethylase, LSD1, in directing specific interchromosomal interaction loci to distinct interchromatin granules, long thought to be "storage" sites for splicing machinery, and demonstrate that these three-dimensional motor-dependent interactions are required to achieve enhanced transcription of specific estrogen-receptor target genes. These findings reveal roles for the modulation of nuclear architecture in orchestrating regulated gene-expression programs in the mammalian nucleus.

LIN28 binds messenger RNAs at GGAGA motifs and regulates splicing factor abundance.

LIN28 is a conserved RNA-binding protein implicated in pluripotency, reprogramming, and oncogenesis. It was previously shown to act primarily by blocking let-7 microRNA (miRNA) biogenesis, but here we elucidate distinct roles of LIN28 regulation via its direct messenger RNA (mRNA) targets. Through crosslinking and immunoprecipitation coupled with high-throughput sequencing (CLIP-seq) in human embryonic stem cells and somatic cells expressing exogenous LIN28, we have defined discrete LIN28-binding sites in a quarter of human transcripts. These sites revealed that LIN28 binds to GGAGA sequences enriched within loop structures in mRNAs, reminiscent of its interaction with let-7 miRNA precursors. Among LIN28 mRNA targets, we found evidence for LIN28 autoregulation and also direct but differing effects on the protein abundance of splicing regulators in somatic and pluripotent stem cells. Splicing-sensitive microarrays demonstrated that exogenous LIN28 expression causes widespread downstream alternative splicing changes. These findings identify important regulatory functions of LIN28 via direct mRNA interactions.

Stability and uniqueness of clonal immunoglobulin CDR3 sequences for MRD tracking in multiple myeloma.

Minimal residual disease (MRD) tracking, by next generation sequencing of immunoglobulin sequences, is moving towards clinical implementation in multiple myeloma. However, there is only sparse information available to address whether clonal sequences remain stable for tracking over time, and to what extent light chain sequences are sufficiently unique for tracking. Here, we analyzed immunoglobulin repertoires from 905 plasma cell myeloma and healthy control samples, focusing on the third complementarity determining region (CDR3). Clonal heavy and/or light chain expression was identified in all patients at baseline, with one or more subclones related to the main clone in 3.2%. In 45 patients with 101 sequential samples, the dominant clonal CDR3 sequences remained identical over time, despite differential clonal evolution by whole exome sequencing in 49% of patients. The low frequency of subclonal CDR3 variants, and absence of evolution over time in active multiple myeloma, indicates that tumor cells at this stage are not under selective pressure to undergo antibody affinity maturation. Next, we establish somatic hypermutation and non-templated insertions as the most important determinants of light chain clonal uniqueness, identifying a potentially trackable sequence in the majority of patients. Taken together, we show that dominant clonal sequences identified at baseline are reliable biomarkers for long-term tracking of the malignant clone, including both IGH and the majority of light chain clones.

Transcriptome-pathology correlation identifies interplay between TDP-43 and the expression of its kinase CK1E in sporadic ALS.

Sporadic amyotrophic lateral sclerosis (sALS) is the most common form of ALS, however, the molecular mechanisms underlying cellular damage and motor neuron degeneration remain elusive. To identify molecular signatures of sALS we performed genome-wide expression profiling in laser capture microdissection-enriched surviving motor neurons (MNs) from lumbar spinal cords of sALS patients with rostral onset and caudal progression. After correcting for immunological background, we discover a highly specific gene expression signature for sALS that is associated with phosphorylated TDP-43 (pTDP-43) pathology. Transcriptome-pathology correlation identified casein kinase 1ε (CSNK1E) mRNA as tightly correlated to levels of pTDP-43 in sALS patients. Enhanced crosslinking and immunoprecipitation in human sALS patient- and healthy control-derived frontal cortex, revealed that TDP-43 binds directly to and regulates the expression of CSNK1E mRNA. Additionally, we were able to show that pTDP-43 itself binds RNA. CK1E, the protein product of CSNK1E, in turn interacts with TDP-43 and promotes cytoplasmic accumulation of pTDP-43 in human stem-cell-derived MNs. Pathological TDP-43 phosphorylation is therefore, reciprocally regulated by CK1E activity and TDP-43 RNA binding. Our framework of transcriptome-pathology correlations identifies candidate genes with relevance to novel mechanisms of neurodegeneration.

Nxf1 natural variant E610G is a semi-dominant suppressor of IAP-induced RNA processing defects.

Endogenous retroviruses and retrotransposons contribute functional genetic variation in animal genomes. In mice, Intracisternal A Particles (IAPs) are a frequent source of both new mutations and polymorphism across laboratory strains. Intronic IAPs can induce alternative RNA processing choices, including alternative splicing. We previously showed IAP I∆1 subfamily insertional mutations are suppressed by a wild-derived allele of the major mRNA export factor, Nxf1. Here we show that a wider diversity of IAP insertions present in the mouse reference sequence induce insertion-dependent alternative processing that is suppressed by Nxf1CAST alleles. These insertions typically show more modest gene expression changes than de novo mutations, suggesting selection or attenuation. Genome-wide splicing-sensitive microarrays and gene-focused assays confirm specificity of Nxf1 genetic modifier activity for IAP insertion alleles. Strikingly, CRISPR/Cas9-mediated genome editing demonstrates that a single amino acid substitution in Nxf1, E610G, is sufficient to recreate a quantitative genetic modifier in a co-isogenic background.

Decoding a signature-based model of transcription cofactor recruitment dictated by cardinal cis-regulatory elements in proximal promoter regions.

Genome-wide maps of DNase I hypersensitive sites (DHSs) reveal that most human promoters contain perpetually active cis-regulatory elements between -150 bp and +50 bp (-150/+50 bp) relative to the transcription start site (TSS). Transcription factors (TFs) recruit cofactors (chromatin remodelers, histone/protein-modifying enzymes, and scaffold proteins) to these elements in order to organize the local chromatin structure and coordinate the balance of post-translational modifications nearby, contributing to the overall regulation of transcription. However, the rules of TF-mediated cofactor recruitment to the -150/+50 bp promoter regions remain poorly understood. Here, we provide evidence for a general model in which a series of cis-regulatory elements (here termed 'cardinal' motifs) prefer acting individually, rather than in fixed combinations, within the -150/+50 bp regions to recruit TFs that dictate cofactor signatures distinctive of specific promoter subsets. Subsequently, human promoters can be subclassified based on the presence of cardinal elements and their associated cofactor signatures. In this study, furthermore, we have focused on promoters containing the nuclear respiratory factor 1 (NRF1) motif as the cardinal cis-regulatory element and have identified the pervasive association of NRF1 with the cofactor lysine-specific demethylase 1 (LSD1/KDM1A). This signature might be distinctive of promoters regulating nuclear-encoded mitochondrial and other particular genes in at least some cells. Together, we propose that decoding a signature-based, expanded model of control at proximal promoter regions should lead to a better understanding of coordinated regulation of gene transcription.

Targeted degradation of sense and antisense C9orf72 RNA foci as therapy for ALS and frontotemporal degeneration.

Expanded hexanucleotide repeats in the chromosome 9 open reading frame 72 (C9orf72) gene are the most common genetic cause of ALS and frontotemporal degeneration (FTD). Here, we identify nuclear RNA foci containing the hexanucleotide expansion (GGGGCC) in patient cells, including white blood cells, fibroblasts, glia, and multiple neuronal cell types (spinal motor, cortical, hippocampal, and cerebellar neurons). RNA foci are not present in sporadic ALS, familial ALS/FTD caused by other mutations (SOD1, TDP-43, or tau), Parkinson disease, or nonneurological controls. Antisense oligonucleotides (ASOs) are identified that reduce GGGGCC-containing nuclear foci without altering overall C9orf72 RNA levels. By contrast, siRNAs fail to reduce nuclear RNA foci despite marked reduction in overall C9orf72 RNAs. Sustained ASO-mediated lowering of C9orf72 RNAs throughout the CNS of mice is demonstrated to be well tolerated, producing no behavioral or pathological features characteristic of ALS/FTD and only limited RNA expression alterations. Genome-wide RNA profiling identifies an RNA signature in fibroblasts from patients with C9orf72 expansion. ASOs targeting sense strand repeat-containing RNAs do not correct this signature, a failure that may be explained, at least in part, by discovery of abundant RNA foci with C9orf72 repeats transcribed in the antisense (GGCCCC) direction, which are not affected by sense strand-targeting ASOs. Taken together, these findings support a therapeutic approach by ASO administration to reduce hexanucleotide repeat-containing RNAs and raise the potential importance of targeting expanded RNAs transcribed in both directions.

Clonal Characterization and Somatic Hypermutation Assessment by Next-Generation Sequencing in Chronic Lymphocytic Leukemia/Small Lymphocytic Lymphoma: A Detailed Description of the Technical Performance, Clinical Utility, and Platform Comparison.

Somatic hypermutation status of the IGHV gene is essential for treating patients with chronic lymphocytic leukemia/small lymphocytic lymphoma. Unlike the conventional low-throughput method, assessment of somatic hypermutation by next-generation sequencing (NGS) has potential for uniformity and scalability. However, it lacks standardization or guidelines for routine clinical use. We critically assessed the performance of an amplicon-based NGS assay across 458 samples. Using a validation cohort (35 samples), the comparison of two platforms (Ion Torrent versus Illumina) and two primer sets [leader versus framework region 1 (FR1)] in their ability to identify clonotypic IGHV rearrangement(s) revealed 97% concordance. The mutation rates were identical by both platforms when using the same primer set (FR1), whereas a slight overestimation bias (+0.326%) was found when comparing FR1 with leader primers. However, for nearly all patients this did not affect the stratification into mutated or unmutated categories, suggesting that use of FR1 may provide comparable results if leader sequencing is not available and allowing for a simpler NGS laboratory workflow. In routine clinical practice (423 samples), the productive rearrangement was successfully detected by either primer set (leader, 97.7%; FR1, 94.7%), and a combination of both in problematic cases reduced the failure rate to 1.2%. Higher sensitivity of the NGS-based analysis also detected a higher frequency of double IGHV rearrangements (19.1%) compared with traditional approaches.

Interlaboratory Analytical Validation of a Next-Generation Sequencing Strategy for Clonotypic Assessment and Minimal Residual Disease Monitoring in Multiple Myeloma.

CONTEXT.­: Minimal residual disease (MRD) is a major prognostic factor in multiple myeloma, although validated technologies are limited. OBJECTIVE.­: To standardize the performance of the LymphoTrack next-generation sequencing (NGS) assays (Invivoscribe), targeting clonal immunoglobulin rearrangements, in order to reproduce the detection of tumor clonotypes and MRD quantitation in myeloma. DESIGN.­: The quantification ability of the assay was evaluated through serial dilution experiments. Paired samples from 101 patients were tested by LymphoTrack, using Sanger sequencing and EuroFlow's next-generation flow (NGF) assay as validated references for diagnostic and follow-up evaluation, respectively. MRD studies using LymphoTrack were performed in parallel at 2 laboratories to evaluate reproducibility. RESULTS.­: Sensitivity was set as 1.3 tumor cells per total number of input cells. Clonality was confirmed in 99% and 100% of cases with Sanger and NGS, respectively, showing great concordance (97.9%), although several samples had minor discordances in the nucleotide sequence of rearrangements. Parallel NGS was performed in 82 follow-up cases, achieving a median sensitivity of 0.001%, while for NGF, median sensitivity was 0.0002%. Reproducibility of LymphoTrack-based MRD studies (85.4%) and correlation with NGF (R2 > 0.800) were high. Bland-Altman tests showed highly significant levels of agreement between flow and sequencing. CONCLUSIONS.­: Taken together, we have shown that LymphoTrack is a suitable strategy for clonality detection and MRD evaluation, with results comparable to gold standard procedures.

Enhancing nuclear receptor-induced transcription requires nuclear motor and LSD1-dependent gene networking in interchromatin granules.

Although the role of liganded nuclear receptors in mediating coactivator/corepressor exchange is well-established, little is known about the potential regulation of chromosomal organization in the 3-dimensional space of the nucleus in achieving integrated transcriptional responses to diverse signaling events. Here, we report that ligand induces rapid interchromosomal interactions among specific subsets of estrogen receptor alpha-bound transcription units, with a dramatic reorganization of nuclear territories, which depends on the actions of nuclear actin/myosin-I machinery and dynein light chain 1. The histone lysine demethylase, LSD1, is required for these ligand-induced interactive loci to associate with distinct interchromatin granules, long thought to serve as "storage" sites for the splicing machinery, some critical transcription elongation factors, and various chromatin remodeling complexes. We demonstrate that this 2-step nuclear rearrangement is essential for achieving enhanced, coordinated transcription of nuclear receptor target genes.

Immunoglobulin gene rearrangement in Koreans with multiple myeloma: Clonality assessment and repertoire analysis using next-generation sequencing.

INTRODUCTION: We assessed the applicability of next-generation sequencing (NGS)-based IGH/IGK clonality testing and analyzed the repertoire of immunoglobulin heavy chain (IGH) or immunoglobulin kappa light chain (IGK) gene usage in Korean patients with multiple myeloma (MM) for the first time. METHODS: Fifty-nine bone marrow samples from 57 Korean patients with MM were analyzed, and NGS-based clonality testing that targeted the IGH and IGK genes was performed using IGH FR1 and IGK primer sets. RESULTS: Clonal IGH and IGK rearrangements were observed in 74.2% and 67.7% of samples from Korean patients with kappa-restricted MM, respectively (90.3% had one or both), and in 60.7% and 95.5% of samples from those with lambda-restricted MM, respectively (85.7% had one or both). In total, 88.1% of samples from Koreans with MM had clonal IGH and/or IGK rearrangement. Clonal rearrangement was not significantly associated with the bone marrow plasma cells as a proportion of all BM lymphoid cells. IGHV3-9 (11.63%) and IGHV4-31 (9.30%) were the most frequently reported IGHV genes and were more common in Koreans with MM than in Western counterparts. IGHD3-10 and IGHD3-3 (13.95% each) were the most frequent IGHD genes; IGHD3-3 was more common in Koreans with MM. No IGK rearrangement was particularly prevalent, but single IGKV-J rearrangements were less common in Koreans with kappa-restricted MM than in Western counterparts. IGKV4-1 was less frequent in Koreans regardless of light chain type. Otherwise, the usages of the IGH V, D, and J genes and of the IGK gene were like those observed in previous Western studies. CONCLUSION: NGS-based IGH/IGK clonality testing ought to be applicable to most Koreans with MM. The overrepresentation of IGHV3-9, IGHV4-31, and IGHD3-3 along with the underrepresentation of IGKV4-1 and the differences in IGK gene rearrangement types suggest the existence of ethnicity-specific variations in this disease.

Routine Evaluation of Minimal Residual Disease in Myeloma Using Next-Generation Sequencing Clonality Testing: Feasibility, Challenges, and Direct Comparison with High-Sensitivity Flow Cytometry.

The 2016 International Myeloma Working Group consensus recommendations emphasize high-sensitivity methods for minimal residual disease (MRD) detection, treatment response assessment, and prognostication. Next-generation sequencing (NGS) of IGH gene rearrangements is highly specific and sensitive, but its description in routine clinical practice and performance comparison with high-sensitivity flow cytometry (hsFC) remain limited. In this large, single-institution study including 438 samples from 251 patients, the use of NGS targeting the IGH and IGK genes for clonal characterization and monitoring, with comparison to hsFC, is described. The index clone characterization success rate was 93.6% (235/251), which depended on plasma cell (PC) cellularity, reaching 98% when PC ≥10% and below 80% when PC <5%. A total of 85% of cases were successfully characterized using leader and FR1 primer sets, and most clones showed high somatic hypermutation rates (median, 8.1%). Among monitoring samples from 124 patients, 78.6% (147/187) had detectable disease by NGS. Concordance with hsFC was 92.9% (170/183). Discordant cases encompassed 8 of 124 hsFC MRD+/NGS MRD- patients (6.5%) and 4 of 124 hsFC MRD-/NGS MRD+ patients (3.2%), all with low-level disease near detection limits for both assays. Among concordant hsFC MRD-/NGS MRD- cases, only 5 of 24 patients (20.8%) showed subsequent overt relapse at 3-year follow-up. HsFC and NGS showed similar operational sensitivity, and the choice of test may depend on practical, rather than test performance, considerations.

Distinct and shared functions of ALS-associated proteins TDP-43, FUS and TAF15 revealed by multisystem analyses.

The RNA-binding protein (RBP) TAF15 is implicated in amyotrophic lateral sclerosis (ALS). To compare TAF15 function to that of two ALS-associated RBPs, FUS and TDP-43, we integrate CLIP-seq and RNA Bind-N-Seq technologies, and show that TAF15 binds to ∼4,900 RNAs enriched for GGUA motifs in adult mouse brains. TAF15 and FUS exhibit similar binding patterns in introns, are enriched in 3' untranslated regions and alter genes distinct from TDP-43. However, unlike FUS and TDP-43, TAF15 has a minimal role in alternative splicing. In human neural progenitors, TAF15 and FUS affect turnover of their RNA targets. In human stem cell-derived motor neurons, the RNA profile associated with concomitant loss of both TAF15 and FUS resembles that observed in the presence of the ALS-associated mutation FUS R521G, but contrasts with late-stage sporadic ALS patients. Taken together, our findings reveal convergent and divergent roles for FUS, TAF15 and TDP-43 in RNA metabolism.

Misregulated RNA processing in amyotrophic lateral sclerosis.

Amyotrophic lateral sclerosis (ALS) research is undergoing an era of unprecedented discoveries with the identification of new genes as major genetic causes of this disease. These discoveries reinforce the genetic, clinical and pathological overlap between ALS and frontotemporal lobar degeneration (FTLD). Common causes of these diseases include mutations in the RNA/DNA-binding proteins, TDP-43 and FUS/TLS and most recently, hexanucleotide expansions in the C9orf72 gene, discoveries that highlight the overlapping pathogenic mechanisms that trigger ALS and FTLD. TDP-43 and FUS/TLS, both of which participate in several steps of RNA processing, are abnormally aggregated and mislocalized in ALS and FTLD, while the expansion in the C9orf72 pre-mRNA strongly suggests sequestration of one or more RNA binding proteins in pathologic RNA foci. Hence, ALS and FTLD converge in pathogenic pathways disrupting the regulation of RNA processing. This article is part of a Special Issue entitled RNA-Binding Proteins.

Divergent roles of ALS-linked proteins FUS/TLS and TDP-43 intersect in processing long pre-mRNAs.

FUS/TLS (fused in sarcoma/translocated in liposarcoma) and TDP-43 are integrally involved in amyotrophic lateral sclerosis (ALS) and frontotemporal dementia. We found that FUS/TLS binds to RNAs from >5,500 genes in mouse and human brain, primarily through a GUGGU-binding motif. We identified a sawtooth-like binding pattern, consistent with co-transcriptional deposition of FUS/TLS. Depletion of FUS/TLS from the adult nervous system altered the levels or splicing of >950 mRNAs, most of which are distinct from RNAs dependent on TDP-43. Abundance of only 45 RNAs was reduced after depletion of either TDP-43 or FUS/TLS from mouse brain, but among these were mRNAs that were transcribed from genes with exceptionally long introns and that encode proteins that are essential for neuronal integrity. Expression levels of a subset of these were lowered after TDP-43 or FUS/TLS depletion in stem cell-derived human neurons and in TDP-43 aggregate-containing motor neurons in sporadic ALS, supporting a common loss-of-function pathway as one component underlying motor neuron death from misregulation of TDP-43 or FUS/TLS.

Long pre-mRNA depletion and RNA missplicing contribute to neuronal vulnerability from loss of TDP-43.

We used cross-linking and immunoprecipitation coupled with high-throughput sequencing to identify binding sites in 6,304 genes as the brain RNA targets for TDP-43, an RNA binding protein that, when mutated, causes amyotrophic lateral sclerosis. Massively parallel sequencing and splicing-sensitive junction arrays revealed that levels of 601 mRNAs were changed (including Fus (Tls), progranulin and other transcripts encoding neurodegenerative disease-associated proteins) and 965 altered splicing events were detected (including in sortilin, the receptor for progranulin) following depletion of TDP-43 from mouse adult brain with antisense oligonucleotides. RNAs whose levels were most depleted by reduction in TDP-43 were derived from genes with very long introns and that encode proteins involved in synaptic activity. Lastly, we found that TDP-43 autoregulates its synthesis, in part by directly binding and enhancing splicing of an intron in the 3' untranslated region of its own transcript, thereby triggering nonsense-mediated RNA degradation.

Comprehensive discovery of endogenous Argonaute binding sites in Caenorhabditis elegans.

MicroRNAs (miRNAs) regulate gene expression by guiding Argonaute proteins to specific target mRNA sequences. Identification of bona fide miRNA target sites in animals is challenging because of uncertainties regarding the base-pairing requirements between miRNA and target as well as the location of functional binding sites within mRNAs. Here we present the results of a comprehensive strategy aimed at isolating endogenous mRNA target sequences bound by the Argonaute protein ALG-1 in C. elegans. Using cross-linking and ALG-1 immunoprecipitation coupled with high-throughput sequencing (CLIP-seq), we identified extensive ALG-1 interactions with specific 3' untranslated region (UTR) and coding exon sequences and discovered features that distinguish miRNA complex binding sites in 3' UTRs from those in other genic regions. Furthermore, our analyses revealed a striking enrichment of Argonaute binding sites in genes important for miRNA function, suggesting an autoregulatory role that may confer robustness to the miRNA pathway.