Genomic analysis of orthologous mouse and human olfactory receptor loci

Bioinformatics Discovery of Putative Enhancers within Mouse Odorant Receptor Gene Clusters

Olfactory neuronal function depends on the expression and proper regulation of odorant receptor (OR) genes. Previous studies have identified 54 putative intergenic enhancers within or flanking 40 mouse OR clusters. At least 2 of these putative enhancers have been shown to regulate the expression of a small subset of proximal OR genes. In recognition of the large size of the mouse OR gene family (~1400 OR genes distributed across multiple chromosomal loci), it is likely that there remain many additional not-as-yet discovered OR enhancers. We utilized 23 of the previously identified enhancers as a training set (TS) and designed an algorithm that combines a broad range of epigenetic criteria (histone-3-lysine-4 monomethylation, histone-3-lysine-79 trimethylation, histone-3-lysine-27 acetylation, and DNase hypersensitivity) and genetic criteria (cross-species sequence conservation and transcription-factor binding site enrichment) to more broadly search OR gene clusters for additional candidates. We identified 181 new candidate enhancers located at 58 (of 68) mouse OR loci, including 25 new candidates identified by stringent search criteria whose signal strengths are not significantly different from the 23 previously characterized OR enhancers used as the TS. Additionally, we compared OR enhancer versus generic enhancer features in order to evaluate likelihoods that new enhancer candidates specifically function in OR regulation. We found that features distinguishing OR-specific function are significantly more evident for enhancer candidates located within OR clusters as compared with those in flanking regions.

A purely bioinformatic pipeline for the prediction of mammalian odorant receptor gene enhancers

Background

In most mammals, a vast array of genes coding for chemosensory receptors mediates olfaction. Odorant receptor (OR) genes generally constitute the largest multifamily (> 1100 intact members in the mouse). From the whole pool, each olfactory neuron expresses a single OR allele following poorly characterized mechanisms termed OR gene choice. OR genes are found in genomic aggregations known as clusters. Nearby enhancers, named elements, are crucial regulators of OR gene choice. Despite their importance, searching for new elements is burdensome. Other chemosensory receptor genes responsible for smell adhere to expression modalities resembling OR gene choice, and are arranged in genomic clusters — often with chromosomal linkage to OR genes. Still, no elements are known for them.

Results

Here we present an inexpensive framework aimed at predicting elements. We redefine cluster identity by focusing on multiple receptor gene families at once, and exemplify thirty — not necessarily OR-exclusive — novel candidate enhancers.

Conclusions

The pipeline we introduce could guide future in vivo work aimed at discovering/validating new elements. In addition, our study provides an updated and comprehensive classification of all genomic loci responsible for the transduction of olfactory signals in mammals.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-3012-1) contains supplementary material, which is available to authorized users.

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.

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.

Changes in Olfactory Receptor Expression Are Correlated With Odor Exposure During Early Development in the zebrafish (Danio rerio)

We have previously shown that exposure to phenyl ethyl alcohol (PEA) causes an increase in the expression of the transcription factor otx2 in the olfactory epithelium (OE) of juvenile zebrafish, and this change is correlated with the formation of an odor memory of PEA. Here, we show that the changes in otx2 expression are specific to βPEA: exposure to αPEA did not affect otx2 expression. We identified 34 olfactory receptors (ORs) representing 16 families on 4 different chromosomes as candidates for direct regulation of OR expression via Otx2. Subsequent in silico analysis uncovered Hnf3b binding sites closely associated with Otx2 binding sites in the regions flanking the ORs. Analysis by quantitative polymerase chain reaction and RNA-seq of OR expression in developing zebrafish exposed to different isoforms of PEA showed that a subset of ORs containing both Otx2/Hnf3b binding sites were downregulated only in βPEA-exposed juveniles and this change persisted through adult life. Localization of OR expression by in situ hybridization indicates the downregulation occurs at the level of RNA and not the number of cells expressing a given receptor. Finally, analysis of immediate early gene expression in the OE did not reveal changes in c-fos expression in response to either αPEA or βPEA.

The Mouse Solitary Odorant Receptor Gene Promoters as Models for the Study of Odorant Receptor Gene Choice

BACKGROUND

In vertebrates, several anatomical regions located within the nasal cavity mediate olfaction. Among these, the main olfactory epithelium detects most conventional odorants. Olfactory sensory neurons, provided with cilia exposed to the air, detect volatile chemicals via an extremely large family of seven-transmembrane chemoreceptors named odorant receptors. Their genes are expressed in a monogenic and monoallelic fashion: a single allele of a single odorant receptor gene is transcribed in a given mature neuron, through a still uncharacterized molecular mechanism known as odorant receptor gene choice.

