Candida spp. isolated from inpatients, the environment, and health practitioners in the Pediatric Unit at the Universitary Hospital of the Jundiai Medical College, State of Sao Paulo, Brazil

Introduction: This study aimed to isolate and identify Candida spp. from the environment, health practitioners, and patients with the presumptive diagnosis of candidiasis in the Pediatric Unit at the Universitary Hospital of the Jundiai Medical College, to verify the production of enzymes regarded as virulence factors, and to determine how susceptible the isolated samples from patients with candidiasis are to antifungal agents. Methods: Between March and November of 2008 a total of 283 samples were taken randomly from the environment and from the hands of health staff, and samples of all the suspected cases of Candida spp. hospital-acquired infection were collected and selected by the Infection Control Committee. The material was processed and the yeast genus Candida was isolated and identified by physiological, microscopic, and macroscopic attributes. Results: The incidence of Candida spp. in the environment and employees was 19.2%. The most frequent species were C. parapsilosis and C. tropicalis among the workers, C. guilliermondii and C. tropicalis in the air, C. lusitanae on the contact surfaces, and C. tropicalis and C. guilliermondii in the climate control equipment. The college hospital had 320 admissions, of which 13 (4%) presented Candida spp. infections; three of them died, two being victims of a C. tropicalis infection and the remaining one of C. albicans. All the Candida spp. in the isolates evidenced sensitivity to amphotericin B, nystatin, and fluconazole. Conclusions: The increase in the rate of hospital-acquired infections caused by Candida spp. indicates the need to take larger measures regarding recurrent control of the environment.

Pollination services at risk: Bee habitats will decrease owing to climate change in Brazil

Native bees are important providers of pollination services, but there are cumulative evidences of their decline. Global changes such as habitat losses, invasions of exotic species and climate change have been suggested as the main causes of the decline of pollinators. In this study, the influence of climate change on the distribution of 10 species of Brazilian bees was estimated with species distribution modelling. We used Maxent algorithm (maximum entropy) and two different scenarios, an optimistic and a pessimistic, to the years 2050 and 2080. We also evaluated the percentage reduction of species habitat based on the future scenarios of climate change through Geographic Information System (GIS). Results showed that the total area of suitable habitats decreased for all species but one under the different future scenarios. The greatest reductions in habitat area were found for Melipona bicolor bicolor and Melipona scutellaris, which occur predominantly in areas related originally to Atlantic Moist Forest. The species analysed have been reported to be pollinators of some regional crops and the consequence of their decrease for these crops needs further clarification. (C) 2012 Elsevier B.V. All rights reserved.

Identifying the areas to preserve passion fruit pollination service in Brazilian Tropical Savannas under climate change

The aim of this study was to identify future distribution areas and propose actions to preserve passion fruit pollination service under a scenario of future climate change. We used four species of Xylocopa bees that are important for passion fruit pollination in Brazilian Tropical Savannas. We also used the known forage plant species (33 species) that are associated with this same area, since passion fruit flowers provide only nectar for bees and only during their blossoming period. We used species distribution modeling to predict the potential areas of occurrence for each bee and plant based on the current day distribution and a future climate scenario (moderate projections of climate change to 2050). We used a geographic information system to classify the models and to analyze the future areas for both groups of species. The current day distribution map showed that Xylocopa and plant species occurred primarily in the southern and central-eastern areas of the Brazilian Tropical Savannas. In the north, Xylocopa species only occurred in a small area between the states of Maranhão and Piauí while forage plant species were only observed in the northern part of the Tocantins State. However, both future scenarios (bees and plants) showed a shift in distribution, with occurrence predominantly detected in the northern areas of Brazilian Tropical Savannas. Possible conservation areas and the use of appropriate agricultural practices were suggested to ensure the maintenance of the bee/plant focal species.

Desafios atuais da modelagem preditiva de distribuição de espécies

A modelagem preditiva tem sido aplicada para analisar a distribuição geográfica de espécies, a partir de extrapolações das características ambientais dos locais conhecidos de ocorrência. O interesse por esse tipo de modelagem deve-se à necessidade de respostas rápidas e fundamentadas para as ameaças que as espécies têm enfrentado, devido à perda de habitat, invasão de espécies exóticas, mudanças climáticas, entre outros. Este artigo oferece uma visão geral dos avanços recentes no campo da modelagem e visa incentivar a discussão e aplicação desse método, que pode auxiliar tanto na aquisição de conhecimento básico sobre a biologia das espécies, quanto na análise e formulação de políticas para sua conservação.

