A new perspective on the evolutionary history of western European Sorex araneus group revealed by paternal and maternal molecular markers

The species of the common shrew (Sorex araneus) group are morphologically very similar, but have undergone a spectacular chromosomal evolution. We investigate here the evolutionary history of the Sorex araneus group distributed in western Europe. In particular, we clarify the position of a difficult species, S. granarius, using sex-specific (mtDNA and Y-chromosome) markers. The karyotype of S. granarius is generally considered similar to the common ancestor of the restricted group considered here. The mtDNA data (1.4 kb) confirms the close relationship between S. granarius and S. araneus sensu stricto (hereafter S. araneus s.s.), but the Y-chromosome (3.4 kb) produces a quite different picture: S. granarius is closely related to another species, S. coronatus. Comparison of mtDNA and Y-chromosome phylogenies suggests that the genetic and chromosomal evolution in this group are disconnected processes. The evolutionary history of the south-western European populations of the S. araneus group can only be understood considering secondary contacts between taxa after their divergence, implying genetic exchanges by means of hybridization and/or introgression.     

Karyotype evolution of eulipotyphla (insectivora): the genome homology of seven sorex species revealed by comparative chromosome painting and banding data

The genus Sorex is one of the most successful genera of Eulipotyphla. Species of this genus are characterized by a striking chromosome variability including XY1Y2 sex chromosome systems and exceptional chromosomal polymorphisms within and between populations. To study chromosomal evolution of the genus in detail, we performed cross-species chromosome painting of 7 Sorex species with S. granarius and S. araneus whole-chromosome probes and found that the tundra shrew S. tundrensis has the most rearranged karyotype among these. We reconstructed robust phylogeny of the genus Sorex based on revealed conserved chromosomal segments and syntenic associations. About 16 rearrangements led to formation of 2 major Palearctic groups after their divergence from the common ancestor: the S. araneus group (10 fusions and 1 fission) and the S. minutus group (5 fusions). Further chromosomal evolution of the 12 species inside the groups, including 5 previously investigated species, was accompanied by multiple reshuffling events: 39 fusions, 20 centromere shifts and 10 fissions. The rate of chromosomal exchanges upon formation of the genus was close to the average rate for eutherians, but increased during recent (about 6-3 million years ago) speciation within Sorex. We propose that a plausible ancestral Sorex karyotype consists of 56 elements. It underwent 20 chromosome rearrangements from the boreoeutherian ancestor, with 14 chromosomes retaining the conserved state. The set of genus-specific chromosome signatures was drawn from the human (HSA)-shrew comparative map (HSA3/12/22, 8/19/3/21, 2/13, 3/18, 11/17, 12/15 and 1/12/22). The syntenic association HSA4/20, that was previously proposed as a common trait of all Eulipotyphla species, is shown here to be an apomorphic trait of S. araneus.     

3D organization of interphase fibroblast nuclei in two closely related shrew species, Sorex granarius and S. araneus, differed by the structure of chromosome termini

To study 3D organization of fibroblast interphase nuclei in two sibling shrew species, Sorex araneus from Cordon race and S. granarius, FISH with probe to telomeric and rDNA repeats, and immunofluorescence with ANA CREST and antibodies to nucleolus protein B23 were used. Karyotypes of studied species are composed of near identical chromosomal arms and differ by the number of metacentrics and the structure of terminal chromosome regions. The large telomeres containing on the average 218 kbp of telomere repeats characterize the short arms in all of 32 S. granarius acrocentrics. Telomere repeats in them alternate with nbosomal repeats. These regions also contain active NORs. In contrast, active NORs in S. araneus are localized at the terminal regions of 8 chromosomal arms (Zhdanova et al., 2005, 2007b). We have shown that telomere associations of chromosomes and contacts of a part of telomere clusters with inner nuclear membrane and nucleolus characterize interphase nuclei of both S. granarius and S. araneus. Moreover, the partial colocalization of telomere and ribosomal clusters, and spatial nearness of centomeric and telomeric regions were revealed in the interphase nuclei of S. granarius. Evidently, only those ribosomal clusters that contain a number of active ribosomal genes display connection with nucleolus. The stripping of nucleolus materials during transition of fibroblasts to mitosis and the role of B23 protein in this process has been studied.     

