Evolutionary Origins and Patterns of Diversification in Animal Brood Parasitism.

AbstractBrood parasitism involves the exploitation of host parental care rather than the extraction of resources directly from hosts. We identify defining characteristics of this strategy and consider its position along continua with adjacent behaviors but focus on canonical brood parasites, where parasitism is obligate and hosts are noneusocial (thereby distinguishing from social parasitism). A systematic literature survey revealed 59 independently derived brood parasitic lineages with most origins (49) in insects, particularly among bees and wasps, and other origins in birds (seven) and fish (three). Insects account for more than 98% of brood parasitic species, with much of that diversity reflecting ancient (≥100-million-year-old) brood parasitic lineages. Brood parasites usually, but not always, evolve from forms that show parental care. In insects, brood parasitism often first evolves through exploitation of a closely related species, following Emery's rule, but this is less typical in birds, which we discuss. We conducted lineage-level comparisons between brood parasitic clades and their sister groups, finding mixed results but an overall neutral to negative effect of brood parasitism on species richness and diversification. Our review of brood parasites reveals many unanswered questions requiring new research, including further modeling of the coevolutionary dynamics of brood parasites and their hosts.

Phylogenetic relationships and the evolution of host preferences in the largest clade of brood parasitic bees (Apidae: Nomadinae).

Brood parasites (also known as cleptoparasites) represent a substantial fraction of global bee diversity. Rather than constructing their own nests, these species instead invade those of host bees to lay their eggs. Larvae then hatch and consume the food provisions intended for the host's offspring. While this life history strategy has evolved numerous times across the phylogeny of bees, the oldest and most speciose parasitic clade is the subfamily Nomadinae (Apidae). However, the phylogenetic relationships among brood parasitic apids both within and outside the Nomadinae have not been fully resolved. Here, we present new findings on the phylogeny of this diverse group of brood parasites based on ultraconserved element (UCE) sequence data and extensive taxon sampling with 114 nomadine species representing all tribes. We suggest a broader definition of the subfamily Nomadinae to describe a clade that includes almost all parasitic members of the family Apidae. The tribe Melectini forms the sister group to all other Nomadinae, while the remainder of the subfamily is composed of two sister clades: a "nomadine line" representing the former Nomadinae sensu stricto, and an "ericrocidine line" that unites several mostly Neotropical lineages. We find the tribe Osirini Handlirsch to be polyphyletic, and divide it into three lineages, including the newly described Parepeolini trib. nov. In addition to our taxonomic findings, we use our phylogeny to explore the evolution of different modes of parasitism, detecting two independent transitions from closed-cell to open-cell parasitism. Finally, we examine how nomadine host-parasite associations have evolved over time. In support of Emery's rule, which suggests close relationships between hosts and parasites, we confirm that the earliest nomadines were parasites of their close free-living relatives within the family Apidae, but that over time their host range broadened to include more distantly related hosts spanning the diversity of bees. This expanded breadth of host taxa may also be associated with the transition to open-cell parasitism.

Gene Tree Estimation Error with Ultraconserved Elements: An Empirical Study on Pseudapis Bees.

Summarizing individual gene trees to species phylogenies using two-step coalescent methods is now a standard strategy in the field of phylogenomics. However, practical implementations of summary methods suffer from gene tree estimation error, which is caused by various biological and analytical factors. Greatly understudied is the choice of gene tree inference method and downstream effects on species tree estimation for empirical data sets. To better understand the impact of this method choice on gene and species tree accuracy, we compare gene trees estimated through four widely used programs under different model-selection criteria: PhyloBayes, MrBayes, IQ-Tree, and RAxML. We study their performance in the phylogenomic framework of $>$800 ultraconserved elements from the bee subfamily Nomiinae (Halictidae). Our taxon sampling focuses on the genus Pseudapis, a distinct lineage with diverse morphological features, but contentious morphology-based taxonomic classifications and no molecular phylogenetic guidance. We approximate topological accuracy of gene trees by assessing their ability to recover two uncontroversial, monophyletic groups, and compare branch lengths of individual trees using the stemminess metric (the relative length of internal branches). We further examine different strategies of removing uninformative loci and the collapsing of weakly supported nodes into polytomies. We then summarize gene trees with ASTRAL and compare resulting species phylogenies, including comparisons to concatenation-based estimates. Gene trees obtained with the reversible jump model search in MrBayes were most concordant on average and all Bayesian methods yielded gene trees with better stemminess values. The only gene tree estimation approach whose ASTRAL summary trees consistently produced the most likely correct topology, however, was IQ-Tree with automated model designation (ModelFinder program). We discuss these findings and provide practical advice on gene tree estimation for summary methods. Lastly, we establish the first phylogeny-informed classification for Pseudapis s. l. and map the distribution of distinct morphological features of the group. [ASTRAL; Bees; concordance; gene tree estimation error; IQ-Tree; MrBayes, Nomiinae; PhyloBayes; RAxML; phylogenomics; stemminess].

