The identification of susceptibility genes for human disease is a major goal of current biomedical research. Both sequence and structural variation have emerged as major genetic sources of phenotypic variability and growing evidence points to copy number variation as a particularly important source of susceptibility for disease. Here we propose and validate a strategy to identify genes in which changes in dosage alter susceptibility to disease-relevant phenotypes in the mouse. Our approach relies on sensitized phenotypic screening of megabase-sized chromosomal deletion and deficiency lines carrying altered copy numbers of ∼30 linked genes. This approach offers several advantages as a method to systematically identify genes involved in disease susceptibility. To examine the feasibility of such a screen, we performed sensitized phenotyping in five therapeutic areas (metabolic syndrome, immune dysfunction, atherosclerosis, cancer and behaviour) of a 0.8 Mb reciprocal chromosomal duplication and deficiency on chromosome 11 containing 27 genes. Gene dosage in the region significantly affected risk for high-fat diet-induced metabolic syndrome, antigen-induced immune hypersensitivity, ApoE-induced atherosclerosis, and home cage activity. Follow up studies on individual gene knockouts for two candidates in the region showed that copy number variation in Stat5 was responsible for the phenotypic variation in antigen-induced immune hypersensitivity and metabolic syndrome. These data demonstrate the power of sensitized phenotypic screening of segmental aneuploidy lines to identify disease susceptibility genes.

Dosage compensation, mediated by the MSL complex, regulates X-chromosomal gene expression in Drosophila. Here we report that the histone H4 lysine 16 (H4K16) specific histone acetyltransferase MOF displays differential binding behavior depending on whether the target gene is located on the X chromosome versus the autosomes. More specifically, on the male X chromosome, where MSL1 and MSL3 are preferentially associated with the 3' end of dosage compensated genes, MOF displays a bimodal distribution binding to promoters and the 3' ends of genes. In contrast, on MSL1/MSL3 independent X-linked genes and autosomal genes in males and females, MOF binds primarily to promoters. Binding of MOF to autosomes is functional, as H4K16 acetylation and the transcription levels of a number of genes are affected upon MOF depletion. Therefore, MOF is not only involved in the onset of dosage compensation, but also acts as a regulator of gene expression in the Drosophila genome.

Coilin protein scaffolds Cajal bodies (CBs)-subnuclear compartments enriched in small nuclear RNAs (snRNAs)-and promotes efficient spliceosomal snRNP assembly. The molecular function of coilin, which is intrinsically disordered with no defined motifs, is poorly understood. We use UV crosslinking and immunoprecipitation (iCLIP) to determine whether mammalian coilin binds RNA in vivo and to identify targets. Robust detection of snRNA transcripts correlated with coilin ChIP-seq peaks on snRNA genes, indicating that coilin binding to nascent snRNAs is a site-specific CB nucleator. Surprisingly, several hundred small nucleolar RNAs (snoRNAs) were identified as coilin interactors, including numerous unannotated mouse and human snoRNAs. We show that all classes of snoRNAs concentrate in CBs. Moreover, snoRNAs lacking specific CB retention signals traffic through CBs en route to nucleoli, consistent with the role of CBs in small RNP assembly. Thus, coilin couples snRNA and snoRNA biogenesis, making CBs the cellular hub of small ncRNA metabolism.

AAUAAA is the most highly conserved motif in eukaryotic mRNA polyadenylation sites and, in mammals, is specifically recognized by the multisubunit CPSF (cleavage and polyadenylation specificity factor) complex. Despite its critical functions in mRNA 3' end formation, the molecular basis for CPSF-AAUAAA interaction remains poorly defined. The CPSF subunit CPSF160 has been implicated in AAUAAA recognition, but direct evidence has been lacking. Using in vitro and in vivo assays, we unexpectedly found that CPSF subunits CPSF30 and Wdr33 directly contact AAUAAA. Importantly, the CPSF30-RNA interaction is essential for mRNA 3' processing and is primarily mediated by its zinc fingers 2 and 3, which are specifically targeted by the influenza protein NS1A to suppress host mRNA 3' processing. Our data suggest that AAUAAA recognition in mammalian mRNA 3' processing is more complex than previously thought and involves multiple protein-RNA interactions.

