Massively targeted evaluation of therapeutic CRISPR off-targets in cells.

Methods for sensitive and high-throughput evaluation of CRISPR RNA-guided nucleases (RGNs) off-targets (OTs) are essential for advancing RGN-based gene therapies. Here we report SURRO-seq for simultaneously evaluating thousands of therapeutic RGN OTs in cells. SURRO-seq captures RGN-induced indels in cells by pooled lentiviral OTs libraries and deep sequencing, an approach comparable and complementary to OTs detection by T7 endonuclease 1, GUIDE-seq, and CIRCLE-seq. Application of SURRO-seq to 8150 OTs from 110 therapeutic RGNs identifies significantly detectable indels in 783 OTs, of which 37 OTs are found in cancer genes and 23 OTs are further validated in five human cell lines by targeted amplicon sequencing. Finally, SURRO-seq reveals that thermodynamically stable wobble base pair (rG•dT) and free binding energy strongly affect RGN specificity. Our study emphasizes the necessity of thoroughly evaluating therapeutic RGN OTs to minimize inevitable off-target effects.

Comparative analysis of clonal evolution among patients with right- and left-sided colon and rectal cancer.

Tumor multiregion sequencing reveals intratumor heterogeneity (ITH) and clonal evolution playing a key role in tumor progression and metastases. Large-scale high-depth multiregional sequencing of colorectal cancer, comparative analysis among patients with right-sided colon cancer (RCC), left-sided colon cancer (LCC), and rectal cancer (RC), as well as the study of lymph node metastasis (LN) with extranodal tumor deposits (ENTDs) from evolutionary perspective remain weakly explored. Here, we recruited 68 patients with RCC (18), LCC (20), and RC (30). We performed high-depth whole-exome sequencing of 206 tumor regions including 176 primary tumors, 19 LN, and 11 ENTD samples. Our results showed ITH with a Darwinian pattern of evolution and the evolution pattern of LCC and RC was more complex and divergent than RCC. Genetic and evolutionary evidences found that both LN and ENTD originated from different clones. Moreover, ENTD was a distinct entity from LN and evolved later.

A porcine brain-wide RNA editing landscape.

Adenosine-to-inosine (A-to-I) RNA editing, catalyzed by ADAR enzymes, is an essential post-transcriptional modification. Although hundreds of thousands of RNA editing sites have been reported in mammals, brain-wide analysis of the RNA editing in the mammalian brain remains rare. Here, a genome-wide RNA-editing investigation is performed in 119 samples, representing 30 anatomically defined subregions in the pig brain. We identify a total of 682,037 A-to-I RNA editing sites of which 97% are not identified before. Within the pig brain, cerebellum and olfactory bulb are regions with most edited transcripts. The editing level of sites residing in protein-coding regions are similar across brain regions, whereas region-distinct editing is observed in repetitive sequences. Highly edited conserved recoding events in pig and human brain are found in neurotransmitter receptors, demonstrating the evolutionary importance of RNA editing in neurotransmission functions. Although potential data biases caused by age, sex or health status are not considered, this study provides a rich resource to better understand the evolutionary importance of post-transcriptional RNA editing.

Single-Cell RNA Sequencing Maps Endothelial Metabolic Plasticity in Pathological Angiogenesis.

Endothelial cell (EC) metabolism is an emerging target for anti-angiogenic therapy in tumor angiogenesis and choroidal neovascularization (CNV), but little is known about individual EC metabolic transcriptomes. By single-cell RNA sequencing 28,337 murine choroidal ECs (CECs) and sprouting CNV-ECs, we constructed a taxonomy to characterize their heterogeneity. Comparison with murine lung tumor ECs (TECs) revealed congruent marker gene expression by distinct EC phenotypes across tissues and diseases, suggesting similar angiogenic mechanisms. Trajectory inference predicted that differentiation of venous to angiogenic ECs was accompanied by metabolic transcriptome plasticity. ECs displayed metabolic transcriptome heterogeneity during cell-cycle progression and in quiescence. Hypothesizing that conserved genes are important, we used an integrated analysis, based on congruent transcriptome analysis, CEC-tailored genome-scale metabolic modeling, and gene expression meta-analysis in cross-species datasets, followed by in vitro and in vivo validation, to identify SQLE and ALDH18A1 as previously unknown metabolic angiogenic targets.

