Color and morphological differentiation in the Sinaloa Wren (Thryophilus sinaloa) in the tropical dry forests of Mexico: The role of environment and geographic isolation.

The role and the degree to which environment and geographic isolation contribute to phenotypic diversity has been widely debated. Here, we studied phenotypic variation (morphology and plumage reflectance) in the Sinaloa Wren, an endemic bird distributed throughout the tropical dry forest (TDF) on the Mexican pacific slope where a pronounced variability in environmental conditions has been reported. In particular, we aimed: 1) to characterize phenotypic variation between subspecies; 2) to analyze the relationship between phenotypic and environmental variation in the context of classic ecogeographic rules, such as Bergmann's, Allen's, Gloger's, and Bogert's, and to quantify the relative roles of environment and geographic isolation and their interaction in shaping phenotypic variation; and 3) to test for niche conservatism between subspecies. Our data revealed significant differences among subspecies morphology and plumage reflectance. The environment explained a higher proportion of the morphological variation, while geography explained a smaller proportion. However, variation in plumage reflectance was mainly explained by the joint effect of geography and environment. Our data did not support for Bergmann´s and Allen´s rule. However, longer tails and wings were positively associated with higher elevations, larger tarsus and culmens were positively related to higher latitudes and to greater tree cover, respectively. Our data partially supported Gloger´s rule, where darker plumages were associated with more humid environments. The effects of temperature on plumage coloration were more consistent with Bogert´s rule. In addition, we found darker plumages related to higher levels of UV-B radiation. Finally, niche divergence was detected between T. s. cinereus and T. s. sinaloa vs. T. s. russeus. In a continuously distributed ecosystem such as the TDF on the pacific slope of Mexico, the environmental conditions and geographic isolation have played an important role in promoting phenotypic differentiation in the Sinaloa Wren.

The geography of evolutionary divergence in the highly endemic avifauna from the Sierra Madre del Sur, Mexico.

BACKGROUND: Mesoamerica is a remarkable region with a high geological and ecological complexity. Within northern Mesoamerica, the biotic province of the Sierra Madre del Sur (SMS) in southwestern Mexico harbors exceptionally high avian endemism and diversity. Herein, we searched for spatially and temporally concordant phylogeographic patterns, in four bird genera from three distinct avian orders co-distributed across Mesoamerica and investigated their causes through hypothesis testing regarding historical processes. Selected species include endemic and differentiated populations across the montane forests of Mesoamerica, and particularly within the SMS. RESULTS: We gathered mitochondrial DNA sequences for at least one locus from 177 individuals across all species. We assessed genetic structure, demographic history, and defined a framework for the coalescent simulations used in biogeographic hypothesis testing temporal and spatial co-variance. Our analyses suggested shared phylogeographic breaks in areas corresponding to the SMS populations, and between the main montane systems in Mesoamerica, with the Central Valley of Oaxaca and the Nicaragua Depression being the most frequently shared breaks among analyzed taxa. Nevertheless, dating analyses and divergence patterns observed were consistent with the hypothesis of broad vicariance across Mesoamerica derived from mechanisms operating at distinct times across taxa in the SMS. CONCLUSIONS: Our study provides a framework for understanding the evolutionary origins and historical factors enhancing speciation in well-defined regions within Mesoamerica, indicating that the evolutionary history of extant biota inhabiting montane forests is complex and often idiosyncratic.

The influence of climatic niche preferences on the population genetic structure of a mistletoe species complex.

The prevalent view on genetic structuring in parasitic plants is that host-race formation is caused by varying degrees of host specificity. However, the relative importance of ecological niche divergence and host specificity to population differentiation remains poorly understood. We evaluated the factors associated with population differentiation in mistletoes of the Psittacanthus schiedeanus complex (Loranthaceae) in Mexico. We used genetic data from chloroplast sequences and nuclear microsatellites to study population genetic structure and tested its association with host preferences and climatic niche variables. Pairwise genetic differentiation was associated with environmental and host preferences, independent of geography. However, environmental predictors appeared to be more important than host preferences to explain genetic structure, supporting the hypothesis that the occurrence of the parasite is largely determined by its own climatic niche and, to a lesser degree, by host specificity. Genetic structure is significant within this mistletoe species complex, but the processes associated with this structure appear to be more complex than previously thought. Although host specificity was not supported as the major determinant of population differentiation, we consider this to be part of a more comprehensive ecological model of mistletoe host-race formation that incorporates the effects of climatic niche evolution.

