Genetic diversity and population divergences of an indigenous tree (Coffea mauritiana) in Reunion Island: role of climatic and geographical factors

Climatic factors, but not geographic distance, promote genetic structure and differentiation of Cleistogenes squarrosa (Trin.) Keng populations

Climate can shape plant genetic diversity and genetic structure, and genetic diversity and genetic structure can reflect the adaptation of plants to climate change. We used rbcl and trnL-trnF sequences to analyze the genetic diversity and genetic structure of C. squarrosa under the influence of different environmental factors in Inner Mongolia grassland. The results showed that the genetic diversity of this species was low. (The trnL-trnF sequences have higher genetic diversity than rbcl sequences.) C. squarrosa had low genetic diversity compared to other prairie plants, but had a more pronounced genetic structure. The haplotype network diagram of the combined sequences could be divided into two categories, and the results of the NJ, MP, and ML trees also showed that the haplotypes were divided into two branches. The results of genetic structure analysis showed that that the populations located in the desert steppe fall into exactly one cluster, and the populations located in the typical steppe fall into exactly another cluster. The neutrality tests were all negative and the mismatch distribution also showed a single peak across the population, suggesting that C. squarrosa had undergone population expansion and was well adapted to the local environment. The results of the mantel test showed that climate had a greater influence on the genetic distance of C. squarrosa, with annual precipitation having a higher influence than mean annual temperature. This study provided basic genetic information on the genetic structure of C. squarrosa and contributes to the study of genetic adaptation mechanisms in grassland plants.

Variations in the root mycobiome and mycorrhizal fungi between different types of Vanilla forest farms on Réunion Island

The mycorrhizal fungi of cultivated Vanilla spp. have mainly been studied in America, while a recent study has investigated them on Réunion Island (Indian Ocean). However, there are many different types of cultivation on Réunion, from shade-house crops to forest farms of endemic or exotic trees. Here we fill a gap in the study of the root mycobiome of Vanilla by sampling vines in forest plantations on recent lava flows in the southeast of Réunion. Specifically, we aimed to characterize the fungal communities between terrestrial and epiphytic roots, between forest farms that differ mainly in the species of trees, and between Vanilla roots and ECM-like roots of nearby trees. By sequencing fungal ITS2, we showed that the Vanilla root mycobiome is diverse and differed between the root types and forest farms. Epiphytic and terrestrial roots host endophytic fungi, while a putative rust with visible urediniospores was abundant in terrestrial roots mainly. Other pathogens were detected in epiphytic roots (Colletotrichum) with no sign of disease. Following sequencing and electron microscopy, Tulasnellaceae, characterized by imperforate parenthesomes and cell wall expansion with an amorphous matrix, were shown to be the main mycorrhizal fungi in both vanilla root types. Interestingly, the dominant Tulasnellaceae OTU was found in ECM-type roots of trees belonging to the ectomycorrhizal family Sapotaceae. Further observations are needed to confirm the ectomycorrhizal association of endemic trees with Tulasnella. Moreover, labeling experiments will be instrumental in investigating the transfer of nutrients between the trees and the Vanilla through the network of mycorrhizal associations in the soil.

Quantifying climatic suitability for tourism in Southwest Indian Ocean Tropical Islands: Applying the Holiday Climate Index to Réunion Island

Climate is a key resource for tourists and tourism providers. Varied approaches to quantifying climate resources for tourism have been developed in the last 40 years, computing indices from a range of meteorological variables to measure the comparative climatic suitability of different destinations through time. This study provides the first application of a tourism climate index in the tropical southwest Indian Ocean, applying the recently developed Holiday Climate Index (HCI) for Réunion Island. The suitability of this index is evaluated for the case of this French department, with a particular focus on air conditioning availability in tourism accommodation establishments as this index excludes night-time thermal comfort. Both iterations of the HCI (HCIBeach and HCIUrban) are computed with meteorological data from Roland Garros Airport for the period 1991-2020, exploring monthly, annual, and seasonal climatic suitability. Mean monthly HCI scores reveal considerable seasonality in climatic suitability for tourism on the island with scores ranging from 89.3 ('excellent') to 36.9 ('marginal') for the HCIBeach and 85.0 ('excellent') to 27.5 ('unacceptable') for the HCIUrban, with more favourable scores calculated for July and August, displaying a clear austral winter peak seasonal classification. Over the 30-year period, there is no statically significant change in mean annual climatic suitability, and at a monthly scale, only one month of the year for each index displays statistically significant trends. These results are important in informing tourism strategies for the island to maximise visitor satisfaction through targeting advertising more deliberately for peak touristic climate suitability during the winter months.