Wild insect diversity increases inter-annual stability in global crop pollinator communities.

While an increasing number of studies indicate that the range, diversity and abundance of many wild pollinators has declined, the global area of pollinator-dependent crops has significantly increased over the last few decades. Crop pollination studies to date have mainly focused on either identifying different guilds pollinating various crops, or on factors driving spatial changes and turnover observed in these communities. The mechanisms driving temporal stability for ecosystem functioning and services, however, remain poorly understood. Our study quantifies temporal variability observed in crop pollinators in 21 different crops across multiple years at a global scale. Using data from 43 studies from six continents, we show that (i) higher pollinator diversity confers greater inter-annual stability in pollinator communities, (ii) temporal variation observed in pollinator abundance is primarily driven by the three-most dominant species, and (iii) crops in tropical regions demonstrate higher inter-annual variability in pollinator species richness than crops in temperate regions. We highlight the importance of recognizing wild pollinator diversity in agricultural landscapes to stabilize pollinator persistence across years to protect both biodiversity and crop pollination services. Short-term agricultural management practices aimed at dominant species for stabilizing pollination services need to be considered alongside longer term conservation goals focussed on maintaining and facilitating biodiversity to confer ecological stability.

Macro-conformality of coatings deposited using high-speed spatial plasma-enhanced atomic layer deposition.

Rotational spatial plasma-enhanced atomic layer deposition has been used to deposit thin films on half-sphere lenses. Non-uniformity of less than ±1% is demonstrated for Nb2O5 deposited at 1.4 Å/s and for SiO2 deposited at 6 Å/s.

Omnivory in Bees: Elevated Trophic Positions among All Major Bee Families.

As pollen and nectar foragers, bees have long been considered strictly herbivorous. Their pollen provisions, however, are host to abundant microbial communities, which feed on the pollen before and/or while it is consumed by bee larvae. In the process, microbes convert pollen into a complex of plant and microbial components. Since microbes are analogous to metazoan consumers within trophic hierarchies, the pollen-eating microbes are, functionally, herbivores. When bee larvae consume a microbe-rich pollen complex, they ingest proteins from plant and microbial sources and thus should register as omnivores on the trophic "ladder." We tested this hypothesis by examining the isotopic compositions of amino acids extracted from native bees collected in North America over multiple years. We measured bee trophic position across the six major bee families. Our findings indicate that bee trophic identity was consistently and significantly higher than that of strict herbivores, providing the first evidence that omnivory is ubiquitous among bee fauna. Such omnivory suggests that pollen-borne microbes represent an important protein source for larval bees, which introduces new questions as to the link between floral fungicide residues and bee development.

Viral transmission in honey bees and native bees, supported by a global black queen cell virus phylogeny.

In recent decades, we have realized that honey bee viruses are not, in fact, exclusive to honey bees. The potential impact of Apis-affiliated viruses on native pollinators is prompting concern. Our research addresses the issue of virus crossover between honey bees and native bees foraging in the same localities. We measured the presence of black queen cell virus (BQCV), deformed wing virus (DWV) and sacbrood virus (SBV) in managed Apis mellifera (honey bees) and native Andrena spp. (subgenus Melandrena) bee populations in five commercial orchards. We identified viral presence across sites and bees and related these data to measures of bee community diversity. All viruses were found in both managed and native bees, and BQCV was the most common virus in each. To establish evidence for viral crossover between taxa, we undertook an additional examination of BQCV where 74 samples were sequenced and placed in a global phylogenic framework of hundreds of BQCV strains. We demonstrate pathogen sharing across managed honey bees and distantly related wild bees. This phylogenetic analysis contributes to growing evidence for host switching and places local incidence patterns in a worldwide context, revealing multispecies viral transmission.

