DNA damage repair proteins across the Tree of Life.

Genome maintenance is orchestrated by a highly regulated DNA damage response with specific DNA repair pathways. Here, we investigate the phylogenetic diversity in the recognition and repair of three well-established DNA lesions, primarily repaired by base excision repair (BER) and ribonucleotide excision repair (RER): (1) 8-oxoguanine, (2) abasic site, and (3) incorporated ribonucleotide in DNA in 11 species: Escherichia coli, Bacillus subtilis, Halobacterium salinarum, Trypanosoma brucei, Tetrahymena thermophila, Saccharomyces cerevisiae, Schizosaccharomyces pombe, Caenorhabditis elegans, Homo sapiens, Arabidopsis thaliana, and Zea mays. Using quantitative mass spectrometry, we identified 337 binding proteins across these species. Of these proteins, 99 were previously characterized to be involved in DNA repair. Through orthology, network, and domain analysis, we linked 44 previously unconnected proteins to DNA repair. Our study presents a resource for future study of the crosstalk and evolutionary conservation of DNA damage repair across all domains of life.

Identification of Pri-miRNA Stem-Loop Interacting Proteins in Plants Using a Modified Version of the Csy4 CRISPR Endonuclease.

Regulation at the RNA level by RNA-binding proteins (RBPs) and microRNAs (miRNAs) is key to coordinating eukaryotic gene expression. In plants, the importance of miRNAs is highlighted by severe developmental defects in mutants impaired in miRNA biogenesis. MiRNAs are processed from long primary-microRNAs (pri-miRNAs) with internal stem-loop structures by endonucleolytic cleavage. The highly structured stem-loops constitute the basis for the extensive regulation of miRNA biogenesis through interaction with RBPs. However, trans-acting regulators of the biogenesis of specific miRNAs are largely unknown in plants. Therefore, we exploit an RNA-centric approach based on modified versions of the conditional CRISPR nuclease Csy4* to pull down interactors of the Arabidopsis pri-miR398b stem-loop (pri-miR398b-SL) in vitro. We designed three epitope-tagged versions of the inactive Csy4* for the immobilization of the protein together with the pri-miR398b-SL bait on high affinity matrices. After incubation with nucleoplasmic extracts from Arabidopsis and extensive washing, pri-miR398b-SL, along with its specifically bound proteins, were released by re-activating the cleavage activity of the Csy4* upon the addition of imidazole. Co-purified proteins were identified via quantitative mass spectrometry and data sets were compared. In total, we identified more than 400 different proteins, of which 180 are co-purified in at least two out of three independent Csy4*-based RNA pulldowns. Among those, the glycine-rich RNA-binding protein AtRZ-1a was identified in all pulldowns. To analyze the role of AtRZ-1a in miRNA biogenesis, we determined the miR398 expression level in the atrz-1a mutant. Indeed, the absence of AtRZ-1a caused a decrease in the steady-state level of mature miR398 with a concomitant reduction in pri-miR398b levels. Overall, we show that our modified Csy4*-based RNA pulldown strategy is suitable to identify new trans-acting regulators of miRNA biogenesis and provides new insights into the post-transcriptional regulation of miRNA processing by plant RBPs.

The nucleolar DExD/H protein Hel66 is involved in ribosome biogenesis in Trypanosoma brucei.

The biosynthesis of ribosomes is a complex cellular process involving ribosomal RNA, ribosomal proteins and several further trans-acting factors. DExD/H box proteins constitute the largest family of trans-acting protein factors involved in this process. Several members of this protein family have been directly implicated in ribosome biogenesis in yeast. In trypanosomes, ribosome biogenesis differs in several features from the process described in yeast. Here, we have identified the DExD/H box helicase Hel66 as being involved in ribosome biogenesis. The protein is unique to Kinetoplastida, localises to the nucleolus and its depletion via RNAi caused a severe growth defect. Loss of the protein resulted in a decrease of global translation and accumulation of rRNA processing intermediates for both the small and large ribosomal subunits. Only a few factors involved in trypanosome rRNA biogenesis have been described so far and our findings contribute to gaining a more comprehensive picture of this essential process.

Ythdf is a N6-methyladenosine reader that modulates Fmr1 target mRNA selection and restricts axonal growth in Drosophila.