AIM

Odorant receptor genes are typically arranged in genomic clusters, but a few are isolated (we call them solitary) from the others within a region broader than 1 Mb upstream and downstream with respect to their transcript's coordinates. The study of clustered genes is problematic, because of redundancy and ambiguities in their regulatory elements: we propose to use the solitary genes as simplified models to understand odorant receptor gene choice.

PROCEDURES

Here we define number and identity of the solitary genes in the mouse genome (C57BL/6J), and assess the conservation of the solitary status in some mammalian orthologs. Furthermore, we locate their putative promoters, predict their homeodomain binding sites (commonly present in the promoters of odorant receptor genes) and compare candidate promoter sequences with those of wild-caught mice. We also provide expression data from histological sections.

RESULTS

In the mouse genome there are eight intact solitary genes: Olfr19 (M12), Olfr49, Olfr266, Olfr267, Olfr370, Olfr371, Olfr466, Olfr1402; five are conserved as solitary in rat. These genes are all expressed in the main olfactory epithelium of three-day-old mice. The C57BL/6J candidate promoter of Olfr370 has considerably varied compared to its wild-type counterpart. Within the putative promoter for Olfr266 a homeodomain binding site is predicted. As a whole, our findings favor Olfr266 as a model gene to investigate odorant receptor gene choice.

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.

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.

Functional interrogation of an odorant receptor locus reveals multiple axes of transcriptional regulation

A transgenic approach in mice allows the functional interrogation of an odorant receptor locus in vivo and reveals characteristics of its monogenic and monoallelic expression.

The odorant receptor (OR) genes constitute the largest mammalian gene family and are expressed in a monogenic and monoallelic fashion, through an unknown mechanism that likely exploits positive and negative regulation. We devised a genetic strategy in mice to examine OR selection by determining the transcriptional activity of an exogenous promoter homologously integrated into an OR locus. Using the tetracycline-dependent transactivator responsive promoter (tet_o), we observed that the OR locus imposes spatial and temporal constraints on tet_o-driven transcription. Conditional expression experiments reveal a developmental change in the permissiveness of the locus. Further, expression of an OR transgene that suppresses endogenous ORs similarly represses the OR-integrated tet_o. Neurons homozygous for the tet_o-modified allele demonstrate predominantly monoallelic expression, despite their potential to express both copies. These data reveal multiple axes of regulation, and support a model of initiation of OR choice limited by nonpermissive chromatin and maintained by repression of nonselected alleles.

Odorant receptor (OR) gene choice is a paradigmatic example of transcriptional regulation in which each olfactory sensory neuron selects a single OR from a repertoire of over 1,000 genes. Two mechanistic models of OR choice have been proposed. One postulates the existence of a specialized transcriptional machinery that selects just one OR allele, while a second, kinetic model proposes that OR chromatin is intrinsically nonpermissive, such that inefficient activation during a critical window of time restricts expression to a single OR allele. Here, we used a transgenic approach in mice in which we inserted a conditionally regulated exogenous promoter into an OR locus by homologous recombination in embryonic stem cells. The resulting novel mouse lines allowed the functional interrogation of the OR locus in vivo during development of the olfactory epithelium, enabling us to directly test models of OR choice. Using this experimental strategy we found that OR loci are indeed slow to activate and that the subsequent phenomenon of spatial restriction of OR expression is accomplished by repression. We also observed a developmental shutdown of OR loci concomitant with expression of the OR repertoire. Together, these experiments provide prima facie evidence for a kinetic model of initiation of OR gene choice, coupled with repression of nonselected OR alleles.

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.

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.

Gene switching and odor induced activity shape expression of the OR37 family of olfactory receptor genes

Olfactory sensory neurons (OSNs) which express distinct odorant receptor (OR) genes are spatially arranged within the mouse olfactory epithelium. Towards an understanding of the mechanisms which determine these patterns, representative OR genes which are typically expressed in the unique central patch of the epithelium were investigated. Inside the patch, numerous OSNs which initially selected a representative gene from this OR group finally expressed another gene from the group, indicating that OSNs inside the patch 'switch' between these genes. If an OSN successively chose genes from the same OR gene cluster, these originated from the same parental chromosome. A deletion of the olfactory cyclic nucleotide-gated ion channel altered the distribution pattern of distinct OSN populations; they were no longer located exclusively inside the patch. Together, the results indicate that OSNs inside the patch initially sample several OR genes for expression; for their correct patterning in the OE, odor-induced activity appears to play a critical role.

