Transcriptomic responses of bat cells to European bat lyssavirus 1 infection under conditions simulating euthermia and hibernation.

BACKGROUND: Coevolution between pathogens and their hosts decreases host morbidity and mortality. Bats host and can tolerate viruses which can be lethal to other vertebrate orders, including humans. Bat adaptations to infection include localized immune response, early pathogen sensing, high interferon expression without pathogen stimulation, and regulated inflammatory response. The immune reaction is costly, and bats suppress high-cost metabolism during torpor. In the temperate zone, bats hibernate in winter, utilizing a specific behavioural adaptation to survive detrimental environmental conditions and lack of energy resources. Hibernation torpor involves major physiological changes that pose an additional challenge to bat-pathogen coexistence. Here, we compared bat cellular reaction to viral challenge under conditions simulating hibernation, evaluating the changes between torpor and euthermia. RESULTS: We infected the olfactory nerve-derived cell culture of Myotis myotis with an endemic bat pathogen, European bat lyssavirus 1 (EBLV-1). After infection, the bat cells were cultivated at two different temperatures, 37 °C and 5 °C, to examine the cell response during conditions simulating euthermia and torpor, respectively. The mRNA isolated from the cells was sequenced and analysed for differential gene expression attributable to the temperature and/or infection treatment. In conditions simulating euthermia, infected bat cells produce an excess signalling by multitude of pathways involved in apoptosis and immune regulation influencing proliferation of regulatory cell types which can, in synergy with other produced cytokines, contribute to viral tolerance. We found no up- or down-regulated genes expressed in infected cells cultivated at conditions simulating torpor compared to non-infected cells cultivated under the same conditions. When studying the reaction of uninfected cells to the temperature treatment, bat cells show an increased production of heat shock proteins (HSPs) with chaperone activity, improving the bat's ability to repair molecular structures damaged due to the stress related to the temperature change. CONCLUSIONS: The lack of bat cell reaction to infection in conditions simulating hibernation may contribute to the virus tolerance or persistence in bats. Together with the cell damage repair mechanisms induced in response to hibernation, the immune regulation may promote bats' ability to act as reservoirs of zoonotic viruses such as lyssaviruses.

Bat population recoveries give insight into clustering strategies during hibernation.

BACKGROUND: Behaviour during hibernation contributes to energy conservation in winter. Hibernating bats select roosts with respect to physiological and environmental stressors, available local microclimate and species-specific requirements. RESULTS: We found that, in the period between 1977 and 2018, hibernating Myotis myotis and Rhinolophus hipposideros bats showed exponential population growth. The growth rates, corrected for local winter seasonal severity and winter duration, were equal to 10 and 13%, respectively. While R. hipposideros only utilised the thermally stable and, at survey time, warmer corridors in the hibernaculum, an increasing proportion of M. myotis roosted in the thermally stable corridors as their abundance increased. About 14% of all hibernating M. myotis displayed solitary roosting, irrespective of other covariates. Those bats that clustered together formed progressively larger clusters with increasing abundance, particularly in cold corridors. We found no statistically significant relationship for clustering behaviour or cluster size with winter severity or winter duration. CONCLUSIONS: Abundance of hibernating bats is increasing in Central Europe. As the number of M. myotis bats increases, thermally unstable corridors become saturated with large clusters and the animals begin to roost deeper underground.

Modelling invasive pathogen load from non-destructive sampling data.

