Phlebotomine sand flies (Diptera: Psychodidae: Phlebotominae) have biological relevance as vectors of several pathogens. To ensure periodic entomological monitoring it is necessary to have efficient and accurate tools for an adequate taxonomic identification. There are only few studies on phylogenetic analyses of phlebotomine sand flies from Neotropics, based mostly on morphological and/or molecular data, which makes the delimitation of intra- and interspecific variability of species challenging. Here we generated new molecular information on sand fly species distributed in endemic areas of leishmaniasis in Mexico, using mitochondrial and ribosomal genes, and incorporating morphological information available. Specifically, we established their phylogenetic relationships, and estimated their divergence time. Our study provides molecular information for 15 phlebotomine sand fly species from different areas of Mexico, contributing to the genetic inventory and phylogenetic relations among Neotropical species of the subfamily Phlebotominae. Mitochondrial genes proved to be suitable markers for the molecular identification of phlebotomine sand flies. However, the incorporation of additional nuclear gene information could increase the significance of phylogenetic inferences. We also provided evidence about a possible divergence time of phlebotomine sand fly species, supporting their presumable origin in the Cretaceous period.