Species-specific transcription in mice carrying human chromosome 21

The androgen receptor induces a distinct transcriptional program in castration-resistant prostate cancer in man

The androgen receptor (AR) regulates prostate cell growth in man, and prostate cancer is the commonest cancer in men in the UK. We present a comprehensive analysis of AR binding sites in human prostate cancer tissues, including castrate-resistant prostate cancer (CRPC). We identified thousands of AR binding sites in CRPC tissue, most of which were not identified in PC cell lines. Many adjacent genes showed AR regulation in xenografts but not in cultured LNCaPs, demonstrating an in-vivo-restricted set of AR-regulated genes. Functional studies support a model of altered signaling in vivo that directs AR binding. We identified a 16 gene signature that outperformed a larger in-vitro-derived signature in clinical data sets, showing the importance of persistent AR signaling in CRPC.

The evolutionary landscape of alternative splicing in vertebrate species

How species with similar repertoires of protein-coding genes differ so markedly at the phenotypic level is poorly understood. By comparing organ transcriptomes from vertebrate species spanning ~350 million years of evolution, we observed significant differences in alternative splicing complexity between vertebrate lineages, with the highest complexity in primates. Within 6 million years, the splicing profiles of physiologically equivalent organs diverged such that they are more strongly related to the identity of a species than they are to organ type. Most vertebrate species-specific splicing patterns are cis-directed. However, a subset of pronounced splicing changes are predicted to remodel protein interactions involving trans-acting regulators. These events likely further contributed to the diversification of splicing and other transcriptomic changes that underlie phenotypic differences among vertebrate species.

Identification and analysis of two novel sites of rat GnRH receptor gene promoter activity: the pineal gland and retina

BACKGROUND AND AIMS

In mammals, activation of pituitary GnRH receptor (GnRHR) by hypothalamic GnRH increases the synthesis and secretion of LH and FSH, which, in turn, regulate gonadal functions. However, GnRHR gene (Gnrhr) expression is not restricted to the pituitary.

METHODS

To gain insight into the extrapituitary expression of Gnrhr, a transgenic mouse model that expresses the human placental alkaline phosphatase reporter gene driven by the rat Gnrhr promoter was created.

RESULTS

This study shows that the rat Gnrhr promoter is operative in two functionally related organs, the pineal gland, as early as embryonic day (E) 13.5, and the retina where activity was only detected at E17.5. Accordingly, Gnrhr mRNA were present in both tissues. Transcription factors known to regulate Gnrhr promoter activity such as the LIM homeodomain factors LHX3 and ISL1 were also detected in the retina. Furthermore, transient transfection studies in CHO and gonadotrope cells revealed that OTX2, a major transcription factor in both pineal and retina cell differentiation, is able to activate the Gnrhr promoter together with either CREB or PROP1, depending on the cell context.

CONCLUSION

Rather than using alternate promoters, Gnrhr expression is directed to diverse cell lineages through specific associations of transcription factors acting on distinct response elements along the same promoter. These data open new avenues regarding GnRH-mediated control of seasonal and circadian rhythms in reproductive physiology.

Differences in enhancer activity in mouse and zebrafish reporter assays are often associated with changes in gene expression

Background

Phenotypic evolution in animals is thought to be driven in large part by differences in gene expression patterns, which can result from sequence changes in cis-regulatory elements (cis-changes) or from changes in the expression pattern or function of transcription factors (trans-changes). While isolated examples of trans-changes have been identified, the scale of their overall contribution to regulatory and phenotypic evolution remains unclear.

Results

Here, we attempt to examine the prevalence of trans-effects and their potential impact on gene expression patterns in vertebrate evolution by comparing the function of identical human tissue-specific enhancer sequences in two highly divergent vertebrate model systems, mouse and zebrafish. Among 47 human conserved non-coding elements (CNEs) tested in transgenic mouse embryos and in stable zebrafish lines, at least one species-specific expression domain was observed in the majority (83%) of cases, and 36% presented dramatically different expression patterns between the two species. Although some of these discrepancies may be due to the use of different transgenesis systems in mouse and zebrafish, in some instances we found an association between differences in enhancer activity and changes in the endogenous gene expression patterns between mouse and zebrafish, suggesting a potential role for trans-changes in the evolution of gene expression.

Conclusions

In total, our results: (i) serve as a cautionary tale for studies investigating the role of human enhancers in different model organisms, and (ii) suggest that changes in the trans environment may play a significant role in the evolution of gene expression in vertebrates.

Improved predictions of transcription factor binding sites using physicochemical features of DNA

Typical approaches for predicting transcription factor binding sites (TFBSs) involve use of a position-specific weight matrix (PWM) to statistically characterize the sequences of the known sites. Recently, an alternative physicochemical approach, called SiteSleuth, was proposed. In this approach, a linear support vector machine (SVM) classifier is trained to distinguish TFBSs from background sequences based on local chemical and structural features of DNA. SiteSleuth appears to generally perform better than PWM-based methods. Here, we improve the SiteSleuth approach by considering both new physicochemical features and algorithmic modifications. New features are derived from Gibbs energies of amino acid-DNA interactions and hydroxyl radical cleavage profiles of DNA. Algorithmic modifications consist of inclusion of a feature selection step, use of a nonlinear kernel in the SVM classifier, and use of a consensus-based post-processing step for predictions. We also considered SVM classification based on letter features alone to distinguish performance gains from use of SVM-based models versus use of physicochemical features. The accuracy of each of the variant methods considered was assessed by cross validation using data available in the RegulonDB database for 54 Escherichia coli TFs, as well as by experimental validation using published ChIP-chip data available for Fis and Lrp.

PPARG binding landscapes in macrophages suggest a genome-wide contribution of PU.1 to divergent PPARG binding in human and mouse

Background

Genome-wide comparisons of transcription factor binding sites in different species can be used to evaluate evolutionary constraints that shape gene regulatory circuits and to understand how the interaction between transcription factors shapes their binding landscapes over evolution.

Results

We have compared the PPARG binding landscapes in macrophages to investigate the evolutionary impact on PPARG binding diversity in mouse and humans for this important nuclear receptor. Of note, only 5% of the PPARG binding sites were shared between the two species. In contrast, at the gene level, PPARG target genes conserved between both species constitute more than 30% of the target genes regulated by PPARG ligand in human macrophages. Moreover, the majority of all PPARG binding sites (55–60%) in macrophages show co-occupancy of the lineage-specification factor PU.1 in both species. Exploring the evolutionary dynamics of PPARG binding sites, we observed that PU.1 co-binding to PPARG sites appears to be important for possible PPARG ancestral functions such as lipid metabolism. Thus we speculate that PU.1 may have guided utilization of these species-specific PPARG conserved binding sites in macrophages during evolution.

