Changing landscapes of Switzerland: assessing effects of historic and future land uses and land covers on the occurrence of selected butterfly species

Summary Landscapes are continuously changing. Natural forces of change such as heavy rainfall and fires can exert lasting influences on their physical form. However, changes related to human activities have often shaped landscapes more distinctly. In Western Europe, especially modern agricultural practices and the expanse of overbuilt land have left their marks in the landscapes since the middle of the 20th century. In the recent years men realised that mare and more changes that were formerly attributed to natural forces might indirectly be the result of their own action. Perhaps the most striking landscape change indirectly driven by human activity we can witness in these days is the large withdrawal of Alpine glaciers. Together with the landscapes also habitats of animal and plant species have undergone vast and sometimes rapid changes that have been hold responsible for the ongoing loss of biodiversity. Thereby, still little knowledge is available about probable effects of the rate of landscape change on species persistence and disappearance. Therefore, the development and speed of land use/land cover in the Swiss communes between the 1950s and 1990s were reconstructed using 10 parameters from agriculture and housing censuses, and were further correlated with changes in butterfly species occurrences. Cluster analyses were used to detect spatial patterns of change on broad spatial scales. Thereby, clusters of communes showing similar changes or transformation rates were identified for single decades and put into a temporally dynamic sequence. The obtained picture on the changes showed a prevalent replacement of non-intensive agriculture by intensive practices, a strong spreading of urban communes around city centres, and transitions towards larger farm sizes in the mountainous areas. Increasing transformation rates toward more intensive agricultural managements were especially found until the 1970s, whereas afterwards the trends were commonly negative. However, transformation rates representing the development of residential buildings showed positive courses at any time. The analyses concerning the butterfly species showed that grassland species reacted sensitively to the density of livestock in the communes. This might indicate the augmented use of dry grasslands as cattle pastures that show altered plant species compositions. Furthermore, these species also decreased in communes where farms with an agricultural area >5ha have disappeared. The species of the wetland habitats were favoured in communes with smaller fractions of agricultural areas and lower densities of large farms (>10ha) but did not show any correlation to transformation rates. It was concluded from these analyses that transformation rates might influence species disappearance to a certain extent but that states of the environmental predictors might generally outweigh the importance of the corresponding rates. Information on the current distribution of species is evident for nature conservation. Planning authorities that define priority areas for species protection or examine and authorise construction projects need to know about the spatial distribution of species. Hence, models that simulate the potential spatial distribution of species have become important decision tools. The underlying statistical analyses such as the widely used generalised linear models (GLM) often rely on binary species presence-absence data. However, often only species presence data have been colleted, especially for vagrant, rare or cryptic species such as butterflies or reptiles. Modellers have thus introduced randomly selected absence data to design distribution models. Yet, selecting false absence data might bias the model results. Therefore, we investigated several strategies to select more reliable absence data to model the distribution of butterfly species based on historical distribution data. The results showed that better models were obtained when historical data from longer time periods were considered. Furthermore, model performance was additionally increased when long-term data of species that show similar habitat requirements as the modelled species were used. This successful methodological approach was further applied to assess consequences of future landscape changes on the occurrence of butterfly species inhabiting dry grasslands or wetlands. These habitat types have been subjected to strong deterioration in the recent decades, what makes their protection a future mission. Four spatially explicit scenarios that described (i) ongoing land use changes as observed between 1985 and 1997, (ii) liberalised agricultural markets, and (iii) slightly and (iv) strongly lowered agricultural production provided probable directions of landscape change. Current species-environment relationships were derived from a statistical model and used to predict future occurrence probabilities in six major biogeographical regions in Switzerland, comprising the Jura Mountains, the Plateau, the Northern and Southern Alps, as well as the Western and Eastern Central Alps. The main results were that dry grasslands species profited from lowered agricultural production, whereas overgrowth of open areas in the liberalisation scenario might impair species occurrence. The wetland species mostly responded with decreases in their occurrence probabilities in the scenarios, due to a loss of their preferred habitat. Further analyses about factors currently influencing species occurrences confirmed anthropogenic causes such as urbanisation, abandonment of open land, and agricultural intensification. Hence, landscape planning should pay more attention to these forces in areas currently inhabited by these butterfly species to enable sustainable species persistence. In this thesis historical data were intensively used to reconstruct past developments and to make them useful for current investigations. Yet, the availability of historical data and the analyses on broader spatial scales has often limited the explanatory power of the conducted analyses. Meaningful descriptors of former habitat characteristics and abundant species distribution data are generally sparse, especially for fine scale analyses. However, this situation can be ameliorated by broadening the extent of the study site and the used grain size, as was done in this thesis by considering the whole of Switzerland with its communes. Nevertheless, current monitoring projects and data recording techniques are promising data sources that might allow more detailed analyses about effects of long-term species reactions on landscape changes in the near future. This work, however, also showed the value of historical species distribution data as for example their potential to locate still unknown species occurrences. The results might therefore contribute to further research activities that investigate current and future species distributions considering the immense richness of historical distribution data. Résumé Les paysages changent continuellement. Des farces naturelles comme des pluies violentes ou des feux peuvent avoir une influence durable sur la forme du paysage. Cependant, les changements attribués aux activités humaines ont souvent modelé les paysages plus profondément. Depuis les années 1950 surtout, les pratiques agricoles modernes ou l'expansion des surfaces d'habitat et d'infrastructure ont caractérisé le développement du paysage en Europe de l'Ouest. Ces dernières années, l'homme a commencé à réaliser que beaucoup de changements «naturels » pourraient indirectement résulter de ses propres activités. Le changement de paysage le plus apparent dont nous sommes témoins de nos jours est probablement l'immense retraite des glaciers alpins. Avec les paysages, les habitats des animaux et des plantes ont aussi été exposés à des changements vastes et quelquefois rapides, tenus pour coresponsable de la continuelle diminution de la biodiversité. Cependant, nous savons peu des effets probables de la rapidité des changements du paysage sur la persistance et la disparition des espèces. Le développement et la rapidité du changement de l'utilisation et de la couverture du sol dans les communes suisses entre les années 50 et 90 ont donc été reconstruits au moyen de 10 variables issues des recensements agricoles et résidentiels et ont été corrélés avec des changements de présence des papillons diurnes. Des analyses de groupes (Cluster analyses) ont été utilisées pour détecter des arrangements spatiaux de changements à l'échelle de la Suisse. Des communes avec des changements ou rapidités comparables ont été délimitées pour des décennies séparées et ont été placées en séquence temporelle, en rendrent une certaine dynamique du changement. Les résultats ont montré un remplacement répandu d'une agriculture extensive des pratiques intensives, une forte expansion des faubourgs urbains autour des grandes cités et des transitions vers de plus grandes surfaces d'exploitation dans les Alpes. Dans le cas des exploitations agricoles, des taux de changement croissants ont été observés jusqu'aux années 70, alors que la tendance a généralement été inversée dans les années suivantes. Par contre, la vitesse de construction des nouvelles maisons a montré des courbes positives pendant les 50 années. Les analyses sur la réaction des papillons diurnes ont montré que les espèces des prairies sèches supportaient une grande densité de bétail. Il est possible que dans ces communes beaucoup des prairies sèches aient été fertilisées et utilisées comme pâturages, qui ont une autre composition floristique. De plus, les espèces ont diminué dans les communes caractérisées par une rapide perte des fermes avec une surface cultivable supérieure à 5 ha. Les espèces des marais ont été favorisées dans des communes avec peu de surface cultivable et peu de grandes fermes, mais n'ont pas réagi aux taux de changement. Il en a donc été conclu que la rapidité des changements pourrait expliquer les disparitions d'espèces dans certains cas, mais que les variables prédictives qui expriment des états pourraient être des descripteurs plus importants. Des informations sur la distribution récente des espèces sont importantes par rapport aux mesures pour la conservation de la nature. Pour des autorités occupées à définir des zones de protection prioritaires ou à autoriser des projets de construction, ces informations sont indispensables. Les modèles de distribution spatiale d'espèces sont donc devenus des moyens de décision importants. Les méthodes statistiques courantes comme les modèles linéaires généralisés (GLM) demandent des données de présence et d'absence des espèces. Cependant, souvent seules les données de présence sont disponibles, surtout pour les animaux migrants, rares ou cryptiques comme des papillons ou des reptiles. C'est pourquoi certains modélisateurs ont choisi des absences au hasard, avec le risque d'influencer le résultat en choisissant des fausses absences. Nous avons établi plusieurs stratégies, basées sur des données de distribution historique des papillons diurnes, pour sélectionner des absences plus fiables. Les résultats ont démontré que de meilleurs modèles pouvaient être obtenus lorsque les données proviennent des périodes de temps plus longues. En plus, la performance des modèles a pu être augmentée en considérant des données de distribution à long terme d'espèces qui occupent des habitats similaires à ceux de l'espèce cible. Vu le succès de cette stratégie, elle a été utilisée pour évaluer les effets potentiels des changements de paysage futurs sur la distribution des papillons des prairies sèches et marais, deux habitats qui ont souffert de graves détériorations. Quatre scénarios spatialement explicites, décrivant (i) l'extrapolation des changements de l'utilisation de sol tels qu'observés entre 1985 et 1997, (ii) la libéralisation des marchés agricoles, et une production agricole (iii) légèrement amoindrie et (iv) fortement diminuée, ont été utilisés pour générer des directions de changement probables. Les relations actuelles entre la distribution des espèces et l'environnement ont été déterminées par le biais des modèles statistiques et ont été utilisées pour calculer des probabilités de présence selon les scénarios dans six régions biogéographiques majeures de la Suisse, comportant le Jura, le Plateau, les Alpes du Nord, du Sud, centrales orientales et centrales occidentales. Les résultats principaux ont montré que les espèces des prairies sèches pourraient profiter d'une diminution de la production agricole, mais qu'elles pourraient aussi disparaître à cause de l'embroussaillement des terres ouvertes dû à la libéralisation des marchés agricoles. La probabilité de présence des espèces de marais a décrû à cause d'une perte générale des habitats favorables. De plus, les analyses ont confirmé que des causes humaines comme l'urbanisation, l'abandon des terres ouvertes et l'intensification de l'agriculture affectent actuellement ces espèces. Ainsi ces forces devraient être mieux prises en compte lors de planifications paysagères, pour que ces papillons diurnes puissent survivre dans leurs habitats actuels. Dans ce travail de thèse, des données historiques ont été intensivement utilisées pour reconstruire des développements anciens et pour les rendre utiles à des recherches contemporaines. Cependant, la disponibilité des données historiques et les analyses à grande échelle ont souvent limité le pouvoir explicatif des analyses. Des descripteurs pertinents pour caractériser les habitats anciens et des données suffisantes sur la distribution des espèces sont généralement rares, spécialement pour des analyses à des échelles fores. Cette situation peut être améliorée en augmentant l'étendue du site d'étude et la résolution, comme il a été fait dans cette thèse en considérant toute la Suisse avec ses communes. Cependant, les récents projets de surveillance et les techniques de collecte de données sont des sources prometteuses, qui pourraient permettre des analyses plus détaillés sur les réactions à long terme des espèces aux changements de paysage dans le futur. Ce travail a aussi montré la valeur des anciennes données de distribution, par exemple leur potentiel pour aider à localiser des' présences d'espèces encore inconnues. Les résultats peuvent contribuer à des activités de recherche à venir, qui étudieraient les distributions récentes ou futures d'espèces en considérant l'immense richesse des données de distribution historiques.