Natural patterns and anthropogenic disturbance in north Adriatic marine benthic assemblages: descriptive and methodological studies

Marine soft bottom systems show a high variability across multiple spatial and temporal scales. Both natural and anthropogenic sources of disturbance act together in affecting benthic sedimentary characteristics and species distribution. The description of such spatial variability is required to understand the ecological processes behind them. However, in order to have a better estimate of spatial patterns, methods that take into account the complexity of the sedimentary system are required. This PhD thesis aims to give a significant contribution both in improving the methodological approaches to the study of biological variability in soft bottom habitats and in increasing the knowledge of the effect that different process (both natural and anthropogenic) could have on the benthic communities of a large area in the North Adriatic Sea. Beta diversity is a measure of the variability in species composition, and Whittaker’s index has become the most widely used measure of beta-diversity. However, application of the Whittaker index to soft bottom assemblages of the Adriatic Sea highlighted its sensitivity to rare species (species recorded in a single sample). This over-weighting of rare species induces biased estimates of the heterogeneity, thus it becomes difficult to compare assemblages containing a high proportion of rare species. In benthic communities, the unusual large number of rare species is frequently attributed to a combination of sampling errors and insufficient sampling effort. In order to reduce the influence of rare species on the measure of beta diversity, I have developed an alternative index based on simple probabilistic considerations. It turns out that this probability index is an ordinary Michaelis-Menten transformation of Whittaker's index but behaves more favourably when species heterogeneity increases. The suggested index therefore seems appropriate when comparing patterns of complexity in marine benthic assemblages. Although the new index makes an important contribution to the study of biodiversity in sedimentary environment, it remains to be seen which processes, and at what scales, influence benthic patterns. The ability to predict the effects of ecological phenomena on benthic fauna highly depends on both spatial and temporal scales of variation. Once defined, implicitly or explicitly, these scales influence the questions asked, the methodological approaches and the interpretation of results. Problem often arise when representative samples are not taken and results are over-generalized, as can happen when results from small-scale experiments are used for resource planning and management. Such issues, although globally recognized, are far from been resolved in the North Adriatic Sea. This area is potentially affected by both natural (e.g. river inflow, eutrophication) and anthropogenic (e.g. gas extraction, fish-trawling) sources of disturbance. Although few studies in this area aimed at understanding which of these processes mainly affect macrobenthos, these have been conducted at a small spatial scale, as they were designated to examine local changes in benthic communities or particular species. However, in order to better describe all the putative processes occurring in the entire area, a high sampling effort performed at a large spatial scale is required. The sedimentary environment of the western part of the Adriatic Sea was extensively studied in this thesis. I have described, in detail, spatial patterns both in terms of sedimentary characteristics and macrobenthic organisms and have suggested putative processes (natural or of human origin) that might affect the benthic environment of the entire area. In particular I have examined the effect of off shore gas platforms on benthic diversity and tested their effect over a background of natural spatial variability. The results obtained suggest that natural processes in the North Adriatic such as river outflow and euthrophication show an inter-annual variability that might have important consequences on benthic assemblages, affecting for example their spatial pattern moving away from the coast and along a North to South gradient. Depth-related factors, such as food supply, light, temperature and salinity play an important role in explaining large scale benthic spatial variability (i.e., affecting both the abundance patterns and beta diversity). Nonetheless, more locally, effects probably related to an organic enrichment or pollution from Po river input has been observed. All these processes, together with few human-induced sources of variability (e.g. fishing disturbance), have a higher effect on macrofauna distribution than any effect related to the presence of gas platforms. The main effect of gas platforms is restricted mainly to small spatial scales and related to a change in habitat complexity due to a natural dislodgement or structure cleaning of mussels that colonize their legs. The accumulation of mussels on the sediment reasonably affects benthic infauna composition. All the components of the study presented in this thesis highlight the need to carefully consider methodological aspects related to the study of sedimentary habitats. With particular regards to the North Adriatic Sea, a multi-scale analysis along natural and anthopogenic gradients was useful for detecting the influence of all the processes affecting the sedimentary environment. In the future, applying a similar approach may lead to an unambiguous assessment of the state of the benthic community in the North Adriatic Sea. Such assessment may be useful in understanding if any anthropogenic source of disturbance has a negative effect on the marine environment, and if so, planning sustainable strategies for a proper management of the affected area.