Conclusions

We propose a model in which PU.1 sites may have served as “anchor” loci for the formation of new and functionally relevant PPARG binding sites throughout evolution. As PU.1 is an essential factor in macrophage biology, such an evolutionary mechanism would allow for the establishment of relevant PPARG regulatory modules in a PU.1-dependent manner and yet permit for nuanced regulatory changes in individual species.

Deciphering the transcriptional cis-regulatory code

Information about developmental gene expression resides in defined regulatory elements, called enhancers, in the non-coding part of the genome. Although cells reliably utilize enhancers to orchestrate gene expression, a cis-regulatory code that would allow their interpretation has remained one of the greatest challenges of modern biology. In this review, we summarize studies from the past three decades that describe progress towards revealing the properties of enhancers and discuss how recent approaches are providing unprecedented insights into regulatory elements in animal genomes. Over the next years, we believe that the functional characterization of regulatory sequences in entire genomes, combined with recent computational methods, will provide a comprehensive view of genomic regulatory elements and their building blocks and will enable researchers to begin to understand the sequence basis of the cis-regulatory code.

Reprogramming and development in nuclear transfer embryos and in interspecific systems

Nuclear transfer (NT) remains the most effective method to reprogram somatic cells to totipotency. Somatic cell nuclear transfer (SCNT) efficiency however remains low, but recurrent problems occurring in partially reprogrammed cloned embryos have recently been identified and some remedied. In particular, the trophectoderm has been identified as a lineage whose reprogramming success has a large influence on SCNT embryo development. Several interspecific hybrid and cybrid reprogramming systems have been developed as they offer various technical advantages and potential applications, and together with SCNT, they have led to the identification of a series of reprogramming events and responsible reprogramming factors. Interspecific incompatibilities hinder full exploitation of cross-species reprogramming systems, yet recent findings suggest that these may not constitute insurmountable obstacles.

Replication-timing boundaries facilitate cell-type and species-specific regulation of a rearranged human chromosome in mouse

In multicellular organisms, developmental changes to replication timing occur in 400-800 kb domains across half the genome. While examples of epigenetic control of replication timing have been described, a role for DNA sequence in mammalian replication-timing regulation has not been substantiated. To assess the role of DNA sequences in directing developmental changes to replication timing, we profiled replication timing in mice carrying a genetically rearranged Human Chromosome 21 (Hsa21). In two distinct mouse cell types, Hsa21 sequences maintained human-specific replication timing, except at points of Hsa21 rearrangement. Changes in replication timing at rearrangements extended up to 900 kb and consistently reconciled with the wild-type replication pattern at developmental boundaries of replication-timing domains. Our results are consistent with DNA sequence-driven regulation of Hsa21 replication timing during development and provide evidence that mammalian chromosomes consist of multiple independent units of replication-timing regulation.

Cistrome plasticity and mechanisms of cistrome reprogramming

Mammalian genomes contain thousands of cis-regulatory elements for each transcription factor (TF), but TFs only occupy a relatively small subset referred to as cistrome. Recent studies demonstrate that a TF cistrome might differ among different organisms, tissue types and individuals. In a cell, a TF cistrome might differ among different physiological states, pathological stages and between physiological and pathological conditions. It is, therefore, remarkable how highly plastic these binding profiles are, and how massively they can be reprogrammed in rapid response to intra/extracellular variations and during cell identity transitions and evolution. Biologically, cistrome reprogramming events tend to be followed by changes in transcriptional outputs, thus serving as transformative mechanisms to synchronically alter the biology of the cell. In this review, we discuss the molecular basis of cistrome plasticity and attempt to integrate the different mechanisms and biological conditions associated with cistrome reprogramming. Emerging data suggest that, when altered, these reprogramming events might be linked to tumor development and/or progression, which is a radical conceptual change in our mechanistic understanding of cancer and, potentially, other diseases.

Comparative epigenomic annotation of regulatory DNA

Despite the explosive growth of genomic data, functional annotation of regulatory sequences remains difficult. Here, we introduce "comparative epigenomics"-interspecies comparison of DNA and histone modifications-as an approach for annotation of the regulatory genome. We measured in human, mouse, and pig pluripotent stem cells the genomic distributions of cytosine methylation, H2A.Z, H3K4me1/2/3, H3K9me3, H3K27me3, H3K27ac, H3K36me3, transcribed RNAs, and P300, TAF1, OCT4, and NANOG binding. We observed that epigenomic conservation was strong in both rapidly evolving and slowly evolving DNA sequences, but not in neutrally evolving sequences. In contrast, evolutionary changes of the epigenome and the transcriptome exhibited a linear correlation. We suggest that the conserved colocalization of different epigenomic marks can be used to discover regulatory sequences. Indeed, seven pairs of epigenomic marks identified exhibited regulatory functions during differentiation of embryonic stem cells into mesendoderm cells. Thus, comparative epigenomics reveals regulatory features of the genome that cannot be discerned from sequence comparisons alone.

Pregnane-x-receptor controls hepatic glucuronidation during pregnancy and neonatal development in humanized UGT1 mice

In humanized UDP glucuronosyltransferase-1 (hUGT1) mice that express the entire UGT1 locus, the maternal hepatic UGT1A genes are dramatically induced 12-14 days after conception. Steroid induction of the UGT1A1 gene indicates that xenobiotic sensors, such as the pregnane X receptor (PXR) and constitutive androstane receptor (CAR), may underlie the induction process. In contrast, neonatal hUGT1 mice display severe hyperbilirubinemia, with limited expression of the UGT1A genes. This study identifies PXR as both a positive and negative regulator of the UGT1A1 gene. Pregnancy hormones, in particular the glucocorticoids, target PXR as a positive regulator of human glucuronidation. Employing reverse genetics, where PXR has been genetically deleted, hUGT1/Pxr(-/-) mice show limited induction of the liver UGT1A genes during pregnancy, whereas the exact opposite occurs in newborn mice. Neonatal hUGT1 mice show delayed expression of hepatic UGT1A1 and are severely hyperbilirubinemic. However, in hUGT1/Pxr(-/-) mice, hyperbilirubinemia is greatly reduced due to induction of hepatic UGT1A1. Thus, PXR serves to repress UGT1A1 gene expression during development. Transcriptional silencing of the UGT1A1 gene was relieved in neonatal hUGT1 hepatocytes through interruption of PXR by small interfering RNA.