Asymmetric and non-uniform evolution of recently duplicated human genes

BACKGROUND

Gene duplications are a source of new genes and protein functions. The innovative role of duplication events makes families of paralogous genes an interesting target for studies in evolutionary biology. Here we study global trends in the evolution of human genes that resulted from recent duplications.

RESULTS

The pressure of negative selection is weaker during a short time immediately after a duplication event. Roughly one fifth of genes in paralogous gene families are evolving asymmetrically: one of the proteins encoded by two closest paralogs accumulates amino acid substitutions significantly faster than its partner. This asymmetry cannot be explained by differences in gene expression levels. In asymmetric gene pairs the number of deleterious mutations is increased in one copy, while decreased in the other copy as compared to genes constituting non-asymmetrically evolving pairs. The asymmetry in the rate of synonymous substitutions is much weaker and not significant.

CONCLUSIONS

The increase of negative selection pressure over time after a duplication event seems to be a major trend in the evolution of human paralogous gene families. The observed asymmetry in the evolution of paralogous genes shows that in many cases one of two gene copies remains practically unchanged, while the other accumulates functional mutations. This supports the hypothesis that slowly evolving gene copies preserve their original functions, while fast evolving copies obtain new specificities or functions.

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.

[Functional analysis of retroviral endogenous inserts in the human genome evolution]

Retroelements, mobile elements produced in DNA by reverse transcription, comprise about 40% of the human genome. A small part of these elements appeared in the genome quite recently after the divergence of humans and chimpanzees had occurred. Evolutionarily young retroelements are represented by the members of four groups, SVA, Alu, L1, and the endogenous HERV-K (HML-2) virus. These retroelements could play a functional role in the course of the molecular evolution of human DNA. We comprehensively studied the contribution of human-specific endogenous viruses (hsERV) to the structural modifications and regulation of the human genome. We found that hsERV presented in 134 copies occupied about 330 000 bp of human DNA. They added to genomic sequences the copies of 50 functional retroviral genes as well as 134 potential promoters and enhancers, 50% of which are located in the regions adjacent to known genes, and 22% in gene introns. At least 67% of these elements are human-specific promoters in vivo. hsERV viruses regulate the activity of known protein-encoding genes by means of RNA interference, function as enhancers, and provide new polyadenylation signals for mRNA.

Mammalian olfactory receptors

Perception of chemical stimuli from the environment is essential to most animals; accordingly, they are equipped with a complex olfactory system capable of receiving a nearly unlimited number of odorous substances and pheromones. This enormous task is accomplished by olfactory sensory neurons (OSNs) arranged in several chemosensory compartments in the nose. The sensitive and selective responsiveness of OSNs to odorous molecules and pheromones is based on distinct receptors in their chemosensory membrane; consequently, olfactory receptors play a key role for a reliable recognition and an accurate processing of chemosensory information. They are therefore considered as key elements for an understanding of the principles and mechanisms underlying the sense of smell. The repertoire of olfactory receptors in mammals encompasses hundreds of different receptor types which are highly diverse and expressed in distinct subcompartments of the nose. Accordingly, they are categorized into several receptor families, including odorant receptors (ORs), vomeronasal receptors (V1Rs and V2Rs), trace amine-associated receptors (TAARs), formyl peptide receptors (FPRs), and the membrane guanylyl cyclase GC-D. This large and complex receptor repertoire is the basis for the enormous chemosensory capacity of the olfactory system.

Mapping of class I and class II odorant receptors to glomerular domains by two distinct types of olfactory sensory neurons in the mouse

The repertoire of approximately 1200 odorant receptors (ORs) is mapped onto the array of approximately 1800 glomeruli in the mouse olfactory bulb (OB). The spatial organization of this array is influenced by the ORs. Here we show that glomerular mapping to broad domains in the dorsal OB is determined by two types of olfactory sensory neurons (OSNs), which reside in the dorsal olfactory epithelium. The OSN types express either class I or class II OR genes. Axons from the two OSN types segregate already within the olfactory nerve and form distinct domains of glomeruli in the OB. These class-specific anatomical domains correlate with known functional odorant response domains. However, axonal segregation and domain formation are not determined by the class of the expressed OR protein. Thus, the two OSN types are determinants of axonal wiring, operate at a higher level than ORs, and contribute to the functional organization of the glomerular array.

Genome-wide patterns of gene flow across a house mouse hybrid zone