Relationships within Bolbitis sinensis Species Complex Using RAD Sequencing

Species identification and phylogenetic relationship clarification are fundamental goals in species delimitation. However, these tasks pose challenges when based on morphologies, geographic distribution, and genomic data. Previously, two species of the fern genus Bolbitis, B. × multipinna and B. longiaurita were described based on morphological traits; they are phylogenetically intertwined with B. sinensis and fail to form monophyletic groups. To address the unclear phylogenetic relationships within the B. sinensis species complex, RAD sequencing was performed on 65 individuals from five populations. Our integrated analysis of phylogenetic trees, neighbor nets, and genetic structures indicate that the B. sinensis species complex should not be considered as separate species. Moreover, our findings reveal differences in the degree of genetic differentiation among the five populations, ranging from low to moderate, which might be influenced by geographical distance and gene flow. The Fst values also confirmed that genetic differentiation intensifies with increasing geographic distance. Collectively, this study clarifies the complex phylogenetic relationships within the B. sinensis species complex, elucidates the genetic diversity and differentiation across the studied populations, and offers valuable genetic insights that contribute to the broader study of evolutionary relationships and population genetics within the Bolbitis species.

Rapid diversification of grey mangroves (Avicennia marina) driven by geographic isolation and extreme environmental conditions in the Arabian Peninsula

Biological systems occurring in ecologically heterogeneous and spatially discontinuous habitats provide an ideal opportunity to investigate the relative roles of neutral and selective factors in driving lineage diversification. The grey mangroves (Avicennia marina) of Arabia occur at the northern edge of the species' range and are subject to variable, often extreme, environmental conditions, as well as historic large fluctuations in habitat availability and connectivity resulting from Quaternary glacial cycles. Here, we analyse fully sequenced genomes sampled from 19 locations across the Red Sea, the Arabian Sea and the Persian/Arabian Gulf (PAG) to reconstruct the evolutionary history of the species in the region and to identify adaptive mechanisms of lineage diversification. Population structure and phylogenetic analyses revealed marked genetic structure correlating with geographic distance and highly supported clades among and within the seas surrounding the Arabian Peninsula. Demographic modelling showed times of divergence consistent with recent periods of geographic isolation and low marine connectivity during glaciations, suggesting the presence of (cryptic) glacial refugia in the Red Sea and the PAG. Significant migration was detected within the Red Sea and the PAG, and across the Strait of Hormuz to the Arabian Sea, suggesting gene flow upon secondary contact among populations. Genetic-environment association analyses revealed high levels of adaptive divergence and detected signs of multi-loci local adaptation driven by temperature extremes and hypersalinity. These results support a process of rapid diversification resulting from the combined effects of historical factors and ecological selection and reveal mangrove peripheral environments as relevant drivers of lineage diversity.

Revisiting the pericarp as a barrier restricting water entry/loss from cotyledons and embryonic axis of temperate desiccation-sensitive Quercus acorns