N6-methyladenosine (m6 A) regulates a variety of physiological processes through modulation of RNA metabolism. This modification is particularly enriched in the nervous system of several species, and its dysregulation has been associated with neurodevelopmental defects and neural dysfunctions. In Drosophila, loss of m6 A alters fly behavior, albeit the underlying molecular mechanism and the role of m6 A during nervous system development have remained elusive. Here we find that impairment of the m6 A pathway leads to axonal overgrowth and misguidance at larval neuromuscular junctions as well as in the adult mushroom bodies. We identify Ythdf as the main m6 A reader in the nervous system, being required to limit axonal growth. Mechanistically, we show that the m6 A reader Ythdf directly interacts with Fmr1, the fly homolog of Fragile X mental retardation RNA binding protein (FMRP), to inhibit the translation of key transcripts involved in axonal growth regulation. Altogether, this study demonstrates that the m6 A pathway controls development of the nervous system and modulates Fmr1 target transcript selection.

PAR-TERRA is the main contributor to telomeric repeat-containing RNA transcripts in normal and cancer mouse cells.

Telomeric repeat-containing RNA (TERRA) molecules play important roles at telomeres, from heterochromatin regulation to telomerase activity control. In human cells, TERRA is transcribed from subtelomeric promoters located on most chromosome ends and associates with telomeres. The origin of mouse TERRA molecules is, however, unclear, as transcription from the pseudoautosomal PAR locus was recently suggested to account for the vast majority of TERRA in embryonic stem cells (ESC). Here, we confirm the production of TERRA from both the chromosome 18q telomere and the PAR locus in mouse embryonic fibroblasts, ESC, and various mouse cancer and immortalized cell lines, and we identify two novel sources of TERRA on mouse chromosome 2 and X. Using various approaches, we show that PAR-TERRA molecules account for the majority of TERRA transcripts, displaying an increase of two to four orders of magnitude compared to the telomeric 18q transcript. Finally, we present a SILAC-based pull-down screen revealing a large overlap between TERRA-interacting proteins in human and mouse cells, including PRC2 complex subunits, chromatin remodeling factors, DNA replication proteins, Aurora kinases, shelterin complex subunits, Bloom helicase, Coilin, and paraspeckle proteins. Hence, despite originating from distinct genomic regions, mouse and human TERRA are likely to play similar functions in cells.

Correction to: Proteotranscriptomics assisted gene annotation and spatial proteomics of Bombyx mori BmN4 cell line.

An amendment to this paper has been published and can be accessed via the original article.

Proteotranscriptomics assisted gene annotation and spatial proteomics of Bombyx mori BmN4 cell line.

BACKGROUND: The process of identifying all coding regions in a genome is crucial for any study at the level of molecular biology, ranging from single-gene cloning to genome-wide measurements using RNA-seq or mass spectrometry. While satisfactory annotation has been made feasible for well-studied model organisms through great efforts of big consortia, for most systems this kind of data is either absent or not adequately precise. RESULTS: Combining in-depth transcriptome sequencing and high resolution mass spectrometry, we here use proteotranscriptomics to improve gene annotation of protein-coding genes in the Bombyx mori cell line BmN4 which is an increasingly used tool for the analysis of piRNA biogenesis and function. Using this approach we provide the exact coding sequence and evidence for more than 6200 genes on the protein level. Furthermore using spatial proteomics, we establish the subcellular localization of thousands of these proteins. We show that our approach outperforms current Bombyx mori annotation attempts in terms of accuracy and coverage. CONCLUSIONS: We show that proteotranscriptomics is an efficient, cost-effective and accurate approach to improve previous annotations or generate new gene models. As this technique is based on de-novo transcriptome assembly, it provides the possibility to study any species also in the absence of genome sequence information for which proteogenomics would be impossible.

The RNA fold interactome of evolutionary conserved RNA structures in S. cerevisiae.

RNA-binding proteins play key roles in regulation of gene expression via recognition of structural features in RNA molecules. Here we apply a quantitative RNA pull-down approach to 186 evolutionary conserved RNA structures and report 162 interacting proteins. Unlike global RNA interactome capture, we associate individual RNA structures within messenger RNA with their interacting proteins. Of our binders 69% are known RNA-binding proteins, whereas some are previously unrelated to RNA binding and do not harbor canonical RNA-binding domains. While current knowledge about RNA-binding proteins relates to their functions at 5' or 3'-UTRs, we report a significant number of them binding to RNA folds in the coding regions of mRNAs. Using an in vivo reporter screen and pulsed SILAC, we characterize a subset of mRNA-RBP pairs and thus connect structural RNA features to functionality. Ultimately, we here present a generic, scalable approach to interrogate the increasing number of RNA structural motifs.

Quantitative Proteomics to Identify Nuclear RNA-Binding Proteins of Malat1.