Effects of aerobic fitness on oxygen uptake kinetics in heavy intensity swimming.

This study aimed to characterise both the [Formula: see text] kinetics within constant heavy-intensity swimming exercise, and to assess the relationships between [Formula: see text] kinetics and other parameters of aerobic fitness, in well-trained swimmers. On separate days, 21 male swimmers completed: (1) an incremental swimming test to determine their maximal oxygen uptake [Formula: see text], first ventilatory threshold (VT), and the velocity associated with [Formula: see text] [Formula: see text] and (2) two square-wave transitions from rest to heavy-intensity exercise, to determine their [Formula: see text] kinetics. All the tests involved breath-by-breath analysis of freestyle swimming using a swimming snorkel. [Formula: see text] kinetics was modelled with two exponential functions. The mean values for the incremental test were 56.0 ± 6.0 ml min(-1) kg(-1), 1.45 ± 0.08 m s(-1); and 42.1 ± 5.7 ml min(-1) kg(-1) for [Formula: see text], [Formula: see text] and VT, respectively. For the square-wave transition, the time constant of the primary phase (τ(p)) averaged 17.3 ± 5.4 s and the relevant slow component (A'(sc)) averaged 4.8 ± 2.9 ml min(-1) kg(-1) [representing 8.9% of the end-exercise [Formula: see text] (%A'(sc))]. τ(p) was correlated with [Formula: see text] (r = -0.55, P = 0.01), but not with either [Formula: see text] (r = 0.05, ns) or VT (r = 0.14, ns). The %A'(sc) did not correlate with either [Formula: see text] (r = -0.14, ns) or [Formula: see text] (r = 0.06, ns), but was inversely related with VT (r = -0.61, P < 0.01). This study was the first to describe the [Formula: see text] kinetics in heavy-intensity swimming using specific swimming exercise and appropriate methods. As has been demonstrated in cycling, faster [Formula: see text] kinetics allow higher aerobic power outputs to be attained. The slow component seems to be reduced in swimmers with higher ventilatory thresholds.