CONCLUSION

PXR is a key regulator of pregnancy induced glucuronidation capacity in addition to modulating the severity of neonatal jaundice.

Cis- and trans-regulatory divergence between progenitor species determines gene-expression novelty in Arabidopsis allopolyploids

Gene-expression divergence between species shapes morphological evolution, but the molecular basis is largely unknown. Here we show cis- and trans-regulatory elements and chromatin modifications on gene-expression diversity in genetically tractable Arabidopsis allotetraploids. In Arabidopsis thaliana and Arabidopsis arenosa, both cis and trans with predominant cis-regulatory effects mediate gene-expression divergence. The majority of genes with both cis- and trans-effects are subjected to compensating interactions and stabilizing selection. Interestingly, cis- and trans-regulation is associated with chromatin modifications. In F1 allotetraploids, Arabidopsis arenosa trans factors predominately affect allelic expression divergence. Arabidopsis arenosa trans factors tend to upregulate Arabidopsis thaliana alleles, whereas Arabidopsis thaliana trans factors up- or down-regulate Arabidopsis arenosa alleles. In resynthesized and natural allotetraploids, trans effects drive expression of both homoeologous loci into the same direction. We provide evidence for natural selection and chromatin regulation in shaping gene-expression diversity during plant evolution and speciation.

Cross-species ChIP-seq studies provide insights into regulatory strategies of PPARγ in adipocytes

Three recent studies have investigated interspecies retention of binding sites of peroxisome proliferator-activated receptor γ (PPARγ), the master regulator of adipocyte differention, between mouse and human adipocytes. Here we discuss the major findings and demonstrate that retention of binding events is highly context-dependent.

Comparative studies of gene expression and the evolution of gene regulation

The hypothesis that differences in gene regulation have an important role in speciation and adaptation is more than 40 years old. With the advent of new sequencing technologies, we are able to characterize and study gene expression levels and associated regulatory mechanisms in a large number of individuals and species at an unprecedented resolution and scale. We have thus gained new insights into the evolutionary pressures that shape gene expression levels and have developed an appreciation for the relative importance of evolutionary changes in different regulatory genetic and epigenetic mechanisms. The current challenge is to link gene regulatory changes to adaptive evolution of complex phenotypes. Here we mainly focus on comparative studies in primates and how they are complemented by studies in model organisms.

Immediate-early gene activation by the MAPK pathways: what do and don't we know?

The study of IE (immediate-early) gene activation mechanisms has provided numerous paradigms for how transcription is controlled in response to extracellular signalling. Many of the findings have been derived from investigating one of the IE genes, FOS, and the models extrapolated to regulatory mechanisms for other IE genes. However, whereas the overall principles of activation appear similar, recent evidence suggests that the underlying mechanistic details may differ depending on cell type, cellular stimulus and IE gene under investigation. In the present paper, we review recent advances in our understanding of IE gene transcription, chiefly focusing on FOS and its activation by ERK (extracellular-signal-regulated kinase) MAPK (mitogen-activated protein kinase) pathway signalling. We highlight important fundamental regulatory principles, but also illustrate the gaps in our current knowledge and the potential danger in making assumptions based on extrapolation from disparate studies.

Allelic expression changes in Medaka (Oryzias latipes) hybrids between inbred strains derived from genetically distant populations

Variations in allele expressions between genetically distant populations are one of the most important factors which affects their morphological and physiological variations. These variations are caused by natural mutations accumulated in their habitats. It has been reported that allelic expression differences in the hybrids of genetically distant populations are different from parental strains. In that case, there is a possibility that allelic expression changes lead to novel phenotypes in hybrids. Based on genomic information of the genetically distant populations, quantification and comparison of allelic expression changes make importance of regulatory sequences (cis-acting factors) or upstream regulatory factors (trans-acting modulators) for these changes clearer. In this study, we focused on two Medaka inbred strains, Hd-rR and HNI, derived from genetically distant populations and their hybrids. They are highly polymorphic and we can utilize whole-genome information. To analyze allelic expression changes, we established a method to quantify and compare allele-specific expressions of 11 genes between the parental strains and their reciprocal hybrids. In intestines of reciprocal hybrids, allelic expression was either similar or different in comparison with the parental strains. Total expressions in Hd-rR and HNI were tissue-dependent in the case of HPRT1, with high up-regulation of Hd-rR allele expression in liver. The proportion of genes with differential allelic expression in Medaka hybrids seems to be the same as that in other animals, despite the high SNP rate in the genomes of the two inbred strains. It is suggested that each tissue of the strain difference in trans-acting modulators is more important than polymorphisms in cis-regulatory sequences in producing the allelic expression changes in reciprocal hybrids.

A selection index for gene expression evolution and its application to the divergence between humans and chimpanzees

The importance of gene regulation in animal evolution is a matter of long-standing interest, but measuring the impact of selection on gene expression has proven a challenge. Here, we propose a selection index of gene expression as a straightforward method for assessing the mode and strength of selection operating on gene expression levels. The index is based on the widely used McDonald-Kreitman test and requires the estimation of four quantities: the within-species and between-species expression variances as well as the sequence heterozygosity and divergence of neutrally evolving sequences. We apply the method to data from human and chimpanzee lymphoblastoid cell lines and show that gene expression is in general under strong stabilizing selection. We also demonstrate how the same framework can be used to estimate the proportion of adaptive gene expression evolution.

Widespread site-dependent buffering of human regulatory polymorphism