An Integrated Gene Expression Landscape Profiling Approach to Identify Lung Tumor Endothelial Cell Heterogeneity and Angiogenic Candidates.

Heterogeneity of lung tumor endothelial cell (TEC) phenotypes across patients, species (human/mouse), and models (in vivo/in vitro) remains poorly inventoried at the single-cell level. We single-cell RNA (scRNA)-sequenced 56,771 endothelial cells from human/mouse (peri)-tumoral lung and cultured human lung TECs, and detected 17 known and 16 previously unrecognized phenotypes, including TECs putatively regulating immune surveillance. We resolved the canonical tip TECs into a known migratory tip and a putative basement-membrane remodeling breach phenotype. Tip TEC signatures correlated with patient survival, and tip/breach TECs were most sensitive to vascular endothelial growth factor blockade. Only tip TECs were congruent across species/models and shared conserved markers. Integrated analysis of the scRNA-sequenced data with orthogonal multi-omics and meta-analysis data across different human tumors, validated by functional analysis, identified collagen modification as a candidate angiogenic pathway.

Haplotyping by CRISPR-mediated DNA circularization (CRISPR-hapC) broadens allele-specific gene editing.

Allele-specific protospacer adjacent motif (asPAM)-positioning SNPs and CRISPRs are valuable resources for gene therapy of dominant disorders. However, one technical hurdle is to identify the haplotype comprising the disease-causing allele and the distal asPAM SNPs. Here, we describe a novel CRISPR-based method (CRISPR-hapC) for haplotyping. Based on the generation (with a pair of CRISPRs) of extrachromosomal circular DNA in cells, the CRISPR-hapC can map haplotypes from a few hundred bases to over 200 Mb. To streamline and demonstrate the applicability of the CRISPR-hapC and asPAM CRISPR for allele-specific gene editing, we reanalyzed the 1000 human pan-genome and generated a high frequency asPAM SNP and CRISPR database (www.crispratlas.com/knockout) for four CRISPR systems (SaCas9, SpCas9, xCas9 and Cas12a). Using the huntingtin (HTT) CAG expansion and transthyretin (TTR) exon 2 mutation as examples, we showed that the asPAM CRISPRs can specifically discriminate active and dead PAMs for all 23 loci tested. Combination of the CRISPR-hapC and asPAM CRISPRs further demonstrated the capability for achieving highly accurate and haplotype-specific deletion of the HTT CAG expansion allele and TTR exon 2 mutation in human cells. Taken together, our study provides a new approach and an important resource for genome research and allele-specific (haplotype-specific) gene therapy.

LION: a simple and rapid method to achieve CRISPR gene editing.

The RNA-guided CRISPR-Cas9 technology has paved the way for rapid and cost-effective gene editing. However, there is still a great need for effective methods for rapid generation and validation of CRISPR/Cas9 gRNAs. Previously, we have demonstrated that highly efficient generation of multiplexed CRISPR guide RNA (gRNA) expression array can be achieved with Golden Gate Assembly (GGA). Here, we present an optimized and rapid method for generation and validation in less than 1 day of CRISPR gene targeting vectors. The method (LION) is based on ligation of double-stranded gRNA oligos into CRISPR vectors with GGA followed by nucleic acid purification. Using a dual-fluorescent reporter vector (C-Check), T7E1 assay, TIDE assay and a traffic light reporter assay, we proved that the LION-based generation of CRISPR vectors are functionally active, and equivalent to CRISPR plasmids generated by traditional methods. We also tested the activity of LION CRISPR vectors in different human cell types. The LION method presented here advances the rapid functional validation and application of CRISPR system for gene editing and simplified the CRISPR gene-editing procedures.

The correlation of copy number variations with longevity in a genome-wide association study of Han Chinese.

Copy number variations (CNVs) have been shown to cause numerous diseases, however, their roles in human lifespan remain elusive. In this study, we investigate the association of CNVs with longevity by comparing the Han Chinese genomes of long-lived individuals from 90 to 117 years of age and the middle-aged from 30 to 65. Our data demonstrate that the numbers of CNVs, especially deletions, increase significantly in a direct correlation with longevity. We identify eleven CNVs that strongly associate with longevity; four of them locate in the chromosome bands, 7p11.2, 20q13.33, 19p12 and 8p23.3 and overlap partially with the CNVs identified in long-lived Danish or U.S. populations, while the other seven have not been reported previously. These CNV regions encode nineteen known genes, and some of which have been shown to affect aging-related phenotypes such as the shortening of telomere length (ZNF208), the risk of cancer (FOXA1, LAMA5, ZNF716), and vascular and immune-related diseases (ARHGEF10, TOR2A, SH2D3C). In addition, we found several pathways enriched in long-lived genomes, including FOXA1 and FOXA transcription factor networks involved in regulating aging or age-dependent diseases such as cancer. Thus, our study has identified longevity-associated CNV regions and their affected genes and pathways. Our results suggest that the human genome structures such as CNVs might play an important role in determining a long life in human.