A hybrid zone with coincident clines for autosomal and sex-specific markers in the Sorex araneus group

In hybrid zones, endogenous counter-selection of hybrids is usually first expressed as reduced fertility or viability in hybrids of the heterogametic sex, a mechanism known as Haldane's rule. This phenomenon often leads to a differential of gene flow between sex-linked markers. Here, we address the possibility of a differential gene flow for Y chromosome, mtDNA and autosomal markers across the hybrid zone between the genetically and chromosomally well-differentiated species Sorex antinorii and Sorex araneus race Vaud. Intermarker comparison clearly revealed coincidental centre and very abrupt clines for all three types of markers. The overall level of genetic differentiation between the two species must be strong enough to hinder asymmetric introgression. Cyto-nuclear mismatches were also observed in the centre of hybrid zone. The significantly lower number of mismatches observed in males than in females possibly results from Y chromosome-mtDNA interactions. Results are compared with those previously reported in another hybrid zone between S. antinori and S. araneus race Cordon.     

Phylogenetic relationships of Caucasian shrew Sorex satunini Ogn. (mammalia) in the superspecies Sorex araneus inferred from the data of karyological analysis and the mtDNA cyt b gene sequencing

Using the data of karyological analysis, the phylogenetic relationships of Caucasian shrew Sorex satunini and the cryptic species of superspecies Sorer araneus were examined. In the population of Sorex satunini from the plain of North Ciscaucasia two deeply radiated cytochrome b genes (A and B) were identified. Genetic distance between haplotype A and B groups constituted 0.0675 +/- 0.008, which is higher than any distance in superspecies S. araneus. Possible introgression of type B haplotypes from the populations of the evolutionary lineage S. subaraneus--S. araneus in Pleistocene and the time of the appearance of the chromosomal polymorphism of S. araneus is discussed. Our results show that the use of only one mitochondrial marker can lead to false conclusions on taxonomic diversity     

Living at the physiological limits: field and maximum metabolic rates of the common shrew (Sorex araneus)

Shrews (genus Sorex, small insectivorous mammals) are well known for their extremely high basal metabolic rates (BMRs) even when corrected for their small body size. We measured energy expenditure of the common shrew (Sorex araneus) under natural conditions (field metabolic rate [FMR]) by doubly labeled water method to test whether FMR is proportional to high BMR in this species. The study was performed in summer in northeastern Poland. In addition to the FMR, we also measured maximum metabolic rates induced by cold exposure and by intense activity (MMRCOLD and MMRRUN, respectively) to evaluate the aerobic reserve (MMR-FMR) in S. araneus. This aerobic reserve was used as an indicator of the potential for metabolic constraints. The FMR averaged 2.31+/-0.32 L CO2 d(-1) (+/-SD) or 58.1+/-8.0 kJ d(-1) in 8.2-g animals. This figure constituted 216%-258% of a value predicted for a "standard" mammal of the same body mass and was the highest mass-specific field metabolic rate in mammals. Because of the high BMR level in S. araneus, the FMR to BMR ratio (2.4) was not far off mammalian standards (median value of 3.1). The rate of water efflux determined in S. araneus (20.2 mL H2O d(-1) or 2.46 mL H2O g(-1) d(-1)) exceeded all figures reported to date in other mammals and was apparently linked to the high FMR level and relatively high water content of shrews' food. Maximal metabolic rates (MMRRUN of 18.1+/-1.6 mL O2 g(-1) h(-1) and MMRCOLD of 23.5+/-1.9 mL O2 g(-1) h(-1)) were not high in proportion to BMR or FMR that resulted in relatively narrow aerobic reserve in S. araneus: 20% when calculated against the MMRRUN and 39% when compared with the MMRCOLD. Our study reveals that S. araneus has high energy costs of living and operates close to its physiological limits.     