Impact of quorum sensing on fitness of Pseudomonas aeruginosa.

In Pseudomonas aeruginosa, cell-cell communication based on N-acyl-homoserine lactone (AHL) signal molecules (termed quorum sensing) is known to control the production of extracellular virulence factors. Hence, in pathogenic interactions with host organisms, the quorum-sensing (QS) machinery can confer a selective advantage on P. aeruginosa. However, as shown by transcriptomic and proteomic studies, many intracellular metabolic functions are also regulated by quorum sensing. Some of these serve to regenerate the AHL precursors methionine and S-adenosyl-methionine and to degrade adenosine via inosine and hypoxanthine. The fact that a significant percentage of clinical and environmental isolates of P. aeruginosa is defective for QS because of mutation in the major QS regulatory gene lasR, raises the question of whether the QS machinery can have a negative impact on the organism's fitness. In vitro, lasR mutants have a higher probability to escape lytic death in stationary phase under alkaline conditions than has the QS-proficient wild type. Similar selective forces might also operate in natural environments.

Inclusive fitness theory and eusociality.

Arising from M. A. Nowak, C. E. Tarnita & E. O. Wilson 466, 1057-1062 (2010); Nowak et al. reply. Nowak et al. argue that inclusive fitness theory has been of little value in explaining the natural world, and that it has led to negligible progress in explaining the evolution of eusociality. However, we believe that their arguments are based upon a misunderstanding of evolutionary theory and a misrepresentation of the empirical literature. We will focus our comments on three general issues.

