Genetically tractable fruit flies have been used for decades to study seizure disorders. However, there is a paucity of data specifically correlating fly and human seizure phenotypes. We have previously shown that mutation of orthologous PRICKLE genes from flies to humans produce seizures. This study aimed to determine whether the prickle-mediated seizure phenotypes in flies closely parallel the epilepsy syndrome found in PRICKLE patients.

Genome-wide association studies (GWAS) have revolutionized the field of cancer genetics, but the causal links between increased genetic risk and onset/progression of disease processes remain to be identified. Here we report the first step in such an endeavor for prostate cancer. We provide a comprehensive annotation of the 77 known risk loci, based upon highly correlated variants in biologically relevant chromatin annotations--we identified 727 such potentially functional SNPs. We also provide a detailed account of possible protein disruption, microRNA target sequence disruption and regulatory response element disruption of all correlated SNPs at r(2) ≥ 0.88%. 88% of the 727 SNPs fall within putative enhancers, and many alter critical residues in the response elements of transcription factors known to be involved in prostate biology. We define as risk enhancers those regions with enhancer chromatin biofeatures in prostate-derived cell lines with prostate-cancer correlated SNPs. To aid the identification of these enhancers, we performed genomewide ChIP-seq for H3K27-acetylation, a mark of actively engaged enhancers, as well as the transcription factor TCF7L2. We analyzed in depth three variants in risk enhancers, two of which show significantly altered androgen sensitivity in LNCaP cells. This includes rs4907792, that is in linkage disequilibrium (r(2) = 0.91) with an eQTL for NUDT11 (on the X chromosome) in prostate tissue, and rs10486567, the index SNP in intron 3 of the JAZF1 gene on chromosome 7. Rs4907792 is within a critical residue of a strong consensus androgen response element that is interrupted in the protective allele, resulting in a 56% decrease in its androgen sensitivity, whereas rs10486567 affects both NKX3-1 and FOXA-AR motifs where the risk allele results in a 39% increase in basal activity and a 28% fold-increase in androgen stimulated enhancer activity. Identification of such enhancer variants and their potential target genes represents a preliminary step in connecting risk to disease process.

Breast Cancer (BCa) genome-wide association studies revealed allelic frequency differences between cases and controls at index single nucleotide polymorphisms (SNPs). To date, 71 loci have thus been identified and replicated. More than 320,000 SNPs at these loci define BCa risk due to linkage disequilibrium (LD). We propose that BCa risk resides in a subgroup of SNPs that functionally affects breast biology. Such a shortlist will aid in framing hypotheses to prioritize a manageable number of likely disease-causing SNPs. We extracted all the SNPs, residing in 1 Mb windows around breast cancer risk index SNP from the 1000 genomes project to find correlated SNPs. We used FunciSNP, an R/Bioconductor package developed in-house, to identify potentially functional SNPs at 71 risk loci by coinciding them with chromatin biofeatures. We identified 1,005 SNPs in LD with the index SNPs (r(2)≥0.5) in three categories; 21 in exons of 18 genes, 76 in transcription start site (TSS) regions of 25 genes, and 921 in enhancers. Thirteen SNPs were found in more than one category. We found two correlated and predicted non-benign coding variants (rs8100241 in exon 2 and rs8108174 in exon 3) of the gene, ANKLE1. Most putative functional LD SNPs, however, were found in either epigenetically defined enhancers or in gene TSS regions. Fifty-five percent of these non-coding SNPs are likely functional, since they affect response element (RE) sequences of transcription factors. Functionality of these SNPs was assessed by expression quantitative trait loci (eQTL) analysis and allele-specific enhancer assays. Unbiased analyses of SNPs at BCa risk loci revealed new and overlooked mechanisms that may affect risk of the disease, thereby providing a valuable resource for follow-up studies.