Where microbes colonizing skin surface may help maintain organism homeostasis, those that invade living skin layers cause disease. In bats, white-nose syndrome is a fungal skin infection that affects animals during hibernation and may lead to mortality in severe cases. Here, we inferred the amount of fungus that had invaded skin tissue of diseased animals. We used simulations to estimate the unobserved disease severity in a non-lethal wing punch biopsy and to relate the simulated pathology to the measured fungal load in paired biopsies. We found that a single white-nose syndrome skin lesion packed with spores and hyphae of the causative agent, Pseudogymnoascus destructans, contains 48.93 pg of the pathogen DNA, which amounts to about 1560 P destructans genomes in one skin lesion. Relating the information to the known UV fluorescence in Nearctic and Palearctic bats shows that Nearctic bats carry about 1.7 µg of fungal DNA per cm2, whereas Palearctic bats have 0.04 µg cm-2 of P. destructans DNA. With the information on the fungal load that had invaded the host skin, the researchers can now calculate disease severity as a function of invasive fungal growth using non-destructive UV light transillumination of each bat's wing membranes. Our results will enable and promote thorough disease severity assessment in protected bat species without the need for extensive animal and laboratory labor sacrifices.

Numerous cold arousals and rare arousal cascades as a hibernation strategy in European Myotis bats.

Hibernating bats optimise the duration of torpor bouts and arousals in relation to hibernaculum microclimatic conditions and fat reserves. Clustering has significant physiological and ecological benefits, promoting successful hibernation of individuals. Such aggregations may help maintain optimal temperatures, allowing better energy utilisation than in solitarily bats. However, aroused bats in a cluster could conceivably disturb those still hibernating, starting an energy-demanding arousal process. Our study was conducted over two winters in two different hibernacula (cave and mine) in the Czech Republic, where Greater mouse-eared bats (Myotis myotis) have previously been diagnosed with white-nose syndrome. In 118 arousal episodes we recorded 193 individual arousals in which a warming phase was observed, 135 (69.9%) being cold arousals, where bats ceased increasing their body temperatures at ≤ 10 °C. The remaining arousals were standard normothermic arousals, where body (fur) surface temperatures reached > 20 °C. Cold arousals occurred during the mid- and late hibernation periods, suggesting they were a response to disturbance by a neighbour in the same cluster. Arousal cascades, where bats aroused in series, were rare (12.7%) and reached a maximum in mid-January. Our data suggest that Myotis bats prolong their torpor bouts using numerous cold arousals but few arousal cascades. Upon arrival of a bat, the clustered bats show tolerance to disturbing by conspecifics.

White-nose syndrome detected in bats over an extensive area of Russia.

BACKGROUND: Spatiotemporal distribution patterns are important infectious disease epidemiological characteristics that improve our understanding of wild animal population health. The skin infection caused by the fungus Pseudogymnoascus destructans emerged as a panzootic disease in bats of the northern hemisphere. However, the infection status of bats over an extensive geographic area of the Russian Federation has remained understudied. RESULTS: We examined bats at the geographic limits of bat hibernation in the Palearctic temperate zone and found bats with white-nose syndrome (WNS) on the European slopes of the Ural Mountains through the Western Siberian Plain, Central Siberia and on to the Far East. We identified the diagnostic symptoms of WNS based on histopathology in the Northern Ural region at 11° (about 1200 km) higher latitude than the current northern limit in the Nearctic. While body surface temperature differed between regions, bats at all study sites hibernated in very cold conditions averaging 3.6 °C. Each region also differed in P. destructans fungal load and the number of UV fluorescent skin lesions indicating skin damage intensity. Myotis bombinus, M. gracilis and Murina hilgendorfi were newly confirmed with histopathological symptoms of WNS. Prevalence of UV-documented WNS ranged between 16 and 76% in species of relevant sample size. CONCLUSIONS: To conclude, the bat pathogen P. destructans is widely present in Russian hibernacula but infection remains at low intensity, despite the high exposure rate.

Deeply torpid bats can change position without elevation of body temperature.

Because body temperature is tightly coupled to physiological function, hibernating animals entering deep torpor are typically immobile. We analysed thermal behaviour and locomotory activity of hibernating greater mouse-eared bats Myotis myotis and found two types of movement behaviour related to body temperature, i.e. movement at high fur temperature and at low fur temperatures (Tflow; <5°C). First Tflow movements appeared at the beginning of March and often occurred during long torpor bouts. In most cases, Tflow events represented slow displacements between clusters of bats. In several cases, however, departure or arrivals from and into clusters was also recorded without any elevation in body temperature. Distance travelled, flight duration and speed of locomotion during Tflow events was lower than in high fur temperature events. Such behaviour could allow bats to save energy long-term and prolong torpor bouts. Tflow movement in torpid bats significantly changes our understanding of basic hibernation principles and we strongly recommend further studies on the subject.

