Exploring the adaptive mechanism of Camellia limonia in karst high calcium environment via an integrative analysis of metabolomics and metagenomics

Exploring the adaptive mechanism of Camellia limonia in karst high calcium environment via an integrative analysis of metabolomics and metagenomics | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Exploring the adaptive mechanism of Camellia limonia in karst high calcium environment via an integrative analysis of metabolomics and metagenomics Chunni Liu, Yang Huang, Feng Wu, Wenjing Liu, Yiqiu Ning, Zhenrong Huang, and 1 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-135639/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract Background: Yellow Camellia is a kind of rare plant with high economic and medicinal value. It is known as the "giant panda" of the botanical world. Camellia group of plants (Camellia Sect. Chrysantha Chang) is the only golden flower of Theaceae . Compared with Camellia nitidissima, Camellia limonia grows in karst areas, where the soil has the characteristic of high calcium content. However, there are few studies about the Camellia limonia in karst soil environment and the adaptation mechanism is no clear. Results: In this study, we found that under high calcium treatment, the chlorophyll content and leaf areas of Camellia limonia increased, while those of Camellia nitidissima decreased. The photosynthetic efficiency of Camellia limonia was more stable and higher than Camellia nitidissima . Compared with Camellia nitidissima , the conductance was larger and the degree of leaf shrinkage was smaller in Camellia limonia . The metabolomics analysis showed that the kaempferol-3-o-rutinoside, tyrosol, 6-o-methyldeacetylisoipecoside and (r)-mandelic acid are the main differently metabolic compounds . The results of karst high calcium soil metagenomics showed that microbacterium-testaceum, intrasporangium-calvum and rubrivivax-gelatinosus significantly changed. Through metabonomics and metagenomics integrative analysis, flavone and flavonol biosynthesis is suggested to be the main regulation pathway, which is regulated by apigenin, kaempferol, astragalin, isoquercitrin metabolites and TT7, UGT78D1 , UGT78D2 genes. This metabolic pathway involves the synthesis of flavonoids. Flavonoids have the functions in drought and salt resistance, which play an important role in the adaptation of Camellia limonia in karst high calcium environment. Conclusion : This omics study identified key regulation metabolites and genes for Camellia and provided important basis for the adaptive mechanism of plants to adapt to the high-calcium environment and the protection of Camellia species. Plant Physiology and Morphology Plant Molecular Biology and Genetics Camellia limonia karst high calcium stress metagenomics metabonomics Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Figure 8 Introduction Camellia is originated from China and have golden flowers, glossy leaves and beautiful appearance with transparent waxy surface [1]. It is a kind of rare plant that has not only high ornamental and medicinal value but also is very useful genetic resource. Camellia is mainly distributed in the southern China and the north of Vietnam [2]. Camellia species are heat-loving tropical plants, belonging to the tropical tree species, which have a high requirement for heat and mainly grow under shady and moist evergreen broad-leaf forests or along streamside in the mountains with the altitudes that is less than 700 meters [3]. Camellia limonia are an evergreen shrub with 1-2m in height. And their elliptic leathery leaves are 5-10cm in length and 2-3.5cm in width with in margin serrate. Their base is broadly wedge-shaped with dark brown glandular spots below. Flowers solitarily have 8 lemon yellow petals [1]. Camellia limonia mainly grow in karst rock mountain area, whereas Camellia nitidissima grow in normal soil mountain. Karst landscape is a typical ecosystem with very fragile ecological environment, mainly formed by carbonate dissolution. Karst regions cover about 20% of the dry ice-free land of the earth [4]. It is approximately 22 million square kilometers, accounting for about 12% of the total land area in the world[5]. China is one of the countries with the extensive karst landscapes and is famous for its diverse karst landforms[6]. Chinese karst landscape covers an area of 3.443 million square kilometers, accounting for one third of the country's land area. In southwest China, karst landforms cover an area of about 540,000 square kilometers, which is the largest and most representative karst fragile environment[6]. Located in the humid subtropical monsoon climate zone, the rainfall is abundant, but the precipitation has seasonal characteristics [7]. However, due to the obvious seasonal rainfall and dual hydrogeological structure, water in the karst area is formed only during the rainy season. As a result, karst soil water storage capacity is very limited, and has significant seasonal variation, which result in the shortage of water resources [4]. Because of the development of rock crevasses and falling-water caves, the situation of soil erosion becomes more serious. In these regions, the soil is shallow and discontinuous. Vegetation grows slowly with low biomass [5]. Karst landforms are mainly distributed in monsoon and valley rain forests, which is beneficial for the growth of Camellia . The soil types, ashes and elemental components (especially calcium) content of the Camellia in the two landform types are significantly different. Studies have shown that the leaf gray matter and calcium contents of Camellia in karst landforms are higher than the Camellia that grows in normal erosion landforms [8]. Karst is a very special environment whose shape is mainly the typical topography of the surface and underground part is formed by dissolving water on soluble rocks, which is the most typical topography in karst landforms [9]. Moreover, the soil in the karst areas is chiefly formed by weathering of bedrock carbonate rock, which causes the soil in the karst area having higher calcium content. Plants in karst areas have been adapted to grow in a drought and high calcium environment for a long time through their own physiological processes and structural adjustments in the evolution. Soil moisture and calcium carbonate content have the greatest impact on vegetation spectral characteristics in karst areas [10]. In the present research, we focused on the adaptation mechanism and physiological changes of Camellia limonia in the karst high calcium environment. In order to systematically obtain the regulation mechanism of Camellia limonia in the karst high calcium environment and find key pathway and important metabolites, integrative analysis of metabonomics and metagenomics was conducted and Camellia nitidissima was used as comparison. It is shown that flavone and flavonol biosynthesis is the main regulation pathway, which involves the synthesis of flavonoids. Flavonoids further play roles in functions in drought and salt stress, which might contribute to the karst adaptation of Camellia limonia . This study can improve the research about adaptive mechanism of plants in karst environments. Also, a multi omics study of Camellia can provide a practical foundation and valuable data resources for endangered plant conservation research. Methods Plant material and high calcium stress treatment Camellia limonia and Camellia nitidissima used in this study were collected from Guangxi Institute of Botany. All the plants grew in outdoor greenhouse (located in Guangxi Normal University, Guilin, China) for shading cultivation. Ten plants of each variety with the same phenotype were selected as an experimental group. Each treatment plant was irrigated with 300mM calcium chloride solution every two days, while the control group was irrigated with pure water. Leaves were selected for metabolomics. Soil ponds with roots were selected as metagenome. The samples were collected on the 15th day after treatment when the phenotype of plants showed significant differences between Camellia limonia and Camellia nitidissima . Morphological measurement The morphological data of plants in high calcium treatment were collected every 5 hours from 7:00 to 18:00 every day. For each plant, at least 3 healthy leaves with the same leaf size and age were randomly selected to measure the photosynthetic parameters, leaf area, temperature, humidity, nitrogen content and chlorophyll content of leaves, respectively. Photosynthetic parameters were measured by Li-6400 photosynthetic apparatus [11]. Plant morphological meter TYS-4N was used to measure nitrogen content, temperature and humidity of leaves. The parameters of leaf length, width and leaf area of plant leaves were measured by a leaf area meter YMJ-B. By repeated freezing and thawing extraction to extract the chlorophyll of leaves were extracted and measured by the methods mentioned in previous relative references [12]. The Camellia nitidissima and Camellia limonia leaves after were selected for scanning electron microscope (SEM) and paraffin sections conservation based on the previous study [13]. Soil moisture, soil temperature and chemical element content were measured daily, including nitrogen (N), phosphorus (P), and potassium (K). The measurements were taken every five hours from 7:00 to 18:00 every day after high calcium treatment. The contents of N, P and K were measured by extracting soil of plants roots [14]. The PH, temperature and humidity of soil were measured by SOIL SENSOR. Metabolomic analysis The leaves under calcium treatment for 15 days were used to conduct metabolomic analysis. All the samples were stored at -80℃ refrigerator and six biological repeat for Camellia limonia and Camellia nitidissima , respectively. Multivariate analysis technologies including Principal Component Analysis (PCA) and Partial Least Squares (PLS) regressions were used at present study to identify several uncorrelated variables (principal components, PCS) from a large set of data [15]. Additionally, methodological advances in untargeted metabolomics and data analysis of tandem mass-spectrometric (MS/MS) data have enabled comprehensive analysis of molecular components in metabolic extracts [16]. These MS/MS spectra can be viewed as barcodes or fingerprints of metabolites and some software tools have been developed to take advantage of this structural information [17]. Metagenomics analysis In order to study the microbial community composition and different express genes in the soil microorganism, we selected plants with obvious different phenotype traits after high calcium treatment to take the soil samples. The soil with the selected plant was taken with a small amount of roots from each plant. Then the samples were stored at -80℃ refrigerator until sequencing and two biological repeat for Camellia limonia and Camellia nitidissima , respectively [18]. Illumina PE sequencing strategy was used to extract total DNA from the soil and then agarose gel electrophoresis was used to detect the quality of the samples and finally the library was obtained. In addition, the Velvet software was used for sequence assembly and finally functional annotation and comparison analysis were performed according to the literature [19, 20]. Data processing and statistics analysis Data analysis was conducted by using the SPSS 20.0 program [4]. Independent-sample was used to detect the indexes of morphology and soil. Leaf morphology indicators include: net photosynthetic rate, stomatal conductance, transpirational rate, leaf area, chlorophyll, nitrogen content, leaf surface moisture and leaf surface temperature; Soil morphology indicators include: pH, temperature and humidity. The methods of mathematical statistics are used to calculate the error value and the significant difference. The results were taken into account as not significant at P > 0.05, significant at 0.01<P<0.05, and extremely significant at P<0.01. PCA and PLS are used to discriminate and analyze the error value and significant difference. A cluster analysis of the identified metabolites was conducted after auto-scaling (mean-centered and divided by the standard deviation) by using MultiExperiment Viewer [21]. Results E ff ects of calcium on the karst adaptive characteristics of Camellia Abiotic stress signals can cause fluctuations of cytoplasmic calcium ions and changes in phosphorylation forms of some proteins. Calcium ions mediated by calcium channels, change the concentration of calcium ions in cytoplasm and amplify the stress signals with intracellular calcium signal receptors [22]. However, high calcium stress can destroy the cytoskeleton and biofilm structure, resulting in intracellular metabolic imbalance. In order to study the adaptability of the karst environment with high calcium content of Camellia limonia , we first processed Camellia limonia with high calcium content and take Camellia nitidissima as the control, the results showed that the adaptability of Camellia limonia to the environment with high calcium content was stronger. Under the calcium stress, the phenotypes of Camellia nitidissima and Camellia limonia became different. The leaves of Camellia nitidissima turned yellow and wilting, while the leaves of Camellia limonia kept relatively healthy (Figure 1). The SEM and microscopic observation also showed that the leaves and roots of Camellia limonia are less wrinkled and dehydrated than Camellia nitidissima (Figure 2). The stomatal conductance and transpiration rate of the two kinds of Camellia showed a downward trend. The stomatal conductance of Camellia nitidissima decreased more than the one of Camellia limonia . While the reduce degree of transpiration rate Camellia limonia is less than in Camellia nitidissima . The chlorophyll content of Camellia limonia increased by 53.12% and Camellia nitidissima decreased by 56.91%. The net photosynthetic rate of Camellia limonia was more stable than Camellia nitidissima . The average net photosynthetic rate of Camellia limonia and Camellia nitidissima was 3.0853 μmol CO 2 m -2 s -1 and 2.4456 μmol CO 2 m -2 s -1 , respectively. Stomatal conductance affects the oxygen and carbon dioxide, indicating that the water might be the more important factor influencing the net photosynthetic rate when gas were equivalently affected [23]. According to the data, the average leaf humidity of Camellia limonia is 21.29% RH, Camellia nitidissima leaf humidity averages of 17.88% RH. The average transpiration rate and leaf temperature of Camellia nitidissima were 0.64 mmolH 2 Om -2 s -1 and 31.94 ℃， Camellia nitidissima were 0.92 mmolH 2 Om -2 s -1 and 30.95 ℃, suggesting the water content of Camellia limonia is more than Camellia nitidissima (Figure 3). According to the results above, Camellia nitidissima has relatively stronger adaptability to the environment with high calcium. The soil moisture and temperature of the two kinds of Camellia did not obviously change. The pH value of the soil of Camellia limonia decreased slightly. In addition, the nitrogen and phosphorus content in the soil of Camellia limonia increased by 58.60% and 53.99%, respectively, while the nitrogen and phosphorus content in the soil of Camellia nitidissima decreased by 13.01% and increased by 49.63%, respectively (Figure S1). Available nitrogen and phosphorus are the main nutrient elements that restrict plant growth, while plant growth is mainly restricted by nitrogen elements under low ratio of nitrogen and phosphorus [24]. Under high calcium stress, nitrogen content of Camellia limonia leaf increased by 14.44% and the available nitrogen in soil of Camellia limonia is increased by 58.60%, while nitrogen content of Camellia nitidissima leaf decreased by 28.57% and the available nitrogen in soil of Camellia nitidissima decreased by 13.01%, suggesting Camellia limonia has better nitrogen use efficiency (Figure 3), which may also help Camellia limonia survived in high calcium environment. Metabonomics of Camellia leaf in karst high calcium environment In order to analysis the different metabolites of Camellia leaf in karst high calcium environment, the metabonomics were conducted. The PCA score chart showed the degree of aggregation and dispersion of Camellia limonia and Camellia nitidissima (Figure S2). The distribution points of Camellia limonia are relatively close, indicating that the composition and concentration of the variables/molecules contained in Camellia limonia are closer. A total of 9133 metabolites were identified through mathematical statistical analysis and screening, of which 6422 were in the positive ion mode and 2711 were in the negative ion mode (Figure S3). 8966 differential metabolites discemed include metabolites of Camellia limonia and Camellia nitidissima . The accurate molecular weight of the metabolites (molecular weight error <15ppm) were confirmed and the fragment information were divided according to the MS/MS mode. We identified 220 different metabolites in HMDB (Human Metabolome Database, HMDB), metlin, massbank, lipiduaps, mzdound databases and obtained the statistical analysis results (Table S1). Under the high calcium stress, Camellia limonia has 53 metabolites up-regulated and 167 metabolites down-regulated, of which the main up-regulated metabolite is kaempferol-3-o-rutinoside and tyrosol. Camellia nitidissima has 167 up-regulated metabolites and 53 down-regulated metabolites, of which the main metabolites up-regulated are n-acetylhistidine, myristic acid, and pelargonin (Figure 4). The five different metabolites mentioned above have the largest significant difference. The content of these metabolites were changed to affect the metabolic process of Camellia , indicating that the significantly different metabolites mainly affect the growth of Camellia limonia and Camellia nitidissima after the treatment of high calcium stress. Among them, kaempferol-3-o-rutinoside, tyrosol, 6-o-methyldeacetylisoipecoside,n-acetylhistidine, myristic acid and pelargonin are positively correlated. But the three metabolites of Camellia limonia were negatively correlated with the three metabolites of Camellia nitidissima (Figure 5). It is inferred that the up-regulated metabolites of Camellia limonia had the opposite trend to the up-regulated metabolites of Camellia nitidissima , indicating that Camellia limonia have different metabolite trends in high calcium stress. MetPA is a part of metaboanalyst (www.metaboanalyst.ca) and is mainly based on the KEGG metabolic pathway. The MetPA related metabolic and KEGG pathway analysis show 42 main metabolic pathways. The most important is biotin metabolism, tyrosine metabolism, flavone and flavonol biosynthesis, inositol phosphate metabolism, alpha-linolenic acid metabolism, ubiquinone and other terpenoid-quinone biosynthesis, phenylpropanoid biosynthesis, lysine biosynthesis, phenylalanine metabolism, glyoxylate and dicarboxylate metabolism, pantothenate and CoA biosynthesis, arginine and proline metabolism and other metabolic (Table S2). The pathways with most significant differences are tyrosine metabolism, flavone and flavonol biosynthesis, phenylpropanoid biosynthesis and biotin metabolism. And more differential metabolites are enriched in these four pathways, indicating that these pathways are more sensitive to high calcium stress for Camellia limonia . Metagenomics of soil in karst high calcium environment Under biotic and abiotic stresses, plants can actively change their own physiological conditions and exudates of their roots. That also has effects on the related microbial communities. Consequently, this directly or indirectly affects the health of plants. The initial central microorganism and exudates of plant root together form a microbial community, including plant central flora, core flora and other microbial communities. Compared with the plants under high calcium treatment, the root microbial populations of the plants in normal condition are different. Therefore, these different populations could affect the growth of plants. Metagenomics analysis of the soil of Camellia limonia and Camellia nitidissima root were conducted to study the effect of microbial communities under high calcium stress. The result showed that tyrosine metabolism, galactose metabolism, inositol phosphate metabolism, biotin metabolism, phenylpropanoid biosynthesis, arginine and proline metabolism were significantly affected. The microorganism regulated these metabolic pathways. At the same time, they regulate the expression of responsive proteins through responsive genes. Finally, these microorganisms influnce the metabolic pathways of Camellia . The software Velvet was used to evaluate the assembly results, including the statistical table of assembly results, the distribution chart of the cumulative length of contigs and the distribution map of GC content (Table 1). We examined nearly 20,000 genes from each sample of rhizosphere microbes. And we found that reduced pathogenicity and loss of toxicity were the most important factors. In order to get a more detailed classification result of species, we match the sequence with the bacterial reference sequence or the virus reference sequence. Moreover, the MetaPhlAn was used to comment the species classification (Table 2). In all samples, rubrivivax is the main member. In root soil of Camellia limonia , rhodanobacter, intrasporangium, microbacterium, ilumatobacter are taking disadvantages. Different from Camellia limonia , the Candidatus_Solibacter was the dominant species of Camellia nitidissima (Figure 6). Further, in order to reveal the distribution and evolution relationship of species, phylogenetic analysis was conducted (Figure S4). The three colors represent the most important four species ( dehalogenimonas lykanthroporepellens , micavibrio aeruginosavorus , polymorphum gilvum and candidatus accumulibacter phosphatis ) and their evolution relationship. The above results about the metagenomic analysis indicated that, compared with Camellia nitidissima , the types of soil microorganisms of Camellia limonia under high calcium stress were different. And these microorganisms are related to the salt tolerance and drought resistance of plants. Camellia limonia adaptive regulation to karst In order to study the karst adaptation of Camellia limonia , we conducted the integrative analysis of metabolomics and metagenomics. The result showed that flavone and flavonol biosynthesis pathway which is one of the most affected pathways can be regulated by apigenin, kaempferol, astragalin, isoquercitrin metabolites and TT7, UGT78D1, UGT78D2 genes. This metabolic pathway involves the synthesis of flavonoids. Flavonoids have functions in drought and salt resistance [25], so flavonoids may also contribute to the karst environment adaptation for Camellia limonia . Flavonoids include kaempferol, rutin, hesperidin, anthocyanin, etc. Among them, kaempferol has the largest difference in content between Camellia limonia and Camellia nitidissima , and the content of Camellia limonia is greater than Camellia nitidissima . Therefore, Camellia limonia grow better than Camellia nitidissima under the same calcium stress, which indicated that kaempferol may play a critical role in the high calcium tolerance. The citrate cycle (TCA cycle, Krebs cycle) is an important aerobic pathway for the final steps of the oxidation of carbohydrates and fatty acids. The citrate cycle is also an important pathway for differential metabolites. In the relative pathway of TCA cycle, the sdhA, frdA genes regulate the synthesis of succinate dehydrogenase/fumarate reductase-flavoprotein subunit. Compared with Camellia nitidissima , the sdhA, frdA were differently expression genes in Camellia limonia , which suggested that this two genes might simultaneously regulate the synthesis of fumarate and succinate and further affect the high calcium tolerance of Camellia limonia . Discussion Changes in gene regulation and metabolic material network lead to Camellia karst adaptation Calcium ion is a necessary mineral element for plant growth and development [26]. It is also an important component of plant cell wall and has the function of maintaining intracellular ion balance and cell membrane stability [27]. Calcium ions in the cytoplasm are generally kept at a low concentration, and high concentration of calcium ions will disturb the physiological functions of plant cells [28]. The absorption of calcium by plants is positively correlated with the concentration of calcium ions in soil. If the concentration of calcium ions in soil is too high, the plants will suffer the high calcium stress. High concentration of calcium ions can seriously affect plant photosynthesis and inhibit plant development [29]. From the microscopic point of view, calcium ions function as a signal molecule of plant cells to environmental changes, but if the concentration of calcium ions in cells is too high, it will lead to cell dysfunction [30]. The phenotypic characteristics of Camellia nitidissima were significantly different from Camellia nitidissima after high calcium treatment. The leaves of Camellia nitidissima are easier to fall off, which is probably caused by increased abscisic acid in plants after high calcium treatment [31]. Moreover, the color of the Camellia leaves has changed obviously. Compared with Camellia nitidissima , leaf temperature, soil moisture and potassium contents of Camellia limonia are basically maintained at the same level after high calcium treatment. Chlorophyll content, phosphorus content, nitrogen content and leaf moisture of Camellia limonia are generally higher than Camellia nitidissima . While the pH of soil, transpiration rate of Camellia limonia are relatively lower than Camellia nitidissima . The Camellia plants lost much water which was mainly caused by higher transpiration under the calcium treatment. For Camellia limonia , its transpiration effect was relatively weak, indicating that this may be one of the results of adapting to the karst high calcium environment. Camellia limonia can properly undermine transpiration, which is beneficial to live in a karst high calcium environment. Transport of Ca 2+ to leaves depends on several factors, such as xylem sap Ca 2+ concentration, water absorption, transpiration and so on. Calcium accumulation may occur in organs with high transpiration. Low transpiration rates can ensure plants well hydrated, which is beneficial to absorb the extra calcium [32]. Therefore, Camellia limonia can adapt better to the high calcium karst environment under a certain transpiration intensity. While Camellia nitidissima had excessive water loss in the plant body due to relatively strong transpiration, which could affect various life activities and could not survive in karst high calcium environment. The metabolites of Camellia have undergone tremendous changes under calcium stress. Some key genes regulate the synthesis and changes of metabolites in the high calcium environment (Figure 7; Figure 8). Among the changed metabolism pathways, biotin metabolism has the highest impact factor. At the same time, flavone and flavonol biosynthesis involve many metabolic pathways and their impact factors are also high. Flavone and flavonol biosynthesis, ubiquinone and other terpenoid-quinone biosynthesis, phenylpropanoid biosynthesis and flavonoid biosynthesis play a regulatory role in metabolite changes and regulation and some of these pathways are interconnected (Figure 7). Special metabolites take part in the defense of plants with both biotic and abiotic environmental factors [33]. Glycosylation is a common modification of the structure of plant specific metabolites. It broadens the structural diversity and improves the substrate molecule solubility, stability and transportability [34]. A variety of glycosylated phenolic metabolites have been accumulated in plant tissues, such as anthocyanins and flavonoidsm [35]. Our results show that flavonoid biosynthesis and anthocyanin biosynthesis are the metabolic pathways of anthocyanins and flavonoids. This suggests that anthocyanins and flavonoids produced by these two pathways could contribute to resisting adverse environment, such as over-acidity, over-alkalinity, high calcium and drought. It is speculated that the up-regulation of these two metabolites may play the key roles for Camellia limonia adaptation to the karst high calcium environment. Moreover, the changes of different metabolites are also related to the overall metabolomic changes. Under high calcium stress, the chlorophyll content of the leaves of Camellia nitidissima decreased significantly and the leaves grew wrinkled (Figure 1; Figure 2). The leaves of Camellia limonia lack green luster on the surface, but their chlorophyll content remained relatively high. Plants with strong drought resistance have higher chlorophyll content and less biomass loss [36]. High calcium in the soil could cause the plants run out of water. Our result shows that the chlorophyll content of Camellia limonia is higher than Camellia nitidissima, which further suggested that enough chlorophyll content could ensure plant stability under high calcium condition. Under the high calcium treatment, some important pathways are very abundant, including energy metabolism associated with TCA Cycle, pathways associated with secondary metabolites, such as flavonol biosynthesis and flavonoids and other important pathways (such as biotin metabolism). All these pathways could also function in the adaptation of Camellia limonia . Response of Camellia to karst high calcium environment and its regulation model In order to further explore the regulation mechanism, we mapped the integrative analysis pathways, including differential metabolites, genes and enzymes (Figure 8). The analysis focused on the relationship of the pathways with higher influencing factors, such as biotin, flavone and flavonol biosynthesis, flavonoid biosynthesis, tyrosine metabolism, inositol phosphate metabolism. The results showed that tyrosine, p-coumaroy-CoA and fumarate are important metabolites. They are the metabolites connecting key metabolic pathways. L-tyrosine and fumarate connected tyrosine metabolism, phenylalanine tyrosine and tryptophan, TCA cycle, phenylpropanoid biosynthesis and other metabolic pathways. P-coumaroy-CoA played a very important role in phenylpropanoid biosynthesis. P-coumaroy-CoA is also an important hub in the regulation network. It links flavone and flavonol biosynthesis, flavonoid biosynthesis, phenylpropanoid biosynthesis. P-coumaroy-CoA was involved in the phenylpropanoid pathway. P-coumaroy-CoA mainly catalyzed the synthesis of 4-coumaroyl-CoA with the substrates of 4-coumaric acid and coenzyme A. The reaction is one of the rate-limiting steps which greatly affect the yields of lignins and flavonoids in phenylpropanoid pathway [37]. In the metabolic regulation processes, enzymes act as biocatalysts in many important bioreactions. In the phenylpropanoid biosynthesis, L-tyrosine is continuously catalyzed by 4-coumarate-CoA ligase, which is regulated by gene ACOS5 . Shikimate O-hydroxycinnamoyl transferase is regulated by gene SHT (AT2G19070). 5-O-(4-coumarayl)-D-quinate 3’-monooxygenase is regulated by gene CYP98A3 . Caffeate O-methyltransferase is regulated by gene OMT1 (AT5G54160). All these genes mentioned above were changed in the present study. These enzymes finally lead to the synthesis of differential metabolite ferulic acid. Phenylpropanoid biosynthesis have an overlapping reaction chain. L-tyrosine can be turned to P-coumaric acid. P-coumaric acid can be transformed to P-coumaroyl-CoA catalyzed by 4-coumarate-CoA ligase (4CL). ACOS5 and CYP98A3 regulated the synthesis of coumarin. Coumarin is involved in the flavonoid synthesis pathway [38]. Studies about Arabidopsis thaliana have shown that flavonoids and flavonoids affect the drought-resistance of plants [39]. OMT (AT5G54160) and SHT (AT2G19070) regulate the synthesis of lignin [40]. Studies about Arabidopsis thaliana have shown that lignin can improve plant tolerance in salt stress [41]. These pathways connect some key differently metabolites in the metabonomics of Camellia . Genes regulate the metabolic processes of Camellia which were also found similar function in other plants [41], suggesting these gene may play an important role in the karst adaption of Camellia . Plants in karst areas have been in the high calcium environment for a long time. And they have also responded to various environmental stresses through their own physiological processes and structural adjustments in the long-term adaptive evolution. But the regulation mechanism of Camellia in karst high calcium environment still remain unclear. The present research explored the regulation mechanism of Camellia in high calcium environment by integrative analysis of metabolomics and metagenomics. These omics results and analysis shown that some key metabolites and genes may affect the Camellia limonia adaption of karst high calcium environment and lead to the difference of phenotype of Camellia limonia and Camellia nitidissima. The comprehensive omics level study on Camellia can provide practical research foundation and valuable data resources for future plant adaptability research in karst environments and conservation research of rare and endangered plants. Declarations Ethics approval and consent to participate Not applicable. Consent for publication Not applicable. Availability of data and materials The datasets generated and analyzed during the present study are available from the corresponding author on reasonable request. Competing interests The authors declare that they have no competing interests. Authors' contributions Y.L., Y.H. and C.L. wrote the main manuscript text. F.W.prepared figures 1-3. W.L. prepared figures 4-5. Y.N. prepared figures 6 and all the tables. Z.H. prepared figure 8 and supplementary materials. All authors reviewed the manuscript. Acknowledgements We acknowledge the assistance from Xiao Wei and Jianmin Tang in Guangxi Institute of Botany, Chinese Academy of Sciences. Funding This work was supported by the Guangxi science and technology base and talent project(Guike AD19245098); Guangxi Natural Science Foundation(2020JJB130026); Research Project of Key Laboratory of Ecology of Rare and Endangered Species and Environmental Protection(ERESEP2020Z09); Guangxi graduate education innovation program(YCSW2020104); Innovation training program for university students(202010602251). References Wei X, Jiang YS, Jiang SY, Qi XX, Xiong ZC, Ye WH, et al. Photosynthetic characteristics of an endangered species Camellia nitidissima and its widespread congener Camellia sinensis . Photosynt. 2008;46:312-314. Tang S, Bin X, Wang L, Zhong Y. Genetic diversity and population structure of yellow Camellia ( Camellia nitidissima) in China as revealed by RAPD and AFLP markers. Biochem Genet. 2006;44:444–456. Song L, Wang X, Zheng X, Huang D. Polyphenolic antioxidant profiles of yellow camellia. Food Chemistry. 2011;129:351–357. Cao M, Wu C, Liu J, Jiang Y. Increasing leaf δ13C values of woody plants in response to water stress induced by tunnel excavation in a karst trough valley: Implication for improving water-use efficiency. Journal of Hydrology. 2020;586:124895. Zhang Z, Huang X, Zhang J. Soil thickness and affecting factors in forestland in a karst basin in Southwest China. Trop Ecol. 2020;61:267–277. Guo F, Jiang G, Yuan D, Polk JS. Evolution of major environmental geological problems in karst areas of Southwestern China. Environ Earth Sci. 2013;69:2427–2435. Zhang Z, Huang X, Zhang J. Soil thickness and affecting factors in forestland in a karst basin in Southwest China. Trop Ecol. 2020;61:267–277. Su Z. A preliminary study on the population ecology of Camellia sect. nitidissima. Guangxi Science. 1994;01:1. Pisano L, Zumpano V, Liso IS, Parise M. Geomorphological and structural characterization of the ‘Canale di Pirro’ polje, Apulia (Southern Italy). Journal of Maps. 2020;16:479–87. Xu M, Li A, Teng Y, Sun Z, Xu M. Exploring the adaptive mechanism of Passiflora edulis in karst areas via an integrative analysis of nutrient elements and transcriptional profiles. BMC Plant Biol. 2019;19:185. Zhang M, Chen Y, Du L, Wu Y, Liu Z, Han L. The potential of P aulownia fortunei seedlings for the phytoremediation of manganese slag amended with spent mushroom compost. Ecotoxicology and Environmental Safety. 2020;196:110538. Wun C, Rho J, Walker RobertW, Litsky W. A simplified method for the simultaneous extraction of phytoplanktonic chlorophyll and fecal sterol from water. Water Air Soil Pollut. 1979;11:1. Buszewska-Forajta M, Patejko M, Macioszek S, Sigorski D, Iżycka-Świeszewska E, Markuszewski MJ. Paraffin-embedded tissue as a novel matrix in metabolomics study: optimization of metabolite extraction method. Chromatographia. 2019;82:1501–1513. Simeonov V, Asenov I, Diadov V. Accelerated extraction of elements from soils by ultrasonic pretreatment. Z Anal Chem. 1977;285:252. Ongo E, Falasconi M, Sberveglieri G, Antonelli A, Montevecchi G, Sberveglieri V, et al. Chemometric discrimination of philippine civet coffee using electronic nose and gas chromatography mass spectrometry. Procedia Engineering. 2012;47:977–980. Shao SY, Ting Y, Wang J, Sun J, Guo XF. Characterization and identification of the major flavonoids in Phyllostachys edulis leaf extract by UPLC–QTOF–MS/MS. Acta Chromatographica. 2020;32:228–237. van der Hooft JJJ, Mohimani H, Bauermeister A, Dorrestein PC, Duncan KR, Medema MH. Linking genomics and metabolomics to chart specialized metabolic diversity. Chem Soc Rev. 2020;49:3297–3314. López-Mondéjar R, Tláskal V, Větrovský T, Štursová M, Toscan R, Nunes da Rocha U, et al. Metagenomics and stable isotope probing reveal the complementary contribution of fungal and bacterial communities in the recycling of dead biomass in forest soil. Soil Biology and Biochemistry. 2020;148:107875. Huang Y. The research of microbial genomics based on high-throughput sequencing. 2013; 12(5): 440-453. Huang X, Liu L, Wen T, Zhu R, Zhang J, Cai Z. Illumina MiSeq investigations on the changes of microbial community in the Fusarium oxysporum f.sp. cubense infected soil during and after reductive soil disinfestation. Microbiological Research. 2015; 181: 33-42. Yu X, Li Y, He C, Zhou J, Chen Y, Yu Z, et al. Nonvolatile metabolism in postharvest tea ( Camellia sinensis L.) leaves: Effects of different withering treatments on nonvolatile metabolites, gene expression levels, and enzyme activity. Food Chemistry. 2020;327:126992. Dietrich P. Physiology and biophysics of plant ligandgated ion channels. Plant Biology. 2010;12:80-93. Hazrati S, Tahmasebi-Sarvestani Z, Modarres-Sanavy S, Mokhtassi-Bidgoli A, Nicola S. Effects of water stress and light intensity on chlorophyll fluorescence parameters and pigments of Aloe vera L. Plant Physiology and Biochemistry. 2016;106:141–148. Hartman WH, Richardson CJ. Differential Nutrient Limitation of Soil Microbial Biomass and Metabolic Quotients (qCO 2 ): Is There a Biological Stoichiometry of Soil Microbes? PLOS ONE. 2013;8:14. Deng J, Qin W, Yang C, Iqbal N, Takpah D, Zhang J, et al. Seed quality deterioration dynamics for isoflavones biosynthesis in soybean ( Glycine max L. Merr.) seeds against field mildew stress. Acta Physiol Plant. 2019;41:57. Yang Q, Wei X, Zeng X, Ye W, Yin X, Zhang W, et al. Seed biology and germination ecophysiology of Camellia nitidissima . Forest Ecology and Management. 2008;255(1):113-118. Al-Whaibi MH, Siddiqui MH, Al-Amri A, Basalah MO. Performance of faba bean under calcium and gibberellic acid application. International Journal of Plant Developmental Biology. 2010; 4(1): 60-63. Silva LRD, Mansfield A. Where does all the calcium go? Evidence of an important regulatory role for trichomes in two calcicoles. Plant, Cell & Environment. 1996; 19(7): 880-886. Broadley MR. Variation in the shoot calcium content of angiosperms. Journal of Experimental Botany. 2003;54:1431–1446. Borer CH, Hamby MN, Hutchinson LH. Plant tolerance of a high calcium environment via foliar partitioning and sequestration. Journal of Arid Environments. 2012;85:128–31. Yoshida T, Fernie AR, Shinozaki K, Takahashi F. Long‐distance stress and developmental signals associated with abscisic acid signaling in environmental responses. The Plant Journal. 2020;15101. Bonomelli C, Alcalde C, Aguilera C, Videla X, Rojas-Silva X, Nario A, et al. Absorption and mobility of radio-labelled calcium in chili pepper plants and sweet cherry trees. Sci agric (Piracicaba, Braz). 2021;78:e20200092. Wu S, Tian L. A new flavone glucoside together with known ellagitannins and flavones with anti-diabetic and anti-obesity activities from the flowers of pomegranate ( Punica granatum ). Natural Product Research. 2019;33:252–7. Wilson AE, Wu S, Tian L. PgUGT95B2 preferentially metabolizes flavones/flavonols and has evolved independently from flavone/flavonol UGTs identified in Arabidopsis thaliana . Phytochemistry. 2019;157:184–93. Wu S, Tian L. Diverse phytochemicals and bioactivities in the ancient fruit and modern functional food pomegranate ( Punica granatum ). Molecules. 2017;22:1606. Gouveia CSS, Ganança JFT, Slaski JJ, Lebot V, Pinheiro de Carvalho MÂA. Involvement of abscisic acid and other stress indicators in taro (Colocasia esculenta (L.) Schott) response to drought conditions. Acta Physiol Plant. 2020;42:173. Liang L, Xiao W, Lei Y, Ling K. Quantitative evaluation of 4-coumarate: coa ligase (4cl) activity and correlated chemical constituents in four plant materials by chromatographic analysis. Chinese Journal of Natural Medicines. 2010; 8(4): 274-279 Siwinska J, Kadzinski L, Banasiuk R, Gwizdek-Wisniewska A, Olry A, Banecki B, et al. Identification of QTLs affecting scopolin and scopoletin biosynthesis in Arabidopsis thaliana. BMC Plant Biol. 2014;14:280. Karamat F, Olry A, Doerper S, Vialart G, Ullmann P, Werck-Reichhart D, et al. CYP98A22, a phenolic ester 3’-hydroxylase specialized in the synthesis of chlorogenic acid, as a new tool for enhancing the furanocoumarin concentration in Ruta graveolens. BMC Plant Biol. 2012;12(1):1-13. Nakatsubo T, Kitamura Y, Sakakibara N, Mizutani M, Hattori T, Sakurai N, et al. At5g54160 gene encodes Arabidopsis thaliana 5-hydroxyconiferaldehyde O-methyltransferase. J Wood Sci. 2008;54:312–317. Chun H J, Baek D, Cho H M, et al. Lignin biosynthesis genes play critical roles in the adaptation of Arabidopsis plants to high-salt stress. Plant signaling & behavior. 2019; 14(8): 1625697. Supplementary Files FigureS1.jpg FigureS2.jpg FigureS3.jpg FigureS4.jpg TableS1.xlsx TableS2.xlsx Supplementaryinformationtitlepagewithlegends.docx Cite Share Download PDF Status: Posted Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-135639","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":7280535,"identity":"614d92ad-72b5-4035-954a-191dbe5bad68","order_by":0,"name":"Chunni Liu","email":"","orcid":"","institution":"Guangxi Normal University","correspondingAuthor":false,"prefix":"","firstName":"Chunni","middleName":"","lastName":"Liu","suffix":""},{"id":7280536,"identity":"fdceae32-1127-4d67-9071-75fa87feedbb","order_by":1,"name":"Yang Huang","email":"","orcid":"","institution":"Guilin University of Electronic Technology","correspondingAuthor":false,"prefix":"","firstName":"Yang","middleName":"","lastName":"Huang","suffix":""},{"id":7280537,"identity":"b07419fa-a151-49ac-8fb5-2ebbec85b44e","order_by":2,"name":"Feng Wu","email":"","orcid":"","institution":"Guangxi Normal University","correspondingAuthor":false,"prefix":"","firstName":"Feng","middleName":"","lastName":"Wu","suffix":""},{"id":7280538,"identity":"febbdc9d-f32b-42ba-9df4-51b3f3618444","order_by":3,"name":"Wenjing Liu","email":"","orcid":"","institution":"Guangxi Normal University","correspondingAuthor":false,"prefix":"","firstName":"Wenjing","middleName":"","lastName":"Liu","suffix":""},{"id":7280539,"identity":"c5775a95-4a2f-40e0-9864-403489cd78dd","order_by":4,"name":"Yiqiu Ning","email":"","orcid":"","institution":"Guangxi Normal University","correspondingAuthor":false,"prefix":"","firstName":"Yiqiu","middleName":"","lastName":"Ning","suffix":""},{"id":7280540,"identity":"59df13ef-1526-4d9c-8fc5-d957866329db","order_by":5,"name":"Zhenrong Huang","email":"","orcid":"","institution":"Guangxi Normal University","correspondingAuthor":false,"prefix":"","firstName":"Zhenrong","middleName":"","lastName":"Huang","suffix":""},{"id":7280541,"identity":"9ef828aa-dbe3-41ff-b3c8-ed03b22f90c9","order_by":6,"name":"Yu Liang","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA0klEQVRIiWNgGAWjYBACA2YILQflMxOvxZgELVA6sYF4LezMzx5+bTucPj8i95gEQ4V1YgP72QMEHMZmbizbdjh345lzaRIMZ9ITG3jyEghoYTCTlgRpae8xk2BsO5zYIMFjQEAL+zeQlnTDZh6gln9EaeExk/zYdjhBnh1kSwNxWsqkGc6lG27gOWNskXAs3biNJwe/Fvv+49skf5RZy8vPyDG88aHGWraf/Qx+LSDAzMsGtO4AkJUAxGwE1QMB448/DAzyDcQoHQWjYBSMghEJAEFOPXIdprtEAAAAAElFTkSuQmCC","orcid":"","institution":"Guangxi Normal University","correspondingAuthor":true,"prefix":"","firstName":"Yu","middleName":"","lastName":"Liang","suffix":""}],"badges":[],"createdAt":"2020-12-24 14:44:06","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-135639/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-135639/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":4824351,"identity":"24cc7763-e403-4830-b474-71a9aa18e79e","added_by":"auto","created_at":"2021-01-08 22:05:05","extension":"jpg","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":2104131,"visible":true,"origin":"","legend":"Phenotype of Camellias in high calcium stress. A: Camellia nitidissima plants, bar=25 cm; B: Camellia limonia plants, bar= 25 cm; Leaves of Camellia nitidissima (lower leaves) and Camellia limonia (upper leaves), bar=15 cm.","description":"","filename":"Figure1.jpg","url":"https://assets-eu.researchsquare.com/files/rs-135639/v1/c3cf627feaf9ea1535dec355.jpg"},{"id":4824554,"identity":"9976d538-9971-4eb8-b6f0-7ed6b8f6ebea","added_by":"auto","created_at":"2021-01-08 22:11:06","extension":"jpg","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":505670,"visible":true,"origin":"","legend":"Microscopic observation of Camellias after high calcium treatment. A-D: Camellia nitidissima plants. A and B: SEM observation of Camellia nitidissima leaves(A) and stoma(B), C: paraffin section of Camellia nitidissima leaf, bar= 100 μm; D: paraffin section of Camellia nitidissima root, bar= 100 μm; E-H: Camellia limonia plants. E and F: SEM observation of Camellia limonia leaves(E) and stoma(F), G: paraffin section of Camellia limonia leaf, bar= 100 μm; H: paraffin section of Camellia limonia root, bar= 100 μm.","description":"","filename":"Figure2.jpg","url":"https://assets-eu.researchsquare.com/files/rs-135639/v1/e751dafee0b156fd20bfab39.jpg"},{"id":4824440,"identity":"24a50123-6407-4f17-bda8-445d2c8dc52f","added_by":"auto","created_at":"2021-01-08 22:08:06","extension":"jpg","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":385771,"visible":true,"origin":"","legend":"Morphological and physiological measurement of plants. A: transpiration rate of leaves. B: stomatal conductance of leaves. C: the net photosynthetic rate of leaves. D: leaf area. E: chlorophyll content. F: nitrogen content of leaves. G: leaf temperature. H: moisture on the leaf. All the results are represented as the mean ± standard deviation (STD; n = 3). Statistically significant differences were determined using a two-tailed paired Student’s t-test compared with WT plants under similar conditions, and the results are indicated by **P \u003c 0.01 and *P \u003c 0.05.","description":"","filename":"Figure3.jpg","url":"https://assets-eu.researchsquare.com/files/rs-135639/v1/5d2ea4df092187966c8b458f.jpg"},{"id":4824363,"identity":"70518db5-fda7-40bd-a501-bd4ad342610e","added_by":"auto","created_at":"2021-01-08 22:05:06","extension":"jpg","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":1483966,"visible":true,"origin":"","legend":"Hierarchical cluster analysis of metabolites. The relative content in the figure is shown by the different colors. The columns represent samples and the rows represent metabolites.","description":"","filename":"Figure4.jpg","url":"https://assets-eu.researchsquare.com/files/rs-135639/v1/38b195ab820ee89eec671353.jpg"},{"id":4824360,"identity":"15470f5f-915d-4eaa-b1ea-1de9bdb82f19","added_by":"auto","created_at":"2021-01-08 22:05:06","extension":"jpg","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":1625764,"visible":true,"origin":"","legend":"Differential metabolite-associated thermogram. Positive correlation tends to 1, and negative correlation tends to-1. At the same time, the cor.test function in R (v3.1.3) is used for statistical test of the metabolite association analysis, and there is a significant correlation with p-value ≤ 0.05.","description":"","filename":"Figure5.jpg","url":"https://assets-eu.researchsquare.com/files/rs-135639/v1/8b31644990d57c627db1fcc3.jpg"},{"id":4824356,"identity":"0451d6bf-83db-4381-b7fd-3b231c323142","added_by":"auto","created_at":"2021-01-08 22:05:06","extension":"jpg","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":627668,"visible":true,"origin":"","legend":"Genus level and species level abundance (A and B). The horizontal axis represents samples and the vertical axis represents species. The color ranges from red to blue and the abundance ranges from low to high. The cluster diagram can only be made when the sample number is ≥3.","description":"","filename":"Figure6.jpg","url":"https://assets-eu.researchsquare.com/files/rs-135639/v1/b33f5ad00bf0d7754a9f553a.jpg"},{"id":4824365,"identity":"8636a44a-91d7-4414-a9ee-65d034f44ad2","added_by":"auto","created_at":"2021-01-08 22:05:06","extension":"jpg","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":785517,"visible":true,"origin":"","legend":"Pathways flavone and flavonol biosynthesis and flavonoid biosynthesis. Pathways flavone and flavonol biosynthesis and Flavonoid biosynthesis are vital. The figure showed their relationship.","description":"","filename":"Figure7.jpg","url":"https://assets-eu.researchsquare.com/files/rs-135639/v1/f8f20849178b92da0376cdb2.jpg"},{"id":4824556,"identity":"6493acc2-7dce-4c4f-87b7-260165365529","added_by":"auto","created_at":"2021-01-08 22:11:06","extension":"jpg","order_by":8,"title":"Figure 8","display":"","copyAsset":false,"role":"figure","size":442150,"visible":true,"origin":"","legend":"Integrative pathway network. The important pathways connect to each other by key metabolites. Genes regulate the production of metabolites by controlling enzyme synthesis.","description":"","filename":"Figure8.jpg","url":"https://assets-eu.researchsquare.com/files/rs-135639/v1/3b8b3a6450ead2a56e690491.jpg"},{"id":15670471,"identity":"337a4b33-cade-489d-8354-8538a9cd5389","added_by":"auto","created_at":"2021-11-18 13:59:34","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":1383672,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-135639/v1/acbb6cb2-01d9-48da-938e-d8f4ec329573.pdf"},{"id":4824439,"identity":"d385583a-3977-4ba1-9773-b6cd58cd6002","added_by":"auto","created_at":"2021-01-08 22:08:05","extension":"jpg","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":268749,"visible":true,"origin":"","legend":"","description":"","filename":"FigureS1.jpg","url":"https://assets-eu.researchsquare.com/files/rs-135639/v1/8e803e34145a796c032eb485.jpg"},{"id":4824353,"identity":"3382ba4d-e491-43af-859f-a3f7779e3c28","added_by":"auto","created_at":"2021-01-08 22:05:06","extension":"jpg","order_by":2,"title":"","display":"","copyAsset":false,"role":"supplement","size":1033811,"visible":true,"origin":"","legend":"","description":"","filename":"FigureS2.jpg","url":"https://assets-eu.researchsquare.com/files/rs-135639/v1/3214742ef8c72d791acafce7.jpg"},{"id":4824441,"identity":"618bfa63-f7fa-4c39-ac1b-ab46532375d2","added_by":"auto","created_at":"2021-01-08 22:08:06","extension":"jpg","order_by":3,"title":"","display":"","copyAsset":false,"role":"supplement","size":401966,"visible":true,"origin":"","legend":"","description":"","filename":"FigureS3.jpg","url":"https://assets-eu.researchsquare.com/files/rs-135639/v1/59da0e3589a0eed3840a1e3e.jpg"},{"id":4824362,"identity":"b655ba6c-4011-40cc-a7e6-7a1d3b1ddbf0","added_by":"auto","created_at":"2021-01-08 22:05:06","extension":"jpg","order_by":4,"title":"","display":"","copyAsset":false,"role":"supplement","size":390649,"visible":true,"origin":"","legend":"","description":"","filename":"FigureS4.jpg","url":"https://assets-eu.researchsquare.com/files/rs-135639/v1/192db59c210d559bb6e3c519.jpg"},{"id":4824445,"identity":"230c9e5a-0753-4ea7-9aab-6efc09dea4df","added_by":"auto","created_at":"2021-01-08 22:08:06","extension":"xlsx","order_by":5,"title":"","display":"","copyAsset":false,"role":"supplement","size":33613,"visible":true,"origin":"","legend":"","description":"","filename":"TableS1.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-135639/v1/95f5362389b603c5ed5d4832.xlsx"},{"id":4824555,"identity":"abe936f5-fed9-4ca1-ba27-219dfe8b7651","added_by":"auto","created_at":"2021-01-08 22:11:06","extension":"xlsx","order_by":6,"title":"","display":"","copyAsset":false,"role":"supplement","size":15570,"visible":true,"origin":"","legend":"","description":"","filename":"TableS2.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-135639/v1/a06b9ff58ea5f376e29aad1f.xlsx"},{"id":4824443,"identity":"3c2e2789-10ba-496c-921e-8f3157aa7677","added_by":"auto","created_at":"2021-01-08 22:08:06","extension":"docx","order_by":7,"title":"","display":"","copyAsset":false,"role":"supplement","size":16832,"visible":true,"origin":"","legend":"","description":"","filename":"Supplementaryinformationtitlepagewithlegends.docx","url":"https://assets-eu.researchsquare.com/files/rs-135639/v1/8baf622e5d1023bb70f2a72b.docx"}],"financialInterests":"","formattedTitle":"\u003cp\u003eExploring the adaptive mechanism of \u003cem\u003eCamellia limonia\u003c/em\u003e in karst high calcium environment via an integrative analysis of metabolomics and metagenomics\u003c/p\u003e","fulltext":[{"header":"Introduction","content":"\u003cp\u003e\u003cem\u003eCamellia\u003c/em\u003e is originated from China and have golden flowers, glossy leaves and beautiful appearance with transparent waxy surface [1]. It is a kind of rare plant that has not only high ornamental and medicinal value but also is very useful genetic resource.\u003cem\u003e Camellia \u003c/em\u003eis mainly distributed in the southern China and the north of Vietnam [2]. \u003cem\u003eCamellia \u003c/em\u003especies are heat-loving tropical plants, belonging to the tropical tree species, which have a high requirement for heat and mainly grow under shady and moist evergreen broad-leaf forests or along streamside in the mountains with the altitudes that is less than 700 meters [3].\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eCamellia limonia\u003c/em\u003e are an evergreen shrub with 1-2m in height. And their elliptic leathery leaves are 5-10cm in length and 2-3.5cm in width with in margin serrate. Their base is broadly wedge-shaped with dark brown glandular spots below. Flowers solitarily have 8 lemon yellow petals [1].\u003cem\u003e Camellia limonia\u003c/em\u003e mainly grow in karst rock mountain area, whereas \u003cem\u003eCamellia nitidissima\u003c/em\u003e grow in normal soil mountain. Karst landscape is a typical ecosystem with very fragile ecological environment, mainly formed by carbonate dissolution. Karst regions cover about 20% of the dry ice-free land of the earth [4]. It is approximately 22 million square kilometers, accounting for about 12% of the total land area in the world[5]. China is one of the countries with the extensive karst landscapes and is famous for its diverse karst landforms[6]. Chinese karst landscape covers an area of 3.443 million square kilometers, accounting for one third of the country's land area. In southwest China, karst landforms cover an area of about 540,000 square kilometers, which is the largest and most representative karst fragile environment[6]. Located in the humid subtropical monsoon climate zone, the rainfall is abundant, but the precipitation has seasonal characteristics [7]. However, due to the obvious seasonal rainfall and dual hydrogeological structure, water in the karst area is formed only during the rainy season. As a result, karst soil water storage capacity is very limited, and has significant seasonal variation, which result in the shortage of water resources [4]. Because of the development of rock crevasses and falling-water caves, the situation of soil erosion becomes more serious. In these regions, the soil is shallow and discontinuous. Vegetation grows slowly with low biomass [5].\u003c/p\u003e\n\u003cp\u003eKarst landforms are mainly distributed in monsoon and valley rain forests, which is beneficial for the growth of \u003cem\u003eCamellia\u003c/em\u003e. The soil types, ashes and elemental components (especially calcium) content of the \u003cem\u003eCamellia \u003c/em\u003ein the two landform types are significantly different. Studies have shown that the leaf gray matter and calcium contents of \u003cem\u003eCamellia\u003c/em\u003e in karst landforms are higher than the \u003cem\u003eCamellia\u003c/em\u003e that grows in normal erosion landforms [8]. Karst is a very special environment whose shape is mainly the typical topography of the surface and underground part is formed by dissolving water on soluble rocks, which is the most typical topography in karst landforms [9]. Moreover, the soil in the karst areas is chiefly formed by weathering of bedrock carbonate rock, which causes the soil in the karst area having higher calcium content. Plants in karst areas have been adapted to grow in a drought and high calcium environment for a long time through their own physiological processes and structural adjustments in the evolution. Soil moisture and calcium carbonate content have the greatest impact on vegetation spectral characteristics in karst areas [10].\u003c/p\u003e\n\u003cp\u003eIn the present research, we focused on the adaptation mechanism and physiological changes of \u003cem\u003eCamellia limonia\u003c/em\u003e in the karst high calcium environment. In order to systematically obtain the regulation mechanism of \u003cem\u003eCamellia limonia\u003c/em\u003e in the karst high calcium environment and find key pathway and important metabolites, integrative analysis of metabonomics and metagenomics was conducted and \u003cem\u003eCamellia nitidissima \u003c/em\u003ewas used as comparison. It is shown that flavone and flavonol biosynthesis is the main regulation pathway, which involves the synthesis of flavonoids. Flavonoids further play roles in functions in drought and salt stress, which might contribute to the karst adaptation of \u003cem\u003eCamellia limonia\u003c/em\u003e. This study can improve the research about adaptive mechanism of plants in karst environments. Also, a multi omics study of \u003cem\u003eCamellia\u003c/em\u003e can provide a practical foundation and valuable data resources for endangered plant conservation research.\u003c/p\u003e"},{"header":"Methods","content":"\u003cp\u003e\u003cstrong\u003ePlant material and high calcium stress treatment\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eCamellia limonia\u003c/em\u003e and \u003cem\u003eCamellia nitidissima\u003c/em\u003e used in this study were collected from Guangxi Institute of Botany. All the plants grew in outdoor greenhouse (located in Guangxi Normal University, Guilin, China) for shading cultivation. Ten plants of each variety with the same phenotype were selected as an experimental group. Each treatment plant was irrigated with 300mM calcium chloride solution every two days, while the control group was irrigated with pure water. Leaves were selected for metabolomics. Soil ponds with roots were selected as metagenome. The samples were collected on the 15th day after treatment when the phenotype of plants showed significant differences between \u003cem\u003eCamellia limonia\u003c/em\u003e and \u003cem\u003eCamellia nitidissima\u003c/em\u003e.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eMorphological measurement \u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe morphological data of plants in high calcium treatment were collected every 5 hours from 7:00 to 18:00 every day. For each plant, at least 3 healthy leaves with the same leaf size and age were randomly selected to measure the photosynthetic parameters, leaf area, temperature, humidity, nitrogen content and chlorophyll content of leaves, respectively. Photosynthetic parameters were measured by Li-6400 photosynthetic apparatus [11]. Plant morphological meter TYS-4N was used to measure nitrogen content, temperature and humidity of leaves. The parameters of leaf length, width and leaf area of plant leaves were measured by a leaf area meter YMJ-B. By repeated freezing and thawing extraction to extract the chlorophyll of leaves were extracted and measured by the methods mentioned in previous relative references [12]. The \u003cem\u003eCamellia nitidissima\u003c/em\u003e and \u003cem\u003eCamellia limonia \u003c/em\u003eleaves after were selected for scanning electron microscope (SEM) and paraffin sections conservation based on the previous study [13].\u003c/p\u003e\n\u003cp\u003eSoil moisture, soil temperature and chemical element content were measured daily, including nitrogen (N), phosphorus (P), and potassium (K). The measurements were taken every five hours from 7:00 to 18:00 every day after high calcium treatment. The contents of N, P and K were measured by extracting soil of plants roots [14]. The PH, temperature and humidity of soil were measured by SOIL SENSOR.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eMetabolomic analysis\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe leaves under calcium treatment for 15 days were used to conduct metabolomic analysis. All the samples were stored at -80℃ refrigerator and six biological repeat for \u003cem\u003eCamellia limonia\u003c/em\u003e and \u003cem\u003eCamellia nitidissima\u003c/em\u003e, respectively.\u003c/p\u003e\n\u003cp\u003eMultivariate analysis technologies including Principal Component Analysis (PCA) and Partial Least Squares (PLS) regressions were used at present study to identify several uncorrelated variables (principal components, PCS) from a large set of data [15].\u003c/p\u003e\n\u003cp\u003eAdditionally, methodological advances in untargeted metabolomics and data analysis of tandem mass-spectrometric (MS/MS) data have enabled comprehensive analysis of molecular components in metabolic extracts [16]. These MS/MS spectra can be viewed as barcodes or fingerprints of metabolites and some software tools have been developed to take advantage of this structural information [17].\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eMetagenomics analysis \u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eIn order to study the microbial community composition and different express genes in the soil microorganism, we selected plants with obvious different phenotype traits after high calcium treatment to take the soil samples. The soil with the selected plant was taken with a small amount of roots from each plant. Then the samples were stored at -80℃ refrigerator until sequencing and two biological repeat for \u003cem\u003eCamellia limonia\u003c/em\u003e and \u003cem\u003eCamellia nitidissima\u003c/em\u003e, respectively [18].\u003c/p\u003e\n\u003cp\u003eIllumina PE sequencing strategy was used to extract total DNA from the soil and then agarose gel electrophoresis was used to detect the quality of the samples and finally the library was obtained. In addition, the Velvet software was used for sequence assembly and finally functional annotation and comparison analysis were performed according to the literature [19, 20].\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eData processing and statistics analysis\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eData analysis was conducted by using the SPSS 20.0 program [4]. Independent-sample was used to detect the indexes of morphology and soil. Leaf morphology indicators include: net photosynthetic rate, stomatal conductance, transpirational rate, leaf area, chlorophyll, nitrogen content, leaf surface moisture and leaf surface temperature; Soil morphology indicators include: pH, temperature and humidity. The methods of mathematical statistics are used to calculate the error value and the significant difference. The results were taken into account as not significant at P \u0026gt; 0.05, significant at 0.01\u0026lt;P\u0026lt;0.05, and extremely significant at P\u0026lt;0.01. PCA and PLS are used to discriminate and analyze the error value and significant difference. A cluster analysis of the identified metabolites was conducted after auto-scaling (mean-centered and divided by the standard deviation) by using MultiExperiment Viewer [21].\u003c/p\u003e"},{"header":"Results","content":"\u003cp\u003e\u003cstrong\u003eE\u003c/strong\u003e\u003cstrong\u003eff\u003c/strong\u003e\u003cstrong\u003eects of calcium on the karst adaptive characteristics of Camellia\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAbiotic stress signals can cause fluctuations of cytoplasmic calcium ions and changes in phosphorylation forms of some proteins. Calcium ions mediated by calcium channels, change the concentration of calcium ions in cytoplasm and amplify the stress signals with intracellular calcium signal receptors [22]. However, high calcium stress can destroy the cytoskeleton and biofilm structure, resulting in intracellular metabolic imbalance.\u003c/p\u003e\n\u003cp\u003eIn order to study the adaptability of the karst environment with high calcium content of \u003cem\u003eCamellia limonia\u003c/em\u003e, we first processed \u003cem\u003eCamellia limonia \u003c/em\u003ewith high calcium content and take \u003cem\u003eCamellia nitidissima \u003c/em\u003eas the control, the results showed that the adaptability of \u003cem\u003eCamellia limonia\u003c/em\u003e to the environment with high calcium content was stronger.\u003c/p\u003e\n\u003cp\u003eUnder the calcium stress, the phenotypes of \u003cem\u003eCamellia nitidissima\u003c/em\u003e and \u003cem\u003eCamellia limonia\u003c/em\u003e became different. The leaves of \u003cem\u003eCamellia nitidissima\u003c/em\u003e turned yellow and wilting, while the leaves of \u003cem\u003eCamellia limonia\u003c/em\u003e kept relatively healthy (Figure 1). The SEM and microscopic observation also showed that the leaves and roots of \u003cem\u003eCamellia limonia\u003c/em\u003e are less wrinkled and dehydrated than \u003cem\u003eCamellia nitidissima \u003c/em\u003e(Figure 2). The stomatal conductance and transpiration rate of the two kinds of \u003cem\u003eCamellia\u003c/em\u003e showed a downward trend. The stomatal conductance of \u003cem\u003eCamellia nitidissima\u003c/em\u003e decreased more than the one of \u003cem\u003eCamellia limonia\u003c/em\u003e. While the reduce degree of transpiration rate \u003cem\u003eCamellia limonia\u003c/em\u003e is less than in \u003cem\u003eCamellia nitidissima\u003c/em\u003e. The chlorophyll content of \u003cem\u003eCamellia limonia \u003c/em\u003eincreased by 53.12% and \u003cem\u003eCamellia nitidissima \u003c/em\u003edecreased by 56.91%. The net photosynthetic rate of \u003cem\u003eCamellia limonia \u003c/em\u003ewas more stable than \u003cem\u003eCamellia nitidissima\u003c/em\u003e. The average net photosynthetic rate of \u003cem\u003eCamellia limonia \u003c/em\u003eand \u003cem\u003eCamellia nitidissima \u003c/em\u003ewas 3.0853 \u0026mu;mol CO\u003csub\u003e2 \u003c/sub\u003em\u003csup\u003e-2\u003c/sup\u003es\u003csup\u003e-1\u003c/sup\u003e and 2.4456 \u0026mu;mol CO\u003csub\u003e2 \u003c/sub\u003em\u003csup\u003e-2\u003c/sup\u003es\u003csup\u003e-1\u003c/sup\u003e, respectively. Stomatal conductance affects the oxygen and carbon dioxide, indicating that the water might be the more important factor influencing the net photosynthetic rate when gas were equivalently affected [23]. According to the data, the average leaf humidity of \u003cem\u003eCamellia limonia \u003c/em\u003eis 21.29% RH, \u003cem\u003eCamellia nitidissima \u003c/em\u003eleaf humidity averages of 17.88% RH. The average transpiration rate and leaf temperature of \u003cem\u003eCamellia nitidissima \u003c/em\u003ewere 0.64 mmolH\u003csub\u003e2\u003c/sub\u003eOm\u003csup\u003e-2\u003c/sup\u003es\u003csup\u003e-1\u003c/sup\u003e and 31.94 ℃，\u003cem\u003eCamellia nitidissima \u003c/em\u003ewere 0.92 mmolH\u003csub\u003e2\u003c/sub\u003eOm\u003csup\u003e-2\u003c/sup\u003es\u003csup\u003e-1\u003c/sup\u003e and 30.95 ℃, suggesting the water content of \u003cem\u003eCamellia limonia \u003c/em\u003eis more than \u003cem\u003eCamellia nitidissima\u003c/em\u003e (Figure 3). According to the results above, \u003cem\u003eCamellia nitidissima \u003c/em\u003ehas relatively stronger adaptability to the environment with high calcium.\u003c/p\u003e\n\u003cp\u003eThe soil moisture and temperature of the two kinds of \u003cem\u003eCamellia \u003c/em\u003edid not obviously change. The pH value of the soil of \u003cem\u003eCamellia limonia \u003c/em\u003edecreased slightly. In addition, the nitrogen and phosphorus content in the soil of \u003cem\u003eCamellia limonia \u003c/em\u003eincreased by 58.60% and 53.99%, respectively, while the nitrogen and phosphorus content in the soil of \u003cem\u003eCamellia nitidissima \u003c/em\u003edecreased by 13.01% and increased by 49.63%, respectively (Figure S1).\u003c/p\u003e\n\u003cp\u003eAvailable nitrogen and phosphorus are the main nutrient elements that restrict plant growth, while plant growth is mainly restricted by nitrogen elements under low ratio of nitrogen and phosphorus [24]. Under high calcium stress, nitrogen content of \u003cem\u003eCamellia limonia \u003c/em\u003eleaf increased by 14.44% and the available nitrogen in soil of \u003cem\u003eCamellia limonia\u003c/em\u003e is increased by 58.60%, while nitrogen content of \u003cem\u003eCamellia nitidissima \u003c/em\u003eleaf decreased by 28.57% and the available nitrogen in soil of \u003cem\u003eCamellia nitidissima\u003c/em\u003e decreased by 13.01%, suggesting \u003cem\u003eCamellia limonia\u003c/em\u003e has better nitrogen use efficiency (Figure 3), which may also help \u003cem\u003eCamellia limonia\u003c/em\u003e survived in high calcium environment.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eMetabonomics of Camellia leaf in karst high calcium environment\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eIn order to analysis the different metabolites of \u003cem\u003eCamellia\u003c/em\u003e leaf in karst high calcium environment, the metabonomics were conducted. The PCA score chart showed the degree of aggregation and dispersion of\u003cem\u003e Camellia limonia\u003c/em\u003e and \u003cem\u003eCamellia nitidissima\u003c/em\u003e (Figure S2). The distribution points of \u003cem\u003eCamellia limonia\u003c/em\u003e are relatively close, indicating that the composition and concentration of the variables/molecules contained in \u003cem\u003eCamellia limonia\u003c/em\u003e are closer.\u003c/p\u003e\n\u003cp\u003eA total of 9133 metabolites were identified through mathematical statistical analysis and screening, of which 6422 were in the positive ion mode and 2711 were in the negative ion mode (Figure S3). 8966 differential metabolites discemed include metabolites of \u003cem\u003eCamellia limonia \u003c/em\u003eand\u003cem\u003e Camellia nitidissima\u003c/em\u003e. The accurate molecular weight of the metabolites (molecular weight error \u0026lt;15ppm) were confirmed and the fragment information were divided according to the MS/MS mode. We identified 220 different metabolites in HMDB (Human Metabolome Database, HMDB), metlin, massbank, lipiduaps, mzdound databases and obtained the statistical analysis results (Table S1). Under the high calcium stress, \u003cem\u003eCamellia limonia\u003c/em\u003e has 53 metabolites up-regulated and 167 metabolites down-regulated, of which the main up-regulated metabolite is kaempferol-3-o-rutinoside and tyrosol. \u003cem\u003eCamellia nitidissima\u003c/em\u003e has 167 up-regulated metabolites and 53 down-regulated metabolites, of which the main metabolites up-regulated are n-acetylhistidine, myristic acid, and pelargonin (Figure 4). The five different metabolites mentioned above have the largest significant difference. The content of these metabolites were changed to affect the metabolic process of \u003cem\u003eCamellia\u003c/em\u003e, indicating that the significantly different metabolites mainly affect the growth of\u003cem\u003e Camellia limonia\u003c/em\u003e and \u003cem\u003eCamellia nitidissima\u003c/em\u003e after the treatment of high calcium stress.\u003c/p\u003e\n\u003cp\u003eAmong them, kaempferol-3-o-rutinoside, tyrosol, 6-o-methyldeacetylisoipecoside,n-acetylhistidine, myristic acid and pelargonin are positively correlated. But the three metabolites of \u003cem\u003eCamellia limonia\u003c/em\u003e were negatively correlated with the three metabolites of \u003cem\u003eCamellia nitidissima\u003c/em\u003e (Figure 5). It is inferred that the up-regulated metabolites of \u003cem\u003eCamellia limonia\u003c/em\u003e had the opposite trend to the up-regulated metabolites of \u003cem\u003eCamellia nitidissima\u003c/em\u003e, indicating that \u003cem\u003eCamellia limonia \u003c/em\u003ehave different metabolite trends in high calcium stress.\u003c/p\u003e\n\u003cp\u003eMetPA is a part of metaboanalyst (www.metaboanalyst.ca) and is mainly based on the KEGG metabolic pathway. The MetPA related metabolic and KEGG pathway analysis show 42 main metabolic pathways. The most important is biotin metabolism, tyrosine metabolism, flavone and flavonol biosynthesis, inositol phosphate metabolism, alpha-linolenic acid metabolism, ubiquinone and other terpenoid-quinone biosynthesis, phenylpropanoid biosynthesis, lysine biosynthesis, phenylalanine metabolism, glyoxylate and dicarboxylate metabolism, pantothenate and CoA biosynthesis, arginine and proline metabolism and other metabolic (Table S2). The pathways with most significant differences are tyrosine metabolism, flavone and flavonol biosynthesis, phenylpropanoid biosynthesis and biotin metabolism. And more differential metabolites are enriched in these four pathways, indicating that these pathways are more sensitive to high calcium stress for \u003cem\u003eCamellia limonia\u003c/em\u003e.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eMetagenomics of soil in karst high calcium environment\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eUnder biotic and abiotic stresses, plants can actively change their own physiological conditions and exudates of their roots. That also has effects on the related microbial communities. Consequently, this directly or indirectly affects the health of plants. The initial central microorganism and exudates of plant root together form a microbial community, including plant central flora, core flora and other microbial communities. Compared with the plants under high calcium treatment, the root microbial populations of the plants in normal condition are different. Therefore, these different populations could affect the growth of plants. Metagenomics analysis of the soil of \u003cem\u003eCamellia limonia\u003c/em\u003e and \u003cem\u003eCamellia nitidissima\u003c/em\u003e root were conducted to study the effect of microbial communities under high calcium stress. The result showed that tyrosine metabolism, galactose metabolism, inositol phosphate metabolism, biotin metabolism, phenylpropanoid biosynthesis, arginine and proline metabolism were significantly affected. The microorganism regulated these metabolic pathways. At the same time, they regulate the expression of responsive proteins through responsive genes. Finally, these microorganisms influnce the metabolic pathways of \u003cem\u003eCamellia\u003c/em\u003e.\u003c/p\u003e\n\u003cp\u003eThe software Velvet was used to evaluate the assembly results, including the statistical table of assembly results, the distribution chart of the cumulative length of contigs and the distribution map of GC content (Table 1).\u003c/p\u003e\n\u003cp\u003e\u003cimg src=\"https://myfiles.space/user_files/58653_1b1c6aeb34a62c68/58653_custom_files/img1609446962.jpg\"\u003e\u003c/p\u003e\n\u003cp\u003eWe examined nearly 20,000 genes from each sample of rhizosphere microbes. And we found that reduced pathogenicity and loss of toxicity were the most important factors. In order to get a more detailed classification result of species, we match the sequence with the bacterial reference sequence or the virus reference sequence. Moreover, the MetaPhlAn was used to comment the species classification (Table 2). In all samples, rubrivivax is the main member. In root soil of \u003cem\u003eCamellia limonia\u003c/em\u003e, rhodanobacter, intrasporangium, microbacterium, ilumatobacter are taking disadvantages. Different from \u003cem\u003eCamellia limonia\u003c/em\u003e, the \u003cem\u003eCandidatus_Solibacter\u003c/em\u003e was the dominant species of \u003cem\u003eCamellia nitidissima\u003c/em\u003e (Figure 6). Further, in order to reveal the distribution and evolution relationship of species, phylogenetic analysis was conducted (Figure S4). The three colors represent the most important four species (\u003cem\u003edehalogenimonas\u003c/em\u003e\u003cem\u003e lykanthroporepellens\u003c/em\u003e, \u003cem\u003emicavibrio aeruginosavorus\u003c/em\u003e, \u003cem\u003epolymorphum gilvum\u003c/em\u003e and \u003cem\u003ecandidatus accumulibacter phosphatis\u003c/em\u003e) and their evolution relationship. The above results about the metagenomic analysis indicated that, compared with \u003cem\u003eCamellia nitidissima\u003c/em\u003e, the types of soil microorganisms of \u003cem\u003eCamellia limonia\u003c/em\u003e under high calcium stress were different. And these microorganisms are related to the salt tolerance and drought resistance of plants.\u003c/p\u003e\n\u003cp\u003e\u003cimg src=\"https://myfiles.space/user_files/58653_1b1c6aeb34a62c68/58653_custom_files/img1609447042.jpg\"\u003e\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003e\u003cem\u003eCamellia limonia\u003c/em\u003e\u003c/strong\u003e\u003cstrong\u003e adaptive regulation to karst\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eIn order to study the karst adaptation of \u003cem\u003eCamellia limonia\u003c/em\u003e, we conducted the integrative analysis of metabolomics and metagenomics. The result showed that flavone and flavonol biosynthesis pathway which is one of the most affected pathways can be regulated by apigenin, kaempferol, astragalin, isoquercitrin metabolites and \u003cem\u003eTT7, UGT78D1, UGT78D2\u003c/em\u003e genes. This metabolic pathway involves the synthesis of flavonoids. Flavonoids have functions in drought and salt resistance [25], so flavonoids may also contribute to the karst environment adaptation for \u003cem\u003eCamellia limonia\u003c/em\u003e. Flavonoids include kaempferol, rutin, hesperidin, anthocyanin, etc. Among them, kaempferol has the largest difference in content between\u003cem\u003e Camellia limonia\u003c/em\u003e and \u003cem\u003eCamellia nitidissima\u003c/em\u003e , and the content of \u003cem\u003eCamellia limonia\u003c/em\u003e is greater than \u003cem\u003eCamellia nitidissima\u003c/em\u003e. Therefore, \u003cem\u003eCamellia limonia\u003c/em\u003e grow better than \u003cem\u003eCamellia nitidissima\u003c/em\u003e under the same calcium stress, which indicated that kaempferol may play a critical role in the high calcium tolerance.\u003c/p\u003e\n\u003cp\u003eThe citrate cycle (TCA cycle, Krebs cycle) is an important aerobic pathway for the final steps of the oxidation of carbohydrates and fatty acids. The citrate cycle is also an important pathway for differential metabolites. In the relative pathway of TCA cycle, the \u003cem\u003esdhA, frdA\u003c/em\u003e genes regulate the synthesis of succinate dehydrogenase/fumarate reductase-flavoprotein subunit. Compared with \u003cem\u003eCamellia nitidissima\u003c/em\u003e, the \u003cem\u003esdhA, frdA\u003c/em\u003e were differently expression genes in \u003cem\u003eCamellia limonia\u003c/em\u003e, which suggested that this two genes might simultaneously regulate the synthesis of fumarate and succinate and further affect the high calcium tolerance of \u003cem\u003eCamellia limonia\u003c/em\u003e.\u003c/p\u003e"},{"header":"Discussion","content":"\u003cp\u003e\u003cstrong\u003eChanges in gene regulation and metabolic material network lead to \u003cem\u003eCamellia\u003c/em\u003e karst adaptation \u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eCalcium ion is a necessary mineral element for plant growth and development [26]. It is also an important component of plant cell wall and has the function of maintaining intracellular ion balance and cell membrane stability [27]. Calcium ions in the cytoplasm are generally kept at a low concentration, and high concentration of calcium ions will disturb the physiological functions of plant cells [28]. The absorption of calcium by plants is positively correlated with the concentration of calcium ions in soil. If the concentration of calcium ions in soil is too high, the plants will suffer the high calcium stress. High concentration of calcium ions can seriously affect plant photosynthesis and inhibit plant development [29]. From the microscopic point of view, calcium ions function as a signal molecule of plant cells to environmental changes, but if the concentration of calcium ions in cells is too high, it will lead to cell dysfunction [30].\u003c/p\u003e\n\u003cp\u003eThe phenotypic characteristics of \u003cem\u003eCamellia nitidissima\u003c/em\u003e were significantly different from \u003cem\u003eCamellia nitidissima\u003c/em\u003e after high calcium treatment. The leaves of \u003cem\u003eCamellia nitidissima\u003c/em\u003e are easier to fall off, which is probably caused by increased abscisic acid in plants after high calcium treatment [31]. Moreover, the color of the \u003cem\u003eCamellia\u003c/em\u003e leaves has changed obviously. Compared with \u003cem\u003eCamellia nitidissima\u003c/em\u003e, leaf temperature, soil moisture and potassium contents of \u003cem\u003eCamellia limonia\u003c/em\u003e are basically maintained at the same level after high calcium treatment. Chlorophyll content, phosphorus content, nitrogen content and leaf moisture of \u003cem\u003eCamellia limonia\u003c/em\u003e are generally higher than\u003cem\u003e Camellia nitidissima\u003c/em\u003e. While the pH of soil, transpiration rate of \u003cem\u003eCamellia limonia\u003c/em\u003e are relatively lower than \u003cem\u003eCamellia nitidissima\u003c/em\u003e. The Camellia plants lost much water which was mainly caused by higher transpiration under the calcium treatment. For \u003cem\u003eCamellia limonia\u003c/em\u003e, its transpiration effect was relatively weak, indicating that this may be one of the results of adapting to the karst high calcium environment. \u003cem\u003eCamellia limonia\u003c/em\u003e can properly undermine transpiration, which is beneficial to live in a karst high calcium environment. Transport of Ca\u003csup\u003e2+\u003c/sup\u003e to leaves depends on several factors, such as xylem sap Ca\u003csup\u003e2+\u003c/sup\u003e concentration, water absorption, transpiration and so on. Calcium accumulation may occur in organs with high transpiration. Low transpiration rates can ensure plants well hydrated, which is beneficial to absorb the extra calcium [32]. Therefore,\u003cem\u003e Camellia limonia \u003c/em\u003ecan adapt better to the high calcium karst environment under a certain transpiration intensity. While\u003cem\u003e Camellia nitidissima\u003c/em\u003e had excessive water loss in the plant body due to relatively strong transpiration, which could affect various life activities and could not survive in karst high calcium environment.\u003c/p\u003e\n\u003cp\u003eThe metabolites of \u003cem\u003eCamellia\u003c/em\u003e have undergone tremendous changes under calcium stress. Some key genes regulate the synthesis and changes of metabolites in the high calcium environment (Figure 7; Figure 8). Among the changed metabolism pathways, biotin metabolism has the highest impact factor. At the same time, flavone and flavonol biosynthesis involve many metabolic pathways and their impact factors are also high. Flavone and flavonol biosynthesis, ubiquinone and other terpenoid-quinone biosynthesis, phenylpropanoid biosynthesis and flavonoid biosynthesis play a regulatory role in metabolite changes and regulation and some of these pathways are interconnected (Figure 7). Special metabolites take part in the defense of plants with both biotic and abiotic environmental factors [33]. Glycosylation is a common modification of the structure of plant specific metabolites. It broadens the structural diversity and improves the substrate molecule solubility, stability and transportability [34]. A variety of glycosylated phenolic metabolites have been accumulated in plant tissues, such as anthocyanins and flavonoidsm [35]. Our results show that flavonoid biosynthesis and anthocyanin biosynthesis are the metabolic pathways of anthocyanins and flavonoids. This suggests that anthocyanins and flavonoids produced by these two pathways could contribute to resisting adverse environment, such as over-acidity, over-alkalinity, high calcium and drought. It is speculated that the up-regulation of these two metabolites may play the key roles for\u003cem\u003e Camellia limonia\u003c/em\u003e adaptation to the karst high calcium environment.\u003c/p\u003e\n\u003cp\u003eMoreover, the changes of different metabolites are also related to the overall metabolomic changes. Under high calcium stress, the chlorophyll content of the leaves of\u003cem\u003e Camellia nitidissima\u003c/em\u003e decreased significantly and the leaves grew wrinkled (Figure 1; Figure 2). The leaves of \u003cem\u003eCamellia limonia\u003c/em\u003e lack green luster on the surface, but their chlorophyll content remained relatively high. Plants with strong drought resistance have higher chlorophyll content and less biomass loss [36]. High calcium in the soil could cause the plants run out of water. Our result shows that the chlorophyll content of \u003cem\u003eCamellia limonia\u003c/em\u003e is higher than \u003cem\u003eCamellia nitidissima,\u003c/em\u003e which further suggested that enough chlorophyll content could ensure plant stability under high calcium condition. Under the high calcium treatment, some important pathways are very abundant, including energy metabolism associated with TCA Cycle, pathways associated with secondary metabolites, such as flavonol biosynthesis and flavonoids and other important pathways (such as biotin metabolism). All these pathways could also function in the adaptation of \u003cem\u003eCamellia limonia\u003c/em\u003e.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eResponse of \u003cem\u003eCamellia\u003c/em\u003e to karst high calcium environment and its regulation model\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eIn order to further explore the regulation mechanism, we mapped the integrative analysis pathways, including differential metabolites, genes and enzymes (Figure 8). The analysis focused on the relationship of the pathways with higher influencing factors, such as biotin, flavone and flavonol biosynthesis, flavonoid biosynthesis, tyrosine metabolism, inositol phosphate metabolism. The results showed that tyrosine, p-coumaroy-CoA and fumarate are important metabolites. They are the metabolites connecting key metabolic pathways.\u003c/p\u003e\n\u003cp\u003eL-tyrosine and fumarate connected tyrosine metabolism, phenylalanine tyrosine and tryptophan, TCA cycle, phenylpropanoid biosynthesis and other metabolic pathways. P-coumaroy-CoA played a very important role in phenylpropanoid biosynthesis. P-coumaroy-CoA is also an important hub in the regulation network. It links flavone and flavonol biosynthesis, flavonoid biosynthesis, phenylpropanoid biosynthesis. P-coumaroy-CoA was involved in the phenylpropanoid pathway. P-coumaroy-CoA mainly catalyzed the synthesis of 4-coumaroyl-CoA with the substrates of 4-coumaric acid and coenzyme A. The reaction is one of the rate-limiting steps which greatly affect the yields of lignins and flavonoids in phenylpropanoid pathway [37]. In the metabolic regulation processes, enzymes act as biocatalysts in many important bioreactions. In the phenylpropanoid biosynthesis, L-tyrosine is continuously catalyzed by 4-coumarate-CoA ligase, which is regulated by gene \u003cem\u003eACOS5\u003c/em\u003e. Shikimate O-hydroxycinnamoyl transferase is regulated by gene \u003cem\u003eSHT\u003c/em\u003e (AT2G19070). 5-O-(4-coumarayl)-D-quinate 3\u0026rsquo;-monooxygenase is regulated by gene \u003cem\u003eCYP98A3\u003c/em\u003e. Caffeate O-methyltransferase is regulated by gene \u003cem\u003eOMT1 \u003c/em\u003e(AT5G54160). All these genes mentioned above were changed in the present study. These enzymes finally lead to the synthesis of differential metabolite ferulic acid. Phenylpropanoid biosynthesis have an overlapping reaction chain. L-tyrosine can be turned to P-coumaric acid. P-coumaric acid can be transformed to P-coumaroyl-CoA catalyzed by 4-coumarate-CoA ligase (4CL). \u003cem\u003eACOS5\u003c/em\u003e and \u003cem\u003eCYP98A3\u003c/em\u003e regulated the synthesis of coumarin. Coumarin is involved in the flavonoid synthesis pathway [38]. Studies about \u003cem\u003eArabidopsis thaliana\u003c/em\u003e have shown that flavonoids and flavonoids affect the drought-resistance of plants [39]. \u003cem\u003eOMT\u003c/em\u003e (AT5G54160) and\u003cem\u003e SHT\u003c/em\u003e (AT2G19070) regulate the synthesis of lignin [40]. Studies about \u003cem\u003eArabidopsis thaliana\u003c/em\u003e have shown that lignin can improve plant tolerance in salt stress [41]. These pathways connect some key differently metabolites in the metabonomics of \u003cem\u003eCamellia\u003c/em\u003e. Genes regulate the metabolic processes of \u003cem\u003eCamellia \u003c/em\u003ewhich were also found similar function in other plants [41], suggesting these gene may play an important role in the karst adaption of \u003cem\u003eCamellia\u003c/em\u003e.\u003c/p\u003e\n\u003cp\u003ePlants in karst areas have been in the high calcium environment for a long time. And they have also responded to various environmental stresses through their own physiological processes and structural adjustments in the long-term adaptive evolution. But the regulation mechanism of \u003cem\u003eCamellia\u003c/em\u003e in karst high calcium environment still remain unclear. The present research explored the regulation mechanism of \u003cem\u003eCamellia\u003c/em\u003e in high calcium environment by integrative analysis of metabolomics and metagenomics. These omics results and analysis shown that some key metabolites and genes may affect the \u003cem\u003eCamellia limonia\u003c/em\u003e adaption of karst high calcium environment and lead to the difference of phenotype of \u003cem\u003eCamellia limonia\u003c/em\u003e and \u003cem\u003eCamellia nitidissima.\u003c/em\u003e The comprehensive omics level study on \u003cem\u003eCamellia \u003c/em\u003ecan provide practical research foundation and valuable data resources for future plant adaptability research in karst environments and conservation research of rare and endangered plants.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eEthics approval and consent to participate\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConsent for publication\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAvailability of data and materials\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe datasets generated and analyzed during the present study are available from the corresponding author on reasonable request.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors declare that they have no competing interests.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthors' contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eY.L., Y.H. and C.L. wrote the main manuscript text. F.W.prepared figures 1-3. W.L. prepared figures 4-5. Y.N. prepared figures 6 and all the tables. Z.H. prepared figure 8 and supplementary materials. All authors reviewed the manuscript.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAcknowledgements\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWe acknowledge the assistance from Xiao Wei and Jianmin Tang in Guangxi Institute of Botany, Chinese Academy of Sciences.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunding\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis work was supported by the Guangxi science and technology base and talent project(Guike AD19245098); Guangxi Natural Science Foundation(2020JJB130026); Research Project of Key Laboratory of Ecology of Rare and Endangered Species and Environmental Protection(ERESEP2020Z09); Guangxi graduate education innovation program(YCSW2020104); Innovation training program for university students(202010602251).\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eWei X, Jiang YS, Jiang SY, Qi XX, Xiong ZC, Ye WH, et al. Photosynthetic characteristics of an endangered species \u003cem\u003eCamellia nitidissima\u003c/em\u003e and its widespread congener \u003cem\u003eCamellia sinensis\u003c/em\u003e. Photosynt. 2008;46:312-314.\u003c/li\u003e\n\u003cli\u003eTang S, Bin X, Wang L, Zhong Y. Genetic diversity and population structure of yellow Camellia (\u003cem\u003eCamellia nitidissima)\u003c/em\u003e in China as revealed by RAPD and AFLP markers. Biochem Genet. 2006;44:444\u0026ndash;456.\u003c/li\u003e\n\u003cli\u003eSong L, Wang X, Zheng X, Huang D. Polyphenolic antioxidant profiles of yellow camellia. Food Chemistry. 2011;129:351\u0026ndash;357.\u003c/li\u003e\n\u003cli\u003eCao M, Wu C, Liu J, Jiang Y. Increasing leaf \u0026delta;13C values of woody plants in response to water stress induced by tunnel excavation in a karst trough valley: Implication for improving water-use efficiency. Journal of Hydrology. 2020;586:124895.\u003c/li\u003e\n\u003cli\u003eZhang Z, Huang X, Zhang J. Soil thickness and affecting factors in forestland in a karst basin in Southwest China. Trop Ecol. 2020;61:267\u0026ndash;277.\u003c/li\u003e\n\u003cli\u003eGuo F, Jiang G, Yuan D, Polk JS. Evolution of major environmental geological problems in karst areas of Southwestern China. Environ Earth Sci. 2013;69:2427\u0026ndash;2435.\u003c/li\u003e\n\u003cli\u003eZhang Z, Huang X, Zhang J. Soil thickness and affecting factors in forestland in a karst basin in Southwest China. Trop Ecol. 2020;61:267\u0026ndash;277.\u003c/li\u003e\n\u003cli\u003eSu Z. A preliminary study on the population ecology of Camellia sect. nitidissima. Guangxi Science. 1994;01:1.\u003c/li\u003e\n\u003cli\u003ePisano L, Zumpano V, Liso IS, Parise M. Geomorphological and structural characterization of the \u0026lsquo;Canale di Pirro\u0026rsquo; polje, Apulia (Southern Italy). Journal of Maps. 2020;16:479\u0026ndash;87.\u003c/li\u003e\n\u003cli\u003eXu M, Li A, Teng Y, Sun Z, Xu M. Exploring the adaptive mechanism of \u003cem\u003ePassiflora edulis\u003c/em\u003e in karst areas via an integrative analysis of nutrient elements and transcriptional profiles. BMC Plant Biol. 2019;19:185.\u003c/li\u003e\n\u003cli\u003eZhang M, Chen Y, Du L, Wu Y, Liu Z, Han L. The potential of P\u003cem\u003eaulownia fortunei\u003c/em\u003e seedlings for the phytoremediation of manganese slag amended with spent mushroom compost. Ecotoxicology and Environmental Safety. 2020;196:110538.\u003c/li\u003e\n\u003cli\u003eWun C, Rho J, Walker RobertW, Litsky W. A simplified method for the simultaneous extraction of phytoplanktonic chlorophyll and fecal sterol from water. Water Air Soil Pollut. 1979;11:1.\u003c/li\u003e\n\u003cli\u003eBuszewska-Forajta M, Patejko M, Macioszek S, Sigorski D, Iżycka-Świeszewska E, Markuszewski MJ. Paraffin-embedded tissue as a novel matrix in metabolomics study: optimization of metabolite extraction method. Chromatographia. 2019;82:1501\u0026ndash;1513.\u003c/li\u003e\n\u003cli\u003eSimeonov V, Asenov I, Diadov V. Accelerated extraction of elements from soils by ultrasonic pretreatment. Z Anal Chem. 1977;285:252.\u003c/li\u003e\n\u003cli\u003eOngo E, Falasconi M, Sberveglieri G, Antonelli A, Montevecchi G, Sberveglieri V, et al. Chemometric discrimination of philippine civet coffee using electronic nose and gas chromatography mass spectrometry. Procedia Engineering. 2012;47:977\u0026ndash;980.\u003c/li\u003e\n\u003cli\u003eShao SY, Ting Y, Wang J, Sun J, Guo XF. Characterization and identification of the major flavonoids in \u003cem\u003ePhyllostachys edulis\u003c/em\u003e leaf extract by UPLC\u0026ndash;QTOF\u0026ndash;MS/MS. Acta Chromatographica. 2020;32:228\u0026ndash;237.\u003c/li\u003e\n\u003cli\u003evan der Hooft JJJ, Mohimani H, Bauermeister A, Dorrestein PC, Duncan KR, Medema MH. Linking genomics and metabolomics to chart specialized metabolic diversity. Chem Soc Rev. 2020;49:3297\u0026ndash;3314.\u003c/li\u003e\n\u003cli\u003eL\u0026oacute;pez-Mond\u0026eacute;jar R, Tl\u0026aacute;skal V, Větrovsk\u0026yacute; T, \u0026Scaron;tursov\u0026aacute; M, Toscan R, Nunes da Rocha U, et al. Metagenomics and stable isotope probing reveal the complementary contribution of fungal and bacterial communities in the recycling of dead biomass in forest soil. Soil Biology and Biochemistry. 2020;148:107875.\u003c/li\u003e\n\u003cli\u003eHuang Y. The research of microbial genomics based on high-throughput sequencing. 2013; 12(5): 440-453.\u003c/li\u003e\n\u003cli\u003eHuang X, Liu L, Wen T, Zhu R, Zhang J, Cai Z. Illumina MiSeq investigations on the changes of microbial community in the \u003cem\u003eFusarium oxysporum\u003c/em\u003e f.sp. \u003cem\u003ecubense\u003c/em\u003e infected soil during and after reductive soil disinfestation. Microbiological Research. 2015; 181: 33-42.\u003c/li\u003e\n\u003cli\u003eYu X, Li Y, He C, Zhou J, Chen Y, Yu Z, et al. Nonvolatile metabolism in postharvest tea (\u003cem\u003eCamellia sinensis\u003c/em\u003e L.) leaves: Effects of different withering treatments on nonvolatile metabolites, gene expression levels, and enzyme activity. Food Chemistry. 2020;327:126992.\u003c/li\u003e\n\u003cli\u003eDietrich P. Physiology and biophysics of plant ligandgated ion channels. Plant Biology. 2010;12:80-93.\u003c/li\u003e\n\u003cli\u003eHazrati S, Tahmasebi-Sarvestani Z, Modarres-Sanavy S, Mokhtassi-Bidgoli A, Nicola S. Effects of water stress and light intensity on chlorophyll fluorescence parameters and pigments of \u003cem\u003eAloe vera\u003c/em\u003e L. Plant Physiology and Biochemistry. 2016;106:141\u0026ndash;148.\u003c/li\u003e\n\u003cli\u003eHartman WH, Richardson CJ. Differential Nutrient Limitation of Soil Microbial Biomass and Metabolic Quotients (qCO\u003csub\u003e2\u003c/sub\u003e): Is There a Biological Stoichiometry of Soil Microbes? PLOS ONE. 2013;8:14.\u003c/li\u003e\n\u003cli\u003eDeng J, Qin W, Yang C, Iqbal N, Takpah D, Zhang J, et al. Seed quality deterioration dynamics for isoflavones biosynthesis in soybean (\u003cem\u003eGlycine max\u003c/em\u003e L. Merr.) seeds against field mildew stress. Acta Physiol Plant. 2019;41:57.\u003c/li\u003e\n\u003cli\u003eYang Q, Wei X, Zeng X, Ye W, Yin X, Zhang W, et al. Seed biology and germination ecophysiology of \u003cem\u003eCamellia nitidissima\u003c/em\u003e. Forest Ecology and Management. 2008;255(1):113-118.\u003c/li\u003e\n\u003cli\u003eAl-Whaibi MH, Siddiqui MH, Al-Amri A, Basalah MO. Performance of faba bean under calcium and gibberellic acid application. International Journal of Plant Developmental Biology. 2010; 4(1): 60-63.\u003c/li\u003e\n\u003cli\u003eSilva LRD, Mansfield A. Where does all the calcium go? Evidence of an important regulatory role for trichomes in two calcicoles. Plant, Cell \u0026amp; Environment. 1996; 19(7): 880-886.\u003c/li\u003e\n\u003cli\u003eBroadley MR. Variation in the shoot calcium content of angiosperms. Journal of Experimental Botany. 2003;54:1431\u0026ndash;1446.\u003c/li\u003e\n\u003cli\u003eBorer CH, Hamby MN, Hutchinson LH. Plant tolerance of a high calcium environment via foliar partitioning and sequestration. Journal of Arid Environments. 2012;85:128\u0026ndash;31.\u003c/li\u003e\n\u003cli\u003eYoshida T, Fernie AR, Shinozaki K, Takahashi F. Long‐distance stress and developmental signals associated with abscisic acid signaling in environmental responses. The Plant Journal. 2020;15101.\u003c/li\u003e\n\u003cli\u003eBonomelli C, Alcalde C, Aguilera C, Videla X, Rojas-Silva X, Nario A, et al. Absorption and mobility of radio-labelled calcium in chili pepper plants and sweet cherry trees. Sci agric (Piracicaba, Braz). 2021;78:e20200092.\u003c/li\u003e\n\u003cli\u003eWu S, Tian L. A new flavone glucoside together with known ellagitannins and flavones with anti-diabetic and anti-obesity activities from the flowers of pomegranate ( \u003cem\u003ePunica granatum\u003c/em\u003e ). Natural Product Research. 2019;33:252\u0026ndash;7.\u003c/li\u003e\n\u003cli\u003eWilson AE, Wu S, Tian L. PgUGT95B2 preferentially metabolizes flavones/flavonols and has evolved independently from flavone/flavonol UGTs identified in \u003cem\u003eArabidopsis thaliana\u003c/em\u003e. Phytochemistry. 2019;157:184\u0026ndash;93.\u003c/li\u003e\n\u003cli\u003eWu S, Tian L. Diverse phytochemicals and bioactivities in the ancient fruit and modern functional food pomegranate (\u003cem\u003ePunica granatum\u003c/em\u003e). Molecules. 2017;22:1606.\u003c/li\u003e\n\u003cli\u003eGouveia CSS, Ganan\u0026ccedil;a JFT, Slaski JJ, Lebot V, Pinheiro de Carvalho M\u0026Acirc;A. Involvement of abscisic acid and other stress indicators in taro (Colocasia esculenta (L.) Schott) response to drought conditions. Acta Physiol Plant. 2020;42:173.\u003c/li\u003e\n\u003cli\u003eLiang L, Xiao W, Lei Y, Ling K. Quantitative evaluation of 4-coumarate: coa ligase (4cl) activity and correlated chemical constituents in four plant materials by chromatographic analysis. Chinese Journal of Natural Medicines. 2010; 8(4): 274-279\u003c/li\u003e\n\u003cli\u003eSiwinska J, Kadzinski L, Banasiuk R, Gwizdek-Wisniewska A, Olry A, Banecki B, et al. Identification of QTLs affecting scopolin and scopoletin biosynthesis in Arabidopsis thaliana. BMC Plant Biol. 2014;14:280.\u003c/li\u003e\n\u003cli\u003eKaramat F, Olry A, Doerper S, Vialart G, Ullmann P, Werck-Reichhart D, et al. CYP98A22, a phenolic ester 3\u0026rsquo;-hydroxylase specialized in the synthesis of chlorogenic acid, as a new tool for enhancing the furanocoumarin concentration in Ruta graveolens. BMC Plant Biol. 2012;12(1):1-13.\u003c/li\u003e\n\u003cli\u003eNakatsubo T, Kitamura Y, Sakakibara N, Mizutani M, Hattori T, Sakurai N, et al. At5g54160 gene encodes \u003cem\u003eArabidopsis thaliana\u003c/em\u003e 5-hydroxyconiferaldehyde O-methyltransferase. J Wood Sci. 2008;54:312\u0026ndash;317.\u003c/li\u003e\n\u003cli\u003eChun H J, Baek D, Cho H M, et al. Lignin biosynthesis genes play critical roles in the adaptation of Arabidopsis plants to high-salt stress. Plant signaling \u0026amp; behavior. 2019; 14(8): 1625697.\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":true,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"Camellia limonia, karst, high calcium stress, metagenomics, metabonomics","lastPublishedDoi":"10.21203/rs.3.rs-135639/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-135639/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003e\u003cstrong\u003eBackground:\u003c/strong\u003e Yellow \u003cem\u003eCamellia \u003c/em\u003eis a kind of rare plant with high economic and medicinal value. It is known as the \"giant panda\" of the botanical world. \u003cem\u003eCamellia \u003c/em\u003egroup of plants (Camellia Sect. Chrysantha Chang) is the only golden flower of \u003cem\u003eTheaceae\u003c/em\u003e. Compared with \u003cem\u003eCamellia nitidissima, Camellia limonia \u003c/em\u003egrows in karst areas, where the soil has the characteristic of high calcium content. However, there are few studies about the \u003cem\u003eCamellia limonia \u003c/em\u003ein karst soil environment and the adaptation mechanism is no clear. \u003c/p\u003e\u003cp\u003e\u003cstrong\u003eResults:\u003c/strong\u003e In this study, we found that under high calcium treatment, the chlorophyll content and leaf areas of \u003cem\u003eCamellia limonia \u003c/em\u003eincreased, while those of \u003cem\u003eCamellia nitidissima \u003c/em\u003edecreased. The photosynthetic efficiency of \u003cem\u003eCamellia limonia \u003c/em\u003ewas more stable and higher than \u003cem\u003eCamellia nitidissima\u003c/em\u003e. Compared with \u003cem\u003eCamellia nitidissima\u003c/em\u003e, the conductance was larger and the degree of leaf shrinkage was smaller\u003cem\u003e \u003c/em\u003ein \u003cem\u003eCamellia limonia\u003c/em\u003e. The metabolomics analysis showed that the kaempferol-3-o-rutinoside, tyrosol, 6-o-methyldeacetylisoipecoside and (r)-mandelic acid are the main differently metabolic compounds . The results of karst high calcium soil metagenomics showed that microbacterium-testaceum, intrasporangium-calvum and rubrivivax-gelatinosus significantly changed. Through metabonomics and metagenomics integrative analysis, flavone and flavonol biosynthesis is suggested to be the main regulation pathway, which is regulated by apigenin, kaempferol, astragalin, isoquercitrin metabolites and \u003cem\u003eTT7, UGT78D1\u003c/em\u003e, \u003cem\u003eUGT78D2\u003c/em\u003e genes. This metabolic pathway involves the synthesis of flavonoids. Flavonoids have the functions in drought and salt resistance, which play an important role in the adaptation of \u003cem\u003eCamellia limonia\u003c/em\u003e in karst high calcium environment.\u003c/p\u003e\u003cp\u003e\u003cstrong\u003eConclusion\u003c/strong\u003e: This omics study identified key regulation metabolites and genes for \u003cem\u003eCamellia\u003c/em\u003e and provided important basis for the adaptive mechanism of plants to adapt to the high-calcium environment and the protection of \u003cem\u003eCamellia \u003c/em\u003especies.\u003c/p\u003e","manuscriptTitle":"Exploring the adaptive mechanism of Camellia limonia in karst high calcium environment via an integrative analysis of metabolomics and metagenomics","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2021-01-08 22:05:03","doi":"10.21203/rs.3.rs-135639/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"67eb9c43-9cd0-4011-aeab-9aa87dc3ff83","owner":[],"postedDate":"January 8th, 2021","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[{"id":1683198,"name":"Plant Physiology and Morphology"},{"id":1683199,"name":"Plant Molecular Biology and Genetics"}],"tags":[],"updatedAt":"2021-01-29T05:29:17+00:00","versionOfRecord":[],"versionCreatedAt":"2021-01-08 22:05:03","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-135639","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-135639","identity":"rs-135639","version":["v1"]},"buildId":"_2-kVJe1T_tPrBINL-cwx","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}