The neotropical genus Varronia P. Browne is one of the most problematic, taxonomically, among the Boraginales. This may be due to the frequency of overlapping the boundaries between intra and interspecific variation of some characters. The purpose of this project was to explore the potential of foliar morphoanatomic characters for the delimitation of the species found within the Varronia bullata complex in Panama: Varronia inermis (Mill.) Borhidi, Varronia bullata L. and Varronia linnaei (Stearn) J.S. Windmill. For this purpose, mature leaves were collected in different populations for each species. Cross sections were made, and superficial sections of the leaf lamina were obtained, which were used to quantify and qualify previously selected morphoanatomic characters. The quantitative information was subjected to a discriminant canonical analysis. The results of the discriminating canonical analysis that detect the density of simple trichomes in the abaxial epidermis and the length of simple trichomes in the adaxial epidermis, are the characters that have the greatest influence on species separation. On the other hand, a descriptive way is reflected that the presence of crystalline sales, the type of trichomes and the distribution of stomata also contribute to the taxonomic separation of species.El género neotropical Varronia P. Browne es uno de los más problemáticos, taxonómicamente, entre las Boraginales. Esto puede deberse a la frecuente superposición de los límites entre la variación intra e interespecífica de algunos caracteres. La finalidad de este trabajo fue explorar el potencial de los caracteres morfoanatómicos foliares para l a delimitación de las especies que se encuentran dentro del complejo Varronia bullata en Panamá: Varronia inermis (Mill.) Borhidi, Varronia bullata L. y Varronia linnaei (Stearn) J.S. Mill. Para esto, se colectaron hojas maduras en diferentes poblaciones de cada especie. Se realizaron cortes transversales y se obtuvieron secciones superficiales de la lámina foliar, los cuales fueron utilizados para cuantificar y cualificar caracteres morfoanatómicos previamente seleccionados. La información cuantitativa fue sometida a un análisis canónico discriminante. Los resultados del análisis canónico discriminante mostraron que la densidad de los tricomas simples en la epidermis abaxial y la longitud de los tricomas simples en la epidermis adaxial, son los caracteres que presentan mayor influencia en la separación de las especies; por otro lado, a manera descriptiva se refleja que la presencia de sales cristalinas, el tipo de tricomas y la distribución de los estomas también contribuyen en la separación taxonómica de las especies