A Skull Might Lie: Modeling Ancestral Ranges and Diet from Genes and Shape of Tree Squirrels.

Tropical forests of Central and South America represent hotspots of biological diversity. Tree squirrels of the tribe Sciurini are an excellent model system for the study of tropical biodiversity as these squirrels disperse exceptional distances, and after colonizing the tropics of the Central and South America, they have diversified rapidly. Here, we compare signals from DNA sequences with morphological signals using pictures of skulls and computational simulations. Phylogenetic analyses reveal step-wise geographic divergence across the Northern Hemisphere. In Central and South America, tree squirrels form two separate clades, which split from a common ancestor. Simulations of ancestral distributions show western Amazonia as the epicenter of speciation in South America. This finding suggests that wet tropical forests on the foothills of Andes possibly served as refugia of squirrel diversification during Pleistocene climatic oscillations. Comparison of phylogeny and morphology reveals one major discrepancy: Microsciurus species are a single clade morphologically but are polyphyletic genetically. Modeling of morphology-diet relationships shows that the only group of species with a direct link between skull shape and diet are the bark-gleaning insectivorous species of Microsciurus. This finding suggests that the current designation of Microsciurus as a genus is based on convergent ecologically driven changes in morphology.

SigHunt: horizontal gene transfer finder optimized for eukaryotic genomes.

MOTIVATION: Genomic islands (GIs) are DNA fragments incorporated into a genome through horizontal gene transfer (also called lateral gene transfer), often with functions novel for a given organism. While methods for their detection are well researched in prokaryotes, the complexity of eukaryotic genomes makes direct utilization of these methods unreliable, and so labour-intensive phylogenetic searches are used instead. RESULTS: We present a surrogate method that investigates nucleotide base composition of the DNA sequence in a eukaryotic genome and identifies putative GIs. We calculate a genomic signature as a vector of tetranucleotide (4-mer) frequencies using a sliding window approach. Extending the neighbourhood of the sliding window, we establish a local kernel density estimate of the 4-mer frequency. We score the number of 4-mer frequencies in the sliding window that deviate from the credibility interval of their local genomic density using a newly developed discrete interval accumulative score (DIAS). To further improve the effectiveness of DIAS, we select informative 4-mers in a range of organisms using the tetranucleotide quality score developed herein. We show that the SigHunt method is computationally efficient and able to detect GIs in eukaryotic genomes that represent non-ameliorated integration. Thus, it is suited to scanning for change in organisms with different DNA composition. AVAILABILITY AND IMPLEMENTATION: Source code and scripts freely available for download at http://www.iba.muni.cz/index-en.php?pg=research-data-analysis-tools-sighunt are implemented in C and R and are platform-independent. CONTACT: 376090@mail.muni.cz or martinkova@ivb.cz.

Tribe-specific satellite DNA in non-domestic Bovidae.

Satellite sequences present in the centromeric and pericentric regions of chromosomes represent useful source of information. Changes in satellite DNA composition may coincide with the speciation and serve as valuable markers of phylogenetic relationships. Here, we examined satellite DNA clones isolated by laser microdissection of centromeric regions of 38 bovid species and categorized them into three types. Sat I sequences from members of Bovini/Tragelaphini/Boselaphini are similar to the well-documented 1.715 sat I DNA family. Sat I DNA from Caprini/Alcelaphini/Hippotragini/Reduncini/Aepycerotini/Cephalophini/Antilopini/Neotragini/Oreotragini form the second group homologous to the common 1.714 sat I DNA. The analysis of sat II DNAs isolated in our study confirmed conservativeness of these sequences within Bovidae. Newly described centromeric clones from Madoqua kirkii and Strepsiceros strepsiceros were similar in length and repetitive tandem arrangement but showed no similarity to any other satellite DNA in the GenBank database. Phylogenetic analysis of sat I sequences isolated in our study from 38 bovid species enabled the description of relationships at the subfamily and tribal levels. The maximum likelihood and Bayesian inference analyses showed a basal position of sequences from Oreotragini in the subfamily Antilopinae. According to the Bayesian inference analysis based on the indels in a partitioned mixed model, Antilopinae satellite DNA split into two groups with those from Neotragini as a basal tribe, followed by a stepwise, successive branching of Cephalophini, Aepycerotini and Antilopini sequences. In the second group, Reduncini sequences were basal followed by Caprini, Alcelaphini and Hippotragini.