BS-virus-finder: virus integration calling using bisulfite sequencing data.

Background: DNA methylation plays a key role in the regulation of gene expression and carcinogenesis. Bisulfite sequencing studies mainly focus on calling single nucleotide polymorphism, different methylation region, and find allele-specific DNA methylation. Until now, only a few software tools have focused on virus integration using bisulfite sequencing data. Findings: We have developed a new and easy-to-use software tool, named BS-virus-finder (BSVF, RRID:SCR_015727), to detect viral integration breakpoints in whole human genomes. The tool is hosted at https://github.com/BGI-SZ/BSVF. Conclusions: BS-virus-finder demonstrates high sensitivity and specificity. It is useful in epigenetic studies and to reveal the relationship between viral integration and DNA methylation. BS-virus-finder is the first software tool to detect virus integration loci by using bisulfite sequencing data.

RNA-Guided Activation of Pluripotency Genes in Human Fibroblasts.

Specific activation of endogenous genes can be achieved by programmable artificial transcription factors (ATFs). In this study, we compared two artificial, programmable, clustered regularly interspaced short palindromic repeats (CRISPR)-based, ubiquitous transcription factors: deficient CRISPR-associated protein 9 (dCas9)-VP64 (CRISPRa) alone, or a combination of dCas9-VP64 and MS2-P65-HSF1 [synergistic activation mediator (SAM) system] mediated activation of five pluripotency genes: KLF4 (K), LIN28 (L), MYC (M), OCT4 (O), and SOX2 (S) in human cells (HEK293T, HeLa, HepG2, and primary fibroblasts). Activation potential was monitored using a luciferase reporter system and we found that both CRISPRa and SAM can efficiently activate the proximal promoter of all five genes. We also observed that the guide RNA (gRNA) target sites and number of gRNAs have a major effect on gRNA-guided activation efficiency. Furthermore, increased activation efficiency (>3-folds) could be achieved by the SAM system compared to CRISPRa. In addition, we discovered that only the SAM system could efficiently activate LIN28, OCT4, and SOX2 expression (up to 100-folds compared to coexpression with a scrambled gRNA) in primary human fibroblasts. This SAM-mediated activation of LOS can be stably maintained for over 20 days in fibroblasts cultured in either fibroblasts or stem cell medium. However, when attempting to use the SAM-LOS activation as an approach for induced pluripotent stem cells-reprogramming, no embryonic stem-like colonies could be obtained from these SAM fibroblasts. In conclusion, our study showed that CRISPR/Cas9-based ATFs are potent to activate and maintain transcription of endogenous human pluripotent genes. However, future improvements of the system are still required to improve activation efficiency and cellular reprogramming using ATFs.

Fusion of SpCas9 to E. coli Rec A protein enhances CRISPR-Cas9 mediated gene knockout in mammalian cells.

Mammalian cells repair double-strand DNA breaks (DSB) by a range of different pathways following DSB induction by the engineered clustered regularly interspaced short palindromic repeats (CRISPR)-associated protein Cas9. While CRISPR-Cas9 thus enables predesigned modifications of the genome, applications of CRISPR-Cas9-mediated genome-editing are frequently hampered by the unpredictable and varying pathways for DSB repair in mammalian cells. Here we present a strategy of fusing Cas9 to recombinant proteins for fine-tuning of the DSB repair preferences in mammalian cells. By fusing Streptococcus Pyogenes Cas9 (SpCas9) to the recombinant protein A (Rec A, NP_417179.1) from Escherichia coli, we create a recombinant Cas9 protein (rSpCas9) which enhances the generation of indel mutations at DSB sites in mammalian cells, increases the frequency of DSB repair by homology-directed single-strand annealing (SSA), and represses homology-directed gene conversion by approximately 33%. Our study thus proves for the first time that fusing SpCas9 to recombinant proteins can influence the balance between DSB repair pathways in mammalian cells. This approach may form the basis for further investigations of the applications of recombinant Cas9 proteins to fine-tuning DSB repair pathways in eukaryotic cells.