The average individual is expected to harbor thousands of variants within non-coding genomic regions involved in gene regulation. However, it is currently not possible to interpret reliably the functional consequences of genetic variation within any given transcription factor recognition sequence. To address this, we comprehensively analyzed heritable genome-wide binding patterns of a major sequence-specific regulator (CTCF) in relation to genetic variability in binding site sequences across a multi-generational pedigree. We localized and quantified CTCF occupancy by ChIP-seq in 12 related and unrelated individuals spanning three generations, followed by comprehensive targeted resequencing of the entire CTCF–binding landscape across all individuals. We identified hundreds of variants with reproducible quantitative effects on CTCF occupancy (both positive and negative). While these effects paralleled protein–DNA recognition energetics when averaged, they were extensively buffered by striking local context dependencies. In the significant majority of cases buffering was complete, resulting in silent variants spanning every position within the DNA recognition interface irrespective of level of binding energy or evolutionary constraint. The prevalence of complex partial or complete buffering effects severely constrained the ability to predict reliably the impact of variation within any given binding site instance. Surprisingly, 40% of variants that increased CTCF occupancy occurred at positions of human–chimp divergence, challenging the expectation that the vast majority of functional regulatory variants should be deleterious. Our results suggest that, even in the presence of “perfect” genetic information afforded by resequencing and parallel studies in multiple related individuals, genomic site-specific prediction of the consequences of individual variation in regulatory DNA will require systematic coupling with empirical functional genomic measurements.

A comprehensive understanding of the contribution of individual genome sequences to disease and quantitative traits will require the general ability to predict consequences of genetic variation in non-protein-coding regions, particularly those involved in gene regulation. Here we tested the power to predict such consequences when presented with “complete” information encompassing the genomic DNA binding site patterns of a well-studied regulatory protein across multiple related individuals, coupled with all individual genome sequences at the binding positions. We find that, while there is reasonable ability to predict the average effects of variation within the consensus recognition sequence of a transcriptional regulator, it is not possible to determine reliably the consequences of variation at any given genomic instance. This suggests that the interpretation of individual genome sequences will require comprehensive complementation with functional genomic studies.

Modelling human regulatory variation in mouse: finding the function in genome-wide association studies and whole-genome sequencing

An increasing body of literature from genome-wide association studies and human whole-genome sequencing highlights the identification of large numbers of candidate regulatory variants of potential therapeutic interest in numerous diseases. Our relatively poor understanding of the functions of non-coding genomic sequence, and the slow and laborious process of experimental validation of the functional significance of human regulatory variants, limits our ability to fully benefit from this information in our efforts to comprehend human disease. Humanized mouse models (HuMMs), in which human genes are introduced into the mouse, suggest an approach to this problem. In the past, HuMMs have been used successfully to study human disease variants; e.g., the complex genetic condition arising from Down syndrome, common monogenic disorders such as Huntington disease and β-thalassemia, and cancer susceptibility genes such as BRCA1. In this commentary, we highlight a novel method for high-throughput single-copy site-specific generation of HuMMs entitled High-throughput Human Genes on the X Chromosome (HuGX). This method can be applied to most human genes for which a bacterial artificial chromosome (BAC) construct can be derived and a mouse-null allele exists. This strategy comprises (1) the use of recombineering technology to create a human variant-harbouring BAC, (2) knock-in of this BAC into the mouse genome using Hprt docking technology, and (3) allele comparison by interspecies complementation. We demonstrate the throughput of the HuGX method by generating a series of seven different alleles for the human NR2E1 gene at Hprt. In future challenges, we consider the current limitations of experimental approaches and call for a concerted effort by the genetics community, for both human and mouse, to solve the challenge of the functional analysis of human regulatory variation.

Intra-nuclear mobility and target search mechanisms of transcription factors: a single-molecule perspective on gene expression

Precise expression of specific genes in time and space is at the basis of cellular viability as well as correct development of organisms. Understanding the mechanisms of gene regulation is fundamental and still one of the great challenges for biology. Gene expression is regulated also by specific transcription factors that recognize and bind to specific DNA sequences. Transcription factors dynamics, and especially the way they sample the nucleoplasmic space during the search for their specific target in the genome, are a key aspect for regulation and it has been puzzling researchers for forty years. The scope of this review is to give a state-of-the-art perspective over the intra-nuclear mobility and the target search mechanisms of specific transcription factors at the molecular level. Going through the seminal biochemical experiments that have raised the first questions about target localization and the theoretical grounds concerning target search processes, we describe the most recent experimental achievements and current challenges in understanding transcription factors dynamics and interactions with DNA using in vitro assays as well as in live prokaryotic and eukaryotic cells. This article is part of a Special Issue entitled: Nuclear Transport and RNA Processing.

Estimation of pairwise sequence similarity of mammalian enhancers with word neighbourhood counts

MOTIVATION

The identity of cells and tissues is to a large degree governed by transcriptional regulation. A major part is accomplished by the combinatorial binding of transcription factors at regulatory sequences, such as enhancers. Even though binding of transcription factors is sequence-specific, estimating the sequence similarity of two functionally similar enhancers is very difficult. However, a similarity measure for regulatory sequences is crucial to detect and understand functional similarities between two enhancers and will facilitate large-scale analyses like clustering, prediction and classification of genome-wide datasets.

RESULTS

We present the standardized alignment-free sequence similarity measure N2, a flexible framework that is defined for word neighbourhoods. We explore the usefulness of adding reverse complement words as well as words including mismatches into the neighbourhood. On simulated enhancer sequences as well as functional enhancers in mouse development, N2 is shown to outperform previous alignment-free measures. N2 is flexible, faster than competing methods and less susceptible to single sequence noise and the occurrence of repetitive sequences. Experiments on the mouse enhancers reveal that enhancers active in different tissues can be separated by pairwise comparison using N2.

CONCLUSION

N2 represents an improvement over previous alignment-free similarity measures without compromising speed, which makes it a good candidate for large-scale sequence comparison of regulatory sequences.

AVAILABILITY

The software is part of the open-source C++ library SeqAn (www.seqan.de) and a compiled version can be downloaded at http://www.seqan.de/projects/alf.html.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Waves of retrotransposon expansion remodel genome organization and CTCF binding in multiple mammalian lineages

CTCF-binding locations represent regulatory sequences that are highly constrained over the course of evolution. To gain insight into how these DNA elements are conserved and spread through the genome, we defined the full spectrum of CTCF-binding sites, including a 33/34-mer motif, and identified over five thousand highly conserved, robust, and tissue-independent CTCF-binding locations by comparing ChIP-seq data from six mammals. Our data indicate that activation of retroelements has produced species-specific expansions of CTCF binding in rodents, dogs, and opossum, which often functionally serve as chromatin and transcriptional insulators. We discovered fossilized repeat elements flanking deeply conserved CTCF-binding regions, indicating that similar retrotransposon expansions occurred hundreds of millions of years ago. Repeat-driven dispersal of CTCF binding is a fundamental, ancient, and still highly active mechanism of genome evolution in mammalian lineages.

PaperClip

