Uncovering Broad Macroecological Patterns by Comparing the Shape of Species' Distributions along Environmental Gradients.

AbstractSpecies' distributions can take many different forms. For example, fat-tailed or skewed distributions are very common in nature, as these can naturally emerge as a result of individual variability and asymmetric environmental tolerances, respectively. Studying the basic shape of distributions can teach us a lot about the ways climatic processes and historical contingencies shape ecological communities. Yet we still lack a general understanding of how their shapes and properties compare to each other along gradients. Here, we use Bayesian nonlinear models to quantify range shape properties in empirical plant distributions. With this approach, we are able to distil the shape of plant distributions and compare them along gradients and across species. Studying the relationship between distribution properties, we revealed the existence of broad macroecological patterns along environmental gradients-such as those expected from Rapoport's rule and the abiotic stress limitation hypothesis. We also find that some aspects of the shape of observed ranges-such as kurtosis and skewness of the distributions-could be intrinsic properties of species or the result of their historical contexts. Overall, our modeling approach and results untangle the general shape of plant distributions and provide a mapping of how this changes along environmental gradients.

Transmission mode and dispersal traits correlate with host specificity in mammalian gut microbes.

Different host species associate with distinct gut microbes in mammals, a pattern sometimes referred to as phylosymbiosis. However, the processes shaping this host specificity are not well understood. One model proposes that barriers to microbial transmission promote specificity by limiting microbial dispersal between hosts. This model predicts that specificity levels measured across microbes is correlated to transmission mode (vertical vs. horizontal) and individual dispersal traits. Here, we leverage two large publicly available gut microbiota data sets (1490 samples from 195 host species) to test this prediction. We found that host specificity varies widely across bacteria (i.e., there are generalist and specialist bacteria) and depends on transmission mode and dispersal ability. Horizontally-like transmitted bacteria equipped with traits that facilitate switches between host (e.g., tolerance to oxygen) were found to be less specific (more generalist) than microbes without those traits, for example, vertically-like inherited bacteria that are intolerant to oxygen. Altogether, our findings are compatible with a model in which limited microbial dispersal abilities foster host specificity.

Taxonomic and functional biogeographies of soil bacterial communities across the Tibet plateau are better explained by abiotic conditions than distance and plant community composition.

The processes governing soil bacteria biogeography are still not fully understood. It remains unknown how the importance of environmental filtering and dispersal differs between bacterial taxonomic and functional biogeography, and whether their importance is scale-dependent. We sampled soils across the Tibet plateau, with distances among plots ranging from 20 m to 1550 km. Taxonomic composition of bacterial community was characterized by 16S amplicon sequencing and functional community composition by qPCR targeting 9 functional groups involved in N dynamics. Factors representing climate, soil and plant community were measured to assess different facets of environmental dissimilarity. Both bacterial taxonomic and functional dissimilarities were more related to abiotic dissimilarity than biotic (vegetation) dissimilarity or distance. Taxonomic dissimilarity was mostly explained by differences in soil pH and mean annual temperature (MAT), while functional dissimilarity was linked to differences in soil N and P availabilities and N:P ratio. Soil pH and MAT remained the main determinants of taxonomic dissimilarity across spatial scales. In contrast, the explanatory variables of N-related functional dissimilarity varied across the scales, with soil moisture and organic matter having the highest role across short distances (<~330 km), and available P, N:P ratio and distance being important over long distances (>~660 km). Our results demonstrate how biodiversity dimension (taxonomic versus functional aspects) and spatial scale influence the factors driving soil bacterial biogeography.

Catchment-based sampling of river eDNA integrates terrestrial and aquatic biodiversity of alpine landscapes.