Systematic in vitro and in vivo characterization of Leukemia-inhibiting factor- and Fibroblast growth factor-derived porcine induced pluripotent stem cells.

Derivation and stable maintenance of porcine induced pluripotent stem cells (piPSCs) is challenging. We herein systematically analyzed two piPSC lines, derived by lentiviral transduction and cultured under either leukemia inhibitory factor (LIF) or fibroblast growth factor (FGF) conditions, to shed more light on the underlying biological mechanisms of porcine pluripotency. LIF-derived piPSCs were more successful than their FGF-derived counterparts in the generation of in vitro chimeras and in teratoma formation. When LIF piPSCs chimeras were transferred into surrogate sows and allowed to develop, only their prescence within the embryonic membranes could be detected. Whole-transcriptome analysis of the piPSCs and porcine neonatal fibroblasts showed that they clustered together, but apart from the two pluripotent cell populations of early porcine embryos, indicating incomplete reprogramming. Indeed, bioinformatic analysis of the pluripotency-related gene network of the LIF- versus FGF-derived piPSCs revealed that ZFP42 (REX1) expression was absent in both piPSC-like cells, whereas it was expressed in the porcine inner cell mass at Day 7/8. A second striking difference was the expression of ATOH1 in piPSC-like cells, which was absent in the inner cell mass. Moreover, our gene expression analyses plus correlation analyses of known pluripotency genes identified unique relationships between pluripotency genes in the inner cell mass, which are to some extent, in the piPSC-like cells. This deficiency in downstream gene activation and divergent gene expression may be underlie the inability to derive germ line-transmitting piPSCs, and provides unique insight into which genes are necessary to achieve fully reprogrammed piPSCs. 84: 229-245, 2017. © 2016 Wiley Periodicals, Inc.

Apolipoprotein E Deficiency Increases Remnant Lipoproteins and Accelerates Progressive Atherosclerosis, But Not Xanthoma Formation, in Gene-Modified Minipigs.

Deficiency of apolipoprotein E (APOE) causes familial dysbetalipoproteinemia in humans resulting in a higher risk of atherosclerotic disease. In mice, APOE deficiency results in a severe atherosclerosis phenotype, but it is unknown to what extent this is unique to mice. In this study, APOE was targeted in Yucatan minipigs. APOE-/- minipigs displayed increased plasma cholesterol and accumulation of apolipoprotein B-48-containing chylomicron remnants on low-fat diet, which was significantly accentuated upon feeding a high-fat, high-cholesterol diet. APOE-/- minipigs displayed accelerated progressive atherosclerosis but not xanthoma formation. This indicates that remnant lipoproteinemia does not induce early lesions but is atherogenic in pre-existing atherosclerosis.

Production of Pigs by Hand-Made Cloning Using Mesenchymal Stem Cells and Fibroblasts.

Mesenchymal stem cells (MSCs) exhibited self-renewal and less differentiation, making the MSCs promising candidates for adult somatic cell nuclear transfer (SCNT). In this article, we tried to produce genome identical pigs through hand-made cloning (HMC), with MSCs and adult skin fibroblasts as donor cells. MSCs were derived from either adipose tissue or peripheral blood (aMSCs and bMSCs, respectively). MSCs usually showed the expression pattern of CD29, CD73, CD90, and CD105 together with lack of expression of the hematopoietic markers CD34and CD45. Flow cytometry results demonstrated high expression of CD29 and CD90 in both MSC lines, while CD73, CD34, and CD45 expression were not detected. In contrary, in reverse transcription-polymerase chain reaction (RT-PCR) analysis, CD73 and CD34 were detected indicating that human antibodies CD73 and CD34 were not suitable to identify porcine cell surface markers and porcine MSC cellular surface markers of CD34 might be different from other species. MSCs also had potential to differentiate successfully into chondrocytes, osteoblasts, and adipocytes. After HMC, embryos reconstructed with aMSCs had higher blastocyst rate on day 5 and 6 than those reconstructed with bMSCs and fibroblasts (29.6% ± 1.3% and 41.1% ± 1.4% for aMSCs vs. 23.9% ± 1.2% and 35.5% ± 1.6% for bMSCs and 22.1% ± 0.9% and 33.3% ± 1.1% for fibroblasts, respectively). Live birth rate per transferred blastocyst achieved with bMSCs (1.59%) was the highest among the three groups. This article was the first report to compare the efficiency among bMSCs, aMSCs, and fibroblasts for boar cloning, which offered a realistic perspective to use the HMC technology for commercial breeding.