Fully mature acorns of Quercus variabilis, Q. aliena, Q. mongolica, and Q. glandulifera are desiccation-sensitive. X-ray computer tomography showed that cotyledons shrink during drying, but embryos are protected. Information available on recalcitrant acorns of tropical and sub-tropical species of Quercus suggests that an impermeable pericarp, which limits the entry and loss of water only through the hilum (scar), is the underlying mechanism that prevents drying of the embryo axis following dispersal until the germination season. However, there is a lack of consensus supporting this proposition across species, and it is not well understood if such mechanisms occur in temperate Quercus species. This study investigated the significance of the acorn pericarp for temperate oak species and presents an ecological framework based on the post-dispersal climatic conditions. Using Quercus variabilis, Q. aliena, Q. mongolica, and Q. glandulifera acorns, the relationship between moisture content (MC) and germination was established, and X-ray computed tomography (X-ray CT) was used to understand the internal structural changes of cotyledons and embryonic axis occurring during desiccation. Water entry and exit routes through the scar, pericarp and apex were determined by imbibition and drying experiments. Climatic data and acorn morphological characteristics and germination were subjected to a principal component analysis (PCA). Freshly dispersed acorns of all species had a moisture content (MC) above 35% fresh weight (FW) basis, but drying to 15-10% MC resulted in complete loss of viability, implying recalcitrance behaviour. X-ray CT images suggested that the pericarp offers some protection to cotyledons and embryonic axis during desiccation, but it is contingent on MC. Extensive drying to a low MC with the scar and apex covered with vaseline resulted in internal tissues shrinkage, corresponding with viability loss. Water could enter or exit through the pericarp, albeit at a much slower rate than through the scar. A combination of factors including acorn anatomy, moisture content at the time of dispersal, microhabitat, the position of acorns in the soil prevent embryo desiccation below the critical MC and thus promotes survival of acorns on/in the soil during winter in temperate regions. Pericarp anatomy, to some extent, prevents excessive drying of the embryonic axis by slowing water movement, but prolonged drying or predatory pressure could result in pericarp cracks, favouring the absorption of water during sporadic rain. In the latter case, the survival of acorns possibly depends extensively on the continuous erratic rainfall, i.e. continuous wet-dry cycle, but in-situ experiments are yet to be performed to test this hypothesis.

How ancient forest fragmentation and riparian connectivity generate high levels of genetic diversity in a microendemic Malagasy tree

Understanding landscape changes is central to predicting evolutionary trajectories and defining conservation practices. While human-driven deforestation is intense throughout Madagascar, exceptions in areas such as the Loky-Manambato region (north) raise questions regarding the causes and age of forest fragmentation. The Loky-Manambato region also harbours a rich and endemic flora, whose evolutionary origin remains poorly understood. We assessed the genetic diversity of an endangered microendemic Malagasy olive species (Noronhia spinifolia Hong-Wa) to better understand the vegetation dynamics in the Loky-Manambato region and its influence on past evolutionary processes. We characterized 72 individuals sampled across eight forests through nuclear and mitochondrial restriction-associated DNA sequencing data and chloroplast microsatellites. Combined population and landscape genetics analyses indicate that N. spinifolia diversity is largely explained by the current forest cover, highlighting a long-standing habitat mosaic in the region. This sustains a major and long-term role of riparian corridors in maintaining connectivity across these antique mosaic habitats, calling for the study of organismal interactions that promote gene flow.

Population genomics study for the conservation management of the endangered shrub Abeliophyllum distichum

Natural monuments are IUCN Category III protected areas that play an important role in biodiversity conservation as they provide species refuge and allow species migration. Despite their status, natural monuments are often confined to cultural and fragmented landscapes due to anthropogenic land-use demands. In this population genomic study, we surveyed 11 populations of the endemic shrub Abeliophyllum distichum Nakai (Oleaceae), including five natural monument habitats, covering its range-wide distribution in South Korea. Using 2,254 SNPs as markers, our results showed a mean expected heterozygosity (He) of 0.319, with populations in the central distribution showing significantly higher He than those at the periphery. There was no significant heterozygote deficiency and inbreeding among studied populations overall (FIS = −0.098), except for a single natural monument population (GS-NM147). Population structure and differentiation was moderate to high (FST = 0.196), while recent gene flow between populations appeared weak, which can be attributed to the fragmented distribution and the outcrossing mating system of the heterostylous plant. Based on these findings, we provide suggestions for the population conservation and management of this endangered species.

Genetic Diversity and Population Structure of Medemia argun (Mart.) Wurttenb. ex H.Wendl. Based on Genome-Wide Markers

