Identification of potential regulatory motifs in odorant receptor genes by analysis of promoter sequences

Establishment and maintenance of random monoallelic expression

This Review elucidates the regulatory principles of random monoallelic expression by focusing on two well-studied examples: the X-chromosome inactivation regulator Xist and the olfactory receptor gene family. Although the choice of a single X chromosome or olfactory receptor occurs in different developmental contexts, common gene regulatory principles guide monoallelic expression in both systems. In both cases, an event breaks the symmetry between genetically and epigenetically identical copies of the gene, leading to the expression of one single random allele, stabilized through negative feedback control. Although many regulatory steps that govern the establishment and maintenance of monoallelic expression have been identified, key pieces of the puzzle are still missing. We provide an overview of the current knowledge and models for the monoallelic expression of Xist and olfactory receptors. We discuss their similarities and differences, and highlight open questions and approaches that could guide the study of other monoallelically expressed genes.

Epigenetic programming of stochastic olfactory receptor choice

The mammalian sense of smell relies upon a vast array of receptor proteins to detect odorant compounds present in the environment. The proper deployment of these receptor proteins in olfactory sensory neurons is orchestrated by a suite of epigenetic processes that remodel the olfactory genes in differentiating neuronal progenitors. The goal of this review is to elucidate the central role of gene regulatory processes acting in neuronal progenitors of olfactory sensory neurons that lead to a singular expression of an odorant receptor in mature olfactory sensory neurons. We begin by describing the principal features of odorant receptor gene expression in mature olfactory sensory neurons. Next, we delineate our current understanding of how these features emerge from multiple gene regulatory mechanisms acting in neuronal progenitors. Finally, we close by discussing the key gaps in our understanding of how these regulatory mechanisms work and how they interact with each other over the course of differentiation.

Emergence and influence of sequence bias in evolutionarily malleable, mammalian tandem arrays

BACKGROUND

The radiation of mammals at the extinction of the dinosaurs produced a plethora of new forms-as diverse as bats, dolphins, and elephants-in only 10-20 million years. Behind the scenes, adaptation to new niches is accompanied by extensive innovation in large families of genes that allow animals to contact the environment, including chemosensors, xenobiotic enzymes, and immune and barrier proteins. Genes in these "outward-looking" families are allelically diverse among humans and exhibit tissue-specific and sometimes stochastic expression.

RESULTS

Here, we show that these tandem arrays of outward-looking genes occupy AT-biased isochores and comprise the "tissue-specific" gene class that lack CpG islands in their promoters. Models of mammalian genome evolution have not incorporated the sharply different functions and transcriptional patterns of genes in AT- versus GC-biased regions. To examine the relationship between gene family expansion, sequence content, and allelic diversity, we use population genetic data and comparative analysis. First, we find that AT bias can emerge during evolutionary expansion of gene families in cis. Second, human genes in AT-biased isochores or with GC-poor promoters experience relatively low rates of de novo point mutation today but are enriched for non-synonymous variants. Finally, we find that isochores containing gene clusters exhibit low rates of recombination.

CONCLUSIONS

Our analyses suggest that tolerance of non-synonymous variation and low recombination are two forces that have produced the depletion of GC bases in outward-facing gene arrays. In turn, high AT content exerts a profound effect on their chromatin organization and transcriptional regulation.

Pseudotime analysis reveals novel regulatory factors for multigenic onset and monogenic transition of odorant receptor expression

During their maturation from horizontal basal stem cells, olfactory sensory neurons (OSNs) are known to select exactly one out of hundreds of olfactory receptors (ORs) and express it on their surface, a process called monogenic selection. Monogenic expression is preceded by a multigenic phase during which several OR genes are expressed in a single OSN. Here, we perform pseudotime analysis of a single cell RNA-Seq dataset of murine olfactory epithelium to precisely align the multigenic and monogenic expression phases with the cell types occurring during OSN differentiation. In combination with motif analysis of OR gene cluster-associated enhancer regions, we identify known and novel transcription (co-)factors (Ebf1, Lhx2, Ldb1, Fos and Ssbp2) and chromatin remodelers (Kdm1a, Eed and Zmynd8) associated with OR expression. The inferred temporal order of their activity suggests novel mechanisms contributing to multigenic OR expression and monogenic selection.

In silico analysis of promoter region and regulatory elements of glucan endo-1,3-beta-glucosidase encoding genes in Solanum tuberosum: cultivar DM 1-3 516 R44

Background

Potato (Solanum tuberosum L.) is one of the most important food crops in the world. Pathogens remain as one of the major constraints limiting potato productivity. Thus, understanding of gene regulation mechanism of pathogenesis-related genes such as glucan endo-1,3-beta-glucosidase is a foundation for genetic engineering of potato for disease resistance and reduces the use of fungicides. In the present study, 19 genes were selected and attempts were made through in silico methods to identify and characterize the promoter regions, regulatory elements, and CpG islands of glucan endo-1,3-beta-glucosidase gene in Solanum tuberosum cultivar DM 1-3 516 R44.

Results

The current analysis revealed that single transcription start sites (TSSs) were present in 12/19 (63.2%) of promoter regions analyzed. The predictive score at a cutoff value of 0.8 for the majority (84.2%) of the promoter regions ranged from 0.90 to 1.00. The locations for 42% of the TSSs were below −500 bp relative to the start codon (ATG). MβGII was identified as the common promoter motif for 94.4% of the genes with an E value of 3.5e−001. The CpG analysis showed low CpG density in the promoter regions of most of the genes except for gene ID102593331 and ID: 102595860. The number of SSRs per gene ranged from 2 to 9 with repeat lengths of 2 to 6 bp. Evolutionary distances ranged from 0.685 to 0.770 (mean = 0.73), demonstrating narrower genetic diversity range. Phylogeny was inferred using the UPGMA method, and gene sequences from different species were found to be clustered together.

Conclusion

In silico identified regulatory elements in promoter regions will contribute to our understanding of the regulatory mechanism of glucan endo-1,3-beta-glucosidase genes and provide a promising target for genetic engineering to improve disease resistance in potatoes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s43141-021-00240-0.

Olfactory expression of trace amine-associated receptors requires cooperative cis-acting enhancers