The stepwise evolution of the exome during acquisition of docetaxel resistance in breast cancer cells.

BACKGROUND: Resistance to taxane-based therapy in breast cancer patients is a major clinical problem that may be addressed through insight of the genomic alterations leading to taxane resistance in breast cancer cells. In the current study we used whole exome sequencing to discover somatic genomic alterations, evolving across evolutionary stages during the acquisition of docetaxel resistance in breast cancer cell lines. RESULTS: Two human breast cancer in vitro models (MCF-7 and MDA-MB-231) of the step-wise acquisition of docetaxel resistance were developed by exposing cells to 18 gradually increasing concentrations of docetaxel. Whole exome sequencing performed at five successive stages during this process was used to identify single point mutational events, insertions/deletions and copy number alterations associated with the acquisition of docetaxel resistance. Acquired coding variation undergoing positive selection and harboring characteristics likely to be functional were further prioritized using network-based approaches. A number of genomic changes were found to be undergoing evolutionary selection, some of which were likely to be functional. Of the five stages of progression toward resistance, most resistance relevant genomic variation appeared to arise midway towards fully resistant cells corresponding to passage 31 (5 nM docetaxel) for MDA-MB-231 and passage 16 (1.2 nM docetaxel) for MCF-7, and where the cells also exhibited a period of reduced growth rate or arrest, respectively. MCF-7 cell acquired several copy number gains on chromosome 7, including ABC transporter genes, including ABCB1 and ABCB4, as well as DMTF1, CLDN12, CROT, and SRI. For MDA-MB-231 numerous copy number losses on chromosome X involving more than 30 genes was observed. Of these genes, CASK, POLA1, PRDX4, MED14 and PIGA were highly prioritized by the applied network-based gene ranking approach. At higher docetaxel concentration MCF-7 subclones exhibited a copy number loss in E2F4, and the gene encoding this important transcription factor was down-regulated in MCF-7 resistant cells. CONCLUSIONS: Our study of the evolution of acquired docetaxel resistance identified several genomic changes that might explain development of docetaxel resistance. Interestingly, the most relevant resistance-associated changes appeared to originate midway through the evolution towards fully resistant cell lines. Our data suggest that no single genomic event sufficiently predicts resistance to docetaxel, but require genomic alterations affecting multiple pathways that in concert establish the final resistance stage.

Enhanced genome editing in mammalian cells with a modified dual-fluorescent surrogate system.

Programmable DNA nucleases such as TALENs and CRISPR/Cas9 are emerging as powerful tools for genome editing. Dual-fluorescent surrogate systems have been demonstrated by several studies to recapitulate DNA nuclease activity and enrich for genetically edited cells. In this study, we created a single-strand annealing-directed, dual-fluorescent surrogate reporter system, referred to as C-Check. We opted for the Golden Gate Cloning strategy to simplify C-Check construction. To demonstrate the utility of the C-Check system, we used the C-Check in combination with TALENs or CRISPR/Cas9 in different scenarios of gene editing experiments. First, we disrupted the endogenous pIAPP gene (3.0 % efficiency) by C-Check-validated TALENs in primary porcine fibroblasts (PPFs). Next, we achieved gene-editing efficiencies of 9.0-20.3 and 4.9 % when performing single- and double-gene targeting (MAPT and SORL1), respectively, in PPFs using C-Check-validated CRISPR/Cas9 vectors. Third, fluorescent tagging of endogenous genes (MYH6 and COL2A1, up to 10.0 % frequency) was achieved in human fibroblasts with C-Check-validated CRISPR/Cas9 vectors. We further demonstrated that the C-Check system could be applied to enrich for IGF1R null HEK293T cells and CBX5 null MCF-7 cells with frequencies of nearly 100.0 and 86.9 %, respectively. Most importantly, we further showed that the C-Check system is compatible with multiplexing and for studying CRISPR/Cas9 sgRNA specificity. The C-Check system may serve as an alternative dual-fluorescent surrogate tool for measuring DNA nuclease activity and enrichment of gene-edited cells, and may thereby aid in streamlining programmable DNA nuclease-mediated genome editing and biological research.

