A Neutralizing Antibody Recognizing Primarily N-Linked Glycan Targets the Silent Face of the HIV Envelope.

Virtually the entire surface of the HIV-1-envelope trimer is recognized by neutralizing antibodies, except for a highly glycosylated region at the center of the "silent face" on the gp120 subunit. From an HIV-1-infected donor, #74, we identified antibody VRC-PG05, which neutralized 27% of HIV-1 strains. The crystal structure of the antigen-binding fragment of VRC-PG05 in complex with gp120 revealed an epitope comprised primarily of N-linked glycans from N262, N295, and N448 at the silent face center. Somatic hypermutation occurred preferentially at antibody residues that interacted with these glycans, suggesting somatic development of glycan recognition. Resistance to VRC-PG05 in donor #74 involved shifting of glycan-N448 to N446 or mutation of glycan-proximal residue E293. HIV-1 neutralization can thus be achieved at the silent face center by glycan-recognizing antibody; along with other known epitopes, the VRC-PG05 epitope completes coverage by neutralizing antibody of all major exposed regions of the prefusion closed trimer.

Convergent antibody responses to SARS-CoV-2 in convalescent individuals.

During the coronavirus disease-2019 (COVID-19) pandemic, severe acute respiratory syndrome-related coronavirus-2 (SARS-CoV-2) has led to the infection of millions of people and has claimed hundreds of thousands of lives. The entry of the virus into cells depends on the receptor-binding domain (RBD) of the spike (S) protein of SARS-CoV-2. Although there is currently no vaccine, it is likely that antibodies will be essential for protection. However, little is known about the human antibody response to SARS-CoV-21-5. Here we report on 149 COVID-19-convalescent individuals. Plasma samples collected an average of 39 days after the onset of symptoms had variable half-maximal pseudovirus neutralizing titres; titres were less than 50 in 33% of samples, below 1,000 in 79% of samples and only 1% of samples had titres above 5,000. Antibody sequencing revealed the expansion of clones of RBD-specific memory B cells that expressed closely related antibodies in different individuals. Despite low plasma titres, antibodies to three distinct epitopes on the RBD neutralized the virus with half-maximal inhibitory concentrations (IC50 values) as low as 2 ng ml-1. In conclusion, most convalescent plasma samples obtained from individuals who recover from COVID-19 do not contain high levels of neutralizing activity. Nevertheless, rare but recurring RBD-specific antibodies with potent antiviral activity were found in all individuals tested, suggesting that a vaccine designed to elicit such antibodies could be broadly effective.

PoMeLo: a systematic computational approach to predicting metabolic loss in pathogen genomes.

BACKGROUND: Genome streamlining, the process by which genomes become smaller and encode fewer genes over time, is a common phenomenon among pathogenic bacteria. This reduction is driven by selection for minimized energy expenditure in a nutrient-rich environment. As pathogens evolve to become more reliant on the host, metabolic genes and resulting capabilities are lost in favor of siphoning metabolites from the host. Characterizing genome streamlining, gene loss, and metabolic pathway degradation can be useful in assessing pathogen dependency on host metabolism and identifying potential targets for host-directed therapeutics. RESULTS: PoMeLo (Predictor of Metabolic Loss) is a novel evolutionary genomics-guided computational approach for identifying metabolic gaps in the genomes of pathogenic bacteria. PoMeLo leverages a centralized public database of high-quality genomes and annotations and allows the user to compare an unlimited number of genomes across individual genes and pathways. PoMeLo runs locally using user-friendly prompts in a matter of minutes and generates tabular and visual outputs for users to compare predicted metabolic capacity between groups of bacteria and individual species. Each pathway is assigned a Predicted Metabolic Loss (PML) score to assess the magnitude of genome streamlining. Optionally, PoMeLo places the results in an evolutionary context by including phylogenetic relationships in visual outputs. It can also initially compute phylogenetically-weighted mean genome sizes to identify genome streamlining events. Here, we describe PoMeLo and demonstrate its use in identifying metabolic gaps in genomes of pathogenic Treponema species. CONCLUSIONS: PoMeLo represents an advance over existing methods for identifying metabolic gaps in genomic data, allowing comparison across large numbers of genomes and placing the resulting data in a phylogenetic context. PoMeLo is freely available for academic and non-academic use at https://github.com/czbiohub-sf/pomelo .