A mistletoe tale: postglacial invasion of Psittacanthus schiedeanus (Loranthaceae) to Mesoamerican cloud forests revealed by molecular data and species distribution modeling.

BACKGROUND: Ecological adaptation to host taxa is thought to result in mistletoe speciation via race formation. However, historical and ecological factors could also contribute to explain genetic structuring particularly when mistletoe host races are distributed allopatrically. Using sequence data from nuclear (ITS) and chloroplast (trnL-F) DNA, we investigate the genetic differentiation of 31 Psittacanthus schiedeanus (Loranthaceae) populations across the Mesoamerican species range. We conducted phylogenetic, population and spatial genetic analyses on 274 individuals of P. schiedeanus to gain insight of the evolutionary history of these populations. Species distribution modeling, isolation with migration and Bayesian inference methods were used to infer the evolutionary transition of mistletoe invasion, in which evolutionary scenarios were compared through posterior probabilities. RESULTS: Our analyses revealed shallow levels of population structure with three genetic groups present across the sample area. Nine haplotypes were identified after sequencing the trnL-F intergenic spacer. These haplotypes showed phylogeographic structure, with three groups with restricted gene flow corresponding to the distribution of individuals/populations separated by habitat (cloud forest localities from San Luis Potosí to northwestern Oaxaca and Chiapas, localities with xeric vegetation in central Oaxaca, and localities with tropical deciduous forests in Chiapas), with post-glacial population expansions and potentially corresponding to post-glacial invasion types. Similarly, 44 ITS ribotypes suggest phylogeographic structure, despite the fact that most frequent ribotypes are widespread indicating effective nuclear gene flow via pollen. Gene flow estimates, a significant genetic signal of demographic expansion, and range shifts under past climatic conditions predicted by species distribution modeling suggest post-glacial invasion of P. schiedeanus mistletoes to cloud forests. However, Approximate Bayesian Computation (ABC) analyses strongly supported a scenario of simultaneous divergence among the three groups isolated recently. CONCLUSIONS: Our results provide support for the predominant role of isolation and environmental factors in driving genetic differentiation of Mesoamerican parrot-flower mistletoes. The ABC results are consistent with a scenario of post-glacial mistletoe invasion, independent of host identity, and that habitat types recently isolated P. schiedeanus populations, accumulating slight phenotypic differences among genetic groups due to recent migration across habitats. Under this scenario, climatic fluctuations throughout the Pleistocene would have altered the distribution of suitable habitat for mistletoes throughout Mesoamerica leading to variation in population continuity and isolation. Our findings add to an understanding of the role of recent isolation and colonization in shaping cloud forest communities in the region.

Molecular and iridescent feather reflectance data reveal recent genetic diversification and phenotypic differentiation in a cloud forest hummingbird.

The present day distribution and spatial genetic diversity of Mesoamerican biota reflects a long history of responses to habitat change. The hummingbird Lampornis amethystinus is distributed in northern Mesoamerica, with geographically disjunct populations. Based on sampling across the species range using mitochondrial DNA (mtDNA) sequences and nuclear microsatellites jointly analysed with phenotypic and climatic data, we (1) test whether the fragmented distribution is correlated with main evolutionary lineages, (2) assess body size and plumage color differentiation of populations in geographic isolation, and (3) evaluate a set of divergence scenarios and demographic patterns of the hummingbird populations. Analysis of genetic variation revealed four main groups: blue-throated populations (Sierra Madre del Sur); two groups of amethyst-throated populations (Trans-Mexican Volcanic Belt and Sierra Madre Oriental); and populations east of the Isthmus of Tehuantepec (IT) with males showing an amethyst throat. The most basal split is estimated to have originated in the Pleistocene, 2.39-0.57 million years ago (MYA), and corresponded to groups of populations separated by the IT. However, the estimated recent divergence time between blue- and amethyst-throated populations does not correspond to the 2-MY needed to be in isolation for substantial plumage divergence, likely because structurally iridescent colors are more malleable than others. Results of species distribution modeling and Approximate Bayesian Computation analysis fit a model of lineage divergence west of the Isthmus after the Last Glacial Maximum (LGM), and that the species' suitable habitat was disjunct during past and current conditions. These results challenge the generality of the contraction/expansion glacial model to cloud forest-interior species and urges management of cloud forest, a highly vulnerable ecosystem to climate change and currently facing destruction, to prevent further loss of genetic diversity or extinction.