Hybrid zones between closely related species or subspecies provide useful settings for studying the genetic architecture of speciation. Using markers distributed throughout the mouse genome, we use a hybrid zone between two recently diverged species of house mice (Mus musculus and Mus domesticus) as a natural mapping experiment to identify genomic regions that may be involved in reproductive isolation. Using cline analysis we document a nearly 50-fold variation in level of introgression among markers. Some markers have extremely narrow cline widths; these genomic regions may contribute to reproductive isolation. Biological processes associated with these narrow clines include physiological and immune responses to the environment as well as physiological and behavioral aspects of reproduction. Other autosomal markers exhibit asymmetrically broad clines, usually with high frequencies of M. domesticus alleles on the M. musculus side of the hybrid zone. These markers identify genome regions likely housing genes with alleles that are spreading from one species to the other. Biological processes associated with these wide clines include cell signaling, olfaction, and pheromone response. These processes play important roles in survival and reproduction, and associated genes are likely targets of selection. Patterns of linkage disequilibrium in the center of the hybrid zone suggest that isolation may be caused by multiple epistatic interactions between sets of genes. These data highlight the complex genetic architecture underlying speciation even at early stages of divergence and point to some of the biological processes that may govern this architecture.

Co-regulation of a large and rapidly evolving repertoire of odorant receptor genes

The olfactory system meets niche- and species-specific demands by an accelerated evolution of its odorant receptor repertoires. In this review, we describe evolutionary processes that have shaped olfactory and vomeronasal receptor gene families in vertebrate genomes. We emphasize three important periods in the evolution of the olfactory system evident by comparative genomics: the adaptation to land in amphibian ancestors, the decline of olfaction in primates, and the delineation of putative pheromone receptors concurrent with rodent speciation. The rapid evolution of odorant receptor genes, the sheer size of the repertoire, as well as their wide distribution in the genome, presents a developmental challenge: how are these ever-changing odorant receptor repertoires coordinated within the olfactory system? A central organizing principle in olfaction is the specialization of sensory neurons resulting from each sensory neuron expressing only ~one odorant receptor allele. In this review, we also discuss this mutually exclusive expression of odorant receptor genes. We have considered several models to account for co-regulation of odorant receptor repertoires, as well as discussed a new hypothesis that invokes important epigenetic properties of the system.

Patterning the developing and regenerating olfactory system

The olfactory system is a remarkable model for investigating the factors that influence the guidance of sensory axon populations to specific targets in the CNS. Since the initial discovery of the vast odorant receptor (ORs) gene family in rodents and the subsequent finding that these molecules directly influence targeting, several additional olfactory axon guidance cues have been identified. Two of these, ephrins and semaphorins, have well-established functions in patterning axon connections in other systems. In addition, lactosamine-containing glycans are also required for proper targeting and maintenance of olfactory axons, and may also function in other sensory regions. It is now apparent that these and likely other additional molecules are required along with ORs to orchestrate the complex pattern of convergence and divergence that is unique to the olfactory system.

Ancient genomic architecture for mammalian olfactory receptor clusters

A new tool for genome-wide definition of genomic gene clusters conserved in multiple species was applied to olfactory receptors in five mammals, demonstrating that most mammalian olfactory receptor clusters have a common ancestry.

Background

Mammalian olfactory receptor (OR) genes reside in numerous genomic clusters of up to several dozen genes. Whole-genome sequence alignment nets of five mammals allow their comprehensive comparison, aimed at reconstructing the ancestral olfactory subgenome.

Results

We developed a new and general tool for genome-wide definition of genomic gene clusters conserved in multiple species. Syntenic orthologs, defined as gene pairs showing conservation of both genomic location and coding sequence, were subjected to a graph theory algorithm for discovering CLICs (clusters in conservation). When applied to ORs in five mammals, including the marsupial opossum, more than 90% of the OR genes were found within a framework of 48 multi-species CLICs, invoking a general conservation of gene order and composition. A detailed analysis of individual CLICs revealed multiple differences among species, interpretable through species-specific genomic rearrangements and reflecting complex mammalian evolutionary dynamics. One significant instance involves CLIC #1, which lacks a human member, implying the human-specific deletion of an OR cluster, whose mouse counterpart has been tentatively associated with isovaleric acid odorant detection.

Conclusion

The identified multi-species CLICs demonstrate that most of the mammalian OR clusters have a common ancestry, preceding the split between marsupials and placental mammals. However, only two of these CLICs were capable of incorporating chicken OR genes, parsimoniously implying that all other CLICs emerged subsequent to the avian-mammalian divergence.

Les récepteurs des molécules odorantes et le codage olfactif