Molecular diversity of Primula apennina and phylogeny of Primula Subsection Euauricula

Primula apennina Widmer is endemic to the North Apennines (Italy). ISSR were used to detect the genetic diversity within and among six populations representative of the species distribution range. High levels of genetic diversity were revealed both at population (PPB = 75.92%, HS = 0.204, Hpop = 0.319) and at species level (PPB = 96.95%, HT = 0.242, Hsp = 0.381). Nei gene diversity statistics (15.7%), Shannon diversity index (16.3%) and AMOVA (14%) detected a moderate level of interpopulation diversity. Principal coordinate and bayesian analyses clustered the populations in three major groups along a geographic gradient. The correlation between genetic and geographic distances was positive (Mantel test, r = 0.232). All together, these analyses revealed a weak but significant spatial genetic structure in P. apennina, with gene flow acting as a homogenizing force that prevents a stronger differentiation of populations. Conservation measures are suggested based on the observed pattern of genetic variability. P. apennina belongs to Primula subsect. Euauricula which includes 15 species distributed on the whole Alps and Apennines. A phylogenetic analysis was carried out using AFLP markers in order both to clarify the relationships among the species of subsection Euauricula that remained unresolved in previous works and to make some hypoteses on their evolutive dynamics. NJ, PCO and BAPS analyses strongly confirmed the monophyly of P. subsect. Euauricula and all the species form strongly supported clades. NJ tree topology suggested a simultaneous fragmentations of ancestral species in a large number of isolated populations that survived in refugia along the unglaciated margins of the Alps in response to the Pleistocene climatic oscillations.

Rekonstruktion der eiszeitlichen Verbreitung und Artbildung vier alpiner Primeln durch Artenverbreitungsmodelle und Phylogeographie

Aim: Previous studies revealed that diversification events in the western clade of the alpine Primula sect. Auricula were concentrated in the Quaternary cold periods. This implies that allopatric speciation in isolated glacial refugia was the most common mode of speciation. In the first part of the present dissertation, this hypothesis is further investigated by locating refugial areas of two sister species, Primula marginata & P. latifolia during the last glacial maximum, 21,000 years ago. In the second part, the glacial and postglacial history of P. hirsuta and P. daonensis is investigated. Location: European Alps. Methods: Glacial refugia were located using species distribution models, which are projected to last glacial maximum climate. These refugia are validated with geographic distribution patterns of intra-specific genetic diversity, rarity and variation. Results 1) Speciation: Glacial refugia of the sister taxa Primula marginata and P. latifolia were largely separated, only a small overlapping zone at the southern margin of the former glacier in the Maritime Alps exists. This overlapping zone is too small to indicate sympatric speciation. The largely separated glacial distribution of both species rather confirms our hypothesis of allopatric speciation in isolated glacial refugia. Results 2) Glacial and postglacial history: Surprizingly, the modelled potential refugia of three out of four Primula species are situated within the former ice-shield, except for P. marginata. This indicates that peripheral and central nunataks played an important role for the glacial survival in P. latifolia, P. hirsuta and P. daonensis, while peripheral refugia outside the maximum extend of the glacier were crucial in P. marginata. In P. hirsuta and P. latifolia SDMs allowed to exclude several hypothetical refugial areas that overlap with today’s distribution as potential refugia for the species. In P. marginata, hypothetical refugial areas at the periphery of the former ice-shield that overlap with today’s distribution were confirmed by the models. The results from the SDMs are confirmed by population genetic patterns in three out of four species. P. daonensis represents an exception, where population genetic data contradict the SDMs. Main conclusions: Species distribution models provide species specific scenarios of glacial distribution and postglacial re-colonization, which can be validated using population genetic analyses. This combined approach is useful and helps to understand the complex processes that have lead to the genetic and floristic patterns of biodiversity that is found today in the Alps.

Landscape ecology and genetics of the wolf in Italy