An approach to costs and energy consumption in private urban Spanish Mediterranean landscapes from a simplified model in sprinkle irrigation

The number of private gardens has increased in recent years, creating a more pleasant urban model, but not without having an environmental impact, including increased energy consumption, which is the focus of this study. The estimation of costs and energy consumption for the generic typology of private urban gardens is based on two simplifying assumptions: square geometry with surface areas from 25 to 500 m2 and hydraulic design with a single pipe. In total, eight sprinkler models have been considered, along with their possible working pressures, and 31 pumping units grouped into 5 series that adequately cover the range of required flow rates and pressures, resultin in 495 hydraulic designs repeated for two climatically different locations in the Spanish Mediterranean area (Girona and Elche). Mean total irrigation costs for the locality with lower water needs (Girona) and greater needs (Elche) were € 2,974 ha-¹ yr-¹ and € 3,383 ha-¹ yr-¹, respectively. Energy costs accounted for 11.4% of the total cost for the first location, and 23.0% for the second. While a suitable choice of the hydraulic elements of the setup is essential, as it may provide average energy savings of 77%, due to the low energy cost in relation to the cost of installation, the potential energy savings do not constitute a significant incentive for the irrigation system design. The low efficiency of the pumping units used in this type of garden is the biggest obstacle and constraint to achieving a high quality energy solution

The stationary distribution of a continuously varying strategy in a class-structured population under mutation-selection-drift balance.

Many traits and/or strategies expressed by organisms are quantitative phenotypes. Because populations are of finite size and genomes are subject to mutations, these continuously varying phenotypes are under the joint pressure of mutation, natural selection and random genetic drift. This article derives the stationary distribution for such a phenotype under a mutation-selection-drift balance in a class-structured population allowing for demographically varying class sizes and/or changing environmental conditions. The salient feature of the stationary distribution is that it can be entirely characterized in terms of the average size of the gene pool and Hamilton's inclusive fitness effect. The exploration of the phenotypic space varies exponentially with the cumulative inclusive fitness effect over state space, which determines an adaptive landscape. The peaks of the landscapes are those phenotypes that are candidate evolutionary stable strategies and can be determined by standard phenotypic selection gradient methods (e.g. evolutionary game theory, kin selection theory, adaptive dynamics). The curvature of the stationary distribution provides a measure of the stability by convergence of candidate evolutionary stable strategies, and it is evaluated explicitly for two biological scenarios: first, a coordination game, which illustrates that, for a multipeaked adaptive landscape, stochastically stable strategies can be singled out by letting the size of the gene pool grow large; second, a sex-allocation game for diploids and haplo-diploids, which suggests that the equilibrium sex ratio follows a Beta distribution with parameters depending on the features of the genetic system.