Combining transcriptomes and ultraconserved elements to illuminate the phylogeny of Apidae.

Two increasingly popular approaches to reconstruct the Tree of Life involve whole transcriptome sequencing and the target capture of ultraconserved elements (UCEs). Both methods can be used to generate large, multigene datasets for analysis of phylogenetic relationships in non-model organisms. While targeted exon sequencing across divergent lineages is now a standard method, it is still not clear if UCE data can be readily combined with published transcriptomes. In this study, we evaluate the combination of UCEs and transcriptomes in a single analysis using genome-, transcriptome-, and UCE data for 79 bees in the largest and most biologically diverse bee family, Apidae. Using existing tools, we first developed a workflow to assemble phylogenomic data from different sources and produced two large nucleotide matrices of combined data. We then reconstructed the phylogeny of the Apidae using concatenation- and coalescent-based methods, and critically evaluated the resulting phylogenies in the context of previously published genetic, genomic, and morphological data sets. Our estimated phylogenetic trees are robustly supported and largely congruent with previous molecular hypotheses, from deep nodes to shallow species-level phylogenies. Moreover, the combined approach allows us to resolve controversial nodes of the apid Tree of Life, by clarifying the relationships among the genera of orchid bees (Euglossini) and the monophyly of the Centridini. Additionally, we present novel phylogenetic evidence supporting the monophyly of the diverse clade of cleptoparasitic Apidae and the placement of two enigmatic, oil-collecting genera (Ctenoplectra and Tetrapedia). Lastly, we propose a revised classification of the family Apidae that reflects our improved understanding of apid higher-level relationships.

Landscape context shifts the balance of costs and benefits from wildflower borders on multiple ecosystem services.

In the face of global biodiversity declines driven by agricultural intensification, local diversification practices are broadly promoted to support farmland biodiversity and multiple ecosystem services. The creation of flower-rich habitats on farmland has been subsidized in both the USA and EU to support biodiversity and promote delivery of ecosystem services. Yet, theory suggests that the landscape context in which local diversification strategies are implemented will influence their success. However, few studies have empirically evaluated this theory or assessed the ability to support multiple ecosystem services simultaneously. Here, we evaluate the impact of creating flower-rich habitats in field margins on pollination, pest control, and crop yield over 3 years using a paired design across a landscape gradient. We find general positive effects of natural habitat cover on fruit weight and that flowering borders increase yields by promoting bee visitation to adjacent crops only in landscapes with intermediate natural habitat cover. Flowering borders had little impact on biological control regardless of landscape context. Thus, knowledge of landscape context can be used to target wildflower border placement in areas where they will have the greatest likelihood for success and least potential for increasing pest populations or yield loss in nearby crops.

Landscape simplification reduces classical biological control and crop yield.

Agricultural intensification resulting in the simplification of agricultural landscapes is known to negatively impact the delivery of key ecosystem services such as the biological control of crop pests. Both conservation and classical biological control may be influenced by the landscape context in which they are deployed; yet studies examining the role of landscape structure in the establishment and success of introduced natural enemies and their interactions with native communities are lacking. In this study, we investigated the relationship between landscape simplification, classical and conservation biological control services and importantly, the outcome of these interactions for crop yield. We showed that agricultural simplification at the landscape scale is associated with an overall reduction in parasitism rates of crop pests. Additionally, only introduced parasitoids were identified, and no native parasitoids were found in crop habitat, irrespective of agricultural landscape simplification. Pest densities in the crop were lower in landscapes with greater proportions of semi-natural habitats. Furthermore, farms with less semi-natural cover in the landscape and consequently, higher pest numbers, had lower yields than farms in less agriculturally dominated landscapes. Our study demonstrates the importance of landscape scale agricultural simplification in mediating the success of biological control programs and highlights the potential risks to native natural enemies in classical biological control programs against native insects. Our results represent an important contribution to an understanding of the landscape-mediated impacts on crop yield that will be essential to implementing effective policies that simultaneously conserve biodiversity and ecosystem services.

Pollinivory and the diversification dynamics of bees.

