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

    Dynamic and adaptive membrane lipid remodeling in leaves of sorghum under salt stress

    2022-12-02 01:00:32ShengnnGeDebinLiuMinChuXinyuLiuYuleiWeiXinyngCheLeiZhuLinHeJingyuXu
    The Crop Journal 2022年6期

    Shengnn Ge,Debin Liu,Min Chu,Xinyu Liu,Yulei Wei,Xinyng Che,Lei Zhu,Lin He,*,Jingyu Xu,b,*

    a Key Laboratory of Modern Agricultural Cultivation and Crop Germplasm Improvement of Heilongjiang Province/Key Laboratory of Low Carbon Green Agriculture in Northeast Plain,Ministry of Agriculture and Rural Affairs,China/Research Center of Saline and Alkali Land Improvement Engineering Technology in Heilongjiang Province/College of Agriculture,Heilongjiang Bayi Agricultural University,Daqing 163319,Heilongjiang,China

    b National Coarse Cereal Engineering Research Center,Daqing 163319,Heilongjiang,China

    Keywords:Sorghum(Sorghum bicolor L.Moench)Membrane lipid remodeling Lipidome Transcriptome Salt stress

    ABSTRACT Salt stress limits plant growth and development.In this study,changes in membrane lipids were investigated in leaves of sorghum seedlings subjected to salt stress(150 mmol L-1 NaCl).Galactolipids(DGDG and MGDG)accounted for more than 65%of the total glycerolipids in sorghum leaves.The predominance of C36 molecular species in MGDG suggested that sorghum is an 18:3 plant.Under NaCl treatment,the content of major phospholipids(PC and PE)increased,accompanied by the activation of their metabolism pathways at the transcriptional level.In contrast,the proportion of MGDG and PG dropped drastically,leading to a decreased ratio of plastidic to non-plastidic lipids.An adjustment of glycerolipid pathway between the cytosolic and plastidic compartments was triggered by salt stress,as reflected by the increased conversion of PC to PA,providing precursors for galactolipid synthesis.The elevated DGDG resulted in increased DGDG/MGDG and bilayer/non-bilayer lipid ratios.The double-bond index of PC,PE,and DGDG increased markedly,evidently owing to the increased expression of FAD3 and FAD8.These findings will be helpful for understanding dynamic membrane lipid changes and adaptive lipid remodeling in sorghum response to salt stress.

    1.Introduction

    Soil salinization is one of the main environmental stresses that limit plant growth and production,and is an increasingly severe global problem[1].More than 800 million hectares of land worldwide are affected by salt stress,and this situation will continue to worsen[2].

    Plant membrane lipids are classified into three main categories:sphingolipids,sterols,and glycerolipids[3].Glycerolipids are the most abundant and the main structural components of cell membranes[4].Membrane glycerolipids bind two fatty acids in the glycerol backbone at the first and second acyl(sn-1 and-2)positions,and condense with phosphorus(phospholipids)or sugar(glycolipid)molecules at the sn-3 site[5].In plant cells,the major phospholipids are phosphatidylcholine(PC)and phosphatidylethanolamine(PE),together with phosphatidylserine(PS),phosphatidylinositol(PI)and phosphatidic acid(PA),constituting the extra-plastidic membrane lipids[6].In plastids,monogalactosediacylglycerol (MGDG), digalactosyldiacylglycerol(DGDG),sulfoquinovosyldiacylglycerol(SQDG),and the only plastid-evolved phospholipid phosphatidylglycerol(PG),form the plastidic bilayer membrane[7].In the green tissues of plants,the dominant membrane lipids are MGDG and DGDG,which strongly influence plant chloroplast(and thylakoid)structure and photosynthetic characteristics[8].

    In plant cells,the synthesis of glycerolipids can be completed via two distinct pathways compartmentalized in cytoplasmic or plastids.The pathway located in the cytoplasmic/endoplasmic reticulum(ER)is called the eukaryotic pathway;the other is implemented in plastids and is referred to as the prokaryotic pathway[9].In the eukaryotic pathway,galactolipids carry C18 fatty acids at the sn-2 position owing the specificity of the acyltransferase,whereas prokaryoticpathway-generated galactolipids carry only C16 fatty acids at the sn-2 position[9].These two pathways work differently among plant species.In Arabidopsis,the synthesis of galactolipids is jointly accomplished by the eukaryotic and prokaryotic pathways,and the major product MGDG contains similar proportions of C34 and C36 fatty acid molecular species,so that this plant is called a 16:3 plant[10].In some plant species,known as 18:3 plants[11],the synthesis of galactolipids depends almost entirely on the eukaryotic pathway,leading to a predominance of C36:6 MGDG molecules.

    Despite the differential spatial distribution of these two pathways,they respond synergistically to various environmental factors[12-14].In the 16:3 plant Arabidopsis thaliana,the chloroplast pathway is activated under low-temperature stress,whereas high temperature promotes the eukaryotic pathway[15].In the 18:3 plants wheat and maize,channeling of glycerolipids from ER to chloroplast is observed under low temperature,as reflected by increased turnover of PC to PA for galactolipid synthesis[13,15].Alkaline stress influences the adjustment of membrane lipid pathways by activating extra-plastidic PC hydrolysis and plastidic galactolipid synthesis[14].

    The altering of the distribution of acyl chains and the exchange of head groups between different membrane lipid classes is also known as lipid remodeling[5,6,15]and occurs frequently in response to external environmental stimuli.Membrane lipid remodeling has been shown in plants exposed to salinity,which affects the permeabilityof the cellmembrane[16].In wheat,PC andPE werenegatively correlated at high temperature,and the higher PC:PE ratio may indicate PE-to-PC conversion[17].Under cold stress,DGDG content increases and MGDG/DGDG ratio decreases,favoring an increase in cell membrane stability[18,19].Changes in membrane lipid composition in plants subjected to stress help maintain the stability and integrity of membrane structure[14,20,21].

    Previous membrane-lipid studies have focused mainly on changes in total lipids of the plasma membrane,and information about the comprehensive lipid remodeling and salt-induced changes at the gene transcriptional level is scarce[5].Sorghum(Sorghum bicolor L.Moench)is an annual C4 cereal providing food,feed and fuel,and is a highly adaptable crop with high water-use efficiency and stress tolerance[22,23].The objective of the present study was to perform a combined lipidomic and transcriptomic analysis to identify dynamic membrane lipid metabolic changes and investigate the regulation of lipid remodeling in sorghum leaves under salt stress.Investigation of the biological response mechanisms of crops to salt stress will help to ensure crop productivity and make effective use of saline land[24].

    2.Materials and methods

    2.1.Plant growth,treatments,and sampling

    Seeds of the sorghum cultivar Hei63027 were disinfected with 10% NaClO solution for 30 min,washed,and soaked in distilled water for 6-8 h.The soaked seeds were placed in germination boxes in an incubator set at 25 °C in the dark for germination.Sprouts with uniform shoot length were planted on perforated foam boards and incubated in 1/2 Hoagland nutrient solution(pH 5.5).The seedlings were placed in an artificial-climate chamber with the temperature set at 25±2°C,16 h/8 h(light/dark)photoperiodic cycle,and relative humidity about 65%.The 2-week-old seedlings were treated with 150 mmol L-1NaCl.The leaf tissues were sampled for RNA sequencing at 0,12,24,and 72 h after NaCl treatment,and three biological replicates were collected for each treatment.For lipidomic analysis,leaf tissues were sampled at 0,3,7,and 14 d after NaCl treatment with six biological replicates for each group of samples.

    2.2.RNA-seq analysis and RT-PCR validation

    RNA from the leaves of sorghum seedlings was extracted with TRIZOL reagent,and a cDNA library was constructed after quality inspection.After the database qualification,RNA-seq sequencing was performed and clean reads were further annotated with the KEGG(Kyoto Encyclopedia of Genes and Genomes),Nr(NCBI non-redundant protein sequences),Nt(NCBI non-redundant nucleotide sequences),COG(Clusters of Orthologous Groups of proteins),and Swiss-port(a manually annotated and reviewed protein sequence database)databases.Identificaion of differentially expressed genes(DEGs)employed a false-discovery rate<0.01 and|log2fold change|≥1.The RNA-seq data have been submitted to the online SRA(sequential read file)database with accession number SRP333181.

    To verify the reliability of the RNA-seq data, real-time quantitative RT-PCR amplification was performed for selected DEGs. The synthesis of cDNA and qRT-PCR was completed in 96-well plates with SYBR qRT-PCR Master Mix (Toyobo, Osaka, Japan) [25]. Sb18SRNA and SbACTIN were used as internal reference genes. The results of qRTPCR (Fig. S1) were quantified by the 2-ΔΔCTmethod [26].

    2.3.Membrane lipid extraction and analysis

    Extraction of polar lipids from sorghum seedlings was performed by a modification of Narayanan’s method[27].Three mL of isopropanol solution containing 0.01% BHT was added to a 50 mL glass tube with a screw cap and preheated to 75 °C,and fresh leaves of sorghum seedlings were quickly placed in the preheated glass tube and heated for 15 min.To the tube were added 1.5 mL of chloroform and 0.6 mL distilled water,followed by shaking at room temperature for 1 h.The extract was transferred to a new glass tube,and 4 mL of chloroform:methanol(2:1)containing 0.01% BHT was added.The mixture was vortexed and then shaken for 30 min.This extraction was repeated 3-4 times until the leaves turned white.Then 1 mL of 1 mol L-1KCl solution was added to the extraction mixture and the supernatant remaining after the vortexing,shaking,and centrifugation was discarded.The lipid extract was washed with 2 mL of distilled water and evaporated in a nitrogen blowing instrument,redissolved in 1 mL of chloroform solution,transferred to a 2-mL brown bottle with a Teflon-lined screw cap,and evaporated again under nitrogen stream.The nitrogen-blown brown vials were stored in a-20 °C refrigerator.Lipid content and composition were determined by electrospray ionization mass spectrometry(ESI-MS/MS)at Kansas Lipidomics Research Center(Manhattan,KS,USA).

    2.4.Principal component analysis

    Principal component analysis(PCA)was performed with Graph-Pad Prism 9.0(https://www.graphpad-prism.cn).Six biological replicates were used for lipidomic data analysis of lipid species:galactolipids(MGDG and DGDG),phospholipids(PC,PE,PA,PI,PS,and PG)and lysophospholipids(LPG,LPC,and LPA).

    2.5.Statistical analysis

    Statistical analysis was performed using GraphPad Prism,Microsoft Excel and SPSS Statistics,and the significance level was calculated by Student’s t-test.

    3.Results

    3.1.Changes of membrane lipids in sorghum leaves under salt stress

    As revealed in Fig.1A,the molar percentages of galactolipids MGDG and DGDG were dominant in the leaves of sorghum seedlings.MGDG and DGDG together accounted for more than 65% of total glycerolipids,with MGDG accounting for 45% and DGDG for 20%.Among the six types of phospholipids,PG and PC were present in higher molar percentages,accounted for respectively 14%and 10%of total glycerolipids.The molar percentages of PI,PS,PA,and the three lysophospholipids were relatively low.Under salt stress,the content of MGDG decreased significantly with the extension of NaCl treatment,and dropped by 47.3% by the 14-d time point;whereas the molar percentage of DGDG increased under long-term NaCl treatment to 29%at 14 d.In the phospholipid classes,the molar percentage of plastidic phospholipid PG decreased by respectively 11.5%and 47.5% after 7 and 14 d treatment.The molar percentages of PC and PE increased drastically at the later stage of NaCl treatment,and rose more than two folds at the 14-d time point.A significant increase of LPC and LPE was observed on the 14-d time point.

    PCA showed that under salt stress,the glycerolipid species in sorghum leaves were clearly separated by PC1 and PC2,which explained respectively 69.5% and 15.8% of total variance.MGDG,PG,PA,and LPG constituted the main components of PC1,while DGDG,PE,and PI constituted the main components of PC2(Fig.1B).

    Spearman’s correlation coefficients(ρ)between the types of glycerolipids were calculated.As shown in Fig.1C,DGDG was negatively correlated with MGDG and PG,with correlations of-0.90 and-0.93,respectively.DGDG was positively correlated with most phospholipid species with correlations ranging from 0.84 to 0.89.MGDG was positively correlated with PG,and strongly negatively correlated with other phospholipids.PG was positively correlated with PA with a coefficient of 0.21,but showed a strong negative correlation with PC,PE,PI,and PS.Most extra-plastidic phospholipids and lysophospholipids were positively correlated,except that PA was negatively correlated with most other phospholipids.

    The proportions of specific glycerolipids also underwent adjustment under salt stress.As shown in Fig.1D,the ratio of DGDG to MGDG increased significantly under salt stress,especially at the 14-d time point;whereas the ratio of PC:PE decreased slightly but not significantly.Bilayer lipids,including DGDG,PC,PG,and PI,are more beneficial to the stability of lipid bilayers than are nonbilayer lipids,including MGDG,PE,PA,and PS.The ratio of bilayer lipid:non-bilayer lipid was significantly increaseed at 14 d,suggesting an increase in membrane stability in response to salt stress.The ratio of plastidic membrane lipids(MGDG,DGDG,and PG)to nonplastidic lipids(phospholipids other than PG)declined under salt stress,particularly at the later stage of NaCl treatment,suggesting that salt stress caused damage to the photosynthetic membrane.

    3.2.Changes in molecular species of phospholipids in sorghum leaves under salt stress

    The fatty acid composition at sn-1&2 acyl positions of individual glycerolipids was determined by lipidomic analysis and marked as Cxx:x(total number of acyl carbon atoms:number of double bonds).As shown in Fig.2A,the six classes of phospholipids were composed of different molecular species dominated by C34 and C36 molecules,ranging from C34:1 to C36:6.Given that PG is the only phospholipid assembled in the plastids,the C34 molecular species was dominant in PG.For other non-plastidic phospholipids,the proportion of C34 and C36 molecules was similar(Fig.2A).Under salt stress,the molecular species of specific phospholipids underwent alterations.As shown in Fig.2B,in PG(plastidic phospholipid),the polyunsaturated molecular species C34:3 and C34:4 decreased significantly under long-term NaCl treatment,particularly at 14 d.In contrast,the C36:6 PA increased drastically at 14 d.The most noteworthy finding was the changes in molecular species of the major phospholipid PC.All C34 and C36 species in PC showed an upward trend in the later stage of salt stress treatment,except for C36:4,the only species that decreased significantly under salt stress(Fig.2B).

    3.3.Molecular species analysis of galactolipids revealed that sorghum is an 18:3 plant species

    As shown in Fig.3A,the molar percentage of C34:3,C36:6,C36:5,and C36:4 MGDG gradually decreased under salt stress and dropped significantly at the 14-d time point.However,the contents of C36:6,C36:5,and C36:3 DGDG gradually increased under salt stress(Fig.3A).

    As revealed in Fig.3B,the proportion of C36 MGDG molecules accounted for most of the total MGDG molecules,with C36:6 MGDG accounting for 91.9% of total MGDGs.C36:6 DGDG accounted for 74.9% of total DGDGs,whereas C34:3 accounted for only 14.8%of total DGDGs(Fig.3B).According to the difference of lipid synthesis pathways[9,28],there is a class of plants called 18:3 plants,in which galactolipid synthesis relies on the ER pathway and produces C36-dominated MGDG and DGDG.In the present study,the predominance of C36 galactolipids implied dependence on ER pathway,suggesting that sorghum is an 18:3 plant.

    C34:2 MGDG and C34:3 MGDG were positively correlated,with a correlation of 0.83,whereas C36:4 MGDG and C36:6 MGDG had a correlation of 0.89.C34:2 DGDG was positively correlated with C34:3 DGDG,and the correlation coefficient was 0.10;C36:4 DGDG was negatively correlated with C36:6 DGDG,and the coefficient was-0.17(Fig.3C).

    3.4.Correlation and unsaturation index analysis of membrane lipids in sorghum leaves under salt stress

    As shown in Fig.4,correlation analysis of phospholipids(PC,PE)and galactolipids(MGDG,DGDG)revealed that phospholipids(PC,PE)were negatively correlated with MGDG and positively correlated with DGDG,but that C36:4 and C34:2 PC were positively correlated with most MGDG and negatively correlated with most DGDG.The correlation between the molecular species of phospholipids PC and PE was calculated and it was found that there was an overall positive correlation among phospholipids,but C34:2 PC and C36:4 PC were negatively correlated with other molecular species of PC and PE,with C36:4 PC showing the largest negative correlation.There was no apparent conversion between PC and PE.There was an overall negative correlation between galactolipids,with stronger negative correlations among C36:2,C36:3,C36:5,and C36:6 DGDG and MGDG molecular species,suggesting a strong conversion relationship.These results are consistent with the previous ratio analysis of classified lipid species(Fig.1B).

    Fatty acid unsaturation index is an indicator of the characteristics of plasma membrane,and strongly influences the adaptation of plants to abiotic stresses.In this study,the double-bond index(DBI)was calculated(Fig.5A),and the DBIs of galactolipids MGDG and DGDG were relatively large.Under salt stress,the DBI of DGDG increased with the prolongation of NaCl treatment,while the DBI of MGDG decreased significantly under salt treatment.In the phospholipids,the DBI of PC and PE increased significantly,especially at the 14-d time point,whereas the unsaturation of PG decreased gradually under salt stress.

    Fig.1.Changes of membrane glycerolipids in sorghum leaves under salt stress.(A)Changes in glycerolipid content;values(mol%)are mean±standard error(n=6);the significance level was calculated by Student’s t-test;*indicates that the value is significantly different from the control(P<0.05);**indicates that the value is extremely different from the control.(B)Principal component analysis of lipid groups.(C)Correlations among lipid species.(D)Ratio comparison among lipid species.The a,b,and c in figure D indicate the significance of the difference between contrasted groups.

    3.5.Transcriptomic analysis of sorghum leaves under salt stress

    Transcriptomic analysis was performed on RNA prepared from sorghum leaf samples collected at 0,12,24,and 72 h after 150 mmol L-1NaCl treatment.A total of 84.44 Gb clean data was obtained,and the clean reads were annotated with KEGG(Kyoto Encyclopedia of Genes and Genomes),Nr(NCBI non-redundant protein sequences),Nt(NCBI non-redundant nucleotide sequences),COG(Clusters of Orthologous Groups of proteins),and Swiss-port(a manually annotated and reviewed protein sequence database)databases.The accession ID of the RNA-seq data in the SRA database is SRP333181.DEGs were identified and a total of 384 co-expressed DEGs(|log2FC|≥1)were retrieved in the‘‘12 h vs.0 h”,‘‘24 h vs.0 h”,and‘‘72 h vs.0 h”comparison groups,of which 140 were up-regulated and 244 down-regulated(Fig.6A).In the‘‘24 h vs.12 h”,‘‘72 h vs.12 h”and‘‘72 h vs.24 h”comparison groups,there were 41 co-expressed genes,of which 30 were up-regulated and 11 down-regulated.There were 41 differentially co-expressed genes in the‘‘24 h vs.12 h”and‘‘72 h vs.24 h”comparison groups.There were 249 differentially co-expressed genes in‘‘72 h vs.12 h”and‘‘72 h vs.24 h”comparison groups,of which 177 were up-regulated and 72 were downregulated(Fig.6A).These genes may function in salt stress response.

    The KEGG annotation analysis of DEGs is shown in Fig.6B,‘‘Metabolism”accounted for the main proportion(around 80%)in all three comparison groups.Partitioning analysis of‘‘Metabolism”showed that‘‘Lipid metabolism”accounted for a large proportion in all the metabolic processes.The‘‘Others”proportion was split into several parts,and‘‘Environmental information processing”accounted for a larger proportion in the‘‘24 h vs.0 h”comparison group;whereas in the‘‘72 h vs.0 h”comparison group,‘‘Genetic information processing”accounted for more than half of that part.

    3.6.Differentially expressed genes associated with lipid metabolism under salt stress

    Genes associated with lipid synthesis were identified in the sorghum transcriptome based on the annotations of the GO and KEGG databases and the Arabidopsis lipid gene database[29,30].As shown in Fig.7A,a total of 395 genes associated with lipid metabolism were identified in the‘‘12 h vs.0 h”comparison group,of which 177 were up-regulated and 218 down-regulated.Respectively 379(207 up-regulated and 172 down-regulated),and 387(198 up-regulated and 189 down-regulated)lipid-related genes were identified in the‘‘24 h vs.0 h”and‘‘72 h vs.0 h”comparison groups.The increase of up-regulated genes in the later stage of salt stress indicates that some lipid metabolism pathways had been activated.The up-regulated genes were enriched mainly in the‘‘Phospholipid signaling”,‘‘Fatty acid elongation & wax biosynthesis”pathways,and the number of up-regulated genes reached 30 in the‘‘24 h vs.0 h”and‘‘72 h vs.0 h”comparison groups.A higher proportion of up-regulated genes were also found in the‘‘Triacylglycerol synthesis”,‘‘Triacylglycerol & fatty acid degradation”,‘‘Oxilipin metabolism”,‘‘Eukaryotic phospholipid synthesis& editing”,and‘‘Eukaryotic galactolipid&sulfolipid synthesis”processes.

    Fig.2.Changes in molecular species of phospholipids in sorghum leaves under salt stress.(A)Changes in molecular species of phospholipids;the significance level was calculated by Student’s t-test;‘‘0 d”is set as control;*indicates that the value is significantly different from the control(P<0.05);**indicates that the value is extremely different from the control.(B)Differential changes of specific molecular species in PC and PG.The a,b,and c in the figure B indicate the significance of the difference between the contrasted groups.

    Fig.3.Changes in galactolipids in sorghum leaves under salt stress.(A)Changes in molecular species of galactolipids.(B)The occupancy of galactolipid molecular species.(C)Correlations among galactolipid molecular species.Values(mol%)are mean value of replicates±standard error(n=6).‘‘0 d”is set as control;*indicates that the value is significantly different from the control(P<0.05);**indicates that the value is extremely different from the control.

    Fig.4.Correlation among molecular species of PC,PE,MGDG,and DGDG.The Spearman correlation was calculated basing on the molar percentage content of C34-C36 PC,PE,MGDG,and DGDG.Spearman’s correlations<0 are negative and represented by green squares and those greater than 0 are positive and represented by red squares.

    A further proportional analysis of specific lipid metabolism pathways also revealed that the‘‘Fatty acid elongation & wax biosynthesis”and‘‘Phospholipid signaling”pathways accounted for a large proportion of total genes associated with lipid synthesis(Fig.7B).These results together indicated that the‘‘Phospholipid signaling”and‘‘Fatty acid elongation & wax biosynthesis”pathways were activated under salt stress treatment.

    The significant DEGs in major lipid metabolic pathways are presented in heat maps in Fig.7C.The expression of genes encoding the enzymes involved in the steps of de novo TAG biosynthesis,which are also the essential steps for PC biosynthesis,including GPAT,LPAAT,PAH,DGAT1,and DGAT3,were markedly upregulated under salt stress treatment.Genes involved in the de novo synthesis of PC were influenced by salt stress as well,including down-regulation of CEK(choline kinase),up-regulation of CCT(choline-phosphate cytidylyltransferase)and AAPT(aminoalcohol phosphotransferase).These findings indicated that the activated DAG-PC and de novo PC assembly steps increased PC synthesis,in agreement with the increase in PC content revealed by the lipidomic analysis.As the central intermediate in lipid metabolism,PC undergoes active degradation to provide DAG and PA as precursors for synthesis of other lipids.In the present study,7 of 8 PLD isoforms and 3 of 4 NPC were significantly up-regulated at most time points under salt stress.The expression of genes involved in the biosynthesis of galactolipids(MGDG and DGDG)was not obviously up-regulated under salt stress.The expression of 2 MGD isoforms was down-regulated,whereas the expression of 2 DGD isoforms was slightly up-regulated at different time points of salt treatment.The expression of DGD1 was up-regulated in the‘‘24 h vs.0 h”comparison group,while DGD2 was up-regulated in the‘‘12 h vs.0 h”comparison group,possibly reflecting an early response of DGD genes to salt stress.The genes encoding FAD3 and FAD8.2 were up-regulated,whereas the other FAD genes were mostly downregulated under salt stress.

    3.7.Regulation of lipid metabolism in sorghum leaves under salt stress

    Based on our findings from the transcriptome and lipidome data,a gene-metabolite regulation network in sorghum leaves under salt stress was constructed(Fig.8).The steps or pathways of transcriptional activation under salt stress are represented by red arrows.As shown in Fig.8,the activated steps or pathways were associated mainly with PC synthesis and turnover,galactolipid synthesis and degradation,and desaturation of fatty acids.

    Fig.5.Changes in double-bond index(DBI)in sorghum leaves under salt stress.(A)Clustered DBI analysis of glycerolipids.(B)DBI of individual glycerolipids.DBI=(∑[N×mol%lipid])/100,N is the total number of double bonds in the two fatty acid chains of each glycerolipid molecule.The significance level was calculated by Student’s ttest;a,b,and c in the figure B indicate the significance of the difference between the contrasted groups.

    In ER compartment,GPAT,LPAAT,PAH,and DGAT are the major catalytic enzymes involved in the biological assembly of TAG,and most of the genes encoding these enzymes were up-regulated under salt stress.The up-regulated expression of PAH and AAPT indicated that the PA-DAG-PC pathway was activated for PC production.The de novo PC synthesis pathway through CCT was also induced under salt stress.These changes together suggest increased PC production in response to salt stress at the transcriptional level,evidenced by elevated PC accumulation at the later stage of salt stress as revealed by lipidomic analysis.PC is a relatively abundant phospholipid and an intermediate in lipid metabolism.The ER-derived PC is converted by PLD to PA and finally by NPC to form DAG,which can be used for synthesis of MGDG and DGDG.Most genes encoding PLD and NPC were up-regulated under salt stress,but this up-regulation did not lead to a decrease in total PC level.However,the accumulation of one specific PC molecule,C36:4 PC,dropped greatly under salt stress,suggesting that it might be the favored substrate for the activated hydrolases PLD and NPC.The increased accumulation of polyunsaturated PCs,C34:3 and C36:6 PC,under salt stress could be attributed to the induced expression of the desaturase FAD3.

    In the chloroplast compartment,although the synthesis of galactolipids seemed to be inhibited at both the transcriptional and metabolic levels under salt stress,the steps governing the synthesis and desaturation of DGDG were induced as illustrated by the up-regulated DGD2 and FAD8 genes,which led to the increased DGDG accumulation.Up-regulation of DGL genes in the chloroplast might result in galactolipid degradation,converting MGDG and DGDG into free fatty acids(FFA).The down-regulation of genes involved in PG synthesis pathway led to lowered PG production under salt stress(Fig.8).The impaired accumulation of plastidic lipids in sorghum leaves may give rise to dysfunction of the photosynthetic membrane of the chloroplast.

    4.Discussion

    The cell membrane is a barrier between plants and external stimuli,owing to its specific structure,typical permeability,and other protective properties[31].Glycerolipids are essential components of the plasma membrane,and large metabolic changes and lipid remodeling occur when plants are subjected to environmental stress[13,25,32,33].

    Phospholipids are the main building blocks of biological membranes,with PC and PE the major membrane phospholipids in plant cells[34].In the present study,the total molar percentage of phospholipids was elevated in sorghum leaves under salt stress,with the molar percentage content of PC and PE increased about twofold at 14 d relative to 0 d.The PC and PE synthesis pathways were activated at the transcriptional level by salt treatment,as reflected by the up-regulation of genes involved in those pathways.In a previous study[35],PC content in chloroplast membranes of leaves and bundle sheaths of maize increased under salt stress.

    Fig.6.Differentially expressed gene(DEG)distribution and Kyoto Encyclopedia of Genes and Genomes(KEGG)classification in the transcriptome of sorghum leaves under salt stress.(A)Venn diagram showing DEGs;colored shapes indicate the expression value.Pink indicates high level,and green indicates low level.(B)KEGG annotation of genes associated with lipid metabolism.

    Membrane galactolipids(MGDG and DGDG)are dominant lipids in photosynthetic membranes of plant leaves and are stronlgy affected by abiotic stresses[36].MGDG is the most abundant galactolipid in the thylakoid membrane and influences chloroplast structure and photosynthetic performance[37].In the present study,the galactolipids in sorghum leaves decreased by 26.6% under salt stress.Although the content of DGDG increased,MGDG decreased by 47.3% at 14 d compared with 0 d.Marked phenotypic changes of sorghum seedlings under salt stress were also observed.The plant growth was severely restricted,and the leaves of sorghum seedlings gradually turned yellow and wilted with increasing duration of salt stress,indicating that the photosynthetic system of the plant was impaired(Fig.S2).The reduction of glycolipids in the leaves of rice and cucumber seedlings results in a decrease in chlorophyll content in the leaves[37,38].Reduction in MGDG seems to be a general response of plants to osmotic stress caused by drought,salt or freezing damage[39].Drought caused a decrease in MGDG and an increase in DGDG in the leaves of two maize varieties[40].Synthesis of DGDG was triggered by the up-regulation of DGD2 under salt-alkali stress[14].The transcriptome analysis in the present study revealed the presence of down-regulated MGD and up-regulated DGD isoforms,possibly accounting for the differing changes in MGDG and DGDG.

    Fig.7.Functional classification of DEGs associated with lipid metabolism in sorghum leaves under salt stress.(A)Functional categorization of lipid-related DEGs involved in various lipid metabolic pathways.In each category,purple columns represent up-regulated genes and green columns down-regulated genes.The number of genes in each category is shown on the x axis.(B)The proportional distribution of lipid-related DEGs.(C)The significant salt response lipid-related DEGs.Heat maps were constructed to illustrate the differential expression profiles of important lipid-associated DEGs.Numbers in each color block represent log2(fold-change)of the corresponding genes,and negative numbers represent down-regulated DEGs.Color scale is provided.

    In plant cells,phospholipids PC,PE,PI,and PS constitute nonplastidic membrane lipids,whereas the glycolipids MGDG,DGDG,SQDG,and PG constitute the plastidic membrane lipids[6,7].In the green tissues of plants,plastidic lipids are predominant and associated closely with plant thylakoid membrane structure and photosynthetic characteristics[8].In the present study,the ratio of plastidic to non-plastidic membrane lipids dropped to a very low level at 14 d salt treatment.It could be speculated that the sharp reduction in the level of galactolipids and the lowered ratio of plastidic to non-plastidic lipids lead to damage to photosynthetic membranes,and result in dysfunction of the chloroplast and photosynthesis in sorghum leaves under salt stress(Fig.9).

    Membrane lipids are usually classified as bilayer or non-bilayer lipids[41].Whereas bilayer lipids ensure membrane stability,nonbilayer lipids function in mediating protein interactions and increasing the plasticity of lipid bilayer structures[42].The bilayer membrane lipids include PC,PG,PI,DGDG,and SQDG;while nonbilayer membrane lipids include PA,PE,PS,and MGDG[6].In the present study,the ratio of bilayer to non-bilayer lipid in the leaves of sorghum seedlings increased by 65.2% under salt stress.The major bilayer lipid PC increased significantly,whereas the ratio of PC/PE remains basically unchanged.In contrast,the ratio of DGDG/MGDG increased drastically at 14 d(145%),suggesting the contribution of bilayer lipid DGDG to the maintenance of chloroplast membrane stability under salt stress(Fig.9).There is evidence[43]that the DGDG/MGDG ratio is necessary for the stability of photosynthetic membranes.In maize,varieties with delayed leaf senescence under drought conditions showed higher DGDG content and DGDG/MGDG ratio[40].

    In the present study,the C36:6 species in MGDG and DGDG in the leaves of sorghum seedlings accounted for the largest proportion,indicating that sorghum is an 18:3 plant.In 18:3 plants,the main synthetic precursors of MGDG and DGDG come from the degradation of eukaryotic-pathway-evolved PC,and the flux of PC hydrolyzates tends to be increased under various stress conditions[15,44].Under salt stress,although the total amount of PC did not show a large decrease,C36:4 PC was the only molecule showing a large reduction among all molecular species,suggesting that C36:4 was the most favorable precursor for glycolipid synthesis.In agreement with this finding,previous studies[13,15]have shown that C36:4 PC is the specific molecular species contributing to lipid flux shared between the cytoplasm and plastid pathways when plants encounter temperature stress.Our findings indicated that under salt stress,C36:4 PC was also the major contributor to plastidic lipid synthesis.Phospholipase D(PLD)and non-specific phospholipase C(NPC)catalyze the degradation of PC into PA and DAG,respectively,and participate in lipid remodeling[45].In the transcriptome data,genes involved in the PLD and PLC/DGK pathways were up-regulated under salt stress,implying activated PC degradation to ensure the flow of lipid intermediates from cytoplasm to chloroplast for photosynthetic membrane lipid assembly.

    Fig.8.The gene-metabolite network of lipid metabolism in sorghum leaves under salt stress.The relative changes of lipid molecular species[(3 d or 7 d or 14 d mol%-0 d mol%)/0 d mol%]and the relative expression levels of related genes(72 h vs.0 h)are marked as heat-map icons.Changes in lipid content are represented by diamond-shaped heat map icons,and transcription levels are represented by square heat-map icons.The steps or pathways of transcriptional activation under salt stress are represented by red arrows.

    Fatty acid unsaturation index is a major indicator of the characteristics of the plasma membrane,and is also closely associated with adaptation of plants to abiotic stresses[46,47].In the present study,the DBI of the bilayer membrane lipids(DGDG,PC,PI)underwent significant elevation.The DBI of DGDG increased by 30.5%,while DBI of PC and PI increased by 144.0% and 171.8%,respectively.Thus,it could be speculated that the fluidity of the bilayer membrane increased under salt stress.In Phanerochaete rose,the DBI of all glycerophospholipids increased after cold treatment[48].In Fabaceae,200 mmol L-1NaCl treatment led to increased unsaturation of MGDG,DGDG,and PG,which might be beneficial to salt stress adaption and protection of the integrity of the chloroplast membrane[37].Fatty acid unsaturation is determined by the activity of fatty acid desaturase(FAD),which introduces double bonds at specific acyl chain positions[32,49].In the ER compartment,FAD2 converts oleic acid(18:1)to linoleic acid(18:2),which is further desaturated by FAD3 to linolenic acid(18:3)[50,51];in plastids,FAD7 and FAD8 are responsible for catalyzing the 18:2 to 18:3 conversion[52].In our transcriptome data,the genes encoding FAD3 and FAD8.2 were up-regulated,while other FAD genes were mostly down-regulated under salt stress,indicating that the 18:2 to 18:3 conversion was promoted under salt stress.In transgenic tobacco plants,overexpression ofω-3 desaturase(FAD3 or FAD8)increased tolerance to salt and drought stress[53].Knockout of FAD8 led to reduced membrane fluidity in rice under cold stress[54].

    Although TAG is not usually produced at high levels in vegetative tissues of plants,its accumulation in these tissues functions in plant response to abiotic stresses[55].TAG accounts for<1% of total lipid content in Arabidopsis leaf tissue,but it is related to the short-term supply of glycerolipid intermediates during membrane lipid remodeling[15,56].TAG regulates the conversion of excessive DAG to TAG through the catalysis of DGAT and DGK,so that plants can resist freezing stress[57,58].In the present study,the expression of genes encoding the enzymes involved in the steps of de novo TAG biosynthesis,including DGAT1 and DGAT3,was up-regulated under salt stress treatment.Although we did not investigate changes in TAG in this study,it can be inferred from the gene transcription that TAG is essential in membrane-lipid metabolic response to salt stress.

    Fig.9.Schematic diagram of membrane lipid changes and adaptive remodeling in response to salt stress in sorghum leaves.A decrease in the proportion of plastidic lipids in sorghum leaves may severely damage the stability of photosynthetic membranes and the function of the chloroplast.Synergistically,several adaptive membrane lipid adjustments also occur,including increased conversion of C36:4 PC into C36:6 MGDG and C36:6 DGDG,increased ratio of bilayer lipid/non-bilayer lipid,and an increased double-bond index(DBI)in the major phospholipids and DGDG,which may increase membrane stability and fluidity in response to salt stress.

    Membrane lipid remolding is a necessary modulation for plant to cope with abiotic stresses.In the present study,active metabolic changes in membrane lipids were observed in sorghum seedlings under salt stress.As summarized in Fig.9,salt stress led to a significant decrease in the proportion of plastidic lipids in sorghum leaves,particularly MGDG and PG,which may greatly damage the stability of photosynthetic membrane and the function of chloroplast.Synergistically,several adaptive membrane lipid adjustments also occur for sorghum leaves to adapt to external stimuli.First,the adjustment of the glycerolipid pathway between the cytosolic and plastidic compartment resulted in an intensive conversion of C36:4 PC into C36:6 MGDG and then into C36:6 DGDG to complement the chloroplast lipids and ensure the integrity of the photosynthetic membrane.Second,the ratio of bilayer lipid to non-bilayer lipid increased,especially the ratio of DGDG/MGDG,which is conducive to stabilizing the photosynthetic membrane.Thirdly,the elevation of the DBI in the major phospholipids and DGDG implied increased membrane fluidity in response to salt stress.These findings are helpful for understanding the dynamic regulation of membrane lipids in plants under salt stress and for identifying some adaptive roles of membrane-lipid remodeling in plant adaptation to salt stress.

    CRediT authorship contribution statement

    Shengnan Ge:Validation,Visualization,Writing-original draft.Debin Liu and Min Chu:Conceptualization,Writing-review &editing.Xinyu Liu and Yulei Wei:Data curation.Xinyang Che and Lei Zhu:Formal analysis.Lin He and Jingyu Xu:Conceptualization,Funding acquisition,Project administration,Writingreview & editing.

    Declaration of competing interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgments

    We gratefully acknowledge Dr.Ruth Welti and Mary Roth(Kansas Lipidomics Research Center)for their help with lipidomic analysis.This work was supported by the Natural Science Foundation of Heilongjiang Province(ZD2020C007,QC2017024),Heilongjiang Bayi Agricultural University Support Program for San Heng San Zong(ZDZX202101),and Special Funds from the Central Finance to Support the Development of Local Universities(TO SXL).

    Appendix A.Supplementary data

    Supplementary data for this article can be found online at https://doi.org/10.1016/j.cj.2022.03.006.

    久久人妻熟女aⅴ| 人人妻人人澡人人看| 国产国语露脸激情在线看| 可以免费在线观看a视频的电影网站 | 三级国产精品片| 午夜日韩欧美国产| 国产片特级美女逼逼视频| 一级黄片播放器| 亚洲精品国产av蜜桃| 精品亚洲成a人片在线观看| 国产亚洲av片在线观看秒播厂| 久久久久久人人人人人| 国产淫语在线视频| 亚洲精品美女久久av网站| 亚洲人成77777在线视频| 国产精品av久久久久免费| 大话2 男鬼变身卡| 99热国产这里只有精品6| 女人精品久久久久毛片| 少妇被粗大猛烈的视频| 高清av免费在线| 国产免费视频播放在线视频| 亚洲一区中文字幕在线| 99久久人妻综合| 久久久久国产网址| 亚洲熟女精品中文字幕| 中文字幕另类日韩欧美亚洲嫩草| 91aial.com中文字幕在线观看| www.熟女人妻精品国产| 97精品久久久久久久久久精品| 精品一区二区三卡| 欧美日韩综合久久久久久| 黄色一级大片看看| 久久久亚洲精品成人影院| 欧美日韩一级在线毛片| 欧美日韩一级在线毛片| 午夜免费鲁丝| 午夜激情av网站| 如何舔出高潮| 精品视频人人做人人爽| 免费看不卡的av| 国产男女超爽视频在线观看| 超色免费av| av在线观看视频网站免费| 国产一区二区在线观看av| 又粗又硬又长又爽又黄的视频| 午夜福利影视在线免费观看| 久久人人97超碰香蕉20202| 男女无遮挡免费网站观看| 丁香六月天网| 国产成人精品福利久久| 亚洲精品视频女| 精品国产乱码久久久久久小说| 日韩中文字幕视频在线看片| 久久99热这里只频精品6学生| 日本欧美国产在线视频| 丰满饥渴人妻一区二区三| 涩涩av久久男人的天堂| 黄色配什么色好看| 久久99精品国语久久久| 熟女少妇亚洲综合色aaa.| 免费观看性生交大片5| 女性被躁到高潮视频| 久久精品夜色国产| 欧美日韩亚洲高清精品| 亚洲精华国产精华液的使用体验| 男女免费视频国产| 哪个播放器可以免费观看大片| 黄色一级大片看看| 亚洲一区二区三区欧美精品| 18禁裸乳无遮挡动漫免费视频| 亚洲欧洲精品一区二区精品久久久 | 一区二区三区激情视频| 美女午夜性视频免费| 1024香蕉在线观看| 国产av码专区亚洲av| 汤姆久久久久久久影院中文字幕| 99久久中文字幕三级久久日本| www.熟女人妻精品国产| 欧美精品亚洲一区二区| 国产日韩欧美在线精品| 大香蕉久久网| 中文欧美无线码| 久久久久久伊人网av| √禁漫天堂资源中文www| 丰满少妇做爰视频| 欧美成人午夜免费资源| 性高湖久久久久久久久免费观看| 只有这里有精品99| av网站在线播放免费| 国精品久久久久久国模美| 欧美精品亚洲一区二区| 亚洲第一av免费看| 久久久久久久久久人人人人人人| 午夜av观看不卡| 大码成人一级视频| 欧美日韩国产mv在线观看视频| 九色亚洲精品在线播放| 日本欧美视频一区| kizo精华| www.av在线官网国产| 亚洲精品aⅴ在线观看| 成人国产麻豆网| 男女国产视频网站| 黄色视频在线播放观看不卡| 国产精品久久久av美女十八| 亚洲中文av在线| 日韩av不卡免费在线播放| 亚洲成国产人片在线观看| 91在线精品国自产拍蜜月| 亚洲欧美一区二区三区国产| 欧美在线黄色| 男女下面插进去视频免费观看| 免费在线观看黄色视频的| 国产日韩一区二区三区精品不卡| 黄网站色视频无遮挡免费观看| 久久精品aⅴ一区二区三区四区 | 国产亚洲精品第一综合不卡| 国产又色又爽无遮挡免| 亚洲精品久久久久久婷婷小说| 大码成人一级视频| 免费黄色在线免费观看| 欧美+日韩+精品| 男的添女的下面高潮视频| 国产精品亚洲av一区麻豆 | 国产免费视频播放在线视频| 欧美成人午夜精品| 亚洲精品美女久久久久99蜜臀 | av不卡在线播放| 男人舔女人的私密视频| 国产探花极品一区二区| 国产精品 国内视频| 中国国产av一级| 午夜精品国产一区二区电影| 大话2 男鬼变身卡| 搡老乐熟女国产| 久久久久久人人人人人| 成人毛片60女人毛片免费| 91久久精品国产一区二区三区| 日韩中文字幕欧美一区二区 | 久久久久久久亚洲中文字幕| 视频在线观看一区二区三区| 久久韩国三级中文字幕| 人妻 亚洲 视频| 欧美精品一区二区免费开放| 一级毛片我不卡| 日韩视频在线欧美| 汤姆久久久久久久影院中文字幕| 一本久久精品| 欧美激情 高清一区二区三区| 久久精品久久久久久噜噜老黄| 久久国产亚洲av麻豆专区| 一个人免费看片子| 久久99蜜桃精品久久| 一二三四中文在线观看免费高清| 精品亚洲成a人片在线观看| 美女大奶头黄色视频| 美女xxoo啪啪120秒动态图| 精品少妇内射三级| 国产xxxxx性猛交| 午夜日本视频在线| 王馨瑶露胸无遮挡在线观看| 日韩电影二区| 成人国产麻豆网| 在现免费观看毛片| 国产探花极品一区二区| 精品久久蜜臀av无| 欧美亚洲日本最大视频资源| 人人妻人人爽人人添夜夜欢视频| 午夜av观看不卡| 丝袜美足系列| 国产无遮挡羞羞视频在线观看| 99国产精品免费福利视频| 极品人妻少妇av视频| 男人爽女人下面视频在线观看| 18禁国产床啪视频网站| av免费在线看不卡| 另类精品久久| 午夜福利在线观看免费完整高清在| 久久久精品94久久精品| 一区二区三区四区激情视频| 极品少妇高潮喷水抽搐| 亚洲色图 男人天堂 中文字幕| videossex国产| 在线观看免费视频网站a站| 久久精品亚洲av国产电影网| 超碰成人久久| 国产精品.久久久| 久久久久国产一级毛片高清牌| 国产亚洲一区二区精品| 免费高清在线观看日韩| 久久ye,这里只有精品| 国产日韩欧美在线精品| 如日韩欧美国产精品一区二区三区| 王馨瑶露胸无遮挡在线观看| 国产毛片在线视频| 亚洲国产欧美在线一区| 老熟女久久久| 亚洲精品,欧美精品| 好男人视频免费观看在线| 一区福利在线观看| 日韩一卡2卡3卡4卡2021年| 国产精品久久久久久精品电影小说| 亚洲成av片中文字幕在线观看 | 男人舔女人的私密视频| 汤姆久久久久久久影院中文字幕| 老司机亚洲免费影院| 国产精品二区激情视频| 亚洲精品国产一区二区精华液| 成人二区视频| 国产女主播在线喷水免费视频网站| 日本爱情动作片www.在线观看| 97在线人人人人妻| 最黄视频免费看| 免费日韩欧美在线观看| 亚洲精品美女久久久久99蜜臀 | 亚洲伊人久久精品综合| 国产黄色视频一区二区在线观看| 午夜影院在线不卡| 丁香六月天网| 欧美人与性动交α欧美软件| 国产色婷婷99| 亚洲欧美精品自产自拍| 国产男人的电影天堂91| 久久久久精品性色| 国产日韩欧美在线精品| 你懂的网址亚洲精品在线观看| 日本猛色少妇xxxxx猛交久久| 亚洲精品日本国产第一区| 午夜福利乱码中文字幕| 精品福利永久在线观看| 十分钟在线观看高清视频www| 啦啦啦在线免费观看视频4| 妹子高潮喷水视频| 不卡视频在线观看欧美| 如日韩欧美国产精品一区二区三区| 制服人妻中文乱码| 老司机亚洲免费影院| 欧美av亚洲av综合av国产av | 久久鲁丝午夜福利片| 国产 一区精品| 精品一区在线观看国产| 中文乱码字字幕精品一区二区三区| 精品亚洲成a人片在线观看| 日韩熟女老妇一区二区性免费视频| 欧美精品一区二区免费开放| 久久久久人妻精品一区果冻| 在线 av 中文字幕| 午夜福利视频在线观看免费| 不卡视频在线观看欧美| 一级片免费观看大全| 亚洲av.av天堂| 黄色怎么调成土黄色| 少妇 在线观看| 国产淫语在线视频| 99国产综合亚洲精品| 日产精品乱码卡一卡2卡三| 少妇人妻 视频| 高清av免费在线| 一级毛片我不卡| 国产无遮挡羞羞视频在线观看| 女人高潮潮喷娇喘18禁视频| 精品人妻熟女毛片av久久网站| 老鸭窝网址在线观看| 99香蕉大伊视频| 欧美少妇被猛烈插入视频| 亚洲色图综合在线观看| 久久国产亚洲av麻豆专区| 777久久人妻少妇嫩草av网站| 国产黄色视频一区二区在线观看| 久久精品国产亚洲av天美| 久久久久久人人人人人| 高清在线视频一区二区三区| 亚洲精品久久成人aⅴ小说| 国产不卡av网站在线观看| 伦理电影大哥的女人| freevideosex欧美| 制服诱惑二区| 久久久久视频综合| 自拍欧美九色日韩亚洲蝌蚪91| 久久97久久精品| 一二三四在线观看免费中文在| 女人高潮潮喷娇喘18禁视频| 亚洲精品美女久久久久99蜜臀 | 国产1区2区3区精品| 三级国产精品片| 亚洲国产精品999| 中文天堂在线官网| 久久久久国产一级毛片高清牌| 我要看黄色一级片免费的| 乱人伦中国视频| 日韩,欧美,国产一区二区三区| av在线观看视频网站免费| 青春草亚洲视频在线观看| 最近的中文字幕免费完整| 欧美bdsm另类| 香蕉国产在线看| 亚洲av男天堂| av一本久久久久| 毛片一级片免费看久久久久| 黄色视频在线播放观看不卡| 亚洲国产精品国产精品| 十八禁网站网址无遮挡| 大陆偷拍与自拍| 国产精品偷伦视频观看了| 夫妻午夜视频| av不卡在线播放| av在线app专区| 捣出白浆h1v1| 国产极品天堂在线| 日韩三级伦理在线观看| 午夜免费鲁丝| www日本在线高清视频| 久久久久久人人人人人| kizo精华| 精品国产露脸久久av麻豆| 少妇被粗大猛烈的视频| 搡老乐熟女国产| 在线 av 中文字幕| 久久婷婷青草| 久久久久久久久久久久大奶| 国产精品亚洲av一区麻豆 | 国产国语露脸激情在线看| 欧美精品国产亚洲| 日韩欧美精品免费久久| 18禁观看日本| 亚洲欧美成人精品一区二区| 欧美精品一区二区免费开放| 丰满乱子伦码专区| 老司机影院成人| 97在线人人人人妻| 黄色配什么色好看| 久久精品亚洲av国产电影网| 色哟哟·www| 亚洲四区av| 午夜福利视频在线观看免费| 国产精品.久久久| 国产精品成人在线| 大码成人一级视频| 99久国产av精品国产电影| 亚洲国产精品国产精品| 国产在线免费精品| 只有这里有精品99| 免费av中文字幕在线| 免费日韩欧美在线观看| 91成人精品电影| 久久鲁丝午夜福利片| 亚洲精品,欧美精品| 中国国产av一级| 亚洲国产精品成人久久小说| 乱人伦中国视频| 捣出白浆h1v1| 欧美精品一区二区免费开放| 黄片播放在线免费| 在现免费观看毛片| 在线观看人妻少妇| 岛国毛片在线播放| 欧美日韩一级在线毛片| 丰满迷人的少妇在线观看| 曰老女人黄片| 秋霞伦理黄片| 亚洲精品视频女| 男女边吃奶边做爰视频| 女性被躁到高潮视频| 一级爰片在线观看| 高清黄色对白视频在线免费看| 最近最新中文字幕免费大全7| av卡一久久| 国产乱人偷精品视频| 夫妻性生交免费视频一级片| 免费播放大片免费观看视频在线观看| 亚洲精品av麻豆狂野| 999久久久国产精品视频| 美女高潮到喷水免费观看| 免费在线观看黄色视频的| 色播在线永久视频| 黑人巨大精品欧美一区二区蜜桃| 熟女电影av网| 十分钟在线观看高清视频www| 亚洲成人手机| av在线app专区| 国产av精品麻豆| xxx大片免费视频| 亚洲精品日本国产第一区| 免费少妇av软件| 久久久亚洲精品成人影院| 制服人妻中文乱码| 婷婷成人精品国产| 深夜精品福利| 伊人亚洲综合成人网| 91久久精品国产一区二区三区| 国产精品一二三区在线看| 黄片小视频在线播放| 黄片无遮挡物在线观看| 亚洲精品久久午夜乱码| 亚洲国产精品国产精品| 日韩人妻精品一区2区三区| 男女国产视频网站| 丝袜在线中文字幕| 亚洲成色77777| 日日撸夜夜添| 丝袜脚勾引网站| 欧美人与善性xxx| 最新中文字幕久久久久| 啦啦啦视频在线资源免费观看| 亚洲美女搞黄在线观看| 汤姆久久久久久久影院中文字幕| 国产片内射在线| 国产一级毛片在线| 久久鲁丝午夜福利片| 国产一区二区在线观看av| 2021少妇久久久久久久久久久| h视频一区二区三区| 亚洲av国产av综合av卡| 国产综合精华液| 如何舔出高潮| 色网站视频免费| 国产精品免费大片| 老汉色∧v一级毛片| 国产精品一区二区在线观看99| 亚洲精品乱久久久久久| 97人妻天天添夜夜摸| 熟女av电影| av免费观看日本| 亚洲综合精品二区| 少妇被粗大的猛进出69影院| 久久99精品国语久久久| 国产1区2区3区精品| 国产一区有黄有色的免费视频| 伊人亚洲综合成人网| 男人爽女人下面视频在线观看| 热99久久久久精品小说推荐| 如日韩欧美国产精品一区二区三区| 国产成人免费无遮挡视频| 精品一区二区三卡| 18在线观看网站| 中文字幕av电影在线播放| 在线观看人妻少妇| 永久免费av网站大全| 黄片播放在线免费| 亚洲国产精品国产精品| 亚洲欧美精品自产自拍| 美女xxoo啪啪120秒动态图| 久久精品国产鲁丝片午夜精品| 制服丝袜香蕉在线| 中国国产av一级| 久久精品国产鲁丝片午夜精品| 在线看a的网站| 亚洲国产精品一区三区| 久久久久久久大尺度免费视频| 亚洲精品aⅴ在线观看| 亚洲,一卡二卡三卡| 欧美老熟妇乱子伦牲交| 亚洲伊人久久精品综合| 狂野欧美激情性bbbbbb| 亚洲,欧美,日韩| 爱豆传媒免费全集在线观看| 亚洲av成人精品一二三区| 久久热在线av| 精品国产超薄肉色丝袜足j| 欧美最新免费一区二区三区| 日本猛色少妇xxxxx猛交久久| 熟妇人妻不卡中文字幕| 国产成人精品无人区| 啦啦啦视频在线资源免费观看| 国产精品香港三级国产av潘金莲 | 欧美日韩精品网址| 亚洲国产精品999| 色婷婷av一区二区三区视频| 下体分泌物呈黄色| 亚洲成色77777| 国产高清国产精品国产三级| 亚洲伊人色综图| 91午夜精品亚洲一区二区三区| 9191精品国产免费久久| 午夜免费观看性视频| 少妇人妻 视频| 午夜福利影视在线免费观看| 国产精品不卡视频一区二区| 久久这里只有精品19| 午夜日韩欧美国产| 久久亚洲国产成人精品v| 国产极品粉嫩免费观看在线| 男人爽女人下面视频在线观看| 欧美老熟妇乱子伦牲交| 一本色道久久久久久精品综合| 九色亚洲精品在线播放| 宅男免费午夜| 国产精品秋霞免费鲁丝片| 99久久综合免费| 国产精品国产av在线观看| 久久久久精品久久久久真实原创| 国产成人aa在线观看| 精品国产露脸久久av麻豆| 国产一区二区 视频在线| 三上悠亚av全集在线观看| 天天躁夜夜躁狠狠久久av| 最近中文字幕高清免费大全6| 国产成人午夜福利电影在线观看| 亚洲伊人久久精品综合| 亚洲欧美中文字幕日韩二区| 久久亚洲国产成人精品v| 欧美精品国产亚洲| 午夜精品国产一区二区电影| 麻豆av在线久日| 亚洲欧美清纯卡通| 精品少妇黑人巨大在线播放| 啦啦啦在线观看免费高清www| 999久久久国产精品视频| 男女边摸边吃奶| 最近手机中文字幕大全| 99re6热这里在线精品视频| 欧美国产精品一级二级三级| 日日啪夜夜爽| 观看av在线不卡| 一级a爱视频在线免费观看| av.在线天堂| 卡戴珊不雅视频在线播放| 宅男免费午夜| 熟女少妇亚洲综合色aaa.| 男女午夜视频在线观看| 99香蕉大伊视频| 久久人人爽人人片av| 肉色欧美久久久久久久蜜桃| 看免费av毛片| 久久精品熟女亚洲av麻豆精品| 成人黄色视频免费在线看| 下体分泌物呈黄色| 波野结衣二区三区在线| 综合色丁香网| 另类亚洲欧美激情| av网站在线播放免费| 午夜影院在线不卡| 亚洲av男天堂| 亚洲av在线观看美女高潮| 国产一区二区三区av在线| 精品少妇一区二区三区视频日本电影 | 成人18禁高潮啪啪吃奶动态图| 成年动漫av网址| 日韩大片免费观看网站| 老汉色∧v一级毛片| 飞空精品影院首页| 久久99一区二区三区| 亚洲成人手机| 欧美+日韩+精品| 国产片内射在线| 日韩电影二区| 桃花免费在线播放| 激情视频va一区二区三区| 成人毛片a级毛片在线播放| 美国免费a级毛片| 岛国毛片在线播放| 久久99热这里只频精品6学生| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲伊人久久精品综合| 三上悠亚av全集在线观看| 国产 一区精品| 性色av一级| 欧美亚洲 丝袜 人妻 在线| 精品少妇一区二区三区视频日本电影 | 国产成人aa在线观看| 午夜福利网站1000一区二区三区| 国产在视频线精品| 人人妻人人澡人人爽人人夜夜| 免费观看性生交大片5| 国产成人午夜福利电影在线观看| 午夜激情久久久久久久| av线在线观看网站| 欧美日韩精品成人综合77777| 婷婷色av中文字幕| av在线观看视频网站免费| 国产精品免费视频内射| 国产av精品麻豆| 咕卡用的链子| 国产精品.久久久| 久久国产亚洲av麻豆专区| 国产一区二区三区综合在线观看| 亚洲情色 制服丝袜| 蜜桃在线观看..| 欧美中文综合在线视频| 欧美av亚洲av综合av国产av | 欧美老熟妇乱子伦牲交| 免费人妻精品一区二区三区视频| 国产人伦9x9x在线观看 | 午夜福利视频在线观看免费| 国产精品一区二区在线观看99| 少妇人妻久久综合中文| 高清在线视频一区二区三区| 国产一区二区三区av在线| 欧美亚洲日本最大视频资源| 男女啪啪激烈高潮av片| 国产在视频线精品| 欧美最新免费一区二区三区| 精品卡一卡二卡四卡免费| 亚洲人成77777在线视频| 亚洲欧美一区二区三区久久| 亚洲av电影在线进入| 咕卡用的链子| 老司机影院毛片| 午夜日韩欧美国产| 叶爱在线成人免费视频播放| 男男h啪啪无遮挡| 99精国产麻豆久久婷婷| 又黄又粗又硬又大视频| 久久精品国产亚洲av高清一级| www日本在线高清视频| 亚洲一码二码三码区别大吗| 丝袜脚勾引网站| 亚洲av国产av综合av卡| 国产精品国产三级国产专区5o| 在线观看免费日韩欧美大片| 香蕉精品网在线| 少妇人妻 视频| 精品亚洲成国产av| 欧美日韩亚洲国产一区二区在线观看 |