Monitoring of terrestrial and aquatic species assemblages at large spatial scales based on environmental DNA (eDNA) has the potential to enable evidence-based environmental policymaking. The spatial coverage of eDNA-based studies varies substantially, and the ability of eDNA metabarcoding to capture regional biodiversity remains to be assessed; thus, questions about best practices in the sampling design of entire landscapes remain open. We tested the extent to which eDNA sampling can capture the diversity of a region with highly heterogeneous habitat patches across a wide elevation gradient for five days through multiple hydrological catchments of the Swiss Alps. Using peristaltic pumps, we filtered 60 L of water at five sites per catchment for a total volume of 1800 L. Using an eDNA metabarcoding approach focusing on vertebrates and plants, we detected 86 vertebrate taxa spanning 41 families and 263 plant taxa spanning 79 families across ten catchments. For mammals, fishes, amphibians and plants, the detected taxa covered some of the most common species in the region according to long-term records while including a few more rare taxa. We found marked turnover among samples from distinct elevational classes indicating that the biological signal in alpine rivers remains relatively localised and is not aggregated downstream. Accordingly, species compositions differed between catchments and correlated with catchment-level forest and grassland cover. Biomonitoring schemes based on capturing eDNA across rivers within biologically integrated catchments may pave the way toward a spatially comprehensive estimation of biodiversity.

Soil protist function varies with elevation in the Swiss Alps.

Protists are abundant and play key trophic functions in soil. Documenting how their trophic contributions vary across large environmental gradients is essential to understand and predict how biogeochemical cycles will be impacted by global changes. Here, using amplicon sequencing of environmental DNA in open habitat soil from 161 locations spanning 2600 m of elevation in the Swiss Alps (from 400 to 3000 m), we found that, over the whole study area, soils are dominated by consumers, followed by parasites and phototrophs. In contrast, the proportion of these groups in local communities shows large variations in relation to elevation. While there is, on average, three times more consumers than parasites at low elevation (400-1000 m), this ratio increases to 12 at high elevation (2000-3000 m). This suggests that the decrease in protist host biomass and diversity toward mountains tops impact protist functional composition. Furthermore, the taxonomic composition of protists that infect animals was related to elevation while that of protists that infect plants or of protist consumers was related to soil pH. This study provides a first step to document and understand how soil protist functions vary along the elevational gradient.

Setting priorities for climate change adaptation of Critical Sites in the Africa-Eurasian waterbird flyways.

Despite their importance for biodiversity and ecosystem services, wetlands are among the most threatened ecosystems globally. The conservation of many migratory waterbirds depends on the conservation of a network of key sites along their flyways. However, the suitability of these sites is changing under climate change, and it is important that management of individual sites in the network adapts to these changes. Using bioclimatic models that also account for changes in inundation, we found that projected climate change will reduce habitat suitability for waterbirds at 57.5% of existing Critical Sites within Africa-Eurasia, varying from 20.1% in Eastern Europe to 87.0% in Africa. African and Middle East sites are particularly threatened, comprising 71 of the 100 most vulnerable sites. By highlighting priority sites for conservation and classifying Critical Sites into Climate Change Adaptation Strategy (CCAS) classes, our results can be used to support the climate change adaptation of both individual sites and the entire site network.

Data integration methods to account for spatial niche truncation effects in regional projections of species distribution.

Many species distribution models (SDMs) are built with precise but geographically restricted presence-absence data sets (e.g., a country) where only a subset of the environmental conditions experienced by a species across its range is considered (i.e., spatial niche truncation). This type of truncation is worrisome because it can lead to incorrect predictions e.g., when projecting to future climatic conditions belonging to the species niche but unavailable in the calibration area. Data from citizen-science programs, species range maps or atlases covering the full species range can be used to capture those parts of the species' niche that are missing regionally. However, these data usually are too coarse or too biased to support regional management. Here, we aim to (1) demonstrate how varying degrees of spatial niche truncation affect SDMs projections when calibrated with climatically truncated regional data sets and (2) test the performance of different methods to harness information from larger-scale data sets presenting different spatial resolutions to solve the spatial niche truncation problem. We used simulations to compare the performance of the different methods, and applied them to a real data set to predict the future distribution of a plant species (Potentilla aurea) in Switzerland. SDMs calibrated with geographically restricted data sets expectedly provided biased predictions when projected outside the calibration area or time period. Approaches integrating information from larger-scale data sets using hierarchical data integration methods usually reduced this bias. However, their performance varied depending on the level of spatial niche truncation and how data were combined. Interestingly, while some methods (e.g., data pooling, downscaling) performed well on both simulated and real data, others (e.g., those based on a Poisson point process) performed better on real data, indicating a dependency of model performance on the simulation process (e.g., shape of simulated response curves). Based on our results, we recommend to use different data integration methods and, whenever possible, to make a choice depending on model performance. In any case, an ensemble modeling approach can be used to account for uncertainty in how niche truncation is accounted for and identify areas where similarities/dissimilarities exist across methods.