The first step of olfactory detection involves interactions between odorant molecules and neuronal protein receptors. Odour coding results from the combinatory activation of a set of receptors and rests on their clonal expression and olfactory neurone connexion, which lead to formation of a specific sensory map in the cortex. This system, sufficient to discriminate myriads of odorants with a mere 350 different receptors, allows humans to smell molecules that are not natural (new cooking flavours, synthetic chemicals...). The extreme olfactory genome diversity explains the absence of odour semantics. Olfactory receptors are also involved in cellular chemotaxis.

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

Mouse odorant receptors (ORs) are encoded by >1000 genes dispersed throughout the genome. Each olfactory neuron expresses one single OR gene, while the rest of the genes remain silent. The mechanisms underlying OR gene expression are poorly understood. Here, we investigated if OR genes share common cis-regulatory sequences in their promoter regions. We carried out a comprehensive analysis in which the upstream regions of a large number of OR genes were compared. First, using RLM-RACE, we generated cDNAs containing the complete 5'-untranslated regions (5'-UTRs) for a total number of 198 mouse OR genes. Then, we aligned these cDNA sequences to the mouse genome so that the 5' structure and transcription start sites (TSSs) of the OR genes could be precisely determined. Sequences upstream of the TSSs were retrieved and browsed for common elements. We found DNA sequence motifs that are overrepresented in the promoter regions of the OR genes. Most motifs resemble O/E-like sites and are preferentially localized within 200 bp upstream of the TSSs. Finally, we show that these motifs specifically interact with proteins extracted from nuclei prepared from the olfactory epithelium, but not from brain or liver. Our results show that the OR genes share common promoter elements. The present strategy should provide information on the role played by cis-regulatory sequences in OR gene regulation.

Promoter motifs of olfactory receptor genes expressed in distinct topographic patterns

Novel olfactory receptor-encoding genes that are expressed in olfactory sensory neurons arranged in a clustered pattern in the nasal epithelium, typical of the mOR262 (approved gene symbol Olfr) family, were identified. The genes share sequence motifs upstream of their transcription start sites that are highly related to those previously identified as characteristic of the mOR262 genes, suggesting that these regulatory elements may contribute to governing their unique expression pattern. Promoter analyses of genes encoding class I receptors that are expressed in the dorsal region of the epithelium revealed a different, but again common set of sequence motifs. A prominent feature of the class I gene promoters are multiple O/E-like binding sites, and O/E-type transcription factors that bind to the putative promoter region of class I OR genes were in fact identified. The findings support the concept that common elements in the promoter region of these OR genes may determine their congenic expression pattern in the epithelium.

Regulation of odorant receptors: one allele at a time

The odorant receptors (ORs) make up the largest gene family in mammals. Each olfactory sensory neuron chooses just one OR from the more than 1000 possibilities encoded in the genome and transcribes it from just one allele. This process generates great neuronal diversity and forms the basis for the development and logic of the olfactory circuit between the nose and the brain. The mechanism behind this monoallelic regulation has been the subject of intense speculation and increasing experimental investigation, yet remains enigmatic. Recent genetic experiments have brought the outlines of the process into sharper relief, identifying a feedback mechanism in which the first odorant receptor expressed, generates a signal that stabilizes its choice, thus maintaining singular selection. In the absence of this signal, the olfactory neuron re-enters the selection process and switches to choose an alternate OR. Irreversible genetic changes in the nuclei of olfactory neurons do not accompany OR selection, which must therefore be initiated by an epigenetic process that may involve a stochastic mechanism.

Comparative evolutionary analysis of olfactory receptor gene clusters between humans and mice

Olfactory receptor (OR) genes form the largest multigene family in mammalian genomes. Humans have approximately 800 OR genes, but >50% of them are pseudogenes. By contrast, mice have approximately 1400 OR genes and pseudogenes are approximately 25%. To understand the evolutionary processes that shaped the difference of OR gene families between humans and mice, we studied the genomic locations of all human and mouse OR genes and conducted a detailed phylogenetic analysis using functional genes and pseudogenes. We identified 40 phylogenetic clades with high bootstrap supports, most of which contain both human and mouse genes. Interestingly, a particular clade contains approximately 100 pseudogenes in humans, whereas the numbers of pseudogenes are <20 for most of the mouse clades. We also found that the organization of OR genomic clusters is well conserved between humans and mice in many chromosomal locations. Despite the difference in the numbers of genes, the numbers of large genomic clusters are nearly the same for humans and mice. These observations suggest that the greater OR gene repertoire in mice has been generated mainly by tandem gene duplication within each genomic cluster.

Detection of odorants through the main olfactory epithelium and vomeronasal organ of mice