This PhD Thesis includes five main parts on diverse topics. The first two parts deal with the trophic ecology of wolves in Italy consequently to a recent increase of wild ungulates abundance. Data on wolf diet across time highlighted how wild ungulates are important food resource for wolves in Italy. Increasing wolf population, increasing numbers of wild ungulates and decreasing livestock consume are mitigating wolf-man conflicts in Italy in the near future. In the third part, non-invasive genetic sampling techniques were used to obtain genotypes and genders of about 400 wolves. Thus, wolf packs were genetically reconstructed using diverse population genetic and parentage software. Combining the results on pack structure and genetic relatedness with sampling locations, home ranges of wolf packs and dispersal patterns were identified. These results, particularly important for the conservation management of wolves in Italy, illustrated detailed information that can be retrieved from genetic identification of individuals. In the fourth part, wolf locations were combined with environmental information obtained as GIS-layers. Modern species distribution models (niche models) were applied to infer potential wolf distribution and predation risk. From the resulting distribution maps, information pastures with the highest risk of depredation were derived. This is particularly relevant as it allows identifying those areas under danger of carnivore attack on livestock. Finally, in the fifth part, habitat suitability models were combined with landscape genetic analysis. On one side landscape genetic analyses on the Italian wolves provided new information on the dynamics and connectivity of the population and, on the other side, a profound analysis of the effects that habitat suitability methods had on the parameterization of landscape genetic analyses was carried out to contributed significantly to landscape genetic theory.

Counterbalancing effects of competition for resources and facilitation against grazing in alpine snowbed communities

Alpine snowbeds are habitats where the major limiting factors for plant growth are herbivory and a small time window for growth due to late snowmelt. Despite these limitations, snowbed vegetation usually forms a dense carpet of palatable plants due to favourable abiotic conditions for plant growth within the short growing season. These environmental characteristics make snowbeds particularly interesting to study the interplay of facilitation and competition. We hypothesised an interplay between resource competition and facilitation against herbivory. Further, we investigated whether these predicted neighbour effects were species-specific and/or dependent on ontogeny, and whether the balance of positive and negative plant–plant interactions shifted along a snowmelt gradient. We determined the neighbour effects by means of neighbour removal experiments along the snowmelt gradient, and linear mixed model analyses. The results showed that the effects of neighbour removal were weak but generally consistent among species and snowmelt dates, and depended on whether biomass production or survival was considered. Higher total biomass and increased fruiting in removal plots indicated that plants competed for nutrients, water, and light, thereby supporting the hypothesis of prevailing competition for resources in snowbeds. However, the presence of neighbours reduced herbivory and thereby also facilitated survival. For plant growth the facilitative effects against herbivores in snowbeds counterbalanced competition for resources, leading to a weak negative net effect. Overall the neighbour effects were not species-specific and did not change with snowmelt date. Our finding of counterbalancing effects of competition and facilitation within a plant community is of special theoretical value for species distribution models and can explain the success of models that give primary importance to abiotic factors and tend to overlook interrelations between biotic and abiotic effects on plants.

Establishment of Regional Herbarium Leads to More Than 200 New Flora Atlas Records for New York State

Loss of small herbaria is an unfortunate global trend, and initiation of new collections at small academic institutions is an increasingly rare occurrence. In 2006, a new herbarium was established at the State University of New York College at Plattsburgh. The PLAT herbarium has since grown to more than 7,000 specimens, many of them representative of the flora of northeastern New York (especially Clinton County). Previous to 2006, this region was without a recognized herbarium, the nearest in-state collections being more than 150 miles away. Although botanists have previously worked in the region, relatively few plant species were recorded for Clinton County by the New York Flora Atlas – a resource providing species distribution records based on specimens accessioned in herbarium collections. Given the dearth of available distribution data for Clinton County (including the eastern Adirondack Mountains and the western Lake Champlain valley), this project sought to provide records of previously unreported species by comparing NY Flora Atlas maps with current holdings. 203 species will now be added to the NY Flora Atlas for Clinton County, roughly half of those considered exotic. This exercise has amplified the importance of supporting and maintaining small regional herbaria as repositories of valuable biodiversity information. Likewise, this project also highlights the enduring value of training in floristics and taxonomy.

Spatial and temporal variation in population genetic structure of wild Nile tilapia (Oreochromis niloticus) across Africa