Consequences of genetic erosion on fitness and phenotypic plasticity in European tree frog populations (Hyla arborea).

The detrimental effects of genetic erosion on small isolated populations are widely recognized contrary to their interactions with environmental changes. The ability of genotypes to plastically respond to variability is probably essential for the persistence of these populations. Genetic erosion impact may be exacerbated if inbreeding affects plastic responses or if their maintenance were at higher phenotypic costs. To understand the interplay 'genetic erosion-fitness-phenotypic plasticity', we experimentally compared, in different environments, the larval performances and plastic responses to predation of European tree frogs (Hyla arborea) from isolated and connected populations. Tadpoles from isolated populations were less performant, but the traits affected were environmental dependant. Heterosis observed in crosses between isolated populations allowed attributing their low fitness to inbreeding. Phenotypic plasticity can be maintained in the face of genetic erosion as inducible defences in response to predator were identical in all populations. However, the higher survival and developmental costs for isolated populations in harsh conditions may lead to an additional fitness loss for isolated populations.

Skeletal muscle mitochondria in the elderly: effects of physical fitness and exercise training.

Context: Sarcopenia is thought to be associated with mitochondrial (M) loss. It is unclear whether the decrease in M content is consequent to aging per se or to decreased physical activity. Objectives: To examine the influence of fitness on M content and function, and to assess whether exercise could improve M function in older adults. Design and subjects: Three distinct studies were conducted: 1) a cross-sectional observation comparing M content and fitness in a large heterogeneous cohort of older adults; 2) a case-control study comparing chronically endurance-trained older adults (A) and sedentary (S) subjects matched for age and gender; 3) a 4-month exercise intervention in S. Setting: University-based clinical research center Outcomes: M volume density (Mv) was assessed by electron microscopy from vastus lateralis biopsies, electron transport chain proteins (ETC) by western blotting, mRNAs for transcription factors involved in M biogenesis by qRT-PCR and in-vivo oxidative capacity (ATPmax) by (31)P-MR spectroscopy. Peak oxygen uptake (VO2peak) was measured by GXT. Results: VO2peak was strongly correlated with Mv in eighty 60-80 yo adults. Comparison of A vs. S revealed differences in Mv, ATPmax and some ETC complexes. Finally, exercise intervention confirmed that S are able to recover Mv, ATPmax and specific transcription factors. Conclusions: These data suggest that 1) aging per se is not the primary culprit leading to M dysfunction, 2) an aerobic exercise program, even at an older age, can ameliorate the loss in skeletal muscle M content and may prevent aging muscle comorbidities and 3) the improvement of M function is all about content.

Increased body fat is independently and negatively related with cardiorespiratory fitness levels in children and adolescents with normal weight.

Body mass index (BMI) is related with cardiorespiratory fitness (CRF), but less is known regarding the combined relationships between BMI and body fat (BF) on CRF. Cross-sectional study included 2361 girls and 2328 boys aged 10–18 years living in the area of Lisbon, Portugal. BMI was calculated by measuring height and weight, and obesity was assessed by international criteria. BF was assessed by bioimpedance. CRF was assessed by the 20-m shuttle run and the participants were classified as normal-to-high or low-CRF level according to Fitness gram criterion-referenced standards. The prevalence of low CRF was 47 and 39% in girls and boys, respectively. The corresponding values for the prevalence of obesity were 4.8 and 5.6% (not significant) and of excess BF of 12.1 and 25.1% (P <0.001), respectively. In both sexes, BMI and BF were inversely related with CRF: r = – 0.53 and – 0.45 for BMI and % BF, respectively, in boys and the corresponding values in girls were – 0.50 and – 0.33 (all P <0.01). When compared with a participant with normal BMI and BF, the odds ratios (95% confidence interval) for low CRF were 1.94 (1.46–2.58) for a participant with normal BMI and high BF, and 6.19 (5.02–7.63) for a participant with high BMI and high BF. The prevalence of low-CRF levels is high in Portuguese youths. BF negatively influences CRF levels among children/adolescents with normal BMI.

Simulations in evolution. II. Relative fitness and the propagation of mutants.

