Mammals sense odors through the gene family of olfactory receptors (ORs). Despite the enormous number of OR genes (∼1,400 in mouse), each olfactory sensory neuron expresses one, and only one, of them. In neurobiology, it remains a long-standing mystery how this singularity can be achieved despite intrinsic stochasticity of gene expression. Recent experiments showed an epigenetic mechanism for maintaining singular OR expression: Once any ORs are activated, their expression inhibits further OR activation by down-regulating a histone demethylase Lsd1 (also known as Aof2 or Kdm1a), an enzyme required for the removal of the repressive histone marker H3K9me3 on OR genes. However, it remains unclear at a quantitative level how singularity can be initiated in the first place. In particular, does a simple activation/ feedback scheme suffice to generate singularity? Here we show theoretically that rare events of histone demethylation can indeed produce robust singularity by separating two timescales: slow OR activation by stepwise H3K9me3 demethylation, and fast feedback to turn off Lsd1. Given a typical 1-h response of transcriptional feedback, to achieve the observed extent of singularity (only 2% of neurons express more than one ORs), we predict that OR activation must be as slow as 5-10 d-a timescale compatible with experiments. Our model further suggests H3K9me3-to-H3K9me2 demethylation as an additional rate-limiting step responsible for OR singularity. Our conclusions may be generally applicable to other systems where monoallelic expression is desired, and provide guidelines for the design of a synthetic system of singular expression.neurogenesis | stochastic gene expression | histone modification | kinetics modeling | negative feedback I n mammals, the ability to sense odors relies on the singular expression of the olfactory receptor (OR) gene family. Despite the enormous family size (∼1,400 genes in the mouse genome), each sensory neuron only expresses one single allele of OR (1-4). It has long been thought that this singularity stems from a negative feedback in which the expression of one OR specifically silences all other ORs (5-7). However, this hypothesis has led to the question of how the one, and only one, active OR can escape its own silencing effect. In other words, it was unclear how such silencing can be biologically feasible.A series of recent studies provoked a different mechanism for OR singularity in light of epigenetic regulation. Contrary to previous belief, Magklara et al. (8) found that silencing of OR genes precedes OR expression: The histone marker H3K9me2 on OR genes become methylated into H3K9me3 as early as in the stage of neuronal progenitors, which do not yet express any ORs. H3K9me2 is a marker for localization of inactive genes into facultative heterochromatin; its further methylation into H3K9me3, however, marks OR genes for constitutive heterochromatinnuclear compartments that normally contain pericentromeric and telomeric regions. With this marker, all OR genes are deeply repressed i...