Mechanisms underlying the tissue-specific and regulated activity of the Gnrhr promoter in mammals

The GnRH receptor (GnRHR) plays a central role in the development and maintenance of reproductive function in mammals. Following stimulation by GnRH originating from the hypothalamus, GnRHR triggers multiple signaling events that ultimately stimulate the synthesis and the periodic release of the gonadotropins, luteinizing-stimulating hormone (LH) and follicle-stimulating hormones (FSH) which, in turn, regulate gonadal functions including steroidogenesis and gametogenesis. The concentration of GnRHR at the cell surface is essential for the amplitude and the specificity of gonadotrope responsiveness. The number of GnRHR is submitted to strong regulatory control during pituitary development, estrous cycle, pregnancy, lactation, or after gonadectomy. These modulations take place, at least in part, at the transcriptional level. To analyze this facet of the reproductive function, the 5' regulatory sequences of the gene encoding the GnRHR have been isolated and characterized through in vitro and in vivo approaches. This review summarizes results obtained with the mouse, rat, human, and ovine promoters either by transient transfection assays or by means of transgenic mice.

A functional selection model explains evolutionary robustness despite plasticity in regulatory networks

Evolutionary rewiring of regulatory networks is an important source of diversity among species. Previous evidence suggested substantial divergence of regulatory networks across species. However, systematically assessing the extent of this plasticity and its functional implications has been challenging due to limited experimental data and the noisy nature of computational predictions. Here, we introduce a novel approach to study cis-regulatory evolution, and use it to trace the regulatory history of 88 DNA motifs of transcription factors across 23 Ascomycota fungi. While motifs are conserved, we find a pervasive gain and loss in the regulation of their target genes. Despite this turnover, the biological processes associated with a motif are generally conserved. We explain these trends using a model with a strong selection to conserve the overall function of a transcription factor, and a much weaker selection over the specific genes it targets. The model also accounts for the turnover of bound targets measured experimentally across species in yeasts and mammals. Thus, selective pressures on regulatory networks mostly tolerate local rewiring, and may allow for subtle fine-tuning of gene regulation during evolution.

Using genomic tools to study regulatory evolution

Differences in gene regulation are thought to play an important role in speciation and adaptation. Comparative genomic studies of gene expression levels have identified a large number of differentially expressed genes among species, and, in a number of cases, also pointed to connections between interspecies differences in gene regulation and differences in ultimate physiological or morphological phenotypes. The mechanisms underlying changes in gene regulation are also being actively studied using comparative genomic approaches. However, the relative importance of different regulatory mechanisms to interspecies differences in gene expression levels is not yet well understood. In particular, it is often difficult to infer causality between apparent differences in regulatory mechanisms and changes in gene expression levels, a challenge that is compounded by the fact that the link between sequence variation and gene regulation is not clear. Indeed, in certain cases, gene regulation can be conserved even when sequences at associated regulatory elements have changed. In this chapter, I examine different genomic approaches to the study of regulatory evolution and the underlying genetic and epigenetic regulatory mechanisms. I try to distinguish between hypothesis-driven and exploratory studies, and argue that the latter class of studies provides valuable information in its own right as well as necessary context for the former. I discuss issues related to study designs and statistical analyses of genomic studies, and review the evidence for natural selection on gene expression levels and associated regulatory mechanisms. Most of the issues that are discussed pertain to the general nature of multivariate genomic data, and thus are often relevant regardless of the technology that is used to collect high-throughput genomic data (for example, microarrays or massively parallel sequencing).

Tumour necrosis factor gene polymorphism and disease prevalence

Tumour necrosis factor (TNF), an important proinflammatory cytokine, plays a role in the regulation of cell differentiation, proliferation and death, as well as in inflammation, innate and adaptive immune responses, and also implicated in a wide variety of human diseases. The presence of DNA sequence variations in regulatory region might interfere with transcription of TNF gene, influencing the circulating level of TNF and thus increases the susceptibility to human diseases (infectious, cancer, autoimmune, neurodegenerative and other diseases). In this review, we have comprehensively analysed various published case-control studies of different types of human diseases, in which TNF gene polymorphism played a role, and computationally predicted several single nucleotide polymorphisms (SNPs) lie in transcription factor-binding sites (TFBS) of transcription factors (TFs). It has been observed that TNF enhancer polymorphism is implicated in several diseases, and TNF rs1800629 and rs361525 SNPs are the most important in human disease susceptibility as these might influence the transcription of TNF gene. Thirty-two SNPs lies in TFBS of 20 TFs have been detected in the TNF upstream region. It has been found that TNF enhancer polymorphism influences the serum level of TNF in different human diseases and thus affects the susceptibility to diseases. The presence of DNA sequence variation in TNF gene causes the modification of transcriptional regulation and thus responsible for association of susceptibility/resistance with human diseases.

Genomically humanized mice: technologies and promises

Mouse models have become an invaluable tool for understanding human health and disease owing to our ability to manipulate the mouse genome exquisitely. Recent progress in genomic analysis has led to an increase in the number and type of disease-causing mutations detected and has also highlighted the importance of non-coding regions. As a result, there is increasing interest in creating 'genomically' humanized mouse models, in which entire human genomic loci are transferred into the mouse genome. The technical challenges towards achieving this aim are large but are starting to be tackled with success.

Deficient induction response in a Xenopus nucleocytoplasmic hybrid

Defects in induction signaling and response underlie the nucleocytoplasmic incompatibility between two evolutionarily distant frog species, while specific treatments partially restore this response in explants and whole embryos.

Incompatibilities between the nucleus and the cytoplasm of sufficiently distant species result in developmental arrest of hybrid and nucleocytoplasmic hybrid (cybrid) embryos. Several hypotheses have been proposed to explain their lethality, including problems in embryonic genome activation (EGA) and/or nucleo-mitochondrial interactions. However, conclusive identification of the causes underlying developmental defects of cybrid embryos is still lacking. We show here that while over 80% of both Xenopus laevis and Xenopus (Silurana) tropicalis same-species androgenetic haploids develop to the swimming tadpole stage, the androgenetic cybrids formed by the combination of X. laevis egg cytoplasm and X. tropicalis sperm nucleus invariably fail to gastrulate properly and never reach the swimming tadpole stage. In spite of this arrest, these cybrids show quantitatively normal EGA and energy levels at the stage where their initial gastrulation defects are manifested. The nucleocytoplasmic incompatibility between these two species instead results from a combination of factors, including a reduced emission of induction signal from the vegetal half, a decreased sensitivity of animal cells to induction signals, and differences in a key embryonic protein (Xbra) concentration between the two species, together leading to inefficient induction and defective convergence-extension during gastrulation. Indeed, increased exposure to induction signals and/or Xbra signalling partially rescues the induction response in animal explants and whole cybrid embryos. Altogether, our study demonstrates that the egg cytoplasm of one species may not support the development promoted by the nucleus of another species, even if this nucleus does not interfere with the cytoplasmic/maternal functions of the egg, while the egg cytoplasm is also capable of activating the genome of that nucleus. Instead, our results provide evidence that inefficient signalling and differences in the concentrations of key proteins between species lead to developmental defects in cybrids. Finally, they show that the incompatibilities of cybrids can be corrected by appropriate treatments.