Acyrthosiphon pisum AQP2: a multifunctional insect aquaglyceroporin.

Annotation of the recently sequenced genome of the pea aphid (Acyrthosiphon pisum) identified a gene ApAQP2 (ACYPI009194, Gene ID: 100168499) with homology to the Major Intrinsic Protein/aquaporin superfamily of membrane channel proteins. Phylogenetic analysis suggests that ApAQP2 is a member of an insect-specific clade of this superfamily. Homology model structures of ApAQP2 showed a novel array of amino acids comprising the substrate selectivity-determining "aromatic/arginine" region of the putative transport pore. Subsequent characterization of the transport properties of ApAQP2 upon expression in Xenopus oocytes supports an unusual substrate selectivity profile. Water permeability analyses show that the ApAQP2 protein exhibits a robust mercury-insensitive aquaporin activity. However unlike the water-specific ApAQP1 protein, ApAQP2 forms a multifunctional transport channel that shows a wide permeability profile to a range of linear polyols, including the potentially biologically relevant substrates glycerol, mannitol and sorbitol. Gene expression analysis indicates that ApAQP2 is highly expressed in the insect bacteriocytes (cells bearing the symbiotic bacteria Buchnera) and the fat body. Overall the results demonstrate that ApAQP2 is a novel insect aquaglyceroporin which may be involved in water and polyol transport in support of the Buchnera symbiosis and aphid osmoregulation.

Divergent evolutionary processes associated with colonization of offshore islands.

Oceanic islands have been a test ground for evolutionary theory, but here, we focus on the possibilities for evolutionary study created by offshore islands. These can be colonized through various means and by a wide range of species, including those with low dispersal capabilities. We use morphology, modern and ancient sequences of cytochrome b (cytb) and microsatellite genotypes to examine colonization history and evolutionary change associated with occupation of the Orkney archipelago by the common vole (Microtus arvalis), a species found in continental Europe but not in Britain. Among possible colonization scenarios, our results are most consistent with human introduction at least 5100 bp (confirmed by radiocarbon dating). We used approximate Bayesian computation of population history to infer the coast of Belgium as the possible source and estimated the evolutionary timescale using a Bayesian coalescent approach. We showed substantial morphological divergence of the island populations, including a size increase presumably driven by selection and reduced microsatellite variation likely reflecting founder events and genetic drift. More surprisingly, our results suggest that a recent and widespread cytb replacement event in the continental source area purged cytb variation there, whereas the ancestral diversity is largely retained in the colonized islands as a genetic 'ark'. The replacement event in the continental M. arvalis was probably triggered by anthropogenic causes (land-use change). Our studies illustrate that small offshore islands can act as field laboratories for studying various evolutionary processes over relatively short timescales, informing about the mainland source area as well as the island.

Complex patterns of host switching in New World arenaviruses.

