ORs and other AT-rich gene families frequently associate with the nuclear lamina (Peric-Hupkes et al., 2010
). However, our DNA FISH analysis with individual BAC probes failed to reveal a significant distribution of OR loci toward the nuclear periphery of OSNs (Lomvar-das et al., 2006
). To obtain a comprehensive view of the distribution of OR loci in OSN nuclei, we sought to generate a DNA FISH probe that would allow the simultaneous detection of most OR loci. First, because OR clusters reside in extremely AT-rich isochores (Clowney et al., 2011
; Glusman et al., 2001
), we digested genomic DNA with restriction enzymes that recognize AT-rich sequences and collected DNA fractions with significant enrichment for ORs. Next, these were amplified and subjected to a second round of purification by sequence capture on a custom tiling array covering OR clusters () (Albert et al., 2007
). This high-density array contains oligonucleotides against the unique sequences within the 46 OR genomic clusters, spanning a total region of 40 MB. Two rounds of capture, elution, and amplification produced a DNA library highly enriched for OR sequences. Quantitative PCR analysis (qPCR) of the final amplicon detects only sequences from OR clusters, suggesting the elimination of unique, non-OR DNA ().
Visualizing the Nuclear Distribution of OR Loci
To further examine the composition of this DNA library, we analyzed its contents by whole-genome microarray hybridization, using a tiling array covering mouse chromosomes 1 to 4. This analysis demonstrates the sensitivity and selectivity of our purification strategy: the probe detects 340 of 346 OR genes located on these 4 chromosomes, 40 of ~80 non-OR genes located cis
to OR clusters (and included on the capturing array), and 6 of ~5,000 non-OR genes (FDR < 0.05, 98.2% sensitivity, 98.4% OR cluster specificity, p ≤10−72
) ( and Figures S1B–S1D
and Table S1
OSN-Specific Aggregation of OR Genes
We used this “panOR” library as a probe for DNA FISH experiments on sections of the MOE. Although there are 92 OR clusters in the diploid nucleus, the panOR probe detects an average of ~5 large foci in OSNs (). This unexpected distribution is specific for OSNs: OR distribution in other cell populations represented in MOE sections (undifferentiated basal cells and sustentacular cells) is diffuse and more consistent with a random arrangement of the 92 OR clusters or ~2,800 alleles. Quantification of the distribution of the DNA FISH signal in the three cell types of the MOE across the same sections in the same experiments supports this conclusion (Figure S1E–S1F
): high-intensity pixels (above 120 in the 8 bit range of 0–255) were found only in OSNs and not in sustentacular or basal cells. To quantify the distribution of panOR signal, we calculated standard deviation of signal intensity across nuclear space. Average standard deviation in OSNs is 42.3, indicating spotty signal distribution, and is 9.3 or 11.3 in basal and sustentacular cells, indicating smoother distribution (n > 100 for each cell type). Finally, DNA FISH with this probe in other neuronal types demonstrates a diffuse distribution of OR loci (data not shown and Figure S2D
), arguing for an OSN-specific nuclear pattern.
The focal nature of the panOR DNA FISH signal suggests that OR alleles from different OR clusters merge in distinct nuclear regions during OSN differentiation. To test this, we pooled 10 OR- or 12 non-OR-BAC probes and performed two-color DNA FISH with the panOR probe. There was extensive colocalization between the panOR probe and OR BAC probes (Pearson’s coefficient r = 0.637, Mander’s coefficient of BAC signal colocalizing with panOR M1 = 0.835, n > 100) and little colocalization between the panOR probe and the non-OR BACs (r = 0.187, M1 = 0.109, n > 100) ( and Table S2
), suggesting selectivity for OR loci in the composition of these aggregates (Figure S1I
). Though the panOR probe includes most OR loci, lack of complete overlap between the panOR and the individual OR BAC probes was expected. The panOR probe is 200-fold more complex than each BAC, and it is outcompeted for binding at ORs targeted by a BAC. Thus, BAC signals colocalized with panOR signal represent OR alleles surrounded by other OR loci labeled by the panOR probe at distances below the optical resolution of confocal microscopy.
The combined OR BACs produce fewer DNA FISH spots in the OSNs (3.94 spots/nucleus/Z stack, n = 38) than in sustentacular (9.1 spots, n = 38) or basal cells (8.52 spots, n = 30), providing an independent verification for the extensive aggregation of these loci: they are optically indiscrete significantly more often in OSNs than in basal and sustentacular cells. Non-OR BAC probes did not appear more aggregated in the OSNs (10.08 spots in OSNs, 6.4 in sustentacular, and 7.1 in basal cells, n = 30 for each cell type).
To explore the contribution of intra- and interchromosomal interactions to the formation of OR foci and the colocalization of ORs, we used two additional pools of OR BACs, one containing seven BACs targeting three clusters on chromosome 2 and the other containing eight BACs, each targeting a cluster from a different chromosome. These pools, when combined with panOR probe, revealed two layers of organization in OSNs: alleles within the same cluster coalesce into optically indiscrete signals, whereas clusters from different chromosomes generate distinguishable signals inside the same panOR focus (). The BACs from the same chromosome produce 5.9 dots in sustentacular cells and 2.3 dots in OSNs (n = 30 for each), whereas BACs from different chromosomes produce equal numbers of dots in both cell types. However, multiple OR BAC dots from different chromosomes were seen in 50% of the panOR foci, and more than two dots per aggregate in 29% of the cells (n = 50) (). Moreover, maternal and paternal alleles of the same OR cluster reside in the same OR aggregate in ~6% of the tested OSNs (Figure S1H
). Finally, panOR foci do not colocalize with large repeat classes, pericentromeric heterochromatin (PH), or other multigene families (Figures S2A–S2C
and data not shown), suggesting that the aggregation of OR clusters produces distinct and selective OR gene territories.
Spatial Segregation between Active and Silent OR Alleles
To reveal the epigenetic signature of OR foci, we combined DNA FISH analysis with immunofluorescence (IF) against the heterochromatic marks found on ORs (Magklara et al., 2011
) or hetero-chromatin-binding protein 1 β (HP1β), the only heterochromatic HP1 member expressed in OSNs (data not shown). This analysis reveals overlap between the OR foci, H3K9me3, H4K20me3, and HP1β ( and Table S2
), but not with Pol II (Figure S2E
), consistent with a heterochromatic nature of these aggregates. This colocalization is differentiation dependent and cell type specific; we do not detect overlap between the two signals in basal cells in the MOE or in retinal neurons (Figure S2D
OR Foci Are Heterochromatic Aggregates from which the Active Allele Escapes
We then performed nascent RNA FISH on sections of the MOE using intronic probes against OR genes MOR28, M50, M71, and P2 combined with IF for H3K9me3, H4K20me3, or HP1β. In contrast with the bulk of the panOR signal, active OR alleles have little overlap with any of the three heterochromatic marks (, S2F, and S2G
and Table S2
). We also combined nascent RNA FISH with IF for Pol II or H3K27-Acetyl and H4K20me3 or HP1β ( and Table S2
). These experiments corroborate that the active OR allele is spatially segregated from the silent ORs and resides in euchromatic territory.
LBR Organizes the Topology of OSN Nuclei
It is intriguing that most OR genes and PH are located near the center of OSN nuclei instead of being distributed toward the nuclear envelope (Figure S1A
). This “inside-out” nuclear morphology is reminiscent of the nuclear architecture reported in homozygous Ichthyosis mice, a spontaneous LBR loss-of-function mutant (Goldowitz and Mullen, 1982
). LBR is a nuclear envelope protein that interacts with HP1 and heterochromatin (Hoffmann et al., 2002
; Okada et al., 2005
; Pyrpasopoulou et al., 1996
). RNA-seq revealed a continuous reduction in LBR mRNA levels during differentiation from HBCs to OSNs, and IF confirmed that whereas LBR is present in the nuclear envelope of basal and sustentacular cells, it is absent in the neuronal lineage of the MOE (, and S3A
LBR Regulates Nuclear Topology in the MOE
PanOR DNA FISH on MOE sections from the Ichthyosis mice revealed no changes in OR aggregation in OSNs, which already lack LBR. However, nuclear architecture and OR organization of Ichthyotic basal and sustentacular cells approach that of wild-type OSNs (, and S3B
). PH forms large, centrally located foci in both cell types, and ORs form aggregates at the periphery of the pericentromeric foci. According to the pooled BAC assay in the Ichthyosis mouse, the number of DNA FISH spots is uniform among the three cell types, and the basal and sustentacular cells have similar numbers of DNA FISH spots to control OSNs (), supporting a role for LBR downregulation in OR aggregation. Because ectopic OR aggregation occurs in two cell types that do not express ORs and likely do not contain the transcription factors responsible for OR activation, an effect of this mutation on OR expression and OSN targeting is unlikely and was not detected (data not shown and Figure S3C
Thus, we sought to perform the opposite experiment: to restore LBR expression to OSNs instead of removing LBR from cells that do not express ORs. We generated a tetO LBR-IRES-GFP transgenic mouse that we crossed to OMPIRES-tTA mice to achieve expression of LBR in OSNs. One transgenic line expresses the transgene in a significant proportion of OSNs (). Like endogenous LBR, transgenic LBR is restricted to the nuclear envelope without diffusing in the nucleoplasm ().
We used this transgenic line to analyze the effects of ectopic LBR expression on the nuclear morphology of mOSNs. DAPI staining becomes less intense, and PH is moved toward the nuclear periphery of LBR+ OSNs ( and S3D
). OSNs in these sections that do not express the transgene have morphology similar to wild-type nuclei (Figure S3D
). IF shows HP1β recruitment to the nuclear envelope in LBR+ OSNs, whereas centrally shifted euchromatin occupies most of the nucleus (, S3E, and S3F
). Thus, ectopic LBR expression in a postmitotic cell is sufficient to reverse the “inside-out” arrangement and to recruit PH to the nuclear periphery.
Ectopic LBR Expression Decondenses OSN Heterochromatin
IF does not provide information about the structural and biophysical changes occurring in OSN chromatin upon ectopic LBR expression. To obtain this information, we imaged control and LBR+ OSNs with soft X-ray tomography (SXT), a high-resolution imaging method that is applied to fully hydrated, unfixed, and unstained cells and measures carbon and nitrogen concentration in biological samples (McDermott et al., 2009
). Orthoslices (computer-generated sections) and three-dimensional (3D) reconstructions of SXT imaging of control OSNs reveal that the more condensed (darker) chromatin is located at the center of the nucleus, in agreement with the morphology seen by IF ( and Movie S1
). We also detect extremely condensed structures at the periphery of this PH core that are specific for this cell type; only sperm nuclei have chromatin particles with higher compaction values (data not shown). Although the arrangement of these dark foci is similar to the arrangement of the OR foci around the PH core of the OSN nucleus, DNA FISH or IF are incompatible with SXT, and therefore it is impossible to prove directly that these are the same structures.
Soft X-Ray Tomography of OSNs Demonstrates Chromatin Decompaction and Nuclear Reorganization upon LBR Expression
SXT imaging of LBR+ OSNs shows the relocation of the most condensed chromatin toward the nuclear membrane ( and Movie S2
). Moreover, LBR expression increases the nuclear volume from 105u3
and induces the folding of the nuclear membrane and an overall change in the nuclear shape ( and Movies S3
). Overall, chromatin decondensation induced by LBR expression in OSNs is quantitatively described by measurements of the linear absorption coefficient (LAC) (McDermott et al., 2009
) of control and LBR+ OSNs (). This measurement, which depicts the concentration of organic material per voxel, corroborates the loss of the densest foci upon LBR expression. Thus, if condensed regions correspond to OR foci, ectopic LBR expression should cause decompaction of OR heterochromatin. DNaseI sensitivity experiments (Magklara et al., 2011
) in nuclei from fluorescence-activated cell-sorted (FAC-sorted) control or LBR+ OSNs confirms a significant decompaction of OR and peri-centromeric heterochromatin upon LBR expression ().
Ectopic LBR Expression Disrupts OR Aggregation
The spatial reorganization of HP1β, the elimination of the dark foci detected by SXT, and the increase in DNase sensitivity of OR chromatin suggest that ectopic LBR expression disrupts the aggregation of OR loci. To test this, we performed DNA FISH with the panOR probe in sections of LBR-expressing trans-genic mice. Low-magnification images show significant effects of ectopic LBR expression on the distribution of OR loci. In the apical LBR+ neuronal layer, the intense OR foci dissolve; in contrast, immature OSNs and progenitors that do not yet express the transgene but have already downregulated the endogenous LBR retain a focal OR arrangement ( and S4A
To investigate whether altering the tertiary organization of OR loci affects the epigenetic characteristics of these genes, we examined association of OR genes with H3K9me3, H4K20me3, and HP1β in LBR+ OSNs. H3K9me3 and H4K20me3 remained enriched on OR loci upon LBR expression by native ChIP-qPCR assays on FAC-sorted OSNs () and FISH-IF (Figures S4B and S4C
and Table S2
). In contrast, association of OR loci with HP1β was reduced as measured by FISH-IF ( and Table S2
). Reduction in overlap between H4K20me3 and HP1β was also observed in the LBR+ OSNs. Thus, despite retaining heterochromatic histone marks, OR loci lose their aggregated arrangement and their nonhistone heterochromatic coat upon LBR expression, which is consistent with the increased DNase sensitivity.
In wild-type OSNs, active OR alleles interact with the H enhancer. To test whether LBR expression also abrogates inter-chromosomal interactions between the active allele and H, we performed circularized chromosome conformation capture (4C) using inverse H PCR primers as previously described (Lomvar-das et al., 2006
) on LBR-expressing or control MOEs. To increase the proportion of LBR-expressing cells in this mixed population, we combined two tTA drivers (OMP-IRES-tTA and Gγ8 tTA). The enrichment of various OR sequences in this 4C library was assayed by qPCR. LBR expression in OSNs results in the loss of most H-OR associations. In LBR transgenics, H retains only its interaction with the linked OR MOR28 (Olfr1507
) located 75 Kb downstream (). Therefore, ectopic LBR expression in OSNs not only prevents heterochromatic OR aggregation, but also disrupts the interaction between the H enhancer and unlinked ORs.
LBR Expression Inhibits OR Transcription
IF and RNA FISH experiments in MOE sections from control and LBR-expressing mice revealed a 3-fold reduction in the numbers of neurons expressing particular ORs in the transgenic mice (). Importantly, most neurons that retain high-level OR expression do not express transgenic LBR (data not shown and ). For more quantitative measure of the effects of LBR in OR expression, we used FAC-sorting to isolate control or LBR+ OSNs and performed quantitative, reverse-transcriptase PCR (qRT-PCR). This analysis supports that LBR expression has significant inhibitory effects on OR expression (). Similarly, whole-mount X-gal staining in MOEs from P2-IRES-τLacZ mice crossed to LBR-expressing transgenics shows reduced X-gal signal, supporting a repressive effect on OR expression (). Neurons that retained high β-gal protein expression often failed to express the LBR transgene, as demonstrated by IF for β-gal and GFP in sections of these mice (). Because OMP drives LBR expression only after OR choice, this result indicates postchoice downregulation of this P2 allele and the rest of the OR repertoire.
LBR Expression Inhibits OR Transcription
To test whether the inhibitory effects of LBR expression apply to genes that do not follow the spatial regulation of endogenous OR genes, we used a transgenic OR that is under the control of the tetO promoter (tetO MOR28-IRES-τLacZ). This transgene also carries H3K9me3 and H4K20me3 (data not shown), but unlike the endogenous ORs, its heterochromatization is not OSN specific and is probably caused by its multicopy (16 tandem copies) insertion (Garrick et al., 1998
). This transgene does not interact with either the endogenous ORs or the H enhancer (Figure S1I
and data not shown). In agreement with the repressive signature of this transgene, its expression is sporadic when crossed to the OMP-IRES-tTA driver but increases in frequency when crossed to LBR-expressing trans-genics (Figure S5B
). This is consistent with a simple mode of gene regulation in which chromatin decompaction allows tTA binding on the tetO promoter in more cells and transcriptional activation at higher frequency; such a linear and straightforward model does not apply to the endogenous ORs.
To determine the genome-wide effects of LBR expression in OSNs, we performed RNA-seq from whole MOE preparations. To increase the proportion of LBR+ OSNs in this mixed population, we combined two tTA drivers (OMP-IRES-tTA and Gγ8 tTA). In agreement with the observations in the FAC-sorted neurons, LBR expression in OSNs induces an ~8-fold downregulation of total OR expression (). Though OR expression is downregulated, most of the genes detected in OSNs (~14,000 genes) are not affected by ectopic LBR expression (367 non-OR transcripts significantly downregulated and 873 transcripts significantly upregulated, Cuffdiff FDR 0.05, genes with at least 10 reads). Interestingly, expression of some markers of immature OSN and progenitor cell populations increases in these animals (Figure S5A
). Overall, these changes suggest a partial transition toward a less-differentiated state, rather than elimination of mature OSNs, as shown by the lack of increased apoptosis in the transgenic MOE (Figure S5C
Ectopic LBR Expression Disrupts the Singularity of OR Transcription
Downregulation of OR transcription in LBR+ neurons is counterintuitive considering that the accessibility of OR chromatin increases and that these loci are stripped from HP1β. The disruption of long-range interactions with activating enhancers, like H, could contribute to this downregulation. Moreover, although global decompaction of OR chromatin might make all of the OR alleles transcriptionally competent, it is possible that OSNs cannot support transcription of ~2,800 OR alleles at the levels of a singularly transcribed OR. To test this, we performed single-cell RT-PCR with degenerate OR primers, followed by restriction enzyme digestion and electrophoresis (Buck and Axel, 1991
; and S6A
). We obtained 10 single-cell cDNA libraries from each genotype (see Extended Experimental Procedures) and examined OR representation by DraI, HinfI, and MseI digestion.
LBR Expression Induces OR Co-expression and Ectopic Targeting to the Olfactory Bulb
The complexity of the degenerate OR amplicons is different between control and LBR+ OSNs. In every LBR+ OSN, the base pair sum of the individual digestion products exceeds the length of the undigested PCR, whereas control amplicons contain only one product (, showing five amplicons from each genotype. Similar results were obtained for the other five amplicons and with MboI digestion [data not shown]). Sequencing 10 clones each from two control and two LBR-expressing amplicons verifies that the OR transcriptome is more complex in LBR+ OSNs. In both control amplicons, all 10 clones were identical, and the sequence of the cloned OR matches the digestion pattern. In the case of the LBR+ neurons, 10/10 and 9/10 clones of each amplicon were unique and different from each other. Moreover, the sequences of these clones did not match their digestion pattern, suggesting that these libraries are very complex and the observed distinct bands represent comigrating bands of similar size digested from a large number of different ORs. To verify this, we sequenced 96 independent colonies from a third LBR+ single-cell library and identified 46 different ORs from 11 chromosomes without any zonal restrictions (Table S3
). Thus, LBR expression in OSNs violates the “one receptor per neuron” rule and induces coexpression of a large number of ORs. To exclude that multiple ORs are expressed in low levels also in control OSNs but are masked by the highly expressed chosen allele, we pooled equal volumes of the single-cell RT reaction from a control and an LBR+ neuron and performed degenerate PCR and digest. Although the OR amplicon from the control neuron dominates the reaction, the ectopically coexpressed ORs from the LBR-expressing neuron are still detectable upon digestion (Figure S6B
For a non-PCR based confirmation for OR coexpression in LBR+ OSNs, we performed IF for ORs M50, M71, and C6 in tetO-LBR/tetO-MOR28 double transgenics. Most OSNs that retain detectable OR levels by IF are LBR negative; thus, in these OSNs, OR coexpression is extremely rare. To overcome this, we exploited the frequent expression of transgenic MOR28 in the LBR-expressing mice. IF in MOE sections from double transgenics revealed double-positive OSNs (), which are not detected in the absence of LBR, as previously shown (Fleischmann et al., 2008
; Nguyen et al., 2007
). Most likely, these OSNs are shutting down the endogenous OR in response to LBR expression while activating the decondensing transgenic OR.
The identity of the expressed OR allele instructs the targeting of like neurons to a single glomerulus (Mombaerts, 2006
). For this reason, we examined the targeting of neurons expressing the P2-IRES-τLacZ allele in LBR transgenics by IF for β-gal in olfactory bulb sections. β-gal protein is stable in the axons long after the cytoplasmic signal has faded (Figure S6D
); thus, we can use this approach to examine the targeting consequences of LBR expression. In control mice, the β-gal-positive fibers coalesce in distinct glomeruli (), with very few axons targeting wrong glomeruli (extreme example shown in Figure S6C
). However, upon LBR expression, β-gal-positive fibers extend to an extraordinary number of glomeruli (~30 per hemisphere per mouse). We detected extra distinct glomeruli and stray fibers both near wild-type P2 glomeruli and in ectopic positions in the bulb ( and S6C