The long non-coding RNA Malat1 has been implicated in several human cancers, while the mechanism of action is not completely understood. As RNAs in cells function together with RNA-binding proteins (RBPs), the composition of their RBP complex can shed light on their functionality. We here performed quantitative interactomics of 14 non-overlapping fragments covering the full length of Malat1 to identify possible nuclear interacting proteins. Overall, we identified 35 candidates including 14 already known binders, which are able to interact with Malat1 in the nucleus. Furthermore, the use of fragments along the full-length RNA allowed us to reveal two hotspots for protein binding, one in the 5'-region and one in the 3'-region of Malat1. Our results provide confirmation on previous RNA-protein interaction studies and suggest new candidates for functional investigations.

A DNA methylation reader complex that enhances gene transcription.

DNA methylation generally functions as a repressive transcriptional signal, but it is also known to activate gene expression. In either case, the downstream factors remain largely unknown. By using comparative interactomics, we isolated proteins in Arabidopsis thaliana that associate with methylated DNA. Two SU(VAR)3-9 homologs, the transcriptional antisilencing factor SUVH1, and SUVH3, were among the methyl reader candidates. SUVH1 and SUVH3 bound methylated DNA in vitro, were associated with euchromatic methylation in vivo, and formed a complex with two DNAJ domain-containing homologs, DNAJ1 and DNAJ2. Ectopic recruitment of DNAJ1 enhanced gene transcription in plants, yeast, and mammals. Thus, the SUVH proteins bind to methylated DNA and recruit the DNAJ proteins to enhance proximal gene expression, thereby counteracting the repressive effects of transposon insertion near genes.

Genome-wide Rules of Nucleosome Phasing in Drosophila.

Regular successions of positioned nucleosomes, or phased nucleosome arrays (PNAs), are predominantly known from transcriptional start sites (TSSs). It is unclear whether PNAs occur elsewhere in the genome. To generate a comprehensive inventory of PNAs for Drosophila, we applied spectral analysis to nucleosome maps and identified thousands of PNAs throughout the genome. About half of them are not near TSSs and are strongly enriched for an uncharacterized sequence motif. Through genome-wide reconstitution of physiological chromatin in Drosophila embryo extracts, we uncovered the molecular basis of PNA formation. We identified Phaser, an unstudied zinc finger protein that positions nucleosomes flanking the motif. It also revealed how the global activity of the chromatin remodelers CHRAC/ACF, together with local barrier elements, generates islands of regular phasing throughout the genome. Our work demonstrates the potential of chromatin assembly by embryo extracts as a powerful tool to reconstitute chromatin features on a global scale in vitro.

Phylointeractomics reconstructs functional evolution of protein binding.

Molecular phylogenomics investigates evolutionary relationships based on genomic data. However, despite genomic sequence conservation, changes in protein interactions can occur relatively rapidly and may cause strong functional diversification. To investigate such functional evolution, we here combine phylogenomics with interaction proteomics. We develop this concept by investigating the molecular evolution of the shelterin complex, which protects telomeres, across 16 vertebrate species from zebrafish to humans covering 450 million years of evolution. Our phylointeractomics screen discovers previously unknown telomere-associated proteins and reveals how homologous proteins undergo functional evolution. For instance, we show that TERF1 evolved as a telomere-binding protein in the common stem lineage of marsupial and placental mammals. Phylointeractomics is a versatile and scalable approach to investigate evolutionary changes in protein function and thus can provide experimental evidence for phylogenomic relationships.

Reader interactome of epigenetic histone marks in birds.

