<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article
  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article" dtd-version="3.0" xml:lang="en">
<front>
<journal-meta>
<journal-id journal-id-type="nlm-ta">PLoS ONE</journal-id>
<journal-id journal-id-type="publisher-id">plos</journal-id>
<journal-id journal-id-type="pmc">plosone</journal-id><journal-title-group>
<journal-title>PLoS ONE</journal-title></journal-title-group>
<issn pub-type="epub">1932-6203</issn>
<publisher>
<publisher-name>Public Library of Science</publisher-name>
<publisher-loc>San Francisco, USA</publisher-loc></publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">PONE-D-14-29808</article-id>
<article-id pub-id-type="doi">10.1371/journal.pone.0111057</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Biology and life sciences</subject><subj-group><subject>Zoology</subject><subj-group><subject>Animal behavior</subject><subj-group><subject>Animal sexual behavior</subject></subj-group></subj-group></subj-group></subj-group></article-categories>
<title-group>
<article-title>Males of a Strongly Polygynous Species Consume More Poisonous Food than Females</article-title>
<alt-title alt-title-type="running-head">Male-Biased Consumption of Poisonous Insects</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Bravo</surname><given-names>Carolina</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib>
<contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Bautista</surname><given-names>Luis Miguel</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib>
<contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>García-París</surname><given-names>Mario</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib>
<contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Blanco</surname><given-names>Guillermo</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib>
<contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Alonso</surname><given-names>Juan Carlos</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib>
</contrib-group>
<aff id="aff1"><label>1</label><addr-line>Department of Evolutionary Ecology, Museo Nacional de Ciencias Naturales, CSIC, Madrid, Madrid, Spain</addr-line></aff>
<aff id="aff2"><label>2</label><addr-line>Department of Biodiversity and Evolutionary Biology, Museo Nacional de Ciencias Naturales, CSIC, Madrid, Madrid, Spain</addr-line></aff>
<contrib-group>
<contrib contrib-type="editor" xlink:type="simple"><name name-style="western"><surname>McGraw</surname><given-names>Kevin</given-names></name>
<role>Editor</role>
<xref ref-type="aff" rid="edit1"/></contrib>
</contrib-group>
<aff id="edit1"><addr-line>Arizona State University, United States of America</addr-line></aff>
<author-notes>
<corresp id="cor1">* E-mail: <email xlink:type="simple">c.bravo@mncn.csic.es</email></corresp>
<fn fn-type="conflict"><p>The authors have declared that no competing interests exist.</p></fn>
<fn fn-type="con"><p>Conceived and designed the experiments: CB JCA LMB GB. Performed the experiments: CB. Analyzed the data: CB LMB. Contributed reagents/materials/analysis tools: CB LMB MGP GB JCA. Wrote the paper: CB JCA LMB MGP GB.</p></fn>
</author-notes>
<pub-date pub-type="collection"><year>2014</year></pub-date>
<pub-date pub-type="epub"><day>22</day><month>10</month><year>2014</year></pub-date>
<volume>9</volume>
<issue>10</issue>
<elocation-id>e111057</elocation-id>
<history>
<date date-type="received"><day>7</day><month>7</month><year>2014</year></date>
<date date-type="accepted"><day>28</day><month>9</month><year>2014</year></date>
</history>
<permissions>
<copyright-year>2014</copyright-year>
<copyright-holder>Bravo et al</copyright-holder><license xlink:type="simple"><license-p>This is an open-access article distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/" xlink:type="simple">Creative Commons Attribution License</ext-link>, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</license-p></license></permissions>
<abstract>
<p>We present evidence of a possible case of self-medication in a lekking bird, the great bustard <italic>Otis tarda</italic>. Great bustards consumed blister beetles (Meloidae), in spite of the fact that they contain cantharidin, a highly toxic compound that is lethal in moderate doses. In addition to anthelminthic properties, cantharidin was effective against gastrointestinal bacteria that cause sexually-transmitted diseases. Although both sexes consumed blister beetles during the mating season, only males selected them among all available insects, and ingested more and larger beetles than females. The male-biased consumption suggests that males could use cantharidin to reduce their parasite load and increase their sexual attractiveness. This plausibly explains the intense cloaca display males perform to approaching females, and the meticulous inspection females conduct of the male's cloaca, a behaviour only observed in this and another similar species of the bustard family. A white, clean cloaca with no infection symptoms (e.g., diarrhoea) is an honest signal of both, resistance to cantharidin and absence of parasites, and represents a reliable indicator of the male quality to the extremely choosy females. Our results do not definitely prove, but certainly strongly suggest that cantharidin, obtained by consumption of blister beetles, acts in great bustards as an oral anti-microbial and pathogen-limiting compound, and that males ingest these poisonous insects to increase their mating success, pointing out that self-medication might have been overlooked as a sexually-selected mechanism enhancing male fitness.</p>
</abstract>
<funding-group><funding-statement>C.B. was supported by a contract HENARSA-CSIC. This study was funded by projects CGL2005-04893, CGL2008-02567, CGL2012-36345 and CGL2009-12753-C02-01 of the Dirección General de Investigación, Spanish Ministry of Economy and Competitiveness. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</funding-statement></funding-group><counts><page-count count="12"/></counts><custom-meta-group><custom-meta id="data-availability" xlink:type="simple"><meta-name>Data Availability</meta-name><meta-value>The authors confirm that all data underlying the findings are fully available without restriction. All relevant data are within the paper and its Supporting Information files.</meta-value></custom-meta></custom-meta-group></article-meta>
</front>
<body><sec id="s1">
<title>Introduction</title>
<p>In polygynous birds, where competition among males for access to females is particularly strong, parasite load of males is assessed by females before mating and hence may affect male breeding success <xref ref-type="bibr" rid="pone.0111057-Ekblom1">[1]</xref>–<xref ref-type="bibr" rid="pone.0111057-Loye1">[6]</xref>. Self-medication could be one of the mechanisms males use to appear as healthy and vigorous as possible, and thus more attractive to females, but this possible function has not been described.</p>
<p>Self-medication is used by some animals to prevent infections, fight against parasites or pathogens, or improve a variety of suboptimal physiological states <xref ref-type="bibr" rid="pone.0111057-deRoode1">[7]</xref>–<xref ref-type="bibr" rid="pone.0111057-Forbey1">[12]</xref>. It involves the consumption of plants, animals or minerals that have a prophylactic or therapeutic effect against disease agents which should lead to an increase in host fitness <xref ref-type="bibr" rid="pone.0111057-deRoode1">[7]</xref>–<xref ref-type="bibr" rid="pone.0111057-Hart1">[10]</xref>.</p>
<p>In this context, we examine here the consumption of poisonous insects by great bustards (<italic>Otis tarda</italic>) during the mating season. The great bustard is a polygynous bird with one of the most strongly skewed male mating success values among birds <xref ref-type="bibr" rid="pone.0111057-Alonso1">[13]</xref>. Males gather each year at traditional arenas (<italic>leks</italic>) where they perform elaborate sexual exhibitions directed towards females in order to express their status and condition <xref ref-type="bibr" rid="pone.0111057-Alonso2">[14]</xref>. Great bustards are among the few birds that feed on blister beetles (<italic>Berberomeloe majalis</italic>, <italic>Physomeloe corallifer</italic>; <xref ref-type="fig" rid="pone-0111057-g001">Fig. 1</xref>) <xref ref-type="bibr" rid="pone.0111057-Lane1">[15]</xref>,<xref ref-type="bibr" rid="pone.0111057-Rocha1">[16]</xref>. These insects are avoided by most animals because they contain cantharidin, a bitter-tasting and highly toxic defensive chemical with high immunogenicity <xref ref-type="bibr" rid="pone.0111057-Dettner1">[17]</xref> that acts in blister beetles as fungicide and nematocide <xref ref-type="bibr" rid="pone.0111057-Carrel1">[18]</xref>. Only a few species, such as the spur-winged goose <italic>Plectropterus gambensis</italic> and the northern leopard frogs <italic>Rana pipiens</italic> consume blister beetles, with a likely side effect of becoming toxic to predators <xref ref-type="bibr" rid="pone.0111057-Bartram1">[19]</xref>,<xref ref-type="bibr" rid="pone.0111057-Eisner1">[20]</xref>. Cantharidin is also a well-known aphrodisiac compound that was obtained in the past from a beetle known as Spanish fly <xref ref-type="bibr" rid="pone.0111057-Sandroni1">[21]</xref>. In humans, it causes priapism in men and pelvic congestion in women <xref ref-type="bibr" rid="pone.0111057-Moed1">[22]</xref>,<xref ref-type="bibr" rid="pone.0111057-Nickolls1">[23]</xref>. Cantharidin-tolerant foragers would benefit from its anti-microbial and anthelminthic properties <xref ref-type="bibr" rid="pone.0111057-Campbell1">[24]</xref>, and thus could enhance their health and their attractiveness to potential mates during the mate selection process.</p>
<fig id="pone-0111057-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0111057.g001</object-id><label>Figure 1</label><caption>
<title>Blister beetles <italic>Berberomeloe majalis</italic> (A) and <italic>Physomeloe corallifer</italic> (D) displaying their characteristic red and black aposematic warning coloration.</title>
<p>Adults only appear in spring. Defensive reactions in <italic>B</italic>. <italic>majalis</italic> (B, C) involve thanatosis and autohaemorrhea, with large droplets of red haemolymph containing cantharidin expelled through thoracic and limb segment joints. Red tegumentary protuberances in the thorax of <italic>P. corallifer</italic> (E) permanently mimic haemolymph droplets. Photographs: M. García-París.</p>
</caption><graphic mimetype="image" xlink:href="info:doi/10.1371/journal.pone.0111057.g001" position="float" xlink:type="simple"/></fig>
<p>We hypothesized that cantharidin, obtained by consumption of blister beetles, has a direct effect as an oral anti-microbial and pathogen-limiting compound. To examine this hypothesis, we combined <italic>in vitro</italic> experiments and field observations. In particular, we (i) evaluated the presence of potential agents of sexually transmitted diseases in male and female faecal samples during the mating season; (ii) tested whether pathogenic microorganisms from bustard faeces were reduced in experimental cultures exposed to dose-dependent cantharidin, and (iii) checked for presence of cantharidin in kidney and liver of great bustards, where it could be excreted or stored. Fulfilment of these evidences would not definitely demonstrate a possible self-medication function of cantharidin in great bustards, but should at least provide partial support for it.</p>
<p>Another appealing hypothesis is that consumption of blister beetles could also enhance the attractiveness of males to females by reducing their parasite load. This would counter-balance the risk of cantharidin poisoning <xref ref-type="bibr" rid="pone.0111057-SnchezBarbudo1">[25]</xref> with reproductive benefits. In addition, males could benefit from cantharidin during the display season, when they may be particularly vulnerable to infections due to their presumably depressed immune system associated to their strenuous investment in sexual display <xref ref-type="bibr" rid="pone.0111057-Folstad1">[26]</xref>,<xref ref-type="bibr" rid="pone.0111057-Roberts1">[27]</xref>. This hypothesis predicts a male-biased blister beetle consumption (i.e., both sexes would consume blister beetles as a healing or prophylactic agent, but males with potential to attract females should consume them also to enhance their body condition and mating opportunities). Previous studies have shown that age, weight, and display effort are the main predictors of male mating success in this species, and that whiskers and neck plumage are reliable indicators of male age and weight <xref ref-type="bibr" rid="pone.0111057-Alonso1">[13]</xref>,<xref ref-type="bibr" rid="pone.0111057-Alonso2">[14]</xref>. However, the function of the meticulous inspection females perform of the male's cloaca prior to copulation is still unknown. A male's sexual display consists of a series of extravagant body postures and movements that end, when an interested female approaches, with a reiterative and almost obstinate exhibition of the cloaca, which is fully surrounded by pure white feathers that allow an easy detection of possible parasites or their remains <xref ref-type="bibr" rid="pone.0111057-Nakamura1">[28]</xref>,<xref ref-type="bibr" rid="pone.0111057-Sheldon1">[29]</xref> (<xref ref-type="fig" rid="pone-0111057-g002">Fig. 2 A–C</xref>). Female bustards are extremely choosy during mating, accepting on average only one of ten males attempting to copulate with them <xref ref-type="bibr" rid="pone.0111057-Alonso1">[13]</xref>. Females closely inspect the male's cloaca, pecking around it as an essential part of mate selection, after which they decide whether to mate or not with the male <xref ref-type="bibr" rid="pone.0111057-Hellmich1">[30]</xref>,<xref ref-type="bibr" rid="pone.0111057-HidalgodeTrucios1">[31]</xref>. Our interpretation of this cloaca pecking is that the female is checking the health status of the male in order to select a copulation partner with genes for tolerance to cantharidin, the consumption of which allows for resistance to pathogens <xref ref-type="bibr" rid="pone.0111057-Hamilton1">[2]</xref>, while avoiding immediate disease transmission during copulations (<xref ref-type="fig" rid="pone-0111057-g002">Fig. 2 D–E</xref>) <xref ref-type="bibr" rid="pone.0111057-Kirkpatrick1">[4]</xref>,<xref ref-type="bibr" rid="pone.0111057-Lombardo1">[32]</xref>,<xref ref-type="bibr" rid="pone.0111057-Birkhead1">[33]</xref>. Similarly, by showing the underwing feather shafts and other parts of their ventral plumage (<xref ref-type="fig" rid="pone-0111057-g002">Fig. 2C</xref>), which are preferred sites for adults and eggs of mites and lice to stay or lay eggs, great bustards might show to visiting females that they are free of these ectoparasites.</p>
<fig id="pone-0111057-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0111057.g002</object-id><label>Figure 2</label><caption>
<title>A great bustard male in courtship display shows its cloaca to prospective females (A–B; Photographs: C. Palacín).</title>
<p>Detail of white feathers of cloaca during display of male great bustard (C; Photographs: F. Kovacs). Detail of other bird with an <italic>Otiditaenia conoides</italic> in extending from the cloaca (D; Photograph: C. Palacín). Detail of a <italic>Otiditaenia conoides</italic> individual separate from the body bird (E; Photograph: A. Lucas). Black arrows show the cloaca position.</p>
</caption><graphic mimetype="image" xlink:href="info:doi/10.1371/journal.pone.0111057.g002" position="float" xlink:type="simple"/></fig>
<p>Following this argument, males should consume an optimal amount of blister beetles to prevent and fight against pathogens that affect their health and/or act as agents of sexually-transmitted diseases. The benefits for health of feeding on blister beetles must be balanced against the toxicity of cantharidin, because an excess of blister beetles in the diet can be lethal for a great bustard <xref ref-type="bibr" rid="pone.0111057-SnchezBarbudo1">[25]</xref>. To test this, we examined the potential of blister beetles to act against pathogens responsible for sexually transmitted diseases, and any preference for blister beetles by great bustard males. These additional analyses involved (i) determination of cantharidin concentration in blister beetles collected in the study area; and (ii) testing whether there is a sex bias in blister beetle consumption by male and female great bustards.</p>
</sec><sec id="s2" sec-type="methods">
<title>Methods</title>
<sec id="s2a">
<title>Ethical statement</title>
<p>Sampling permits for insect collection and research permits for dead bustard transfer to the Museo Nacional de Ciencias Naturales were provided by the environmental conservation departments of Madrid and Castilla-La Mancha regions (Consejería de Medio Ambiente y Ordenación del Territorio de la Comunidad de Madrid, and Consejería de Agricultura y Medio Ambiente de la Junta de Comunidades de Castilla-La Mancha). The research was carried out on private lands with landowner permission and on the SPA 139 ‘Estepas Cerealistas de los Ríos Jarama y Henares’ with permit from the Consejería de Medio Ambiente y Ordenación del Territorio de la Comunidad de Madrid. Sampling protocols were approved by National R+D Agencies (Ministerio de Economía y Competitividad).</p>
</sec><sec id="s2b">
<title>Study area</title>
<p>The study was conducted in Madrid province, central Spain (40°43′N-40°00′N, 3°29′W-3°06′W; for details see <xref ref-type="bibr" rid="pone.0111057-Alonso3">[34]</xref>), where a population of ca. 1500 great bustards has been intensively studied over two decades <xref ref-type="bibr" rid="pone.0111057-Martn1">[35]</xref>. Great bustards live in dry cereal farmland, where blister beetles <italic>Berberomeloe majalis</italic> and <italic>Physomeloe corallifer</italic> are also found. Great bustards are omnivorous, their diet consisting of green plant material, arthropods and seeds. Arthropods represent a much larger component in summer than in spring (respectively, ca. 40% and 7%) <xref ref-type="bibr" rid="pone.0111057-Lane1">[15]</xref>,<xref ref-type="bibr" rid="pone.0111057-Lucio1">[36]</xref>,<xref ref-type="bibr" rid="pone.0111057-Palacios1">[37]</xref>. Great bustards are distributed in fragmented populations through the Palaearctic from the Iberian Peninsula and Morocco to Eastern China <xref ref-type="bibr" rid="pone.0111057-Alonso4">[38]</xref>, whilst <italic>B. majalis</italic> is present in north-western Africa, Iberian Peninsula and southern France <xref ref-type="bibr" rid="pone.0111057-Bologna1">[39]</xref>,<xref ref-type="bibr" rid="pone.0111057-GarcaParis1">[40]</xref>, and <italic>P. corallifer</italic> is endemic to the Iberian Peninsula <xref ref-type="bibr" rid="pone.0111057-Bologna1">[39]</xref>.The season when adult blister beetles are present coincides with the display period of great bustards. Blister beetles of genus <italic>Berberomeloe</italic> are diurnal, widespread, but not particularly abundant <xref ref-type="bibr" rid="pone.0111057-GarcaParis1">[40]</xref>,<xref ref-type="bibr" rid="pone.0111057-GarcaPars1">[41]</xref>. These two blister beetle species are among the largest of the family Meloidae (<xref ref-type="fig" rid="pone-0111057-g001">Fig. 1</xref>) <xref ref-type="bibr" rid="pone.0111057-GarcaParis1">[40]</xref>. Furthermore, the average body size of these species is higher than that of all other beetles in the study area (20–65 mm).</p>
<p>Cereal (mainly wheat and barley) is usually grown in a traditional two-year rotation system. Some fields are left as fallows for two or more years. Ploughed and sown fields, stubbles, fallows, and a few fields with olives, legumes, sunflower and grape vines create a dynamic mosaic which, complemented with occasional sheep grazing in stubbles and fallows, constitutes an optimal habitat for great bustards <xref ref-type="bibr" rid="pone.0111057-Lane2">[42]</xref>. A thorough description of the study area is available in previous studies <xref ref-type="bibr" rid="pone.0111057-Lane2">[42]</xref>,<xref ref-type="bibr" rid="pone.0111057-Magaa1">[43]</xref>.</p>
</sec><sec id="s2c">
<title>Determination of microflora, parasites and pathogens in great bustard faeces</title>
<p>We used a sample of 53 faeces (25 males, 28 females) collected at dawn in April 2012 to search for organisms of the common digestive flora, as well as for known parasites and pathogens. These included protozoa, helminths and bacteria. Fresh faeces were sampled with sterile microbiological swabs that were subsequently inserted into tubes containing Amies transport medium. The tubes were transported to the laboratory in a container with ice within the same day of collection and processed within the two following hours. Swaps were cultured on 5% sheep blood agar, Vogel Johnson agar, COBA, and MacConkey agar (Oxoid, UK). Plates were incubated at 37°C, either in normal atmospheric or anaerobic environment (just in the case of sheep blood agar plates) (bioMérieux, France) conditions respectively, for 24 h. All potential colonies were then subcultured on appropriate medium and identified to species level by using multi-substrate identification stripes (API 20 E, API 20NE, API STAPH; bioMérieux, France), Lancefield grouping kit (Oxoid, UK) and other biochemical tests (catalase; bioMérieux, and oxidase; BD). Criteria for the presumptive identification of bacterial pathogens were: colonial morphology and colour; presence or absence of haemolysis on blood agar; appearance when stained by the gram method; ability to grow on MacConkey agar; and reactions in catalase and oxidase tests <xref ref-type="bibr" rid="pone.0111057-Quinn1">[44]</xref>.</p>
<p>Coccidian species (Protozooa) were examined by oocyst sporulation with 2.5% potassium dichromate during fourteen days followed by zinc sulfate flotation. For the detection of helminth eggs in faeces (trematodes, acantocephalans, cestodes and nematodes) we used the flotation method with zinc sulphate solution as well as slide direct examination procedure.</p>
</sec><sec id="s2d">
<title>Test of the bactericidal activity of cantharidin on great bustard faecal microflora</title>
<p>To test the bactericidal activity of cantharidin, we firstly determined the faecal microflora in great bustards and secondly, the susceptibility of bacteria species to cantharidin by exposing them to different cantharidin concentrations in cultures. We also determined the cantharidin content (mg/g) in 63 individuals of the two blister beetle species present in our study area (<italic>Berberomeloe majalis</italic> and <italic>Physomeloe corallifer</italic>) by gas chromatography.</p>
<sec id="s2d1">
<title>Susceptibility of great bustard's faecal microflora to cantharidin</title>
<p>We isolated six bacteria species (<italic>Escherichia coli, Staphylococcus aureus, S. lentus, Bacillus</italic> sp, <italic>Clostridium</italic> sp, <italic>Kocuria</italic> sp) with distinctive colony morphology, re-streaked them on nutrient agar, a generalized medium, followed by incubation at 37°C for 48–72 hours until we could determine the purity of our cultures as evidenced by unique colony morphology characteristics (gram staining, Panreac).</p>
<p>The bactericidal activity of cantharidin over the isolated species was tested according to Kirby-Bauer disc diffusion method <xref ref-type="bibr" rid="pone.0111057-Institute1">[45]</xref> on Mueller Hinton agar (Oxoid). Cantharidin (Sigma–Aldrich Laboratories, St. Louis, MO) was dissolved in distilled water at five concentrations (1∶10, 1∶100, 1∶1000, 1∶10000, 1∶100000). Six sterilized discs (bioMérieux, France) were impregnated with these five concentrations and a control with only distilled water. Between 5 and 10 culture plates with six discs were performed per each bacteria species, totalling 45 cultures. After 24 hours, the plates with the discs were removed from the incubator and growth-inhibition halos were measured (mm). Sensitive percentage per cantharidin concentration was calculated as percentage of plates with an inhibition effect (&gt; 6 mm) of cantharidin.</p>
</sec><sec id="s2d2">
<title>Cantharidin concentration in blister beetles</title>
<p>Cantharidin content was analysed in dried individuals of <italic>B. majalis</italic> (<italic>n</italic> = 58) and <italic>P. corallifer</italic> (<italic>n</italic> = 5). Body dry weight was determined after 48 hours of freeze drying in a VIRTIS 25 LE-53 GENESIS SQ. The lyophilized individuals (0.057–0.725 g) were transferred to a teflon digestion reactor with hydrochloric acid (3 ml) and were heated to 120°C during three hours. Once the reactor was cooled, its content was added to chloroform (5 ml), followed by mixing in a vortex for 15 seconds. Once phases were separated, lower layer of chloroform was collected by Pasteur pipette. This extraction was repeated twice to ensure complete extraction. Samples were placed in 2 ml vials for analysis by gas chromatography coupled to mass spectrometry (GC–MS). The GC–MS system consisted of a Thermo Finnigan Trace GC 2000 coupled with a Trace MS mass selective detector. The chromatographic conditions were controlled using the Xcalibur software version 1.2 (Thermo Finnigan; San José, CA). The GC column was a SLB-5 ms (30 m×0.32 mm, 0.25 µm, Supelco Analytical; Bellefonte, PA). The flow rate of helium was 0.8 ml/min. The injection volume was 1 µl in splitless mode for 2 min. Injector conditions were 250°C in constant flow mode. The column oven had an initial temperature of 50°C for two minutes. The subsequent temperature was programmed at heating rate of 10°C/min to 310°C. The final temperature was held isothermally for 5 min. Total run time was 30 min. Cantharidin detection was performed by selected ion monitoring (SIM), registering m/z = 128 that is the majority ion of cantharidin's mass spectra. Cantharidin identification was performed by comparison with mass spectra available in NIST MS search 2.0 library. Limit of detection (LOD) was 1.5 µg/mg. Confirmation and quantification was achieved with the retention time and calibration curves (range: 0.015–48 µg/ml, slope: 476.401, r<sup>2</sup> = 0.999) obtained from the injection of cantharidin standard purchased from Sigma–Aldrich (St. Louis, MO).</p>
</sec></sec><sec id="s2e">
<title>Cantharidin in kidney and liver of great bustards: post-mortem analyses</title>
<p>We searched for cantharidin in 18 livers and 19 kidneys of 25 great bustard carcasses, collected during 1999–2013 (winter: n = 8; spring: n = 10; summer: n = 3, fall: n = 4). All individuals were adults, we distinguished males (<italic>n</italic> = 14) and females (<italic>n</italic> = 11). In addition, the stomach content of spring season carcasses (<italic>n</italic> = 15) was analysed to check blister beetle ingestion. As in faeces analysis, the number and biomass of ingested blister beetles was estimated.</p>
<p>Cantharidin could be found in kidneys because it is eliminated by renal excretion <xref ref-type="bibr" rid="pone.0111057-SnchezBarbudo1">[25]</xref>,<xref ref-type="bibr" rid="pone.0111057-Ray1">[46]</xref>. On the other hand, great bustards may prolong the pathogenic activity of cantharidin over periods longer than the few hours after blister beetle ingestion if cantharidin is accumulated in the liver, because from the liver it may be slowly disseminated <xref ref-type="bibr" rid="pone.0111057-Bailey1">[47]</xref>. Cantharidin is a lipophilic substance that may cross the gastrointestinal epithelium and therefore accumulate in the liver, although to our knowledge this has not been reported in great bustards. It must be acknowledged that the sole presence of cantharidin in liver or kidney does not definitely demonstrate that great bustards manage the use of this chemical at their will, because presence in these organs may be an intermediate step of a detoxification procedure. A long-lasting therapeutic use of cantharidin by great bustards will only be demonstrated if a body reservoir for this compound is discovered.</p>
<p>We used a modified chromatography method described in Sánchez-Barbudo et al. <xref ref-type="bibr" rid="pone.0111057-SnchezBarbudo1">[25]</xref>. For calibrating and validating the method, five chicken liver samples spiked with 0.023 mg/g of cantharidin were previously analyzed. The calculated recovery was 52.24±7.15%, so the detected concentrations were corrected by this value.</p>
<p>Liver and kidney samples (1 g frozen sample or 200 mg lyophilized sample) were transferred to a mortar and homogenised with 9 g of anhydrous sodium sulphate. The homogenate was placed in a beaker, then 15 ml of dichloromethane were added and the samples were shaken for 10 min in a horizontal shaker (IKA KS 250), followed by 5 min of sonication (Elma transonic 700). The extracts were filtered and the remaining sample homogenates were extracted twice again with 5 ml of dichloromethane. The extracts were pooled, evaporated in a helium stream, adjusted to 1 ml and placed in a 2 ml vial for GC analysis.</p>
<p>The GC–MS system consisted of a Hewlett Packard 5890A II GC coupled with a Hewlett Packard 5972 mass selective detector. The chromatographic conditions were controlled using the MS Chemstation software version C.01.05 (Hewlett Packard; Palo Alto, CA).</p>
<p>The GC column was a 007-5MS (30 m×0.32 mm, 0.25 µm, Quadrex Corporation; Woodbridge, CT). The flow rate of helium was 0.8 ml/min. The injection volume was 1 µl in splitless mode for 2 min. Injector conditions were 250°C in constant flow mode. The column oven had an initial temperature of 50°C for two minutes. The subsequent temperature was programmed at heating rate of 10°C/min to 310°C. The final temperature was held isothermally for 5 min. Total run time was 30 min. Cantharidin detection was performed by selected ion monitoring (SIM), registering m/z  =  128 that is the majority ion of cantharidin's mass spectra.</p>
<p>Cantharidin identification was performed by comparison with mass spectra available in NIST/EPA/NIH 75k Mass Spectral Database. Limit of detection (LOD) was 1.44 µg/mg. Confirmation and quantification was achieved with the retention time and calibration curves (range: 0.015–48 µg/ml, slope: 476.401, r<sup>2</sup> = 0.999) obtained from the injection of cantharidin standard purchased from Sigma–Aldrich (St. Louis, MO). When cantharidin concentration was below LOD, we divided LOD by two, as a valid way of dealing with samples reported to contain non-detectable values of cantharidin <xref ref-type="bibr" rid="pone.0111057-Finkelstein1">[48]</xref>.</p>
</sec><sec id="s2f">
<title>Selection of blister beetles by great bustards</title>
<sec id="s2f1">
<title>Abundance of blister beetles in the field</title>
<p>Blister beetles were sampled in 2006 and 2007 between March and May. An observer walked slowly at constant speed (0.04 ms<sup>−1</sup>) and counted the blister beetles along 30×1 m transects. Abundance was calculated as the number of beetles per hectare (n/ha) in each of the main six substrate types: sprouted cereal (wheat <italic>Triticum aestivum</italic> and barley <italic>Hordeum vulgare</italic>), legumes (vetch <italic>Vicia sativa</italic>, <italic>Vicia ervilia</italic>, lentils <italic>Lens squlenta</italic>, peas <italic>Pisum sativum</italic>, and chickpeas <italic>Cicer arietinum</italic>), ploughed grounds, stubbles, fallows, and borders between fields. Other invertebrates consumed by great bustards were also counted to calculate the relative abundance of blister beetles. Between 14 and 33 transects were performed in each substrate type (194 transects in total), covering most areas with great bustards in Madrid.</p>
<p>The abundance of blister beetles in the study area was calculated multiplying the mean abundance in each substrate type by its relative surface. Relative surfaces of all substrate types were calculated from a sample of 3222 fields. The surface of borders was calculated by multiplying the total length of the border by 0.3 m (average width of borders in all areas).</p>
<p>Blister beetle dry biomass was calculated through published length-weight equations <xref ref-type="bibr" rid="pone.0111057-Hdar1">[49]</xref>. Maximum body length was measured in all adult arthropods with a digital calliper (0.01 mm precision). Dry weights were estimated for each arthropod family in each transect and substrate type by means of linear regressions of body weight on body length. Dry biomasses of other invertebrates consumed by great bustards were also calculated to determine the relative biomass of blister beetles in the field. Available biomass was computed as grams of beetle per hectare (g/ha) and grams of other invertebrates per hectare (g/ha), for each substrate type.</p>
</sec><sec id="s2f2">
<title>Abundance of blister beetles in great bustard faeces</title>
<p>Two hundred and twelve fresh faeces were collected between March and May (2006–2007), with a total dry weight of 556.6 g. Faeces were collected at roosting sites at dawn (70.8% of all droppings, 82.8% g dry weight), at midday resting sites (26%, 16% g dry weight), and also throughout the day on fields where great bustard flocks had been previously located. Since males and females of this species live in separate flocks, faeces from both sexes could be easily distinguished <xref ref-type="bibr" rid="pone.0111057-Bautista1">[50]</xref>. Sexual segregation does not imply habitat segregation in this species (pers. obs.). Faeces collected at roosting and resting sites accounted for the total intake, respectively in the previous evenings and mornings. Only complete faeces were collected, and each faeces was stored in a single plastic bag, recording the date, coordinates and substrate type of the collection site, and sex of the flock. The faeces were weighted (accuracy: 0.001 g) after drying them during 48 h at 60°C. The number of blister beetles in faeces was estimated by counting the remains of characteristic exoskeleton components such as pronotum, head or elytra. The biomass of blister beetles consumed was calculated using published length-weight equations <xref ref-type="bibr" rid="pone.0111057-Hdar2">[51]</xref>. The same method was used to estimate the abundance (n) and biomass (g) of other insect taxa consumed (Acrididae, Formicidae, Hemiptera, Curculionidae, Scarabeidae and Tenebrionidae).</p>
</sec></sec><sec id="s2g">
<title>Data analyses</title>
<p>Determination of faecal microflora, parasites and pathogens in great bustard faeces was performed with multiple comparisons (Fisher's exact probability tests and Odds ratio) on the prevalence of the different organisms between male and female great bustards. Therefore, the sequential Bonferroni correction was used to adjust the significance level to control for type I errors.</p>
<p>Differences in bactericide activity with increasing cantharidin concentration (control, 1∶100000, 1∶10000, 1∶1000, 1∶100 and 1∶10) were analysed by generalized linear models (GLZs) with logit link function on six isolated bacteria species (<italic>Escherichia coli, Staphylococcus aureus, S. lentus, Bacillus</italic> sp, <italic>Clostridium</italic> sp, <italic>Kocuria</italic> sp).</p>
<p>The presence of cantharidin in liver and kidneys was analyzed by GLZs with logit link function, where presence/absence of cantharidin was the dependent variable and sex was a fixed factor. Statistical significance of sex differences in concentration of cantharidin in liver and in kidneys was analysed with a two-tailed Wilcoxon Kruskall-Wallis Rank test.</p>
<p>Sex differences in the frequency occurrence of faeces and stomachs content with and without blister beetles were calculated with Chi-square test. Sex differences in mean abundance (n) and biomass (g) of blister beetles per faeces and per gram of faeces were calculated with GLZ (Poisson distribution). The proportion of blister beetles with respect to all invertebrates identified in faeces was calculated as both, relative abundance (n%) and relative biomass (g%). The effects of sex and year on the mean proportions (relative abundance and relative biomass) of blister beetles in faeces were calculated with GLZ (binomial distribution and a logit link function). Statistics calculated with pooled data from males and females' faeces were weighted by the relative population size of each sex, since there is a high adult sex-ratio bias in this population: 2.42 females per male <xref ref-type="bibr" rid="pone.0111057-Alonso3">[34]</xref>.</p>
<p>The abundance of invertebrates in the habitat of great bustards was calculated using transects, and compared to the proportion of blister beetles in excrements. Differences between proportions available and consumed by sex and year were tested first with a simple comparison of proportions, and subsequently with a GLZ (binomial distribution and a logit link function). Differences in blister beetle body sizes (g) available and consumed by sex and year were also analysed with Student's t-test. All analyses above were carried out in R version 2.15.1 <xref ref-type="bibr" rid="pone.0111057-R1">[52]</xref>.</p>
</sec></sec><sec id="s3">
<title>Results</title>
<sec id="s3a">
<title>Microflora, parasites and pathogens in great bustard faeces</title>
<p>Eight species of bacteria were found in faecal samples (<xref ref-type="table" rid="pone-0111057-t001">Table 1</xref>). There were no differences in the prevalence of cantharidin-sensitive bacteria in faeces of males and females (<xref ref-type="table" rid="pone-0111057-t001">Table 1</xref>), with the exception of <italic>Kocuria</italic> sp and <italic>S. intermedius</italic>. Prevalence of Coccidia was 3.77% in faecal samples, whilst 71.7% of faecal samples presented intestinal parasitism by nematodes (64.2%) and cestodes (20.8%, <xref ref-type="table" rid="pone-0111057-t001">Table 1</xref>). These nematodes were identified as <italic>Trichostrongylus tennuis</italic>, <italic>Ascaridia</italic> sp, <italic>Heterakis isolenche</italic> and <italic>Capillaria</italic> sp; and cestodes as <italic>Otiditaenia conoides</italic>. There were no sexual differences in the prevalence of coccidia, nematodes and cestodes (<xref ref-type="table" rid="pone-0111057-t001">Table 1</xref>).</p>
<table-wrap id="pone-0111057-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0111057.t001</object-id><label>Table 1</label><caption>
<title>Prevalence of each pathogenic organism in male and female great bustards, estimated as percentage of faeces with pathogen presence.</title>
</caption><alternatives><graphic id="pone-0111057-t001-1" position="float" mimetype="image" xlink:href="info:doi/10.1371/journal.pone.0111057.t001" xlink:type="simple"/>
<table><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup>
<thead>
<tr>
<td align="left" rowspan="1" colspan="1">Pathogens</td>
<td align="left" rowspan="1" colspan="1">Males (<italic>n</italic> = 25)</td>
<td align="left" rowspan="1" colspan="1">Females (<italic>n</italic> = 28)</td>
<td align="left" rowspan="1" colspan="1">χ<sup>2</sup></td>
<td align="left" rowspan="1" colspan="1"><italic>P</italic></td>
<td align="left" rowspan="1" colspan="1">Odds ratio</td>
<td align="left" rowspan="1" colspan="1">Upper CI</td>
<td align="left" rowspan="1" colspan="1">Lower CI</td>
</tr>
</thead>
<tbody>
<tr>
<td align="left" rowspan="1" colspan="1">Bacteria</td>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"><italic>E. coli</italic></td>
<td align="left" rowspan="1" colspan="1">100.0</td>
<td align="left" rowspan="1" colspan="1">100.0</td>
<td align="left" rowspan="1" colspan="1">0.00</td>
<td align="left" rowspan="1" colspan="1">1.00</td>
<td align="left" rowspan="1" colspan="1">1.000</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"><italic>Staphylococcus intermedius</italic></td>
<td align="left" rowspan="1" colspan="1">12.0</td>
<td align="left" rowspan="1" colspan="1">60.7</td>
<td align="left" rowspan="1" colspan="1">13.34</td>
<td align="left" rowspan="1" colspan="1">0.00</td>
<td align="left" rowspan="1" colspan="1">0.088</td>
<td align="left" rowspan="1" colspan="1">0.021</td>
<td align="left" rowspan="1" colspan="1">0.367</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"><italic>Staphylcoccus aureus</italic></td>
<td align="left" rowspan="1" colspan="1">0.0</td>
<td align="left" rowspan="1" colspan="1">7.1</td>
<td align="left" rowspan="1" colspan="1">1.86</td>
<td align="left" rowspan="1" colspan="1">0.27</td>
<td align="left" rowspan="1" colspan="1">0.929</td>
<td align="left" rowspan="1" colspan="1">0.838</td>
<td align="left" rowspan="1" colspan="1">1.029</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"><italic>Staphylcoccus lentus</italic></td>
<td align="left" rowspan="1" colspan="1">20.0</td>
<td align="left" rowspan="1" colspan="1">0.0</td>
<td align="left" rowspan="1" colspan="1">6.18</td>
<td align="left" rowspan="1" colspan="1">0.02</td>
<td align="left" rowspan="1" colspan="1">1.250</td>
<td align="left" rowspan="1" colspan="1">1.028</td>
<td align="left" rowspan="1" colspan="1">1.521</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"><italic>Bacillus</italic> sp</td>
<td align="left" rowspan="1" colspan="1">8.0</td>
<td align="left" rowspan="1" colspan="1">3.6</td>
<td align="left" rowspan="1" colspan="1">0.49</td>
<td align="left" rowspan="1" colspan="1">0.46</td>
<td align="left" rowspan="1" colspan="1">2.348</td>
<td align="left" rowspan="1" colspan="1">0.200</td>
<td align="left" rowspan="1" colspan="1">27.592</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"><italic>Clostridium</italic> sp</td>
<td align="left" rowspan="1" colspan="1">0.0</td>
<td align="left" rowspan="1" colspan="1">10.7</td>
<td align="left" rowspan="1" colspan="1">2.84</td>
<td align="left" rowspan="1" colspan="1">0.14</td>
<td align="left" rowspan="1" colspan="1">0.893</td>
<td align="left" rowspan="1" colspan="1">0.785</td>
<td align="left" rowspan="1" colspan="1">1.015</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"><italic>Kocuria</italic> sp</td>
<td align="left" rowspan="1" colspan="1">56.0</td>
<td align="left" rowspan="1" colspan="1">0.0</td>
<td align="left" rowspan="1" colspan="1">21.40</td>
<td align="left" rowspan="1" colspan="1">0.00</td>
<td align="left" rowspan="1" colspan="1">2.273</td>
<td align="left" rowspan="1" colspan="1">1.460</td>
<td align="left" rowspan="1" colspan="1">3.537</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"><italic>Acinetobacter</italic> sp</td>
<td align="left" rowspan="1" colspan="1">4.0</td>
<td align="left" rowspan="1" colspan="1">0.0</td>
<td align="left" rowspan="1" colspan="1">1.14</td>
<td align="left" rowspan="1" colspan="1">0.47</td>
<td align="left" rowspan="1" colspan="1">1.042</td>
<td align="left" rowspan="1" colspan="1">0.962</td>
<td align="left" rowspan="1" colspan="1">1.128</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Coccidia</td>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"><italic>Caryospora</italic> sp</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">3.6</td>
<td align="left" rowspan="1" colspan="1">0.91</td>
<td align="left" rowspan="1" colspan="1">0.53</td>
<td align="left" rowspan="1" colspan="1">0.964</td>
<td align="left" rowspan="1" colspan="1">0.898</td>
<td align="left" rowspan="1" colspan="1">1.036</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"><italic>Eimeria</italic> sp</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">3.6</td>
<td align="left" rowspan="1" colspan="1">0.91</td>
<td align="left" rowspan="1" colspan="1">0.53</td>
<td align="left" rowspan="1" colspan="1">0.964</td>
<td align="left" rowspan="1" colspan="1">0.898</td>
<td align="left" rowspan="1" colspan="1">1.036</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Nematode</td>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"><italic>Trichostrongylus tennuis</italic></td>
<td align="left" rowspan="1" colspan="1">68.0</td>
<td align="left" rowspan="1" colspan="1">53.6</td>
<td align="left" rowspan="1" colspan="1">1.15</td>
<td align="left" rowspan="1" colspan="1">0.22</td>
<td align="left" rowspan="1" colspan="1">1.842</td>
<td align="left" rowspan="1" colspan="1">0.600</td>
<td align="left" rowspan="1" colspan="1">5.653</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"><italic>Ascaridia</italic> sp</td>
<td align="left" rowspan="1" colspan="1">0.0</td>
<td align="left" rowspan="1" colspan="1">3.6</td>
<td align="left" rowspan="1" colspan="1">0.91</td>
<td align="left" rowspan="1" colspan="1">0.53</td>
<td align="left" rowspan="1" colspan="1">0.964</td>
<td align="left" rowspan="1" colspan="1">0.898</td>
<td align="left" rowspan="1" colspan="1">1.036</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Morphotype 1<xref ref-type="table-fn" rid="nt101">†</xref></td>
<td align="left" rowspan="1" colspan="1">36.0</td>
<td align="left" rowspan="1" colspan="1">17.9</td>
<td align="left" rowspan="1" colspan="1">2.24</td>
<td align="left" rowspan="1" colspan="1">0.12</td>
<td align="left" rowspan="1" colspan="1">2.588</td>
<td align="left" rowspan="1" colspan="1">0.730</td>
<td align="left" rowspan="1" colspan="1">9.175</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Morphotype 2<xref ref-type="table-fn" rid="nt101">†</xref></td>
<td align="left" rowspan="1" colspan="1">4.0</td>
<td align="left" rowspan="1" colspan="1">3.6</td>
<td align="left" rowspan="1" colspan="1">0.01</td>
<td align="left" rowspan="1" colspan="1">0.73</td>
<td align="left" rowspan="1" colspan="1">1.125</td>
<td align="left" rowspan="1" colspan="1">0.067</td>
<td align="left" rowspan="1" colspan="1">18.984</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Cestoda</td>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"><italic>Otiditaenia conoides</italic></td>
<td align="left" rowspan="1" colspan="1">0.0</td>
<td align="left" rowspan="1" colspan="1">7.1</td>
<td align="left" rowspan="1" colspan="1">1.86</td>
<td align="left" rowspan="1" colspan="1">0.27</td>
<td align="left" rowspan="1" colspan="1">0.929</td>
<td align="left" rowspan="1" colspan="1">0.838</td>
<td align="left" rowspan="1" colspan="1">1.029</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Undetermined<xref ref-type="table-fn" rid="nt102">††</xref></td>
<td align="left" rowspan="1" colspan="1">24.0</td>
<td align="left" rowspan="1" colspan="1">10.7</td>
<td align="left" rowspan="1" colspan="1">1.65</td>
<td align="left" rowspan="1" colspan="1">0.18</td>
<td align="left" rowspan="1" colspan="1">1.654</td>
<td align="left" rowspan="1" colspan="1">0.582</td>
<td align="left" rowspan="1" colspan="1">11.899</td>
</tr>
</tbody>
</table>
</alternatives><table-wrap-foot><fn id="nt101"><label>†</label><p><italic>Heterakis isolenche</italic> and <italic>Capillaria sp</italic>. Prevalence of these species was calculated, but cannot be with certainty assigned.</p></fn><fn id="nt102"><label>††</label><p>Hispaniolepis villosa or Idiogenes otidis.</p></fn></table-wrap-foot></table-wrap></sec><sec id="s3b">
<title>Bactericidal activity of cantharidin on great bustard faecal microflora</title>
<p>Bactericide activity of cantharidin varied among species (<italic>F</italic><sub>5,228</sub> = 17.63, <italic>P</italic>&lt;0.001) and increased with concentration (<italic>F</italic><sub>5,258</sub> = 17.37, <italic>P</italic>&lt;0.001, <xref ref-type="fig" rid="pone-0111057-g003">Fig. 3</xref>). Sensitivity (% of samples with an inhibition effect &gt; 6 mm of cantharidin) and inhibition halo were higher in <italic>Bacillus</italic> sp, <italic>Clostridium</italic> sp and <italic>Kocuria</italic> sp. <italic>Bacillus</italic> sp was the most sensitive species at any concentration, and <italic>Staphylococcus lentus</italic> the most resistant at all concentrations (<xref ref-type="fig" rid="pone-0111057-g003">Fig. 3</xref>). There were significant differences in bactericide activity at different concentrations for all species (<italic>P</italic>&lt;0.05), except for <italic>S. lentus</italic> (χ<sup>2</sup> = 4.69, d.f. = 5, <italic>P</italic> = 0.455).</p>
<fig id="pone-0111057-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0111057.g003</object-id><label>Figure 3</label><caption>
<title>Sensitivity (%; A) and inhibition halo (mm, mean ± SD; B) of bacteria species isolated from faeces after 24 h exposure to cantharidin concentration (control without cantharidin, 1∶100000, 1∶10000, 1∶1000, 1∶100 and 1∶10).</title>
</caption><graphic mimetype="image" xlink:href="info:doi/10.1371/journal.pone.0111057.g003" position="float" xlink:type="simple"/></fig></sec><sec id="s3c">
<title>Cantharidin in kidney and liver of great bustards</title>
<p>Cantharidin was found in kidneys and livers of 28% of 25 available great bustard carcasses. Cantharidin was found in 3 carcasses collected in spring (March-April), but also in 4 carcasses collected in winter (January-February), suggesting that in the latter individuals cantharidin was stored since the previous spring, as no blister beetles are found in winter. Presence and concentration of cantharidin did not differ between sexes (21% males, 36% females with cantharidin, χ<sup>2</sup> = 28.97, <italic>P</italic> = 0.41; male liver: 6.46±4.50 µg/g, female liver: 8.02±4.62 µg/g; χ<sup>2</sup> = 0.05, <italic>P</italic> = 0.83; male kidney: 0.76±0.23 µg/g, female kidney: 0.55±0.15 µg/g; χ<sup>2</sup> = 0.33, <italic>P</italic> = 0.56).</p>
</sec><sec id="s3d">
<title>Cantharidin concentration in blister beetles</title>
<p>Average cantharidin concentrations in blister beetles were 39.51±47.58 mg/g (range: 1.5–156.7 mg/g) and 19.32±10.96mg/g (range: 7.4–30.9 mg/g), respectively in <italic>Berberomeloe majalis</italic> (<italic>n</italic> = 58) and <italic>Physomeloe corallifer</italic> (<italic>n</italic> = 5).</p>
</sec><sec id="s3e">
<title>Selection of blister beetles by great bustards</title>
<p>The frequency distribution of blister beetles in faeces did not differ between years (χ<sup>2</sup> = 1.89, d.f. = 3, <italic>P</italic> = 0.595). At least one blister beetle was found in 34% of the samples (<italic>n</italic> = 72), with 62 faeces containing just one individual (29%), 6 faeces two individuals (3%) and 4 faeces three individuals (2%). Blister beetle consumption by great bustards was clearly male-biased. First, blister beetles were significantly more frequent in male than in female faeces (<xref ref-type="table" rid="pone-0111057-t002">Table 2</xref>). Second, the biomass and abundance of blister beetles were significantly greater in male than in female faeces. Third, males selected blister beetles among other invertebrates, both in number of individuals and biomass, compared to their abundance in the field, whereas females did not (<xref ref-type="fig" rid="pone-0111057-g004">Fig. 4 A, B</xref>). And finally, males selected the largest blister beetles among those available, whereas females did not (<xref ref-type="fig" rid="pone-0111057-g004">Fig. 4 C</xref>). Moreover, among seven insect taxa consumed by great bustards (Acrididae, Formicidae, Hemiptera, Curculionidae, Meloidae, Scarabeidae and Tenebrionidae), there was only a selection of larger prey sizes by males compared to females in Tenebrionidae (F<sub>1,40</sub> = 9.61 p = 0.003) and Meloidae (F<sub>1,22</sub> = 4.37 p = 0.048), and the largest size differences were found in the latter (Meloidae: 1.6 times larger prey in great bustard males; 595.3 mg average size in males vs 370.7 mg in females; Tenebrionidae: 1.3 times larger prey in males; 144.8 mg average size in males <italic>vs</italic> 109.5 mg in females).</p>
<fig id="pone-0111057-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0111057.g004</object-id><label>Figure 4</label><caption>
<title>Blister beetle relative abundance (A) and biomass (B) in relation to all invertebrates sampled in field transects (column ‘available’), and found in faeces of females and males (columns ‘consumed’) in two springs (2006 and 2007).</title>
<p>(C) Mean body size (± SD, g) of blister beetles sampled in field transects (column ‘available’) and found in faeces of females and males (columns ‘consumed’) in two springs (2006 and 2007). Differences in A and B were compared with the test of differences between percentages, and those in C with a Student's t-test, in all cases with data from both years pooled (horizontal lines; ns: <italic>P</italic>&gt;0.05, *: <italic>P</italic>&lt;0.05, **: <italic>P</italic>&lt;0.01).</p>
</caption><graphic mimetype="image" xlink:href="info:doi/10.1371/journal.pone.0111057.g004" position="float" xlink:type="simple"/></fig><table-wrap id="pone-0111057-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0111057.t002</object-id><label>Table 2</label><caption>
<title>Differences in blister beetle consumption between male and female great bustards.</title>
</caption><alternatives><graphic id="pone-0111057-t002-2" position="float" mimetype="image" xlink:href="info:doi/10.1371/journal.pone.0111057.t002" xlink:type="simple"/>
<table><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup>
<thead>
<tr>
<td align="left" rowspan="1" colspan="1">Consumption of blister beetles</td>
<td align="left" rowspan="1" colspan="1">♂♂</td>
<td align="left" rowspan="1" colspan="1">♀♀</td>
<td align="left" rowspan="1" colspan="1">test</td>
<td align="left" rowspan="1" colspan="1">df</td>
<td align="left" rowspan="1" colspan="1"><italic>P</italic></td>
</tr>
</thead>
<tbody>
<tr>
<td align="left" rowspan="1" colspan="1">Faeces (<italic>n</italic> total)</td>
<td align="left" rowspan="1" colspan="1">131</td>
<td align="left" rowspan="1" colspan="1">81</td>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Occurrence (% of faeces)</td>
<td align="left" rowspan="1" colspan="1">42.74</td>
<td align="left" rowspan="1" colspan="1">19.75</td>
<td align="left" rowspan="1" colspan="1">11.800<xref ref-type="table-fn" rid="nt103">†</xref></td>
<td align="left" rowspan="1" colspan="1">1</td>
<td align="left" rowspan="1" colspan="1">0.001</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Abundance (n/faeces)</td>
<td align="left" rowspan="1" colspan="1">0.52±0.69</td>
<td align="left" rowspan="1" colspan="1">0.22±0.50</td>
<td align="left" rowspan="1" colspan="1">3.201<xref ref-type="table-fn" rid="nt104">††</xref></td>
<td align="left" rowspan="1" colspan="1">1</td>
<td align="left" rowspan="1" colspan="1">0.001</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Weighted abundance (n/g of faeces)</td>
<td align="left" rowspan="1" colspan="1">0.35±0.72</td>
<td align="left" rowspan="1" colspan="1">0.19±0.47</td>
<td align="left" rowspan="1" colspan="1">2.041<xref ref-type="table-fn" rid="nt104">††</xref></td>
<td align="left" rowspan="1" colspan="1">1</td>
<td align="left" rowspan="1" colspan="1">0.041</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Biomass (g/faeces)</td>
<td align="left" rowspan="1" colspan="1">0.28±0.38</td>
<td align="left" rowspan="1" colspan="1">0.06±0.14</td>
<td align="left" rowspan="1" colspan="1">3.181<xref ref-type="table-fn" rid="nt104">††</xref></td>
<td align="left" rowspan="1" colspan="1">1</td>
<td align="left" rowspan="1" colspan="1">0.001</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Weighted biomass (g/g of faeces)</td>
<td align="left" rowspan="1" colspan="1">0.19±0.39</td>
<td align="left" rowspan="1" colspan="1">0.06±0.14</td>
<td align="left" rowspan="1" colspan="1">2.420<xref ref-type="table-fn" rid="nt104">††</xref></td>
<td align="left" rowspan="1" colspan="1">1</td>
<td align="left" rowspan="1" colspan="1">0.016</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Stomach content (<italic>n</italic> total)</td>
<td align="left" rowspan="1" colspan="1">12</td>
<td align="left" rowspan="1" colspan="1">5</td>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Occurrence (% of stomachs)</td>
<td align="left" rowspan="1" colspan="1">41.66</td>
<td align="left" rowspan="1" colspan="1">20</td>
<td align="left" rowspan="1" colspan="1">0.726<xref ref-type="table-fn" rid="nt103">†</xref></td>
<td align="left" rowspan="1" colspan="1">1</td>
<td align="left" rowspan="1" colspan="1">0.600</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Abundance (n/stomach content)</td>
<td align="left" rowspan="1" colspan="1">1.80±0.80</td>
<td align="left" rowspan="1" colspan="1">1±0</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"><xref ref-type="table-fn" rid="nt105">†††</xref></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Weighted abundance (n/g of stomach content)</td>
<td align="left" rowspan="1" colspan="1">0.15±0.14</td>
<td align="left" rowspan="1" colspan="1">0.06±0</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"><xref ref-type="table-fn" rid="nt105">†††</xref></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Biomass (g/stomach content)</td>
<td align="left" rowspan="1" colspan="1">1.07±0.51</td>
<td align="left" rowspan="1" colspan="1">0.37±0</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"><xref ref-type="table-fn" rid="nt105">†††</xref></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Weighted biomass (g/g of stomach content)</td>
<td align="left" rowspan="1" colspan="1">0.11±0.09</td>
<td align="left" rowspan="1" colspan="1">0.02±0</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"><xref ref-type="table-fn" rid="nt105">†††</xref></td>
</tr>
</tbody>
</table>
</alternatives><table-wrap-foot><fn id="nt103"><label>†</label><p>χ<sup>2</sup> test.</p></fn><fn id="nt104"><label>††</label><p>GLM test with Poisson distribution.</p></fn><fn id="nt105"><label>†††</label><p>Significance values cannot be determined because there was only one female stomach with a blister beetle.</p></fn><fn id="nt106"><label/><p>Occurrence (%) is the percentage of faeces or stomachs where blister beetles were found. Abundance is the number (n) of blister beetles per excrement and stomach. Weighted abundance (n/g) takes into account the larger size of male faeces and stomachs, and expresses the number of blister beetles per gram of faeces and stomach content (n/g). Biomass is expressed as dry weight (g) of blister beetles per excrement and stomach, and weighted biomass is defined as dry weight (g) of blister beetle per gram of excrement and stomach content (g/g) (<italic>n</italic> = 212 faeces and <italic>n</italic> = 17 stomach contents). Data are means ± SD.</p></fn></table-wrap-foot></table-wrap></sec></sec><sec id="s4">
<title>Discussion</title>
<p>Great bustards consumed blister beetles, which contain cantharidin, a highly toxic compound that is lethal in high doses <xref ref-type="bibr" rid="pone.0111057-SnchezBarbudo1">[25]</xref>, but has a demonstrated nematocide and bactericide activity even at low concentrations <xref ref-type="bibr" rid="pone.0111057-Campbell2">[53]</xref>–<xref ref-type="bibr" rid="pone.0111057-Lckmann1">[55]</xref>. There are various hypotheses that could explain this behaviour, the predictions of which are summarized in <xref ref-type="supplementary-material" rid="pone.0111057.s001">Table S1</xref> and discussed in detail below.</p>
<p>The <italic>nutritional value hypothesis</italic> would explain blister beetle consumption based on their large size, low mobility and short handling time (which make them a highly profitable prey). In spite of the clumped distribution of blister beetles, which would allow bustards to eat more than one at a time, the birds usually consumed only one individual (a maximum of three, <xref ref-type="table" rid="pone-0111057-t002">Table 2</xref>). This limited intake suggests that great bustards are not fully resistant to cantharidin, and indeed the intake of more than a few blister beetles is lethal <xref ref-type="bibr" rid="pone.0111057-SnchezBarbudo1">[25]</xref>. Blister beetles may be consumed for nutritional purposes to the limits of their toxicity, however, male-biased blister beetle consumption is not explained through this hypothesis.</p>
<p>The <italic>predator deterrence hypothesis</italic> suggests that cantharidin would make great bustards appear toxic to their predators, thus helping them to avoid predator attacks, as suggested for the spur-winged goose <italic>Plectropterus gambensis</italic>, an African bird that also eats blister beetles <xref ref-type="bibr" rid="pone.0111057-Bartram1">[19]</xref>. In support of this hypothesis, we found cantharidin in the kidneys and livers of great bustards. Against this hypothesis, intoxication of bustard predators like golden eagles <italic>Aquila chrysaetus</italic> or foxes <italic>Vulpes vulpes</italic> has never been reported. In several Spanish localities, local people report that great bustards were once hunted for their meat, and this species was sometimes rated as a delicacy <xref ref-type="bibr" rid="pone.0111057-Rotherham1">[56]</xref>, with no reported cases of intoxication in humans. Thus, we discard the predator deterrence hypothesis as an explanation of the presence of cantharidin in the kidneys or livers of great bustards.</p>
<p>The <italic>self-medication hypothesis</italic> proposes a health benefit to the consumer. Great bustards may use cantharidin to fight against disease agents. Following de Roode et al. <xref ref-type="bibr" rid="pone.0111057-deRoode1">[7]</xref>, five conditions must be met for food intake behaviour to qualify as medication:</p>
<list list-type="roman-lower"><list-item>
<p><bold>Ingestion or external application of a third species or chemical.</bold> Great bustards fulfil this condition because they consume blister beetles that contain cantharidin, a toxic chemical with high immunogenicity which is avoided by most animals <xref ref-type="bibr" rid="pone.0111057-Dettner1">[17]</xref>,<xref ref-type="bibr" rid="pone.0111057-Eisner1">[20]</xref>. Cantharidin could be excreted through the uropygial gland as a defence against ectoparasites and bacterial infection, but great bustards lack an uropygial gland <xref ref-type="bibr" rid="pone.0111057-King1">[57]</xref>. Cantharidin may be accumulated in kidneys and liver, and slowly disseminated for periods longer than few hours after ingestion.</p>
</list-item><list-item>
<p><bold>The chemical ingested should enhance host fitness by increasing tolerance to infection or reducing parasite load or parasite fitness.</bold> Although the conservation status of our study species did not allow us to perform manipulative experiments to show this (see below), it is known that cantharidin acts in blister beetles as fungicide <xref ref-type="bibr" rid="pone.0111057-Carrel1">[18]</xref> and nematocide <xref ref-type="bibr" rid="pone.0111057-Lckmann1">[55]</xref>, and that it has a demonstrated nematocide capacity <xref ref-type="bibr" rid="pone.0111057-Campbell2">[53]</xref>,<xref ref-type="bibr" rid="pone.0111057-Hoste1">[54]</xref>. We showed that cantharidin also demostrates bactericidal activity against potentially pathogenic bacteria affecting great bustards (<xref ref-type="fig" rid="pone-0111057-g003">Fig. 3</xref>).</p>
</list-item><list-item>
<p><bold>It should entail an antinutritional cost due to its toxicity.</bold> The negative effects of cantharidin for most animals are well known <xref ref-type="bibr" rid="pone.0111057-Bartram1">[19]</xref>,<xref ref-type="bibr" rid="pone.0111057-Eisner1">[20]</xref>,<xref ref-type="bibr" rid="pone.0111057-Meynier1">[58]</xref>,<xref ref-type="bibr" rid="pone.0111057-Vzien1">[59]</xref>. Great bustards tolerate cantharidin, but are not fully resistant to it, as suggested by the limited number of blister beetles found in faeces and stomachs (either none or just one blister beetle in most faeces; <xref ref-type="table" rid="pone-0111057-t002">Table 2</xref>; see also <xref ref-type="bibr" rid="pone.0111057-SnchezBarbudo1">[25]</xref>)</p>
</list-item><list-item>
<p><bold>The medication behaviour must be shown to operate in the natural environment of the host; the effect of an artificial diet does not demonstrate its relevance in nature.</bold> In our case, all samples of pathogens and tissues were taken from birds in their natural environment.</p>
</list-item><list-item>
<p><bold>The behaviour must be initiated by parasite infection (therapeutic medication) or in response to parasite risk (prophylactic medication).</bold> Great bustards show high prevalence of several parasites and pathogens <xref ref-type="bibr" rid="pone.0111057-GarcaMontijano1">[60]</xref> (<xref ref-type="table" rid="pone-0111057-t001">Table 1</xref>), which suggests blister beetle consumption could be both, therapeutic and prophylactic against most bacteria and parasites identified in faeces. As in ii), due to ethical and logistic reasons explained below it was not possible to test predictions in field conditions or through manipulative experiments.</p>
</list-item></list>
<p>Our results were consistent with most of Roode <italic>et al.</italic>'s conditions specified above. Thus, consumption of blister beetles by great bustards could be interpreted as a self-medication behaviour. Moreover, the higher consumption by males suggests an additional function only in this sex. As discussed below, the most plausible function we can think of is its use as a mechanism to appear more healthy to females and thus increase their mating success. Nevertheless, the suggested link between self-medication and sexual selection cannot be definitely tested through manipulative experiments using great bustards for various ethical and logistic reasons: (1) the species is globally endangered and strictly protected; (2) the risk of poisoning adult males in captivity, by allowing them to feed ad libitum on blister beetles; (3) capturing and recapturing wild individuals during courtship to take measurements and samples aiming at testing this prediction would imply increasing their stress level and exposing them to serious risk of death; (4) correlating mating success or pathogen prevalence to cantharidin consumption in individual males is practically impossible due to the low number of effective copulations <xref ref-type="bibr" rid="pone.0111057-Alonso1">[13]</xref>,<xref ref-type="bibr" rid="pone.0111057-Hellmich1">[30]</xref>,<xref ref-type="bibr" rid="pone.0111057-HidalgodeTrucios1">[31]</xref>,<xref ref-type="bibr" rid="pone.0111057-Gewalt1">[61]</xref>,<xref ref-type="bibr" rid="pone.0111057-Morales1">[62]</xref>, and to the difficulties mentioned in (3).</p>
<p>As for the higher consumption of blister beetles by males, it was clear from our results. First, the occurrence, abundance and biomass of blister beetles were much higher in the diet of male than female great bustards (<xref ref-type="table" rid="pone-0111057-t002">Table 2</xref>), even after controlling for the sexual difference in body size (<xref ref-type="table" rid="pone-0111057-t002">Table 2</xref>, weighted abundance, weighted biomass). Second, males but not females showed a clear preference for this prey among all invertebrates available (<xref ref-type="fig" rid="pone-0111057-g004">Fig. 4</xref>). And third, males selected those of largest size, showing significant large-size selection compared to females only in this (Meloidae) and another family (Tenebrionidae) of their prey spectrum, with much higher differences in average prey size in Meloidae. The absence of a sexual difference in cantharidin content in livers and kidneys of carcasses collected throughout the year did not necessarily contradict the male-biased consumption established through faeces analysis, because the presence of cantharidin in winter carcasses means that it was stored since the previous spring, and during this long time males and females could have made a differential use of the compound. For example, although males ingested more cantharidin than females, it is likely that males need more cantharidin to deal with a higher pathogen and parasite burden than females, since during courtship males surely have a reduced immune resistance associated to their strenuous investment in sexual display <xref ref-type="bibr" rid="pone.0111057-Alonso1">[13]</xref>,<xref ref-type="bibr" rid="pone.0111057-Folstad1">[26]</xref>,<xref ref-type="bibr" rid="pone.0111057-Roberts1">[27]</xref>.</p>
<p>The absence of a sex difference in prevalence of cantharidin-sensitive pathogens in faeces during the breeding season (with the exception of <italic>Kocuria</italic> sp, which has not been cited as a bacterium producing STDs <xref ref-type="bibr" rid="pone.0111057-Poiani1">[63]</xref>), supports this conclusion; otherwise the observed male-biased cantharidin intake would have resulted in a higher abundance of pathogens in females.</p>
<p>If males use cantharidin to reduce their intestinal parasite and pathogen load, a clean cloaca could be an honest signal of a low probability of being a carrier of STD agents. Infections with some of the identified pathogens in bustard faeces such as <italic>E. coli</italic> and <italic>Sthaphylococcus</italic> sp, or parasite helminths, usually produce diarrhoea symptoms <xref ref-type="bibr" rid="pone.0111057-Bailey2">[64]</xref>, which may be detected by females via visual inspection of the males' cloaca (<xref ref-type="fig" rid="pone-0111057-g002">Fig. 2E</xref>). Females probably accept mating with a male only after checking that helminths are not visible (<xref ref-type="fig" rid="pone-0111057-g002">Fig. 2D</xref>), and that diarrhoea and odour effects caused by these and other gastrointestinal parasites are absent around its cloaca <xref ref-type="bibr" rid="pone.0111057-Poiani1">[63]</xref>. Female pecking may induce the protrusion of the male's cloaca to better show the presence or absence of disease agents or their signs. We suggest that it is the importance of selecting the healthiest male that explains the exhaustive inspection of the male's cloaca performed by the female prior to mating, and contributes to the strongly skewed male mating success and extremely choosy mate selection shown by females in this species compared to other lekking birds <xref ref-type="bibr" rid="pone.0111057-Alonso1">[13]</xref>. Behaviours facilitating male health evaluation are particularly important in strongly polygynous species, in which a female is likely to mate with a male that has already copulated with other females, increasing the risk of disease transmission <xref ref-type="bibr" rid="pone.0111057-Benskin1">[65]</xref>. The ability of females to recognize and avoid infected males has been observed in other lekking birds <xref ref-type="bibr" rid="pone.0111057-Boyce1">[66]</xref>,<xref ref-type="bibr" rid="pone.0111057-Johnsgard1">[67]</xref>. It is not unlikely that self-medication has also evolved in great bustards as a sexual selection mechanism capable of transmitting to females a signal of good resistance to a poisonous compound, in a similar way as other costly secondary sexual traits are exhibited by males of many species, and correlated with their mating success.</p>
<p>Exposure of the cloaca is widespread in pre-copulatory displays of birds and has been related with high probability of transmission of sexual diseases <xref ref-type="bibr" rid="pone.0111057-Sheldon1">[29]</xref>. However, cloaca pecking has only been observed in three species <xref ref-type="bibr" rid="pone.0111057-Nakamura1">[28]</xref>,<xref ref-type="bibr" rid="pone.0111057-Davies1">[68]</xref>,<xref ref-type="bibr" rid="pone.0111057-Mller1">[69]</xref>. In these species, it is males who inspect the cloaca of females, contrasting with great bustard, and the behaviour has been explained within a context of high sperm competition, extreme degree of polygynandry, or a high level of extra-pair copulations <xref ref-type="bibr" rid="pone.0111057-Sheldon1">[29]</xref>,<xref ref-type="bibr" rid="pone.0111057-Birkhead1">[33]</xref>. Long-lived birds with mating systems characterized by multiple matings are the most likely candidates to have coevolved with virulent STDs <xref ref-type="bibr" rid="pone.0111057-Sheldon1">[29]</xref>,<xref ref-type="bibr" rid="pone.0111057-Lombardo1">[32]</xref>,<xref ref-type="bibr" rid="pone.0111057-Benskin1">[65]</xref>. This condition is fulfilled in the great bustard, and also in the great Indian bustard (<italic>Ardeotis nigriceps</italic>), a species that also consumes blister beetles (<italic>Cylindrothorax tenuicollis</italic>) and exhibits a similar cloaca inspection behaviour during courtship <xref ref-type="bibr" rid="pone.0111057-Rahmani1">[70]</xref>.</p>
<p>Another indirect argument supporting a relationship between self-medication and sexual selection is derived from the temporal scenario for the evolution of blister beetle consumption in great bustards. The intake of blister beetles in other species of the bustard family (Otididae) has been occasionally reported <xref ref-type="bibr" rid="pone.0111057-Rahmani1">[70]</xref>,<xref ref-type="bibr" rid="pone.0111057-Kok1">[71]</xref> in non-sister branches along their phylogeny <xref ref-type="bibr" rid="pone.0111057-Pitra1">[72]</xref>, suggesting that the potential to tolerate low levels of cantharidin could have evolved at early stages of the evolution of the family. Elaborate courting displays, including exposure of the cloaca, are a common feature for all known populations of great bustards <xref ref-type="bibr" rid="pone.0111057-Hellmich1">[30]</xref>,<xref ref-type="bibr" rid="pone.0111057-HidalgodeTrucios1">[31]</xref>,<xref ref-type="bibr" rid="pone.0111057-Gewalt1">[61]</xref>, and therefore, sexual selection through female choice was developed in <italic>O. tarda</italic> before subsequent geographic isolation of populations occurred <xref ref-type="bibr" rid="pone.0111057-Pitra2">[73]</xref>. On the other hand, consumption of blister beetles by great bustards has been only reported for the Iberian populations <xref ref-type="bibr" rid="pone.0111057-Lane1">[15]</xref>,<xref ref-type="bibr" rid="pone.0111057-Rocha1">[16]</xref>,<xref ref-type="bibr" rid="pone.0111057-Lucio1">[36]</xref>,<xref ref-type="bibr" rid="pone.0111057-Palacios1">[37]</xref>, suggesting that self-medicating behaviour might have evolved independently in these populations. An independent origin for self-medication in Iberian great bustards is fully congruent with the species' phylogeographic patterns and fossil record, which suggest genetic isolation of the Iberian populations with respect to other European nuclei for at least 200,000 years <xref ref-type="bibr" rid="pone.0111057-Pitra2">[73]</xref>–<xref ref-type="bibr" rid="pone.0111057-Snchez1">[75]</xref>. It is also congruent with the geographic distribution of blister beetles of the genera <italic>Berberomeloe</italic> and <italic>Physomeloe</italic>, which are present only in the Iberian Peninsula. Beetles of the family Meloidae co-occur with great bustards all over their geographic range <xref ref-type="bibr" rid="pone.0111057-Bologna1">[39]</xref>, but large species of the genus <italic>Berberomeloe</italic>, are actively and conspicuously wandering at full sunlight in the steppe areas where great bustards perform their courtship display <xref ref-type="bibr" rid="pone.0111057-GarcaParis1">[40]</xref>,<xref ref-type="bibr" rid="pone.0111057-PercinoDaniel1">[76]</xref>. Under this scenario, the rapid and recent evolution of blister beetle consumption, so far restricted at the population level, could be a direct consequence of the action of strong sexual selection at a local geographic scale.</p>
<p>In conclusion, the hypothesis that self-medication by male great bustards might have evolved in a sexual selection context still needs a definitive demonstration, but its proposal here is based on reasoned arguments and relevant literature (<xref ref-type="supplementary-material" rid="pone.0111057.s001">Table S1</xref>). We suggest that males use self-medication not only to overcome the pathogen burden during a strenuous display, but also to withstand the close scrutiny of choosy females, showing a clean cloaca as honest signal of health, and so increasing their mating success. The expected relationship between self-medication and subsequent greater mating success, necessary to demonstrate a direct link with sexual selection, is extremely difficult, if not impossible to measure in great bustards at present, but it could be explored in other species. In spite of the universality of sexual selection theory, self-medication has probably been overlooked as a sexually-selected mechanism enhancing male fitness.</p>
</sec><sec id="s5">
<title>Supporting Information</title>
<supplementary-material id="pone.0111057.s001" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document" xlink:href="info:doi/10.1371/journal.pone.0111057.s001" position="float" xlink:type="simple"><label>Table S1</label><caption>
<p>Hypotheses about the potential function of blister beetle consumption in great bustards. Hypotheses were grouped by selective processes: <italic>natural selection</italic> and <italic>sexual selection</italic>. The <italic>nutritional value hypothesis</italic> explains consumption of blister beetles based on their large size, low mobility and short handling time (i.e., high rate of net energy intake). The <italic>predator deterrence hypothesis</italic> suggests that great bustards may incorporate cantharidin in their tissues to become toxic to predators. The <italic>self-medication hypothesis</italic> poses a health benefit to the consumer. The <italic>sexual selection mechanism hypothesis</italic> suggests an indicator function based on the observed male-biased blister beetle consumption and male cloaca inspection by females. All hypotheses predict blister beetles are ingested, but differ in the amount ingested.</p>
<p>(DOCX)</p>
</caption></supplementary-material></sec></body>
<back>
<ack>
<p>We thank C. Ponce and A. Lucas for their cooperation in field samplings of invertebrates, and A. Lucas, N. Panizo, L. Hernández, A. Alcalde, M. de Simón, P. Mas-Peinado, I. Noguerón, A. Mora, J. Morente and B. Casares for their help in invertebrate identification in the laboratory. We also thank J. A. Díaz De Tuesta for identifying faecal microflora and carrying out antibiograms (Laboratorio Regional de Sanidad Animal de la Comunidad de Madrid), I. Luaces for identifying protozoa and helminths (GIR-Diagnostics), L. Cuadra and E. Fernández for carrying out chromatography analysis (Laboratorio de cromatografía MNCN-CSIC), and C. Palacín and F. Kovacs for their pictures of great bustard males. Necropsies were performed with dead great bustards submitted from rehabilitation centres or found in the field and frozen in the Museo Nacional de Ciencias Naturales-CSIC between 1991 and 2013 by J.C. Alonso and co-workers of his research group (J.A. Alonso, C. Palacín, C.A. Martín, C. Ponce, M. Magaña, S. Lane, M.B. Morales, A. García, R. Barrientos). We thank E. Crespo and P. Cervera (CERI Castilla-La Mancha), A. López (Dirección General de Áreas Protegidas y Biodiversidad, Castilla-La Mancha), J. A. Arranz (Dirección General Medio Natural, Castilla y León), M.J. Martos and I. Zorrilla (Centro de Análisis y Diagnóstico de la Fauna Silvestre, Andalucía), M.M. Elósegui, M. Lamuela (Sección de Hábitats, Gobierno de Navarra) and A. Onrubia for submitting bustard carcasses. We thank Daphne J. Fairbairn and three anonymous reviewers for their comments.</p>
<sec>
<title>Contributions</title>
<p>Conceived and designed the experiments: CB JCA LMB GB. Performed the experiments: CB. Analyzed the data: CB LMB. Contributed reagents/materials/analysis tools: CB LMB MGP GB JCA. Wrote the paper: CB JCA LMB MGP GB.</p>
</sec></ack>
<ref-list>
<title>References</title>
<ref id="pone.0111057-Ekblom1"><label>1</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Ekblom</surname><given-names>R</given-names></name>, <name name-style="western"><surname>Saether</surname><given-names>SA</given-names></name>, <name name-style="western"><surname>Hasselquist</surname><given-names>D</given-names></name>, <name name-style="western"><surname>Hannersjo</surname><given-names>D</given-names></name>, <name name-style="western"><surname>Fiske</surname><given-names>P</given-names></name>, <etal>et al</etal>. (<year>2005</year>) <article-title>Female choice and male humoral immune response in the lekking great snipe (Gallinago media)</article-title>. <source>Behavioral Ecology</source> <volume>16</volume>: <fpage>346</fpage>–<lpage>351</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Hamilton1"><label>2</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Hamilton</surname><given-names>WD</given-names></name>, <name name-style="western"><surname>Zuk</surname><given-names>M</given-names></name> (<year>1982</year>) <article-title>Heritable true fitness and bright birds: a role for parasites?</article-title> <source>Science</source> <volume>218</volume>: <fpage>384</fpage>–<lpage>387</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Hglund1"><label>3</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Höglund</surname><given-names>J</given-names></name>, <name name-style="western"><surname>Alatalo</surname><given-names>RV</given-names></name>, <name name-style="western"><surname>Lundberg</surname><given-names>A</given-names></name> (<year>1992</year>) <article-title>The effects of parasites on male ornaments and female choice in the lek-breeding black grouse (Tetrao tetrix)</article-title>. <source>Behavioral Ecology and Sociobiology</source> <volume>30</volume>: <fpage>71</fpage>–<lpage>76</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Kirkpatrick1"><label>4</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Kirkpatrick</surname><given-names>M</given-names></name>, <name name-style="western"><surname>Ryan</surname><given-names>M</given-names></name> (<year>1991</year>) <article-title>The evolution of mating preferences and the paradox of the lek</article-title>. <source>Nature</source> <volume>350</volume>: <fpage>33</fpage>–<lpage>38</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Clayton1"><label>5</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Clayton</surname><given-names>DH</given-names></name> (<year>1991</year>) <article-title>The influence of parasites on host sexual selection</article-title>. <source>Parasitology Today</source> <volume>7</volume>: <fpage>329</fpage>–<lpage>334</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Loye1"><label>6</label>
<mixed-citation publication-type="other" xlink:type="simple">Loye JE, Zuk M (1991) Bird-parasite interactions. Ecology, evolution, and behaviour. Oxford Ornithology Series: i-xv, 1–406.</mixed-citation>
</ref>
<ref id="pone.0111057-deRoode1"><label>7</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>de Roode</surname><given-names>JC</given-names></name>, <name name-style="western"><surname>Lefèvre</surname><given-names>T</given-names></name>, <name name-style="western"><surname>Hunter</surname><given-names>MD</given-names></name> (<year>2013</year>) <article-title>Self-medication in animals</article-title>. <source>Science</source> <volume>340</volume>: <fpage>150</fpage>–<lpage>151</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Lozano1"><label>8</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Lozano</surname><given-names>GA</given-names></name> (<year>1998</year>) <article-title>Parasitic stress and self-medication in wild animals</article-title>. <source>Stress and Behavior</source> <volume>27</volume>: <fpage>291</fpage>–<lpage>317</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Rodrguez1"><label>9</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Rodríguez</surname><given-names>E</given-names></name>, <name name-style="western"><surname>Wrangham</surname><given-names>R</given-names></name> (<year>1993</year>) <article-title>Zoopharmacognosy - the use of medicinal-plants by animals</article-title>. <source>Recent Advances Phytochemical</source> <volume>27</volume>: <fpage>89</fpage>–<lpage>105</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Hart1"><label>10</label>
<mixed-citation publication-type="other" xlink:type="simple">Hart BL (1997) Behavioral defence. In: Clayton DH, Moore J, editors. Host-parasite evolution: general principles and avian models. Oxford: Oxford University Press.</mixed-citation>
</ref>
<ref id="pone.0111057-Beaulieu1"><label>11</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Beaulieu</surname><given-names>M</given-names></name>, <name name-style="western"><surname>Schaefer</surname><given-names>HM</given-names></name> (<year>2013</year>) <article-title>Rethinking the role of dietary antioxidants through the lens of self-medication</article-title>. <source>Animal Behaviour</source> <volume>86</volume>: <fpage>17</fpage>–<lpage>24</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Forbey1"><label>12</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Forbey</surname><given-names>JS</given-names></name>, <name name-style="western"><surname>Harvey</surname><given-names>AL</given-names></name>, <name name-style="western"><surname>Huffman</surname><given-names>Ma</given-names></name>, <name name-style="western"><surname>Provenza</surname><given-names>FD</given-names></name>, <name name-style="western"><surname>Sullivan</surname><given-names>R</given-names></name>, <etal>et al</etal>. (<year>2009</year>) <article-title>Exploitation of secondary metabolites by animals: A response to homeostatic challenges</article-title>. <source>Integrative and comparative biology</source> <volume>49</volume>: <fpage>314</fpage>–<lpage>328</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Alonso1"><label>13</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Alonso</surname><given-names>JC</given-names></name>, <name name-style="western"><surname>Magaña</surname><given-names>M</given-names></name>, <name name-style="western"><surname>Palacín</surname><given-names>C</given-names></name>, <name name-style="western"><surname>Martín</surname><given-names>CA</given-names></name> (<year>2010</year>) <article-title>Correlates of male mating success in great bustard leks: the effects of age, weight, and display effort</article-title>. <source>Behavioral Ecology and Sociobiology</source> <volume>64</volume>: <fpage>1589</fpage>–<lpage>1600</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Alonso2"><label>14</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Alonso</surname><given-names>JC</given-names></name>, <name name-style="western"><surname>Magaña</surname><given-names>M</given-names></name>, <name name-style="western"><surname>Martín</surname><given-names>CA</given-names></name>, <name name-style="western"><surname>Palacín</surname><given-names>C</given-names></name> (<year>2010</year>) <article-title>Sexual traits as quality indicators in lekking male Great Bustards</article-title>. <source>Ethology</source> <volume>116</volume>: <fpage>1084</fpage>–<lpage>1098</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Lane1"><label>15</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Lane</surname><given-names>SJ</given-names></name>, <name name-style="western"><surname>Alonso</surname><given-names>JC</given-names></name>, <name name-style="western"><surname>Alonso</surname><given-names>JA</given-names></name>, <name name-style="western"><surname>Naveso</surname><given-names>MA</given-names></name> (<year>1999</year>) <article-title>Seasonal changes in diet and diet selection of great bustards (<italic>Otis tarda</italic>) in north-west Spain</article-title>. <source>Journal of Zoology</source> <volume>247</volume>: <fpage>201</fpage>–<lpage>214</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Rocha1"><label>16</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Rocha</surname><given-names>P</given-names></name>, <name name-style="western"><surname>Marques</surname><given-names>AT</given-names></name>, <name name-style="western"><surname>Moreira</surname><given-names>F</given-names></name> (<year>2005</year>) <article-title>Seasonal variation in Great Bustard <italic>Otis tarda</italic> diet in south Portugal with a focus on the animal component</article-title>. <source>Ardeola</source> <volume>52</volume>: <fpage>371</fpage>–<lpage>376</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Dettner1"><label>17</label>
<mixed-citation publication-type="other" xlink:type="simple">Dettner K (1997) Inter- and intraspecific transfer of toxic insect compound cantharidin. In: Dettner K, Bauer G, Völkl W, editors. Vertical Food Web Interactions. Heidelberg Berlin Springer pp. 115–145.</mixed-citation>
</ref>
<ref id="pone.0111057-Carrel1"><label>18</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Carrel</surname><given-names>JE</given-names></name>, <name name-style="western"><surname>Eisner</surname><given-names>T</given-names></name> (<year>1974</year>) <article-title>Cantharidin: potent feeding deterrent to insects</article-title>. <source>Science</source> <volume>183</volume>: <fpage>755</fpage>–<lpage>757</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Bartram1"><label>19</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Bartram</surname><given-names>S</given-names></name>, <name name-style="western"><surname>Boland</surname><given-names>W</given-names></name> (<year>2001</year>) <article-title>Chemistry and ecology of toxic birds</article-title>. <source>ChemBioChem</source> <volume>2</volume>: <fpage>809</fpage>–<lpage>811</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Eisner1"><label>20</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Eisner</surname><given-names>T</given-names></name>, <name name-style="western"><surname>Conner</surname><given-names>J</given-names></name>, <name name-style="western"><surname>Carrel</surname><given-names>JE</given-names></name>, <name name-style="western"><surname>McCormick</surname><given-names>JP</given-names></name>, <name name-style="western"><surname>Slagle</surname><given-names>AJ</given-names></name>, <etal>et al</etal>. (<year>1990</year>) <article-title>Systemic retention of ingested cantharidin by frogs</article-title>. <source>Chemoecology</source> <volume>1</volume>: <fpage>57</fpage>–<lpage>62</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Sandroni1"><label>21</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Sandroni</surname><given-names>P</given-names></name> (<year>2001</year>) <article-title>Aphrodisiacs past and present: a historical review</article-title>. <source>Clinical Autonomic Research</source> <volume>11</volume>: <fpage>303</fpage>–<lpage>307</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Moed1"><label>22</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Moed</surname><given-names>L</given-names></name>, <name name-style="western"><surname>Shwayder</surname><given-names>TA</given-names></name>, <name name-style="western"><surname>Chang</surname><given-names>MW</given-names></name> (<year>2001</year>) <article-title>A blistering defense of an ancient medicine</article-title>. <source>Archives of Dermatology</source> <volume>137</volume>: <fpage>1357</fpage>–<lpage>1360</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Nickolls1"><label>23</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Nickolls</surname><given-names>LC</given-names></name>, <name name-style="western"><surname>Teare</surname><given-names>D</given-names></name> (<year>1954</year>) <article-title>Poisoning by cantharidin</article-title>. <source>British Medical Journal</source> <volume>2</volume>: <fpage>1384</fpage>–<lpage>1386</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Campbell1"><label>24</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Campbell</surname><given-names>BE</given-names></name>, <name name-style="western"><surname>Hofmann</surname><given-names>A</given-names></name>, <name name-style="western"><surname>McCluskey</surname><given-names>A</given-names></name>, <name name-style="western"><surname>Gasser</surname><given-names>RB</given-names></name> (<year>2011</year>) <article-title>Serine/threonine phosphatases in socioeconomically important parasitic nematodes–Prospects as novel drug targets?</article-title> <source>Biotechnology Advances</source> <volume>29</volume>: <fpage>28</fpage>–<lpage>39</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-SnchezBarbudo1"><label>25</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Sánchez-Barbudo</surname><given-names>IS</given-names></name>, <name name-style="western"><surname>Camarero</surname><given-names>P</given-names></name>, <name name-style="western"><surname>García Montijano</surname><given-names>M</given-names></name>, <name name-style="western"><surname>Mateo</surname><given-names>R</given-names></name> (<year>2012</year>) <article-title>Possible cantharidin poisoning of a great bustard (<italic>Otis tarda</italic>)</article-title>. <source>Toxicon</source> <volume>59</volume>: <fpage>100</fpage>–<lpage>103</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Folstad1"><label>26</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Folstad</surname><given-names>I</given-names></name>, <name name-style="western"><surname>Karter</surname><given-names>AJ</given-names></name> (<year>1992</year>) <article-title>Parasites, bright males, and the immunocompetence handicap</article-title>. <source>American Naturalist</source> <volume>139</volume>: <fpage>603</fpage>–<lpage>622</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Roberts1"><label>27</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Roberts</surname><given-names>ML</given-names></name>, <name name-style="western"><surname>Buchanan</surname><given-names>KL</given-names></name>, <name name-style="western"><surname>Evans</surname><given-names>MR</given-names></name> (<year>2004</year>) <article-title>Testing the immunocompetence handicap hypothesis: a review of the evidence</article-title>. <source>Animal Behaviour</source> <volume>68</volume>: <fpage>227</fpage>–<lpage>239</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Nakamura1"><label>28</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Nakamura</surname><given-names>M</given-names></name> (<year>1990</year>) <article-title>Cloacal protuberance and copulatory-behavior of the alpine accentor (Prunella collaris)</article-title>. <source>Auk</source> <volume>107</volume>: <fpage>284</fpage>–<lpage>295</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Sheldon1"><label>29</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Sheldon</surname><given-names>BC</given-names></name> (<year>1993</year>) <article-title>Sexually transmitted disease in birds: occurrence and evolutionary significance</article-title>. <source>Philosophical Transactions - Royal Society Biological Sciences</source> <volume>339</volume>: <fpage>491</fpage>–<lpage>497</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Hellmich1"><label>30</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Hellmich</surname><given-names>J</given-names></name> (<year>1991</year>) <article-title>El display de cortejo de la avutarda (<italic>Otis tarda</italic> L.)</article-title>. <source>Alytes</source> <volume>2</volume>: <fpage>127</fpage>–<lpage>150</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-HidalgodeTrucios1"><label>31</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Hidalgo de Trucios</surname><given-names>SJ</given-names></name>, <name name-style="western"><surname>Carranza</surname><given-names>J</given-names></name> (<year>1991</year>) <article-title>Timing, structure and functions of the courtship display in male great bustard</article-title>. <source>Ornis Scandinavica</source> <volume>22</volume>: <fpage>360</fpage>–<lpage>366</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Lombardo1"><label>32</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Lombardo</surname><given-names>M</given-names></name> (<year>1998</year>) <article-title>On the evolution of sexually transmitted diseases in birds</article-title>. <source>Journal of Avian Biology</source> <volume>29</volume>: <fpage>314</fpage>–<lpage>321</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Birkhead1"><label>33</label>
<mixed-citation publication-type="other" xlink:type="simple">Birkhead T (2000) Promiscuity: an evolutionary history of sperm competition and sexual conflict: Harvard University Press.</mixed-citation>
</ref>
<ref id="pone.0111057-Alonso3"><label>34</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Alonso</surname><given-names>JC</given-names></name>, <name name-style="western"><surname>Martín</surname><given-names>CA</given-names></name>, <name name-style="western"><surname>Palacín</surname><given-names>C</given-names></name>, <name name-style="western"><surname>Magaña</surname><given-names>M</given-names></name>, <name name-style="western"><surname>Martín</surname><given-names>B</given-names></name> (<year>2003</year>) <article-title>Distribution, size and recent trends of the great bustard <italic>Otis tarda</italic> population in Madrid region, Spain</article-title>. <source>Ardeola</source> <volume>50</volume>: <fpage>21</fpage>–<lpage>29</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Martn1"><label>35</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Martín</surname><given-names>B</given-names></name>, <name name-style="western"><surname>Alonso</surname><given-names>JC</given-names></name>, <name name-style="western"><surname>Martín</surname><given-names>CA</given-names></name>, <name name-style="western"><surname>Palacín</surname><given-names>C</given-names></name>, <name name-style="western"><surname>Magaña</surname><given-names>M</given-names></name>, <etal>et al</etal>. (<year>2012</year>) <article-title>Influence of spatial heterogeneity and temporal variability in habitat selection: A case study on a great bustard metapopulation</article-title>. <source>Ecological Modelling</source> <volume>228</volume>: <fpage>39</fpage>–<lpage>48</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Lucio1"><label>36</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Lucio</surname><given-names>AJ</given-names></name> (<year>1985</year>) <article-title>Datos sobre la alimentación de la avutarda (<italic>Otis tarda</italic> L. 1758) en la Cuenca del Duero</article-title>. <source>Alytes</source> <volume>3</volume>: <fpage>69</fpage>–<lpage>86</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Palacios1"><label>37</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Palacios</surname><given-names>F</given-names></name>, <name name-style="western"><surname>Garzón</surname><given-names>J</given-names></name>, <name name-style="western"><surname>Castroviejo</surname><given-names>J</given-names></name> (<year>1975</year>) <article-title>La alimentación de la avutarda (<italic>Otis tarda</italic>) en España, especialmente en primavera</article-title>. <source>Ardeola</source> <volume>21</volume>: <fpage>347</fpage>–<lpage>406</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Alonso4"><label>38</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Alonso</surname><given-names>JC</given-names></name>, <name name-style="western"><surname>Palacín</surname><given-names>C</given-names></name> (<year>2010</year>) <article-title>The world status and population trends of the great bustard (<italic>Otis tarda</italic>): 2010 update</article-title>. <source>Chinese Birds</source> <volume>1</volume>: <fpage>141</fpage>–<lpage>147</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Bologna1"><label>39</label>
<mixed-citation publication-type="other" xlink:type="simple">Bologna MA (1991) Coleoptera Meloidae. In: Calderini, editor. Fauna de Italia. Bologna.</mixed-citation>
</ref>
<ref id="pone.0111057-GarcaParis1"><label>40</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>García-Paris</surname><given-names>M</given-names></name> (<year>1998</year>) <article-title>Revisión sistemática del género <italic>Berberomeloe</italic> Bologna, 1988 (Coleoptera, Meloidae) y diagnosis de un endemismo ibérico olvidado</article-title>. <source>Graellsia</source> <volume>54</volume>: <fpage>97</fpage>–<lpage>109</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-GarcaPars1"><label>41</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>García-París</surname><given-names>M</given-names></name>, <name name-style="western"><surname>Trotta-Moreu</surname><given-names>N</given-names></name>, <name name-style="western"><surname>Capote</surname><given-names>L</given-names></name> (<year>2006</year>) <article-title>Estado de conocimiento actual y problemas de conservación de los Meloidae (Coleoptera) de la Comunidad de Madrid</article-title>. <source>Graellsia</source> <volume>62</volume>: <fpage>333</fpage>–<lpage>370</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Lane2"><label>42</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Lane</surname><given-names>SJ</given-names></name>, <name name-style="western"><surname>Alonso</surname><given-names>JC</given-names></name>, <name name-style="western"><surname>Martín</surname><given-names>CA</given-names></name> (<year>2001</year>) <article-title>Habitat preferences of great bustard <italic>Otis tarda</italic> flocks in the arable steppes of central Spain: are potentially suitable areas unoccupied?</article-title> <source>Journal of Applied Ecology</source> <volume>38</volume>: <fpage>193</fpage>–<lpage>203</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Magaa1"><label>43</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Magaña</surname><given-names>M</given-names></name>, <name name-style="western"><surname>Alonso</surname><given-names>JC</given-names></name>, <name name-style="western"><surname>Martín</surname><given-names>CA</given-names></name>, <name name-style="western"><surname>Bautista</surname><given-names>LM</given-names></name>, <name name-style="western"><surname>Martín</surname><given-names>B</given-names></name> (<year>2010</year>) <article-title>Nest-site selection by Great Bustards <italic>Otis tarda</italic> suggests a trade-off between concealment and visibility</article-title>. <source>Ibis</source> <volume>152</volume>: <fpage>77</fpage>–<lpage>89</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Quinn1"><label>44</label>
<mixed-citation publication-type="book" xlink:type="simple">Quinn PJ, Markey B.K., Carter W.J.C., Donnelly F.C. and Leonhard F.C. (2002) Veterinary Microbiology and Microbial Diseases. Oxford: Blackwell Science Ltd.</mixed-citation>
</ref>
<ref id="pone.0111057-Institute1"><label>45</label>
<mixed-citation publication-type="other" xlink:type="simple">Institute CLS (2006) Performance standards for antimicrobial disk susceptibility tests; Approved standard. In: Institute CLS, editor. CLSI document M2-A9 <volume>26</volume>: :1. Wayne, PA.</mixed-citation>
</ref>
<ref id="pone.0111057-Ray1"><label>46</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Ray</surname><given-names>AC</given-names></name>, <name name-style="western"><surname>Post</surname><given-names>LO</given-names></name>, <name name-style="western"><surname>Hurst</surname><given-names>JM</given-names></name>, <name name-style="western"><surname>Edwards</surname><given-names>WC</given-names></name>, <name name-style="western"><surname>Reagor</surname><given-names>JC</given-names></name> (<year>1980</year>) <article-title>Evaluation of an analytical method for the diagnosis of cantharidin toxicosis due to ingestion of blister beetles (<italic>Epicauta lemniscata</italic>) by horses and sheep</article-title>. <source>American Journal of Veterinary Research</source> <volume>41</volume>: <fpage>932</fpage>–<lpage>933</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Bailey1"><label>47</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Bailey</surname><given-names>T</given-names></name>, <name name-style="western"><surname>John</surname><given-names>A</given-names></name>, <name name-style="western"><surname>Mensah-Brown</surname><given-names>EP</given-names></name>, <name name-style="western"><surname>Garner</surname><given-names>A</given-names></name>, <name name-style="western"><surname>Samour</surname><given-names>J</given-names></name>, <etal>et al</etal>. (<year>1998</year>) <article-title>Drug metabolizing enzyme systems in the houbara bustard (<italic>Chlamydotis undulata</italic>)</article-title>. <source>Comparative biochemistry and physiology C Comparative pharmacology and toxicology</source> <volume>120</volume>: <fpage>365</fpage>–<lpage>372</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Finkelstein1"><label>48</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Finkelstein</surname><given-names>MM</given-names></name>, <name name-style="western"><surname>Verma</surname><given-names>DK</given-names></name> (<year>2001</year>) <article-title>Exposure estimation in the presence of nondetectable values: Another look</article-title>. <source>Aihaj</source> <volume>62</volume>: <fpage>195</fpage>–<lpage>198</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Hdar1"><label>49</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Hódar</surname><given-names>JA</given-names></name> (<year>1996</year>) <article-title>The use of regression equations for estimation of arthropod biomass in ecological studies</article-title>. <source>Acta Oecologica-International Journal of Ecology</source> <volume>17</volume>: <fpage>421</fpage>–<lpage>433</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Bautista1"><label>50</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Bautista</surname><given-names>LM</given-names></name>, <name name-style="western"><surname>Silván</surname><given-names>G</given-names></name>, <name name-style="western"><surname>Cáceres</surname><given-names>S</given-names></name>, <name name-style="western"><surname>Martínez-Fernández</surname><given-names>L</given-names></name>, <name name-style="western"><surname>Bravo</surname><given-names>C</given-names></name>, <etal>et al</etal>. (<year>2013</year>) <article-title>Faecal sexual steroids in sex typing and endocrine status of great bustards</article-title>. <source>European Journal of Wildlife Research</source> <volume>59</volume>: <fpage>815</fpage>–<lpage>822</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Hdar2"><label>51</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Hódar</surname><given-names>JA</given-names></name> (<year>1997</year>) <article-title>The use of regression equations for the estimation of prey length and biomass in diet studies of insectivore vertebrates</article-title>. <source>Miscelània Zoològica</source> <volume>20</volume>: <fpage>1</fpage>–<lpage>10</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-R1"><label>52</label>
<mixed-citation publication-type="other" xlink:type="simple">R Development Team Core (2012) R: A Language and environment for statistical computing. In: R Foundation for Statistical Computing, editor. 2.15.1 ed. Vienna, Austria.</mixed-citation>
</ref>
<ref id="pone.0111057-Campbell2"><label>53</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Campbell</surname><given-names>BE</given-names></name>, <name name-style="western"><surname>Tarleton</surname><given-names>M</given-names></name>, <name name-style="western"><surname>Gordon</surname><given-names>CP</given-names></name>, <name name-style="western"><surname>Sakoff</surname><given-names>JA</given-names></name>, <name name-style="western"><surname>Gilbert</surname><given-names>J</given-names></name>, <etal>et al</etal>. (<year>2011</year>) <article-title>Norcantharidin analogues with nematocidal activity in <italic>Haemonchus contortus</italic></article-title>. <source>Bioorganic and Medicinal Chemistry Letters</source> <volume>21</volume>: <fpage>3277</fpage>–<lpage>3281</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Hoste1"><label>54</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Hoste</surname><given-names>H</given-names></name>, <name name-style="western"><surname>Sotiraki</surname><given-names>S</given-names></name>, <name name-style="western"><surname>Landau</surname><given-names>S</given-names></name>, <name name-style="western"><surname>Jackson</surname><given-names>F</given-names></name>, <name name-style="western"><surname>Beveridge</surname><given-names>I</given-names></name> (<year>2010</year>) <article-title>Goat-nematode interactions: think differently</article-title>. <source>Trends in Parasitology</source> <volume>26</volume>: <fpage>376</fpage>–<lpage>381</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Lckmann1"><label>55</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Lückmann</surname><given-names>J</given-names></name>, <name name-style="western"><surname>Poinar</surname><given-names>G</given-names></name> (<year>2003</year>) <article-title>First record of a Mermithidae (Nematoda) from the meloid beetle <italic>Meloe violaceus</italic> Marsham, 1802 (Coleoptera: Meloidae)</article-title>. <source>Parasitology Research</source> <volume>90</volume>: <fpage>82</fpage>–<lpage>83</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Rotherham1"><label>56</label>
<mixed-citation publication-type="other" xlink:type="simple">Rotherham ID (2008) In: Hall CM, Sharples L, editors. Food and Wine Festivals and Events Around the World: Development, Management and Markets: Routledge. pp. 47–61.</mixed-citation>
</ref>
<ref id="pone.0111057-King1"><label>57</label>
<mixed-citation publication-type="other" xlink:type="simple">King AS, McLelland J (1984) Birds, their structure and function. Philadelphia: Baillière Tindall.</mixed-citation>
</ref>
<ref id="pone.0111057-Meynier1"><label>58</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Meynier</surname><given-names>J</given-names></name> (<year>1893</year>) <article-title>Empoisonnement par la chair de grenovilles infestées par des insectes du genre Mylabris de la familie des méloides</article-title>. <source>Archiv de Medecine et de Pharmacie Militaires</source> <volume>22</volume>: <fpage>53</fpage>–<lpage>56</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Vzien1"><label>59</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Vézien</surname><given-names>M</given-names></name> (<year>1861</year>) <article-title>Note sur la cystide cantharidienne par l'ingestion de grenouilles qui sont nourries de coléoptères vésicants Recueil de Mémoires de Medecine de Chirurgie et de Pharmacie Militaires</article-title>. <volume>4</volume>: <fpage>457</fpage>–<lpage>460</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-GarcaMontijano1"><label>60</label>
<mixed-citation publication-type="other" xlink:type="simple">García-Montijano M, Tébar AM, Barreiro B, Rodríguez P, Alonso JC, et al. Postmortem findings in wild great bustards (Otis tarda) from Spain: a clinical approach. In: (EAZWV) EAoZ-aWV, editor; 2002; Heidelberg, Germany.</mixed-citation>
</ref>
<ref id="pone.0111057-Gewalt1"><label>61</label>
<mixed-citation publication-type="other" xlink:type="simple">Gewalt W (1959) Die Großtrappe (Otis tarda L.). A Ziemsen Verlag, Wittenberg-Lutherstadt.</mixed-citation>
</ref>
<ref id="pone.0111057-Morales1"><label>62</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Morales</surname><given-names>MB</given-names></name>, <name name-style="western"><surname>Alonso</surname><given-names>JC</given-names></name>, <name name-style="western"><surname>Martín</surname><given-names>C</given-names></name>, <name name-style="western"><surname>Martín</surname><given-names>E</given-names></name>, <name name-style="western"><surname>Alonso</surname><given-names>J</given-names></name> (<year>2003</year>) <article-title>Male sexual display and attractiveness in the great bustard Otis tarda: the role of body condition</article-title>. <source>Journal of Ethology</source> <volume>21</volume>: <fpage>51</fpage>–<lpage>56</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Poiani1"><label>63</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Poiani</surname><given-names>A</given-names></name> (<year>2010</year>) <article-title>Do cloacal pathogenic microbes behave as sexually transmitted parasites in birds</article-title>. <source>The Open Ornithology Journal</source> <volume>3</volume>: <fpage>72</fpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Bailey2"><label>64</label>
<mixed-citation publication-type="other" xlink:type="simple">Bailey TA (2008) Disease and medical management of Houbara Bustards and other Otididae. Dubai: Emirates Printing Press.</mixed-citation>
</ref>
<ref id="pone.0111057-Benskin1"><label>65</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Benskin</surname><given-names>CMH</given-names></name>, <name name-style="western"><surname>Wilson</surname><given-names>K</given-names></name>, <name name-style="western"><surname>Jones</surname><given-names>K</given-names></name>, <name name-style="western"><surname>Hartley</surname><given-names>I</given-names></name> (<year>2009</year>) <article-title>Bacterial pathogens in wild birds: a review of the frequency and effects of infection</article-title>. <source>Biological Reviews</source> <volume>84</volume>: <fpage>349</fpage>–<lpage>373</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Boyce1"><label>66</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Boyce</surname><given-names>MS</given-names></name> (<year>1990</year>) <article-title>The red queen visits sage grouse leks</article-title>. <source>American Zoologist</source> <volume>30</volume>: <fpage>263</fpage>–<lpage>270</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Johnsgard1"><label>67</label>
<mixed-citation publication-type="other" xlink:type="simple">Johnsgard PA (1994) Arena birds: sexual selection and behaviour. Washington: Smithsonian Institution Press.</mixed-citation>
</ref>
<ref id="pone.0111057-Davies1"><label>68</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Davies</surname><given-names>NB</given-names></name> (<year>1983</year>) <article-title>Polyandry, cloaca-pecking and sperm competition in dunnocks</article-title>. <source>Nature</source> <volume>302</volume>: <fpage>334</fpage>–<lpage>336</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Mller1"><label>69</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Møller</surname><given-names>A</given-names></name> (<year>1987</year>) <article-title>House sparrow, <italic>Passer domesticus</italic>, communal displays</article-title>. <source>Animal Behaviour</source> <volume>35</volume>: <fpage>203</fpage>–<lpage>210</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Rahmani1"><label>70</label>
<mixed-citation publication-type="other" xlink:type="simple">Rahmani AR (1989) The Great Indian Bustard. Bombay Natural History Society.</mixed-citation>
</ref>
<ref id="pone.0111057-Kok1"><label>71</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Kok</surname><given-names>OB</given-names></name>, <name name-style="western"><surname>Earle</surname><given-names>RA</given-names></name> (<year>1990</year>) <article-title>Diet of the black korhaamn <italic>Eupodotis afra</italic> in the orange free state and northwest cape</article-title>. <source>Ostrich</source> <volume>61</volume>: <fpage>107</fpage>–<lpage>110</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Pitra1"><label>72</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Pitra</surname><given-names>C</given-names></name>, <name name-style="western"><surname>Lieckfeldt</surname><given-names>D</given-names></name>, <name name-style="western"><surname>Frahnert</surname><given-names>S</given-names></name>, <name name-style="western"><surname>Fickel</surname><given-names>J</given-names></name> (<year>2002</year>) <article-title>Phylogenetic relationships and ancestral areas of the bustards (Gruiformes: Otididae), inferred from mitochondrial DNA and nuclear intron sequences</article-title>. <source>Molecular Phylogenetics and Evolution</source> <volume>23</volume>: <fpage>63</fpage>–<lpage>74</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Pitra2"><label>73</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Pitra</surname><given-names>C</given-names></name>, <name name-style="western"><surname>Lieckfeldt</surname><given-names>D</given-names></name>, <name name-style="western"><surname>Alonso</surname><given-names>JC</given-names></name> (<year>2000</year>) <article-title>Population subdivision in Europe's great bustard inferred from mitochondrial and nuclear DNA sequence variation</article-title>. <source>Molecular Ecology</source> <volume>9</volume>: <fpage>1165</fpage>–<lpage>1170</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Garca1"><label>74</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>García</surname><given-names>JT</given-names></name>, <name name-style="western"><surname>Manosa</surname><given-names>S</given-names></name>, <name name-style="western"><surname>Morales</surname><given-names>MB</given-names></name>, <name name-style="western"><surname>Ponjoan</surname><given-names>A</given-names></name>, <name name-style="western"><surname>García de la Morena</surname><given-names>EL</given-names></name>, <etal>et al</etal>. (<year>2011</year>) <article-title>Genetic consequences of interglacial isolation in a steppe bird</article-title>. <source>Molecular Phylogenetics and Evolution</source> <volume>61</volume>: <fpage>671</fpage>–<lpage>676</lpage>.</mixed-citation>
</ref>
<ref id="pone.0111057-Snchez1"><label>75</label>
<mixed-citation publication-type="other" xlink:type="simple">Sánchez AM. Birds of the stratigraphic unit TG-11 at la Galería (Sierra de Atapuerca, Burgos, Spain). In: León JdCy, editor; 1995; Valladolid, Spain. pp.137–146.</mixed-citation>
</ref>
<ref id="pone.0111057-PercinoDaniel1"><label>76</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Percino-Daniel</surname><given-names>N</given-names></name>, <name name-style="western"><surname>Buckley</surname><given-names>D</given-names></name>, <name name-style="western"><surname>Garcia-Paris</surname><given-names>M</given-names></name> (<year>2013</year>) <article-title>Pharmacological properties of blister beetles (Coleoptera: Meloidae) promoted their integration into the cultural heritage of native rural Spain as inferred by vernacular names diversity, traditions, and mitochondrial DNA</article-title>. <source>Journal of Ethnopharmacology</source> <volume>147</volume>: <fpage>570</fpage>–<lpage>583</lpage>.</mixed-citation>
</ref>
</ref-list></back>
</article>