RNA-binding proteins (RBPs) regulate splicing according to position-dependent principles, which can be exploited for analysis of regulatory motifs. Here we present RNAmotifs, a method that evaluates the sequence around differentially regulated alternative exons to identify clusters of short and degenerate sequences, referred to as multivalent RNA motifs. We show that diverse RBPs share basic positional principles, but differ in their propensity to enhance or repress exon inclusion. We assess exons differentially spliced between brain and heart, identifying known and new regulatory motifs, and predict the expression pattern of RBPs that bind these motifs. RNAmotifs is available at https://bitbucket.org/rogrro/rna_motifs.

Heterogeneous nuclear ribonucleoprotein C1/C2 (hnRNP C) is a core component of 40S ribonucleoprotein particles that bind pre-mRNAs and influence their processing, stability and export. Breast cancer tumor suppressors BRCA1, BRCA2 and PALB2 form a complex and play key roles in homologous recombination (HR), DNA double strand break (DSB) repair and cell cycle regulation following DNA damage.

The evolution of pediatric solid tumors is poorly understood. There is conflicting evidence of intra-tumor genetic homogeneity vs. heterogeneity (ITGH) in a small number of studies in pediatric solid tumors. A number of copy number aberrations (CNA) are proposed as prognostic biomarkers to stratify patients, for example 1q+ in Wilms tumor (WT); current clinical trials use only one sample per tumor to profile this genetic biomarker. We multisampled 20 WT cases and assessed genome-wide allele-specific CNA and loss of heterozygosity, and inferred tumor evolution, using Illumina CytoSNP12v2.1 arrays, a custom analysis pipeline, and the MEDICC algorithm. We found remarkable diversity of ITGH and evolutionary trajectories in WT. 1q+ is heterogeneous in the majority of tumors with this change, with variable evolutionary timing. We estimate that at least three samples per tumor are needed to detect >95% of cases with 1q+. In contrast, somatic 11p15 LOH is uniformly an early event in WT development. We find evidence of two separate tumor origins in unilateral disease with divergent histology, and in bilateral WT. We also show subclonal changes related to differential response to chemotherapy. Rational trial design to include biomarkers in risk stratification requires tumor multisampling and reliable delineation of ITGH and tumor evolution.

Although the proteins that read the gene regulatory code, transcription factors (TFs), have been largely identified, it is not well known which sequences TFs can recognize. We have analyzed the sequence-specific binding of human TFs using high-throughput SELEX and ChIP sequencing. A total of 830 binding profiles were obtained, describing 239 distinctly different binding specificities. The models represent the majority of human TFs, approximately doubling the coverage compared to existing systematic studies. Our results reveal additional specificity determinants for a large number of factors for which a partial specificity was known, including a commonly observed A- or T-rich stretch that flanks the core motifs. Global analysis of the data revealed that homodimer orientation and spacing preferences, and base-stacking interactions, have a larger role in TF-DNA binding than previously appreciated. We further describe a binding model incorporating these features that is required to understand binding of TFs to DNA.

IHF and HU are two heterodimeric nucleoid-associated proteins (NAP) that belong to the same protein family but interact differently with the DNA. IHF is a sequence-specific DNA-binding protein that bends the DNA by over 160°. HU is the most conserved NAP, which binds non-specifically to duplex DNA with a particular preference for targeting nicked and bent DNA. Despite their importance, the in vivo interactions of the two proteins to the DNA remain to be described at a high resolution and on a genome-wide scale. Further, the effects of these proteins on gene expression on a global scale remain contentious. Finally, the contrast between the functions of the homo- and heterodimeric forms of proteins deserves the attention of further study. Here we present a genome-scale study of HU- and IHF binding to the Escherichia coli K12 chromosome using ChIP-seq. We also perform microarray analysis of gene expression in single- and double-deletion mutants of each protein to identify their regulons. The sequence-specific binding profile of IHF encompasses ∼30% of all operons, though the expression of <10% of these is affected by its deletion suggesting combinatorial control or a molecular backup. The binding profile for HU is reflective of relatively non-specific binding to the chromosome, however, with a preference for A/T-rich DNA. The HU regulon comprises highly conserved genes including those that are essential and possibly supercoiling sensitive. Finally, by performing ChIP-seq experiments, where possible, of each subunit of IHF and HU in the absence of the other subunit, we define genome-wide maps of DNA binding of the proteins in their hetero- and homodimeric forms.

