Transcriptome analysis in mouse skin after exposure to ultraviolet radiation from a canopy sunbed

Exposure to ultraviolet radiation (UV-R), from both natural and artificial tanning, heightens the risk of skin cancer by inducing molecular changes in cells and tissues. Despite established transcriptional alterations at a molecular level due to UV-R exposure, uncertainties persist regarding UV radiation characterization and subsequent genomic changes. Our study aimed to mechanistically explore dose- and time-dependent gene expression changes, that may drive short-term (e.g., sunburn) and long-term actinic (e.g., skin cancer) consequences. Using C57BL/6N mouse skin, we analyzed transcriptomic expression following exposure to five erythemally weighted UV-R doses (0, 5, 10, 20, and 40 mJ/cm2 ) emitted by a UV-R tanning device. At 96 h post-exposure, 5 mJ/cm2 induced 116 statistically significant differentially expressed genes (DEGs) associated with structural changes from UV-R damage. The highest number of significant gene expression changes occurred at 6 and 48 h post-exposure in the 20 and 40 mJ/cm2 dose groups. Notably, at 40 mJ/cm2 , 13 DEGs related to skin barrier homeostasis were consistently perturbed across all timepoints. UV-R exposure activated pathways involving oxidative stress, P53 signaling, inflammation, biotransformation, skin barrier maintenance, and innate immunity. This in vivo study's transcriptional data offers mechanistic insights into both short-term and potential non-threshold-dependent long-term health effects of UV-R tanning.

Measuring significant changes in chromatin conformation with ACCOST

Chromatin conformation assays such as Hi-C cannot directly measure differences in 3D architecture between cell types or cell states. For this purpose, two or more Hi-C experiments must be carried out, but direct comparison of the resulting Hi-C matrices is confounded by several features of Hi-C data. Most notably, the genomic distance effect, whereby contacts between pairs of genomic loci that are proximal along the chromosome exhibit many more Hi-C contacts that distal pairs of loci, dominates every Hi-C matrix. Furthermore, the form that this distance effect takes often varies between different Hi-C experiments, even between replicate experiments. Thus, a statistical confidence measure designed to identify differential Hi-C contacts must accurately account for the genomic distance effect or risk being misled by large-scale but artifactual differences. ACCOST (Altered Chromatin COnformation STatistics) accomplishes this goal by extending the statistical model employed by DEseq, re-purposing the ‘size factors,’ which were originally developed to account for differences in read depth between samples, to instead model the genomic distance effect. We show via analysis of simulated and real data that ACCOST provides unbiased statistical confidence estimates that compare favorably with competing methods such as diffHiC, FIND and HiCcompare. ACCOST is freely available with an Apache license at https://bitbucket.org/noblelab/accost.

Changes in genome organization of parasite-specific gene families during the Plasmodium transmission stages

The development of malaria parasites throughout their various life cycle stages is coordinated by changes in gene expression. We previously showed that the three-dimensional organization of the Plasmodium falciparum genome is strongly associated with gene expression during its replication cycle inside red blood cells. Here, we analyze genome organization in the P. falciparum and P. vivax transmission stages. Major changes occur in the localization and interactions of genes involved in pathogenesis and immune evasion, host cell invasion, sexual differentiation, and master regulation of gene expression. Furthermore, we observe reorganization of subtelomeric heterochromatin around genes involved in host cell remodeling. Depletion of heterochromatin protein 1 (PfHP1) resulted in loss of interactions between virulence genes, confirming that PfHP1 is essential for maintenance of the repressive center. Our results suggest that the three-dimensional genome structure of human malaria parasites is strongly connected with transcriptional activity of specific gene families throughout the life cycle.

RNAcompete-S: Combined RNA sequence/structure preferences for RNA binding proteins derived from a single-step in vitro selection

RNA-binding proteins recognize RNA sequences and structures, but there is currently no systematic and accurate method to derive large (>12base) motifs de novo that reflect a combination of intrinsic preference to both sequence and structure. To address this absence, we introduce RNAcompete-S, which couples a single-step competitive binding reaction with an excess of random RNA 40-mers to a custom computational pipeline for interrogation of the bound RNA sequences and derivation of SSMs (Sequence and Structure Models). RNAcompete-S confirms that HuR, QKI, and SRSF1 prefer binding sites that are single stranded, and recapitulates known 8-10bp sequence and structure preferences for Vts1p and RBMY. We also derive an 18-base long SSM for Drosophila SLBP, which to our knowledge has not been previously determined by selections from pure random sequence, and accurately discriminates human replication-dependent histone mRNAs. Thus, RNAcompete-S enables accurate identification of large, intrinsic sequence-structure specificities with a uniform assay.