Previous research has indicated that volatile odorants are detected through the main olfactory epithelium (MOE), whereas pheromones are detected via the vomeronasal organ (VNO). Gene disruption studies have established that olfactory signaling through the MOE is mediated through receptor stimulation of type 3 adenylyl cyclase (AC3). Mice lacking AC3 cannot detect odorants through the MOE. Recently, it was discovered using olfactory-based behavioral assays that AC3 mutant mice can detect some volatile odorants. An analysis of these mutant mice led to the surprising discovery that some odorants are detected through the VNO.

Using shared genomic synteny and shared protein functions to enhance the identification of orthologous gene pairs

MOTIVATION

The identification of orthologous gene pairs is generally based on sequence similarity. Gene pairs that are mutually 'best hits' between the genomes being compared are asserted to be orthologs. Although this method identifies most orthologous gene pairs with high confidence, it will miss a fraction of them, especially genes in duplicated gene families. In addition, the approach depends heavily on the completeness and quality of gene annotation. When the gene sequences are not correctly represented the approach is unlikely to find the correct ortholog. To overcome these limitations, we have developed an approach to identify orthologous gene pairs using shared chromosomal synteny and the annotation of protein function.

RESULTS

Assembled mouse and human genomes were used to identify the regions of conserved synteny between these genomes. 'Syntenic anchors' are conserved non-repetitive locations between mouse and human genomes. Using these anchors, we identified blocks of sequences that contain consistently ordered anchors between the two genomes (syntenic blocks). The synteny information has been used to help us identify orthologous gene pairs between mouse and human genomes. The approach combines the mutual selection of the best tBlastX hits between human and mouse transcripts, and inferring gene orthologous relationships based on sharing syntenic anchors, collocating in the same syntenic blocks and sharing the same annotated protein function. Using this approach, we were able to find 19,357 orthologous gene pairs between human and mouse genomes, a 20% increase in the number of orthologs identified by conventional approaches.

A contextual model for axonal sorting into glomeruli in the mouse olfactory system

No models fully account for how odorant receptors (ORs) function in the guidance of axons of olfactory sensory neurons (OSNs) to glomeruli in the olfactory bulb. Here, we use gene targeting in mice to demonstrate that the OR amino acid sequence imparts OSN axons with an identity that allows them to coalesce into glomeruli. Replacements between the coding regions of the M71 and M72 OR genes reroute axons to their respective glomeruli. A series of M71-M72 hybrid ORs uncover a spectrum of glomerular phenotypes, leading to the concept that the identity of OSN axons is revealed depending on what other axons are present. Naturally occurring amino acid polymorphisms in other ORs also produce distinct axonal identities. These critical amino acid residues are distributed throughout the protein and reside predominantly within transmembrane domains. We propose a contextual model for axon guidance in which ORs mediate homotypic interactions between like axons.

The LIM-homeodomain protein Lhx2 is required for complete development of mouse olfactory sensory neurons

In mice, approximately 1,000 odorant receptor (OR) genes are expressed in olfactory sensory neurons (OSNs). Homeodomain sites can be recognized in the promoter and upstream regions of several OR genes. Here, using the yeast one-hybrid system and electrophoretic mobility shift assay, we report that Lhx2, a LIM-homeodomain protein, binds to the homeodomain site in the mouse M71 OR promoter region. In Lhx2-deficient mice, the morphology of the olfactory epithelium is grossly normal. However, expression of OMP is abolished and that of GAP43 is severely reduced, indicating that no mature and few immature OSNs are produced. M71 and other OR genes also are not expressed. OSN development appears to be arrested between the terminal differentiation into neurons and the transition to immature neurons. Thus, Lhx2 is required for complete development of OSNs in mice.

Evolution and comparative genomics of odorant- and pheromone-associated genes in rodents

Chemical cues influence a range of behavioral responses in rodents. The involvement of protein odorants and odorant receptors in mediating reproductive behavior, foraging, and predator avoidance suggests that their genes may have been subject to adaptive evolution. We have estimated the consequences of selection on rodent pheromones, their receptors, and olfactory receptors. These families were chosen on the basis of multiple gene duplications since the common ancestor of rat and mouse. For each family, codons were identified that are likely to have been subject to adaptive evolution. The majority of such sites are situated on the solvent-accessible surfaces of putative pheromones and the lumenal portions of their likely receptors. We predict that these contribute to physicochemical and functional diversity within pheromone-receptor interaction sites.

Species specificity in rodent pheromone receptor repertoires