In past years, much attention has focused on the gene networks that regulate early developmental processes, but less attention has been paid to how multiple networks and processes are temporally coordinated. Recently the discovery of the transcriptional activator Zelda (Zld), which binds to CAGGTAG and related sequences present in the enhancers of many early-activated genes in Drosophila, hinted at a mechanism for how batteries of genes could be simultaneously activated. Here we use genome-wide binding and expression assays to identify Zld target genes in the early embryo with the goal of unraveling the gene circuitry regulated by Zld. We found that Zld binds to genes involved in early developmental processes such as cellularization, sex determination, neurogenesis, and pattern formation. In the absence of Zld, many target genes failed to be activated, while others, particularly the patterning genes, exhibited delayed transcriptional activation, some of which also showed weak and/or sporadic expression. These effects disrupted the normal sequence of patterning-gene interactions and resulted in highly altered spatial expression patterns, demonstrating the significance of a timing mechanism in early development. In addition, we observed prevalent overlap between Zld-bound regions and genomic "hotspot" regions, which are bound by many developmental transcription factors, especially the patterning factors. This, along with the finding that the most over-represented motif in hotspots, CAGGTA, is the Zld binding site, implicates Zld in promoting hotspot formation. We propose that Zld promotes timely and robust transcriptional activation of early-gene networks so that developmental events are coordinated and cell fates are established properly in the cellular blastoderm embryo.

Expression profiling of T24 cells revealed that 17 out of 313 human miRNAs were upregulated more than 3-fold by simultaneous treatment with the chromatin-modifying drugs 5-aza-2'-deoxycytidine and 4-phenylbutyric acid. One of these, miR-127, is embedded in a CpG island and is highly induced from its own promoter after treatment. miR-127 is usually expressed as part of a miRNA cluster in normal cells but not in cancer cells, suggesting that it is subject to epigenetic silencing. In addition, the proto-oncogene BCL6, a potential target of miR-127, was translationally downregulated after treatment. These results suggest that DNA demethylation and histone deacetylase inhibition can activate expression of miRNAs that may act as tumor suppressors.

Multiple discrete regions at 8q24 were recently shown to contain alleles that predispose to many cancers including prostate, breast, and colon. These regions are far from any annotated gene and their biological activities have been unknown. Here we profiled a 5-megabase chromatin segment encompassing all the risk regions for RNA expression, histone modifications, and locations occupied by RNA polymerase II and androgen receptor (AR). This led to the identification of several transcriptional enhancers, which were verified using reporter assays. Two enhancers in one risk region were occupied by AR and responded to androgen treatment; one contained a single nucleotide polymorphism (rs11986220) that resides within a FoxA1 binding site, with the prostate cancer risk allele facilitating both stronger FoxA1 binding and stronger androgen responsiveness. The study reported here exemplifies an approach that may be applied to any risk-associated allele in non-protein coding regions as it emerges from genome-wide association studies to better understand the genetic predisposition of complex diseases.

Alpha-synuclein accumulation in dopaminergic neurons is one of the primary features of Parkinson's disease (PD). Despite its toxic properties during PD, alpha-synuclein has some important physiological functions. Although the activity of the protein has been extensively studied in neurons, the protein is also expressed in other cell types including immune cells and glia. Genetic studies show that mutations in synuclein alpha (SNCA), the gene that encodes alpha-synuclein, and alterations in its expression levels are a significant risk factor for PD, which likely impact the functions of a broad range of cell types. The consequences of altered SNCA expression in other cell types is beginning to be explored. Microglia, the primary macrophage population in the Central Nervous System (CNS), for example, are affected by variations in alpha-synuclein levels and functions. Studies suggest that deviations of alpha-synuclein's normal activity influence hematopoiesis, the process that gives rise to microglia, and microglia's immune functions. Alpha-synuclein levels also dictate the efficiency of SNARE-mediated vesicle formation, which could influence autophagy and cytokine release in microglia. Starting from the time of conception, these effects could impact one's risk for developing PD. Further studies are needed to determine the physiological role of alpha-synuclein and how the protein is affected during PD in non-neuronal cells such as microglia. In this review we will discuss the known roles of alpha-synuclein in differentiation, immune responses, and vesicle formation, with insights into how abnormal alpha-synuclein expression and activity are linked to altered functions of microglia during PD.