Medemia argun is a wild, dioecious palm, adapted to the harsh arid environment of the Nubian Desert in Sudan and southern Egypt. There is a concern about its conservation status, since little is known about its distribution, abundance, and genetic variation. M. argun grows on the floodplains of seasonal rivers (wadis). The continuing loss of suitable habitats in the Nubian Desert is threatening the survival of this species. We analyzed the genetic diversity, population genetic structure, and occurrence of M. argun populations to foster the development of conservation strategies for M. argun. Genotyping-by-sequencing (GBS) analyses were performed using a whole-genome profiling service. We found an overall low genetic diversity and moderate genetic structuring based on 40 single-nucleotide polymorphisms (SNPs) and 9,866 SilicoDArT markers. The expected heterozygosity of the total population (HT) equaled 0.036 and 0.127, and genetic differentiation among populations/groups (FST) was 0.052 and 0.092, based on SNP and SilicoDArT markers, respectively. Bayesian clustering analyses defined five genetic clusters that did not display any ancestral gene flow among each other. Based on SilicoDArT markers, the results of the analysis of molecular variance (AMOVA) confirmed the previously observed genetic differentiation among generation groups (23%; p < 0.01). Pairwise FST values indicated a genetic gap between old and young individuals. The observed low genetic diversity and its loss among generation groups, even under the detected high gene flow, show genetically vulnerable M. argun populations in the Nubian Desert in Sudan. To enrich and maintain genetic variability in these populations, conservation plans are required, including collection of seed material from genetically diverse populations and development of ex situ gene banks.

Multi-Approach Analysis Reveals Local Adaptation in a Widespread Forest Tree of Reunion Island

Detecting processes of local adaptation in forest trees and identifying environmental selective drivers are of primary importance for forest management and conservation. Transplant experiments, functional genomics and population genomics are complementary tools to efficiently characterize heritable phenotypic traits and to decipher the genetic bases of adaptive traits. Using an integrative approach combining phenotypic assessment in common garden, transcriptomics and landscape genomics, we investigated leaf adaptive traits in Coffea mauritiana, a forest tree endemic to Reunion Island. Eight populations of C. mauritiana originating from sites with contrasted environmental conditions were sampled in common garden to assess several leaf morphological traits, to analyze the leaf transcriptome and leaf cuticular wax composition. The relative alkane content of cuticular waxes was significantly correlated with major climatic gradients, paving the way for further transcriptome-based analyses. The expression pattern of cuticle biosynthetic genes was consistent with a modulation of alkane accumulation across the population studied, supporting the hypothesis that the composition of cuticular wax is involved in the local adaptation of C. mauritiana. Association tests in landscape genomics performed using RNA-seq-derived single-nucleotide polymorphisms revealed that genes associated with cell wall remodeling also likely play an adaptive role. By combining these different approaches, this study efficiently identified local adaptation processes in a non-model species. Our results provide the first evidence for local adaptation in trees endemic to Reunion Island and highlight the importance of cuticle composition for the adaptation of trees to the high evaporative demand in warm climates.

Genetic diversity and structure of Saussurea polylepis (Asteraceae) on continental islands of Korea: Implications for conservation strategies and management

Saussurea polylepis Nakai is an herbaceous perennial endemic to Korea and is highly restricted to several continental islands in the southwestern part of the Korean Peninsula. Given its very narrow geographical distribution, it is more vulnerable to anthropogenic activities and global climate changes than more widely distributed species. Despite the need for comprehensive genetic information for conservation and management, no such population genetic studies of S. polylepis have been conducted. In this study, genetic diversity and population structure were evaluated for 97 individuals from 5 populations (Gwanmaedo, Gageodo, Hongdo, Heusando, and Uido) using 19 polymorphic microsatellites. The populations were separated by a distance of 20–90 km. We found moderate levels of genetic diversity in S. polylepis (Ho = 0.42, He = 0.43). This may be due to long lifespans, outcrossing, and gene flow, despite its narrow range. High levels of gene flow (Nm = 1.76, mean Fst = 0.09), especially from wind-dispersed seeds, would contribute to low levels of genetic differentiation among populations. However, the small population size and reduced number of individuals in the reproductive phase of S. polylepis can be a major threat leading to inbreeding depression and genetic diversity loss. Bayesian cluster analysis revealed three significant structures at K = 3, consistent with DAPC and UPGMA. It is thought that sea level rise after the last glacial maximum may have acted as a geographical barrier, limiting the gene flow that would lead to distinct population structures. We proposed the Heuksando population, which is the largest island inhabited by S. polylepis, as a source population because of its large population size and high genetic diversity. Four management units (Gwanmaedo, Gageodo, Hongdo-Heuksando, and Uido) were suggested for conservation considering population size, genetic diversity, population structure, unique alleles, and geographical location (e.g., proximity).

Drivers of population divergence and species differentiation in a recent group of indigenous orchids (Vanilla spp.) in Madagascar

