PRAS: Predicting functional targets of RNA binding proteins based on CLIP-seq peaks.

RNA-protein interaction plays important roles in post-transcriptional regulation. Recent advancements in cross-linking and immunoprecipitation followed by sequencing (CLIP-seq) technologies make it possible to detect the binding peaks of a given RNA binding protein (RBP) at transcriptome scale. However, it is still challenging to predict the functional consequences of RBP binding peaks. In this study, we propose the Protein-RNA Association Strength (PRAS), which integrates the intensities and positions of the binding peaks of RBPs for functional mRNA targets prediction. We illustrate the superiority of PRAS over existing approaches on predicting the functional targets of two related but divergent CELF (CUGBP, ELAV-like factor) RBPs in mouse brain and muscle. We also demonstrate the potential of PRAS for wide adoption by applying it to the enhanced CLIP-seq (eCLIP) datasets of 37 RNA decay related RBPs in two human cell lines. PRAS can be utilized to investigate any RBPs with available CLIP-seq peaks. PRAS is freely available at http://ouyanglab.jax.org/pras/

Temporal Regulation of Foregut Development by HTZ-1/H2A.Z and PHA-4/FoxA

The histone variant H2A.Z is evolutionarily conserved and plays an essential role in mice, Drosophila, and Tetrahymena. The essential function of H2A.Z is unknown, with some studies suggesting a role in transcriptional repression and others in activation. Here we show that Caenorhabditis elegans HTZ-1/H2A.Z and the remodeling complex MYS-1/ESA1–SSL-1/SWR1 synergize with the FoxA transcription factor PHA-4 to coordinate temporal gene expression during foregut development. We observe dramatic genetic interactions between pha-4 and htz-1, mys-1, and ssl-1. A survey of transcription factors reveals that this interaction is specific, and thus pha-4 is acutely sensitive to reductions in these three proteins. Using a nuclear spot assay to visualize HTZ-1 in living embryos as organogenesis proceeds, we show that HTZ-1 is recruited to foregut promoters at the time of transcriptional onset, and this recruitment requires PHA-4. Loss of htz-1 by RNAi is lethal and leads to delayed expression of a subset of foregut genes. Thus, the effects of PHA-4 on temporal regulation can be explained in part by recruitment of HTZ-1 to target promoters. We suggest PHA-4 and HTZ-1 coordinate temporal gene expression by modulating the chromatin environment

Local Regulatory Variation in Saccharomyces cerevisiae

Naturally occurring sequence variation that affects gene expression is an important source of phenotypic differences among individuals within a species. We and others have previously shown that such regulatory variation can occur both at the same locus as the gene whose expression it affects (local regulatory variation) and elsewhere in the genome at trans-acting factors. Here we present a detailed analysis of genome-wide local regulatory variation in Saccharomyces cerevisiae. We used genetic linkage analysis to show that nearly a quarter of all yeast genes contain local regulatory variation between two divergent strains. We measured allele-specific expression in a diploid hybrid of the two strains for 77 genes showing strong self-linkage and found that in 52%–78% of these genes, local regulatory variation acts directly in cis. We also experimentally confirmed one example in which local regulatory variation in the gene AMN1 acts in trans through a feedback loop. Genome-wide sequence analysis revealed that genes subject to local regulatory variation show increased polymorphism in the promoter regions, and that some but not all of this increase is due to polymorphisms in predicted transcription factor binding sites. Increased polymorphism was also found in the 3′ untranslated regions of these genes. These findings point to the importance of cis-acting variation, but also suggest that there is a diverse set of mechanisms through which local variation can affect gene expression levels

Hsp90 Selectively Modulates Phenotype in Vertebrate Development