Proper gene regulation is critical for both neuronal development and maintenance as the brain matures. We previously demonstrated that Akirin2, an essential nuclear protein that interacts with transcription factors and chromatin remodeling complexes, is required for the embryonic formation of the cerebral cortex. Here we show that Akirin2 plays a mechanistically distinct role in maintaining healthy neurons during cortical maturation. Restricting Akirin2 loss to excitatory cortical neurons resulted in progressive neurodegeneration via necroptosis and severe cortical atrophy with age. Comparing transcriptomes from Akirin2-null postnatal neurons and cortical progenitors revealed that targets of the tumor suppressor p53, a regulator of both proliferation and cell death encoded by Trp53, were consistently upregulated. Reduction of Trp53 rescued neurodegeneration in Akirin2-null neurons. These data: (1) implicate Akirin2 as a critical neuronal maintenance protein, (2) identify p53 pathways as mediators of Akirin2 functions, and (3) suggest Akirin2 dysfunction may be relevant to neurodegenerative diseases.

We used paired-end next generation sequencing (NGS) to characterize the classic isoform-specific pk pk1 and pk sple1 mutations of the prickle gene in Drosophila melanogaster . Here we provide evidence that these previously reported null mutations are caused by either a tirant transposon insertion into the 5' UTR of pk pk1 or a premature stop codon in the second exon of pk sple1 . Additional likely benign missense mutations were identified in both mutant isoforms.

Vesicoureteral reflux (VUR) is a complex, heritable disorder. Genome-wide linkage analyses of families affected by VUR have revealed multiple genomic loci linked to VUR. These loci normally harbor a number of genes whose biologically functional variant is yet to be identified. DNA copy number variations (CNVs) have not been extensively studied at high resolution in VUR patients. In this study, we performed array comparative genomic hybridization (aCGH) on a cohort of patients with a history of both VUR and urinary tract infection (UTI) with the objective of identifying genetic variations responsible for VUR and/or UTI susceptibility. UTI/VUR-associated CNVs were identified by aCGH results from the 192 Randomized Intervention for Children With Vesicoureteral Reflux (RIVUR) patients compared to 683 controls. Rare, large CNVs that are likely pathogenic and lead to VUR development were identified using stringent analysis criteria. Because UTI is a common affliction with multiple risk factors, we utilized standard analysis to identify potential disease-modifying CNVs that can contribute to UTI risk. Gene ontology analysis identified that CNVs in innate immunity and development genes were enriched in RIVUR patients. CNVs affecting innate immune genes may contribute to UTI susceptibility in VUR patients and may provide the first step in assisting clinical medicine in determining adverse outcome risk in children with VUR.

Single nucleotide polymorphisms (SNPs) are increasingly used to tag genetic loci associated with phenotypes such as risk of complex diseases. Technically, this is done genome-wide without prior restriction or knowledge of biological feasibility in scans referred to as genome-wide association studies (GWAS). Depending on the linkage disequilibrium (LD) structure at a particular locus, such tagSNPs may be surrogates for many thousands of other SNPs, and it is difficult to distinguish those that may play a functional role in the phenotype from those simply genetically linked. Because a large proportion of tagSNPs have been identified within non-coding regions of the genome, distinguishing functional from non-functional SNPs has been an even greater challenge. A strategy was recently proposed that prioritizes surrogate SNPs based on non-coding chromatin and epigenomic mapping techniques that have become feasible with the advent of massively parallel sequencing. Here, we introduce an R/Bioconductor software package that enables the identification of candidate functional SNPs by integrating information from tagSNP locations, lists of linked SNPs from the 1000 genomes project and locations of chromatin features which may have functional significance.