Converting non-neutralizing SARS-CoV-2 antibodies into broad-spectrum inhibitors.

Omicron and its subvariants have rendered most authorized monoclonal antibody-based treatments for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) ineffective, highlighting the need for biologics capable of overcoming SARS-CoV-2 evolution. These mostly ineffective antibodies target variable epitopes. Here we describe broad-spectrum SARS-CoV-2 inhibitors developed by tethering the SARS-CoV-2 receptor, angiotensin-converting enzyme 2 (ACE2), to known non-neutralizing antibodies that target highly conserved epitopes in the viral spike protein. These inhibitors, called receptor-blocking conserved non-neutralizing antibodies (ReconnAbs), potently neutralize all SARS-CoV-2 variants of concern (VOCs), including Omicron. Neutralization potency is lost when the linker joining the binding and inhibitory ReconnAb components is severed. In addition, a bi-functional ReconnAb, made by linking ACE2 to a bi-specific antibody targeting two non-overlapping conserved epitopes, defined here, shows sub-nanomolar neutralizing activity against all VOCs, including Omicron and BA.2. Given their conserved targets and modular nature, ReconnAbs have the potential to act as broad-spectrum therapeutics against SARS-CoV-2 and other emerging pandemic diseases.

AIRRscape: An interactive tool for exploring B-cell receptor repertoires and antibody responses.

The sequencing of antibody repertoires of B-cells at increasing coverage and depth has led to the identification of vast numbers of immunoglobulin heavy and light chains. However, the size and complexity of these Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) datasets makes it difficult to perform exploratory analyses. To aid in data exploration, we have developed AIRRscape, an R Shiny-based interactive web browser application that enables B-cell receptor (BCR) and antibody feature discovery through comparisons among multiple repertoires. Using AIRR-seq data as input, AIRRscape starts by aggregating and sorting repertoires into interactive and explorable bins of germline V-gene, germline J-gene, and CDR3 length, providing a high-level view of the entire repertoire. Interesting subsets of repertoires can be quickly identified and selected, and then network topologies of CDR3 motifs can be generated for further exploration. Here we demonstrate AIRRscape using patient BCR repertoires and sequences of published monoclonal antibodies to investigate patterns of humoral immunity to three viral pathogens: SARS-CoV-2, HIV-1, and DENV (dengue virus). AIRRscape reveals convergent antibody sequences among datasets for all three pathogens, although HIV-1 antibody datasets display limited convergence and idiosyncratic responses. We have made AIRRscape available as a web-based Shiny application, along with code on GitHub to encourage its open development and use by immuno-informaticians, virologists, immunologists, vaccine developers, and other scientists that are interested in exploring and comparing multiple immune receptor repertoires.

Impact of past climate warming on genomic diversity and demographic history of collared lemmings across the Eurasian Arctic.

The Arctic climate was warmer than today at the last interglacial and the Holocene thermal optimum. To reveal the impact of past climate-warming events on the demographic history of an Arctic specialist, we examined both mitochondrial and nuclear genomic variation in the collared lemming (Dicrostonyx torquatus, Pallas), a keystone species in tundra communities, across its entire distribution in northern Eurasia. The ancestral phylogenetic position of the West Beringian group and divergence time estimates support the hypothesis of continental range contraction to a single refugial area located in West Beringia during high-magnitude warming of the last interglacial, followed by westward recolonization of northern Eurasia in the last glacial period. The West Beringian group harbors the highest mitogenome diversity and its inferred demography indicates a constantly large effective population size over the Late Pleistocene to Holocene. This suggests that northward forest expansion during recent warming of the Holocene thermal optimum did not affect the gene pool of the collared lemming in West Beringia but reduced genomic diversity and effective population size in all other regions of the Eurasian Arctic. Demographic inference from genomic diversity was corroborated by species distribution modeling showing reduction in species distribution during past climate warming. These conclusions are supported by recent paleoecological evidence suggesting smaller temperature increases and moderate northward forest advances in the extreme northeast of Eurasia during the Late Pleistocene-to-Holocene warming events. This study emphasizes the importance of West Beringia as a potential refugium for cold-adapted Arctic species under ongoing climate warming.

Prediction of phylogeographic endemism in an environmentally complex biome.