The mouse V1R putative pheromone receptor gene family consists of at least 137 intact genes clustered at multiple chromosomal locations in the genome. Species-specific pheromone receptor repertoires may partly explain species-specific social behavior. We conducted a genomic analysis of an orthologous pair of mouse and rat V1R gene clusters to test for species specificity in rodent pheromone systems. Mouse and rat have lineage-specific V1R repertoires in each of three major subfamilies at these loci as a result of postspeciation duplications, gene loss, and gene conversions. The onset of this diversification roughly coincides with a wave of Line1 (L1) retrotranspositions into the two loci. We propose that L1 activity has facilitated postspeciation V1R duplications and gene conversions. In addition, we find extensive homology among putative V1R promoter regions in both species. We propose a regulatory model in which promoter homogenization could ensure that V1R genes are equally competitive for a limiting transcriptional structure to account for mutually exclusive V1R expression in vomeronasal neurons.

After the holy grail: establishing a molecular basis for Mammalian olfaction

The quest to identify mammalian odorant receptors was a triumph of molecular biology. The characterization of these molecules has provided extraordinary insight into the strategy used by one neuronal system to organize sensory structures and code complex information. The odorant receptor genes have also served as powerful tools in understanding genomic organization and gene regulation.

Functional bias and spatial organization of genes in mutational hot and cold regions in the human genome

The neutral mutation rate is known to vary widely along human chromosomes, leading to mutational hot and cold regions. We provide evidence that categories of functionally related genes reside preferentially in mutationally hot or cold regions, the size of which we have measured. Genes in hot regions are biased toward extracellular communication (surface receptors, cell adhesion, immune response, etc.), while those in cold regions are biased toward essential cellular processes (gene regulation, RNA processing, protein modification, etc.). From a selective perspective, this organization of genes could minimize the mutational load on genes that need to be conserved and allow fast evolution for genes that must frequently adapt. We also analyze the effect of gene duplication and chromosomal recombination, which contribute significantly to these biases for certain categories of hot genes. Overall, our results show that genes are located nonrandomly with respect to hot and cold regions, offering the possibility that selection acts at the level of gene location in the human genome.

Functionally related genes tend to be found in regions of the genome with high or low mutation rates, which implies that natural selection can also operate at the level of gene location

The clustered olfactory receptor gene family 262: genomic organization, promotor elements, and interacting transcription factors

For six mouse olfactory receptor genes from family 262 which are expressed in clustered populations of olfactory sensory neurons, the genomic as well as cDNA structures were deciphered. All genes contained several exons which in some cases were alternatively spliced. Immediately upstream of the transcription start sites, sequence motif blocks were identified that are highly conserved among olfactory receptor (OR) genes which are expressed in clustered neuronal populations. By means of electrophoretic mobility shift assays, it was demonstrated that segments of the motif block region interact with proteins extracted from nuclear fractions of the olfactory epithelium. Yeast one-hybrid screenings of an olfactory cDNA library led to the identification of a set of transcription factors that specifically bind to particular elements of the motif block region. The identified factors can be categorized into two types: One group is known to be involved in transcriptional initiation, and the second group represents factors involved in pattern formations. The identified components may contribute to govern the precise topographic expression pattern of olfactory receptor genes.

Gene structure prediction in syntenic DNA segments

The accurate prediction of higher eukaryotic gene structures and regulatory elements directly from genomic sequences is an important early step in the understanding of newly assembled contigs and finished genomes. As more new genomes are sequenced, comparative approaches are becoming increasingly practical and valuable for predicting genes and regulatory elements. We demonstrate the effectiveness of a comparative method called pattern filtering; it utilizes synteny between two or more genomic segments for the annotation of genomic sequences. Pattern filtering optimally detects the signatures of conserved functional elements despite the stochastic noise inherent in evolutionary processes, allowing more accurate annotation of gene models. We anticipate that pattern filtering will facilitate sequence annotation and the discovery of new functional elements by the genetics and genomics communities.

Microarray-based discovery of highly expressed olfactory mucosal genes: potential roles in the various functions of the olfactory system

We sought to gain a global view of tissue-specific gene expression in the olfactory mucosa (OM), the major site of neurogenesis and neuroregeneration in adult vertebrates, by examination of its overexpressed genes relative to that in 81 other developing and adult mouse tissues. We used a combination of statistical and fold-difference criteria to identify the top 269 cloned cDNAs from an array of 8,734 mouse cDNA elements on the Incyte Mouse GEM1 array. These clones, representing known and poorly characterized gene transcripts, were grouped according to their relative expression patterns across the other tissues and then further examined with respect to gene ontology categories. Approximately one-third of the 269 genes were also highly expressed in developing and/or adult central nervous system tissues. Several of these have been suggested or demonstrated to play roles in neurogenesis, neuronal differentiation, and/or neuronal migration, further suggesting that many of the unknown genes that share this expression pattern may play similar roles. Highly OM-specific genes included a palate, lung, and nasal epithelium carcinoma-associated gene (Plunc); sphingosine phosphate lyase (Sgpl1), and paraoxonase 1 (Pon1). Cell-type-specific expression within OM was established using in situ hybridization for several representative expression pattern clusters. Using the ENSEMBL-assembled mouse genome and comparative genomics analyses to the human genome, we assigned many of the unknown expressed sequence tags (ESTs) and poorly characterized genes to either novel or known gene products and provided predictive classification. Further exploration of this database will provide additional insights into genes and pathways critical for olfactory neurogenesis, neuronal differentiation, olfaction, and mucosal defense.

