The three-dimensional (3D) architecture of the mammalian nucleus is now being unraveled thanks to the recent development of chromatin conformation capture (3C) technologies. Here we report the results of a combined multiscale analysis of genome-wide mean replication timing and chromatin conformation data that reveal some intimate relationships between chromatin folding and human DNA replication. We previously described megabase replication N/U-domains as mammalian multiorigin replication units, and showed that their borders are 'master' replication initiation zones that likely initiate cascades of origin firing responsible for the stereotypic replication of these domains. Here, we demonstrate that replication N/U-domains correspond to the structural domains of self-interacting chromatin, and that their borders act as insulating regions both in high-throughput 3C (Hi-C) data and high-resolution 3C (4C) experiments. Further analyses of Hi-C data using a graph-theoretical approach reveal that N/U-domain borders are long-distance, interconnected hubs of the chromatin interaction network. Overall, these results and the observation that a well-defined ordering of chromatin states exists from N/U-domain borders to centers suggest that 'master' replication initiation zones are at the heart of a high-order, epigenetically controlled 3D organization of the human genome.