With over 25,000 species, orchids are among families with remarkable high rate of diversification. Since Darwin's time, major advances attributed the exceptional diversity of orchids to plant-pollinator interactions. However, unraveling the processes and factors that determine the phenotypic and genotypic variation of natural orchid populations remains a challenge. Here, we assessed genetic population structure and floral differentiation in recently diverged leafless Vanilla species in a world biodiversity hotspot, Madagascar, using seven microsatellite loci and 26 morphometric variables. Additionally, analyses were performed to test for the occurrence of any patterns of isolation by distance, isolation by environment, and isolation by adaptation and to detect possible physical barriers that might have caused genetic discontinuities between populations. Positive inbreeding coefficients detected in 22 populations were probably due to the presence of null alleles, geitonogamy and/or some admixture (sympatric species). In contrast, the only high-altitude population showed an important rate of clonality leading to heterozygote excess. Genetic diversity was maximum in western populations, suggesting a postglacial colonization to the north and south. Clustering analyses identified seven genetic groups characterized by specific floral traits that matched five botanical descriptions in the literature. A contribution of montane refugia and river barriers on population differentiation was detected. We also detected combined effects of IBD/IBE and IBE/IBA on genetic differentiation and suggested this pattern is more likely determined by ecological isolation, although pollinator-mediated divergent selection could not be ruled out for some of the species. Overall, this study provides further insights on speciation in orchids, a group for which Madagascar shows one of the world's highest level of endemism and confirms the importance of the peculiar biogeography of the island in shaping species differentiation.

Terminalia bentzoë, a Mascarene Endemic Plant, Inhibits Human Hepatocellular Carcinoma Cells Growth In Vitro via G0/G1 Phase Cell Cycle Arrest

Tropical forests constitute a prolific sanctuary of unique floral diversity and potential medicinal sources, however, many of them remain unexplored. The scarcity of rigorous scientific data on the surviving Mascarene endemic taxa renders bioprospecting of this untapped resource of utmost importance. Thus, in view of valorizing the native resource, this study has as its objective to investigate the bioactivities of endemic leaf extracts. Herein, seven Mascarene endemic plants leaves were extracted and evaluated for their in vitro antioxidant properties and antiproliferative effects on a panel of cancer cell lines, using methyl thiazolyl diphenyl-tetrazolium bromide (MTT) and clonogenic cell survival assays. Flow cytometry and comet assay were used to investigate the cell cycle and DNA damaging effects, respectively. Bioassay guided-fractionation coupled with liquid chromatography mass spectrometry (MS), gas chromatography-MS, and nuclear magnetic resonance spectroscopic analysis were used to identify the bioactive compounds. Among the seven plants tested, Terminaliabentzoë was comparatively the most potent antioxidant extract, with significantly (p < 0.05) higher cytotoxic activities. T. bentzoë extract further selectively suppressed the growth of human hepatocellular carcinoma cells and significantly halted the cell cycle progression in the G0/G1 phase, decreased the cells' replicative potential and induced significant DNA damage. In total, 10 phenolic compounds, including punicalagin and ellagic acid, were identified and likely contributed to the extract's potent antioxidant and cytotoxic activities. These results established a promising basis for further in-depth investigations into the potential use of T. bentzoë as a supportive therapy in cancer management.

The East Asian Winter Monsoon Acts as a Major Selective Factor in the Intraspecific Differentiation of Drought-Tolerant Nitraria tangutorum in Northwest China

The influence of Quaternary climate fluctuation on the geographical structure and genetic diversity of species distributed in the regions of the Qinghai–Tibet Plateau (QTP) has been well established. However, the underlying role of the East Asian monsoon system (EAMS) in shaping the genetic structure of the population and the demography of plants located in the arid northwest of China has not been explored. In the present study, Nitraria tangutorum, a drought-tolerant desert shrub that is distributed in the EAMS zone and has substantial ecological and economic value, was profiled to better understand the influence of EAMS evolution on its biogeographical patterns and demographic history. Thus, the phylogeographical structure and historical dynamics of this plant species were elucidated using its five chloroplast DNA (cpDNA) fragments. Hierarchical structure analysis revealed three distinct, divergent lineages: West, East-A, and East-B. The molecular dating was carried out using a Bayesian approach to estimate the time of intraspecies divergence. Notably, the eastern region, which included East-A and East-B lineages, was revealed to be the original center of distribution and was characterized by a high level of genetic diversity, with the intraspecific divergence time dated to be around 2.53 million years ago (Ma). These findings, combined with the data obtained by ecological niche modeling analysis, indicated that the East lineages have undergone population expansion and differentiation, which were closely correlated with the development of the EAMS, especially the East Asian winter monsoon (EAWM). The West lineage appears to have originated from the migration of N. tangutorum across the Hexi corridor at around 1.85 Ma, and subsequent colonization of the western region. These results suggest that the EAWM accelerated the population expansion of N. tangutorum and subsequent intraspecific differentiation. These findings collectively provide new information on the impact of the evolution of the EAMS on intraspecific diversification and population demography of drought-tolerant plant species in northwest China.