Development and characterization of microsatellite loci in the mistletoe Psittacanthus schiedeanus (Loranthaceae).

PREMISE OF THE STUDY: Microsatellite primers were developed for the parasitic Psittacanthus schiedeanus, a common mistletoe species on cloud forest-adapted tree hosts in Mesoamerica, to investigate intraspecific genetic patterns of diversity and genetic structure. • METHODS AND RESULTS: Using an enriched library, 10 polymorphic microsatellite loci were developed in P. schiedeanus. All loci consisted of dinucleotide repeats. Average alleles per locus were 12 (4-17), and a total of 120 alleles were recorded across 39 individuals from four populations in Mexico. Primers were tested in 11 additional species, but only amplified successfully in P. calyculatus and P. angustifolius. • CONCLUSIONS: The polymorphic loci described will be useful in studies of genetic diversity and genetic population differentiation in natural populations of these parasitic plants, and will provide valuable information to understand the importance of host distribution.

Acoustic divergence with gene flow in a lekking hummingbird with complex songs.

Hummingbirds have developed a remarkable diversity of learned vocalizations, from single-note songs to phonologically and syntactically complex songs. In this study we evaluated if geographic song variation of wedge-tailed sabrewings (Campylopterus curvipennis) is correlated with genetic divergence, and examined processes that explain best the origin of intraspecific song variation. We contrasted estimates of genetic differentiation, genetic structure, and gene flow across leks from microsatellite loci of wedge-tailed sabrewings with measures for acoustic signals involved in mating derived from recordings of males singing at leks throughout eastern Mexico. We found a strong acoustic structure across leks and geography, where lek members had an exclusive assemblage of syllable types, differed in spectral and temporal measurements of song, and song sharing decreased with geographic distance. However, neutral genetic and song divergence were not correlated, and measures of genetic differentiation and migration estimates indicated gene flow across leks. The persistence of acoustic structuring in wedge-tailed sabrewings may thus best be explained by stochastic processes across leks, in which intraspecific vocal variation is maintained in the absence of genetic differentiation by postdispersal learning and social conditions, and by geographical isolation due to the accumulation of small differences, producing most dramatic changes between populations further apart.

Interglacial genetic diversification of Moussonia deppeana (Gesneriaceae), a hummingbird-pollinated, cloud forest shrub in northern Mesoamerica.

Recent empirical work on cloud forest-adapted species supports the role of both old divergences across major geographical areas and more recent divergences attributed to Pleistocene climate changes. The shrub Moussonia deppeana is distributed in northern Mesoamerica, with geographically disjunct populations. Based on sampling throughout the species range and employing plastid and nuclear markers, we (i) test whether the fragmented distribution is correlated with main evolutionary lineages, (ii) reconstruct its phylogeographical history to infer the history of cloud forest in northern Mesoamerica and (iii) evaluate a set of refugia/vicariance scenarios for the region and demographic patterns of the populations whose ranges expanded and tracked cloud forest conditions during the Last Glacial Maximum. We found a deep evolutionary split in M. deppeana about 6-3 Ma, which could be consistent with a Pliocene divergence. Comparison of variation in plastid and nuclear markers revealed several lineages mostly congruent with their isolated geographical distribution and restricted gene flow among groups. Results of species distribution modelling and coalescent simulations fit a model of multiple refugia diverging during interglacial cycles. The demographic history of M. deppeana is not consistent with an expanding-contracting cloud forest archipelago model during the Last Glacial Maximum. Instead, our data suggest that populations persisted across the geographical range throughout the glacial cycles, and experienced isolation and divergence during interglacial periods.

Permanent genetic resources added to Molecular Ecology Resources Database 1 February 2013-31 March 2013.

This article documents the addition of 142 microsatellite marker loci to the Molecular Ecology Resources database. Loci were developed for the following species: Agriophyllum squarrosum, Amazilia cyanocephala, Batillaria attramentaria, Fungal strain CTeY1 (Ascomycota), Gadopsis marmoratus, Juniperus phoenicea subsp. turbinata, Liriomyza sativae, Lupinus polyphyllus, Metschnikowia reukaufii, Puccinia striiformis and Xylocopa grisescens. These loci were cross-tested on the following species: Amazilia beryllina, Amazilia candida, Amazilia rutila, Amazilia tzacatl, Amazilia violiceps, Amazilia yucatanensis, Campylopterus curvipennis, Cynanthus sordidus, Hylocharis leucotis, Juniperus brevifolia, Juniperus cedrus, Juniperus osteosperma, Juniperus oxycedrus, Juniperus thurifera, Liriomyza bryoniae, Liriomyza chinensis, Liriomyza huidobrensis and Liriomyza trifolii.

