Like this?
using GMT
city_coords = [12.50 41.88; -9.13 38.72; -3.68 40.44; 2.32 48.85; 23.74 37.97; 16.36 48.21; 10.19 36.79];
city_names = ["Rome", "Lisbon", "Madrid", "Lutecia", "Athens", "Vienna", "Tunis"];
cities = mat2ds(city_coords, city_names);
# The image will be downloaded automatically
grdimage("@earth_day_03m", region=(-10,28,34,51), proj=:guess)
plot!(cities, connection="p", lw=1, lc=:yellow, marker=:circ, ms=0.2, mc=:white)
text!(cities, offset=0.3, font=(10,:red), show=true)
