Histone modifications play an important role in chromatin organization and transcriptional regulation, but despite the large amount of genome-wide histone modification data collected in different cells and tissues, little is known about co-occurrence of modifications on the same nucleosome. Here we present a genome-wide quantitative method for combinatorial indexed chromatin immunoprecipitation (co-ChIP) to characterize co-occurrence of histone modifications on nucleosomes. Using co-ChIP, we study the genome-wide co-occurrence of 14 chromatin marks (70 pairwise combinations), and find previously undescribed co-occurrence patterns, including the co-occurrence of H3K9me1 and H3K27ac in super-enhancers. Finally, we apply co-ChIP to measure the distribution of the bivalent H3K4me3-H3K27me3 domains in two distinct mouse embryonic stem cell (mESC) states and in four adult tissues. We observe dynamic changes in 5,786 regions and discover both loss and de novo gain of bivalency in key tissue-specific regulatory genes, suggesting a functional role for bivalent domains during different stages of development. These results show that co-ChIP can reveal the complex interactions between histone modifications.