Comparative phylogeographic analyses illustrate the complex evolutionary history of threatened cloud forests of northern Mesoamerica.

Comparative phylogeography can elucidate the influence of historical events on current patterns of biodiversity and can identify patterns of co-vicariance among unrelated taxa that span the same geographic areas. Here we analyze temporal and spatial divergence patterns of cloud forest plant and animal species and relate them to the evolutionary history of naturally fragmented cloud forests--among the most threatened vegetation types in northern Mesoamerica. We used comparative phylogeographic analyses to identify patterns of co-vicariance in taxa that share geographic ranges across cloud forest habitats and to elucidate the influence of historical events on current patterns of biodiversity. We document temporal and spatial genetic divergence of 15 species (including seed plants, birds and rodents), and relate them to the evolutionary history of the naturally fragmented cloud forests. We used fossil-calibrated genealogies, coalescent-based divergence time inference, and estimates of gene flow to assess the permeability of putative barriers to gene flow. We also used the hierarchical Approximate Bayesian Computation (HABC) method implemented in the program msBayes to test simultaneous versus non-simultaneous divergence of the cloud forest lineages. Our results show shared phylogeographic breaks that correspond to the Isthmus of Tehuantepec, Los Tuxtlas, and the Chiapas Central Depression, with the Isthmus representing the most frequently shared break among taxa. However, dating analyses suggest that the phylogeographic breaks corresponding to the Isthmus occurred at different times in different taxa. Current divergence patterns are therefore consistent with the hypothesis of broad vicariance across the Isthmus of Tehuantepec derived from different mechanisms operating at different times. This study, coupled with existing data on divergence cloud forest species, indicates that the evolutionary history of contemporary cloud forest lineages is complex and often lineage-specific, and thus difficult to capture in a simple conservation strategy.

Isolation and characterization of microsatellite markers in distylous Palicourea padifolia (Rubiaceae).

PREMISE OF THE STUDY: Eighteen microsatellite loci of distylous Palicourea padifolia were isolated and characterized for population genetics studies. METHODS AND RESULTS: Following a microsatellite enrichment protocol, 18 primer pairs amplified successfully, and the polymorphism of the loci was initially evaluated in 15 individuals from three populations in Mexico. Seven loci were polymorphic, and their variability was further assessed in 60 individuals from three populations. The total number of alleles per locus, combining samples from all populations, ranged from 7 to 12. Nei's genetic diversity average across loci was 0.681, 0.714, and 0.703 for El Riscal, Montebello, and Ruiz Cortínez, respectively. Palicourea padifolia individuals had one to four alleles, confirming the polyploidy of this species. CONCLUSIONS: These markers will facilitate genetic diversity studies of P. padifolia across its distribution range and facilitate investigations on the evolution of its breeding system.

Selection and geographic isolation influence hummingbird speciation: genetic, acoustic and morphological divergence in the wedge-tailed sabrewing (Campylopterus curvipennis).

BACKGROUND: Mesoamerica is one of the most threatened biodiversity hotspots in the world, yet we are far from understanding the geologic history and the processes driving population divergence and speciation for most endemic taxa. In species with highly differentiated populations selective and/or neutral factors can induce rapid changes to traits involved in mate choice, promoting reproductive isolation between allopatric populations that can eventually lead to speciation. We present the results of genetic differentiation, and explore drift and selection effects in promoting acoustic and morphological divergence among populations of Campylopterus curvipennis, a lekking hummingbird with an extraordinary vocal variability across Mesoamerica. RESULTS: Analyses of two mitochondrial genes and ten microsatellite loci genotyped for 160 individuals revealed the presence of three lineages with no contemporary gene flow: C. c. curvipennis, C. c. excellens, and C. c. pampa disjunctly distributed in the Sierra Madre Oriental, the Tuxtlas region and the Yucatan Peninsula, respectively. Sequence mtDNA and microsatellite data were congruent with two diversification events: an old vicariance event at the Isthmus of Tehuantepec (c. 1.4 Ma), and a more recent Pleistocene split, isolating populations in the Tuxtlas region. Hummingbirds of the excellens group were larger, and those of the pampa group had shorter bills, and lineages that have been isolated the longest shared fewer syllables and differed in spectral and temporal traits of a shared syllable. Coalescent simulations showed that fixation of song types has occurred faster than expected under neutrality but the null hypothesis that morphological divergence resulted from drift was not rejected. CONCLUSIONS: Our phylogeographic analyses uncovered the presence of three Mesoamerican wedge-tailed sabrewing lineages, which diverged at different time scales. These results highlight the importance of the Isthmus of Tehuantepec and more recent Pleistocene climatic events in driving isolation and population divergence. Coalescent analyses of the evolution of phenotypic traits suggest that selection is driving song evolution in wedge-tailed sabrewings but drift could not be rejected as a possibility for morphological divergence.