In Neo-Darwinism, variation and natural selection are the two evolutionary mechanisms which propel biological evolution. Our previous article presented a histogram model [1] consisting in populations of individuals whose number changed under the influence of variation and/or fitness, the total population remaining constant. Individuals are classified into bins, and the content of each bin is calculated generation after generation by an Excel spreadsheet. Here, we apply the histogram model to a stable population with fitness F(1)=1.00 in which one or two fitter mutants emerge. In a first scenario, a single mutant emerged in the population whose fitness was greater than 1.00. The simulations ended when the original population was reduced to a single individual. The histogram model was validated by excellent agreement between its predictions and those of a classical continuous function (Eqn. 1) which predicts the number of generations needed for a favorable mutation to spread throughout a population. But in contrast to Eqn. 1, our histogram model is adaptable to more complex scenarios, as demonstrated here. In the second and third scenarios, the original population was present at time zero together with two mutants which differed from the original population by two higher and distinct fitness values. In the fourth scenario, the large original population was present at time zero together with one fitter mutant. After a number of generations, when the mutant offspring had multiplied, a second mutant was introduced whose fitness was even greater. The histogram model also allows Shannon entropy (SE) to be monitored continuously as the information content of the total population decreases or increases. The results of these simulations illustrate, in a graphically didactic manner, the influence of natural selection, operating through relative fitness, in the emergence and dominance of a fitter mutant.

Long maximal incremental tests accurately assess aerobic fitness in class II and III obese men.

This study aimed to compare two different maximal incremental tests with different time durations [a maximal incremental ramp test with a short time duration (8-12 min) (STest) and a maximal incremental test with a longer time duration (20-25 min) (LTest)] to investigate whether an LTest accurately assesses aerobic fitness in class II and III obese men. Twenty obese men (BMI≥35 kg.m-2) without secondary pathologies (mean±SE; 36.7±1.9 yr; 41.8±0.7 kg*m-2) completed an STest (warm-up: 40 W; increment: 20 W*min-1) and an LTest [warm-up: 20% of the peak power output (PPO) reached during the STest; increment: 10% PPO every 5 min until 70% PPO was reached or until the respiratory exchange ratio reached 1.0, followed by 15 W.min-1 until exhaustion] on a cycle-ergometer to assess the peak oxygen uptake [Formula: see text] and peak heart rate (HRpeak) of each test. There were no significant differences in [Formula: see text] (STest: 3.1±0.1 L*min-1; LTest: 3.0±0.1 L*min-1) and HRpeak (STest: 174±4 bpm; LTest: 173±4 bpm) between the two tests. Bland-Altman plot analyses showed good agreement and Pearson product-moment and intra-class correlation coefficients showed a strong correlation between [Formula: see text] (r=0.81 for both; p≤0.001) and HRpeak (r=0.95 for both; p≤0.001) during both tests. [Formula: see text] and HRpeak assessments were not compromised by test duration in class II and III obese men. Therefore, we suggest that the LTest is a feasible test that accurately assesses aerobic fitness and may allow for the exercise intensity prescription and individualization that will lead to improved therapeutic approaches in treating obesity and severe obesity.

In vivo survival of teicoplanin-resistant Staphylococcus aureus and fitness cost of teicoplanin resistance.

Glycopeptide resistance, in a set of in vitro step-selected teicoplanin-resistant mutants derived from susceptible Staphylococcus aureus SA113, was associated with slower growth, thickening of the bacterial cell wall, increased N-acetylglucosamine incorporation, and decreased hemolysis. Differential transcriptome analysis showed that as resistance increased, some virulence-associated genes became downregulated. In a mouse tissue cage infection model, an inoculum of 10(4) CFU of strain SA113 rapidly produced a high-bacterial-load infection, which triggered MIP-2 release, leukocyte infiltration, and reduced leukocyte viability. In contrast, with the same inoculum of the isogenic glycopeptide-resistant derivative NM67, CFU initially decreased, resulting in the elimination of the mutant in three out of seven cages. In the four cages in which NM67 survived, it partially regained wild-type characteristics, including thinning of the cell wall, reduced N-acetylglucosamine uptake, and increased hemolysis; however, the survivors also became teicoplanin hypersusceptible. The elimination of the teicoplanin-resistant mutants and selection of teicoplanin-hypersusceptible survivors in the tissue cages indicated that glycopeptide resistance imposes a fitness burden on S. aureus and is selected against in vivo, with restoration of fitness incurring the price of resistance loss.