Phylogeographic endemism, the degree to which the history of recently evolved lineages is spatially restricted, reflects fundamental evolutionary processes such as cryptic divergence, adaptation and biological responses to environmental heterogeneity. Attempts to explain the extraordinary diversity of the tropics, which often includes deep phylogeographic structure, frequently invoke interactions of climate variability across space, time and topography. To evaluate historical versus contemporary drivers of phylogeographic endemism in a tropical system, we analyse the effects of current and past climatic variation on the genetic diversity of 25 vertebrates in the Brazilian Atlantic rainforest. We identify two divergent bioclimatic domains within the forest and high turnover around the Rio Doce. Independent modelling of these domains demonstrates that endemism patterns are subject to different climatic drivers. Past climate dynamics, specifically areas of relative stability, predict phylogeographic endemism in the north. Conversely, contemporary climatic heterogeneity better explains endemism in the south. These results accord with recent speleothem and fossil pollen studies, suggesting that climatic variability through the last 250 kyr impacted the northern and the southern forests differently. Incorporating sub-regional differences in climate dynamics will enhance our ability to understand those processes shaping high phylogeographic and species endemism, in the Neotropics and beyond.

Whole-genome sequencing rule-out of suspected hospital-onset Rhizopus outbreaks.

Two independent temporal-spatial clusters of hospital-onset Rhizopus infections were evaluated using whole-genome sequencing (WGS). Phylogenetic analysis confirmed that isolates within each cluster were unrelated despite epidemiological suspicion of outbreaks. The ITS1 region alone was insufficient for accurate analysis. WGS has utility for rapid rule-out of suspected nosocomial Rhizopus outbreaks.

MultiSero: An Open-Source Multiplex-ELISA Platform for Measuring Antibody Responses to Infection.

A multiplexed enzyme-linked immunosorbent assay (ELISA) that simultaneously measures antibody binding to multiple antigens can extend the impact of serosurveillance studies, particularly if the assay approaches the simplicity, robustness, and accuracy of a conventional single-antigen ELISA. Here, we report on the development of multiSero, an open-source multiplex ELISA platform for measuring antibody responses to viral infection. Our assay consists of three parts: (1) an ELISA against an array of proteins in a 96-well format; (2) automated imaging of each well of the ELISA array using an open-source plate reader; and (3) automated measurement of optical densities for each protein within the array using an open-source analysis pipeline. We validated the platform by comparing antibody binding to Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) antigens in 217 human sera samples, showing high sensitivity (0.978), specificity (0.977), positive predictive value (0.978), and negative predictive value (0.977) for classifying seropositivity, a high correlation of multiSero determined antibody titers with commercially available SARS-CoV-2 antibody tests, and antigen-specific changes in antibody titer dynamics upon vaccination. The open-source format and accessibility of our multiSero platform can contribute to the adoption of multiplexed ELISA arrays for serosurveillance studies, for SARS-CoV-2 and other pathogens of significance.

Persistence and diversification of the Holarctic shrew, Sorex tundrensis (Family Soricidae), in response to climate change.

Environmental processes govern demography, species movements, community turnover and diversification and yet in many respects these dynamics are still poorly understood at high latitudes. We investigate the combined effects of climate change and geography through time for a widespread Holarctic shrew, Sorex tundrensis. We include a comprehensive suite of closely related outgroup taxa and three independent loci to explore phylogeographic structure and historical demography. We then explore the implications of these findings for other members of boreal communities. The tundra shrew and its sister species, the Tien Shan shrew (Sorex asper), exhibit strong geographic population structure across Siberia and into Beringia illustrating local centres of endemism that correspond to Late Pleistocene refugia. Ecological niche predictions for both current and historical distributions indicate a model of persistence through time despite dramatic climate change. Species tree estimation under a coalescent process suggests that isolation between populations has been maintained across timeframes deeper than the periodicity of Pleistocene glacial cycling. That some species such as the tundra shrew have a history of persistence largely independent of changing climate, whereas other boreal species shifted their ranges in response to climate change, highlights the dynamic processes of community assembly at high latitudes.

multiSero: open multiplex-ELISA platform for analyzing antibody responses to SARS-CoV-2 infection.