Permanent Genetic Resources added to Molecular Ecology Resources Database 1 August 2009-30 September 2009.

This article documents the addition of 238 microsatellite marker loci and 72 pairs of Single Nucleotide Polymorphism (SNP) sequencing primers to the Molecular Ecology Resources Database. Loci were developed for the following species: Adelges tsugae, Artemisia tridentata, Astroides calycularis, Azorella selago, Botryllus schlosseri, Botrylloides violaceus, Cardiocrinum cordatum var. glehnii, Campylopterus curvipennis, Colocasia esculenta, Cynomys ludovicianus, Cynomys leucurus, Cynomys gunnisoni, Epinephelus coioides, Eunicella singularis, Gammarus pulex, Homoeosoma nebulella, Hyla squirella, Lateolabrax japonicus, Mastomys erythroleucus, Pararge aegeria, Pardosa sierra, Phoenicopterus ruber ruber and Silene latifolia. These loci were cross-tested on the following species: Adelges abietis, Adelges cooleyi, Adelges piceae, Pineus pini, Pineus strobi, Tubastrea micrantha, three other Tubastrea species, Botrylloides fuscus, Botrylloides simodensis, Campylopterus hemileucurus, Campylopterus rufus, Campylopterus largipennis, Campylopterus villaviscensio, Phaethornis longuemareus, Florisuga mellivora, Lampornis amethystinus, Amazilia cyanocephala, Archilochus colubris, Epinephelus lanceolatus, Epinephelus fuscoguttatus, Symbiodinium temperate-A clade, Gammarus fossarum, Gammarus roeselii, Dikerogammarus villosus and Limnomysis benedeni. This article also documents the addition of 72 sequencing primer pairs and 52 allele specific primers for Neophocaena phocaenoides.

Between-year changes in functional gender expression of Palicourea padifolia (Rubiaceae), a distylous, hummingbird-pollinated shrub.

BACKGROUND AND AIMS: Because distylous species have two hermaphroditic style-length floral morphs, they face two sex allocation problems: the equilibrium morph ratio and the optimal allocation to pollen and seed production in each floral morph. Gender specialization is expected among distylous species when floral morphs differ in reproductive output. However, spatio-temporal variability in female reproductive output between morphs needs to be investigated to fully understand patterns of sexual expression and gender specialization in distylous plants. Between-year variation in flower and fruit production of hummingbird-pollinated Palicourea padifolia (Rubiaceae) was examined, focusing on functional gender expression of long- and short-styled morphs and comparing their reproductive performance in five consecutive years (1998-2002). METHODS: Between-year variation in inflorescence, floral bud and fruit production was monitored and quantified. These traits were then used as parameters to determine functional gender differences between floral morphs through time. KEY RESULTS: Inflorescence production varied among years but no significant differences were found between floral morphs. Long-styled plants initiated more floral buds per inflorescence every year than short-styled plants, suggesting higher allocation to pollinator attraction and, potentially, an increase in male fitness through pollen donation. Although fruit production was similar between morphs, their functional gender shifted across years. CONCLUSIONS: The gender expression inconsistency across years is surprising because a number of floral characters and attributes that contribute to differently attracting and rewarding effective pollinators in P. padifolia suggest gender specialization. The evidence that morphs of distylous species might specialize in functional gender mostly comes from differences among populations in seed production and non-equilibrium morph ratios based on 1-year field population surveys. The results suggest that more sampling through time is needed to detect gender specialization among distylous species with a perennial habit.