Olfactory sensory neurons express a large family of odorant receptors (ORs) and a small family of trace amine-associated receptors (TAARs). While both families are subject to so-called singular expression (expression of one allele of one gene), the mechanisms underlying TAAR gene choice remain obscure. Here, we report the identification of two conserved sequence elements in the mouse TAAR cluster (T-elements) that are required for TAAR gene expression. We observed that cell-type-specific expression of a TAAR-derived transgene required either T-element. Moreover, deleting either element reduced or abolished expression of a subset of TAAR genes, while deleting both elements abolished olfactory expression of all TAARs in cis with the mutation. The T-elements exhibit several features of known OR enhancers but also contain highly conserved, unique sequence motifs. Our data demonstrate that TAAR gene expression requires two cooperative cis-acting enhancers and suggest that ORs and TAARs share similar mechanisms of singular expression.

Genetic Background Effects on the Expression of an Odorant Receptor Gene

There are more than 1000 odorant receptor (OR) genes in the mouse genome. Each olfactory sensory neuron expresses only one of these genes, in a monoallelic fashion. The transcript abundance of homologous OR genes vary between distinct mouse strains. Here we analyzed the expression of the OR gene Olfr17 (also named P2) in different genomic contexts. Olfr17 is expressed at higher levels in the olfactory epithelium from 129 mice than from C57BL/6 (B6) mice. However, we found that in P2-IRES-tauGFP knock-in mice, the transcript levels of the 129 Olfr17 allele are highly reduced when compared to the B6 Olfr17 allele. To address the mechanisms involved in this variation we compared the 5′ region sequence and DNA methylation patterns of the B6 and 129 Olfr17 alleles. Our results show that genetic variations in cis regulatory regions can lead to differential DNA methylation frequencies in these OR gene alleles. They also show that expression of the Olfr17 alleles is largely affected by the genetic background, and suggest that in knock-in mice, expression can be affected by epigenetic modifications in the region of the targeted locus.

Maturation of the Olfactory Sensory Neuron and Its Cilia

Olfactory sensory neurons (OSNs) are bipolar neurons, unusual because they turn over continuously and have a multiciliated dendrite. The extensive changes in gene expression accompanying OSN differentiation in mice are largely known, especially the transcriptional regulators responsible for altering gene expression, revealing much about how differentiation proceeds. Basal progenitor cells of the olfactory epithelium transition into nascent OSNs marked by Cxcr4 expression and the initial extension of basal and apical neurites. Nascent OSNs become immature OSNs within 24-48 h. Immature OSN differentiation requires about a week and at least 2 stages. Early-stage immature OSNs initiate expression of genes encoding key transcriptional regulators and structural proteins necessary for further neuritogenesis. Late-stage immature OSNs begin expressing genes encoding proteins important for energy production and neuronal homeostasis that carry over into mature OSNs. The transition to maturity depends on massive expression of one allele of one odorant receptor gene, and this results in expression of the last 8% of genes expressed by mature OSNs. Many of these genes encode proteins necessary for mature function of axons and synapses or for completing the elaboration of non-motile cilia, which began extending from the newly formed dendritic knobs of immature OSNs. The cilia from adjoining OSNs form a meshwork in the olfactory mucus and are the site of olfactory transduction. Immature OSNs also have a primary cilium, but its role is unknown, unlike the critical role in proliferation and differentiation played by the primary cilium of the olfactory epithelium's horizontal basal cell.

The role of the NMD factor UPF3B in olfactory sensory neurons

The UPF3B-dependent branch of the nonsense-mediated RNA decay (NMD) pathway is critical for human cognition. Here, we examined the role of UPF3B in the olfactory system. Single-cell RNA-sequencing (scRNA-seq) analysis demonstrated considerable heterogeneity of olfactory sensory neuron (OSN) cell populations in wild-type (WT) mice, and revealed that UPF3B loss influences specific subsets of these cell populations. UPF3B also regulates the expression of a large cadre of antimicrobial genes in OSNs, and promotes the selection of specific olfactory receptor (Olfr) genes for expression in mature OSNs (mOSNs). RNA-seq and Ribotag analyses identified classes of mRNAs expressed and translated at different levels in WT and Upf3b-null mOSNs. Integrating multiple computational approaches, UPF3B-dependent NMD target transcripts that are candidates to mediate the functions of NMD in mOSNs were identified in vivo. Together, our data provides a valuable resource for the olfactory field and insights into the roles of NMD in vivo.

Depletion of Ric-8B leads to reduced mTORC2 activity

mTOR, a serine/threonine protein kinase that is involved in a series of critical cellular processes, can be found in two functionally distinct complexes, mTORC1 and mTORC2. In contrast to mTORC1, little is known about the mechanisms that regulate mTORC2. Here we show that mTORC2 activity is reduced in mice with a hypomorphic mutation of the Ric-8B gene. Ric-8B is a highly conserved protein that acts as a non-canonical guanine nucleotide exchange factor (GEF) for heterotrimeric Gαs/olf type subunits. We found that Ric-8B hypomorph embryos are smaller than their wild type littermates, fail to close the neural tube in the cephalic region and die during mid-embryogenesis. Comparative transcriptome analysis revealed that signaling pathways involving GPCRs and G proteins are dysregulated in the Ric-8B mutant embryos. Interestingly, this analysis also revealed an unexpected impairment of the mTOR signaling pathway. Phosphorylation of Akt at Ser473 is downregulated in the Ric-8B mutant embryos, indicating a decreased activity of mTORC2. Knockdown of the endogenous Ric-8B gene in cultured cell lines leads to reduced phosphorylation levels of Akt (Ser473), further supporting the involvement of Ric-8B in mTORC2 activity. Our results reveal a crucial role for Ric-8B in development and provide novel insights into the signals that regulate mTORC2.

Gene inactivation in mice can be used to identify genes that are involved in important biological processes and that may contribute to disease. We used this approach to study the Ric-8B gene, which is highly conserved in mammals, including humans. We found that Ric-8B is essential for embryogenesis and for the proper development of the nervous system. Ric-8B mutant mouse embryos are smaller than their wild type littermates and show neural tube defects at the cranial region. This approach also allowed us to identify the biological pathways that potentially contribute to the observed phenotypes, and uncover a novel role for Ric-8B in the mTORC2 signaling pathway. mTORC2 plays particular important roles in the adult brain, and has been implicated in neurological disorders. Our mutant mice provide a model to study the complex molecular and cellular processes underlying the interplay between Ric-8B and mTORC2 in neuronal function.

In silico analysis of the promoter region of olfactory receptors in cattle (Bos indicus) to understand its gene regulation