Pollinivory-the consumption of pollen rather than arthropod prey-is a defining feature of bees (Anthophila; the flower lovers). In virtually all bee species, larvae consume a diet composed of pollen mixed with nectar or floral oils. Bees arose from within a group of solitary, carnivorous, apoid wasps in the Early to Mid-Cretaceous, coincident with the rapid rise of flowering plants. It is assumed that the switch from carnivory to pollen-feeding was a key innovation that led to the rapid diversification of bees, but this has never been examined empirically. Here, we explore the hypothesis that pollinivory led to the increased diversification of bees. In contrast to common perception, we find that the switch to pollen-feeding per se does not explain their extensive diversification. Rather, our results indicate that pollinivory was a necessary but not sufficient condition for diversification, and that other complementary innovations, such as a broadening of host-plant diet, allowed the diversification of the major bee lineages. Our results have broad implications for understanding tempo and mode of bee diversification dynamics in light of their floral resources.

Long-legged bees make adaptive leaps: linking adaptation to coevolution in a plant-pollinator network.

Adaptation is evolution in response to natural selection. Hence, an adaptation is expected to originate simultaneously with the acquisition of a particular selective environment. Here we test whether long legs evolve in oil-collecting Rediviva bees when they come under selection by long-spurred, oil-secreting flowers. To quantify the selective environment, we drew a large network of the interactions between Rediviva species and oil-secreting plant species. The selective environment of each bee species was summarized as the average spur length of the interacting plant species weighted by interaction frequency. Using phylogenetically independent contrasts, we calculated divergence in selective environment and evolutionary divergence in leg length between sister species (and sister clades) of Rediviva We found that change in the selective environment explained 80% of evolutionary change in leg length, with change in body size contributing an additional 6% of uniquely explained variance. The result is one of four proposed steps in testing for plant-pollinator coevolution.

The impact of GC bias on phylogenetic accuracy using targeted enrichment phylogenomic data.

The field of sequence based phylogenetic analyses is currently being transformed by novel hybrid-based targeted enrichment methods, such as the use of ultraconserved elements (UCEs). Rather than analyzing relationships among organisms using a small number of genes, these methods now allow us to evaluate relationships with many hundreds to thousands of individual gene loci. However, the inclusion of thousands of loci does not necessarily overcome the long-standing challenge of incongruence among phylogenetic trees derived from different genes or gene regions. One factor that impacts the level of incongruence in phylogenomic data sets is the level of GC bias. GC rich gene regions are prone to higher recombination rates than AT rich regions, driven by a process referred to as "GC biased gene conversion". As a result, high GC content can be negatively associated with phylogenetic accuracy, but the extent to which this impacts incongruence among UCEs is currently unstudied. We investigated the impact of GC content on phylogeny reconstruction using in silico captured UCE data for the corbiculate bees (Hymenoptera: Apidae). The phylogeny of this group has been the subject of extensive study, and incongruence among gene trees is thought to be a source of phylogenetic error. We conducted coalescent- and concatenation-based analyses of 810 individual gene loci from all 13 currently available bee genomes, including 8 corbiculate taxa. Both coalescent- and concatenation-based methods converged on a single topology for the corbiculate tribes. In contrast to concatenation, the coalescent-based methods revealed significant topological conflict at nodes involving the orchid bees (Euglossini) and honeybees (Apini). Partitioning the loci by GC content reveals decreasing support for the inferred topology with increasing GC bias. Based on the results of this study, we report the first evidence that GC biased gene conversion may contribute to topological incongruence in studies based on ultraconserved elements.

Playing with extremes: Origins and evolution of exaggerated female forelegs in South African Rediviva bees.

Despite close ecological interactions between plants and their pollinators, only some highly specialised pollinators adapt to a specific host plant trait by evolving a bizarre morphology. Here we investigated the evolution of extremely elongated forelegs in females of the South African bee genus Rediviva (Hymenoptera: Melittidae), in which long forelegs are hypothesised to be an adaptation for collecting oils from the extended spurs of their Diascia host flowers. We first reconstructed the phylogeny of the genus Rediviva using seven genes and inferred an origin of Rediviva at around 29MYA (95% HPD=19.2-40.5), concurrent with the origin and radiation of the Succulent Karoo flora. The common ancestor of Rediviva was inferred to be a short-legged species that did not visit Diascia. Interestingly, all our analyses strongly supported at least two independent origins of long legs within Rediviva. Leg length was not correlated with any variable we tested (ecological specialisation, Diascia visitation, geographic distribution, pilosity type) but seems to have evolved very rapidly. Overall, our results indicate that foreleg length is an evolutionary highly labile, rapidly evolving trait that might enable Rediviva bees to respond quickly to changing floral resource availability.