The SR proteins comprise a family of essential, structurally related RNA binding proteins. The complexity of their RNA targets and specificity of RNA recognition in vivo is not well understood. Here we use iCLIP to globally analyze and compare the RNA binding properties of two SR proteins, SRSF3 and SRSF4, in murine cells.

RNA-binding proteins have emerged as causal agents of complex neurological diseases. Mice deficient for neuronal RNA-binding protein CELF4 have a complex neurological disorder with epilepsy as a prominent feature. Human CELF4 has recently been associated with clinical features similar to those seen in mutant mice. CELF4 is expressed primarily in excitatory neurons, including large pyramidal cells of the cerebral cortex and hippocampus, and it regulates excitatory but not inhibitory neurotransmission. We examined mechanisms underlying neuronal hyperexcitability in Celf4 mutants by identifying CELF4 target mRNAs and assessing their fate in the absence of CELF4 in view of their known functions. CELF4 binds to at least 15%-20% of the transcriptome, with striking specificity for the mRNA 3' untranslated region. CELF4 mRNA targets encode a variety of proteins, many of which are well established in neuron development and function. While the overall abundance of these mRNA targets is often dysregulated in Celf4 deficient mice, the actual expression changes are modest at the steady-state level. In contrast, by examining the transcriptome of polysome fractions and the mRNA distribution along the neuronal cell body-neuropil axis, we found that CELF4 is critical for maintaining mRNA stability and availability for translation. Among biological processes associated with CELF4 targets that accumulate in neuropil of mutants, regulation of synaptic plasticity and transmission are the most prominent. Together with a related study of the impact of CELF4 loss on sodium channel Na(v)1.6 function, we suggest that CELF4 deficiency leads to abnormal neuronal function by combining a specific effect on neuronal excitation with a general impairment of synaptic transmission. These results also expand our understanding of the vital roles RNA-binding proteins play in regulating and shaping the activity of neural circuits.

Transcriptomic studies routinely measure expression levels across numerous conditions. These datasets allow identification of genes that are specifically expressed in a small number of conditions. However, there are currently no statistically robust methods for identifying such genes. Here we present SpeCond, a method to detect condition-specific genes that outperforms alternative approaches. We apply the method to a dataset of 32 human tissues to determine 2,673 specifically expressed genes. An implementation of SpeCond is freely available as a Bioconductor package at http://www.bioconductor.org/packages/release/bioc/html/SpeCond.html.

A report on the European Science Foundation workshop 'Transcription Networks: A Global View', Madrid, Spain, 26-28 May 2005.

We established a modified iCLIP protocol, called 'read-through marking', which facilitates the detection of cDNAs that have not been truncated upon encountering the RNA-peptide complex during reverse transcription (read-through cDNAs). A large proportion of these cDNAs would be undesirable in an iCLIP library, as it could affect the resolution of the method. To this end, we added an oligonucleotide to the 5'-end of RNA fragments-a 5'-marker-to mark the read-through cDNAs. By applying this modified iCLIP protocol to PTBP1 and eIF4A3, we found that the start sites of read-through cDNAs are enriched in adenosines, while the remaining cDNAs have a markedly different sequence content at their starts, preferentially containing thymidines. This finding in turn indicates that most of the reads in our iCLIP libraries are a product of truncation with valuable information regarding the proteins' RNA-binding sites. Thus, cDNA start sites confidently identify a protein's RNA-crosslink sites and we can account for the impact of read-through cDNAs by commonly adding a 5'-marker.

Recursive splicing (RS) starts by defining an "RS-exon," which is then spliced to the preceding exon, thus creating a recursive 5' splice site (RS-5ss). Previous studies focused on cryptic RS-exons, and now we find that the exon junction complex (EJC) represses RS of hundreds of annotated, mainly constitutive RS-exons. The core EJC factors, and the peripheral factors PNN and RNPS1, maintain RS-exon inclusion by repressing spliceosomal assembly on RS-5ss. The EJC also blocks 5ss located near exon-exon junctions, thus repressing inclusion of cryptic microexons. The prevalence of annotated RS-exons is high in deuterostomes, while the cryptic RS-exons are more prevalent in Drosophila, where EJC appears less capable of repressing RS. Notably, incomplete repression of RS also contributes to physiological alternative splicing of several human RS-exons. Finally, haploinsufficiency of the EJC factor Magoh in mice is associated with skipping of RS-exons in the brain, with relevance to the microcephaly phenotype and human diseases.