Identifications of transcription start sites (TSSs) and promoter regions are first step to understand the regulation mechanisms of gene expression and association with genetic variations in the regions. This analysis was conducted with the objectives to identify TSSs, determine the promoter regions, identify common candidate motifs and transcription factors (TFs), and search for CpG islands (CGIs) in cattle olfactory receptors (ORs) genes promoter regions. In the analysis, TSSs of cattle olfactory genes were first identified. The locations for 60% of the TSSs were below -500 bp relative to the start codon and five candidate motifs (MOR1, MOR2, MOR3, MOR4, and MOR5) were identified that are shared by at least 50% of the cattle ORs promoter input sequences from both strands. Among the five candidate motifs, MOR4 was revealed as the common promoter motif for 85.71% of cattle ORs genes that serves as binding sites for TFs involved in the expression regulation of these genes. MOR4 was also compared to registered motifs in publically available databases to see if they are similar to known regulatory motifs for TF by using the TOMTOM web application. Hence, it was revealed that MOR4 may serve as the binding site mainly for the Zinc finger (ZNF) TF gene family to regulate expression of cattle ORs genes. Further gene ontology analysis for MOR4 demonstrated ORs belong to the G-protein-coupled receptor superfamily and MOR4 tend to be located near the genes involved in the detection of chemical stimulus involved in sensory perception and in innate immune responses such as cytokine-mediated signaling. In silico digestion of cattle OR sequences was performed using restriction enzyme MspI. CGIs from OR10K1 and OR2L13 gene was found. In the present analysis, the poor CGIs observed might suggest their gene expression regulation pattern is in tissue specific manner.

Alternative polyadenylation produces multiple 3' untranslated regions of odorant receptor mRNAs in mouse olfactory sensory neurons

BACKGROUND

Odorant receptor genes constitute the largest gene family in mammalian genomes and this family has been extensively studied in several species, but to date far less attention has been paid to the characterization of their mRNA 3' untranslated regions (3'UTRs). Given the increasing importance of UTRs in the understanding of RNA metabolism, and the growing interest in alternative polyadenylation especially in the nervous system, we aimed at identifying the alternative isoforms of odorant receptor mRNAs generated through 3'UTR variation.

RESULTS

We implemented a dedicated pipeline using IsoSCM instead of Cufflinks to analyze RNA-Seq data from whole olfactory mucosa of adult mice and obtained an extensive description of the 3'UTR isoforms of odorant receptor mRNAs. To validate our bioinformatics approach, we exhaustively analyzed the 3'UTR isoforms produced from 2 pilot genes, using molecular approaches including northern blot and RNA ligation mediated polyadenylation test. Comparison between datasets further validated the pipeline and confirmed the alternative polyadenylation patterns of odorant receptors. Qualitative and quantitative analyses of the annotated 3' regions demonstrate that 1) Odorant receptor 3'UTRs are longer than previously described in the literature; 2) More than 77% of odorant receptor mRNAs are subject to alternative polyadenylation, hence generating at least 2 detectable 3'UTR isoforms; 3) Splicing events in 3'UTRs are restricted to a limited subset of odorant receptor genes; and 4) Comparison between male and female data shows no sex-specific differences in odorant receptor 3'UTR isoforms.

CONCLUSIONS

We demonstrated for the first time that odorant receptor genes are extensively subject to alternative polyadenylation. This ground-breaking change to the landscape of 3'UTR isoforms of Olfr mRNAs opens new avenues for investigating their respective functions, especially during the differentiation of olfactory sensory neurons.

Analysis of Pig Vomeronasal Receptor Type 1 (V1R) Promoter Region Reveals a Common Promoter Motif but Poor CpG Islands

Promoters are, generally, located immediately upstream of a transcription start site (TSS) and have a variety of regulatory motifs, such as transcription factors (TFs) and CpG islands (CGIs), that participate in the regulation of gene expression. Here analysis of the promoter region for pig vomeronasal receptor type 1 (V1R) was described. In the analysis, TSSs for pig V1R genes was first identified and five motifs (MV1, MV2, MV3, MV4, and MV5) were found that are shared by at least 50% of the pig V1R promoter input sequences from both strands. Among the five motifs, MV2 was identified as a common promoter motif shared by all (100%) pig V1R promoters. For further analysis, to better characterize and get deeper biological insight associated with MV2, TOMTOM web application was used. MV2 was compared to the known motif databases (such as JASPAR) to see if they are similar to a known regulatory motif (transcription factor). Hence, it was revealed that MV2 serves as the binding site mainly for the BetaBetaAlpha-zinc finger (BTB-ZF) transcription factor gene family to regulate expression of pig V1R genes. Moreover, it was shown that pig V1R promoters are CpG poor, suggesting that their gene expression regulation pattern is in tissue specific manner.

Cooperative interactions enable singular olfactory receptor expression in mouse olfactory neurons

The monogenic and monoallelic expression of only one out of >1000 mouse olfactory receptor (ORs) genes requires the formation of large heterochromatic chromatin domains that sequester the OR gene clusters. Within these domains, intergenic transcriptional enhancers evade heterochromatic silencing and converge into interchromosomal hubs that assemble over the transcriptionally active OR. The significance of this nuclear organization in OR choice remains elusive. Here, we show that transcription factors Lhx2 and Ebf specify OR enhancers by binding in a functionally cooperative fashion to stereotypically spaced motifs that defy heterochromatin. Specific displacement of Lhx2 and Ebf from OR enhancers resulted in pervasive, long-range, and trans downregulation of OR transcription, whereas pre-assembly of a multi-enhancer hub increased the frequency of OR choice in cis. Our data provide genetic support for the requirement and sufficiency of interchromosomal interactions in singular OR choice and generate general regulatory principles for stochastic, mutually exclusive gene expression programs.

Lysine-specific demethylase-1 (LSD1) depletion disrupts monogenic and monoallelic odorant receptor (OR) expression in an olfactory neuronal cell line

Function of the mammalian olfactory system depends on specialized olfactory sensory neurons (OSNs) that each express only one allele ("monoallelic") of one odorant receptor (OR) gene ("monogenic"). The lysine-specific demethylase-1 (LSD1) protein removes activating H3K4 or silencing H3K9 methylation marks in a variety of developmental contexts, and is thought to be important for proper OR regulation. Most of the focus in the field has been on a potential "activating" function for LSD1; e.g., in the demethylation of H3K9 associated with the expressed OR allele. Here we show that depletion of LSD1 in an immortalized olfactory-placode-derived cell line (OP6) results in multigenic and multiallelic OR transcription per cell, while not seemingly disrupting the ability of these cells to activate new OR genes during clonal expansion. These results are consistent with LSD1 having a role in silencing additional OR alleles, as opposed to being required for the activation of OR alleles, within the OP6 cellular context.