A global synthesis of the effects of diversified farming systems on arthropod diversity within fields and across agricultural landscapes.

Agricultural intensification is a leading cause of global biodiversity loss, which can reduce the provisioning of ecosystem services in managed ecosystems. Organic farming and plant diversification are farm management schemes that may mitigate potential ecological harm by increasing species richness and boosting related ecosystem services to agroecosystems. What remains unclear is the extent to which farm management schemes affect biodiversity components other than species richness, and whether impacts differ across spatial scales and landscape contexts. Using a global metadataset, we quantified the effects of organic farming and plant diversification on abundance, local diversity (communities within fields), and regional diversity (communities across fields) of arthropod pollinators, predators, herbivores, and detritivores. Both organic farming and higher in-field plant diversity enhanced arthropod abundance, particularly for rare taxa. This resulted in increased richness but decreased evenness. While these responses were stronger at local relative to regional scales, richness and abundance increased at both scales, and richness on farms embedded in complex relative to simple landscapes. Overall, both organic farming and in-field plant diversification exerted the strongest effects on pollinators and predators, suggesting these management schemes can facilitate ecosystem service providers without augmenting herbivore (pest) populations. Our results suggest that organic farming and plant diversification promote diverse arthropod metacommunities that may provide temporal and spatial stability of ecosystem service provisioning. Conserving diverse plant and arthropod communities in farming systems therefore requires sustainable practices that operate both within fields and across landscapes.

Crop domestication facilitated rapid geographical expansion of a specialist pollinator, the squash bee Peponapis pruinosa.

Squash was first domesticated in Mexico and is now found throughout North America (NA) along with Peponapis pruinosa, a pollen specialist bee species of the squash genus Cucurbita The origin and spread of squash cultivation is well-studied archaeologically and phylogenetically; however, no study has documented how cultivation of this or any other crop has influenced species in mutualistic interactions. We used molecular markers to reconstruct the demographic range expansion and colonization routes of P. pruinosa from its native range into temperate NA. Populations east of the Rocky Mountains expanded from the wild host plant's range in Mexico and were established by a series of founder events. Eastern North America was most likely colonized from squash bee populations in the present-day continental Midwest USA and not from routes that followed the Gulf and Atlantic coasts from Mexico. Populations of P. pruinosa west of the Rockies spread north from the warm deserts much more recently, showing two genetically differentiated populations with no admixture: one in California and the other one in eastern Great Basin. These bees have repeatedly endured severe bottlenecks as they colonized NA, following human spread of their Cucurbita pollen hosts during the Holocene.

Phylogenetic systematics and a revised generic classification of anthidiine bees (Hymenoptera: Megachilidae).

The bee tribe Anthidiini (Hymenoptera: Megachilidae) is a large, cosmopolitan group of solitary bees that exhibit intriguing nesting behavior. We present the first molecular-based phylogenetic analysis of relationships within Anthidiini using model-based methods and a large, multi-locus dataset (five nuclear genes, 5081 base pairs), as well as a combined analysis using our molecular dataset in conjunction with a previously published morphological matrix. We discuss the evolution of nesting behavior in Anthidiini and the relationship between nesting material and female mandibular morphology. Following an examination of the morphological characters historically used to recognize anthidiine genera, we recommend the use of a molecular-based phylogenetic backbone to define taxonomic groups prior to the assignment of diagnostic morphological characters for these groups. Finally, our results reveal the paraphyly of numerous genera and have significant consequences for anthidiine classification. In order to promote a classification system based on stable, monophyletic clades, we hereby make the following changes to Michener's (2007) classification: The subgenera Afranthidium (Zosteranthidium) Michener and Griswold, 1994, Afranthidium (Branthidium) Pasteels, 1969 and Afranthidium (Immanthidium) Pasteels, 1969 are moved into the genus Pseudoanthidium, thus forming the new combinations Pseudoanthidium (Zosteranthidium), Pseudoanthidium (Branthidium), and Pseudoanthidium (Immanthidium). The genus Neanthidium Pasteels, 1969 is also moved into the genus Pseudoanthidium, thus forming the new combination Pseudoanthidium (Neanthidium). Based on morphological characters shared with our new definition of the genus Pseudoanthidium, the subgenus Afranthidium (Mesanthidiellum) Pasteels, 1969 and the genus Gnathanthidium Pasteels, 1969 are also moved into the genus Pseudoanthidium, thus forming the new combinations Pseudoanthidium (Mesanthidiellum) and Pseudoanthidium (Gnathanthidium).

