Discovery of directional and nondirectional pioneer transcription factors by modeling DNase profile magnitude and shape

We describe protein interaction quantitation (PIQ), a computational method for modeling the magnitude and shape of genome-wide DNase I hypersensitivity profiles to identify transcription factor (TF) binding sites. Through the use of machine-learning techniques, PIQ identified binding sites for >700 TFs from one DNase I hypersensitivity analysis followed by sequencing (DNase-seq) experiment with accuracy comparable to that of chromatin immunoprecipitation followed by sequencing (ChIP-seq). We applied PIQ to analyze DNase-seq data from mouse embryonic stem cells differentiating into prepancreatic and intestinal endoderm. We identified 120 and experimentally validated eight 'pioneer' TF families that dynamically open chromatin. Four pioneer TF families only opened chromatin in one direction from their motifs. Furthermore, we identified 'settler' TFs whose genomic binding is principally governed by proximity to open chromatin. Our results support a model of hierarchical TF binding in which directional and nondirectional pioneer activity shapes the chromatin landscape for population by settler TFs.National Institutes of Health (U.S.) (Common Fund 5UL1DE019581)National Institutes of Health (U.S.) (Common Fund RL1DE019021)National Institutes of Health (U.S.) (Common Fund 5TL1EB008540)National Institutes of Health (U.S.) (Grant 1U01HG007037)National Institutes of Health (U.S.) (Grant 5P01NS055923

Genome-wide association study identifies loci associated with liability to alcohol and drug dependence that is associated with variability in reward-related ventral striatum activity in African- and European-Americans.

Genetic influences on alcohol and drug dependence partially overlap, however, specific loci underlying this overlap remain unclear. We conducted a genome-wide association study (GWAS) of a phenotype representing alcohol or illicit drug dependence (ANYDEP) among 7291 European-Americans (EA; 2927 cases) and 3132 African-Americans (AA: 1315 cases) participating in the family-based Collaborative Study on the Genetics of Alcoholism. ANYDEP was heritable (h 2 in EA = 0.60, AA = 0.37). The AA GWAS identified three regions with genome-wide significant (GWS; P < 5E-08) single nucleotide polymorphisms (SNPs) on chromosomes 3 (rs34066662, rs58801820) and 13 (rs75168521, rs78886294), and an insertion-deletion on chromosome 5 (chr5:141988181). No polymorphisms reached GWS in the EA. One GWS region (chromosome 1: rs1890881) emerged from a trans-ancestral meta-analysis (EA + AA) of ANYDEP, and was attributable to alcohol dependence in both samples. Four genes (AA: CRKL, DZIP3, SBK3; EA: P2RX6) and four sets of genes were significantly enriched within biological pathways for hemostasis and signal transduction. GWS signals did not replicate in two independent samples but there was weak evidence for association between rs1890881 and alcohol intake in the UK Biobank. Among 118 AA and 481 EA individuals from the Duke Neurogenetics Study, rs75168521 and rs1890881 genotypes were associated with variability in reward-related ventral striatum activation. This study identified novel loci for substance dependence and provides preliminary evidence that these variants are also associated with individual differences in neural reward reactivity. Gene discovery efforts in non-European samples with distinct patterns of substance use may lead to the identification of novel ancestry-specific genetic markers of risk

High Energy Neutrinos from Quasars

We review and clarify the assumptions of our basic model for neutrino production in the cores of quasars, as well as those modifications to the model subsequently made by other workers. We also present a revised estimate of the neutrino background flux and spectrum obtained using more recent empirical studies of quasars and their evolution. We compare our results with other thoeretical calculations and experimental upper limits on the AGN neutrino background flux. We also estimate possible neutrino fluxes from the jets of blazars detected recently by the EGRET experiment on the Compton Gamma Ray Observatory. We discuss the theoretical implications of these estimates.Comment: 14 pg., ps file (includes figures), To be published in Space Science Review

A Model-Based Analysis of GC-Biased Gene Conversion in the Human and Chimpanzee Genomes

GC-biased gene conversion (gBGC) is a recombination-associated process that favors the fixation of G/C alleles over A/T alleles. In mammals, gBGC is hypothesized to contribute to variation in GC content, rapidly evolving sequences, and the fixation of deleterious mutations, but its prevalence and general functional consequences remain poorly understood. gBGC is difficult to incorporate into models of molecular evolution and so far has primarily been studied using summary statistics from genomic comparisons. Here, we introduce a new probabilistic model that captures the joint effects of natural selection and gBGC on nucleotide substitution patterns, while allowing for correlations along the genome in these effects. We implemented our model in a computer program, called phastBias, that can accurately detect gBGC tracts about 1 kilobase or longer in simulated sequence alignments. When applied to real primate genome sequences, phastBias predicts gBGC tracts that cover roughly 0.3% of the human and chimpanzee genomes and account for 1.2% of human-chimpanzee nucleotide differences. These tracts fall in clusters, particularly in subtelomeric regions; they are enriched for recombination hotspots and fast-evolving sequences; and they display an ongoing fixation preference for G and C alleles. They are also significantly enriched for disease-associated polymorphisms, suggesting that they contribute to the fixation of deleterious alleles. The gBGC tracts provide a unique window into historical recombination processes along the human and chimpanzee lineages. They supply additional evidence of long-term conservation of megabase-scale recombination rates accompanied by rapid turnover of hotspots. Together, these findings shed new light on the evolutionary, functional, and disease implications of gBGC. The phastBias program and our predicted tracts are freely available. © 2013 Capra et al

The clinical effectiveness of transurethral incision of the prostate : a systematic review of randomised controlled trials

The original publication is available at www.springerlink.com.Peer reviewedPostprin

Cancer somatic mutations cluster in a subset of regulatory sites predicted from the ENCODE data

Background: Transcriptional regulation of gene expression is essential for cellular differentiation and function, and defects in the process are associated with cancer. The ENCODE project has mapped potential regulatory sites across the complete genome in many cell types, and these regions have been shown to harbour many of the somatic mutations that occur in cancer cells, suggesting that their effects may drive cancer initiation and development. The ENCODE data suggests a very large number of regulatory sites, and methods are needed to identify those that are most relevant and to connect them to the genes that they control. Methods: Predictive models of gene expression were developed by integrating the ENCODE data for regulation, including transcription factor binding and DNase1 hypersensitivity, with RNA-seq data for gene expression. A penalized regression method was used to identify the most predictive potential regulatory sites for each transcript. Known cancer somatic mutations from the COSMIC database were mapped to potential regulatory sites, and we examined differences in the mapping frequencies associated with sites chosen in regulatory models and other (rejected) sites. The effects of potential confounders, for example replication timing, were considered. Results: Cancer somatic mutations preferentially occupy those regulatory regions chosen in our models as most predictive of gene expression. Conclusion: Our methods have identified a significantly reduced set of regulatory sites that are enriched in cancer somatic mutations and are more predictive of gene expression. This has significance for the mechanistic interpretation of cancer mutations, and the understanding of genetic regulation

Toxoplasma Effector MAF1 Mediates Recruitment of Host Mitochondria and Impacts the Host Response

Recent information has revealed the functional diversity and importance of mitochondria in many cellular processes including orchestrating the innate immune response. Intriguingly, several infectious agents, such as Toxoplasma, Legionella, and Chlamydia, have been reported to grow within vacuoles surrounded by host mitochondria. Although many hypotheses have been proposed for the existence of host mitochondrial association (HMA), the causes and biological consequences of HMA have remained unanswered. Here we show that HMA is present in type I and III strains of Toxoplasma but missing in type II strains, both in vitro and in vivo. Analysis of F1 progeny from a type II×III cross revealed that HMA is a Mendelian trait that we could map. We use bioinformatics to select potential candidates and experimentally identify the polymorphic parasite protein involved, mitochondrial association factor 1 (MAF1). We show that introducing the type I (HMA+) MAF1 allele into type II (HMA-) parasites results in conversion to HMA+ and deletion of MAF1 in type I parasites results in a loss of HMA. We observe that the loss and gain of HMA are associated with alterations in the transcription of host cell immune genes and the in vivo cytokine response during murine infection. Lastly, we use exogenous expression of MAF1 to show that it binds host mitochondria and thus MAF1 is the parasite protein directly responsible for HMA. Our findings suggest that association with host mitochondria may represent a novel means by which Toxoplasma tachyzoites manipulate the host. The existence of naturally occurring HMA+ and HMA- strains of Toxoplasma, Legionella, and Chlamydia indicates the existence of evolutionary niches where HMA is either advantageous or disadvantageous, likely reflecting tradeoffs in metabolism, immune regulation, and other functions of mitochondria. © 2014 Pernas et al

Circuit dissection of the role of somatostatin in itch and pain

Stimuli that elicit itch are detected by sensory neurons that innervate the skin. This information is processed by the spinal cord; however, the way in which this occurs is still poorly understood. Here we investigated the neuronal pathways for itch neurotransmission, particularly the contribution of the neuropeptide somatostatin. We find that in the periphery, somatostatin is exclusively expressed in Nppb+ neurons, and we demonstrate that Nppb+somatostatin+ cells function as pruriceptors. Employing chemogenetics, pharmacology and cell-specific ablation methods, we demonstrate that somatostatin potentiates itch by inhibiting inhibitory dynorphin neurons, which results in disinhibition of GRPR+ neurons. Furthermore, elimination of somatostatin from primary afferents and/or from spinal interneurons demonstrates differential involvement of the peptide released from these sources in itch and pain. Our results define the neural circuit underlying somatostatin-induced itch and characterize a contrasting antinociceptive role for the peptide

Connecting the dots: Potential of data integration to identify regulatory snps in late-onset alzheimer's disease GWAS findings

Late-onset Alzheimer's disease (LOAD) is a multifactorial disorder with over twenty loci associated with disease risk. Given the number of genome-wide significant variants that fall outside of coding regions, it is possible that some of these variants alter some function of gene expression rather than tagging coding variants that alter protein structure and/or function. RegulomeDB is a database that annotates regulatory functions of genetic variants. In this study, we utilized RegulomeDB to investigate potential regulatory functions of lead single nucleotide polymorphisms (SNPs) identified in five genome-wide association studies (GWAS) of risk and age-at onset (AAO) of LOAD, as well as SNPs in LD (r2≥0.80) with the lead GWAS SNPs. Of a total 614 SNPs examined, 394 returned RegulomeDB scores of 1-6. Of those 394 variants, 34 showed strong evidence of regulatory function (RegulomeDB score ,3), and only 3 of them were genome-wide significant SNPs (ZCWPW1/ rs1476679, CLU/rs1532278 and ABCA7/rs3764650). This study further supports the assumption that some of the non-coding GWAS SNPs are true associations rather than tagged associations and demonstrates the application of RegulomeDB to GWAS data.©2014 Rosenthal et al

Discovery and characterization of chromatin states for systematic annotation of the human genome

A plethora of epigenetic modifications have been described in the human genome and shown to play diverse roles in gene regulation, cellular differentiation and the onset of disease. Although individual modifications have been linked to the activity levels of various genetic functional elements, their combinatorial patterns are still unresolved and their potential for systematic de novo genome annotation remains untapped. Here, we use a multivariate Hidden Markov Model to reveal 'chromatin states' in human T cells, based on recurrent and spatially coherent combinations of chromatin marks. We define 51 distinct chromatin states, including promoter-associated, transcription-associated, active intergenic, large-scale repressed and repeat-associated states. Each chromatin state shows specific enrichments in functional annotations, sequence motifs and specific experimentally observed characteristics, suggesting distinct biological roles. This approach provides a complementary functional annotation of the human genome that reveals the genome-wide locations of diverse classes of epigenetic function.National Science Foundation (U.S.). (Award 0905968)National Human Genome Research Institute (U.S.) (Award U54-HG004570)National Human Genome Research Institute (U.S.) (Award RC1-HG005334