Regulatory Features for Odorant Receptor Genes in the Mouse Genome

The odorant receptor genes, seven transmembrane receptor genes constituting the vastest mammalian gene multifamily, are expressed monogenically and monoallelicaly in each sensory neuron in the olfactory epithelium. This characteristic, often referred to as the one neuron-one receptor rule, is driven by mostly uncharacterized molecular dynamics, generally named odorant receptor gene choice. Much attention has been paid by the scientific community to the identification of sequences regulating the expression of odorant receptor genes within their loci, where related genes are usually arranged in genomic clusters. A number of studies identified transcription factor binding sites on odorant receptor promoter sequences. Similar binding sites were also found on a number of enhancers that regulate in cis their transcription, but have been proposed to form interchromosomal networks. Odorant receptor gene choice seems to occur via the local removal of strongly repressive epigenetic markings, put in place during the maturation of the sensory neuron on each odorant receptor locus. Here we review the fast-changing state of art for the study of regulatory features for odorant receptor genes.

Lhx2 Determines Odorant Receptor Expression Frequency in Mature Olfactory Sensory Neurons

A developmental program of epigenetic repression prepares each mammalian olfactory sensory neuron (OSN) to strongly express one allele from just one of hundreds of odorant receptor (OR) genes, but what completes this process of OR gene choice by driving the expression of this allele is incompletely understood. Conditional deletion experiments in mice demonstrate that Lhx2 is necessary for normal expression frequencies of nearly all ORs and all trace amine-associated receptors, irrespective of whether the deletion of Lhx2 is initiated in immature or mature OSNs. Given previous evidence that Lhx2 binds OR gene control elements, these findings indicate that Lhx2 is directly involved in driving OR expression. The data also support the conclusion that OR expression is necessary to allow immature OSNs to complete differentiation and become mature. In contrast to the robust effects of conditional deletion of Lhx2, the loss of Emx2 has much smaller effects and more often causes increased expression frequencies. Lhx2:Emx2 double mutants show opposing effects on Olfr15 expression that reveal independent effects of these two transcription factors. While Lhx2 is necessary for OR expression that supports OR gene choice, Emx2 can act differently; perhaps by helping to control the availability of OR genes for expression.

Monogenic and Monoallelic Expression of Odorant Receptors

Odorant receptors (ORs) belong to a large gene family of rhodopsin-like G protein-coupled receptors (GPCRs). The mouse OR gene family is composed of ∼1000 OR genes, and the human OR gene family is composed of ∼400 OR genes. The OR genes are spread throughout the genome, and can be found in clusters or as solitary genes in almost all chromosomes. These chemosensory GPCRs are expressed in highly specialized cells, the olfactory sensory neurons of the nose. Each one of these neurons expresses a single OR gene out of the complete repertoire of genes. In addition, only one of the two homologous alleles of the chosen OR gene, the maternal or the paternal, is expressed per neuron. Here we review recent findings that help to elucidate the mechanisms underlying monogenic and monoallelic expression of OR genes.

Expression of homeobox genes in the mouse olfactory epithelium

Homeobox genes constitute a large family of genes widely studied because of their role in the establishment of the body pattern. However, they are also involved in many other events during development and adulthood. The main olfactory epithelium (MOE) is an excellent model to study neurogenesis in the adult nervous system. Analyses of homeobox genes during development show that some of these genes are involved in the formation and establishment of cell diversity in the MOE. Moreover, the mechanisms of expression of odorant receptors (ORs) constitute one of the biggest enigmas in the field. Analyses of OR promoters revealed the presence of homeodomain binding sites in their sequences. Here we characterize the expression patterns of a set of 49 homeobox genes in the MOE with in situ hybridization. We found that seven of them (Dlx3, Dlx5, Dlx6, Msx1, Meis1, Isl1, and Pitx1) are zonally expressed. The homeobox gene Emx1 is expressed in three guanylate cyclase(+) populations, two located in the MOE and the third one in an olfactory subsystem known as Grüneberg ganglion located at the entrance of the nasal cavity. The homeobox gene Tshz1 is expressed in a unique patchy pattern across the MOE. Our findings provide new insights to guide functional studies that aim to understand the complexity of transcription factor expression and gene regulation in the MOE. J. Comp. Neurol. 524:2713-2739, 2016. © 2016 Wiley Periodicals, Inc.

Monoallelic expression of olfactory receptors

The sense of smell collects vital information about the environment by detecting a multitude of chemical odorants. Breadth and sensitivity are provided by a huge number of chemosensory receptor proteins, including more than 1,400 olfactory receptors (ORs). Organizing the sensory information generated by these receptors so that it can be processed and evaluated by the central nervous system is a major challenge. This challenge is overcome by monogenic and monoallelic expression of OR genes. The single OR expressed by each olfactory sensory neuron determines the neuron's odor sensitivity and the axonal connections it will make to downstream neurons in the olfactory bulb. The expression of a single OR per neuron is accomplished by coupling a slow chromatin-mediated activation process to a fast negative-feedback signal that prevents activation of additional ORs. Singular OR activation is likely orchestrated by a network of interchromosomal enhancer interactions and large-scale changes in nuclear architecture.

Nuclear architecture and gene silencing in olfactory sensory neurons

Odorants are discriminated by hundreds of odorant receptor (OR) genes, which are dispersed throughout the mammalian genome. The OR genes are expressed in a highly specialized type of cell, the olfactory sensory neuron. Each one of these neurons expresses one of the 2 alleles from one single OR gene type. The mechanisms underlying OR gene expression are unclear. Here we describe recent work demonstrating that the olfactory sensory neuron shows a particular nuclear architecture, and that the genomic OR loci are colocalized in silencing heterochromatin compartments within the nucleus. These discoveries highlight the important role played by epigenetic modifications and nuclear genome organization in the regulation of OR gene expression.

Identification of novel post-transcriptional features in olfactory receptor family mRNAs