Serology has provided valuable diagnostic and epidemiological data on antibody responses to SARS-CoV-2 in diverse patient cohorts. Deployment of high content, multiplex serology platforms across the world, including in low and medium income countries, can accelerate longitudinal epidemiological surveys. Here we report multiSero, an open platform to enable multiplex serology with up to 48 antigens in a 96-well format. The platform consists of three components: ELISA-array of printed proteins, a commercial or home-built plate reader, and modular python software for automated analysis (pysero). We validate the platform by comparing antibody titers against the SARS-CoV-2 Spike, receptor binding domain (RBD), and nucleocapsid (N) in 114 sera from COVID-19 positive individuals and 87 pre-pandemic COVID-19 negative sera. We report data with both a commercial plate reader and an inexpensive, open plate reader (nautilus). Receiver operating characteristic (ROC) analysis of classification with single antigens shows that Spike and RBD classify positive and negative sera with the highest sensitivity at a given specificity. The platform distinguished positive sera from negative sera when the reactivity of the sera was equivalent to the binding of 1 ng mL âˆ'1 RBD-specific monoclonal antibody. We developed normalization and classification methods to pool antibody responses from multiple antigens and multiple experiments. Our results demonstrate a performant and accessible pipeline for multiplexed ELISA ready for multiple applications, including serosurveillance, identification of viral proteins that elicit antibody responses, differential diagnosis of circulating pathogens, and immune responses to vaccines.

Adaption of a conventional ELISA to a 96-well ELISA-Array for measuring the antibody responses to influenza virus proteins and vaccines.

We describe an adaptation of conventional ELISA methods to an ELISA-Array format using non-contact Piezo printing of up to 30 spots of purified recombinant viral fusion proteins and vaccine on 96 well high-protein binding plates. Antigens were printed in 1 nanoliter volumes of protein stabilizing buffer using as little as 0.25 nanograms of protein, 2000-fold less than conventional ELISA. The performance of the ELISA-Array was demonstrated by serially diluting n = 9 human post-flu vaccination plasma samples starting at a 1/1000 dilution and measuring binding to the array of Influenza antigens. Plasma polyclonal antibody levels were detected using a cocktail of biotinylated anti-human kappa and lambda light chain antibodies, followed by a Streptavidin-horseradish peroxidase conjugate and the dose-dependent signal was developed with a precipitable TMB substrate. Intra- and inter-assay precision of absorbance units among the eight donor samples showed mean CVs of 4.8% and 10.8%, respectively. The plasma could be differentiated by donor and antigen with titer sensitivities ranging from 1 × 103 to 4 × 106, IC50 values from 1 × 104 to 9 × 106, and monoclonal antibody sensitivities in the ng/mL range. Equivalent sensitivities of ELISA versus ELISA-Array, compared using plasma and an H1N1 HA trimer, were achieved on the ELISA-Array printed at 0.25 ng per 200um spot and 1000 ng per ELISA 96-well. Vacuum-sealed array plates were shown to be stable when stored for at least 2 days at ambient temperature and up to 1 month at 4-8 °C. By the use of any set of printed antigens and analyte matrices the methods of this multiplexed ELISA-Array format can be broadly applied in translational research.

VSV-Displayed HIV-1 Envelope Identifies Broadly Neutralizing Antibodies Class-Switched to IgG and IgA.

The HIV-1 envelope (Env) undergoes conformational changes during infection. Broadly neutralizing antibodies (bNAbs) are typically isolated by using soluble Env trimers, which do not capture all Env states. To address these limitations, we devised a vesicular stomatitis virus (VSV)-based probe to display membrane-embedded Env trimers and isolated five bNAbs from two chronically infected donors, M4008 and M1214. Donor B cell receptor (BCR) repertoires identified two bNAb lineages, M4008_N1 and M1214_N1, that class-switched to immunoglobulin G (IgG) and IgA. Variants of these bNAbs reconstituted as IgA demonstrated broadly neutralizing activity, and the IgA fraction of M1214 plasma conferred neutralization. M4008_N1 epitope mapping revealed a glycan-independent V3 epitope conferring tier 2 virus neutralization. A 4.86-Å-resolution cryogenic electron microscopy (cryo-EM) structure of M1214_N1 complexed with CH505 SOSIP revealed another elongated epitope, the V2V5 corridor, extending from V2 to V5. Overall, the VSVENV probe identified bNAb lineages with neutralizing IgG and IgA members targeting distinct sites of HIV-1 Env vulnerability.

Convergent Antibody Responses to SARS-CoV-2 Infection in Convalescent Individuals.