Polycomb-group (PcG) and Trithorax-group proteins together form a maintenance machinery that is responsible for stable heritable states of gene activity. While the best-studied target genes are the Hox genes of the Antennapedia and Bithorax complexes, a large number of key developmental genes are also Polycomb (Pc) targets, indicating a widespread role for this maintenance machinery in cell fate determination. We have studied the linkage between the binding of PcG proteins and the developmental regulation of gene expression using whole-genome mapping to identify sites bound by the PcG proteins, Pc and Pleiohomeotic (Pho), in the Drosophila embryo and in a more restricted tissue, the imaginal discs of the third thoracic segment. Our data provide support for the idea that Pho is a general component of the maintenance machinery, since the majority of Pc targets are also associated with Pho binding. We find, in general, considerable developmental stability of Pc and Pho binding at target genes and observe that Pc/Pho binding can be associated with both expressed and inactive genes. In particular, at the Hox complexes, both active and inactive genes have significant Pc and Pho binding. However, in comparison to inactive genes, the active Hox genes show reduced and altered binding profiles. During development, Pc target genes are not simply constantly associated with Pc/Pho binding, and we identify sets of genes with clear differential binding between embryo and imaginal disc. Using existing datasets, we show that for specific fate-determining genes of the haemocyte lineage, the active state is characterised by lack of Pc binding. Overall, our analysis suggests a dynamic relationship between Pc/Pho binding and gene transcription. Pc/Pho binding does not preclude transcription, but levels of Pc/Pho binding change during development, and loss of Pc/Pho binding can be associated with both stable gene activity and inactivity.

Recent genome-wide association studies (GWAS) have identified more than 100 loci associated with increased risk of prostate cancer, most of which are in non-coding regions of the genome. Understanding the function of these non-coding risk loci is critical to elucidate the genetic susceptibility to prostate cancer.

Epilepsy is a common disabling disease with complex, multifactorial genetic and environmental etiology. The small fraction of epilepsies subject to Mendelian inheritance offers key insight into epilepsy disease mechanisms; and pathologies brought on by mutations in a single gene can point the way to generalizable therapeutic strategies. Mutations in the PRICKLE genes can cause seizures in humans, zebrafish, mice, and flies, suggesting the seizure-suppression pathway is evolutionarily conserved. This pathway has never been targeted for novel anti-seizure treatments. Here, the mammalian PRICKLE-interactome was defined, identifying prickle-interacting proteins that localize to synapses and a novel interacting partner, USP9X, a substrate-specific de-ubiquitinase. PRICKLE and USP9X interact through their carboxy-termini; and USP9X de-ubiquitinates PRICKLE, protecting it from proteasomal degradation. In forebrain neurons of mice, USP9X deficiency reduced levels of Prickle2 protein. Genetic analysis suggests the same pathway regulates Prickle-mediated seizures. The seizure phenotype was suppressed in prickle mutant flies by the small-molecule USP9X inhibitor, Degrasyn/WP1130, or by reducing the dose of fat facets a USP9X orthologue. USP9X mutations were identified by resequencing a cohort of patients with epileptic encephalopathy, one patient harbored a de novo missense mutation and another a novel coding mutation. Both USP9X variants were outside the PRICKLE-interacting domain. These findings demonstrate that USP9X inhibition can suppress prickle-mediated seizure activity, and that USP9X variants may predispose to seizures. These studies point to a new target for anti-seizure therapy and illustrate the translational power of studying diseases in species across the evolutionary spectrum.