Lysine methylation is part of the posttranscriptional histone code employed to recruit modification specific readers to chromatin. Unbiased, quantitative mass spectrometry approaches combined with peptide pull-downs have been used to study histone methylation-dependent binders in mammalian cells. Here, we extend the study to birds by investigating the interaction partners for H3K4me3, H3K9me3, H3K27me3 and H3K36me3 in chicken (Gallus gallus) and zebra finch (Taeniopygia guttata) using label-free quantitative proteomics. In general, we find very strong overlap in interaction partners for the trimethyl marks in birds compared to mammals, underscoring the known conserved function of these modifications. In agreement with their epigenetic role, we find binding of PHF2 and members of the TFIID, SAGA, SET1 and NURF complex to the activation mark H3K4me3. Our data furthermore supports the existence of a LID complex in vertebrates recruited to the H3K4me3 mark. The repressive marks are bound by the HP1 proteins and the EED subunit of the PRC2 complex as well as by WIZ. Like reported in the previous mammalian screens, we found ZNF462, ZNF828 and POGZ enriched at H3K9me3. However, we noted some unexpected differences. N-PAC (also known as GLYR1), an H3K36me3 interactor in mammals, is reproducible not enriched at this modification in our screen in birds. This initial finding suggests that despite strong conservation of the histone tail sequence, a few species-specific differences in epigenetic readers may have evolved between birds and mammals. All MS data have been deposited in the ProteomeXchange with identifier PXD002282 (http://proteomecentral.proteomexchange.org/dataset/PXD002282).

Identification of TTAGGG-binding proteins in Neurospora crassa, a fungus with vertebrate-like telomere repeats.

BACKGROUND: To date, telomere research in fungi has mainly focused on Saccharomyces cerevisiae and Schizosaccharomyces pombe, despite the fact that both yeasts have degenerated telomeric repeats in contrast to the canonical TTAGGG motif found in vertebrates and also several other fungi. RESULTS: Using label-free quantitative proteomics, we here investigate the telosome of Neurospora crassa, a fungus with canonical telomeric repeats. We show that at least six of the candidates detected in our screen are direct TTAGGG-repeat binding proteins. While three of the direct interactors (NCU03416 [ncTbf1], NCU01991 [ncTbf2] and NCU02182 [ncTay1]) feature the known myb/homeobox DNA interaction domain also found in the vertebrate telomeric factors, we additionally show that a zinc-finger protein (NCU07846) and two proteins without any annotated DNA-binding domain (NCU02644 and NCU05718) are also direct double-strand TTAGGG binders. We further find two single-strand binders (NCU02404 [ncGbp2] and NCU07735 [ncTcg1]). CONCLUSION: By quantitative label-free interactomics we identify TTAGGG-binding proteins in Neurospora crassa, suggesting candidates for telomeric factors that are supported by phylogenomic comparison with yeast species. Intriguingly, homologs in yeast species with degenerated telomeric repeats are also TTAGGG-binding proteins, e.g. in S. cerevisiae Tbf1 recognizes the TTAGGG motif found in its subtelomeres. However, there is also a subset of proteins that is not conserved. While a rudimentary core TTAGGG-recognition machinery may be conserved across yeast species, our data suggests Neurospora as an emerging model organism with unique features.

Quantitative interaction screen of telomeric repeat-containing RNA reveals novel TERRA regulators.

Telomeres are actively transcribed into telomeric repeat-containing RNA (TERRA), which has been implicated in the regulation of telomere length and heterochromatin formation. Here, we applied quantitative mass spectrometry (MS)-based proteomics to obtain a high-confidence interactome of TERRA. Using SILAC-labeled nuclear cell lysates in an RNA pull-down experiment and two different salt conditions, we distinguished 115 proteins binding specifically to TERRA out of a large set of background binders. While TERRA binders identified in two previous studies showed little overlap, using quantitative mass spectrometry we obtained many candidates reported in these two studies. To test whether novel candidates found here are involved in TERRA regulation, we performed an esiRNA-based interference analysis for 15 of them. Knockdown of 10 genes encoding candidate proteins significantly affected total cellular levels of TERRA, and RNAi of five candidates perturbed TERRA recruitment to telomeres. Notably, depletion of SRRT/ARS2, involved in miRNA processing, up-regulated both total and telomere-bound TERRA. Conversely, knockdown of MORF4L2, a component of the NuA4 histone acetyltransferase complex, reduced TERRA levels both globally and for telomere-bound TERRA. We thus identified new proteins involved in the homeostasis and telomeric abundance of TERRA, extending our knowledge of TERRA regulation.

Quantitative proteomic analysis reveals concurrent RNA-protein interactions and identifies new RNA-binding proteins in Saccharomyces cerevisiae.

A growing body of evidence supports the existence of an extensive network of RNA-binding proteins (RBPs) whose combinatorial binding affects the post-transcriptional fate of every mRNA in the cell-yet we still do not have a complete understanding of which proteins bind to mRNA, which of these bind concurrently, and when and where in the cell they bind. We describe here a method to identify the proteins that bind to RNA concurrently with an RBP of interest, using quantitative mass spectrometry combined with RNase treatment of affinity-purified RNA-protein complexes. We applied this method to the known RBPs Pab1, Nab2, and Puf3. Our method significantly enriched for known RBPs and is a clear improvement upon previous approaches in yeast. Our data reveal that some reported protein-protein interactions may instead reflect simultaneous binding to shared RNA targets. We also discovered more than 100 candidate RBPs, and we independently confirmed that 77% (23/30) bind directly to RNA. The previously recognized functions of the confirmed novel RBPs were remarkably diverse, and we mapped the RNA-binding region of one of these proteins, the transcriptional coactivator Mbf1, to a region distinct from its DNA-binding domain. Our results also provided new insights into the roles of Nab2 and Puf3 in post-transcriptional regulation by identifying other RBPs that bind simultaneously to the same mRNAs. While existing methods can identify sets of RBPs that interact with common RNA targets, our approach can determine which of those interactions are concurrent-a crucial distinction for understanding post-transcriptional regulation.

HOT1 is a mammalian direct telomere repeat-binding protein contributing to telomerase recruitment.

Telomeres are repetitive DNA structures that, together with the shelterin and the CST complex, protect the ends of chromosomes. Telomere shortening is mitigated in stem and cancer cells through the de novo addition of telomeric repeats by telomerase. Telomere elongation requires the delivery of the telomerase complex to telomeres through a not yet fully understood mechanism. Factors promoting telomerase-telomere interaction are expected to directly bind telomeres and physically interact with the telomerase complex. In search for such a factor we carried out a SILAC-based DNA-protein interaction screen and identified HMBOX1, hereafter referred to as homeobox telomere-binding protein 1 (HOT1). HOT1 directly and specifically binds double-stranded telomere repeats, with the in vivo association correlating with binding to actively processed telomeres. Depletion and overexpression experiments classify HOT1 as a positive regulator of telomere length. Furthermore, immunoprecipitation and cell fractionation analyses show that HOT1 associates with the active telomerase complex and promotes chromatin association of telomerase. Collectively, these findings suggest that HOT1 supports telomerase-dependent telomere elongation.

Proteome-wide analysis of disease-associated SNPs that show allele-specific transcription factor binding.

A causative role for single nucleotide polymorphisms (SNPs) in many genetic disorders has become evident through numerous genome-wide association studies. However, identification of these common causal variants and the molecular mechanisms underlying these associations remains a major challenge. Differential transcription factor binding at a SNP resulting in altered gene expression is one possible mechanism. Here we apply PWAS ("proteome-wide analysis of SNPs"), a methodology based on quantitative mass spectrometry that enables rapid screening of SNPs for differential transcription factor binding, to 12 SNPs that are highly associated with type 1 diabetes at the IL2RA locus, encoding the interleukin-2 receptor CD25. We report differential, allele-specific binding of the transcription factors RUNX1, LEF1, CREB, and TFAP4 to IL2RA SNPs rs12722508*A, rs12722522*C, rs41295061*A, and rs2104286*A and demonstrate the functional influence of RUNX1 at rs12722508 by reporter gene assay. Thus, PWAS may be able to contribute to our understanding of the molecular consequences of human genetic variability underpinning susceptibility to multi-factorial disease.

Quantitative mass spectrometry and PAR-CLIP to identify RNA-protein interactions.

Systematic analysis of the RNA-protein interactome requires robust and scalable methods. We here show the combination of two completely orthogonal, generic techniques to identify RNA-protein interactions: PAR-CLIP reveals a collection of RNAs bound to a protein whereas SILAC-based RNA pull-downs identify a group of proteins bound to an RNA. We investigated binding sites for five different proteins (IGF2BP1-3, QKI and PUM2) exhibiting different binding patterns. We report near perfect agreement between the two approaches. Nevertheless, they are non-redundant, and ideally complement each other to map the RNA-protein interaction network.

Unbiased RNA-protein interaction screen by quantitative proteomics.

Mass spectrometry (MS)-based quantitative interaction proteomics has successfully elucidated specific protein-protein, DNA-protein, and small molecule-protein interactions. Here, we developed a gel-free, sensitive, and scalable technology that addresses the important area of RNA-protein interactions. Using aptamer-tagged RNA as bait, we captured RNA-interacting proteins from stable isotope labeling by amino acids in cell culture (SILAC)-labeled mammalian cell extracts and analyzed them by high-resolution, quantitative MS. Binders specific to the RNA sequence were distinguished from background by their isotope ratios between bait and control. We demonstrated the approach by retrieving known and novel interaction partners for the HuR interaction motif, H4 stem loop, "zipcode" sequence, tRNA, and a bioinformatically-predicted RNA fold in DGCR-8/Pasha mRNA. In all experiments we unambiguously identified known interaction partners by a single affinity purification step. The 5' region of the mRNA of DGCR-8/Pasha, a component of the microprocessor complex, specifically interacts with components of the translational machinery, suggesting that it contains an internal ribosome entry site.