Background: Reconstructing the evolutionary history of a species is challenging. It often depends not only on the past biogeographic and climatic events but also the contemporary and ecological factors, such as current connectivity and habitat heterogeneity. In fact, these factors might interact with each other and shape the current species distribution. However, to what extent the current population genetic structure reflects the past and the contemporary factors is largely unknown. Here we investigated spatio-temporal genetic structures of Nile tilapia (Oreochromis niloticus) populations, across their natural distribution in Africa. While its large biogeographic distribution can cause genetic differentiation at the paleo-biogeographic scales, its restricted dispersal capacity might induce a strong genetic structure at micro-geographic scales. Results: Using nine microsatellite loci and 350 samples from ten natural populations, we found the highest genetic differentiation among the three ichthyofaunal provinces and regions (Ethiopian, Nilotic and Sudano-Sahelian) (R(ST) = 0.38 - 0.69). This result suggests the predominant effect of paleo-geographic events at macro-geographic scale. In addition, intermediate divergences were found between rivers and lakes within the regions, presumably reflecting relatively recent interruptions of gene flow between hydrographic basins (R(ST) = 0.24 - 0.32). The lowest differentiations were observed among connected populations within a basin (R(ST) = 0.015 in the Volta basin). Comparison of temporal sample series revealed subtle changes in the gene pools in a few generations (F = 0 - 0.053). The estimated effective population sizes were 23 - 143 and the estimated migration rate was moderate (m similar to 0.094 - 0.097) in the Volta populations. Conclusions: This study revealed clear hierarchical patterns of the population genetic structuring of O. niloticus in Africa. The effects of paleo-geographic and climatic events were predominant at macro-geographic scale, and the significant effect of geographic connectivity was detected at micro-geographic scale. The estimated effective population size, the moderate level of dispersal and the rapid temporal change in genetic composition might reflect a potential effect of life history strategy on population dynamics. This hypothesis deserves further investigation. The dynamic pattern revealed at micro-geographic and temporal scales appears important from a genetic resource management as well as from a biodiversity conservation point of view.

Growth response to climate and drought change along an aridity gradient in the southernmost Pinus nigra relict forests

Tree populations at the rear edge of species distribution are sensitive to climate stress and drought. However, growth responses of these tree populations to those stressors may vary along climatic gradients. To analyze growth responses to climate and drought using dendrochronology in rear-edge Pinus nigra populations located along an aridity gradient. Tree-ring width chronologies were built for the twentieth century and related to monthly climatic variables, a drought index (Standardized Precipitation-Evapotranspiration Index), and two atmospheric circulation patterns (North Atlantic and Western Mediterranean Oscillations). Growth was enhanced by wet and cold previous autumns and warm late winters before tree-ring formation. The influence of the previous year conditions on growth increased during the past century. Growth was significantly related to North Atlantic and Western Mediterranean Oscillations in two out of five sites. The strongest responses of growth to the drought index were observed in the most xeric sites. Dry conditions before tree-ring formation constrain growth in rear-edge P. nigra populations. The comparisons of climate-growth responses along aridity gradients allow characterizing the sensitivity of relict stands to climate warming.

Isolation and characterization of nine polymorphic microsatellite markers for the deep-sea shrimp Nematocarcinus lanceopes (Crustacea: Decapoda: Caridea)

Background: The shrimp Nematocarcinus lanceopes Bate, 1888 is found in the deep sea around Antarctica and sub-Antarctic islands. Previous studies on mitochondrial data and species distribution models provided evidence for a homogenous circum-Antarctic population of N. lanceopes. However, to analyze the fine-scale population genetic structure and to examine influences of abiotic environmental conditions on population composition and genetic diversity, a set of fast evolving nuclear microsatellite markers is required. Findings: We report the isolation and characterization of nine polymorphic microsatellite markers from the Antarctic deep-sea shrimp species Nematocarcinus lanceopes (Crustacea: Decapoda: Caridea). Microsatellite markers were screened in 55 individuals from different locations around the Antarctic continent. All markers were polymorphic with 9 to 25 alleles per locus. The observed heterozygosity ranged from 0.545 to 0.927 and the expected heterozygosity from 0.549 to 0.934. Conclusions: The reported markers provide a novel tool to study genetic structure and diversity in Nematocarcinus lanceopes populations in the Southern Ocean and monitor effects of ongoing climate change in the region on the populations inhabiting these.