Eco-genetic additivity of diploids in allopolyploid wild wheats.

Underpinnings of the distribution of allopolyploid species (hybrids with duplicated genome) along spatial and ecological gradients are elusive. As allopolyploid speciation combines the range of genetic and ecological characteristics of divergent diploids, allopolyploids initially show their additivity and are predicted to evolve differentiated ecological niches to establish in face of their competition. Here, we use four diploid wild wheats that differentially combined into four independent allopolyploid species to test for such additivity and assess the impact of ecological constraints on species ranges. Divergent genetic variation from diploids being fixed in heterozygote allopolyploids supports their genetic additivity. Spatial integration of comparative phylogeography and modelling of climatic niches supports ecological additivity of locally adapted diploid progenitors into allopolyploid species which subsequently colonised wide ranges. Allopolyploids fill suitable range to a larger extent than diploids and conservative evolution following the combination of divergent species appears to support their expansion under environmental changes.

Greater topoclimatic control of above- versus below-ground communities.

Assessing the degree to which climate explains the spatial distributions of different taxonomic and functional groups is essential for anticipating the effects of climate change on ecosystems. Most effort so far has focused on above-ground organisms, which offer only a partial view on the response of biodiversity to environmental gradients. Here including both above- and below-ground organisms, we quantified the degree of topoclimatic control on the occurrence patterns of >1,500 taxa and phylotypes along a c. 3,000 m elevation gradient, by fitting species distribution models. Higher model performances for animals and plants than for soil microbes (fungi, bacteria and protists) suggest that the direct influence of topoclimate is stronger on above-ground species than on below-ground microorganisms. Accordingly, direct climate change effects are predicted to be stronger for above-ground than for below-ground taxa, whereas factors expressing local soil microclimate and geochemistry are likely more important to explain and forecast the occurrence patterns of soil microbiota. Detailed mapping and future scenarios of soil microclimate and microhabitats, together with comparative studies of interacting and ecologically dependent above- and below-ground biota, are thus needed to understand and realistically forecast the future distribution of ecosystems.

Integrated Methods for Monitoring the Invasive Potential and Management of Heracleum mantegazzianum (giant hogweed) in Switzerland.

Biological invasions are a major driver of human-induced global environmental change. This makes monitoring of potential spread, population changes and control measures necessary for guiding management. We illustrate the value of integrated methods (species distribution modelling (SDM), plant population monitoring and questionnaires) for monitoring and assessing invasions of Heracleum mantegazzianum (giant hogweed) over time in Switzerland. SDMs highlighted the potential spread of the species, uncovered ecological mechanisms underlying invasions and guided monitoring at a regional level. We used adaptive and repeat plant sampling to monitor invasive population status and changes, and assess the effectiveness of H. mantegazzianum management over three periods (2005, 2013 and 2018) within the pre-Alps, Vaud. We also conducted questionnaire surveys with managers and the public. Multiscale modelling, and integrating global and regional SDMs, provided the best predictions, showing that H. mantegazzianum can potentially invade large parts of Switzerland, especially below 2 000 m a.s.l. Over time, populations of invasive H. mantegazzianum in the Vaud pre-Alps have declined, which is most likely due to a sharp rise in management uptake post 2007 (7% of municipalities before 2007 to 86% in 2018). The level of known invasive populations has decreased by 54% over time. Some municipalities have even successfully eradicated H. mantegazzianum within their borders. However, a few areas, particularly in the rural, higher-altitude municipalities, where management was not implemented effectively, populations have expanded, which could hamper control efforts at lower altitudes. We provide encouraging evidence that control measures can be effective in reducing plant invasions with long-term commitment, as well as a good template for using integrated methodological approaches to better study and monitor invasive alien species.

Genetic diversity in frogs linked to past and future climate changes on the roof of the world.