Reproductive ecology of distylous Palicourea Padifolia (Rubiaceae) in a tropical montane cloud forest. I. Hummingbirds' effectiveness as pollen vectors.

The adaptiveness of distyly has been typically investigated in terms of its female function, specifically pollen receipt. However, pollen loads on stigmas can only provide moderate support for Darwin's hypothesis of the promotion of legitimate crosses. To determine the effectiveness of hummingbirds as pollen vectors between floral morphs and the consequences in terms of male (pollen transfer) and female function (pollen receipt) in Palicourea padifolia (Rubiaceae), floral visitors, their foraging modes, and temporal patterns of floral visitation were observed and documented. Differences in pollen and stigma morphology, pollen flow, rates of pollen deposition, and/or stigmatic pollen loads were then evaluated for their contribution toward differences in reproductive output between floral morphs. A pollination experiment with stuffed hummingbirds that varied in bill size was done to evaluate the contribution of bill variation toward differences between floral morphs in pollen receipt and pollen transfer and female reproductive output. Anthers of long-styled flowers contained significantly more and smaller pollen grains than those of short-styled flowers, independently of corolla and anther lengths. The shape and orientation of the stigma lobes differed between morphs and were significantly longer among short-styled flowers. Hummingbird visitation rates did not differ significantly between floral morphs, and foraging movements from focal plants towards neighboring plants were independent of floral morph. Stigmatic pollen loads under field conditions and those after controlled hummingbird visitation, along with rates of pollen accumulation through the day indicated that stigmas of short-styled flowers receive proportionately more legitimate (intermorph) pollen grains than did those of long-styled flowers. However, the species of hummingbird was marginally significant in explaining variation in pollen deposition on stigmas. Lastly, intermorph pollinations of P. padifolia resulted in significant differences in fruit production between floral morphs, independent of pollination treatment and pollinator species; short-styled flowers proportionately developed almost twice the number of fruits developed by long-styled flowers.

Reproductive ecology of distylous Palicourea padifolia (Rubiaceae) in a tropical montane cloud forest. II. Attracting and rewarding mutualistic and antagonistic visitors.

By definition, the floral morphs of distylous plants differ in floral architecture. Yet, because cross-pollination is necessary for reproductive success in both morphs, they should not differ in attributes that contribute to attracting and rewarding floral visitors. Floral and vegetative attributes that function in distylous polymorphism in hummingbird-pollinated Palicourea padifolia (Rubiaceae) and the responses of pollinators and insect herbivores to the resources offered by both morphs were investigated. The performance of each morph along multiple stages of the reproductive cycle, from inflorescence and nectar production to fruit production, was surveyed, and pollinator behavior and nectar standing crops were then observed. Costs associated with such attractiveness were also evaluated in terms of herbivore attack and of plant reproductive fitness (female function) as a function of leaf herbivory. The number of inflorescences, floral buds, open flowers, and ripe fruits offered by either floral morph were similar, but short-styled plants almost doubled the number of developing fruits of long-styled plants. Long-styled flowers produced higher nectar volumes and accumulated more nectar over time than short-styled flowers. Measures of nectar standing crop and data on pollinator behavior suggest that hummingbirds respond to this morph-specific scheduling of nectar production. Lastly, long-styled plants suffered a higher herbivore attack and lost more leaf area over time than those with short-styled flowers. Herbivory was negatively correlated with fruit number and fruit mass, and long-styled plants set significantly less fruit mass than short-styled plants. The results suggest that pollinators and herbivores may exert selective pressures on floral and vegetative traits that could also influence gender function.

High-performance thin-layer chromatography-bioautography for multiple antibiotic residues in cow's milk.

An analytical method to identify and quantify multiple antibiotic residues (chloramphenicol, ampicillin, benzylpenicillin, dicloxacillin and erythromycin) in cow's milk by high-performance thin-layer chromatography (HPTLC) combined with bioautography was developed. The test microorganism used for bioautography was Bacillus subtilis ATCC 6633. Antibiotic residues were extracted with acetonitrile, fat eliminated with petroleum ether and residues isolated with dichloromethane The sensitivity of the method guarantees the detection of the above-mentioned antibiotics at levels below maximum residue limits (MRL) allowed for milk. Percentage recoveries ranged between 90 and 100%, with coefficients of variation between 7.2 and 21.3%. Some advantages of this methodology over thin-layer chromatography (TLC)/bioautography are also discussed.