Co-duplication of olfactory receptor and MHC class I genes in the mouse major histocompatibility complex

We report the 897 kb sequence of a cluster of olfactory receptor (OR) genes located at the distal end of the major histocompatibility complex (MHC) class I region on mouse chromosome 17 of strain 129/SvJ (H2bc). With additional information from the mouse genome draft sequence, we identified 59 OR loci (approximately 20% pseudogenes) in contrast to only 25 OR loci (approximately 50% pseudogenes) in the corresponding centromeric OR cluster that is part of the 'extended MHC class I region' on human chromosome 6. Comparative analysis leads to three major observations: (i) most of the OR subfamilies have evolved independently in the two species, expanding more in the mouse, and resulting in co-orthologs--subfamilies of highly similar paralogs that keep orthologous relationships with their human counterparts; (ii) three of the mouse OR subfamilies have no orthologs in humans; and (iii) MHC class I loci are interspersed in the OR cluster in mouse but not in human, and were subjected to co-duplication with OR genes. Screening of our sequence against the available sequences of other strains/haplotypes revealed that most of the OR loci are polymorphic and that the number of OR loci may vary among strains/haplotypes. Our findings that MHC-linked OR loci share duplication with MHC class I loci, have duplicated extensively and are polymorphic revives questions about potential reciprocal influences acting on the dynamics and evolution of the H2 region and the H2-linked OR loci.

The potential role of gene duplications in the evolution of imprinting mechanisms

Using the completed genomic sequences of mouse and human we performed a comparative analyses of imprinted genes and gene clusters. For many imprinted genes we could detect imprinted as well as non-imprinted paralogues. The inter- and intrachromosomal similarities between paralogues and their linkage to imprinting clusters suggests that imprinted genes were dispersed throughout the genome by gene duplications as well as translocation and transposition events. Our findings indicate that imprinting clusters may have been linked together on one (or a few) ancestral pre-imprinted chromosome(s), arguing for a common mechanistic origin of imprinting control. Imprinting may originally have evolved on a simple basis of dosage compensation required for some duplicated genes (chromosomes) followed by selection of sex-biased expression control.

Odorant receptor expressed sequence tags demonstrate olfactory expression of over 400 genes, extensive alternate splicing and unequal expression levels

BACKGROUND

The olfactory receptor gene family is one of the largest in the mammalian genome. Previous computational analyses have identified approximately 1,500 mouse olfactory receptors, but experimental evidence confirming olfactory function is available for very few olfactory receptors. We therefore screened a mouse olfactory epithelium cDNA library to obtain olfactory receptor expressed sequence tags, providing evidence of olfactory function for many additional olfactory receptors, as well as identifying gene structure and putative promoter regions.

RESULTS

We identified more than 1,200 odorant receptor cDNAs representing more than 400 genes. Using real-time PCR to confirm expression level differences suggested by our screen, we find that transcript levels in the olfactory epithelium can differ between olfactory receptors by up to 300-fold. Differences for one gene pair are apparently due to both unequal numbers of expressing cells and unequal transcript levels per expressing cell. At least two-thirds of olfactory receptors exhibit multiple transcriptional variants, with alternative isoforms of both 5' and 3' untranslated regions. Some transcripts (5%) utilize splice sites within the coding region, contrary to the stereotyped olfactory receptor gene structure. Most atypical transcripts encode nonfunctional olfactory receptors, but can occasionally increase receptor diversity.

CONCLUSIONS

Our cDNA collection confirms olfactory function of over one-third of the intact mouse olfactory receptors. Most of these genes were previously annotated as olfactory receptors based solely on sequence similarity. Our finding that different olfactory receptors have different expression levels is intriguing given the one-neuron, one-gene expression regime of olfactory receptors. We provide 5' untranslated region sequences and candidate promoter regions for more than 300 olfactory receptors, valuable resources for computational regulatory motif searches and for designing olfactory receptor microarrays and other experimental probes.