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

    Genomic analysis of isopentenyltransferase genes and functional characterization of TaIPT8 indicates positive effects of cytokinins on drought tolerance in wheat

    2023-01-30 04:47:10NaWangJunChenYuanGaoYonginZhouMingChenZhaoshiXuZhengwuFangYouzhi
    The Crop Journal 2023年1期

    Na Wang,Jun Chen,Yuan Gao,Yongin Zhou,Ming Chen,Zhaoshi Xu,Zhengwu Fang,Youzhi M

    a Hubei Collaborative Innovation Center for Grain Industry/College of Agriculture,Yangtze University,Jingzhou 434000,Hubei,China

    b Institute of Crop Sciences,Chinese Academy of Agricultural Sciences/National Key Facility for Crop Gene Resources and Genetic Improvement,Key Laboratory of Biology and Genetic Improvement of Triticeae Crops,Ministry of Agriculture,Beijing 100081,China

    Keywords:TaIPT Phytohormone Abiotic stress Cytokinin biosynthesis Wheat

    ABSTRACT Cytokinins(CKs)function in plant development and during stress responses,but their role in drought tolerance in wheat is unknown.In the present study,24 isopentenyltransferase(IPT)genes,encoding ratelimiting enzymes in CK biosynthesis were identified in the wheat genome.The chromosomal locations and structures of the genes,protein properties,and phylogenetic relationships were characterized.ATP/ADP TaIPT genes showed tissue-specific expression.TaIPT2,TaIPT7,and TaIPT8 expression was rapidly induced by 0.5–1 h drought treatments,which decreased to low levels after 2 h drought treatment,as did most other TaIPT genes.TaIPT8-5a/5b/5d triple mutants showed decreased levels of tZtype CK under normal and drought conditions and reduced drought tolerance,which,however,did not manifest as phenotype alterations.By contrast,transgenic wheat plants with drought-induced TaIPT8 showed increased drought tolerance.Our study provides a foundation for further investigation of TaIPT genes and novel insights into the role of CKs in the drought response of wheat.

    1.Introduction

    Wheat is a staple crop,particularly in arid and semi-arid areas where drought limits plant growth and productivity[1].Evolving in rapidly changing environments,plants have developed sensitive protection systems to respond and adapt to abiotic stressors.Plant hormones,even at minute concentrations,transmit local or longdistance signals in response to stress[2–4].

    Although abscisic acid(ABA)has long been recognized as the primary stress-induced hormone,regulating biological activities including stomatal closure and induction of antioxidant defense.Increasing evidence[5–7]suggests that cytokinins(CK)too are involved in responses to abiotic stressors.CKs are a class of phytohormones first characterized as inducers of cell division[8,9],and they are synthesized primarily in growing root tips,from where they are transported to shoots[10,11].Along with auxins,CKs function in multiple biotic processes including inhibition of lateral root initiation,leaf expansion,shoot branching,photosynthesis,floral transition,seed germination,and leaf senescence[12].However,the levels of endogenous CKs and the molecular mechanisms of CK effects during stress are not well understood.Generally,production and transport of endogenous CKs decrease in response to extended periods of stress and result in a high ABA:CK ratio,which influences root-to-shoot signaling under abiotic stress[13–18].But CK levels can be rapidly elevated at the initial stage of stress and continue to increase in response to severe stress[19,20].In plants,CK biosynthesis and degradation are catalyzed by respectively phosphate-isopentenyl transferases(IPTs)and CK oxidase/dehydrogenases(CKXs),which are responsible for the fine-tuning of endogenous CK levels[21–23].Lonely Guy(LOG)converts inactive cytokinin nucleotides directly to the active free bases,while conjugation of CKs with glucose mediated by O-glucosyltransferases(ZOG)leads to rapid CK deactivation.The conjugation is reversible and O-glucosides can be converted back to free bases and ribosides by β-glucosidase(GLU),which is necessary for re-establishment or maintenance of CK homeostasis[21,24,25].The expression patterns of individual members of the IPT,CKX,LOG,ZOG,and GLU multigene families are tissue-and development-specific in wheat[25,26].IPTs act as rate-limiting enzymes in the biosynthesis of CKs,and were first identified in Agrobacterium tumefaciens[27,28].In plants,IPTs mediate two CK synthesis pathways:the tRNA and ATP/ADP pathways.The Arabidopsis genome encodes seven ATP/ADP IPT genes(AtIPT1,3,and 4–8)that control the biosynthesis of isopentenyladenine(iP)-and trans-zeatin(tZ)-type CKs and of two tRNA IPTs(AtIPT2 and 9),which are responsible for the synthesis of cis-zeatin(cZ)-type CKs[29–32].Their orthologs have been annotated or functionally verified in several crop species.Ten OsIPTs,ten ZmIPTs,and six TaIPTs have been isolated from rice,maize and wheat,respectively[26].Expression analyses of plant IPT genes under abiotic stress or during ABA treatment showed that the majority of abiotic responsive IPT genes in Arabidopsis were repressed;however,the contrasting pattern was observed in soybean,cabbage,and apple,in which they were up-regulated before returning to pretreatment levels,and a few were maintained at high expression levels under severe stress conditions[16,33–35].Upregulation of CKXs,which reduces CK levels,may be a common mechanism in plants continuously exposed to a stressor[33–37].

    CK and ABA signaling are generally considered antagonistic in regulating seed germination and in response to stressors.In fact,both negative and positive effects of CKs on stress tolerance have been reported in numerous functional studies.In Arabidopsis and tobacco,CK downregulation via overexpression of the CKX gene led to higher drought and salt tolerance accompanied by lower leaf water loss,while under drought conditions,long-term CK elevation by overexpression of IPT or application of exogenous CK resulted in higher water loss,than observed in CKX transformants[16,36,38,39].Transgenic plants expressing an IPT gene through stress-or senescence-induced promoters such as PSARK and PSAG12 were generated[40–45]to examine the effects of increasing CK concentrations on stress tolerance,particularly drought tolerance,and these transgenic plants showed better protection of the photosynthetic system,delayed senescence,and increased redox potential,improving their drought resistance.

    The concentrations of endogenous plant CKs are strongly regulated in response to abiotic stressors,and the roles of IPT genes in regulating CK levels have been identified in Arabidopsis and other plants.However,in wheat,the relationship between CKs and stress tolerance remains unclear.The objective of this study was to perform genome-wide identification of IPT genes in wheat and use phylogenetic and expression analyses to characterize the evolution of the wheat IPT gene family and the functions of its members in drought stress responses.Production of transgenic wheat plants with TaIPT8 homologs induced by CRISPR/Cas9 and droughtinducible expression of PNAC48::TaIPT8 was undertaken to investigate the regulatory roles of CK in wheat drought tolerance.

    2.Materials and methods

    2.1.In silico analysis of IPT genes in wheat

    To identify all latent IPT genes in wheat,IPT protein sequences of Arabidopsis and rice were retrieved from the TAIR(https://www.arabidopsis.org)and RGAP(https://rice.uga.edu/)databases and were used as query sequences to search the wheat genome sequence database in Ensembl Plants(https://plants.ensembl.org/Triticum_aestivum/).A BLASTP search was performed,and all proteins with an E-value<1e-5 were selected.The protein sequences of all putative wheat IPT genes were submitted to the Pfam(https://pfam.janelia.org)and SMART(https://smart.embl-heidelberg.de/)webtools to confirm the presence of a conserved IPPTbinding domain.Sequences lacking the IPPT domain were excluded.In total,24 wheat IPT gene family members were identified.Given that wheat is a hexaploid with three subgenomes(2n=6x=42;AABBDD),we further divided the 24 wheat IPT genes into nine homologous groups,which were designated TaIPTX-YA,TaIPTX-YB,or TaIPTX-YD,where X denotes the homologous group number and Y the wheat chromosome on which it was located.The molecular weights and isoelectric points of the TaIPT proteins were predicted using ExPASy(https://web.expasy.org/compute_pi).

    2.2.Chromosomal locations and phylogenetic analyses

    TBtools[46]was used to visualize the positions of IPT genes on chromosomes based on the physical location of each gene and the length of the respective chromosomes.The location information of each TaIPT gene was obtained from Ensembl Plants(https://plants.ensembl.org/).All TaIPT genes were mapped to 15 wheat chromosomes.For phylogenetic analysis,multiple sequence alignments were produced using ClustalX[47].The neighbor-joining method was used to construct a phylogenetic tree of 50 IPT proteins in wheat,Arabidopsis,rice and barley using MEGA 7.0[48].The confidence limits of the monophyletic groups in the phylogenetic tree were assessed using bootstrap analyses with 1000 replications.

    2.3.Construction of a CRISPR/Cas9-based gene editing and overexpression vector

    To generate TaIPT8-5A/5B/5D (TraesCS5A02G460000,TraesCS5B02G469600,and TraesCS5D02G471100)mutants,two sgRNA target sequences were designed according to the conserved exon sequences of TaIPT8-5A/5B/5D using CRISPR-P 2.0(https://crispr.hzau.edu.cn/CRISPR2/).The segment containing the wheat U3 promoter was amplified from the pCBC-MT1T2 vector using a pair of primers containing the two designed sgRNAs and was then cloned into the BsaI site of the CRISPR/Cas9 vector pBUE411[49]using T4 DNA ligase(TransGen Biotech,Beijing,China).

    To generate overexpressing plants,the coding regions of TaIPT8-5A were amplified from wild-type(WT)cDNA using the DNA polymerase KOD FX(Toyobo,Osaka,Japan)and were cloned into a pMWB110 vector[50]using an In-Fusion Advantage PCR cloning kit(Clontech Laboratories,Takara,Tokyo,Japan).The recombinant plasmids were transformed into the spring wheat cultivar Fielder using Agrobacterium-mediated gene transformation[51].To generate transgenic plants with drought-inducible TaIPT8 overexpression,the ubiquitin promoter in the pMWB110 vector was replaced with a 1.4-Kb TaNAC48 promoter[52]using two restriction enzymes,HindIII and BamHI.All primer sequences are listed in Table S1.

    2.4.Growth conditions and stress treatments

    WT and transgenic plants were planted in pots containing peat soil and were grown in a greenhouse with a 16/8 h light/dark photoperiod,70% relative humidity,and 25/23 °C day/night temperature.Fifteen-day-old WT seedlings were dehydrated by placing on filter paper to induce rapid drought,and the seedlings were then sampled at 0,0.5,1,2,4,8,and 12 h.The samples were immediately frozen at-80 °C.

    To allow WT and transgenic plants to grow under the same soil moisture conditions,seeds of both plant types were sown in the same pot.For drought-tolerance tests,seedlings were grown in the greenhouse under normal conditions for three weeks after which the drought treatment was applied by ceasing irrigation,and when marked differences in wilting were observed between transgenic and WT plants,all plants were rehydrated.All drought treatments were independently repeated three times.

    2.5.RNA extraction and RT-PCR analyses

    Total RNA was extracted from tissues using an RNeasy Plant Mini Kit(Qiagen,Valencia,CA,USA)according to the manufacturer’s instructions,and genomic DNA contamination was removed by digestion with DNase I(Takara).The first strand of cDNA was synthesized using the PrimeScript First Strand cDNA Synthesis Kit(Takara).For semi-quantitative RT-PCR,the wheat Actin gene was used as an internal control.Specific primers for TaIPT genes were designed using Premier 5.0(Premier Biosoft International,Palo Alto,CA,USA),based on cDNA sequences.All primer sequences are listed in Table S1.The PCR conditions included denaturation at 94 °C for 4 min,followed by predetermined number of cycles of denaturation at 94°C for 30 s,annealing at 58°C for 30 s,extension at 72 °C for 30 s,and a final extension at 72 °C for 5 min.The numbers of PCR cycles were chosen to ensure an exponential range of amplification.For the amplification of TaActin sequences,18 to 30 cycles were initially tested,and 24 cycles were used consistently thereafter.For TaIPT genes,the numbers of PCR cycles varied from 29 to 35,according to their expression levels.PCR products were separated by electrophoresis on a 1.2%agarose gel and stained with ethidium bromide.qRT-PCR was performed using gene-specific primers and SYBR Premix Ex Taq reagent(Takara)on an ABI 7500 Fast Real-Time PCR system(Applied Biosystems,Foster City,CA,USA).The TaActin gene was used as an internal control.PCR was performed in triplicate for each sample,using two independent biological replicates.

    2.6.Detection of phytohormones

    Fresh third-leaf samples of WT and transgenic seedlings were harvested 30 days after germination before and after drought treatment(0.5 h);they were immediately frozen in liquid nitrogen,ground to powder,and stored at-80°C until analysis.CK concentrations were detected using MetWare(https://www.metware.cn/)based on the AB Sciex QTRAP4500 LC-MS/MS platform.Three replicates were used per assay.

    2.7.Physiological characteristics

    To measure malondialdehyde(MDA),H2O2,O2-,and proline concentrations,wheat seedlings at the three-leaf stage were grown for approximately-two weeks under drought conditions,and 0.1 g third-leaf blade was sampled before and after drought treatment.Proline and MDA concentrations were measured using a commercial assay kit(Comin,Beijing,China)according to the manufacturer’s instructions.O2-and H2O2concentrations were measured using a superoxide anion content detection kit(BC1295;Solarbio Life Science,Beijing,China)and an H2O2content detection kit(BC3595;Solarbio Life Science,Beijing,China),respectively.Dry and fresh weights of WT and transgenic plants were recorded after one week of rehydration,and water loss was measured as previously described[53].Measurements were performed using three biological replicates,and ANOVA was fitted using SPSS 19.0(IBM Corp.,Armonk,NY,USA).

    3.Results

    3.1.In silico identification of IPT genes in wheat

    To identify wheat IPT genes,a BLAST search using the amino acid sequences of Arabidopsis,rice,and barley IPT genes as queries was performed in the wheat DNA database.Twenty-four nonredundant putative wheat IPT sequences were identified and were designated according to their chromosomal order and genetic relationships(Figs.1A,S1).The 24 TaIPT genes were distributed widely and unevenly on all wheat chromosomes(Fig.S1)except 4A,4B,4D,6A,6B,and 6D.Sequence and phylogenetic analyses revealed that the 24 TaIPTs constituted nine sets,with each set including three homoeologous genes in the A,B,and D subgenomes,except for TaIPT1 and TaIPT2.TaIPT1-1A/1B had no homologs in the D subgenomes,and TaIPT2-1B had no homologs in either the A or D subgenomes(Figs.1A,S1).Each predicted gene was submitted to the Pfam database(https://pfam.janelia.org)to confirm the presence of domain signatures,and all 24 TaIPT proteins were predicted to contain one or two IPP transferase domains(PF01715;Fig.1B;Table 1).Amino acid sequence lengths of the 24 TaIPT proteins ranged from 280 to 466,the isoelectric points ranged from 5.10 to 10.32,and the molecular weights ranged from 30.64 to 52.18 kDa(Table 1).Among these 24 TaIPTs,TaIPT4-2A,TaIPT4-2B,and TaIPT4-2D contained many exons(11,11,and 10,respectively);TaIPT1-1A contained two exons,and the remaining 20 TaIPT members lacked introns(Fig.1C).

    Table 1 IPT genes identified by in silico analysis in wheat.

    Fig.1.Phylogenetic relationships,conserved domains,and gene structures of TaIPTs.(A)The unrooted phylogenetic tree was constructed using MEGA 7.0 by the neighborjoining method with 1000 bootstrap replications.(B)TaIPTs contain one or two IPPT domains.(C)Exon–intron structure of TaIPTs.Blue boxes indicate untranslated 5′-and 3′-regions;yellow boxes indicate exons;black lines indicate introns.

    Fig.2.Phylogenetic relationship of IPTs.The phylogenetic tree was generated using MEGA 7.0 software based on multiple alignments of amino acid sequences from Arabidopsis,rice,barley,and wheat.

    3.2.Phylogenetics of the IPT gene family in Arabidopsis,rice,and wheat

    A phylogenetic analysis of 24 wheat TaIPTs,nine Arabidopsis AtIPTs[29],ten rice OsIPTs[54],and seven barley HvIPTs[55]was performed by generating a phylogenetic tree that classified all IPTs into four major groups.TaIPT4-2A/2B/2D and TaIPT9-7A7/B/7D grouped in groups I and group II with AtIPT2,OsIPT9 and AtIPT9,OsIPT10 respectively,suggesting their involvement in tRNA prenylation.The other 38 IPT genes presumably encoding ATP/ADP IPTs clustered in groups III and IV and were further divided into small subgroups(Fig.2).TaIPT2-1B and TaIPT8-5A/B/D showed the closest relationship with AtIPT3,5,and 7,which affect predominantly the formation of iP-and tZ-types in Arabidopsis.

    3.3.Expression profiles of TaIPTs in tissues of wheat

    Tissue-specific expression of the CK metabolic genes indicated their functional specification,which is useful for the genetic modification of specific traits.Semi-quantitative RT-PCR analysis revealed that TaIPT1 was expressed specifically in roots(Fig.3).By contrast,TaIPT6 and TaIPT9 were expressed mainly in aerial tissues.TaIPT2,TaIPT5,and TaIPT7 were expressed in relatively mature tissues,such as sheaths,stems,and spikes,whereas TaIPT8 was expressed in young roots,shoots,and leaves(Fig.3).Two TaIPTs,TaIPT3,and TaIPT4,were ubiquitously expressed in all examined tissues(Fig.3).

    3.4.Expression analysis of TaIPTs in response to drought stress

    To investigate the response of TaIPTs to the drought stress,the expression patterns of TaIPT genes under drought stress treatment were further studied.Transcription of TaIPT2 and TaIPT8 was induced by 0.5 h drought treatment,while TaIPT7 was highly expressed at 1 and 2 h drought treatments,after which they decreased to a low level.TaIPT3,TaIPT4,and TaIPT9 began to be suppressed at 1 h drought treatment and TaIPT6 at 2 h drought treatment(Fig.4).No significant changes were observed in TaIPT1 or TaIPT5(Fig.4).To further verify drought-induced transcription of TaIPT8,qRT-PCR analyses were performed using TaIPT8 universal primers and TaIPT8-5A,TaIPT8-5B,and TaIPT8-5D specific primers,and stronger upregulation of TaIPT8-5A and TaIPT8-5D was observed than TaIPT8-5B after 0.5 h drought treatment(Fig.S2).

    3.5.Targeted mutagenesis of TaIPT8 induced by CRISPR/Cas9

    To investigate the roles of CK in drought tolerance,CRISPR/Cas9 system was used to specifically disrupt the TaIPT8 gene which was closely related to the crucial ATP/ADP type Arabidopsis IPTs AtIPT3,5,7,and showed rapid responses to drought stress.

    Fig.3.Tissue-specific expression patterns of TaIPT genes.The wheat actin gene was used as the control for equal cDNA concentrations.Roots and shoots were collected at seedling stage.Leaf blade,sheath,stem,and spike tissue was collected at heading stage.

    Fig.5.Homozygous targeted mutagenesis of TaIPT8.(A)Schematic map of gRNA target site selection in the conserved regions of TaIPT8-5A,TaIPT8-5B,and TaIPT8-5D.(B)The sequences and sequencing chromatograms corresponding to the wild type and seven representative editing lines at target sites are shown;‘‘-”indicates deletion of nucleotides;PAM motifs are highlighted in red;Red arrow indicates the start of translation.

    In order to simultaneously knock out TaIPT8 in all subgenomes(A,B,and D),we designed an sgRNA targeting the conserved region near the start codons of TaIPT8-5A,TaIPT8-5B,and TaIPT8-5D(Fig.5A),and a corresponding sgRNA/Cas9 vector was transformed into the wheat cultivar Fielder using Agrobacterium tumefaciens-mediated transformation.After DNA sequencing analysis of T2plants derived from five independent transgenic lines,we obtained four homozygous Taipt8-5a/5b/5d triple mutants designated Taipt8-M1 to Taipt8-M4 and Taipt8-M1;Taipt8-M2 plants were chosen for further functional analysis(Fig.5B).To avoid the influence of potential off-target sites on phenotype and functional analyses,six most likely off-target sites with 2–3 mismatches were selected for amplification and analysis using genome-specific PCR and Sanger sequencing.Sequencing revealed no mutations in the examined potential off-target loci of all homozygous Taipt8 mutants(Table S2).

    3.6.Mutation of TaIPT8 results in reduced drought tolerance

    There were no visible phenotypes in Taipt8-M1 and Taipt8-M2 under normal growth conditions.However,after drought treatment,the Taipt8 mutants showed significantly accelerated leaf rolling,compared to WT plants.After rehydration,the WT plants quickly recovered from water stress,whereas the Taipt8 mutants continuously withered and showed lower fresh and dry weights(Fig.6A–C).Although the water loss rate was comparable between mutant and WT plants(Fig.S3),the proline content of Taipt8 mutants was lower,while the concentrations of MDA,H2O2,and O2-were significantly higher than those of WT plants(Fig.6D–G).

    Fig.6.Knockout of TaIPT8 leads to decreased drought tolerance in wheat.(A)Phenotypes of TaIPT8 mutants and control plants under normal growth and drought conditions.(B,C)Fresh weight(B)and dry weight(C)of TaIPT8 mutants and control plants under drought-stress conditions.(D–G)MDA content(D),proline content(E),O2-content(F),and H2O2 content(G)of TaIPT8 mutants and wild-type plants under normal and drought stress conditions.Shown are means±standard deviation(SD)of three independent biological replications.Differences were evaluated by Student’s t-test(*,P<0.05;**,P<0.01).

    Fig.7.Cytokinin levels in Taipt8 mutant before and after drought treatment.(A)tZ,(B)IP,(C)IPA,and(D)cZ concentrations were measured in three-week-old seedlings before and after 0.5 h drought treatment.Shown are the means±SD of three independent biological replicates.

    Fig.8.Quantitative RT-PCR analysis showing the relative expression levels of TaIPT8 in five drought-induced TaIPT8-overexpressing lines under drought stress.The x-axes shows treatment times of 0,0.5,1,2,4,and 8 h.The wheat actin gene was used as an internal reference.Shown are the means±SD of three independent biological replicates.Different letters in bar graphs indicate differences at P<0.05.

    3.7.Concentrations of CKs in WT and TaIPT8 mutant plants

    As the expression of some TaIPT genes was upregulated after 0.5–1 h drought treatments,we first examined the effects of drought stress on CK metabolism in wheat plants by comparing CK concentrations in WT plants before and after the 0.5 h drought treatment.The tZ content was increased and iP and iPA content was decreased in WT leaves after drought treatment(Fig.7A–C).When the CK content was compared between WT and Taipt8-M1,the bioactive tZ content in the Taipt8 mutants was significantly reduced before and after drought treatment(Fig.7A).Levels of iP and iPA were not affected before treatment,whereas the iP content was decreased in Taipt8 mutants after drought treatment(Fig.7B,C).cZ levels were significantly increased in Taipt8 mutants with and without drought treatment(Fig.7D).

    3.8.Drought-inducible expression of TaIPT8 improves wheat drought tolerance

    To assess the effects of CK overproduction on wheat development and drought tolerance,we overexpressed TaIPT8 in transgenic wheat under the control of a maize ubiquitin promoter.Root formation was completely inhibited in all transgenic lines,which eventually withered without forming seeds(Fig.S4).We induced the ectopic expression of TaIPT8 with the previously identified drought-inducible promoter PNAC48[52],and based on the qRT-PCR results,two independent transgenic lines,DI-OE2 and DI-OE5,which maintained high expression levels of TaIPT8 under drought conditions,were selected from five positive transgenic lines for further analysis(Fig.8).No growth-restraining effect was observed in PNAC48::TaIPT8 transgenic wheat plants grown under well-watered conditions.After drought treatment,both PNAC48::TaIPT8 transgenic lines delayed wilting,and more transgenic plants were recovered with higher fresh weight and dry weight after seven days of rehydration(Fig.9).

    Fig.9.Drought-induced overexpression of TaIPT8 confers improved drought tolerance in wheat.(A)Phenotypes of drought-induced overexpressing lines and wild-type plants under normal growth and drought conditions.(B,C)Fresh weight(B)and dry weight(C)of TaNAC48-overexpressing plants and wild-type plants under drought-stress conditions.Shown are the means±SD of three independent biological replications.Differences were evaluated by Student’s t-test(*,P<0.05;**,P<0.01).

    4.Discussion

    The IPT gene family has been characterized in many plant species because of its important biological functions in catalyzing the rate-limiting step of CK biosynthesis[32,54,56].In addition to the regulation of developmental processes,increasing evidence has shown that CK is essential for stress responses in plants.However,observations of endogenous CK levels under various stress conditions suggested[19]that CK metabolism is differentially regulated in response to abiotic stress in different plant species.To investigate the relationship between CK and drought tolerance in wheat,24 TaIPT genes were identified by in silico analysis.Not all TaIPT genes were represented by three homoeologous genes in the A,B,and D subgenomes,such as TaIPT2-1B,which had no homologs in the A and D subgenomes but is closely related to TaIPT8-5A/5B/5D.We predicted that TaIPT2-1B originated from gene duplication in the distant evolutionary history of the B genome donors.Phylogenetic analysis indicated that TaIPT4 and TaIPT9 would also be involved in the synthesis of cZ-type CKs,as they clustered with two Arabidopsis tRNA-dependent IPTs,whereas TaIPT1,2,3,5,6,7,and 8 may act as ATP/ADP-dependent TaIPTs that are responsible for bioactive CK biosynthesis.Previous study[26]have identified the six wheat IPT genes TaIPT2,3,5,6,7,and 8,which were named by their similarity to ZmIPT2,3,5,6,7,and 8 in maize.We found that TaIPT5 and TaIP6(previously identified)were two homoeologous alleles of TaIPT8(newly named)from 5A and 5B subgenomes and renamed the six previously identified TaIPT as TaIPT1-1D,TaIPT6-3D,TaIPT8-5A,TaIPT8-5B,TaIPT5-3D,and TaIPT3-1D according to the results of amino acid sequence alignment.

    Although CKs have long been considered to be synthesized primarily in roots and transported acropetally to the shoots,recent study[57]confirmed that CKs synthesized locally rather than in roots are even more critical for axillary bud growth.The expression patterns of individual TaIPT genes,especially ATP/ADP-dependent TaIPTs,are differential.TaIPT1 was expressed specifically in roots,and TaIPT8 was expressed mainly in shoots,where it may be involved in the regulation of vegetative growth,whereas TaIPT2,TaIPT5,and TaIPT7 were expressed predominantly in sheath,stem,and spike tissue,where they may be involved in the regulation of reproductive growth.Similarly,considering their specificity of expression,15 wheat TaCKX genes were assigned to four groups specific to leaves,developing spikes,roots,and all the organs tested respectively.These assignments indicated that CK content is subtly controlled in diverse organs and development stages[23].The expression of IPTs has been shown[33–35]to be a response to abiotic and biotic stressors in many plant species.In maize,salt and osmotic stresses induce the expression of most CK biosynthetic genes in seedlings of maize,leading to a moderate increase of CKs lasting several days during acclimatization to stress[58].In the present study,expression of most TaIPTs was gradually reduced after 1 h drought treatment;however,TaIPT2,TaIPT7,and TaIPT8 were up-regulated within 0.5–1 h,which may have induced the elevation of tZ but not of iP or iPA in WT plants at the early stages of drought treatment.In TaIPT8-5a/5b/5d triple mutants,levels of tZ-type CK rather than iP-type CK were decreased under normal and drought conditions,indicating that TaIPT8 homologs are more critical for the formation of tZ-type CK in wheat.

    In Arabidopsis,neither single-nor double-atipt mutants show visible changes in their phenotypes due to gene redundancy[29].Similarly,the decrease in CK content did not change the morphology of the TaIPT8-5a/5b/5d triple mutants in the present study;however,it reduced the drought tolerance of mutant plants.In contrast,downregulation of CK content via mutation of multiple IPT or overexpression of CKX genes has been reported[16,38]to increase drought tolerance in Arabidopsis and tobacco.These multi-gene mutants or overexpression plants show altered phenotypes,including slow growth rates,dwarf shoots,and small leaves,but enlarged root systems,which may contribute to lower transpiration rates and higher relative water content[29,38,39].The CKdeficient phenotype with enlarged root system and suppressed aerial part could also be observed in monocots.Constitutive overexpression of CKX in barley plants led to a substantially enlarged root system,but shoot growth was retarded and the plants died without reaching flowering[55].It thus remains uncertain whether the decreased CK levels achieved by overproduction of the CKX or mutation of IPT has a same effect on abiotic stress tolerance in other cereal crops.When the CK effect on drought tolerance in wheat was evaluated,the water loss rate was comparable between TaIPT8 mutant and WT plants,while the contents of MDA,H2O2,and O2-in transgenic plants were higher than those in WT plants under drought,owing to a decrease in active oxygen-scavenging ability of their antioxidant systems.Thus,although these TaIPTs may function redundantly during plant organ development,their spatially specific expression pattern probably influences plant stress response.

    Given that employment of a constitutive promoter for IPT overexpression resulted in diminished root growth,genetic manipulation of crop plants using inducible promoters to increase stress tolerance is inevitable.Conditionally activated expression of IPT under stress-or senescence-induced promoters stimulated extreme drought tolerance without impairing plant morphology in cotton[40],maize[41],peanut[42],rice[44],and cassava[45].Most transgenic plants showed improved drought resistance by increased protection of the photosynthetic system,delayed senescence,and increased redox ability.A similar strategy was used in the present study,where TaIPT8 expression was induced by the wheat native drought-induced-promoter PNAC48.TaIPT8 maintained a high expression level in the resulting transgenic lines,increasing drought tolerance without change in plant development.Our findings of water loss rates in wheat Taipt8 mutants and WT plants also suggested that CKs exerted their beneficial effects on drought tolerance via regulation of cellular redox status.

    In summary,we performed an in silico analysis of 24 TaIPTs in the wheat genome and characterized their chromosomal locations,phylogenetic relationships,gene structures,protein properties,and gene expression patterns.TaIPT8 mutants and drought-induced overexpression plants were produced,and functional analyses revealed beneficial effects of CKs on wheat drought tolerance.Further research may reveal the mechanisms and signaling pathways involved in the interactions of CK metabolism and stress tolerance in wheat.

    CRediT authorship contribution statement

    Na Wang:Data curation,Formal analysis,Visualization,Writing-review&editing.Jun Chen:Data curation,Formal analysis,Visualization,Writing-original draft.Yuan Gao:Investigation,Methodology,Validation.Yongbin Zhou:Resources,Software.Ming Chen:Resources,Software.Zhaoshi Xu:Supervision,Validation.Zhengwu Fang:Data curation,Project administration,Resources,Writing-review & editing.Youzhi Ma:Conceptualization,Data curation,Funding acquisition,Project administration,Resources,Writing-review & 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 thank Dr.Genying Li from the Shandong Academy of Agricultural Sciences for help with the wheat genetic transformation.This work was financially supported by the National Key Research and Development Program of China(2020YFE0202300)and the Agricultural Science and Technology Innovation Program(CAASZDRW202109 and CAAS-ZDRW202002).

    Appendix A.Supplementary data

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

    黄色配什么色好看| 人人妻人人爽人人添夜夜欢视频 | 国产探花在线观看一区二区| 又黄又爽又刺激的免费视频.| 一级爰片在线观看| 80岁老熟妇乱子伦牲交| 亚洲美女搞黄在线观看| 久久午夜福利片| 一区二区三区乱码不卡18| 亚洲高清免费不卡视频| 亚洲最大成人手机在线| 国产精品精品国产色婷婷| 国产精品爽爽va在线观看网站| 国产男人的电影天堂91| 国产国拍精品亚洲av在线观看| 精品久久久久久电影网| eeuss影院久久| 两个人的视频大全免费| 日本与韩国留学比较| 亚洲精品中文字幕在线视频 | 久久久久久久大尺度免费视频| 欧美成人精品欧美一级黄| 特大巨黑吊av在线直播| 高清毛片免费看| 欧美成人午夜免费资源| 人人妻人人看人人澡| 亚洲在久久综合| 国产色爽女视频免费观看| 简卡轻食公司| 国产高清有码在线观看视频| 免费电影在线观看免费观看| 高清av免费在线| 精华霜和精华液先用哪个| 人体艺术视频欧美日本| 日韩一区二区视频免费看| 国产69精品久久久久777片| 欧美成人午夜免费资源| 久久久久精品性色| 嫩草影院新地址| 免费电影在线观看免费观看| 亚洲精品亚洲一区二区| 又大又黄又爽视频免费| 久久久久久国产a免费观看| 视频中文字幕在线观看| 成人国产av品久久久| 久久久久久久精品精品| 菩萨蛮人人尽说江南好唐韦庄| 麻豆久久精品国产亚洲av| 免费av不卡在线播放| 亚洲精华国产精华液的使用体验| 亚洲精品国产av蜜桃| 中文乱码字字幕精品一区二区三区| 我要看日韩黄色一级片| 日本爱情动作片www.在线观看| 国产欧美另类精品又又久久亚洲欧美| 亚洲精品日韩在线中文字幕| 精品酒店卫生间| 日韩视频在线欧美| 日日摸夜夜添夜夜添av毛片| 久久精品国产鲁丝片午夜精品| 国产美女午夜福利| 久久韩国三级中文字幕| 一边亲一边摸免费视频| 日韩成人av中文字幕在线观看| 18禁裸乳无遮挡动漫免费视频 | 好男人在线观看高清免费视频| 综合色丁香网| 一本色道久久久久久精品综合| av又黄又爽大尺度在线免费看| 街头女战士在线观看网站| 欧美最新免费一区二区三区| 亚洲av中文字字幕乱码综合| 亚洲不卡免费看| 欧美极品一区二区三区四区| 丝袜脚勾引网站| 欧美xxxx黑人xx丫x性爽| 久久ye,这里只有精品| 免费不卡的大黄色大毛片视频在线观看| 亚洲国产精品专区欧美| 亚洲精品,欧美精品| 九九在线视频观看精品| 国产伦在线观看视频一区| 亚洲人成网站高清观看| 国产一区二区三区综合在线观看 | 一级a做视频免费观看| 日韩精品有码人妻一区| 26uuu在线亚洲综合色| 蜜臀久久99精品久久宅男| 你懂的网址亚洲精品在线观看| 内射极品少妇av片p| 国产色婷婷99| 成人美女网站在线观看视频| 亚洲精品乱码久久久v下载方式| 在线免费观看不下载黄p国产| 一区二区三区乱码不卡18| 校园人妻丝袜中文字幕| 综合色丁香网| 嫩草影院精品99| 亚洲精品色激情综合| 26uuu在线亚洲综合色| 婷婷色麻豆天堂久久| av在线蜜桃| 亚洲精品亚洲一区二区| 免费观看在线日韩| 日韩,欧美,国产一区二区三区| 午夜视频国产福利| 国产 一区 欧美 日韩| 日韩成人伦理影院| 在线观看一区二区三区激情| 中文字幕免费在线视频6| 99热这里只有是精品在线观看| 在线 av 中文字幕| 日本三级黄在线观看| 九九爱精品视频在线观看| 欧美丝袜亚洲另类| 水蜜桃什么品种好| 97热精品久久久久久| 国产成人免费无遮挡视频| 色视频www国产| 亚洲精品第二区| 在线免费十八禁| 插阴视频在线观看视频| 日本av手机在线免费观看| 狂野欧美激情性xxxx在线观看| 最后的刺客免费高清国语| 91在线精品国自产拍蜜月| 国产国拍精品亚洲av在线观看| 免费电影在线观看免费观看| 99热国产这里只有精品6| 97在线视频观看| 在线观看三级黄色| 国产精品国产三级国产av玫瑰| 18禁在线无遮挡免费观看视频| 美女xxoo啪啪120秒动态图| 一级a做视频免费观看| 性色avwww在线观看| 亚洲,一卡二卡三卡| 青春草亚洲视频在线观看| 久久人人爽人人爽人人片va| 亚洲综合色惰| 久久综合国产亚洲精品| 在线观看美女被高潮喷水网站| 日本猛色少妇xxxxx猛交久久| 亚洲精品乱久久久久久| 中国三级夫妇交换| 九九久久精品国产亚洲av麻豆| 免费观看性生交大片5| 噜噜噜噜噜久久久久久91| 国产男女内射视频| 大片电影免费在线观看免费| 日韩大片免费观看网站| 最近最新中文字幕免费大全7| 视频中文字幕在线观看| 久久久精品欧美日韩精品| 男女边吃奶边做爰视频| 色哟哟·www| 久久久久久国产a免费观看| 卡戴珊不雅视频在线播放| 好男人视频免费观看在线| 亚洲国产精品成人综合色| 日本av手机在线免费观看| 久久久久久久精品精品| 黄色怎么调成土黄色| 国产中年淑女户外野战色| 午夜激情福利司机影院| 亚洲欧美一区二区三区国产| 国产美女午夜福利| 国产黄色视频一区二区在线观看| 简卡轻食公司| 人妻 亚洲 视频| 99久久精品一区二区三区| 99热这里只有是精品50| 国产日韩欧美在线精品| 国产探花在线观看一区二区| 欧美一级a爱片免费观看看| 伦精品一区二区三区| 伦精品一区二区三区| 日韩国内少妇激情av| 日韩,欧美,国产一区二区三区| 久久久久久久久久成人| 1000部很黄的大片| 国产女主播在线喷水免费视频网站| av.在线天堂| 欧美xxⅹ黑人| 免费观看在线日韩| 午夜福利网站1000一区二区三区| 国产男人的电影天堂91| 欧美国产精品一级二级三级 | 欧美一级a爱片免费观看看| 亚洲人成网站高清观看| 国产又色又爽无遮挡免| 制服丝袜香蕉在线| 国产午夜福利久久久久久| 自拍偷自拍亚洲精品老妇| 男女边吃奶边做爰视频| 国产精品三级大全| 国产高清有码在线观看视频| 好男人在线观看高清免费视频| 免费看av在线观看网站| 免费黄频网站在线观看国产| 亚洲av中文字字幕乱码综合| 五月开心婷婷网| 亚洲激情五月婷婷啪啪| 精品国产三级普通话版| 欧美国产精品一级二级三级 | 99视频精品全部免费 在线| 亚洲精品日韩在线中文字幕| 亚洲av成人精品一二三区| 国产成人福利小说| 卡戴珊不雅视频在线播放| 国产 精品1| 别揉我奶头 嗯啊视频| 久久久欧美国产精品| 在线a可以看的网站| 欧美亚洲 丝袜 人妻 在线| 美女国产视频在线观看| 久久ye,这里只有精品| 欧美日韩国产mv在线观看视频 | 久久人人爽人人片av| 久久精品国产亚洲av天美| 另类亚洲欧美激情| 大陆偷拍与自拍| 国产真实伦视频高清在线观看| 日韩三级伦理在线观看| 日本爱情动作片www.在线观看| 成年人午夜在线观看视频| 蜜桃久久精品国产亚洲av| 赤兔流量卡办理| 国产老妇女一区| 亚洲成人精品中文字幕电影| 国产欧美另类精品又又久久亚洲欧美| 欧美另类一区| 王馨瑶露胸无遮挡在线观看| 亚洲久久久久久中文字幕| 另类亚洲欧美激情| 日韩欧美 国产精品| 国产白丝娇喘喷水9色精品| 肉色欧美久久久久久久蜜桃 | 99热这里只有是精品在线观看| 国产伦在线观看视频一区| 亚洲伊人久久精品综合| 一区二区三区乱码不卡18| 女的被弄到高潮叫床怎么办| 久久久久久久久久成人| 高清午夜精品一区二区三区| 国产亚洲5aaaaa淫片| 精品国产三级普通话版| 在线观看国产h片| 久久久久网色| 99热这里只有精品一区| 小蜜桃在线观看免费完整版高清| 久久女婷五月综合色啪小说 | 97超视频在线观看视频| 一级二级三级毛片免费看| 精品久久久精品久久久| 干丝袜人妻中文字幕| 国产精品麻豆人妻色哟哟久久| 久久鲁丝午夜福利片| 汤姆久久久久久久影院中文字幕| 亚洲国产精品专区欧美| 午夜亚洲福利在线播放| 国产精品99久久99久久久不卡 | 制服诱惑二区| 在线观看三级黄色| 亚洲 欧美一区二区三区| 下体分泌物呈黄色| 亚洲欧美色中文字幕在线| 亚洲人成电影观看| 黄色毛片三级朝国网站| 男女午夜视频在线观看| 电影成人av| 欧美另类一区| 一级黄片播放器| 国产在线视频一区二区| 日日撸夜夜添| 日韩一区二区三区影片| 性高湖久久久久久久久免费观看| 少妇的丰满在线观看| 午夜福利影视在线免费观看| 欧美成人精品欧美一级黄| 国产一区二区激情短视频 | 一级黄片播放器| tube8黄色片| 看十八女毛片水多多多| 街头女战士在线观看网站| 最近中文字幕高清免费大全6| 丝瓜视频免费看黄片| av有码第一页| 最近最新中文字幕免费大全7| 高清欧美精品videossex| 亚洲国产成人一精品久久久| 国产不卡av网站在线观看| 蜜桃在线观看..| 亚洲av男天堂| 亚洲欧美精品综合一区二区三区| 在线观看免费视频网站a站| 亚洲国产欧美一区二区综合| 国产一区二区 视频在线| 在线观看人妻少妇| 最新的欧美精品一区二区| 99久久综合免费| 91成人精品电影| 永久免费av网站大全| 国产97色在线日韩免费| 少妇猛男粗大的猛烈进出视频| 日韩中文字幕欧美一区二区 | 国产一区有黄有色的免费视频| 51午夜福利影视在线观看| 免费av中文字幕在线| av女优亚洲男人天堂| 日韩精品免费视频一区二区三区| 99国产综合亚洲精品| 日韩不卡一区二区三区视频在线| 丝袜美足系列| 国产精品久久久av美女十八| 在线看a的网站| 亚洲情色 制服丝袜| 精品亚洲成a人片在线观看| 最近最新中文字幕大全免费视频 | 日本午夜av视频| 一级毛片电影观看| 国产日韩欧美视频二区| 成年人免费黄色播放视频| 亚洲欧洲国产日韩| 久久精品国产综合久久久| 久久久久久久国产电影| 国产在线免费精品| 亚洲欧洲国产日韩| 天天躁日日躁夜夜躁夜夜| 亚洲第一青青草原| 国产一区二区在线观看av| 精品久久久精品久久久| 一个人免费看片子| 国产黄频视频在线观看| 91精品国产国语对白视频| 91成人精品电影| 国产一区二区激情短视频 | 久久精品亚洲熟妇少妇任你| 纯流量卡能插随身wifi吗| av视频免费观看在线观看| 成人亚洲欧美一区二区av| 青草久久国产| 欧美中文综合在线视频| 亚洲欧美中文字幕日韩二区| 亚洲在久久综合| 中文精品一卡2卡3卡4更新| 一级毛片我不卡| 纯流量卡能插随身wifi吗| 男女床上黄色一级片免费看| 日韩人妻精品一区2区三区| 国产一区亚洲一区在线观看| 大片免费播放器 马上看| 国产精品蜜桃在线观看| 欧美黄色片欧美黄色片| 国产xxxxx性猛交| 亚洲在久久综合| 久久性视频一级片| 免费久久久久久久精品成人欧美视频| 亚洲自偷自拍图片 自拍| 男女午夜视频在线观看| av有码第一页| 哪个播放器可以免费观看大片| 夫妻午夜视频| 老司机亚洲免费影院| 国产乱人偷精品视频| 国产在线免费精品| 香蕉国产在线看| 国产精品一区二区在线不卡| 国产精品久久久久久精品电影小说| 国产乱人偷精品视频| 女性生殖器流出的白浆| 波多野结衣一区麻豆| 天美传媒精品一区二区| 欧美成人午夜精品| 亚洲人成网站在线观看播放| 满18在线观看网站| 成人国产av品久久久| 一级毛片 在线播放| 久久精品国产综合久久久| 日韩制服丝袜自拍偷拍| 欧美激情极品国产一区二区三区| 青草久久国产| 午夜免费鲁丝| 麻豆乱淫一区二区| 秋霞伦理黄片| 国产成人精品久久久久久| 亚洲成国产人片在线观看| 精品国产国语对白av| 男人舔女人的私密视频| 男女床上黄色一级片免费看| 制服诱惑二区| 最近中文字幕高清免费大全6| 欧美成人精品欧美一级黄| 国产欧美亚洲国产| 美女主播在线视频| av国产精品久久久久影院| 18禁动态无遮挡网站| 九草在线视频观看| 看非洲黑人一级黄片| 高清不卡的av网站| 亚洲伊人色综图| 欧美激情高清一区二区三区 | 国产一级毛片在线| 亚洲精品美女久久久久99蜜臀 | 1024视频免费在线观看| 男女高潮啪啪啪动态图| 九色亚洲精品在线播放| 午夜福利一区二区在线看| 亚洲精品日韩在线中文字幕| 中文字幕最新亚洲高清| 日韩中文字幕欧美一区二区 | 亚洲精品久久成人aⅴ小说| 男女边摸边吃奶| 久久人人爽av亚洲精品天堂| 老司机靠b影院| 90打野战视频偷拍视频| 最近最新中文字幕免费大全7| 青春草国产在线视频| 国产成人精品在线电影| 精品久久蜜臀av无| 精品国产超薄肉色丝袜足j| 国产一区二区激情短视频 | 高清在线视频一区二区三区| 日韩人妻精品一区2区三区| 老司机影院成人| 亚洲欧美精品自产自拍| 一区二区三区乱码不卡18| 亚洲男人天堂网一区| 夫妻午夜视频| 黄色 视频免费看| 免费女性裸体啪啪无遮挡网站| 丁香六月欧美| 伊人久久国产一区二区| 国产成人精品久久久久久| 亚洲欧洲日产国产| 亚洲熟女毛片儿| 黄色视频不卡| 日本vs欧美在线观看视频| 高清欧美精品videossex| 成人午夜精彩视频在线观看| 黄片无遮挡物在线观看| 中文字幕高清在线视频| 巨乳人妻的诱惑在线观看| 视频在线观看一区二区三区| 最近中文字幕高清免费大全6| 精品人妻一区二区三区麻豆| 久久精品国产亚洲av高清一级| 国产亚洲最大av| 国产成人欧美| 成人亚洲欧美一区二区av| 国产亚洲欧美精品永久| 99热全是精品| 19禁男女啪啪无遮挡网站| 亚洲伊人久久精品综合| 国产精品久久久久成人av| 久热这里只有精品99| 这个男人来自地球电影免费观看 | 亚洲人成77777在线视频| av不卡在线播放| 亚洲,欧美精品.| 国产熟女午夜一区二区三区| 毛片一级片免费看久久久久| 亚洲成人免费av在线播放| 一级片免费观看大全| av电影中文网址| 夫妻性生交免费视频一级片| 我的亚洲天堂| 日韩视频在线欧美| 中国三级夫妇交换| 婷婷成人精品国产| 亚洲av欧美aⅴ国产| 美女福利国产在线| 18禁国产床啪视频网站| 亚洲精品国产色婷婷电影| 欧美变态另类bdsm刘玥| 午夜激情久久久久久久| 天天躁狠狠躁夜夜躁狠狠躁| 国产成人免费观看mmmm| 丁香六月天网| av免费观看日本| 婷婷成人精品国产| 久久久亚洲精品成人影院| 欧美人与性动交α欧美软件| 夫妻性生交免费视频一级片| 国产1区2区3区精品| 中文字幕人妻丝袜一区二区 | 不卡av一区二区三区| 午夜福利免费观看在线| 熟女少妇亚洲综合色aaa.| 热re99久久精品国产66热6| 国产有黄有色有爽视频| 夫妻午夜视频| 日韩视频在线欧美| 国产精品成人在线| 97精品久久久久久久久久精品| 国产熟女欧美一区二区| 啦啦啦在线免费观看视频4| 国产毛片在线视频| 成人漫画全彩无遮挡| 成年av动漫网址| 如日韩欧美国产精品一区二区三区| 777久久人妻少妇嫩草av网站| 国产精品.久久久| 人人澡人人妻人| 国产成人一区二区在线| 国产视频首页在线观看| 国产又色又爽无遮挡免| 国产一区亚洲一区在线观看| 少妇精品久久久久久久| xxx大片免费视频| 国产深夜福利视频在线观看| 欧美日韩视频精品一区| 久久人人爽人人片av| 免费高清在线观看视频在线观看| 免费黄色在线免费观看| 午夜日本视频在线| 一边摸一边做爽爽视频免费| 婷婷色麻豆天堂久久| 一级黄片播放器| 十八禁人妻一区二区| 建设人人有责人人尽责人人享有的| 母亲3免费完整高清在线观看| a 毛片基地| 视频在线观看一区二区三区| 久久精品久久久久久久性| 欧美人与性动交α欧美精品济南到| 国产亚洲av高清不卡| 成人黄色视频免费在线看| 久久人人97超碰香蕉20202| 日韩一卡2卡3卡4卡2021年| 纯流量卡能插随身wifi吗| 日本一区二区免费在线视频| 日日啪夜夜爽| 国产成人免费观看mmmm| 伊人久久国产一区二区| 色吧在线观看| 国产xxxxx性猛交| 一级毛片电影观看| 熟女av电影| 啦啦啦中文免费视频观看日本| 19禁男女啪啪无遮挡网站| 大片免费播放器 马上看| 久久女婷五月综合色啪小说| av电影中文网址| 亚洲第一区二区三区不卡| 久久韩国三级中文字幕| 午夜福利在线免费观看网站| 啦啦啦在线免费观看视频4| 亚洲成人av在线免费| 超碰成人久久| 男女下面插进去视频免费观看| 女性被躁到高潮视频| 成人免费观看视频高清| 国产免费一区二区三区四区乱码| 丝袜美腿诱惑在线| 色婷婷久久久亚洲欧美| 天天操日日干夜夜撸| 亚洲三区欧美一区| 中文欧美无线码| 伊人久久大香线蕉亚洲五| 七月丁香在线播放| 国产野战对白在线观看| av片东京热男人的天堂| 国产淫语在线视频| 午夜久久久在线观看| 1024香蕉在线观看| 国产成人av激情在线播放| 伦理电影大哥的女人| 亚洲成人手机| 久久精品国产亚洲av高清一级| a 毛片基地| 91国产中文字幕| 在线观看www视频免费| 69精品国产乱码久久久| 飞空精品影院首页| 日本色播在线视频| 国产无遮挡羞羞视频在线观看| 99香蕉大伊视频| 亚洲图色成人| 狂野欧美激情性xxxx| 精品视频人人做人人爽| 免费av中文字幕在线| 高清视频免费观看一区二区| 天美传媒精品一区二区| 国产免费现黄频在线看| 不卡av一区二区三区| av网站免费在线观看视频| 亚洲熟女毛片儿| av免费观看日本| 精品人妻一区二区三区麻豆| 自拍欧美九色日韩亚洲蝌蚪91| 日韩欧美一区视频在线观看| 激情五月婷婷亚洲| 亚洲第一青青草原| 日本爱情动作片www.在线观看| 免费高清在线观看日韩| a级毛片黄视频| 午夜av观看不卡| 18禁观看日本| 免费不卡黄色视频| tube8黄色片| 国产男人的电影天堂91| 欧美黑人精品巨大| 国产亚洲午夜精品一区二区久久| 赤兔流量卡办理| 国产人伦9x9x在线观看| 91精品三级在线观看| 男男h啪啪无遮挡| 久久女婷五月综合色啪小说| 午夜福利影视在线免费观看| av网站免费在线观看视频| 亚洲精品第二区| 亚洲天堂av无毛| 91成人精品电影| 久久热在线av|