We empirically tested the long-standing hypothesis of codivergence of New World arenaviruses (NWA) with their hosts. We constructed phylogenies for NWA and all known hosts and used them in reconciliation analyses. We also constructed a phylogenetic tree of all Sigmodontinae and Neotominae rodents and tested whether viral-host associations were phylogenetically clustered. We determined host geographical overlap to determine to what extent opportunity to switch hosts was limited by host relatedness or physical proximity. With the exception of viruses from North America, no phylogenetically codivergent pattern between NWA and their hosts was found. We found that different virus clades were clustered differently and that Clade B with members pathogenic to humans was randomly distributed across the rodent phylogeny. Furthermore, viral relatedness within Clade B was significantly explained by the geographic overlap of their hosts' ranges rather than host relatedness, indicating that they are capable of host switching opportunistically. This has important bearings on their potential to become panzootic. Together, these analyses suggest that NWA have not codiverged with their hosts and instead have evolved predominantly via host switching.

High predation risk decimates survival during the reproduction season.

Predators attack conspicuous prey phenotypes that are present in the environment. Male display behavior of conspicuous nuptial coloration becomes risky in the presence of a predator, and adult males face higher predation risk. High predation risk in one sex will lead to low survival and sex ratio bias in adult cohorts, unless the increased predation risk is compensated by higher escape rate.Here, we tested the hypothesis that sand lizards (Lacerta agilis) have sex-specific predation risk and escape rate. We expected the differences to manifest in changes in sex ratio with age, differences in frequency of tail autotomy, and in sex-specific survival rate.We developed a statistical model to estimate predation risk and escape rate, combining the observed sex ratio and frequency of tail autotomy with likelihood-based survival rate. Using Bayesian framework, we estimated the model parameters. We projected the date of the tail autotomy events from growth rates derived from capture-recapture data measurements.We found statistically stable sex ratio in age groups, equal frequency of tail regenerates between sexes, and similar survival rate. Predation risk is similar between sexes, and escape rate increases survival by about 5%. We found low survival rate and a low number of tail autotomy events in females during months when sand lizards mate and lay eggs, indicating high predator pressure throughout reproduction. Our data show that gravid females fail to escape predation.The risks of reproduction season in an ectotherm are a convolution of morphological changes (conspicuous coloration in males and body allometry changes in gravid females), behavior (nuptial displays), and environmental conditions which challenge lizard thermal performance. Performance of endotherm predators in cold spring months endangers gravid females more than displaying males in bright nuptial coloration.

Satellite DNA in Neotropical Deer Species.

The taxonomy and phylogenetics of Neotropical deer have been mostly based on morphological criteria and needs a critical revision on the basis of new molecular and cytogenetic markers. In this study, we used the variation in the sequence, copy number, and chromosome localization of satellite I-IV DNA to evaluate evolutionary relationships among eight Neotropical deer species. Using FISH with satI-IV probes derived from Mazama gouazoubira, we proved the presence of satellite DNA blocks in peri/centromeric regions of all analyzed deer. Satellite DNA was also detected in the interstitial chromosome regions of species of the genus Mazama with highly reduced chromosome numbers. In contrast to Blastocerus dichotomus, Ozotoceros bezoarticus, and Odocoileus virginianus, Mazama species showed high abundance of satIV DNA by FISH. The phylogenetic analysis of the satellite DNA showed close relationships between O. bezoarticus and B. dichotomus. Furthermore, the Neotropical and Nearctic populations of O. virginianus formed a single clade. However, the satellite DNA phylogeny did not allow resolving the relationships within the genus Mazama. The high abundance of the satellite DNA in centromeres probably contributes to the formation of chromosomal rearrangements, thus leading to a fast and ongoing speciation in this genus, which has not yet been reflected in the satellite DNA sequence diversification.

Sequence Analysis and FISH Mapping of Four Satellite DNA Families among Cervidae.