Mountains, representing storehouses of biodiversity, endemism and threatened species, are biodiversity hotspots of great conservation importance. However, increasing evidence indicates that mountain species throughout the world are responding to climate change, past or contemporary, by shifting their geographic distributions and patterns of genetic diversity, potentially affecting their adaptive capacity and increasing risk of extinction. Using the iconic high-elevation frog Nanorana parkeri as indicator, we showed how spatial analyses of climatic stability combined with genetic data allow unravelling amphibian responses to past and future climate changes on 'the roof of the world'-the Qinghai-Tibetan Plateau. We found that areas along the Yarlung Tsangpo Valley were climatically more stable relative to other regions, apparently serving as a large climatic refugium during Quaternary glaciations, but that these areas will likely be affected by future climate change. As populations closer to Quaternary refugia usually had higher genetic diversity, current genetic diversity can be explained in the largest part by distance to historically stable areas, outweighing other historical and contemporary factors. Along with the dynamics of suitable range, a fluctuating habitat fragmentation supported the pattern of historical changes in genetic diversity (Ne ) over time. Our results emphasize strong relationships between amphibian genetic diversity, past range dynamics and where to preserve suitable habitats in the face of future climate changes. More generally, our findings highlighted a central role of refugia during Quaternary climatic fluctuations, and how isolation from refugia may have modulated amphibian genetic diversity across the Qinghai-Tibetan Plateau.

Reconstructing geographical parthenogenesis: effects of niche differentiation and reproductive mode on Holocene range expansion of an alpine plant.

Asexual taxa often have larger ranges than their sexual progenitors, particularly in areas affected by Pleistocene glaciations. The reasons given for this 'geographical parthenogenesis' are contentious, with expansion of the ecological niche or colonisation advantages of uniparental reproduction assumed most important in case of plants. Here, we parameterized a spread model for the alpine buttercup Ranunculus kuepferi and reconstructed the joint Holocene range expansion of its sexual and apomictic cytotype across the European Alps under different simulation settings. We found that, rather than niche broadening or a higher migration rate, a shift of the apomict's niche towards colder conditions per se was crucial as it facilitated overcoming of topographical barriers, a factor likely relevant for many alpine apomicts. More generally, our simulations suggest potentially strong interacting effects of niche differentiation and reproductive modes on range formation of related sexual and asexual taxa arising from their differential sensitivity to minority cytotype disadvantage.

Less favourable climates constrain demographic strategies in plants.

Correlative species distribution models are based on the observed relationship between species' occurrence and macroclimate or other environmental variables. In climates predicted less favourable populations are expected to decline, and in favourable climates they are expected to persist. However, little comparative empirical support exists for a relationship between predicted climate suitability and population performance. We found that the performance of 93 populations of 34 plant species worldwide - as measured by in situ population growth rate, its temporal variation and extinction risk - was not correlated with climate suitability. However, correlations of demographic processes underpinning population performance with climate suitability indicated both resistance and vulnerability pathways of population responses to climate: in less suitable climates, plants experienced greater retrogression (resistance pathway) and greater variability in some demographic rates (vulnerability pathway). While a range of demographic strategies occur within species' climatic niches, demographic strategies are more constrained in climates predicted to be less suitable.

Genetic consequences of Quaternary climatic oscillations in the Himalayas: Primula tibetica as a case study based on restriction site-associated DNA sequencing.

The effects of Quaternary climatic oscillations on the demography of organisms vary across regions and continents. In taxa distributed in Europe and North America, several paradigms regarding the distribution of refugia have been identified. By contrast, less is known about the processes that shaped the species' spatial genetic structure in areas such as the Himalayas, which is considered a biodiversity hotspot. Here, we investigated the phylogeographic structure and population dynamics of Primula tibetica by combining genomic phylogeography and species distribution models (SDMs). Genomic data were obtained for 293 samples of P. tibetica using restriction site-associated DNA sequencing (RADseq). Ensemble SDMs were carried out to predict potential present and past distribution ranges. Four distinct lineages were identified. Approximate Bayesian computation analyses showed that each of them have experienced both expansions and bottlenecks since their divergence, which occurred during or across the Quaternary glacial cycles. The two lineages at both edges of the distribution were found to be more vulnerable and responded in different ways to past climatic changes. These results illustrate how past climatic changes affected the demographic history of Himalayan organisms. Our findings highlight the significance of combining genomic approaches with environmental data when evaluating the effects of past climatic changes.

Numerical ragweed pollen forecasts using different source maps: a comparison for France.