Towards Personalized Regenerative Cell Therapy: Mesenchymal Stem Cells Derived from Human Induced Pluripotent Stem Cells.

Mesenchymal stem cells (MSCs) are adult stem cells with the capacity of self-renewal and multilineage differentiation, and can be isolated from several adult tissues. However, isolating MSCs from adult tissues for cell therapy is hampered by the invasive procedure, the rarity of the cells and their attenuated proliferation capacity when cultivated and expanded in vitro. Human MSCs derived from induced pluripotent stem cells (iPSC-MSCs) have now evolved as a promising alternative cell source for MSCs and regenerative medicine. Several groups, including ours, have reported successful derivation of functional iPSC-MSCs and applied these cells in MSC-based therapeutic testing. Still, the current experience and understanding of iPSC-MSCs with respect to production methods, safety and efficacy are primitive. In this review, we highlight the methodological progress in iPSC-MSC research, describing the importance of choosing the right sources of iPSCs, iPSC reprogramming methods, iPSC culture systems, embryoid body intermediates, pathway inhibitors, basal medium, serum, growth factors and culture surface coating. We also highlight some progress in the application of iPSC-MSCs in direct cell therapy, tissue engineering and gene therapy.

Interaction Between the FOXO1A-209 Genotype and Tea Drinking Is Significantly Associated with Reduced Mortality at Advanced Ages.

On the basis of the genotypic/phenotypic data from Chinese Longitudinal Healthy Longevity Survey (CLHLS) and Cox proportional hazard model, the present study demonstrates that interactions between carrying FOXO1A-209 genotypes and tea drinking are significantly associated with lower risk of mortality at advanced ages. Such a significant association is replicated in two independent Han Chinese CLHLS cohorts (p = 0.028-0.048 in the discovery and replication cohorts, and p = 0.003-0.016 in the combined dataset). We found the associations between tea drinking and reduced mortality are much stronger among carriers of the FOXO1A-209 genotype compared to non-carriers, and drinking tea is associated with a reversal of the negative effects of carrying FOXO1A-209 minor alleles, that is, from a substantially increased mortality risk to substantially reduced mortality risk at advanced ages. The impacts are considerably stronger among those who carry two copies of the FOXO1A minor allele than those who carry one copy. On the basis of previously reported experiments on human cell models concerning FOXO1A-by-tea-compounds interactions, we speculate that results in the present study indicate that tea drinking may inhibit FOXO1A-209 gene expression and its biological functions, which reduces the negative impacts of FOXO1A-209 gene on longevity (as reported in the literature) and offers protection against mortality risk at oldest-old ages. Our empirical findings imply that the health outcomes of particular nutritional interventions, including tea drinking, may, in part, depend upon individual genetic profiles, and the research on the effects of nutrigenomics interactions could potentially be useful for rejuvenation therapies in the clinic or associated healthy aging intervention programs.

EWS and FUS bind a subset of transcribed genes encoding proteins enriched in RNA regulatory functions.

BACKGROUND: FUS (TLS) and EWS (EWSR1) belong to the FET-protein family of RNA and DNA binding proteins. FUS and EWS are structurally and functionally related and participate in transcriptional regulation and RNA processing. FUS and EWS are identified in translocation generated cancer fusion proteins and involved in the human neurological diseases amyotrophic lateral sclerosis and fronto-temporal lobar degeneration. RESULTS: To determine the gene regulatory functions of FUS and EWS at the level of chromatin, we have performed chromatin immunoprecipitation followed by next generation sequencing (ChIP-seq). Our results show that FUS and EWS bind to a subset of actively transcribed genes, that binding often is downstream the poly(A)-signal, and that binding overlaps with RNA polymerase II. Functional examinations of selected target genes identified that FUS and EWS can regulate gene expression at different levels. Gene Ontology analyses showed that FUS and EWS target genes preferentially encode proteins involved in regulatory processes at the RNA level. CONCLUSIONS: The presented results yield new insights into gene interactions of EWS and FUS and have identified a set of FUS and EWS target genes involved in pathways at the RNA regulatory level with potential to mediate normal and disease-associated functions of the FUS and EWS proteins.