Brain tumor is a sequence-specific RNA-binding protein that directs maternal mRNA clearance during the Drosophila maternal-to-zygotic transition

Background

Brain tumor (BRAT) is a Drosophila member of the TRIM-NHL protein family. This family is conserved among metazoans and its members function as post-transcriptional regulators. BRAT was thought to be recruited to mRNAs indirectly through interaction with the RNA-binding protein Pumilio (PUM). However, it has recently been demonstrated that BRAT directly binds to RNA. The precise sequence recognized by BRAT, the extent of BRAT-mediated regulation, and the exact roles of PUM and BRAT in post-transcriptional regulation are unknown.

Results

Genome-wide identification of transcripts associated with BRAT or with PUM in Drosophila embryos shows that they bind largely non-overlapping sets of mRNAs. BRAT binds mRNAs that encode proteins associated with a variety of functions, many of which are distinct from those implemented by PUM-associated transcripts. Computational analysis of in vitro and in vivo data identified a novel RNA motif recognized by BRAT that confers BRAT-mediated regulation in tissue culture cells. The regulatory status of BRAT-associated mRNAs suggests a prominent role for BRAT in post-transcriptional regulation, including a previously unidentified role in transcript degradation. Transcriptomic analysis of embryos lacking functional BRAT reveals an important role in mediating the decay of hundreds of maternal mRNAs during the maternal-to-zygotic transition.

Conclusions

Our results represent the first genome-wide analysis of the mRNAs associated with a TRIM-NHL protein and the first identification of an RNA motif bound by this protein family. BRAT is a prominent post-transcriptional regulator in the early embryo through mechanisms that are largely independent of PUM.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0659-4) contains supplementary material, which is available to authorized users.

High-throughput characterization of protein-RNA interactions

RNA-binding proteins (RBPs) are important regulators of eukaryotic gene expression. Genomes typically encode dozens to hundreds of proteins containing RNA-binding domains, which collectively recognize diverse RNA sequences and structures. Recent advances in high-throughput methods for assaying the targets of RBPs in vitro and in vivo allow large-scale derivation of RNA-binding motifs as well as determination of RNA–protein interactions in living cells. In parallel, many computational methods have been developed to analyze and interpret these data. The interplay between RNA secondary structure and RBP binding has also been a growing theme. Integrating RNA–protein interaction data with observations of post-transcriptional regulation will enhance our understanding of the roles of these important proteins.

Paralogs hnRNP L and hnRNP LL exhibit overlapping but distinct RNA binding constraints

HnRNP (heterogeneous nuclear ribonucleoprotein) proteins are a large family of RNA-binding proteins that regulate numerous aspects of RNA processing. Interestingly, several paralogous pairs of hnRNPs exist that exhibit similar RNA-binding specificity to one another, yet have non-redundant functional targets in vivo. In this study we systematically investigate the possibility that the paralogs hnRNP L and hnRNP LL have distinct RNA binding determinants that may underlie their lack of functional redundancy. Using a combination of RNAcompete and native gel analysis we find that while both hnRNP L and hnRNP LL preferentially bind sequences that contain repeated CA dinucleotides, these proteins differ in their requirement for the spacing of the CAs. Specifically, hnRNP LL has a more stringent requirement for a two nucleotide space between CA repeats than does hnRNP L, resulting in hnRNP L binding more promiscuously than does hnRNP LL. Importantly, this differential requirement for the spacing of CA dinucleotides explains the previously observed differences in the sensitivity of hnRNP L and LL to mutations within the CD45 gene. We suggest that overlapping but divergent RNA-binding preferences, as we show here for hnRNP L and hnRNP LL, may be commonplace among other hnRNP paralogs.

RBPDB: a database of RNA-binding specificities

The RNA-Binding Protein DataBase (RBPDB) is a collection of experimental observations of RNA-binding sites, both in vitro and in vivo, manually curated from primary literature. To build RBPDB, we performed a literature search for experimental binding data for all RNA-binding proteins (RBPs) with known RNA-binding domains in four metazoan species (human, mouse, fly and worm). In total, RPBDB contains binding data on 272 RBPs, including 71 that have motifs in position weight matrix format, and 36 sets of sequences of in vivo-bound transcripts from immunoprecipitation experiments. The database is accessible by a web interface which allows browsing by domain or by organism, searching and export of records, and bulk data downloads. Users can also use RBPDB to scan sequences for RBP-binding sites. RBPDB is freely available, without registration at http://rbpdb.ccbr.utoronto.ca/.