When two species evolve separately for several million years, their respective genomes accumulate many small changes that together are responsible for the differences in their characters. Some of these affect the way eggs are prepared inside the germline, and/or how embryos develop, such that the egg cytoplasm of a given species can only support development promoted by its own genome or nucleus. Thus, developmental incompatibility arises between the cytoplasm and the nucleus of distant species during evolution and we don't know its mechanism. We have studied this phenomenon in an advantageous system using two evolutionarily distant frog species (Xenopus laevis and Xenopus tropicalis). We found that hybrid frog embryos with X. laevis cytoplasm and X. tropicalis nuclei are always defective in an important process that is necessary to generate morphogenetic cell movements during development. Through a series of experiments in which we dissect out and/or recombine parts of such hybrid embryos and observe their behaviour in culture, we show that this phenomenon occurs because of malfunctions in the signalling cascade that is responsible for generating these cell movements. Thus, we postulate that inefficient molecular signalling contributes to the death of such hybrids.

An interspecies analysis reveals a key role for unmethylated CpG dinucleotides in vertebrate Polycomb complex recruitment

The role of DNA sequence in determining chromatin state is incompletely understood. We have previously demonstrated that large chromosomal segments from human cells recapitulate their native chromatin state in mouse cells, but the relative contribution of local sequences versus their genomic context remains unknown. In this study, we compare orthologous chromosomal regions for which the human locus establishes prominent sites of Polycomb complex recruitment in pluripotent stem cells, whereas the corresponding mouse locus does not. Using recombination-mediated cassette exchange at the mouse locus, we establish the primacy of local sequences in the encoding of chromatin state. We show that the signal for chromatin bivalency is redundantly encoded across a bivalent domain and that this reflects competition between Polycomb complex recruitment and transcriptional activation. Furthermore, our results suggest that a high density of unmethylated CpG dinucleotides is sufficient for vertebrate Polycomb recruitment. This model is supported by analysis of DNA methyltransferase-deficient embryonic stem cells.

Rapid turnover of functional sequence in human and other genomes

The amount of a genome's sequence that is functional has been surprisingly difficult to estimate accurately. This has severely hindered analyses asking whether the amount of functional genomic sequence correlates with organismal complexity. Most studies estimate these amounts by considering nucleotide substitution rates within aligned sequences. These approaches show reduced power to identify sequence that is aligned, functional, and constrained only within narrowly defined phyla. The neutral indel model exploits insertions or deletions (indels) rather than substitutions in predicting functional sequence. Surprisingly, this method indicates that half of all functional sequence is specific to individual eutherian lineages. This review considers the rates at which coding or noncoding and functional or nonfunctional sequence changes among mammalian genomes. In contrast to the slow rate at which protein-coding sequence changes, functional noncoding sequence appears to change or be turned over at rapid rates in mammals.

Inferring regulatory mechanisms from patterns of evolutionary divergence

The number of sequenced species is increasing at a staggering rate, calling for new approaches for incorporating evolutionary information in the study of biological mechanisms. Evolutionary conservation is widely used for assigning a function to new proteins and for predicting functional coding or non-coding sequences. Here, we argue for a complementary approach that focuses on the divergence of regulatory programs. Regulatory mechanisms can be learned from patterns of evolutionary divergence in regulatory properties such as gene expression, transcription factor binding or nucleosome positioning. We review examples of this concept using yeast as a model system, and highlight a hybrid-based approach that is highly instrumental in this analysis.

Discriminative prediction of mammalian enhancers from DNA sequence

Accurately predicting regulatory sequences and enhancers in entire genomes is an important but difficult problem, especially in large vertebrate genomes. With the advent of ChIP-seq technology, experimental detection of genome-wide EP300/CREBBP bound regions provides a powerful platform to develop predictive tools for regulatory sequences and to study their sequence properties. Here, we develop a support vector machine (SVM) framework which can accurately identify EP300-bound enhancers using only genomic sequence and an unbiased set of general sequence features. Moreover, we find that the predictive sequence features identified by the SVM classifier reveal biologically relevant sequence elements enriched in the enhancers, but we also identify other features that are significantly depleted in enhancers. The predictive sequence features are evolutionarily conserved and spatially clustered, providing further support of their functional significance. Although our SVM is trained on experimental data, we also predict novel enhancers and show that these putative enhancers are significantly enriched in both ChIP-seq signal and DNase I hypersensitivity signal in the mouse brain and are located near relevant genes. Finally, we present results of comparisons between other EP300/CREBBP data sets using our SVM and uncover sequence elements enriched and/or depleted in the different classes of enhancers. Many of these sequence features play a role in specifying tissue-specific or developmental-stage-specific enhancer activity, but our results indicate that some features operate in a general or tissue-independent manner. In addition to providing a high confidence list of enhancer targets for subsequent experimental investigation, these results contribute to our understanding of the general sequence structure of vertebrate enhancers.

Transcription factor interactions in genomic nuclear receptor function

Transcription factors (TF) regulate gene expression acting as DNA sequence-specific binding factors, orchestrating cofactor recruitment and assembly of the transcriptional machinery. Nuclear receptors, a ligand-inducible TF class, regulate a large proportion of the genome, yet achieve highly cell-specific and context-dependent transcription, despite their widespread expression. High-throughput genome-wide profiling of TF binding reveals a startling proportion of colocalized cell- and context-specific TF-binding patterns, implying TF interactions play a critical role in transcription. These interactions depend on the chromatin architecture, that predominantly acts to predetermine accessibility of TF-binding sites at regulatory elements. Here, we summarize recent findings that highlight the importance of combinatorial TF interactions in determining diverse nuclear receptor-mediated transcriptional responses, emphasizing the significance of chromatin structure in directing TF and nuclear receptor recruitment. Interactions between TFs are likely to be a general mechanism of regulatory factors, contributing to transcriptional control in health and disease.