Olfactory receptor (Olfr) genes comprise the largest gene family in mice. Despite their importance in olfaction, how most Olfr mRNAs are regulated remains unexplored. Using RNA-seq analysis coupled with analysis of pre-existing databases, we found that Olfr mRNAs have several atypical features suggesting that post-transcriptional regulation impacts their expression. First, Olfr mRNAs, as a group, have dramatically higher average AU-content and lower predicted secondary structure than do control mRNAs. Second, Olfr mRNAs have a higher density of AU-rich elements (AREs) in their 3'UTR and upstream open reading frames (uORFs) in their 5 UTR than do control mRNAs. Third, Olfr mRNAs have shorter 3' UTR regions and with fewer predicted miRNA-binding sites. All of these novel properties correlated with higher Olfr expression. We also identified striking differences in the post-transcriptional features of the mRNAs from the two major classes of Olfr genes, a finding consistent with their independent evolutionary origin. Together, our results suggest that the Olfr gene family has encountered unusual selective forces in neural cells that have driven them to acquire unique post-transcriptional regulatory features. In support of this possibility, we found that while Olfr mRNAs are degraded by a deadenylation-dependent mechanism, they are largely protected from this decay in neural lineage cells.

Chemosensory receptor specificity and regulation

The senses provide a means by which data on the physical and chemical properties of the environment may be collected and meaningfully interpreted. Sensation begins at the periphery, where a multitude of different sensory cell types are activated by environmental stimuli as different as photons and odorant molecules. Stimulus sensitivity is due to expression of different cell surface sensory receptors, and therefore the receptive field of each sense is defined by the aggregate of expressed receptors in each sensory tissue. Here, we review current understanding on patterns of expression and modes of regulation of sensory receptors.

Developmental regulation of neural map formation in the mouse olfactory system

In the mouse olfactory system, various odorants are detected by approximately 1000 different odorant receptors (ORs) expressed in the olfactory sensory neurons (OSNs). It is well established that each OSN expresses only one functional OR gene in a monoallelic manner. Furthermore, OSN axons expressing the same OR converge to a set of glomeruli in the olfactory bulb (OB). During embryonic development, a coarse map is formed by the combination of two genetically programmed processes. One is OR-independent axonal projection along the dorsal-ventral (D-V) axis, and the other is OR-dependent projection along the anterior-posterior (A-P) axis. D-V projection is regulated by the anatomical location of OSNs within the olfactory epithelium (OE), whereas A-P projection is instructed by expressed OR molecules using cyclic adenosine monophosphate (cAMP) signals. After birth, the map is further refined in an activity-dependent manner by its conversion from a continuous to a discrete map through segregation of glomerular structures. Here, we summarize recent progress from our laboratory in understanding neural map formation in the mouse olfactory system.

Mammalian olfactory receptors: molecular mechanisms of odorant detection, 3D-modeling, and structure-activity relationships

This chapter describes the main characteristics of olfactory receptor (OR) genes of vertebrates, including generation of this large multigenic family and pseudogenization. OR genes are compared in relation to evolution and among species. OR gene structure and selection of a given gene for expression in an olfactory sensory neuron (OSN) are tackled. The specificities of OR proteins, their expression, and their function are presented. The expression of OR proteins in locations other than the nasal cavity is regulated by different mechanisms, and ORs display various additional functions. A conventional olfactory signal transduction cascade is observed in OSNs, but individual ORs can also mediate different signaling pathways, through the involvement of other molecular partners and depending on the odorant ligand encountered. ORs are engaged in constitutive dimers. Ligand binding induces conformational changes in the ORs that regulate their level of activity depending on odorant dose. When present, odorant binding proteins induce an allosteric modulation of OR activity. Since no 3D structure of an OR has been yet resolved, modeling has to be performed using the closest G-protein-coupled receptor 3D structures available, to facilitate virtual ligand screening using the models. The study of odorant binding modes and affinities may infer best-bet OR ligands, to be subsequently checked experimentally. The relationship between spatial and steric features of odorants and their activity in terms of perceived odor quality are also fields of research that development of computing tools may enhance.

Nuclear compartmentalization of odorant receptor genes

Odorants are detected by odorant receptors, which are located on olfactory sensory neurons of the nose. Each olfactory sensory neuron expresses one single odorant receptor gene allele from a large family of odorant receptor genes. To gain insight into the mechanisms underlying this monogenic and monoallelic expression, we examined the 3D nuclear organization of olfactory sensory neurons and determined the positions of homologous odorant receptor gene alleles in relation to different nuclear compartments. Our results show that olfactory neurons exhibit a singular nuclear architecture that is characterized by a large centrally localized constitutive heterochromatin block and by the presence of prominent facultative heterochromatin domains that are localized around this constitutive heterochromatin block. We also found that the two homologous alleles of a given odorant receptor gene are frequently segregated to separate compartments in the nucleus, with one of the alleles localized to the constitutive heterochromatin block and the other one localized to the more plastic facultative heterochromatin, or next to it. Our findings suggest that this nuclear compartmentalization may play a critical role in the expression of odorant receptor genes.

Extrasensory perception: odorant and taste receptors beyond the nose and mouth

G protein-coupled receptors (GPCRs) represent the largest family of transmembrane receptors and are prime therapeutic targets. The odorant and taste receptors account for over half of the GPCR repertoire, yet they are generally excluded from large-scale, drug candidate analyses. Accumulating molecular evidence indicates that the odorant and taste receptors are widely expressed throughout the body and functional beyond the oronasal cavity - with roles including nutrient sensing, autophagy, muscle regeneration, regulation of gut motility, protective airway reflexes, bronchodilation, and respiratory disease. Given this expanding array of actions, the restricted perception of these GPCRs as mere mediators of smell and taste is outdated. Moreover, delineation of the precise actions of odorant and taste GPCRs continues to be hampered by the relative paucity of selective and specific experimental tools, as well as the lack of defined receptor pharmacology. In this review, we summarize the evidence for expression and function of odorant and taste receptors in tissues beyond the nose and mouth, and we highlight their broad potential in physiology and pathophysiology.

Single cell RT-PCR identification of odorant receptors expressed by olfactory neurons

Mammals have between 400 and 1,300 functional odorant receptor (OR) genes in their genomes. Each olfactory sensory neuron in the nose expresses only one single type of OR out of this vast repertoire. The OR expressed by an olfactory sensory neuron determines its functional activity and wiring to the olfactory bulb. Therefore, the identification of the OR type expressed by individual neurons contributes to the understanding of the mechanisms of odorant perception. Here we describe the protocol that combines single cell RT-PCR and degenerate PCR to identify the OR expressed by single olfactory neurons. In a primary PCR reaction, cDNAs corresponding to all mRNAs present in a single neuron are synthesized. In a secondary PCR reaction, the cDNAs corresponding to the OR expressed by this neuron are amplified using degenerate primers that match OR sequences.