During the COVID-19 pandemic, SARS-CoV-2 infected millions of people and claimed hundreds of thousands of lives. Virus entry into cells depends on the receptor binding domain (RBD) of the SARS-CoV-2 spike protein (S). Although there is no vaccine, it is likely that antibodies will be essential for protection. However, little is known about the human antibody response to SARS-CoV-21-5. Here we report on 149 COVID-19 convalescent individuals. Plasmas collected an average of 39 days after the onset of symptoms had variable half-maximal neutralizing titers ranging from undetectable in 33% to below 1:1000 in 79%, while only 1% showed titers >1:5000. Antibody cloning revealed expanded clones of RBD-specific memory B cells expressing closely related antibodies in different individuals. Despite low plasma titers, antibodies to three distinct epitopes on RBD neutralized at half-maximal inhibitory concentrations (IC50s) as low as single digit ng/mL. Thus, most convalescent plasmas obtained from individuals who recover from COVID-19 do not contain high levels of neutralizing activity. Nevertheless, rare but recurring RBD-specific antibodies with potent antiviral activity were found in all individuals tested, suggesting that a vaccine designed to elicit such antibodies could be broadly effective.

Functional Enrichment and Analysis of Antigen-Specific Memory B Cell Antibody Repertoires in PBMCs.

Phenotypic screening of antigen-specific antibodies in human blood is a common diagnostic test for infectious agents and a correlate of protection after vaccination. In addition to long-lived antibody secreting plasma cells residing in the bone marrow, memory B cells are a latent source of antigen-experienced, long-term immunity that can be found at low frequencies in circulating peripheral blood mononuclear cells (PBMCs). Assessing the genotype, clonal frequency, quality, and function of antibodies resulting from an individual's persistent memory B cell repertoire can help inform the success or failure of immune protection. Using in vitro polyclonal stimulation, we functionally expand the memory repertoire from PBMCs and clonally map monoclonal antibodies from this population. We show that combining deep sequencing of stimulated memory B cell repertoires with retrieving single antigen-specific cells is a promising approach in evaluating the latent, functional B cell memory in PBMCs.

Broadly neutralizing human antibodies against dengue virus identified by single B cell transcriptomics.

Eliciting broadly neutralizing antibodies (bNAbs) against the four dengue virus serotypes (DENV1-4) that are spreading into new territories is an important goal of vaccine design. To define bNAb targets, we characterized 28 antibodies belonging to expanded and hypermutated clonal families identified by transcriptomic analysis of single plasmablasts from DENV-infected individuals. Among these, we identified J9 and J8, two somatically related bNAbs that potently neutralized DENV1-4. Mutagenesis studies showed that the major recognition determinants of these bNAbs are in E protein domain I, distinct from the only known class of human bNAbs against DENV with a well-defined epitope. B cell repertoire analysis from acute-phase peripheral blood suggested that J9 and J8 followed divergent somatic hypermutation pathways, and that a limited number of mutations was sufficient for neutralizing activity. Our study suggests multiple B cell evolutionary pathways leading to DENV bNAbs targeting a new epitope that can be exploited for vaccine design.

5' Rapid Amplification of cDNA Ends and Illumina MiSeq Reveals B Cell Receptor Features in Healthy Adults, Adults With Chronic HIV-1 Infection, Cord Blood, and Humanized Mice.

Using 5' rapid amplification of cDNA ends, Illumina MiSeq, and basic flow cytometry, we systematically analyzed the expressed B cell receptor (BCR) repertoire in 14 healthy adult PBMCs, 5 HIV-1+ adult PBMCs, 5 cord blood samples, and 3 HIS-CD4/B mice, examining the full-length variable region of µ, γ, α, κ, and λ chains for V-gene usage, somatic hypermutation (SHM), and CDR3 length. Adding to the known repertoire of healthy adults, Illumina MiSeq consistently detected small fractions of reads with high mutation frequencies including hypermutated µ reads, and reads with long CDR3s. Additionally, the less studied IgA repertoire displayed similar characteristics to that of IgG. Compared to healthy adults, the five HIV-1 chronically infected adults displayed elevated mutation frequencies for all µ, γ, α, κ, and λ chains examined and slightly longer CDR3 lengths for γ, α, and λ. To evaluate the reconstituted human BCR sequences in a humanized mouse model, we analyzed cord blood and HIS-CD4/B mice, which all lacked the typical SHM seen in the adult reference. Furthermore, MiSeq revealed identical unmutated IgM sequences derived from separate cell aliquots, thus for the first time demonstrating rare clonal members of unmutated IgM B cells by sequencing.