Climate and vegetation changes during the Lateglacial and early–middle Holocene at Lake Ledro (southern Alps, Italy)

Adding to the on-going debate regarding vegetation recolonisation (more particularly the timing) in Europe and climate change since the Lateglacial, this study investigates a long sediment core (LL081) from Lake Ledro (652ma.s.l., southern Alps, Italy). Environmental changes were reconstructed using multiproxy analysis (pollen-based vegetation and climate reconstruction, lake levels, magnetic susceptibility and X-ray fluorescence (XRF) measurements) recorded climate and land-use changes during the Lateglacial and early-middle Holocene. The well-dated and high-resolution pollen record of Lake Ledro is compared with vegetation records from the southern and northern Alps to trace the history of tree species distribution. An altitudedependent progressive time delay of the first continuous occurrence of Abies (fir) and of the Larix (larch) development has been observed since the Lateglacial in the southern Alps. This pattern suggests that the mid-altitude Lake Ledro area was not a refuge and that trees originated from lowlands or hilly areas (e.g. Euganean Hills) in northern Italy. Preboreal oscillations (ca. 11 000 cal BP), Boreal oscillations (ca. 10 200, 9300 cal BP) and the 8.2 kyr cold event suggest a centennial-scale climate forcing in the studied area. Picea (spruce) expansion occurred preferentially around 10 200 and 8200 cal BP in the south-eastern Alps, and therefore reflects the long-lasting cumulative effects of successive boreal and the 8.2 kyr cold event. The extension of Abies is contemporaneous with the 8.2 kyr event, but its development in the southern Alps benefits from the wettest interval 8200-7300 cal BP evidenced in high lake levels, flood activity and pollen-based climate reconstructions. Since ca. 7500 cal BP, a weak signal of pollen-based anthropogenic activities suggest weak human impact. The period between ca. 5700 and ca. 4100 cal BP is considered as a transition period to colder and wetter conditions (particularly during summers) that favoured a dense beech (Fagus) forest development which in return caused a distinctive yew (Taxus) decline.We conclude that climate was the dominant factor controlling vegetation changes and erosion processes during the early and middle Holocene (up to ca. 4100 cal BP).

Orbital changes, variation in solar activity and increased anthropogenic activities: controls on the Holocene flood frequency in the Lake Ledro area, Northern Italy

Adding to the on-going debate regarding vegetation recolonisation (more particularly the timing) in Europe and climate change since the Lateglacial, this study investigates a long sediment core (LL081) from Lake Ledro (652ma.s.l., southern Alps, Italy). Environmental changes were reconstructed using multiproxy analysis (pollen-based vegetation and climate reconstruction, lake levels, magnetic susceptibility and X-ray fluorescence (XRF) measurements) recorded climate and land-use changes during the Lateglacial and early-middle Holocene. The well-dated and high-resolution pollen record of Lake Ledro is compared with vegetation records from the southern and northern Alps to trace the history of tree species distribution. An altitudedependent progressive time delay of the first continuous occurrence of Abies (fir) and of the Larix (larch) development has been observed since the Lateglacial in the southern Alps. This pattern suggests that the mid-altitude Lake Ledro area was not a refuge and that trees originated from lowlands or hilly areas (e.g. Euganean Hills) in northern Italy. Preboreal oscillations (ca. 11 000 cal BP), Boreal oscillations (ca. 10 200, 9300 cal BP) and the 8.2 kyr cold event suggest a centennial-scale climate forcing in the studied area. Picea (spruce) expansion occurred preferentially around 10 200 and 8200 cal BP in the south-eastern Alps, and therefore reflects the long-lasting cumulative effects of successive boreal and the 8.2 kyr cold event. The extension of Abies is contemporaneous with the 8.2 kyr event, but its development in the southern Alps benefits from the wettest interval 8200-7300 cal BP evidenced in high lake levels, flood activity and pollen-based climate reconstructions. Since ca. 7500 cal BP, a weak signal of pollen-based anthropogenic activities suggest weak human impact. The period between ca. 5700 and ca. 4100 cal BP is considered as a transition period to colder and wetter conditions (particularly during summers) that favoured a dense beech (Fagus) forest development which in return caused a distinctive yew (Taxus) decline.We conclude that climate was the dominant factor controlling vegetation changes and erosion processes during the early and middle Holocene (up to ca. 4100 cal BP).