• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Effects of long-term fencing on soil microbial community structure and function in the desert steppe, China

    2024-03-19 02:29:32PANYaqingKANGPengQUXuanRANYichaoLIXinrong
    Journal of Arid Land 2024年3期

    PAN Yaqing, KANG Peng, QU Xuan, RAN Yichao, LI Xinrong*

    1 Shapotou Desert Research and Experiment Station, Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Lanzhou 730000, China;

    2 School of Biological Science and Engineering, North Minzu University, Yinchuan 750021, China

    Abstract: One of the goals of grazing management in the desert steppe is to improve its ecosystem.However, relatively little is known about soil microbe communities in the desert steppe ecosystem under grazing management.In this study, we investigated the diversity and aboveground biomass of Caragana korshinskii Kom.shrub communities in long-term fencing and grazing areas, combined with an analysis of soil physical-chemical properties and genomics, with the aim of understanding how fence management affects plant-soil-microbial inter-relationships in the desert steppe, China.The results showed that fence management (exclosure) increased plant diversity and aboveground biomass in C. korshinskii shrub area and effectively enhanced soil organic carbon (233.94%), available nitrogen (87.77%), and available phosphorus (53.67%) contents.As well, the Shannon indices of soil bacteria and fungi were greater in the fenced plot.Plant-soil changes profoundly affected the alpha- and beta-diversity of soil bacteria.Fence management also altered the soil microbial community structure, significantly increasing the relative abundances of Acidobacteriota (5.31%-8.99%), Chloroflexi (3.99%-5.58%), and Glomeromycota(1.37%-3.28%).The soil bacterial-fungal co-occurrence networks under fence management had higher complexity and connectivity.Based on functional predictions, fence management significantly increased the relative abundance of bacteria with nitrification and nitrate reduction functions and decreased the relative abundance of bacteria with nitrate and nitrite respiration functions.The relative abundances of ecologically functional fungi with arbuscular mycorrhizal fungi, ectomycorrhizal fungi, and saprotrophs also significantly increased under fence management.In addition, the differential functional groups of bacteria and fungi were closely related to plant-soil changes.The results of this study have significant positive implications for the ecological restoration and reconstruction of dry desert steppe and similar areas.

    Keywords: desert steppe; fence management; Caragana korshinskii; soil physical-chemical property; soil microorganism

    1 Introduction

    Desert steppe (the transition zone from a typical steppe to a desert) is an important ecosystem,unique to Eurasia, and has a long history of livestock grazing (Kang et al., 2007; Zhang et al.,2018).In the desert steppe ecosystem with relatively simple plant species richness and community structure, high-intensity grazing exacerbates the environmental pressure on plant community (Herrero-Jáuregui and Oesterheld, 2018; Gao and Carmel, 2020).In order to alleviate this threat and mitigate the negative impacts on grassland ecosystem, government in China has initiated a series of ecological protection policies and restoration projects since the late 1990s and early 2000s.For example, in 2003, the Returning Pasture to Grassland project was implemented,which uses a system of grazing bans and rotational grazing to reduce grazing pressure on natural grasslands (Xiong et al., 2016; Xu et al., 2016).Subsequently, the Grassland Ecological Conservation project was implemented in 2011 to mitigate grassland degradation and increase grassland vegetation biomass (Hao et al., 2014).The implementation of these ecological engineering projects has had a positive impact on the structure of plant communities and the physical and chemical properties of soil, contributing to the restoration of the grassland ecosystem(Li et al., 2017; Sun et al., 2023).

    Researchers generally agree that improving grazing management system is the principal strategy for the revegetation and reconstruction of arid and semi-arid grassland ecosystem (Xu et al., 2018).Limiting selective foraging by herbivores directly alters the availability of light to subcanopy species, which, in turn, affects the structure and composition of plant community (Sun et al., 2023).However, studies have shown that the response of species richness to fencing is not consistent.For example, long-term fencing of grasslands in the Mu Us Sandy Land, China increased vegetation cover and the Shannon-Wiener diversity index, but decreased plant species richness (Liu et al., 2019).In shrubbery grassland ecosystem of Inner Mongolia Autonomous Region, China, long-term fencing was found to increase Shannon-Wiener diversity and species richness (Hao et al., 2019).The responses of plant community diversity and richness to fencing in desert grassland ecosystem dominated by xerophytic shrubs deserve further investigation.

    In the desert steppe ecosystems, the physical and chemical properties of soil have long been important factors influencing vegetation cover.The trampling of vegetation by livestock has been shown to destroy the soil surface crust and aggregates and increase soil erodibility (Munkhtsetseg et al., 2017).The loss of soil nutrients is further exacerbated by external forces, such as wind(Aubault et al., 2017).In addition, the degradation of vegetation induced by grazing reduces the quantity and quality of aboveground plant litter input, which adversely affects the physical and chemical properties of soil (Penner and Frank, 2019).The implementation of fencing usually increases vegetation cover in arid and semi-arid grassland ecosystem, reduces evapotranspiration,and facilitates water storage in soil (Mekuria and Veldkamp, 2012).The removal of animal feces and urine can also effectively reduce soil pH (Fan et al., 2020).Furthermore, fencing protects the aboveground and belowground parts of plants and further enhances plant-soil interactions, which not only facilitates the accumulation of soil organic carbon but also effectively improves the accumulation of soil nutrients, such as available nitrogen and phosphorus (Zheng et al., 2023).

    It is well known that soil microorganisms are an important part of ecosystem and an important medium for the transfer of nutrients and information between plants and environment (Faghihinia et al., 2020).Studies have shown that changes in plant community structure and soil physical-chemical properties as a result of fencing also affect the structure and function of soil microbial community (Aubault et al., 2015; Zhang et al., 2016).For example, fencing significantly increases the abundance of soil bacteria and fungi in desert grasslands (Fan et al.,2020), and alters the relative abundances of Actinobacteria, Proteobacteria, and Ascomycota (Zhu et al., 2021).Additionally, the effects of fencing on soil moisture, pH, and nutrients create ecological niches for underground microbial community (Yao et al., 2018).Microbial network relationships can further explain this change.More stable and connected network relationships are conducive to the occupation of ecological niches by bacteria and fungi and play a crucial role in maintaining the stability of microbial communities (da Costa et al., 2021; Fu et al., 2022).For example, fencing increases the efficiency of soil microbial nitrogen utilization, reduces soil N2O loss (Zhang et al., 2021), and increases the relative abundance of functional fungal groups(Eldridge and Delgado-Baquerizo, 2018).However, the responses of soil bacterial and fungal functional groups to fencing in desert grasslands have not been sufficiently investigated.

    In many studies of grassland management, the focus has been on the synergistic changes and responses of plant-soil microorganisms under different grazing regimes.In view of the importance of soil microbial communities in ecosystem restoration, we focused our attention on theCaragana korshinkiiKom.shrub area of the Ordos Plateau, northern China, an area between the Helan Mountains to the west and Mu Us Sandy Land to the east.C.korshinkiiis the dominant species for the ecological restoration of desert grasslands and has the characteristics of wind prevention and drought resistance (Chen et al., 2022).Through the investigation of plant communities, analysis of soil physical and chemical properties, and results of bacterial 16S ribosomal ribonucleic acid (rRNA) gene and internal transcribed spacer (ITS) high-throughput sequencing, the following questions were addressed: (1) whether changes in plant community structure and soil physical-chemical properties under fence management positively affect ecological niche competition among soil microorganisms; and (2) how potential functions of soil bacteria and fungi respond to fence management.

    2 Materials and methods

    2.1 Study area and experimental design

    The study area was located in the Otog Front Banner, Inner Mongolia Autonomous Region, China(37°43′-37°55′N, 107°21′-107°38′E; 1300-1350 m a.s.l.; Fig.1) with a typical temperate continental climate.The sampling time was in late June 2021, the planting time ofC.korshinskiiwas approximately 2004, and the fences were installed in 2014.Nine study plots (5 m×5 m) were placed inside and outside the fenced enclosures, and the samples were spaced more than 50 m apart.The diversity and richness of plants were assessed using the per-tree survey method.The plants studied includedC.korshinskii,Achnatherum splendens(Trin.) Nevski,Stipa glareosaP.A.Smirn.,Lespedeza potaniniiVass.,Cynanchum komaroviiAl.,andSophora alopecuroidesL.

    Fig.1 Study area (a) and study plots (b and c) in the Otog Front Banner, Inner Mongolia Autonomous Region,China

    2.2 Soil sample collection

    Soil samples were collected in August 2022.Within the study area, nine grazing plots (5 m×5 m)and nine fencing plots (5 m×5 m) were divided, of which we further divided five sub-samples in each quadrat based on five-point sampling method, and mixed as one soil sample (Fig.S1).A total of 18 soil samples were obtained.The sampling depth was 0-20 cm, and the samples were taken after removing litter from the surface and avoiding plant roots (Pan et al., 2021).These collected soil samples were immediately placed in an icebox and returned to the laboratory (Kang et al., 2022).These fresh soil samples were sent to Novogene Bioinformatics Technology Co.,Ltd., Beijing, China and stored at low temperature for the determination of soil bacteria and fungi.

    2.3 Analysis of soil physical-chemical properties

    Air-dried soil was used to determine soil physical-chemical properties.The soil samples were sieved and weighed to determine soil water content (SWC).Ten grams of soil was weighed in a beaker, 50 mL of deionized water was added, a conductivity meter (DDS-307A, LEICI Scientific Instruments Co., Ltd., Shanghai, China) was used to determine the soil electrical conductivity(EC), and a pH meter (PHS-3C, LEICI Scientific Instruments Co., Ltd., Shanghai, China) was used to determine soil pH (Bao, 2000).Soil samples were air-dried and an Elementar Vario Macro cube was used to determine total carbon (TC) and nitrogen (TN) contents in the soil.Total phosphorus (TP) was measured using the alkaline potassium persulfate digestion method(Bettinelli et al., 2000).Soil organic carbon (SOC) content was determined using the dichromate oxidation method, available nitrogen (AN) content was determined using the alkali hydrolysis method, and available phosphorus (AP) content was determined using the ammonium molybdate colorimetric method (Bao, 2000; Bettinelli et al., 2000).

    2.4 Polymerase chain reaction (PCR) amplification and sequence processing

    Soil deoxyribonucleic acid (DNA) was extracted from 18 soil samples using the hexadecyl trimethyl ammonium bromide (CTAB) method.The soil bacterial V3-V4 gene region and fungal internal transcribed spacer-1 (ITS-1) were amplified using primers (341F, 806R, ITS1F, and ITS2)(Claesson et al., 2009; Cheng et al., 2022).After PCR amplification, the products were extracted from a 2% agarose gel and purified for quantification (Pan et al., 2021; Kang et al., 2022).We used the NEXTflex? Rapid DNA-Seq Kit (Bioo Scientific, Austin, USA) to construct the PE(paired-end) library and sequenced it using a NovaSeq PE250 platform.Trimmomatic method was used for quality control of the raw MiSeq data, which was then merged with fast length adjustment of short reads (FLASH) v.1.2.11 (Wang et al., 2007).UPARSE v.7.0.1090 was used to perform operational taxonomic unit (OTU) cluster analysis on the sequence (similarity=97%)(Edgar, 2013).The representative sequence of OTU was obtained by removing chimeras during the clustering process.The species classification annotation of each sequence was performed using ribosomal database project (RDP) classifier and compared with the Silva database (SSU128)using a confidence threshold of 0.7 (Quast et al., 2013).The 16S rRNA and ITS gene sequences obtained in this study were submitted to NCBI (National Center for Biotechnology Information)SRA (Sequence Read Archive) database with serial numbers PRJNA1026058 and PRJNA1026303, respectively.

    2.5 Data analysis

    Plant community indices, including Shannon-Wiener (H), Simpson (H'), Pielou (E), and Patrick(R) indices (Shannon, 1948; Zhang et al., 2003) were calculated using the following equations:

    wherePiis the relative importance value of theithspecies; andSis the number of species.

    The above-ground biomass (AGB; kg) estimation in this study was performed using the method described by Conti et al.(2019).

    whereBDis the ground diameter (cm);CDis the crown breadth (cm); andHPis the plant height(cm).

    Mantel's test was used to describe the correlation between plant community indices and soil physical-chemical properties (Sunagawa et al., 2015).Quantitative insights into microbial ecology (QIIME) v.1.9.1 software was used to analyze the microbiome alpha-diversity (OTUs,Shannon, and richness) (Caporaso et al., 2012).Non-metric multidimensional scaling analysis(NMDS) was used to analyze the beta-diversity of bacterial and fungal communities (Bell et al.,2014).The phyla (relative abundance>1%) of the microbial community were analyzed using the"ggalluvial" package in the R software program (Pan et al., 2021).The relative abundances of bacterial and fungal phyla were analyzed using Spearman's test (P<0.05) (Kang et al., 2022).We calculated Spearman's coefficients for the OTUs of bacteria and fungi under different plots(|r|>0.9,P<0.001) and then displayed the co-occurrence network with established data (e.g.,showing the number of nodes in Actinobacteria, Proteobacteria, Acidobacteria, Chloroflexi, and Ascomycota) using Cytoscape software v.3.7.1, and data on nodes and edges were obtained (Kang et al., 2022).Additionally, we calculated the network parameters for each OTU for subsequent analyses.We used redundancy analysis (RDA) to assess the effects of fence management on soil physical-chemical properties and soil microbial communities.Subsequently, the biogeochemical functions of bacterial communities were predicted using the Functional Annotation of Prokaryotic Taxa (FAPROTAX) (Louca et al., 2016).Fungi+Functional+Guild (FUNGuild) has also been used to predict the ecological functions of fungal communities (Nguyen et al., 2016).The bacterial and fungal functional groups between the two plots were analyzed using one-way analysis of variance (ANOVA) and Duncan's multiple range tests (P<0.050).Finally, Spearman's analysis showed the correlation of plant community indices, soil physical-chemical properties,and differential functional taxa (Pan et al., 2022).

    3 Results

    3.1 Effects of fencing on plant community structure and soil physical-chemical properties

    Fig.2 Plant community indices (a) of Caragana korshinskii shrub and their correlation with soil physical-chemical properties of the GRA (b) and FEN (c) plots in the desert steppe.GRA, grazing; FEN, fencing.In Figure 1a, different lowercase letters indicate significant differences between the two plots at P<0.050 level.Lines in the box are median values.Bars are standard errors.*, P<0.050 level; **, P<0.010 level; ***, P<0.001 level.AGB, above-ground biomass; SWC, soil water content; EC, electrical conductivity; TC, total carbon; TN,total nitrogen; TP, total phosphorus; SOC, soil organic carbon; AN, available nitrogen; AP, available phosphorus.The abbreviations are the same in the following figures.

    Fencing increased Shannon-Wiener, Simpson, and Pielou indices ofC.korshinskiicommunities in the desert steppe and significantly increased plant AGB, but had no significant effect on the Patrick index (Fig.2a).Fencing also decreased soil pH and EC, but had little effect on SWC.Further analysis revealed that fencing significantly increased soil TC, TN, and TP contents, with more pronounced increases in SOC, AN, and AP contents (Table 1).In addition, fencing significantly increased SOC/TN and SOC/TP ratios.Mantel's test showed that Shannon-Wiener,Simpson, and Pielou indices had significant correlation with SOC/TP ratio in the grazing plot,respectively, whereas AGB was significantly related to EC (P<0.050; Fig.2b).In the fencing plot,Shannon-Wiener, Simpson, and Pielou indices were significantly correlated with SOC/TN ratio as well as SOC/TP ratio (P<0.050; Fig.2c).

    Table 1 Soil physical-chemical properties of the GRA and FEN plots in the desert steppe

    3.2 Effects of fencing on soil microbial diversity

    The number of OTUs and Shannon indices of bacteria and fungi in the soil increased under fencing, but the richness index was not significantly different (P<0.050; Fig.3a and b).NMDS analyses showed that the stress values of bacterial and fungal communities were less than 0.2000,indicating that there were between-group differences in both bacterial and fungal communities in both the grazing and fencing plots, and that the differences in soil fungal communities were more pronounced in the two plots (Fig.3c and d).

    3.3 Effects of fencing on soil microbial community structure

    In the fencing and grazing plots, 84.36% of the OTUs were annotated at the phylum level, with Actinobacteriota, Proteobacteria, Acidobacteriota, and Firmicutes as the dominant phyla in all samples.Among the fungi, 70.33% of the OTUs were annotated at the phylum level, with Ascomycota, Basidiomycota, Mortierellomycota, and Glomeromycota being the dominant phyla in all samples (Fig.4a and b).Fencing significantly increased the relative abundances of Acidobacteriota (5.31%-8.99%) and Chloroflexi (3.99%-5.58%), whereas it decreased the relative abundance of Proteobacteria (25.34%-15.24%).However, fencing also significantly increased the relative abundance of Glomeromycota (1.37%-3.28%) (Fig.4c and d).

    3.4 Effects of fencing on soil bacterial-fungal network relationships

    In both two plots, the OTUs of Actinobacteriota, Proteobacteria, Acidobacteriota, Chloroflexi, and Ascomycota exhibited large ecological niches in the network (>5% of the nodes).In the grazing plot, the number of nodes of Actinobacteria and Proteobacteria in the network was higher than in the fencing plot (Fig.5a and b).The grazing plot had 1311 nodes and 2926 edges, with 2216 positive and 710 negative edges.The fencing plot contained 1715 nodes and 3706 edges, with 1965 positive and 1741 negative edges.In addition, the average path length and graph diameter parameters of the bacterial-fungal co-occurrence network were also higher in the fencing plot than in the grazing plot (Table 2).

    Fig.3 Alpha- and beta-diversity of soil bacterial (a1-a3 and c) and fungal (b1-b3 and d) communities of the GRA and FEN plots in the desert steppe.In Figure 3a and b, different lowercase letters indicate significant differences between the two plots at P<0.050 level.Lines in the box are median values.Bars are standard errors.OTUs, operational taxonomic units; ACE, abundance-based coverage estimator; NMDS, non-metric multidimensional scaling.

    Fig.4 Relative abundance of bacteria (a) and fungi (b) greater than 1% at the phylum level, and the difference of relative abundance of bacteria (c) and fungi (d) of the GRA and FEN plots in the desert steppe.Bars are standard errors.

    Fig.5 Co-occurrence network of soil bacteria and fungi of the GRA (a) and FEN (b) plots in the desert steppe

    Table 2 Network topological features of the GRA and FEN plots in the desert steppe

    3.5 Effects of fencing on soil bacterial and fungal functions

    RDA of plant community indices and soil physical-chemical properties revealed that RDA1 explained 82.96% of the soil physical-chemical properties, and AGB had a greater effect on soil physical-chemical properties (P=0.005; Fig.6).Further RDA of soil microbial community by soil physical-chemical properties revealed that RDA1 explained 45.34% of the soil microbial community, and all the factors except SWC and TN/TP ratio, had a strong correlation on soil microbial community (P<0.050; Fig.6b).

    Functional predictions of bacterial (FAPROTAX) and fungal (FUNGuild) communities showed significant differences in nitrogen cycle functional taxa (e.g., aerobic ammonia oxidation,nitrification, denitrification, nitrite respiration, and fermentation) between grazing and fencing bacterial communities; whereas saprotroph fungi and ectomycorrhizal were the dominant functional taxa of fungal communities (Table S1).Through a FAPROTAX prediction of soil bacterial function, we found that fencing significantly increased the relative abundance of bacteria with nitrification and nitrate reduction functions, while decreasing the relative abundance of bacteria with nitrate and nitrite respiration functions.FUNGuild predictions showed that the relative abundance of fungi with ecological functions, such as arbuscular mycorrhizal and ectomycorrhizal fungi, increased under fencing.In addition, we analyzed the differential functions of bacteria and fungi with plant community indices and soil physical-chemical properties in the two plots, and in the grazing plot, SOC, SOC/TN and SOC/TP ratios were negatively correlated with pathogen-saprotroph taxa (Fig.7a).In the fencing plot, TN/TP ratio showed negative correlations with aerobic chemoheterotrophy and chemoheterotrophy (P<0.050;Fig.7b).

    Fig.6 RDA (redundancy analysis) of plant community indices (a) and soil physical-chemical properties (b) of the GRA and FEN plots in the desert steppe

    Fig.7 Spearman's correlation of differential functioning of soil bacteria and fungi with soil physical-chemical properties of the GRA (a) and FEN (b) plots in the desert steppe.*, P<0.050 level, **, P<0.010 level.

    4 Discussion

    4.1 Soil nutrient accumulation in the desert steppe under fencing

    Fencing is a major driver of plant community structure and function in the desert steppe (Ye et al.,2023).The results of our study showed that fencing significantly increased Shannon-Wiener,Simpson, and Pielou indices, and AGB compared with those under grazing (Fig.2).It has been shown that fencing increases plant Shannon-Wiener and Simpson indices and effectively improves the biodiversity of desert areas (Tang et al., 2016).Zhang et al.(2018) pointed out that grazing intensity was negatively correlated with AGB, and helped to maintain higher species diversity and productivity in the desert steppe of northern China.The results of the present study confirmed the previous findings and also showed that pH and EC were significantly lower in the fencing plot compared with those of the grazing plot (Table 1).Li et al.(2014) noted that reduced AGB due to grazing increased surface water evaporation (resulting in surface salt accumulation),which in turn increased soil pH and EC.In addition, the accumulation of herbivorous urine and feces has been reported to increase soil pH (Li et al., 2014).Notably, Shannon-Wiener, Simpson,and Pielou indices were significantly correlated with SOC/TP ratio in both grazing and fencing plots (Fig.2b and c).Under the background of limited soil phosphorus, the input of plant litter causes the increase of SOC content (Deng et al., 2014; Bai et al., 2020).In this study, the effect of fencing on the biomass of plant species in the desert steppes is one of the most important factors affecting the physical properties of soil.

    Findings from several studies demonstrate that an increase in plant AGB after fencing, as well as in litter and root exudates, contributes to soil carbon and nitrogen accumulation (Wu et al.,2014).Similarly, in the present study, soil in the fencing plot had higher TC and TN contents compared with those in the grazing plot (Table 1).Wang et al.(2016) found that soil carbon and nitrogen stocks changed most significantly in the surface layer (0-30 cm) after fencing and pointed out that fencing is an effective measure for improving carbon and nitrogen pools of the desert steppe ecosystems.It is well known that grassland productivity in northern China is constrained by soil phosphorous content (Lu et al., 2023).In our study, fencing was found to increase soil TP content (Table 1).A meta-analysis also indicated that fencing could increase TP content of grassland in arid and semi-arid areas of northern China (Wang et al., 2018).Fencing can improve carbon, nitrogen, and phosphorous contents in the surface layer of desert steppe soils.

    SOC is a major source of soil nutrients, and fencing can directly affect SOC content by influencing the accumulation of AGB (Wilson et al., 2018).These findings are further supported by the fact that SOC content was more significantly associated with AGB in our study and that fencing significantly increased SOC content (Table 1).In addition, fencing also increased the contents of AN and AP, and it has been noted that soil TN and AN contents were not only significantly correlated after long-term fence management but also effectively increased soil TN and AN contents (Zhu et al., 2016).There are conflicting reports in the literature; in some cases,short-term fencing had less of an effect on soil AN.However, in the present study, we report results from a long-term (20 a) exclosure experiment, which may explain the more pronounced effect on AN.Another explanation is thatC.korshinskiibelongs to the leguminous family and thus has a stronger biological N fixation capacity than common ordinary forage.As a result, plant growth was further promoted in the absence of grazing, which improved the amount of AN in the soil (Doll, 2023).Overall, fencing is an effective measure for increasing grassland productivity and improving soil fertility, particularly in the arid and semi-arid desert steppes of northern China.

    4.2 Soil microbial community structure in the desert steppe under fencing

    In our study, fencing significantly increased the diversity of soil bacterial and fungal communities(Fig.3).It has been demonstrated that ecological restoration of grasslands by fencing increases the diversity of soil microbial bacteria and fungi along with plant diversity and biomass (Hu et al.,2016).In addition, there were differences between soil bacterial and fungal communities in both the two plots, which have been observed in other study (Wang et al., 2022).Further analysis found that the dominant groups of soil microorganisms in the desert steppe are Proteobacteria,Acidobacteria, Bacteroidetes, Firmicutes, Ascomycota and Basidiomycota (Fan et al., 2021; Kang et al., 2023).Our study not only revealed changes in the relative abundance of these dominant phyla but also found that Acidobacteriota has a high explanatory rate in the desert steppe.Zhang et al.(2021) and Li et al.(2022) pointed out that fencing could increase the relative abundance of Acidobacteria in the soil and was positively correlated with SOC, TN, AK, and AP.

    The response characteristics of soil microbial communities to environmental changes have been visualized using co-occurrence networks (de Vries and Wallenstein, 2017).We constructed soil bacterial-fungal co-occurrence networks for the grazing and fencing plots and found that fencing had more complex network relationships.The number of nodes and edges in the fencing plot was higher than in the grazing plot, and the average path length and graph diameter were also higher than those in the grazing plot.Previous studies have found that fencing increases the complexity of soil bacterial and fungal networks (Jing et al., 2023), which plays a crucial role in maintaining the stability of microbial communities (Fu et al., 2022).We further analyzed and found that the nodes and edges of the bacterial-fungal co-occurrence network were higher than under fencing than under grazing, with higher network complexity and connectivity.More studies have indicated that increased grazing intensity reduces the complexity of soil bacterial-fungal co-occurrence networks; more stable and connected network relationships favor improved plant production performance (da Costa et al., 2022).Therefore, soil bacterial-fungal network relationships can also be an important basis for evaluating fence management in the desert steppe.

    4.3 Function of soil microbial community in the desert steppe under fencing

    Soil microorganisms play key roles in energy flow and material circulation in ecosystems(Veresoglou et al., 2015).In our study, FAPROTAX predictions showed that fencing significantly increased the relative abundance of bacteria with nitrification and nitrate reduction functions,while decreasing the relative abundance of bacteria with nitrate and nitrite respiration functions.Similarly, a recent study by Zhang et al.(2023) found that grazing increased the relative abundance of soil bacteria with denitrification functions.It is worth noting that there was a correlation between SOC/TN ratio, SOC/TP ratio, and nitrogen cycle taxa in the fencing plot,which may be caused by the effect of increased SOC content on denitrification processes (Zhang et al., 2021).Additional studies have also pointed out the effects of grazing on the relative abundance of bacteria with nitrogen cycling functions, which may lead to a loss in soil nutrient availability (Olsen et al., 2011; Pan et al., 2018).Combined with plant communities and soil physical-chemical properties, it is clear that fencing has a positive effect on soil nitrogen cycling and accumulation.

    In this study, FUNGuild predicted that the relative abundance of fungi with ecological functions, such as arbuscular mycorrhizal and ectomycorrhizal fungi, could be enhanced by fencing (Table S1).Arbuscular mycorrhizal fungi are ubiquitous and specialized symbiotic fungi belonging to Glomeromycota (?milauer et al., 2021).According to Figure 4, we found that the relative abundance of Glomeromycota was significantly higher in the fencing plot than in the grazing plot, which may be an important reason for the increase of arbuscular mycorrhizal and ectomycorrhizal fungi groups.In view of the important role of mycorrhizal fungi in above- and below-ground multi-nutrient cycling, it can be further inferred that the relative abundance of mycorrhizal fungi groups increases by fencing, thus indirectly enhancing soil nutrient cycling (Fei et al., 2022; Chen et al., 2023).In addition, the relative abundance of saprotrophic fungi increased significantly by fencing in this study.Previous studies have reported that an increase in SOC increases competition by saprophytic fungi groups (Kang et al., 2023).Our study further demonstrated the effect of fencing on the potential relationship between saprotrophic fungi groups and soil physical-chemical properties.

    5 Conclusions

    In the context of increased disturbances from grazing, excluding livestock through fencing can restore and protect plant community structure and soil physical-chemical properties in the desert steep ecosystems.As shown by the results, soil microbial community structure and function respond positively to fencing.The study demonstrated that fencing improved soil by increasing AGB levels in the soil ofC.korshinskii.In addition, long-term fencing resulted in marked enhancement of SOC, AN, and AP contents, as well as increased Shannon indices of bacteria and fungi in the soil.Plant-soil changes profoundly affect the alpha- and beta-diversity of soil bacteria.The exclusion of grazing also altered the soil microbial community structure, significantly increasing the relative abundances of Acidobacteriota, Chloroflexi, and Glomeromycota.The soil bacterial-fungal co-occurrence networks under fencing had higher complexity and connectivity.RDA demonstrated the role of AGB in influencing the desert steppe soil environment and the importance of SOC, AN, and AP contents for soil microbial communities.Based on functional predictions, fencing significantly increased the relative abundance of bacteria with nitrification and nitrate reduction functions and decreased the relative abundance of bacteria with nitrate and nitrite respiration functions.The relative abundances of ecologically functional fungi with arbuscular mycorrhizal fungi, ectomycorrhizal fungi, and saprotrophs also significantly increased under fencing.In addition, the differential functional groups of bacteria and fungi are closely related to plant-soil changes.The comprehensive results of this long-term study offer a solid basis for the ecological restoration and reconstruction of the dry desert steppes.

    Conflict of interest

    LI Xinrong is an editorial board member of Journal of Arid Land and was not involved in the editorial review or the decision to publish this article.All authors declare that there are no competing interests.

    Acknowledgements

    This research was funded by the National Natural Science Foundation of China (32061123006, 32360426).

    Author contributions

    Conceptualization: PAN Yaqing, KANG Peng, LI Xinrong; Methodology: PAN Yaqing, KANG Peng, QU Xuan;Investigation: PAN Yaqing, KANG Peng, QU Xuan, RAN Yichao; Formal analysis: KANG Peng, RAN Yichao;Writing - original draft preparation: PAN Yaqing, KANG Peng; Writing - review and editing: PAN Yaqing, LI Xinrong; Funding acquisition: KANG Peng, LI Xinrong.All authors approved the manuscript.

    Appendix

    Fig.S1 Sampling diagram of Caragana korshinskii Kom.shrub in the desert steppe

    Table S1 Differential functional taxa of bacterial and fungal communities of the GRA and FEN plots in the desert steppe

    青春草亚洲视频在线观看| 午夜日韩欧美国产| 亚洲五月婷婷丁香| 亚洲,欧美,日韩| 精品国产一区二区三区四区第35| 一级黄片播放器| 亚洲,一卡二卡三卡| 久久久久视频综合| 麻豆av在线久日| www.av在线官网国产| 人妻一区二区av| avwww免费| 岛国毛片在线播放| 国产亚洲一区二区精品| netflix在线观看网站| 国产精品一区二区免费欧美 | 国产一卡二卡三卡精品| 另类亚洲欧美激情| 91字幕亚洲| 国产成人免费无遮挡视频| 国产精品欧美亚洲77777| 国产麻豆69| 天堂中文最新版在线下载| 18禁观看日本| 国产99久久九九免费精品| 国产精品秋霞免费鲁丝片| 亚洲色图 男人天堂 中文字幕| 日韩人妻精品一区2区三区| 国产精品av久久久久免费| 久久久精品免费免费高清| 精品亚洲成a人片在线观看| 成年人午夜在线观看视频| 99久久99久久久精品蜜桃| 国产精品一区二区在线不卡| 成人国产av品久久久| 亚洲国产日韩一区二区| 少妇精品久久久久久久| 婷婷成人精品国产| 老司机在亚洲福利影院| 啦啦啦在线免费观看视频4| 色综合欧美亚洲国产小说| 国产成人一区二区三区免费视频网站 | 99精品久久久久人妻精品| 欧美精品人与动牲交sv欧美| 人人妻人人澡人人爽人人夜夜| 免费在线观看完整版高清| 国产精品国产av在线观看| 人人妻,人人澡人人爽秒播 | 久久久久久人人人人人| 亚洲av美国av| 亚洲精品成人av观看孕妇| 一本大道久久a久久精品| 久久综合国产亚洲精品| 高清av免费在线| 中文精品一卡2卡3卡4更新| 国产一区二区三区av在线| 久久天躁狠狠躁夜夜2o2o | 色婷婷av一区二区三区视频| 免费看十八禁软件| 国产精品成人在线| 好男人电影高清在线观看| 不卡av一区二区三区| 亚洲国产欧美网| 国产成人a∨麻豆精品| 丰满迷人的少妇在线观看| 国产精品 欧美亚洲| 爱豆传媒免费全集在线观看| 亚洲伊人久久精品综合| 久久天躁狠狠躁夜夜2o2o | 波野结衣二区三区在线| 男女床上黄色一级片免费看| 午夜免费观看性视频| 美女午夜性视频免费| 日韩制服丝袜自拍偷拍| 欧美日韩黄片免| 九色亚洲精品在线播放| 99国产精品一区二区三区| 男男h啪啪无遮挡| 成人三级做爰电影| 一本—道久久a久久精品蜜桃钙片| 麻豆国产av国片精品| 久久免费观看电影| 日韩,欧美,国产一区二区三区| 赤兔流量卡办理| 久久狼人影院| 久久久久久久国产电影| 久久久久久免费高清国产稀缺| 麻豆乱淫一区二区| 久久国产精品影院| 老司机深夜福利视频在线观看 | 满18在线观看网站| 90打野战视频偷拍视频| 欧美黑人欧美精品刺激| 一区福利在线观看| 亚洲精品久久成人aⅴ小说| 另类亚洲欧美激情| 成人午夜精彩视频在线观看| 久久久久久久大尺度免费视频| a 毛片基地| 美国免费a级毛片| 人妻人人澡人人爽人人| 一区二区三区乱码不卡18| 国产高清videossex| 一区在线观看完整版| 精品人妻熟女毛片av久久网站| 久久久欧美国产精品| 在线观看一区二区三区激情| 飞空精品影院首页| 亚洲精品国产av成人精品| 无限看片的www在线观看| 人妻人人澡人人爽人人| 99国产精品一区二区三区| 老司机深夜福利视频在线观看 | 久久鲁丝午夜福利片| 赤兔流量卡办理| 国产精品久久久人人做人人爽| 国产又色又爽无遮挡免| 久久亚洲精品不卡| 无遮挡黄片免费观看| 十八禁高潮呻吟视频| 欧美精品高潮呻吟av久久| 啦啦啦啦在线视频资源| 一边摸一边做爽爽视频免费| 下体分泌物呈黄色| 午夜激情av网站| 亚洲三区欧美一区| 激情视频va一区二区三区| 国产无遮挡羞羞视频在线观看| 又紧又爽又黄一区二区| 手机成人av网站| 久久国产精品影院| 又紧又爽又黄一区二区| 99热全是精品| 日韩一区二区三区影片| 久久久久国产一级毛片高清牌| 亚洲av电影在线进入| 亚洲av电影在线进入| 国产亚洲欧美在线一区二区| 欧美黑人欧美精品刺激| 青春草视频在线免费观看| av视频免费观看在线观看| 午夜影院在线不卡| 老司机午夜十八禁免费视频| 国产国语露脸激情在线看| 亚洲精品日本国产第一区| 亚洲色图综合在线观看| 亚洲成人手机| 丝袜美足系列| 中文乱码字字幕精品一区二区三区| 国产精品免费视频内射| 日韩 亚洲 欧美在线| 精品国产超薄肉色丝袜足j| 在线观看www视频免费| 伊人久久大香线蕉亚洲五| 日韩精品免费视频一区二区三区| 国产主播在线观看一区二区 | 国产无遮挡羞羞视频在线观看| 亚洲国产成人一精品久久久| 天天躁狠狠躁夜夜躁狠狠躁| avwww免费| 免费看不卡的av| 免费在线观看完整版高清| 一本一本久久a久久精品综合妖精| 国产成人精品在线电影| 伊人亚洲综合成人网| 18在线观看网站| 久9热在线精品视频| 亚洲国产av影院在线观看| 狂野欧美激情性xxxx| 丰满迷人的少妇在线观看| 少妇猛男粗大的猛烈进出视频| 免费在线观看视频国产中文字幕亚洲 | 国产亚洲av高清不卡| 又紧又爽又黄一区二区| 老司机靠b影院| tube8黄色片| 国产精品久久久人人做人人爽| 久久女婷五月综合色啪小说| cao死你这个sao货| 91国产中文字幕| 七月丁香在线播放| 精品一区在线观看国产| 午夜福利视频精品| 日韩免费高清中文字幕av| 曰老女人黄片| 亚洲精品在线美女| 日本91视频免费播放| 男人添女人高潮全过程视频| 看免费成人av毛片| 国产成人精品久久二区二区91| 中文字幕人妻熟女乱码| 一级片'在线观看视频| 中文字幕av电影在线播放| 在现免费观看毛片| 亚洲三区欧美一区| 久久精品国产a三级三级三级| 午夜影院在线不卡| 日韩欧美一区视频在线观看| 亚洲成人国产一区在线观看 | 久久久久久久国产电影| 91成人精品电影| 深夜精品福利| 国产成人精品久久久久久| 人人妻人人澡人人看| 人人澡人人妻人| 精品少妇一区二区三区视频日本电影| 久久久久久久国产电影| kizo精华| 午夜福利免费观看在线| 亚洲五月婷婷丁香| 各种免费的搞黄视频| 激情五月婷婷亚洲| 午夜日韩欧美国产| 国产视频一区二区在线看| 亚洲伊人色综图| 18禁观看日本| 国产淫语在线视频| 日韩,欧美,国产一区二区三区| 交换朋友夫妻互换小说| 亚洲五月婷婷丁香| 爱豆传媒免费全集在线观看| 中国国产av一级| netflix在线观看网站| 自线自在国产av| 亚洲久久久国产精品| 久久久久国产精品人妻一区二区| 久久人人爽人人片av| 亚洲图色成人| 国产淫语在线视频| 国产在线视频一区二区| 99国产精品一区二区蜜桃av | 成年人午夜在线观看视频| 女人被躁到高潮嗷嗷叫费观| 男女无遮挡免费网站观看| 亚洲精品美女久久av网站| 大型av网站在线播放| 交换朋友夫妻互换小说| 免费在线观看黄色视频的| 可以免费在线观看a视频的电影网站| 黄片播放在线免费| 99国产精品一区二区三区| 啦啦啦啦在线视频资源| 国产成人系列免费观看| 黄色一级大片看看| 免费观看人在逋| 欧美亚洲 丝袜 人妻 在线| 女人精品久久久久毛片| 各种免费的搞黄视频| 高清不卡的av网站| 国产亚洲av高清不卡| 国产日韩欧美在线精品| 一区二区日韩欧美中文字幕| 你懂的网址亚洲精品在线观看| 亚洲天堂av无毛| 高清av免费在线| 王馨瑶露胸无遮挡在线观看| 亚洲中文字幕日韩| 久久狼人影院| 中文字幕另类日韩欧美亚洲嫩草| 亚洲色图 男人天堂 中文字幕| 国产亚洲午夜精品一区二区久久| 亚洲中文日韩欧美视频| 秋霞在线观看毛片| 午夜免费观看性视频| 你懂的网址亚洲精品在线观看| 女警被强在线播放| 黑人欧美特级aaaaaa片| 色婷婷久久久亚洲欧美| 999久久久国产精品视频| 精品人妻1区二区| 国产三级黄色录像| 日本av免费视频播放| 激情五月婷婷亚洲| 国产精品熟女久久久久浪| 精品人妻1区二区| 九草在线视频观看| 丝袜人妻中文字幕| 午夜久久久在线观看| 91九色精品人成在线观看| 婷婷成人精品国产| 91麻豆精品激情在线观看国产 | 十八禁网站网址无遮挡| 亚洲一码二码三码区别大吗| e午夜精品久久久久久久| 国产成人免费无遮挡视频| a级毛片黄视频| 国产亚洲av片在线观看秒播厂| www.av在线官网国产| 男的添女的下面高潮视频| 亚洲国产av新网站| 国产免费视频播放在线视频| 免费在线观看影片大全网站 | 欧美激情高清一区二区三区| 亚洲色图 男人天堂 中文字幕| √禁漫天堂资源中文www| 亚洲专区国产一区二区| 国产免费又黄又爽又色| 美女高潮到喷水免费观看| 亚洲人成电影观看| 欧美中文综合在线视频| 不卡av一区二区三区| 嫩草影视91久久| 欧美激情高清一区二区三区| 中文精品一卡2卡3卡4更新| 国产精品久久久久久精品电影小说| kizo精华| 国产熟女午夜一区二区三区| 国产成人精品久久二区二区免费| 人人妻,人人澡人人爽秒播 | 青草久久国产| 一区二区三区四区激情视频| 999久久久国产精品视频| 日韩熟女老妇一区二区性免费视频| 亚洲国产av影院在线观看| 五月天丁香电影| 80岁老熟妇乱子伦牲交| 十八禁人妻一区二区| 2018国产大陆天天弄谢| 免费在线观看日本一区| 久久国产精品影院| 日韩精品免费视频一区二区三区| 成年美女黄网站色视频大全免费| 日本一区二区免费在线视频| 国精品久久久久久国模美| 久久精品久久久久久久性| 国产精品国产av在线观看| 交换朋友夫妻互换小说| 亚洲综合色网址| 91九色精品人成在线观看| www.自偷自拍.com| 国产欧美亚洲国产| 交换朋友夫妻互换小说| 后天国语完整版免费观看| 国产片内射在线| 真人做人爱边吃奶动态| 中文字幕制服av| 欧美日本中文国产一区发布| 青草久久国产| 色视频在线一区二区三区| 热re99久久精品国产66热6| 国产伦理片在线播放av一区| 国产成人精品久久久久久| 啦啦啦在线观看免费高清www| 丝袜美足系列| 亚洲专区中文字幕在线| 1024香蕉在线观看| 中国美女看黄片| 夫妻午夜视频| 曰老女人黄片| 美女中出高潮动态图| 国产成人精品久久久久久| 精品久久蜜臀av无| 久久精品亚洲熟妇少妇任你| 日本91视频免费播放| 国产伦人伦偷精品视频| 一级毛片 在线播放| 日本欧美国产在线视频| 国产精品久久久av美女十八| 一级毛片黄色毛片免费观看视频| 亚洲av男天堂| 夜夜骑夜夜射夜夜干| 日本91视频免费播放| 男女下面插进去视频免费观看| 免费av中文字幕在线| 亚洲专区国产一区二区| 午夜影院在线不卡| 久久精品国产亚洲av高清一级| 国产成人啪精品午夜网站| 女性被躁到高潮视频| 精品亚洲乱码少妇综合久久| 夫妻性生交免费视频一级片| 免费不卡黄色视频| 一级毛片我不卡| tube8黄色片| 国产高清不卡午夜福利| 日韩一卡2卡3卡4卡2021年| 国产1区2区3区精品| 欧美+亚洲+日韩+国产| 满18在线观看网站| 中国国产av一级| 黑人巨大精品欧美一区二区蜜桃| 亚洲成色77777| 99热网站在线观看| 国产一区二区 视频在线| 激情视频va一区二区三区| 国产欧美日韩一区二区三 | 欧美国产精品va在线观看不卡| 亚洲av在线观看美女高潮| 日本a在线网址| 亚洲欧洲国产日韩| 国产成人精品久久二区二区免费| 亚洲国产精品国产精品| 黑丝袜美女国产一区| 精品久久蜜臀av无| 亚洲欧美一区二区三区国产| 人妻一区二区av| 少妇人妻 视频| 一本久久精品| 欧美精品高潮呻吟av久久| 精品国产乱码久久久久久男人| 色综合欧美亚洲国产小说| 亚洲精品日本国产第一区| 91麻豆精品激情在线观看国产 | 人成视频在线观看免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 午夜日韩欧美国产| 久久久久久久大尺度免费视频| 久久女婷五月综合色啪小说| 首页视频小说图片口味搜索 | 亚洲欧美激情在线| 热99久久久久精品小说推荐| 男人操女人黄网站| 久久青草综合色| 欧美亚洲日本最大视频资源| 亚洲中文av在线| 国产亚洲午夜精品一区二区久久| 欧美+亚洲+日韩+国产| 后天国语完整版免费观看| 亚洲免费av在线视频| av天堂久久9| 精品一区二区三区av网在线观看 | 啦啦啦在线免费观看视频4| 麻豆乱淫一区二区| 国产伦理片在线播放av一区| 一级,二级,三级黄色视频| 伦理电影免费视频| 黄网站色视频无遮挡免费观看| 国产不卡av网站在线观看| 香蕉国产在线看| 日本av手机在线免费观看| 久久久精品区二区三区| 国产欧美日韩精品亚洲av| 视频区图区小说| 又黄又粗又硬又大视频| 日韩 欧美 亚洲 中文字幕| 日韩av在线免费看完整版不卡| 青青草视频在线视频观看| 国产爽快片一区二区三区| 国产av精品麻豆| 国产精品免费大片| 激情视频va一区二区三区| 成年人免费黄色播放视频| 欧美亚洲 丝袜 人妻 在线| 亚洲色图综合在线观看| 亚洲少妇的诱惑av| 纵有疾风起免费观看全集完整版| 国产人伦9x9x在线观看| 欧美成人精品欧美一级黄| 成人亚洲欧美一区二区av| 久久久国产一区二区| 国产日韩欧美视频二区| 国产在线视频一区二区| 亚洲免费av在线视频| 亚洲精品日本国产第一区| a级片在线免费高清观看视频| 免费在线观看黄色视频的| 高清欧美精品videossex| 欧美变态另类bdsm刘玥| 波多野结衣一区麻豆| 国语对白做爰xxxⅹ性视频网站| 黑人巨大精品欧美一区二区蜜桃| 国产亚洲欧美在线一区二区| 狂野欧美激情性xxxx| 麻豆乱淫一区二区| 成人亚洲精品一区在线观看| 国产一区亚洲一区在线观看| 日韩中文字幕欧美一区二区 | 丰满人妻熟妇乱又伦精品不卡| 欧美亚洲日本最大视频资源| 精品国产乱码久久久久久男人| 亚洲欧洲精品一区二区精品久久久| 欧美日韩亚洲高清精品| 国产精品av久久久久免费| 色综合欧美亚洲国产小说| 性色av一级| 香蕉丝袜av| 国产一区二区 视频在线| 99久久人妻综合| 麻豆av在线久日| 在线观看免费高清a一片| 日韩电影二区| 中文字幕人妻熟女乱码| 亚洲精品国产色婷婷电影| 操美女的视频在线观看| 欧美日韩av久久| 黄色一级大片看看| 免费久久久久久久精品成人欧美视频| 亚洲av综合色区一区| 亚洲欧美精品自产自拍| 黑人欧美特级aaaaaa片| 免费黄频网站在线观看国产| 极品人妻少妇av视频| 自线自在国产av| 99香蕉大伊视频| 丝袜脚勾引网站| 亚洲欧美一区二区三区黑人| 日日摸夜夜添夜夜爱| 精品人妻熟女毛片av久久网站| 青春草亚洲视频在线观看| 热99久久久久精品小说推荐| 国产人伦9x9x在线观看| 在线观看免费视频网站a站| 中文精品一卡2卡3卡4更新| 欧美日韩亚洲高清精品| 色精品久久人妻99蜜桃| 亚洲欧美一区二区三区国产| 1024视频免费在线观看| a 毛片基地| 99国产综合亚洲精品| 久久久久久久精品精品| 少妇被粗大的猛进出69影院| 爱豆传媒免费全集在线观看| 中文字幕制服av| 久久99精品国语久久久| 在线 av 中文字幕| 大陆偷拍与自拍| 啦啦啦中文免费视频观看日本| 国产又色又爽无遮挡免| 中文精品一卡2卡3卡4更新| 久久精品国产a三级三级三级| 如日韩欧美国产精品一区二区三区| 成年人午夜在线观看视频| 亚洲人成77777在线视频| 久久热在线av| av视频免费观看在线观看| 狂野欧美激情性xxxx| 波多野结衣av一区二区av| e午夜精品久久久久久久| 18禁观看日本| 国产熟女午夜一区二区三区| 婷婷色av中文字幕| 日本一区二区免费在线视频| 丝袜在线中文字幕| 只有这里有精品99| 亚洲欧美精品综合一区二区三区| 免费观看av网站的网址| 少妇猛男粗大的猛烈进出视频| www.999成人在线观看| 制服人妻中文乱码| 亚洲国产精品999| 亚洲 欧美一区二区三区| 国产在线视频一区二区| 亚洲精品av麻豆狂野| 人妻一区二区av| 国产视频首页在线观看| 国产精品香港三级国产av潘金莲 | 午夜福利影视在线免费观看| 日韩熟女老妇一区二区性免费视频| 高潮久久久久久久久久久不卡| 宅男免费午夜| 欧美日韩亚洲国产一区二区在线观看 | 精品一区二区三区四区五区乱码 | 男女边摸边吃奶| 黄色视频在线播放观看不卡| 汤姆久久久久久久影院中文字幕| 亚洲av综合色区一区| 日本wwww免费看| 天天影视国产精品| 国产麻豆69| 少妇精品久久久久久久| 黄色片一级片一级黄色片| 精品亚洲成a人片在线观看| 十分钟在线观看高清视频www| av欧美777| 国产日韩一区二区三区精品不卡| 国产欧美亚洲国产| 欧美人与性动交α欧美软件| 国产精品 国内视频| 激情视频va一区二区三区| 日韩伦理黄色片| 午夜影院在线不卡| 成人18禁高潮啪啪吃奶动态图| 国产在线免费精品| 婷婷色综合www| 黄频高清免费视频| 亚洲精品日本国产第一区| 久久久亚洲精品成人影院| 又紧又爽又黄一区二区| 欧美亚洲 丝袜 人妻 在线| 99热全是精品| 欧美日韩视频精品一区| 免费在线观看影片大全网站 | 精品久久久久久久毛片微露脸 | 制服人妻中文乱码| 亚洲国产日韩一区二区| 国语对白做爰xxxⅹ性视频网站| 亚洲av欧美aⅴ国产| 性色av一级| 精品一区二区三区四区五区乱码 | 婷婷色麻豆天堂久久| 在线亚洲精品国产二区图片欧美| 七月丁香在线播放| 美女午夜性视频免费| 久久午夜综合久久蜜桃| 欧美日本中文国产一区发布| 王馨瑶露胸无遮挡在线观看| 国产免费现黄频在线看| 一级毛片我不卡| 亚洲av日韩在线播放| 女性被躁到高潮视频| 日韩av不卡免费在线播放| 国产黄色视频一区二区在线观看| 日韩一本色道免费dvd| 精品国产超薄肉色丝袜足j| 99热国产这里只有精品6| 国产成人系列免费观看| 男女高潮啪啪啪动态图| 亚洲精品美女久久久久99蜜臀 | 久久女婷五月综合色啪小说| 热99国产精品久久久久久7| 久久精品成人免费网站| 男人舔女人的私密视频| 久久久久国产一级毛片高清牌| 亚洲,欧美精品.| 真人做人爱边吃奶动态|