Historical changes in northeastern US bee pollinators related to shared ecological traits.

Pollinators such as bees are essential to the functioning of terrestrial ecosystems. However, despite concerns about a global pollinator crisis, long-term data on the status of bee species are limited. We present a long-term study of relative rates of change for an entire regional bee fauna in the northeastern United States, based on >30,000 museum records representing 438 species. Over a 140-y period, aggregate native species richness weakly decreased, but richness declines were significant only for the genus Bombus. Of 187 native species analyzed individually, only three declined steeply, all of these in the genus Bombus. However, there were large shifts in community composition, as indicated by 56% of species showing significant changes in relative abundance over time. Traits associated with a declining relative abundance include small dietary and phenological breadth and large body size. In addition, species with lower latitudinal range boundaries are increasing in relative abundance, a finding that may represent a response to climate change. We show that despite marked increases in human population density and large changes in anthropogenic land use, aggregate native species richness declines were modest outside of the genus Bombus. At the same time, we find that certain ecological traits are associated with declines in relative abundance. These results should help target conservation efforts focused on maintaining native bee abundance and diversity and therefore the important ecosystems services that they provide.

Negative effects of pesticides on wild bee communities can be buffered by landscape context.

Wild bee communities provide underappreciated but critical agricultural pollination services. Given predicted global shortages in pollination services, managing agroecosystems to support thriving wild bee communities is, therefore, central to ensuring sustainable food production. Benefits of natural (including semi-natural) habitat for wild bee abundance and diversity on farms are well documented. By contrast, few studies have examined toxicity of pesticides on wild bees, let alone effects of farm-level pesticide exposure on entire bee communities. Whether beneficial natural areas could mediate effects of harmful pesticides on wild bees is also unknown. Here, we assess the effect of conventional pesticide use on the wild bee community visiting apple (Malus domestica) within a gradient of percentage natural area in the landscape. Wild bee community abundance and species richness decreased linearly with increasing pesticide use in orchards one year after application; however, pesticide effects on wild bees were buffered by increasing proportion of natural habitat in the surrounding landscape. A significant contribution of fungicides to observed pesticide effects suggests deleterious properties of a class of pesticides that was, until recently, considered benign to bees. Our results demonstrate extended benefits of natural areas for wild pollinators and highlight the importance of considering the landscape context when weighing up the costs of pest management on crop pollination services.

Climate, physiological tolerance and sex-biased dispersal shape genetic structure of Neotropical orchid bees.

Understanding the impact of past climatic events on the demographic history of extant species is critical for predicting species' responses to future climate change. Palaeoclimatic instability is a major mechanism of lineage diversification in taxa with low dispersal and small geographical ranges in tropical ecosystems. However, the impact of these climatic events remains questionable for the diversification of species with high levels of gene flow and large geographical distributions. In this study, we investigate the impact of Pleistocene climate change on three Neotropical orchid bee species (Eulaema bombiformis, E. meriana and E. cingulata) with transcontinental distributions and different physiological tolerances. We first generated ecological niche models to identify species-specific climatically stable areas during Pleistocene climatic oscillations. Using a combination of mitochondrial and nuclear markers, we inferred calibrated phylogenies and estimated historical demographic parameters to reconstruct the phylogeographical history of each species. Our results indicate species with narrower physiological tolerance experienced less suitable habitat during glaciations and currently exhibit strong population structure in the mitochondrial genome. However, nuclear markers with low and high mutation rates show lack of association with geography. These results combined with lower migration rate estimates from the mitochondrial than the nuclear genome suggest male-biased dispersal. We conclude that despite large effective population sizes and capacity for long-distance dispersal, climatic instability is an important mechanism of maternal lineage diversification in orchid bees. Thus, these Neotropical pollinators are susceptible to disruption of genetic connectivity in the event of large-scale climatic changes.