The occurrence of Centrorhynchus (Acanthocephala) in shrews (Sorex araneus and Sorex minutus) in the United Kingdom

Encysted acanthocephalans belonging to the genus Centrorhynchus were found in the body cavities of Sorex araneus (common shrew) and Sorex minutus (pygmy shrew) from Boxworth, Cambridgeshire, U.K. Fifty percent of the male S. araneus and 67% of the male S. minutus examined were found to be infected, with the mean intensity (+/-SD) being 54.3 +/- 91.3 and 14.7 +/- 18.4, respectively. The species of Centrorhynchus in the shrews may be Centrorhynchus aluconis, which is distributed widely in tawny owls, Strix aluco, in the United Kingdom. Shrews appear to serve as paratenic hosts for C. aluconis.     

Genetic and karyotypic structure in the shrews of the Sorex araneus group: are they independent?

The species of the common shrew (Sorex araneus) group are morphologically very similar but exhibit high levels of karyotypic variation. Here we used genetic variation at 10 microsatellite markers in a data set of 212 individuals mostly sampled in the western Alps and composed of five karyotypic taxa (Sorex coronatus, Sorex antinorii and the S. araneus chromosome races Cordon, Bretolet and Vaud) to investigate the concordance between genetic and karyotypic structure. Bayesian analysis confirmed the taxonomic status of the three sampled species since individuals consistently grouped according to their taxonomical status. However, introgression can still be detected between S. antinorii and the race Cordon of S. araneus. This observation is consistent with the expected low karyotypic complexity of hybrids between these two taxa. Geographically based cryptic substructure was discovered within S. antinorii, a pattern consistent with the different postglaciation recolonization routes of this species. Additionally, we detected two genetic groups within S. araneus notwithstanding the presence of three chromosome races. This pattern can be explained by the probable hybrid status of the Bretolet race but also suggests a relatively low impact of chromosomal differences on genetic structure compared to historical factors. Finally, we propose that the current data set (available at http://www.unil.ch/dee/page7010_en.html#1) could be used as a reference by those wanting to identify Sorex individuals sampled in the western Alps.     

Telomeric DNA allocation in chromosomes of common shrew (Sorex araneus, eulipotyphla)

Recently, we displayed an Iberian shrew species (Sorex granarius) with telomere structures unusual for mammals. Long telomeres on the short acrocentric arms contain an average of 213 kb of telomere repeats, whereas the other chromosomal ends have only 3.8 kb (Zhdanova et al., 2005; 2007). However, it is not clear whether these telomeres are typical for all shrew species or only for S. granarius. S. granarius and common shrew Sorex araneus are sibling species. In this study, using modified Q-FISH we demonstrated that telomeres in S. araneus from various chromosomal races distinguished by their number of metacentrics contain 6.8–15.2 kb of telomeric tracts. The S. araneus telomere lengths appear to correspond to telomere lengths in the majority of both shrew species and wild mammals, whereas S. granarius has telomeres with unique or rare structures. Using DNA and RNA high-specific modified probes to telomeric repeats (PNA and LNA), we showed that interstitial telomeric sites in S. araneus chromosomes contain mainly telomeric DNA and that their localization coincide with some evolutionary breakpoints. Interstitial telomeric DNA in S. granarius chromosomes was not revealed. Thus, the distribution of telomeric DNA may be significantly different, even in closely related species whose chromosomes are composed of almost identical chromosomal arms.     

Interspecific competition in the shrews Sorex araneus and Sorex minutus (Soricidae, Insectivora): a population study of the Irish pygmy shrew

Sorex araneus and S. minutus , which are likely to show a temporary vertical habitat segregation, are largely sympatric, except in Ireland, where S. minutus occurs alone. Two mechanisms for this segregation may be suggested: a direct interspecific interaction or an evolutionary fixed relationship.
Population densities and surface activity were measured, both in areas where the two species occur together (The Netherlands) and in areas in Ireland, in (order to investigate the possibility of a habitat shift and its influence on population density. No such effect was found in the present study.     

Chromosomal rearrangements do not seem to affect the gene flow in hybrid zones between karyotypic races of the common shrew (Sorex araneus)

Chromosomal rearrangements are proposed to promote genetic differentiation between chromosomally differentiated taxa and therefore promote speciation. Due to their remarkable karyotypic polymorphism, the shrews of the Sorex araneus group were used to investigate the impact of chromosomal rearrangements on gene flow. Five intraspecific chromosomal hybrid zones characterized by different levels of karyotypic complexity were studied using 16 microsatellites markers. We observed low levels of genetic differentiation even in the hybrid zones with the highest karyotypic complexity. No evidence of restricted gene flow between differently rearranged chromosomes was observed. Contrary to what was observed at the interspecific level, the effect of chromosomal rearrangements on gene flow was undetectable within the S. araneus species.     

Phylogenetic relationships among Nearctic shrews of the genus Sorex (Insectivora, Soricidae) inferred from combined cytochrome b and inter-SINE fingerprint data using Bayesian analysis

The field of molecular systematics has relied heavily on mitochondrial DNA (mtDNA) analysis since its inception. Despite the obvious utility of mtDNA, such data inevitably only presents a limited (i.e., single genome) perspective on species evolution. A combination of mitochondrial and nuclear markers is essential for reconstructing more robust phylogenetic trees. To evaluate the utility of one category of nuclear marker (short interspersed elements or SINEs) for resolving phylogenetic relationships, we constructed an inter-SINE fingerprint for nine putative species of the genus Sorex. In addition, we analyzed 1011 nucleotides of the cytochrome b gene. Traditional neighbor-joining and maximum parsimony analyses were applied to the individual cytochrome b and inter-SINE fingerprint data sets, along with Bayesian analysis to the combined data sets. We found inter-SINE fingerprinting to be an effective species level marker; however, we were unable to reconstruct deeper branching patterns within the Sorex genus using these data. The combined data analyzed under a Bayesian analysis showed higher levels of structuring within the Otisorex subgenus, most notably recognizing a monophyletic group consisting of sister-taxa S. palustris and S. monticolus, S. cinereus and S. haydeni, and S. hoyi. An additional noteworthy result was the detection of an historic mitochondrial introgression event between S. monticolus and S. palustris. When combining disparate data sets, we emphasize researcher diligence as certain types of data and processes may overly influence the analysis. However, there is considerable phylogenetic potential stemming from inter-SINE fingerprinting.     

Pneumocystis carinii of the Common Shrew, Sorex araneus, Shows a Discrete Phenotype

ABSTRACT. We carried out an immunohistological and morphological study on Pneumocystis carinii originating from the common shrew, Sorex araneus . Immunologic properties were studied by applying two commercially available immunofluorescence staining kits with differing developmental form specificity to a lung homogenate. The cyst form-specific staining kit reacted with cysts originating from S. araneus . Ultrastructurally this particular antigen epitope specifically deposited on the electron-lucent middle layer of the cyst pellicle. The immunohistochemical staining kit reacting with both cyst and trophozoite forms from human and rat origin did not react with any developmental forms of P. carinii originating from S. araneus . Both kits demonstrated P. carinii in the lung homogenate of a field vole, Microtus agrestis. In morphologic examination, the methenamine silver-stained cyst forms of P. carinii from S. araneus and from M. agrestis differed in size from each other and from those originating from laboratory rats. Ultrastructurally P. carinii from S. aruneus did not differ from organisms of rat origin.     

The evolutionary history of the two karyotypic groups of the common shrew, Sorex araneus, in Poland

Genetic variability within and among two karyotypic groups and five chromosome races of the common shrew (Sorex araneus) in Poland was assayed by sequencing a 1023 bp part of the cytochrome b gene (mtDNA) from 28 individuals. Thirty-four variable positions defined 21 distinct haplotypes with a maximum sequence divergence of 0.88%. No significant differentiation in the cytochrome b gene between Western and Eastern Karyotypic groups was found. Haplotype diversity estimates within the races and groups sampled were high (h = 0.800-0.928), while nucleotide diversity estimates were low (pi = 0.0034-0.0053). The distribution of pairwise nucleotide differences fits well with expectations of a "sudden expansion" model. High haplotype diversity was accompanied by relatively high expected heterozygosity (H(E)) values in nuclear genes (calculated over 47 enzyme loci: H(E) = 0.031 - 0.049), giving no evidence for a recent bottleneck after the process of post-Pleistocene recolonization of Poland by the shrews. Thus, for S. araneus chromosome races in Poland, the data on the cytochrome b gene variability support the hypothesis assuming the Robertsonian fusions having spread into an ancestral acrocentric distribution.     

Effects of habitat size and quality on equilibrium density and extinction time of Sorex araneus populations

1. The effects of changes in habitat size and quality on the expected population density and the expected time to extinction of Sorex araneus are studied by means of mathematical models that incorporate demographic stochasticity.
2. Habitat size is characterized by the number of territories, while habitat quality is represented by the expected number of offspring produced during the lifetime of an individual.
3. The expected population density of S. araneus is shown to be mainly influenced by the habitat size. The expected time to extinction of S. araneus populations due to demographic stochasticity, on the other hand, is much more affected by the habitat quality.
4. In a more general setting we demonstrate that, irrespective of the actual species under consideration, the likelihood of extinction as a consequence of demographic stochasticity is more effectively countered by increasing the reproductive success and survival of individuals then by increasing total population size.     

CHROMOSOMAL VERSUS MITOCHONDRIAL DNA EVOLUTION: TRACKING THE EVOLUTIONARY HISTORY OF THE SOUTHWESTERN EUROPEAN POPULATIONS OF THE SOREX ARANEUS GROUP (MAMMALIA,INSECTIVORA)

The shrews of the Sorex araneus group have undergone a spectacular chromosome evolution. The karyotype of Sorex granarius is generally considered ancestral to those of Sorex coronatus and S. araneus. However, a sequence of 777 base pairs of the cytochrome b gene of the mitochondrial DNA (mtDNA) produces a quite different picture: S. granarius is closely related to the populations of S. araneus from the Pyrenees and from the northwestern Alps, whereas S. coronatus and S. araneus from Italy and the southern Alps represent two well-separated lineages. It is suggested that mtDNA and chromosomal evolution are in this case largely independant processes. Whereas mtDNA haplotypes are closely linked to the geographical history of the populations, chromosomal mutations were probably transmitted from one population to another. Available data suggest that the impressive chromosome polymorphism of this group is quite a recent phenomenon.     

Comparative analysis of karyotypes in European shrew species. I. The sibling species Sorex araneus and S. gemellus: Q-bands, G-bands, and position of NORs

The karyotypes of two closely related species of the genus Sorex (Mammalia, Insectivora) were compared with each other by G- and Q-banding techniques and by Ag-AS staining (GOODPASTURE and BLOOM, 1975). By comparing the G-banded karyotypes, it could be ascertained that the basic differences in karyotype between the two species lie in three pericentric inversions, three paracentric inversions, and one reciprocal translocation. This is in near agreement with FORD and HAMERTON (1970), who assumed that both species differ by three pericentric inversions and one tandem translocation. Furthermore, the karyotype of S. araneus (race C) presented by HALKKA et al. (1974) has been compared with the S. araneus of the present report. Considering the species with respect to karyotypic evolution, it is supposed that S. araneus and S. gemellus derive from a common ancestor.     

Three-dimensional organization of interphase fibroblast nuclei in two closely related shrew species (<Emphasis Type="Italic">Sorex granarius</Emphasis> and <Emphasis Type="Italic">Sorex araneus</Emphasis>) differing in the structures of their chromosome termini

A FISH with a probe for telomeric and rDNA repeats and immunofluorescence with ANA CREST and antibodies to nucleolae protein B23 were used to study the three-dimensional (3D) organization of fibroblast interphase nuclei in two shrew twin species, Sorex granarius and Sorex araneus, of the Cordon race. Karyotypes of these species are composed of nearly identical chromosomal arms and differ in the number of their metacentrics and the structures of their terminal chromosome regions. In the short arms of S. granarius, 32 of the acrocentrics have large telomeres that contain an average of 218 kb telomere repeats, which alternate with ribosomal repeats. These regions also contain active nucleolar organizing regions (NORs). In contrast, in active NORs in S. araneus are localized at the terminal regions of 8 chromosomal arms (Zhdanova et al., 2005; 2007b). Here, we show that associations of chromosomes by telomeres and the contact of a part of the telomere clusters with the inner nuclear membrane and nucleolus characterize the interphase nuclei of both Sorex granarius and Sorex araneus. We also reveal the partial colocalization of telomere and ribosomal clusters and the spatial proximity of centomeric and telomeric regions in the interphase nuclei of S. granarius. It appears that only ribosomal clusters containing a sufficient number of active ribosomal genes exhibit a connection with the nucleolus. Nucleolus disassembly during the fibroblastís transition to mitosis and the role of the B23 protein in this process have been studied.     

Molecular variability in the commom shrew Sorex araneus L. from European Russia and Siberia inferred from the length polymorphism of DNA regions flanked by short interspersed elements (Inter-SINE PCR) and the relationships between the Moscow and Seliger chromosome races

Genetic exchange among chromosomal races of the common shrew Sorex araneus and the problem of reproductive barriers have been extensively studied by means of such molecular markers as mtDNA, microsatellites, and allozymes. In the present study, the interpopulation and interracial polymorphism in the common shrew was derived, using fingerprints generated by amplified DNA regions flanked by short interspersed repeats (SINEs)-interSINE PCR (IS-PCR). We used primers, complementary to consensus sequences of two short retroposons: mammalian element MIR and the SOR element from the genome of Sorex araneus. Genetic differentiation among eleven populations of the common shrew from eight chromosome races was estimated. The NP and MJ analyses, as well as multidimensional scaling showed that all samples examined grouped into two main clusters, corresponding to European Russia and Siberia. The bootstrap support of the European Russia cluster in the NJ and MP analyses was respectively 76 and 61%. The bootstrap index for the Siberian cluster was 100% in both analyses; the Tomsk race, included into this cluster, was separated with the bootstrap support of NJ/MP 92/95%.     

Origin and evolution of homologous repeated sequences in the mitochondrial DNA control region of shrews

The complete mitochondrial DNA (mtDNA) control region was amplified and directly sequenced in two species of shrew, Crocidura russula and Sorex araneus (Insectivora, Mammalia). The general organization is similar to that found in other mammals: a central conserved region surrounded by two more variable domains. However, we have found in shrews the simultaneous presence of arrays of tandem repeats in potential locations where repeats tend to occur separately in other mammalian species. These locations correspond to regions which are associated with a possible interruption of the replication processes, either at the end of the three-stranded D-loop structure or toward the end of the heavy-strand replication. In the left domain the repeated sequences (R1 repeats) are 78 bp long, whereas in the right domain the repeats are 12 bp long in C. russula and 14 bp long in S. araneus (R2 repeats). Variation in the copy number of these repeated sequences results in mtDNA control region length differences. Southern blot analysis indicates that level of heteroplasmy (more than one mtDNA form within an individual) differs between species. A comparative study of the R2 repeats in 12 additional species representing three shrew subfamilies provides useful indications for the understanding of the origin and the evolution of these homologous tandemly repeated sequences. An asymmetry in the distribution of variants within the arrays, as well as the constant occurrence of shorter repeated sequences flanking only one side of the R2 arrays, could be related to asymmetry in the replication of each strand of the mtDNA molecule. The pattern of sequence and length variation within and between species, together with the capability of the arrays to form stable secondary structures, suggests that the dominant mechanism involved in the evolution of these arrays in unidirectional replication slippage.