Genome-wide integration on transcription factors, histone acetylation and gene expression reveals genes co-regulated by histone modification patterns

N-terminal tails of H2A, H2B, H3 and H4 histone families are subjected to posttranslational modifications that take part in transcriptional regulation mechanisms, such as transcription factor binding and gene expression. Regulation mechanisms under control of histone modification are important but remain largely unclear, despite of emerging datasets for comprehensive analysis of histone modification. In this paper, we focus on what we call genetic harmonious units (GHUs), which are co-occurring patterns among transcription factor binding, gene expression and histone modification. We present the first genome-wide approach that captures GHUs by combining ChIP-chip with microarray datasets from Saccharomyces cerevisiae. Our approach employs noise-robust soft clustering to select patterns which share the same preferences in transcription factor-binding, histone modification and gene expression, which are all currently implied to be closely correlated. The detected patterns are a well-studied acetylation of lysine 16 of H4 in glucose depletion as well as co-acetylation of five lysine residues of H3 with H4 Lys12 and H2A Lys7 responsible for ribosome biogenesis. Furthermore, our method further suggested the recognition of acetylated H4 Lys16 being crucial to histone acetyltransferase ESA1, whose essential role is still under controversy, from a microarray dataset on ESA1 and its bypass suppressor mutants. These results demonstrate that our approach allows us to provide clearer principles behind gene regulation mechanisms under histone modifications and detect GHUs further by applying to other microarray and ChIP-chip datasets. The source code of our method, which was implemented in MATLAB (http://www.mathworks.com/), is available from the supporting page for this paper: http://www.bic.kyoto-u.ac.jp/pathway/natsume/hm_detector.htm.

Genome-wide identification of conserved regulatory function in diverged sequences

Plasticity of gene regulatory encryption can permit DNA sequence divergence without loss of function. Functional information is preserved through conservation of the composition of transcription factor binding sites (TFBS) in a regulatory element. We have developed a method that can accurately identify pairs of functional noncoding orthologs at evolutionarily diverged loci by searching for conserved TFBS arrangements. With an estimated 5% false-positive rate (FPR) in approximately 3000 human and zebrafish syntenic loci, we detected approximately 300 pairs of diverged elements that are likely to share common ancestry and have similar regulatory activity. By analyzing a pool of experimentally validated human enhancers, we demonstrated that 7/8 (88%) of their predicted functional orthologs retained in vivo regulatory control. Moreover, in 5/7 (71%) of assayed enhancer pairs, we observed concordant expression patterns. We argue that TFBS composition is often necessary to retain and sufficient to predict regulatory function in the absence of overt sequence conservation, revealing an entire class of functionally conserved, evolutionarily diverged regulatory elements that we term "covert."

The androgen receptor fuels prostate cancer by regulating central metabolism and biosynthesis

The androgen receptor (AR) is a key regulator of prostate growth and the principal drug target for the treatment of prostate cancer. Previous studies have mapped AR targets and identified some candidates which may contribute to cancer progression, but did not characterize AR biology in an integrated manner. In this study, we took an interdisciplinary approach, integrating detailed genomic studies with metabolomic profiling and identify an anabolic transcriptional network involving AR as the core regulator. Restricting flux through anabolic pathways is an attractive approach to deprive tumours of the building blocks needed to sustain tumour growth. Therefore, we searched for targets of the AR that may contribute to these anabolic processes and could be amenable to therapeutic intervention by virtue of differential expression in prostate tumours. This highlighted calcium/calmodulin-dependent protein kinase kinase 2, which we show is overexpressed in prostate cancer and regulates cancer cell growth via its unexpected role as a hormone-dependent modulator of anabolic metabolism. In conclusion, it is possible to progress from transcriptional studies to a promising therapeutic target by taking an unbiased interdisciplinary approach.

Identification of cis-elements and transcription factors regulating neuronal activity-dependent transcription of human BDNF gene

Brain-derived neurotrophic factor (BDNF) is an important mediator of activity-dependent functions of the nervous system and its expression is dysregulated in several neuropsychiatric disorders. Regulation of rodent BDNF neuronal activity-dependent transcription has been relatively well characterized. Here, we have studied regulation of human BDNF (hBDNF) transcription by membrane depolarization of cultured mouse or rat primary cortical neurons expressing hBDNF gene or transfected with hBDNF promoter constructs, respectively. We identified an asymmetric E-box-like element, PasRE [basic helix-loop-helix (bHLH)-PAS transcription factor response element], in hBDNF promoter I and demonstrate that binding of this element by bHLH-PAS transcription factors ARNT2 (aryl hydrocarbon receptor nuclear translocator 2) and NPAS4 (neuronal PAS domain protein 4) is crucial for neuronal activity-dependent transcription from promoter I. We show that binding of CREB (cAMP response element-binding protein) to the cAMP/Ca(2+)-response element (CRE) in hBDNF promoter IV is critical for activity-dependent transcription from this promoter and that upstream stimulatory factor (USF) transcription factors also contribute to the activation by binding to the upstream stimulatory factor binding element (UBE) in hBDNF promoter IV. However, we report that full induction of hBDNF exon IV mRNA transcription is dependent on ARNT2 and NPAS4 binding to a PasRE in promoter IV. Finally, we demonstrate that CRE and PasRE elements in hBDNF promoter IX are required for the induction of this promoter by neuronal activity. Together, the results of this study have identified the cis-elements and transcription factors regulating neuronal activity-dependent transcription of human BDNF gene.

Genome-wide approaches to the study of adaptive gene expression evolution: systematic studies of evolutionary adaptations involving gene expression will allow many fundamental questions in evolutionary biology to be addressed

The role of gene expression in evolutionary adaptation has been a subject of debate for over 40 years. cis-regulation of transcription has been proposed to be the primary source of morphological novelty in evolution, though this is based on only a handful of examples. Recently the first genome-wide studies of gene expression adaptation have been published, giving us an initial global view of this process. Systematic studies such as these will allow a number of key questions currently facing the field of gene expression evolution to be addressed.

Genetic architecture of regulatory variation in Arabidopsis thaliana

Studying the genetic regulation of expression variation is a key method to dissect complex phenotypic traits. To examine the genetic architecture of regulatory variation in Arabidopsis thaliana, we performed genome-wide association (GWA) mapping of gene expression in an F(1) hybrid diversity panel. At a genome-wide false discovery rate (FDR) of 0.2, an associated single nucleotide polymorphism (SNP) explains >38% of trait variation. In comparison with SNPs that are distant from the genes to which they were associated, locally associated SNPs are preferentially found in regions with extended linkage disequilibrium (LD) and have distinct population frequencies of the derived alleles (where Arabidopsis lyrata has the ancestral allele), suggesting that different selective forces are acting. Locally associated SNPs tend to have additive inheritance, whereas distantly associated SNPs are primarily dominant. In contrast to results from mapping of expression quantitative trait loci (eQTL) in linkage studies, we observe extensive allelic heterogeneity for local regulatory loci in our diversity panel. By association mapping of allele-specific expression (ASE), we detect a significant enrichment for cis-acting variation in local regulatory variation. In addition to gene expression variation, association mapping of splicing variation reveals both local and distant genetic regulation for intron and exon level traits. Finally, we identify candidate genes for 59 diverse phenotypic traits that were mapped to eQTL.

Cross species comparison of C/EBPα and PPARγ profiles in mouse and human adipocytes reveals interdependent retention of binding sites

Background

The transcription factors peroxisome proliferator activated receptor γ (PPARγ) and CCAAT/enhancer binding protein α (C/EBPα) are key transcriptional regulators of adipocyte differentiation and function. We and others have previously shown that binding sites of these two transcription factors show a high degree of overlap and are associated with the majority of genes upregulated during differentiation of murine 3T3-L1 adipocytes.

Results

Here we have mapped all binding sites of C/EBPα and PPARγ in human SGBS adipocytes and compared these with the genome-wide profiles from mouse adipocytes to systematically investigate what biological features correlate with retention of sites in orthologous regions between mouse and human. Despite a limited interspecies retention of binding sites, several biological features make sites more likely to be retained. First, co-binding of PPARγ and C/EBPα in mouse is the most powerful predictor of retention of the corresponding binding sites in human. Second, vicinity to genes highly upregulated during adipogenesis significantly increases retention. Third, the presence of C/EBPα consensus sites correlate with retention of both factors, indicating that C/EBPα facilitates recruitment of PPARγ. Fourth, retention correlates with overall sequence conservation within the binding regions independent of C/EBPα and PPARγ sequence patterns, indicating that other transcription factors work cooperatively with these two key transcription factors.

Conclusions

This study provides a comprehensive and systematic analysis of what biological features impact on retention of binding sites between human and mouse. Specifically, we show that the binding of C/EBPα and PPARγ in adipocytes have evolved in a highly interdependent manner, indicating a significant cooperativity between these two transcription factors.

Functional and mechanistic diversity of distal transcription enhancers

Biological differences among metazoans and between cell types in a given organism arise in large part due to differences in gene expression patterns. Gene-distal enhancers are key contributors to these expression patterns, exhibiting both sequence diversity and cell type specificity. Studies of long-range interactions indicate that enhancers are often important determinants of nuclear organization, contributing to a general model for enhancer function that involves direct enhancer-promoter contact. However, mechanisms for enhancer function are emerging that do not fit solely within such a model, suggesting that enhancers as a class of DNA regulatory element may be functionally and mechanistically diverse.

Gene expression differences among primates are associated with changes in a histone epigenetic modification

Changes in gene regulation are thought to play an important role in speciation and adaptation, especially in primates. However, we still know relatively little about the mechanisms underlying regulatory evolution. In particular, the extent to which epigenetic modifications underlie gene expression differences between primates is not yet known. Our study focuses on an epigenetic histone modification, H3K4me3, which is thought to promote transcription. To investigate the contribution of H3K4me3 to regulatory differences between species, we collected gene expression data and identified H3K4me3-associated genomic regions in lymphoblastoid cell lines (LCLs) from humans, chimpanzees, and rhesus macaques, using three cell lines from each species. We found strong evidence for conservation of H3K4me3 localization in primates. Moreover, regardless of species, H3K4me3 is consistently enriched near annotated transcription start sites (TSS), and highly expressed genes are more likely than lowly expressed genes to have the histone modification near their TSS. Interestingly, we observed an enrichment of interspecies differences in H3K4me3 at the TSS of genes that are differentially expressed between species. We estimate that as much as 7% of gene expression differences between the LCLs of humans, chimpanzees, and rhesus macaques may be explained, at least in part, by changes in the status of H3K4me3 histone modifications. Our results suggest a modest, yet important role for epigenetic changes in gene expression differences between primates.

GnRH receptor gene expression in the developing rat hippocampus: transcriptional regulation and potential roles in neuronal plasticity

In the pituitary of mammals, the GnRH receptor (GnRHR) plays a primary role in the control of reproductive function. It is further expressed in the hippocampus, where its function, however, is not well defined. By quantitative RT-PCR analyses, we demonstrate herein that the onset of GnRHR gene (Gnrhr) expression in the rat hippocampus was unexpectedly delayed as compared to the pituitary and only occurred after birth. Using a previously described transgenic mouse model bearing the human placental alkaline phosphatase reporter gene under the control of the rat Gnrhr promoter, we established a positive correlation between the temporal pattern of Gnrhr mRNA levels and promoter activity in the hippocampal formation. The gradual appearance of human placental alkaline phosphatase transgene expression occurred simultaneously in the hippocampus and interconnected structures such as the lateral septum and the amygdala, coinciding with the establishment of hippocampo-septal projections. Analysis of transcription factors together with transient transfection assays in hippocampal neurons indicated that the combinatorial code governing the hippocampus-specific expression of the Gnrhr is distinct from the pituitary, likely involving transactivating factors such as NUR77, cyclic AMP response element binding protein, and Finkel-Biskis-Jinkins murine osteosarcoma virus oncogene homolog. A silencing transcription factor acting via the -3255/-1135 promoter region of the Gnrhr may be responsible for the transcriptional repression observed around birth. Finally, GnRH directly stimulated via activation of its receptor the expression of several marker genes of neuronal plasticity such as Egr1, synaptophysin, and spinophilin in hippocampal primary cultures, suggesting a role for GnRHR in neuronal plasticity. Further characterization of these mechanisms may help unravel important functions of GnRH/GnRHR signaling in the brain.