Compromised heat shock protein 90 (Hsp90) function reveals cryptic phenotypes in flies and plants. These observations were interpreted to suggest that this molecular stress-response chaperone has a capacity to buffer underlying genetic variation. Conversely, the protective role of Hsp90 could account for the variable penetrance or severity of some heritable developmental malformations in vertebrates. Using zebrafish as a model, we defined Hsp90 inhibitor levels that did not induce a heat shock response or perturb phenotype in wild-type strains. Under these conditions the severity of the recessive eye phenotype in sunrise, caused by a pax6b mutation, was increased, while in dreumes, caused by a sufu mutation, it was decreased. In another strain, a previously unobserved spectrum of severe structural eye malformations, reminiscent of anophthalmia, microphthalmia, and nanophthalmia complex in humans, was uncovered by this limited inhibition of Hsp90 function. Inbreeding of offspring from selected unaffected carrier parents led to significantly elevated malformation frequencies and revealed the oligogenic nature of this phenotype. Unlike in Drosophila, Hsp90 inhibition can decrease developmental stability in zebrafish, as indicated by increased asymmetric presentation of anophthalmia, microphthalmia, and nanophthalmia and sunrise phenotypes. Analysis of the sunrise pax6b mutation suggests a molecular mechanism for the buffering of mutations by Hsp90. The zebrafish studies imply that mild perturbation of Hsp90 function at critical developmental stages may underpin the variable penetrance and expressivity of many developmental anomalies where the interaction between genotype and environment plays a major role

Quantitative Genomics of Aggressive Behavior in Drosophila melanogaster

Aggressive behavior is important for animal survival and reproduction, and excessive aggression is an enormous social and economic burden for human society. Although the role of biogenic amines in modulating aggressive behavior is well characterized, other genetic mechanisms affecting this complex behavior remain elusive. Here, we developed an assay to rapidly quantify aggressive behavior in Drosophila melanogaster, and generated replicate selection lines with divergent levels of aggression. The realized heritability of aggressive behavior was approximately 0.10, and the phenotypic response to selection specifically affected aggression. We used whole-genome expression analysis to identify 1,539 probe sets with different expression levels between the selection lines when pooled across replicates, at a false discovery rate of 0.001. We quantified the aggressive behavior of 19 mutations in candidate genes that were generated in a common co-isogenic background, and identified 15 novel genes affecting aggressive behavior. Expression profiling of genetically divergent lines is an effective strategy for identifying genes affecting complex traits

Role of APOBEC3 in Genetic Diversity among Endogenous Murine Leukemia Viruses

The ability of human and murine APOBECs (specifically, APOBEC3) to inhibit infecting retroviruses and retrotransposition of some mobile elements is becoming established. Less clear is the effect that they have had on the establishment of the endogenous proviruses resident in the human and mouse genomes. We used the mouse genome sequence to study diversity and genetic traits of nonecotropic murine leukemia viruses (polytropic [Pmv], modified polytropic [Mpmv], and xenotropic [Xmv] subgroups), the best-characterized large set of recently integrated proviruses. We identified 49 proviruses. In phylogenetic analyses, Pmvs and Mpmvs were monophyletic, whereas Xmvs were divided into several clades, implying a greater number of replication cycles between the integration events. Four distinct primer binding site types (Pro, Gln1, Gln2 and Thr) were dispersed within the phylogeny, indicating frequent mispriming. We analyzed the frequency and context of G-to-A mutations for the role of mA3 in formation of these proviruses. In the Pmv and Mpmv (but not Xmv) groups, mutations attributable to mA3 constituted a large fraction of the total. A significant number of nonsense mutations suggests the absence of purifying selection following mutation. A strong bias of G-to-A relative to C-to-T changes was seen, implying a strand specificity that can only have occurred prior to integration. The optimal sequence context of G-to-A mutations, TTC, was consistent with mA3. At least in the Pmv group, a significant 5′ to 3′ gradient of G-to-A mutations was consistent with mA3 editing. Altogether, our results for the first time suggest mA3 editing immediately preceding the integration event that led to retroviral endogenization, contributing to inactivation of infectivity

Elevated Rates of Sister Chromatid Exchange at Chromosome Ends

Chromosome ends are known hotspots of meiotic recombination and double-strand breaks. We monitored mitotic sister chromatid exchange (SCE) in telomeres and subtelomeres and found that 17% of all SCE occurs in the terminal 0.1% of the chromosome. Telomeres and subtelomeres are significantly enriched for SCEs, exhibiting rates of SCE per basepair that are at least 1,600 and 160 times greater, respectively, than elsewhere in the genome