Genome-wide association studies (GWAS) have identified numerous common prostate cancer (PrCa) susceptibility loci. We have fine-mapped 64 GWAS regions known at the conclusion of the iCOGS study using large-scale genotyping and imputation in 25 723 PrCa cases and 26 274 controls of European ancestry. We detected evidence for multiple independent signals at 16 regions, 12 of which contained additional newly identified significant associations. A single signal comprising a spectrum of correlated variation was observed at 39 regions; 35 of which are now described by a novel more significantly associated lead SNP, while the originally reported variant remained as the lead SNP only in 4 regions. We also confirmed two association signals in Europeans that had been previously reported only in East-Asian GWAS. Based on statistical evidence and linkage disequilibrium (LD) structure, we have curated and narrowed down the list of the most likely candidate causal variants for each region. Functional annotation using data from ENCODE filtered for PrCa cell lines and eQTL analysis demonstrated significant enrichment for overlap with bio-features within this set. By incorporating the novel risk variants identified here alongside the refined data for existing association signals, we estimate that these loci now explain ∼38.9% of the familial relative risk of PrCa, an 8.9% improvement over the previously reported GWAS tag SNPs. This suggests that a significant fraction of the heritability of PrCa may have been hidden during the discovery phase of GWAS, in particular due to the presence of multiple independent signals within the same region.

Although technological advances now allow increased tumor profiling, a detailed understanding of the mechanisms leading to the development of different cancers remains elusive. Our approach toward understanding the molecular events that lead to cancer is to characterize changes in transcriptional regulatory networks between normal and tumor tissue. Because enhancer activity is thought to be critical in regulating cell fate decisions, we have focused our studies on distal regulatory elements and transcription factors that bind to these elements.

Mutations in the gene encoding transcription factor TFAP2A result in pigmentation anomalies in model organisms and premature hair graying in humans. However, the pleiotropic functions of TFAP2A and its redundantly-acting paralogs have made the precise contribution of TFAP2-type activity to melanocyte differentiation unclear. Defining this contribution may help to explain why TFAP2A expression is reduced in advanced-stage melanoma compared to benign nevi. To identify genes with TFAP2A-dependent expression in melanocytes, we profile zebrafish tissue and mouse melanocytes deficient in Tfap2a, and find that expression of a small subset of genes underlying pigmentation phenotypes is TFAP2A-dependent, including Dct, Mc1r, Mlph, and Pmel. We then conduct TFAP2A ChIP-seq in mouse and human melanocytes and find that a much larger subset of pigmentation genes is associated with active regulatory elements bound by TFAP2A. These elements are also frequently bound by MITF, which is considered the "master regulator" of melanocyte development. For example, the promoter of TRPM1 is bound by both TFAP2A and MITF, and we show that the activity of a minimal TRPM1 promoter is lost upon deletion of the TFAP2A binding sites. However, the expression of Trpm1 is not TFAP2A-dependent, implying that additional TFAP2 paralogs function redundantly to drive melanocyte differentiation, which is consistent with previous results from zebrafish. Paralogs Tfap2a and Tfap2b are both expressed in mouse melanocytes, and we show that mouse embryos with Wnt1-Cre-mediated deletion of Tfap2a and Tfap2b in the neural crest almost completely lack melanocytes but retain neural crest-derived sensory ganglia. These results suggest that TFAP2 paralogs, like MITF, are also necessary for induction of the melanocyte lineage. Finally, we observe a genetic interaction between tfap2a and mitfa in zebrafish, but find that artificially elevating expression of tfap2a does not increase levels of melanin in mitfa hypomorphic or loss-of-function mutants. Collectively, these results show that TFAP2 paralogs, operating alongside lineage-specific transcription factors such as MITF, directly regulate effectors of terminal differentiation in melanocytes. In addition, they suggest that TFAP2A activity, like MITF activity, has the potential to modulate the phenotype of melanoma cells.