Modelling and simulating Chikungunya spread with an unstructured triangular cellular automata

In this work we propose a mathematical model to simulate Chikungunya spread; the spread model is implemented in a C++ cellular automata code defined on unstructured triangular grids and space visualizations are performed with Python. In order to simulate the time space spread of the Chikungunya diseases we include assumptions such as: heterogeneous human and vector densities, population mobility, geographically localized points of infection using geographical information systems, changes in the probabilities of infection, extrinsic incubation and mosquito death rate due to environmental variables. Numerical experiments reproduce the qualitative behavior of diseases spread and provide an insight to develop strategies to prevent the diseases spread.

Genetic Diversity of Phyllanthus emblica From Two Different Climate Type Areas

Phyllanthus emblica L. is a well-known medicinal and edible plant species. Various medicinal compounds in the fruit make it an important medicinal and promising economic material. The plant is widely distributed in Southwestern and Southern China. However, due to massive deforestation and land reclamation as well as deterioration of its natural habitat in recent years, the wild resources of this species have been sharply reduced, and it is rare to see large-scale wild P. emblica forests so far. In order to effectively protect and rationally utilize this species, we investigated the genetic diversity, genetic structure, and population dynamics of 260 individuals from 10 populations of P. emblica sampled from the dry climate area in Yunnan and wet climate area in Guangxi using 20 polymorphic EST-SSR markers. We found high genetic diversity at the species level (He = 0.796) and within populations (He = 0.792), but low genetic differentiation among populations (F _ST = 0.084). In addition, most genetic variation existed within populations (92.44%) compared with variation among the populations (7.56%). Meanwhile, the NJ tree, STRUCTURE, and hierarchical analysis suggested that the sampled individuals were clustered into two distinct genetic groups. In contrast, the genetic diversity of the dry climate group (He = 0.786, Na = 11.790, I = 1.962) was higher than that of the wet climate group (He = 0.673, Na = 9.060, I = 1.555), which might be attributed to the combined effects of altitude, precipitation, and geographic distance. Interestingly, only altitude and precipitation had significant pure effects on the genetic diversity, and the former was slightly stronger. In addition, DIYABC analysis suggested the effective population size of P. emblica might have contracted in the beginning of the Last Glacial Maximum. These genetic features provided vital information for the conservation and sustainable development of genetic resources of P. emblica, and they also provided new insights and guidelines for ecological restoration and economic development in dry-hot valleys of Yunnan and karst areas in Guangxi.

Plant population dynamics on oceanic islands during the Late Quaternary climate changes: genetic evidence from a tree species (Coffea mauritiana) in Reunion Island

Past climatic fluctuations have played a major role in shaping the current plant biodiversity. Although harbouring an exceptional biota, oceanic islands have received little attention in studies on species demographic history and past vegetation patterns. We investigated the impact of past climatic changes on the effective population size of a tree (Coffea mauritiana) that is endemic to Reunion Island, located in the south-western Indian Ocean (SWIO). Demographic changes were inferred using summary statistics calculated from genomic data. Using ecological niche modelling and the current distribution of genetic diversity, the paleodistribution of the species was also assessed. A reduction in the effective population size of C. mauritiana during the last glaciation maximum was inferred. The distribution of the species was reduced on the western side of the island, due to low rainfall. It appeared that a major reduction in rainfall and a slight temperature decrease prevailed in the SWIO. Our findings indicated that analyses on the current patterns of intraspecific genetic variations can efficiently contribute to past climatic changes characterisation in remote islands. Identifying area with higher resilience in oceanic islands could provide guidance in forest management and conservation faced to the global climate change.