The authors propose a fully nonparametric algorithm that accurately recovers the causal polytree structure from data, even when the number of variables is much larger than the sample size.