Correlations in gene expression are used to infer functional and regulatory relationships between genes. However, correlations are often calculated across different cell types or perturbations, causing genes with unrelated functions to be correlated. Here, we demonstrate that correlated modules can be better captured by measuring correlations of steady-state gene expression fluctuations in single cells. We report a high-precision single-cell RNA-seq method called MALBAC-DT to measure the correlation between any pair of genes in a homogenous cell population. Using this method, we were able to identify numerous cell-type specific and functionally enriched correlated gene modules. We confirmed through knockdown that a module enriched for p53 signaling predicted p53 regulatory targets more accurately than a consensus of ChIP-seq studies and that steady-state correlations were predictive of transcriptome-wide response patterns to perturbations. This approach provides a powerful way to advance our functional understanding of the genome.

Olfactory sensory neurons (OSNs) are functionally defined by their expression of a unique odorant receptor (OR). Mechanisms underlying singular OR expression are well studied, and involve a massive cross-chromosomal enhancer interaction network. Trace amine-associated receptors (TAARs) form a distinct family of olfactory receptors, and here we find that mechanisms regulating Taar gene choice display many unique features. The epigenetic signature of Taar genes in TAAR OSNs is different from that in OR OSNs. We further identify that two TAAR enhancers conserved across placental mammals are absolutely required for expression of the entire Taar gene repertoire. Deletion of either enhancer dramatically decreases the expression probabilities of different Taar genes, while deletion of both enhancers completely eliminates the TAAR OSN populations. In addition, both of the enhancers are sufficient to drive transgene expression in the partially overlapped TAAR OSNs. We also show that the TAAR enhancers operate in cis to regulate Taar gene expression. Our findings reveal a coordinated control of Taar gene choice in OSNs by two remote enhancers, and provide an excellent model to study molecular mechanisms underlying formation of an olfactory subsystem. In our nose, some neuron subpopulations express a family of trace amine associated receptors (TAARs, smelling e.g., rotten fish). Fei et al. identify two conserved enhancers across placental mammals named TAAR enhancer 1 and 2 that coordinately regulate expression of the entire Taar gene repertoire.

Three-dimensional genome structures play a key role in gene regulation and cell functions. Characterization of genome structures necessitates single-cell measurements. This has been achieved for haploid cells but has remained a challenge for diploid cells. We developed a single-cell chromatin conformation capture method, termed Dip-C, that combines a transposon-based whole-genome amplification method to detect many chromatin contacts, called META (multiplex end-tagging amplification), and an algorithm to impute the two chromosome haplotypes linked by each contact. We reconstructed the genome structures of single diploid human cells from a lymphoblastoid cell line and from primary blood cells with high spatial resolution, locating specific single-nucleotide and copy number variations in the nucleus. The two alleles of imprinted loci and the two X chromosomes were structurally different. Cells of different types displayed statistically distinct genome structures. Such structural cell typing is crucial for understanding cell functions.

Different regions of the mammalian nose smell different odors. In the mouse olfactory system, spatially regulated expression of >1000 olfactory receptors (ORs) along the dorsomedial-ventrolateral (DV) axis forms a topological map in the main olfactory epithelium (MOE). However, the locations of most ORs along the DV axis are currently unknown. By sequencing mRNA of 12 isolated MOE pieces, we mapped out the DV locations-as quantified by "zone indices" on a scale of 1-5-of 1033 OR genes with an estimated error of 0.3 zone indices. Our map covered 81% of all intact OR genes and 99.4% of the total OR mRNA abundance. Spatial regulation tended to vary gradually along chromosomes. We further identified putative non-OR genes that may exhibit spatial expression along the DV axis.