One of the key input parameters for numerical pollen forecasts is the distribution of pollen sources. Generally, three different methodologies exist to assemble such distribution maps: (1) plant inventories, (2) land use data in combination with annual pollen counts, and (3) ecological modeling. We have used six exemplary maps for all of these methodologies to study their applicability and usefulness in numerical pollen forecasts. The ragweed pollen season of 2012 in France has been simulated with the numerical weather prediction model COSMO-ART using each of the distribution maps in turn. The simulated pollen concentrations were statistically compared to measured values to derive a ranking of the maps with respect to their performance. Overall, approach (2) resulted in the best correspondence between observed and simulated pollen concentrations for the year 2012. It is shown that maps resulting from ecological modeling that does not include a sophisticated estimation of the plant density have a very low predictive skill. For inventory maps and the maps based on land use data and pollen counts, the results depend very much on the observational site. The use of pollen counts to calibrate the map enhances the performance of the model considerably.

Local Environmental Factors Drive Divergent Grassland Soil Bacterial Communities in the Western Swiss Alps.

Mountain ecosystems are characterized by a diverse range of climatic and topographic conditions over short distances and are known to shelter a high biodiversity. Despite important progress, still little is known on bacterial diversity in mountain areas. Here, we investigated soil bacterial biogeography at more than 100 sampling sites randomly stratified across a 700-km2 area with 2,200-m elevation gradient in the western Swiss Alps. Bacterial grassland communities were highly diverse, with 12,741 total operational taxonomic units (OTUs) across 100 sites and an average of 2,918 OTUs per site. Bacterial community structure was correlated with local climatic, topographic, and soil physicochemical parameters with high statistical significance. We found pH (correlated with % CaO and % mineral carbon), hydrogen index (correlated with bulk gravimetric water content), and annual average number of frost days during the growing season to be among the groups of the most important environmental drivers of bacterial community structure. In contrast, bacterial community structure was only weakly stratified as a function of elevation. Contrasting patterns were discovered for individual bacterial taxa. Acidobacteria responded both positively and negatively to pH extremes. Various families within the Bacteroidetes responded to available phosphorus levels. Different verrucomicrobial groups responded to electrical conductivity, total organic carbon, water content, and mineral carbon contents. Alpine grassland bacterial communities are thus highly diverse, which is likely due to the large variety of different environmental conditions. These results shed new light on the biodiversity of mountain ecosystems, which were already identified as potentially fragile to anthropogenic influences and climate change. IMPORTANCE: This article addresses the question of how microbial communities in alpine regions are dependent on local climatic and soil physicochemical variables. We benefit from a unique 700-km2 study region in the western Swiss Alps region, which has been exhaustively studied for macro-organismal and fungal ecology, and for topoclimatic modeling of future ecological trends, but without taking into account soil bacterial diversity. Here, we present an in-depth biogeographical characterization of the bacterial community diversity in this alpine region across 100 randomly stratified sites, using 56 environmental variables. Our exhaustive sampling ensured the detection of ecological trends with high statistical robustness. Our data both confirm previously observed general trends and show many new detailed trends for a wide range of bacterial taxonomic groups and environmental parameters.

Imprints of multiple glacial refugia in the Pyrenees revealed by phylogeography and palaeodistribution modelling of an endemic spider.

Mediterranean mountain ranges harbour highly endemic biota in islandlike habitats. Their topographic diversity offered the opportunity for mountain species to persist in refugial areas during episodes of major climatic change. We investigate the role of Quaternary climatic oscillations in shaping the demographic history and distribution ranges in the spider Harpactocrates ravastellus, endemic to the Pyrenees. Gene trees and multispecies coalescent analyses on mitochondrial and nuclear DNA sequences unveiled two distinct lineages with a hybrid zone around the northwestern area of the Catalan Pyrenees. The lineages were further supported by morphological differences. Climatic niche-based species distribution models (SDMs) identified two lowland refugia at the western and eastern extremes of the mountain range, which would suggest secondary contact following postglacial expansion of populations from both refugia. Neutrality test and approximate Bayesian computation (ABC) analyses indicated that several local populations underwent severe bottlenecks followed by population expansions, which in combination with the deep population differentiation provided evidence for population survival during glacial periods in microrefugia across the mountain range, in addition to the main Atlantic and Mediterranean (western and eastern) refugia. This study sheds light on the complexities of Quaternary climatic oscillations in building up genetic diversity and local endemicity in the southern Europe mountain ranges.

Will climate change increase the risk of plant invasions into mountains?

