The partition scheme and evolutionary models selected were then used for maximum likelihood phylogenetic reconstruction using Garli 2.0.1 (Zwickl, 2006) and for Bayesian phylogenetic analysis using MrBayes 3.2.2 (Ronquist and Huelsenbeck, 2003). For our Garli analyses, we ran 50 bootstrap replications using the default parameter values recommended for datasets with small numbers of taxa, with 2 searches per run and 4 individuals in the population for the genetic algorithm. For our Bayesian analyses, we ran 4 chains and 2 runs of 11 million generations each, with trees sampled every 1000 generations, discarding the first 1 million generations and 1000 trees as burnin. All Garli and MrBayes runs were executed on the CIPRES Science Gateway v 3.3 (Miller et al., 2010). As we were able to collect less sequence data overall on Ateles marginatus than on any of the other taxa (and as we were unable to sequence any portion of the NAD6 gene for that sample), we ran our phylogenetic analyses in three ways: (a) including data from Ateles marginatus and analyzing all three coding regions, (b) including A. marginatus, but only considering the NAD5 and cyto- chrome B loci, and (c) excluding A. marginatus completely from the alignment. Once the phylogenetic trees were inferred, we replaced the sample names with the species name assigned according to Groves (2001) to check if the samples putatively assigned to a specific taxon formed monophyletic groups.