Centromeric and pericentromeric chromosome regions are occupied by satellite DNA. Satellite DNAs play essential roles in chromosome segregation, and, thanks to their extensive sequence variability, to some extent, they can also be used as phylogenetic markers. In this paper, we isolated and sequenced satellite DNA I-IV in 11 species of Cervidae. The obtained satellite DNA sequences and their chromosomal distribution were compared among the analysed representatives of cervid subfamilies Cervinae and Capreolinae. Only satI and satII sequences are probably present in all analysed species with high abundance. On the other hand, fluorescence in situ hybridisation (FISH) with satIII and satIV probes showed signals only in a part of the analysed species, indicating interspecies copy number variations. Several indices, including FISH patterns, the high guanine and cytosine (GC) content, and the presence of centromere protein B (CENP-B) binding motif, suggest that the satII DNA may represent the most important satellite DNA family that might be involved in the centromeric function in Cervidae. The absence or low intensity of satellite DNA FISH signals on biarmed chromosomes probably reflects the evolutionary reduction of heterochromatin following the formation of chromosome fusions. The phylogenetic trees constructed on the basis of the satellite I-IV DNA relationships generally support the present cervid taxonomy.

Comparative eco-physiology revealed extensive enzymatic curtailment, lipases production and strong conidial resilience of the bat pathogenic fungus Pseudogymnoascus destructans.

The genus Pseudogymnoascus encompasses soil psychrophilic fungi living also in caves. Some are opportunistic pathogens; nevertheless, they do not cause outbreaks. Pseudogymnoascus destructans is the causative agent of the white-nose syndrome, which is decimating cave-hibernating bats. We used comparative eco-physiology to contrast the enzymatic potential and conidial resilience of P. destructans with that of phylogenetically diverse cave fungi, including Pseudogymnoascus spp., dermatophytes and outdoor saprotrophs. Enzymatic potential was assessed by Biolog MicroArray and by growth on labelled substrates and conidial viability was detected by flow cytometry. Pseudogymnoascus destructans was specific by extensive losses of metabolic variability and by ability of lipid degradation. We suppose that lipases are important enzymes allowing fungal hyphae to digest and invade the skin. Pseudogymnoascus destructans prefers nitrogenous substrates occurring in bat skin and lipids. Additionally, P. destructans alkalizes growth medium, which points to another possible virulence mechanism. Temperature above 30 °C substantially decreases conidial viability of cave fungi including P. destructans. Nevertheless, survival of P. destructans conidia prolongs by the temperature regime simulating beginning of the flight season, what suggests that conidia could persist on the body surface of bats and contribute to disease spreading during bats active season.

A common partitivirus infection in United States and Czech Republic isolates of bat white-nose syndrome fungal pathogen Pseudogymnoascus destructans.

The psychrophilic (cold-loving) fungus Pseudogymnoascus destructans was discovered more than a decade ago to be the pathogen responsible for white-nose syndrome, an emerging disease of North American bats causing unprecedented population declines. The same species of fungus is found in Europe but without associated mortality in bats. We found P. destructans was infected with a mycovirus [named Pseudogymnoascus destructans partitivirus 1 (PdPV-1)]. The virus is bipartite, containing two double-stranded RNA (dsRNA) segments designated as dsRNA1 and dsRNA2. The cDNA sequences revealed that dsRNA1 dsRNA is 1,683 bp in length with an open reading frame (ORF) that encodes 539 amino acids (molecular mass of 62.7 kDa); dsRNA2 dsRNA is 1,524 bp in length with an ORF that encodes 434 amino acids (molecular mass of 46.9 kDa). The dsRNA1 ORF contains motifs representative of RNA-dependent RNA polymerase (RdRp), whereas the dsRNA2 ORF sequence showed homology with the putative capsid proteins (CPs) of mycoviruses. Phylogenetic analyses with PdPV-1 RdRp and CP sequences indicated that both segments constitute the genome of a novel virus in the family Partitiviridae. The purified virions were isometric with an estimated diameter of 33 nm. Reverse transcription PCR (RT-PCR) and sequencing revealed that all US isolates and a subset of Czech Republic isolates of P. destructans were infected with PdPV-1. However, PdPV-1 appears to be not widely dispersed in the fungal genus Pseudogymnoascus, as non-pathogenic fungi P. appendiculatus (1 isolate) and P. roseus (6 isolates) tested negative. P. destructans PdPV-1 could be a valuable tool to investigate fungal biogeography and the host-pathogen interactions in bat WNS.