An epigenetic trap stabilizes singular olfactory receptor expression

The molecular mechanisms regulating olfactory receptor (OR) expression in the mammalian nose are not yet understood. Here, we identify the transient expression of histone demethylase LSD1 and the OR-dependent expression of adenylyl cyclase 3 (Adcy3) as requirements for initiation and stabilization of OR expression. As a transcriptional coactivator, LSD1 is necessary for desilencing and initiating OR transcription, but as a transcriptional corepressor, it is incompatible with maintenance of OR expression, and its downregulation is imperative for stable OR choice. Adcy3, a sensor of OR expression and a transmitter of an OR-elicited feedback, mediates the downregulation of LSD1 and promotes the differentiation of olfactory sensory neurons (OSNs). This novel, three-node signaling cascade locks the epigenetic state of the chosen OR, stabilizes its singular expression, and prevents the transcriptional activation of additional OR alleles for the life of the neuron.

HORDE: comprehensive resource for olfactory receptor genomics

Olfactory receptors (ORs) constitute the largest gene family in the mammalian genome. The existence of these proteins underlies the nature of, and variability in, odorant perception. The Human Olfactory Receptor Data Explorer (HORDE, http://genome.weizmann.ac.il/horde/ ) is a free online resource, which presents a complete compendium of all OR genes and pseudogenes in the genome of human and four other vertebrates. HORDE includes three parts: (1) an automated pipeline, which mines OR gene and pseudogene sequences out of complete genomes, and generates gene symbols based on sequence similarity; (2) a card generator that obtains and displays annotative information on individual ORs retrieved from external databases and relevant studies; and (3) a search engine that allows user retrieval of OR information. For human ORs, HORDE specifically addresses the universe of interindividual variation, as obtained from several sources, including whole genome sequences made possible by next-generation sequencing. This encompasses single nucleotide polymorphisms (SNP) and copy number variation (CNV), including deleterious mutational events. HORDE also hosts a number of tools designed specifically to assist in the study of OR evolution and function. In this chapter, we describe the status of HORDE (build #43). We also discuss plans for future enhancements and a road map for HORDE to become a better community-based bioinformatics tool. We highlight HORDE's role as a major research tool in the study of an expanding cohort of OR repertoires.

Common promoter elements in odorant and vomeronasal receptor genes

In mammals, odorants and pheromones are detected by hundreds of odorant receptors (ORs) and vomeronasal receptors (V1Rs and V2Rs) expressed by sensory neurons that are respectively located in the main olfactory epithelium and in the vomeronasal organ. Even though these two olfactory systems are functionally and anatomically separate, their sensory neurons show a common mechanism of receptor gene regulation: each neuron expresses a single receptor gene from a single allele. The mechanisms underlying OR and VR gene expression remain unclear. Here we investigated if OR and V1R genes share common sequences in their promoter regions.

We conducted a comparative analysis of promoter regions of 39 mouse V1R genes and found motifs that are common to a large number of promoters. We then searched mouse OR promoter regions for motifs that resemble the ones found in the V1R promoters. We identified motifs that are present in both the V1R and OR promoter regions. Some of these motifs correspond to the known O/E like binding sites while others resemble binding sites for transcriptional repressors. We show that one of these motifs specifically interacts with proteins extracted from both nuclei from olfactory and vomeronasal neurons. Our study is the first to identify motifs that resemble binding sites for repressors in the promoters of OR and V1R genes. Analysis of these motifs and of the proteins that bind to these motifs should reveal important aspects of the mechanisms of OR/V1R gene regulation.

Promoter architecture of mouse olfactory receptor genes

Odorous chemicals are detected by the mouse main olfactory epithelium (MOE) by about 1100 types of olfactory receptors (OR) expressed by olfactory sensory neurons (OSNs). Each mature OSN is thought to express only one allele of a single OR gene. Major impediments to understand the transcriptional control of OR gene expression are the lack of a proper characterization of OR transcription start sites (TSSs) and promoters, and of regulatory transcripts at OR loci. We have applied the nanoCAGE technology to profile the transcriptome and the active promoters in the MOE. nanoCAGE analysis revealed the map and architecture of promoters for 87.5% of the mouse OR genes, as well as the expression of many novel noncoding RNAs including antisense transcripts. We identified candidate transcription factors for OR gene expression and among them confirmed by chromatin immunoprecipitation the binding of TBP, EBF1 (OLF1), and MEF2A to OR promoters. Finally, we showed that a short genomic fragment flanking the major TSS of the OR gene Olfr160 (M72) can drive OSN-specific expression in transgenic mice.

Rigorous and thorough bioinformatic analyses of olfactory receptor promoters confirm enrichment of O/E and homeodomain binding sites but reveal no new common motifs

Background

Mammalian olfactory receptors (ORs) are subject to a remarkable but poorly understood regime of transcriptional regulation, whereby individual olfactory neurons each express only one allele of a single member of the large OR gene family.

Results

We performed a rigorous search for enriched sequence motifs in the largest dataset of OR promoter regions analyzed to date. We combined measures of cross-species conservation with databases of known transcription factor binding sites and ab initio motif-finding algorithms. We found strong enrichment of binding sites for the O/E family of transcription factors and for homeodomain factors, both already known to be involved in the transcriptional control of ORs, but did not identify any novel enriched sequences. We also found that TATA-boxes are present in at least a subset of OR promoters.

Conclusions

Our rigorous approach provides a template for the analysis of the regulation of large gene families and demonstrates some of the difficulties and pitfalls of such analyses. Although currently available bioinformatics methods cannot detect all transcriptional regulatory elements, our thorough analysis of OR promoters shows that in the case of this gene family, experimental approaches have probably already identified all the binding factors common to large fractions of OR promoters.

How is the olfactory map formed and interpreted in the mammalian brain?

Odor signals received by odorant receptors (ORs) expressed by olfactory sensory neurons (OSNs) in the olfactory epithelium (OE) are represented as an odor map in the olfactory bulb (OB). In the mouse, there are ~1,000 different OR species, and each OSN expresses only one functional OR gene in a monoallelic manner. Furthermore, OSN axons expressing the same type of OR converge on a specific target site in the OB, forming a glomerular structure. Because each glomerulus represents a single OR species, and a single odorant can interact with multiple OR species, odor signals received in the OE are converted into a topographic map of multiple glomeruli activated with varying magnitudes. Here we review recent progress in the study of the mammalian olfactory system, focusing on the formation of the olfactory map and the transmission of topographical information in the OB to the olfactory cortex to elicit various behaviors.

Simple Shared Motifs (SSM) in conserved region of promoters: a new approach to identify co-regulation patterns

Background

Regulation of gene expression plays a pivotal role in cellular functions. However, understanding the dynamics of transcription remains a challenging task. A host of computational approaches have been developed to identify regulatory motifs, mainly based on the recognition of DNA sequences for transcription factor binding sites. Recent integration of additional data from genomic analyses or phylogenetic footprinting has significantly improved these methods.

Results

Here, we propose a different approach based on the compilation of Simple Shared Motifs (SSM), groups of sequences defined by their length and similarity and present in conserved sequences of gene promoters. We developed an original algorithm to search and count SSM in pairs of genes. An exceptional number of SSM is considered as a common regulatory pattern. The SSM approach is applied to a sample set of genes and validated using functional gene-set enrichment analyses. We demonstrate that the SSM approach selects genes that are over-represented in specific biological categories (Ontology and Pathways) and are enriched in co-expressed genes. Finally we show that genes co-expressed in the same tissue or involved in the same biological pathway have increased SSM values.

Conclusions

Using unbiased clustering of genes, Simple Shared Motifs analysis constitutes an original contribution to provide a clearer definition of expression networks.

High-throughput mapping of the promoters of the mouse olfactory receptor genes reveals a new type of mammalian promoter and provides insight into olfactory receptor gene regulation

The olfactory receptor (OR) genes are the largest mammalian gene family and are expressed in a monogenic and monoallelic fashion in olfactory neurons. Using a high-throughput approach, we mapped the transcription start sites of 1085 of the 1400 murine OR genes and performed computational analysis that revealed potential transcription factor binding sites shared by the majority of these promoters. Our analysis produced a hierarchical model for OR promoter recognition in which unusually high AT content, a unique epigenetic signature, and a stereotypically positioned O/E site distinguish OR promoters from the rest of the murine promoters. Our computations revealed an intriguing correlation between promoter AT content and evolutionary plasticity, as the most AT-rich promoters regulate rapidly evolving gene families. Within the AT-rich promoter category the position of the TATA-box does not correlate with the transcription start site. Instead, a spike in GC composition might define the exact location of the TSS, introducing the concept of "genomic contrast" in transcriptional regulation. Finally, our experiments show that genomic neighborhood rather than promoter sequence correlates with the probability of different OR genes to be expressed in the same olfactory cell.

Expression of the IP3R1 promoter-driven nls-lacZ transgene in Purkinje cell parasagittal arrays of developing mouse cerebellum

The cerebellar Purkinje cell monolayer is organized into heterogeneous Purkinje cell compartments that have different molecular compositions. Here we describe a transgenic mouse line, 1NM13, that shows heterogeneous transgene expression in parasagittal Purkinje cell arrays. The transgene consists of a nuclear localization signal (nls) fused to the beta-galactosidase (lacZ) composite gene driven by the type 1 inositol 1,4,5-trisphosphate receptor (IP(3)R1) gene promoter. IP(3)R1-nls-lacZ transgene expression was detected at a single Purkinje cell level over the surface of a whole-mount X-gal-stained cerebellum because of nuclear accumulation of the nls-lacZ activity. Developing cerebella of 1NM13 mice showed stripe-like X-gal staining patterns of parasagittal Purkinje cell subsets. The X-gal stripe pattern was likely determined by an intrinsic property as early as E15 and showed increasing complexity with cerebellar development. The X-gal stripe pattern was reminiscent of, but not identical to, the stripe pattern of zebrin II immunoreactivity. We designated the symmetrical X-gal-positive (transgene-positive, Tg(+)) Purkinje cell stripes about the midline as vermal Tg1(+), Tg2(a, b)(+) and Tg3(a, b)(+) stripes and hemispheric Tg4(a, b)(+), Tg5(a, b)(+), Tg6(a, b, c)(+), and Tg7(a, b)(+) stripes, where a, b, and c indicate substripes. We also assigned three parafloccular substripes Tg8(a, b, c)(+). The boundaries of X-gal stripes at P5 were consistent with raphes in the Purkinje cell layer through which granule cells migrate, suggesting a possible association of the X-gal stripes with raphe formation. Our results indicate that 1NM13 is a good mouse model with a reproducible and clear marker for the compartmentalization of Purkinje cell arrays.

Homeodomain binding motifs modulate the probability of odorant receptor gene choice in transgenic mice

Odorant receptor (OR) genes constitute with 1200 members the largest gene family in the mouse genome. A mature olfactory sensory neuron (OSN) is thought to express just one OR gene, and from one allele. The cell bodies of OSNs that express a given OR gene display a mosaic pattern within a particular region of the main olfactory epithelium. The mechanisms and cis-acting DNA elements that regulate the expression of one OR gene per OSN - OR gene choice - remain poorly understood. Here, we describe a reporter assay to identify minimal promoters for OR genes in transgenic mice, which are produced by the conventional method of pronuclear injection of DNA. The promoter transgenes are devoid of an OR coding sequence, and instead drive expression of the axonal marker tau-β-galactosidase. For four mouse OR genes (M71, M72, MOR23, and P3) and one human OR gene (hM72), a mosaic, OSN-specific pattern of reporter expression can be obtained in transgenic mice with contiguous DNA segments of only ~300 bp that are centered around the transcription start site (TSS). The ~150bp region upstream of the TSS contains three conserved sequence motifs, including homeodomain (HD) binding sites. Such HD binding sites are also present in the H and P elements, DNA sequences that are known to strongly influence OR gene expression. When a 19mer encompassing a HD binding site from the P element is multimerized nine times and added upstream of a MOR23 minigene that contains the MOR23 coding region, we observe a dramatic increase in the number of transgene-expressing founders and lines and in the number of labeled OSNs. By contrast, a nine times multimerized 19mer with a mutant HD binding site does not have these effects. We hypothesize that HD binding sites in the H and P elements and in OR promoters modulate the probability of OR gene choice.

Neural map formation in the mouse olfactory system

Odor signals received by odorant receptors (ORs) in the olfactory epithelium are represented as an odor map of activated glomeruli in the olfactory bulb. In the mouse olfactory system, it appears that much of axon pathfinding and sorting occurs autonomously by olfactory neuron axons. Here, we review the recent progress on the study of olfactory map formation in rodents. We will discuss how neuronal identity is represented at axon termini and how the OR-instructed axonal projection is regulated.

Achieving singularity in mammalian odorant receptor gene choice

Odor discrimination requires differential expression of odor detectors. In fact, olfactory input to the brain is organized in units (glomeruli) innervated only by olfactory sensory neurons that express the same odorant receptor (OR). Therefore, discriminatory capacity is maximized if each sensory neuron expresses only one allele of a single OR gene, a postulate sometimes canonized as the "one neuron-one receptor rule." OR gene choice appears to result from a hierarchy of processes: differential availability of the alleles of each OR gene, zonal exclusion (or selection), OR gene switching during the initiation of OR gene transcription, and OR-dependent feedback to solidify the choice of one OR gene. The mechanisms underlying these processes are poorly understood, though a few elements are known or suspected. For example, the mechanism of activation of OR gene transcription appears to work in part through a few homeobox transcription factors (Emx2, and perhaps Lhx2) and the Ebf family of transcription factors. Further insights will probably come from several directions, but a promising hypothesis is that epigenetic mechanisms contribute to all levels of the hierarchical control of OR gene expression, especially the repressive events that seem to be necessary to achieve the singularity of OR gene choice.

Mechanisms of odorant receptor gene choice in Drosophila and vertebrates

Odorant receptors are encoded by extremely large and divergent families of genes. Each receptor is expressed in a small proportion of neurons in the olfactory organs, and each neuron in turn expresses just one odorant receptor gene. This fundamental property of the peripheral olfactory system is widely conserved across evolution, and observed in vertebrates, like mice, and invertebrates, like Drosophila, despite their olfactory receptor gene families being evolutionarily unrelated. Here we review the progress that has been made in these two systems to understand the intriguing and elusive question: how does a single neuron choose to express just one of many possible odorant receptors and exclude expression of all others?

Emx2 stimulates odorant receptor gene expression

The mechanisms selecting a single odorant receptor (OR) gene for expression in each olfactory sensory neuron (OSN) establish an OR expression pattern critical for odor discrimination. These mechanisms are largely unknown, but putative OR promoters contain homeodomain-like sites, implicating homeobox transcription factors such as Emx2. At embryonic day 18.5, expression of 49-76% of ORs was decreased in mice lacking Emx2, depending on the metric used. The decreases were due to fewer OSNs expressing each OR. Affected ORs showed changes that were disproportionately greater than the 42% reduction in mature neurons and similar decreases in unrelated olfactory neuron-enriched messenger RNAs in Emx2(-/-) mice. Both Class I and Class II ORs decreased, as did ORs expressed in both the dorsal and ventral regions of the epithelium. Conversely, 7% of Class II ORs tested were expressed more frequently, suggesting that some ORs are independent of Emx2. Emx2 helps stimulate transcription for many OR genes, which we hypothesize is through direct action at OR promoters, but Emx2 appears to have no significant role in regulating other aspects of OR gene expression, including the zonal patterns, OR gene cluster selection mechanisms, and singularity of OR gene choice.

Deletion of the core-H region in mice abolishes the expression of three proximal odorant receptor genes in cis

We have previously reported that a 2.1-kb homology (H) sequence, conserved between mouse and human, regulates the odorant receptor (OR) gene MOR28 in transgenic mice. Here, we narrowed down the essential sequences of the H to a core of 124 bp by using a transient expression system in zebrafish embryos. Transgenic experiments in mice demonstrated that the core-H sequence is sufficient to endow expression of the MOR28 minigene. Deletion and mutation analyses of the core-H region revealed two homeodomain sequences to be essential for the H enhancer activity. Targeted deletion of the core-H abolished expression of three proximal OR genes, MOR28, MOR10, and MOR83, in cis, indicating the presence of another locus control region/enhancer in the downstream region, that regulates four distal OR genes in the same MOR28 cluster. In the heterozygous mice, the H(-) phenotype of the mutant allele was not rescued by the wild-type H(+) allele in trans.

Odorant and pheromone receptor gene regulation in vertebrates

The largest mammalian gene family codes for odorant receptors and is exclusively devoted to the perception of the outside world. Its expression is very peculiar, since olfactory sensory neurons are only allowed to express a single of its numerous members, from a single parental allele. How this is achieved is unknown, but recent work points to multiple regulatory mechanisms, possibly shared by pheromone receptor genes, acting at (a) a general level, via the expression of the chemoreceptor itself and (b) a more restricted level, defined by activator elements.

Promotor elements governing the clustered expression pattern of odorant receptor genes

Odorant receptor (OR) genes of family mOR262 are only expressed in olfactory sensory neurons (OSNs) segregated in a central patch of the nasal turbinates; they comprise conserved DNA elements upstream of their transcription start sites that are proposed to govern the distinct expression pattern. In mouse lines with a transgene containing the coding sequence and a short upstream region of the mOR262-12 gene, expression was restricted to OSNs that were segregated in the characteristic central patch, although the number of cells varied considerably. Only in one line, the transgene was also expressed in OSNs ectopically positioned outside the patch. The axons of transgene-expressing OSNs co-converged with those expressing the endogenous gene. The transgene was found to be expressed in a mutually exclusive manner and from only one allele indicating that the conserved upstream DNA elements play a critical role in controlling the specific expression pattern of these genes.

Transcription factor modularity in a gene-centered C. elegans core neuronal protein-DNA interaction network

Transcription regulatory networks play a pivotal role in the development, function, and pathology of metazoan organisms. Such networks are comprised of protein-DNA interactions between transcription factors (TFs) and their target genes. An important question pertains to how the architecture of such networks relates to network functionality. Here, we show that a Caenorhabditis elegans core neuronal protein-DNA interaction network is organized into two TF modules. These modules contain TFs that bind to a relatively small number of target genes and are more systems specific than the TF hubs that connect the modules. Each module relates to different functional aspects of the network. One module contains TFs involved in reproduction and target genes that are expressed in neurons as well as in other tissues. The second module is enriched for paired homeodomain TFs and connects to target genes that are often exclusively neuronal. We find that paired homeodomain TFs are specifically expressed in C. elegans and mouse neurons, indicating that the neuronal function of paired homeodomains is evolutionarily conserved. Taken together, we show that a core neuronal C. elegans protein-DNA interaction network possesses TF modules that relate to different functional aspects of the complete network.