The CRISPR-Cas9 system has successfully been adapted to edit the genome of various organisms. However, our ability to predict the editing outcome at specific sites is limited. Here, we examined indel profiles at over 1,000 genomic sites in human cells and uncovered general principles guiding CRISPR-mediated DNA editing. We find that precision of DNA editing (i.e., recurrence of a specific indel) varies considerably among sites, with some targets showing one highly preferred indel and others displaying numerous infrequent indels. Editing precision correlates with editing efficiency and a preference for single-nucleotide homologous insertions. Precise targets and editing outcome can be predicted based on simple rules that mainly depend on the fourth nucleotide upstream of the protospacer adjacent motif (PAM). Indel profiles are robust, but they can be influenced by chromatin features. Our findings have important implications for clinical applications of CRISPR technology and reveal general patterns of broken end joining that can provide insights into DNA repair mechanisms.

Canonical Wnt/beta-catenin signalling regulates self-renewal and lineage selection within the mammalian epidermis. Although the transcriptional response of keratinocytes that receive a Wnt signal is well characterized, little is known about the mechanism by which keratinocytes in proximity to the Wnt-receiving cell are co-opted to undergo a change in cell fate.

Mutations causing amyotrophic lateral sclerosis (ALS) strongly implicate ubiquitously expressed regulators of RNA processing. To understand the molecular impact of ALS-causing mutations on neuronal development and disease, we analysed transcriptomes during in vitro differentiation of motor neurons (MNs) from human control and patient-specific VCP mutant induced-pluripotent stem cells (iPSCs). We identify increased intron retention (IR) as a dominant feature of the splicing programme during early neural differentiation. Importantly, IR occurs prematurely in VCP mutant cultures compared with control counterparts. These aberrant IR events are also seen in independent RNAseq data sets from SOD1- and FUS-mutant MNs. The most significant IR is seen in the SFPQ transcript. The SFPQ protein binds extensively to its retained intron, exhibits lower nuclear abundance in VCP mutant cultures and is lost from nuclei of MNs in mouse models and human sporadic ALS. Collectively, we demonstrate SFPQ IR and nuclear loss as molecular hallmarks of familial and sporadic ALS.

Methyl-5-uridine (m5U) is one the most abundant non-canonical bases present in cellular RNA, and in yeast is found at position U54 of tRNAs where modification is catalysed by the methyltransferase Trm2. Although the mammalian enzymes that catalyse m5U formation are yet to be identified via experimental evidence, based on sequence homology to Trm2, two candidates currently exist, TRMT2A and TRMT2B. Here we developed a genome-wide single-nucleotide resolution mapping method, Fluorouracil-Induced-Catalytic-Crosslinking-Sequencing (FICC-Seq), in order to identify the relevant enzymatic targets. We demonstrate that TRMT2A is responsible for the majority of m5U present in human RNA, and that it commonly targets U54 of cytosolic tRNAs. By comparison to current methods, we show that FICC-Seq is a particularly robust method for accurate and reliable detection of relevant enzymatic target sites. Our associated finding of extensive irreversible TRMT2A-tRNA crosslinking in vivo following 5-Fluorouracil exposure is also intriguing, as it suggests a tangible mechanism for a previously suspected RNA-dependent route of Fluorouracil-mediated cytotoxicity.

Nonsense-mediated mRNA decay (NMD) controls gene expression by eliminating mRNAs with premature or aberrant translation termination. Degradation of NMD substrates is initiated by the central NMD factor UPF1, which recruits the endonuclease SMG6 and the deadenylation-promoting SMG5/7 complex. The extent to which SMG5/7 and SMG6 contribute to the degradation of individual substrates and their regulation by UPF1 remains elusive. Here we map transcriptome-wide sites of SMG6-mediated endocleavage via 3' fragment capture and degradome sequencing. This reveals that endogenous transcripts can have NMD-eliciting features at various positions, including upstream open reading frames (uORFs), premature termination codons (PTCs), and long 3' UTRs. We find that NMD substrates with PTCs undergo constitutive SMG6-dependent endocleavage, rather than SMG7-dependent exonucleolytic decay. In contrast, the turnover of NMD substrates containing uORFs and long 3' UTRs involves both SMG6- and SMG7-dependent endo- and exonucleolytic decay, respectively. This suggests that the extent to which SMG6 and SMG7 degrade NMD substrates is determined by the mRNA architecture.