Mountain ecosystems have been less adversely affected by invasions of non-native plants than most other ecosystems, partially because most invasive plants in the lowlands are limited by climate and cannot grow under harsher high-elevation conditions. However, with ongoing climate change, invasive species may rapidly move upwards and threaten mid-, and then high-elevation mountain ecosystems. We evaluated this threat by modeling the current and future habitat suitability for 48 invasive plant species in Switzerland and New South Wales, Australia. Both regions had contrasting climate interactions with elevation, resulting in possible different responses of species distributions to climate change. Using a species distribution modeling approach that combines data from two spatial scales, we built high-resolution species distribution models (≤ 250 m) that account for the global climatic niche of species and also finer variables depicting local climate and disturbances. We found that different environmental drivers limit the elevation range of invasive species in each of the two regions, leading to region-specific species responses to climate change. The optimal suitability for plant invaders is predicted to markedly shift from the lowland to the montane or subalpine zone in Switzerland, whereas the upward shift is far less pronounced in New South Wales where montane and subalpine elevations are already suitable. The results suggest that species most likely to invade high elevations in Switzerland will be cold-tolerant, whereas species with an affinity to moist soils are most likely to invade higher elevations in Australia. Other plant traits were only marginally associated with elevation limits. These results demonstrate that a more systematic consideration of future distributions of invasive species is required in conservation plans of not yet invaded mountainous ecosystems.

Tree cover at fine and coarse spatial grains interacts with shade tolerance to shape plant species distributions across the Alps.

The role of competition for light among plants has long been recognised at local scales, but its importance for plant species distributions at larger spatial scales has generally been ignored. Tree cover modifies the local abiotic conditions below the canopy, notably by reducing light availability, and thus, also the performance of species that are not adapted to low-light conditions. However, this local effect may propagate to coarser spatial grains, by affecting colonisation probabilities and local extinction risks of herbs and shrubs. To assess the effect of tree cover at both the plot- and landscape-grain sizes (approximately 10-m and 1-km), we fit Generalised Linear Models (GLMs) for the plot-level distributions of 960 species of herbs and shrubs using 6,935 vegetation plots across the European Alps. We ran four models with different combinations of variables (climate, soil and tree cover) at both spatial grains for each species. We used partial regressions to evaluate the independent effects of plot- and landscape-grain tree cover on plot-level plant communities. Finally, the effects on species-specific elevational range limits were assessed by simulating a removal experiment comparing the species distributions under high and low tree cover. Accounting for tree cover improved the model performance, with the probability of the presence of shade-tolerant species increasing with increasing tree cover, whereas shade-intolerant species showed the opposite pattern. The tree cover effect occurred consistently at both the plot and landscape spatial grains, albeit most strongly at the former. Importantly, tree cover at the two grain sizes had partially independent effects on plot-level plant communities. With high tree cover, shade-intolerant species exhibited narrower elevational ranges than with low tree cover whereas shade-tolerant species showed wider elevational ranges at both limits. These findings suggest that forecasts of climate-related range shifts for herb and shrub species may be modified by tree cover dynamics.

Density-based hierarchical clustering of pyro-sequences on a large scale--the case of fungal ITS1.

MOTIVATION: Analysis of millions of pyro-sequences is currently playing a crucial role in the advance of environmental microbiology. Taxonomy-independent, i.e. unsupervised, clustering of these sequences is essential for the definition of Operational Taxonomic Units. For this application, reproducibility and robustness should be the most sought after qualities, but have thus far largely been overlooked. RESULTS: More than 1 million hyper-variable internal transcribed spacer 1 (ITS1) sequences of fungal origin have been analyzed. The ITS1 sequences were first properly extracted from 454 reads using generalized profiles. Then, otupipe, cd-hit-454, ESPRIT-Tree and DBC454, a new algorithm presented here, were used to analyze the sequences. A numerical assay was developed to measure the reproducibility and robustness of these algorithms. DBC454 was the most robust, closely followed by ESPRIT-Tree. DBC454 features density-based hierarchical clustering, which complements the other methods by providing insights into the structure of the data. AVAILABILITY: An executable is freely available for non-commercial users at ftp://ftp.vital-it.ch/tools/dbc454. It is designed to run under MPI on a cluster of 64-bit Linux machines running Red Hat 4.x, or on a multi-core OSX system. CONTACT: dbc454@vital-it.ch or nicolas.guex@isb-sib.ch.