Transcriptional host-pathogen responses of Pseudogymnoascus destructans and three species of bats with white-nose syndrome.

Understanding how context (e.g., host species, environmental conditions) drives disease susceptibility is an essential goal of disease ecology. We hypothesized that in bat white-nose syndrome (WNS), species-specific host-pathogen interactions may partly explain varying disease outcomes among host species. We characterized bat and pathogen transcriptomes in paired samples of lesion-positive and lesion-negative wing tissue from bats infected with Pseudogymnoascus destructans in three parallel experiments. The first two experiments analyzed samples collected from the susceptible Nearctic Myotis lucifugus and the less-susceptible Nearctic Eptesicus fuscus, following experimental infection and hibernation in captivity under controlled conditions. The third experiment applied the same analyses to paired samples from infected, free-ranging Myotis myotis, a less susceptible, Palearctic species, following natural infection and hibernation (n = 8 sample pairs/species). Gene expression by P. destructans was similar among the three host species despite varying environmental conditions among the three experiments and was similar within each host species between saprophytic contexts (superficial growth on wings) and pathogenic contexts (growth in lesions on the same wings). In contrast, we observed qualitative variation in host response: M. lucifugus and M. myotis exhibited systemic responses to infection, while E. fuscus up-regulated a remarkably localized response. Our results suggest potential phylogenetic determinants of response to WNS and can inform further studies of context-dependent host-pathogen interactions.

Cytochrome b phylogeography of chamois (Rupicapra spp.). Population contractions, expansions and hybridizations governed the diversification of the genus.

The chamois provides an excellent model for exploring the effect of historical and evolutionary events on diversification. We investigate cytochrome b (cytb) sequences in the 10 recognized subspecies of Rupicapra classified within 2 species: Rupicapra pyrenaica, with the subspecies parva, pyrenaica, and ornata, and Rupicapra rupicapra, with cartusiana, rupicapra, tatrica, carpatica, balcanica, asiatica, and caucasica. A fragment of 349 bp of the cytb was sequenced in 189 individuals. We identified 3 cytb lineages: Clade West in Iberia and Western Alps; Clade Central in the Apennines and the Massif of Chartreuse; and Clade East present in populations to the east of the Alps. The 2 proposed species were polyphyletic; the clades West and Central are represented in both, whereas the Clade East is restricted to R. rupicapra. In contrast to the current systematic, cytb phylogenies suggest the classification of the 10 subspecies of chamois into a single species, R. rupicapra. Phylogeny and geographical distribution of the 3 lineages show the effects of limited latitudinal range expansions, contractions, and hybridizations among highly divergent lineages, along with a major role of the glacial ice sheets of the Alps and the Pyrenees as barriers to gene flow, on the diversification of extant taxa.

The genetic impact of translocations and habitat fragmentation in chamois (Rupicapra) spp.

The chamois is a useful species with which to investigate the combined genetic impact of habitat fragmentation, over hunting, and translocations. Genetic variation within and between chamois (genus Rupicapra) populations was analyzed in 259 individuals from 16 sampling sites located in Italy, Spain, Slovakia, and the Czech Republic. Two mitochondrial DNA markers (control region and cytochrome b) and 11 nuclear microsatellites were typed. The principal results of this study can be summarized as follows: 1) high and significant differentiation between almost all chamois populations is observed even on a microgeographical scale, probably caused by the patchy distribution of this species, sharp geographical barriers to gene flow, and drift effects related to recent bottlenecks; 2) historical translocation events have left a clear genetic signature, including interspecific hybridization in some Alpine localities; 3) the Apennine subspecies of chamois, Rupicapra pyrenaica ornata, shows a high and similar level of divergence (about 1.5 My) from the Pyrenean (Rupicapra pyrenaica pyrenaica) and the Alpine (Rupicapra rupicapra) chamois; therefore, the specific status of these taxa should be revised. These results confirm the potential of population genetic analyses to dissect and interpret complex patterns of diversity in order to define factors important to conservation and management.