Evolutionary dynamics of intron size, genome size, and physiological correlates in archosaurs.

It has been proposed that intron and genome sizes in birds are reduced in comparison with mammals because of the metabolic demands of flight. To test this hypothesis, we examined the sizes of 14 introns in a nonflying relative of birds, the American alligator (Alligator mississippiensis), and in 19 flighted and flightless birds in 12 taxonomic orders. Our results indicate that a substantial fraction (66%) of the reduction in intron size as well as in genome size had already occurred in nonflying archosaurs. Using phylogenetically independent contrasts, we found that the proposed inverse correlation of genome size and basal metabolic rate (BMR) is significant among amniotes and archosaurs, whereas intron and genome size variation within birds showed no significant correlation with BMR. We show statistically that the distribution of genome sizes in birds and mammals is underdispersed compared with the Brownian motion model and consistent with strong stabilizing selection; that genome size differences between vertebrate clades are overdispersed and punctuational; and that evolution of BMR and avian intron size is consistent with Brownian motion. These results suggest that the contrast between genome size/BMR and intron size/BMR correlations may be a consequence of different intensities of selection for these traits and that we should not expect changes in intron size to be significantly associated with metabolically costly behaviors such as flight.

Losing ground: past history and future fate of Arctic small mammals in a changing climate.

According to the IPCC, the global average temperature is likely to increase by 1.4-5.8 °C over the period from 1990 to 2100. In Polar regions, the magnitude of such climatic changes is even larger than in temperate and tropical biomes. This amplified response is particularly worrisome given that the so-far moderate warming is already impacting Arctic ecosystems. Predicting species responses to rapid warming in the near future can be informed by investigating past responses, as, like the rest of the planet, the Arctic experienced recurrent cycles of temperature increase and decrease (glacial-interglacial changes) in the past. In this study, we compare the response of two important prey species of the Arctic ecosystem, the collared lemming and the narrow-skulled vole, to Late Quaternary climate change. Using ancient DNA and Ecological Niche Modeling (ENM), we show that the two species, which occupy similar, but not identical ecological niches, show markedly different responses to climatic and environmental changes within broadly similar habitats. We empirically demonstrate, utilizing coalescent model-testing approaches, that collared lemming populations decreased substantially after the Last Glacial Maximum; a result consistent with distributional loss over the same period based on ENM results. Given this strong association, we projected the current niche onto future climate conditions based on IPCC 4.0 scenarios, and forecast accelerating loss of habitat along southern range boundaries with likely associated demographic consequences. Narrow-skulled vole distribution and demography, by contrast, was only moderately impacted by past climatic changes, but predicted future changes may begin to affect their current western range boundaries. Our work, founded on multiple lines of evidence suggests a future of rapidly geographically shifting Arctic small mammal prey communities, some of whom are on the edge of existence, and whose fate may have ramifications for the whole Arctic food web and ecosystem.

Effects of late quaternary climate change on Palearctic shrews.

The Late Quaternary was a time of rapid climatic oscillations and drastic environmental changes. In general, species can respond to such changes by behavioral accommodation, distributional shifts, ecophenotypic modifications (nongenetic), evolution (genetic) or ultimately face local extinction. How those responses manifested in the past is essential for properly predicting future ones especially as the current warm phase is further intensified by rising levels of atmospheric carbon dioxide. Here, we use ancient DNA (aDNA) and morphological features in combination with ecological niche modeling (ENM) to investigate genetic and nongenetic responses of Central European Palearctic shrews to past climatic change. We show that a giant form of shrew, previously described as an extinct Pleistocene Sorex species, represents a large ecomorph of the common shrew (Sorex araneus), which was replaced by populations from a different gene-pool and with different morphology after the Pleistocene Holocene transition. We also report the presence of the cold-adapted tundra shrew (S. tundrensis) in Central Europe. This species is currently restricted to Siberia and was hitherto unknown as an element of the Pleistocene fauna of Europe. Finally, we show that there is no clear correlation between climatic oscillations within the last 50 000 years and body size in shrews and conclude that a special nonanalogous situation with regard to biodiversity and food supply in the Late Glacial may have caused the observed large body size.