Constraint-Based Learning Of Interventional Markov Equivalence Classes On High-Dimensional Data