Prostate cancer (PCa) is the second-leading cause of cancer-related mortality, after lung cancer, in men from developed countries. In its early stages, primary tumor growth is dependent on androgens, thus generally can be controlled by androgen deprivation therapy (ADT). Eventually however, the disease progresses to castration-resistant prostate cancer (CRPC), a lethal form in need of more effective treatments. G-protein coupled receptors (GPCRs) comprise a large clan of cell surface proteins that have been implicated as therapeutic targets in PCa growth and progression. The findings reported here provide intriguing evidence of a role for the newly characterized glutamate family member GPR158 in PCa growth and progression. We found that GPR158 promotes PCa cell proliferation independent of androgen receptor (AR) functionality and that this requires its localization in the nucleus of the cell. This suggests that GPR158 acts by mechanisms different from other GPCRs. GPR158 expression is stimulated by androgens and GPR158 stimulates AR expression, implying a potential to sensitize tumors to low androgen conditions during ADT via a positive feedback loop. Further, we found GPR158 expression correlates with a neuroendocrine (NE) differentiation phenotype and promotes anchorage-independent colony formation implying a role for GPR158 in therapeutic progression and tumor formation. GPR158 expression was increased at the invading front of prostate tumors that formed in the genetically defined conditional Pten knockout mouse model, and co-localized with elevated AR expression in the cell nucleus. Kaplan-Meier analysis on a dataset from the Memorial Sloan Kettering cancer genome portal showed that increased GPR158 expression in tumors is associated with lower disease-free survival. Our findings strongly suggest that pharmaceuticals targeting GPR158 activities could represent a novel and innovative approach to the prevention and management of CRPC.

We report the development of OikoBase (http://oikoarrays.biology.uiowa.edu/Oiko/), a tiling array-based genome browser resource for Oikopleura dioica, a metazoan belonging to the urochordates, the closest extant group to vertebrates. OikoBase facilitates retrieval and mining of a variety of useful genomics information. First, it includes a genome browser which interrogates 1260 genomic sequence scaffolds and features gene, transcript and CDS annotation tracks. Second, we annotated gene models with gene ontology (GO) terms and InterPro domains which are directly accessible in the browser with links to their entries in the GO (http://www.geneontology.org/) and InterPro (http://www.ebi.ac.uk/interpro/) databases, and we provide transcript and peptide links for sequence downloads. Third, we introduce the transcriptomics of a comprehensive set of developmental stages of O. dioica at high resolution and provide downloadable gene expression data for all developmental stages. Fourth, we incorporate a BLAST tool to identify homologs of genes and proteins. Finally, we include a tutorial that describes how to use OikoBase as well as a link to detailed methods, explaining the data generation and analysis pipeline. OikoBase will provide a valuable resource for research in chordate development, genome evolution and plasticity and the molecular ecology of this important marine planktonic organism.

Paraoxonases (PON) are a family of proteins (PON1, 2 and 3) with multiple enzymatic activities. PON1 interferes with homoserine lactone-mediated quorum sensing in bacteria and with reactive oxygen species (ROS) in humans and mice. PON1 gene mutations have been linked to multiple traits, including aging, and diseases of the cardiovascular, nervous and gastrointestinal system. The overlapping enzymatic activities in the PON family members and high linkage disequilibrium rates within their polymorphisms confound animal and human studies of PON1 function. In contrast, arthropods such as Drosophila melanogaster have no PON homologs, resulting in an ideal model to study interactions between PON genotype and host phenotypes. We hypothesized that expression of PON1 in D. melanogaster would alter ROS. We found that PON1 alters expression of multiple oxidative stress genes and decreases superoxide anion levels in normal and germ-free D. melanogaster. We also found differences in the composition of the gut microbiota, with a remarkable increase in levels of Lactobacillus plantarum and associated changes in expression of antimicrobial and cuticle-related genes. PON1 expression directly decreased superoxide anion levels and altered bacterial colonization of the gut and its gene expression profile, highlighting the complex nature of the interaction between host genotype and gut microbiota. We speculate that the interaction between some genotypes and human diseases may be mediated by the presence of certain gut bacteria that can induce specific immune responses in the gut and other host tissues.