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

    ldentification and expression of the CEP gene family in apple(Malus×domestica)

    2018-02-05 07:10:46LlRuiANJianpingYOUChunxiangSHUJingWANGXiaofeiHAOYujin
    Journal of Integrative Agriculture 2018年2期

    Ll Rui, AN Jian-ping, YOU Chun-xiang, SHU Jing, WANG Xiao-fei, HAO Yu-jin

    1 National Key Laboratory of Crop Biology, National Research Center for Apple Engineering and Technology/College of Horticulture Science and Engineering, Shandong Agricultural University, Tai’an 271018, P.R.China

    2 Shandong Agriculture and Engineering University, Jinan 251100, P.R.China

    1. lntroduction

    The majority of signaling peptides are small cleavage products of precursor peptides. Several of these precursors must be post-translationally modified to form mature peptides, which usually contain 20 amino acids (Butenkoet al.2009; Murphyet al. 2012). In theArabidopsisgenome, over 1 000 putative small signaling peptides have been predicted,and some of them have been functionally identified, including C-TERMINALLY ENCODED PEPTIDE1 (CEP1), CLAVATA3(CLV3), CLV3/EMBRYO SURROUNDING REGIONRELATED (CLE), and RAPID ALKALINIZATION FACTOR (RALF)(Czyzewiczet al. 2013). Small and secreted regulatory peptides are a growing class of signaling molecules that are involved in regulating plant developmental programs and adapting to extreme environmentviacell-to-cell communica-tion (Katsiret al. 2011; Murphyet al. 2012). Several families of regulatory peptides have been functionally identified and have been found to participate in this process (Butenkoet al. 2009;Czyzewiczet al. 2013; Leeet al. 2015; Songet al. 2016).

    Among the peptide families, the C-TERMINALLYENCODED PEPTIDE (CEP) genes contain a conserved 15-amino acid peptide domain at or near the C-terminus and are characterized inArabidopsis(Ohyamaet al. 2008). The post-translationally modified CEP family members contain an N-terminal secretion signal (NSS) and one or more conserved CEP domains (Mohd-Radzmanet al. 2015).CEPgenes are widely present among gymnosperm and angiosperm plants, but absent in land plants that lack true branching roots or root vasculature, indicating that their emergence coincides with the evolution of seed plants(Delayet al. 2013; Ogilvieet al. 2014). Members of theCEPfamily have already shown to regulate plant lateral root and root nodule development as well as root/shoot growth(Delayet al. 2013; Robertset al. 2013; Mohd-Radzman and Laffont 2016). For example, overexpression ofAtCEP1or treatment with chemically synthesized CEP1 peptide inArabidopsisresults in a reduction in the number of emerged lateral roots and the inhibition of primary root growth (Ohyamaet al. 2008). InMedicago, overexpression ofMtCEP1increases nodulation by promoting rhizobial infections and exhibits repression in lateral root development (Iminet al.2013; Mohd-Radzman and Laffont 2016). In addition,CEPgenes are reported to be negative regulators that mediate environmental influences on plant development (Delayet al. 2013). TheAtCEP3loss-of-function mutant enhances root development under adverse environmental conditions(Delayet al. 2013).

    Currently, an enormous number of possible peptide ligand-receptors (kinases) has been identified and several receptor-like proteins such as CLV1/2, FER, and CEPR1/2 have been found to be involved with root development through interaction with peptides (Kondoet al. 2011; Duet al. 2016; Robertset al. 2016). The leucine-rich repeat(LRR) receptor kinases CEP RECEPTOR 1 (CEPR1;At5g49660) and CEP RECEPTOR 2 (CEPR2; At1g72180)have been shown to be the receptor for CEP1 and other CEPs (Bryanet al. 2012; Tabataet al. 2014; Robertset al.2016). Further studies may reveal that CEPs and CEPR1 participate in the N-dependent responses in long-distance systemic signaling pathways (Okamotoet al. 2016).

    Apple (Malus×domestica) is one of the most widely cultivated fruit crops worldwide and is the most economically important woody plant in temperate regions (Dimick and Hoskin 1983; Leeet al. 2007). The draft genome sequence of apple has been completed, which allows genome-wide analyses of specific gene families (Velascoet al. 2010). Genome-wide analyses of theCEPgenes have been reported inArabidopsis thaliana(Delayet al. 2013) andOryza sativa(Suiet al. 2016), and several gene families have been identified in apple. However, there is no genome wide data regarding the appleCEPgenes. In brief, small signaling peptides play an essential role in all stages of plant growth and development.The characterization of apple CEP peptides provides insight into the molecular mechanism of apple root growth and the responses to different environmental factors.

    In this study, a genome-wide analysis of theCEPgene family was conducted using the apple genome database,and the chromosome locations and gene structures of the putativeCEPgenes were analyzed. Next, the motif compositions of theMdCEPswere obtained using the MEME Program. Subsequently, the expression patterns ofMdCEPgenes in different tissues and in response to abiotic stresses were analyzed. This study provides a foundation for future research into the functional roles ofMdCEPgenes.

    2. Materials and methods

    2.1. ldentification and annotation MdCEP genes in apple

    To identify members of theCEPgene family, all knownArabidopsisCEP protein sequences were retrieved from the database of Institute for Genomic Research (TIGR) and used as queries in BLASTP searches against the Genome Database for Rosaceae (GDR) (http://www.rosaceae.org/).Stand-alone versions of BLASP (http://blast.ncbi.nlm.nih.gov), which are available from NCBI, were used with the e-value cutoff of 1e-003. Then, the predictedCEPgene family sequences were downloaded from the GDR database. All of the protein sequences that were derived from the selectedMdCEPcandidate genes were examined with the domain analysis programs Pfam (http://pfam.sanger.ac.uk/) and Simple Modular Architecture Research Tool(SMART; http://smart.embl-heidelberg.de/), with the default cutoff parameters.

    2.2. Chromosomal locations and gene structures of MdCEP genes

    The chromosomal locations and gene structures were retrieved from the apple genome data that were downloaded from the GDR database. The chromosomal map showing the physical location of all of theMdCEPgenes was generated with the MapDraw Software and the gene structures of theMdCEPgenes were generated with GSDS (http://gsds.cbi.pku.edu.cn/). The isoelectric point (pI) and molecular weight of MdCEPs were obtained with the assistance of proteomics and sequence analysis tools on the ExPASy Proteomics Server (http://expasy.org/). All putative MdCEPs were analyzed by Multiple Em for Motif Elicitation (MEME, ver.4.11.2), a motif search tool to identify conserved motifs shared among CEP proteins with the following parameters, optimum motif width ranges from 6 to 20 and a maximum number of motifs of 10. All other parameters were set at default.

    2.3. Sequence alignment and phylogenetic analysis

    Multiple sequence alignments of the MdCEP sequences(MdCEP1-27) and AtCEP sequences (AtCEP1-15) were performed using Clustal X (ver. 2.1). The MUSCLE Program(ver. 3.52) was also used to perform multiple sequence alignments to confirm the Clustal X result (Edgar 2004).Phylogenetic trees of the MdCEPs were constructed based on the protein sequences according to the neighbor-joining(NJ) method of the MEGA5 Program (ver. 5.1, http://www.megasoftware.net/) (Wanget al. 2012; Hall 2013), with p-distance and the complete deletion option parameters engaged. The reliability of the trees was tested using bootstrapping with 1 000 replicates. The images of the phylogenetic treeos were drawn using the MEGA5 Program.

    2.4. Plant materials, treatments and gene expression analysis

    To investigate the expression ofMdCEPgenes under different apple tissues, the samples (roots, stems, leaves, flowers and fruits) from a five-year-old apple tree were collected.In vitroshoot cultures of apple Gala were maintained on Murashige and Skoog (MS) medium containing 0.5 mg L–16-BA, 0.2 mg L–1NAA, and 0.1 mg L–1gibberellin (GA)under long-day conditions (16 h light/8 h dark cycle) and then subcultured with a 4 weeks interval (Anet al. 2016).Two-week-old apple shoot cultures were treated separately with either 0.5 mmol L–1nitrates, 200 mmol L–1NaCl, 300 mmol L–1mannitol, 10% PEG, 4°C temperature, or 100 μmol L–1ABA. Untreated apple tissue cultures were used as controls. The samples were treated for 6 h before the expression analysis. All the samples were stored at ?80°C.RNA was extracted from triplicate biological replicates of the above samples using the Trizol method (Caoet al. 2013).All of the primers are shown in Appendix A. The results were based on the average of three parallel experiments.

    3. Results

    3.1. ldentification of CEP genes in apple

    In total, 27 typicalCEPgenes were identified using theArabidopsisCEP domain as the unique identifier after manually removing the repeated and incomplete sequences from our data. These genes contained full length open reading frames (ORFs). Each gene was named based on its location on the chromosome (Table 1). The amino acids from the 89th(MdCEP26) to the 1 125th (MdCEP1) were encoded by the identifiedMdCEPgenes, which predicted that the molecular masses were 9.1-125.4 kDa and the pI of these proteins was between 5.58 (MdCEP3) and 10.62 (MdCEP25).

    3.2. Phylogenetic tree of the apple CEP proteins

    To evaluate the evolutionary relationships among the MdCEP proteins, the phylogenetic analysis was performed for theArabidopsis(15 members) and apple (27 members)CEP proteins using NJ method of the MEGA5 Program(Fig. 1). Next, the unrooted phylogenetic tree was constructed using the multiple sequence alignment file. The phylogenetic tree of the MdCEPs was divided into three classes (I, II and III) of monophyletic clades and 7 sister pairs of paralogousMdCEPgenes were found with strong bootstrap support (>90%). In addition, Class I contained sixMdCEPgenes. Class III contained 10MdCEPgenes, and class II constituted the largest clade, with 11MdCEPgenes.The CEP protein in the same group has a similar amino acid sequence of the CEP domain (Appendix B).

    3.3. Chromosomal localizations of CEPs in the apple genomes

    The chromosomal location for eachMdCEPgene was analyzed to visualize the relationships among genes. The results showed that the genes were distributed unevenly. Among these 27MdCEPgenes, threeMdCEPgenes (MdCEP25,MdCEP26andMdCEP27) were localized to unassembled genomic sequence scaffolds and could not be conclusively mapped to any chromosome. The remaining 24MdCEPgenes were present on 10 of 17 chromosomes (Fig. 2). Notably, chromosome 1 encompassed the largest quantity with sixMdCEPgenes, whereas chromosomes 10, 11 and 12 contained only oneMdCEPgene. For the remaining genes,threeMdCEPgenes were distributed on chromosomes 1,5 and 16, and two genes were located on chromosomes 2,6 and 8. Interestingly, threeMdCEPgenes (MdCEP1,Md-CEP2, andMdCEP3) and three pairs ofMdCEPs(MdCEP11andMdCEP12, MdCEP17andMdCEP18,MdCEP22andMdCEP23) were tightly colocalized on the apple genome. In this study, the large-scale segmental duplication events were investigated, which were thought to have occurred during the evolutionary process. Among theMdCEPgenes, a total of 7 genes (MdCEP6,MdCEP13,MdCEP14,MdCEP15,MdCEP16,MdCEP20, andMdCEP21) were found in the segmental duplication blocks. We also noted thatMdCEP6andMdCEP13were on the same pair of homologous chromosomes (Chr. 5 and Chr. 10) and thatMdCEP11,MdCEP12,MdCEP17,MdCEP18, andMdCEP19were on another pair of homologous chromosomes (Chr. 8 and Chr. 15). Additionally,they have been classified into the class II phylogenetic group.This confirms that genome duplications have contributed to gene family expansions throughout plant evolution (Cannonet al. 2004; Liet al. 2015).

    Table 1 Information of C-TERMINALLYENCODED PEPTIDE (CEP) gene family in apple

    3.4. Gene structure of CEP genes

    To analyze the structural diversity and understand the possible structural evolution ofMdCEPgenes, the exon/intron organization in the coding sequences of the individualMdCEPgenes were analyzed. As illustrated in Fig. 3, the number of introns varied from 0 to 21 in theMdCEPgene family. Among the 27MdCEPgenes, there were 16 genes with no introns,and four genes (MdCEP11,MdCEP14,MdCEP20 andMd-CEP27) with only one intron. Unlike otherMdCEPs,MdCEP1,MdCEP2, andMdCEP3had many more introns. It should be noted that closely related members in the same subfamilies share similar gene length and exon/intron structure.

    3.5. Protein motif analysis of the apple CEP gene family

    To reveal the diversification of theMdCEPgenes, all of the MdCEP protein sequences were predicted by the MEME Program and 10 types of motifs were identified (Fig. 4). As expected, a common motif composition of the closely related members in the phylogenetic tree was revealed, suggesting functional similarities among the MdCEP proteins within the same subgroup (Fig. 4-A). The results showed that the CEP domain (motif 1) contained 15 amino acids and was shared by all members. In addition, three MdCEPs(MdCEP11, MdCEP12 and MdCEP19) were composed of more than five CEP domains; MdCEP14 and MdCEP17 had two CEP domains and the remaining MdCEPs had just one CEP domain (Fig. 4-B). Moreover, motif 5 existed in mostCEPgenes and it was predicted by the SignalP4.1 (http://www.cbs.dtu.dk/services/SignalP) to be the signal peptide which possibly contributed to the functions of CEP (Delayet al. 2013; Suiet al. 2016). Next, the WebLogo plots were assembled, and revealed that motifs 1 and 5 were highly conserved (Fig. 4-C), indicating that these sequences might have important biological functions.

    3.6. Expression profiles of MdCEP genes in different tissues and in response to different abiotic conditions

    Fig. 1 Phylogenetic relationship between Arabidopsis and apple C-TERMINALLYENCODED PEPTIDE (CEP) proteins. The protein sequences of the CEPs domain, including 27 from apple (MdCEPs, ★) and 15 from Arabidopsis (AtCEPs, ☆), were aligned by ClustalX (ver. 2.1) and were classified into 3 distinct classes (I, II and III). The phylogenetic tree was constructed using the neighborjoining (NJ) method by the MEGA (ver. 5.1) Software. The bootstrap values of 1 000 replicates were calculated at each node.

    Fig. 2 Chromosome positions of C-TERMINALLYENCODED PEPTIDE (CEP) genes in apple (Malus×domestica, MdCEP). The chromosome number is indicated at the 10 of each chromosome representation. To simplify the presentation, we named the putative CEP family genes from MdCEP1 to MdCEP27 according to the gene order on chromosomes (MdCEP25, MdCEP26 and MdCEP27 were not found). Segmented duplicate homologous blocks are indicated with shadows.

    Fig. 3 Exon/intron structure of C-TERMINALLYENCODED PEPTIDE (CEP) genes in apple (Malus×domestica, MdCEP).

    Gene expression patterns can provide important clues for gene functions (Kidokoroet al. 2015). The expression patterns ofMdCEPgenes were detected by quantitative real-time polymerase chain reaction (qRT-PCR) in five different organs (roots, stems, leaves, flowers, and fruits).The results showed that the tissue expression patterns of theCEPgenes were specific and expressed divergently(Fig. 5-A). Compared to other genes,MdCEP1,MdCEP6,MdCEP8andMdCEP16had a relatively high expression levels in stems, whereas,MdCEP9andMdCEP10were highly expressed in flowers and had lower levels of expression in fruits. It was noteworthy thatMdCEP14was mainly expressed in the leaves and thatMdCEP25was chiefly expressed in the fruits. Conversely, the transcript levels of the remainingMdCEPgenes were significantly higher in roots, withMdCEP22having the highest levels. Mounting evidence concludes thatCEPsexpression is regulated by environmental cues inArabidopsis(Delayet al. 2013). To examine whether the expression of appleMdCEPgenes were also induced by environmental cues and ABA treatment, the transcript abundance of all of theMdCEPgenes were tested under nitrogen limitation, salinity (NaCl) osmotic(Mannitol, PEG), cold (4°C) and exogenous ABA stress for 6 h. The results showed that expression ofMdCEPgenes were regulated by environmental factors and following ABA treatment (Fig. 5-B). The transcript levels ofMdCEP5,MdCEP6,MdCEP8andMdCEP10increased in response to the nitrogen limitation, while the expression of one-third ofMdCEPsgenes was decreased. The levels ofMdCEP1,MdCEP9,MdCEP13,MdCEP16andMdCEP27were enhanced in response to NaCl treatment. The expression ofMdCEP3,MdCEP5,MdCEP6,MdCEP8andMdCEP25was primarily induced by mannitol and PEG. MoreoverMdCEP15,MdCEP19andMdCEP20were significantly up-regulated under cold treatment. Notably,MdCEP17had the highest expression after ABA treatment. In addition, the genes within the same classes showed similar expression patterns. For example, members of the class II were induced by the ABA treatment. Most of theMdCEPgenes that were significantly up-regulated under low nitrogen conditions were classified as class III, while class I genes, includingMdCEP14,MdCEP21,MdCEP22, andMdCEP23,did not show significant changes. Interestingly, the expressions ofMdCEP21,MdCEP22,MdCEP23andMdCEP24were up-regulated by PEG treatment, but inhibited by the other treatments.

    Fig. 4 Phylogenetic tree (A), protein motif analysis (B) and conserved domains (C) of C-TERMINALLYENCODED PEPTIDE (CEP)genes in apple (Malus×domestica, MdCEP). The amino acid sequences of the MdCEP proteins were aligned using Clustal X,and the phylogenetic tree was constructed using the neighbor-joining (NJ) method in the MEGA5 (ver. 5.1) Software. Schematic representation of the conserved motifs in the CEP proteins, which were elucidated using the MEME (ver. 4.11.4) Software. The black lines represent the unconserved sequences. The scale bar represents 200 aa.

    4. Discussion

    Fig. 5 Expression patterns of C-TERMINALLYENCODED PEPTIDE (CEP) genes of apple (Malus×domestica, MdCEP) in different organs (A) and under various abiotic conditions (B). Quantitative real-time polymerase chain reaction (qRT-PCR) was used to analyze the expression patterns of all the MdCEP genes in seedlings under the following conditions, nitrates (0.5 mmol L-1), salinity(NaCl, 200 mmol L–1), mannitol (300 mmol L–1), drought (PEG 6 000, 10% ), cold (4°C), and exogenous ABA (100 μmol L–1) for 6 h. The data were normalized to the apple actin expression level.

    The CEPs, a class of small post-translationally modified signaling peptides, have been shown to be involved in regulating root architecture inArabidopsisand cell size in rice seeds (Robertset al. 2016; Suiet al. 2016). The genome-wide analyses of theCEPgene family have been reported inArabidopsisand rice. TheCEPfamily members displayed distinct expression patterns inArabidopsisdevelopment and were expressed divergently and exhibited tissue-specific expression patterns inOryza sativa(Delayet al. 2013; Suiet al. 2016). However, the biological functions for most CEPs remain largely unknown (Lee and De Smet 2016; Okamotoet al. 2016). In this study, twenty-sevenCEPgenes that contained the conserved CEP peptide domain were identified in apple. Then, a phylogenetic tree ofMdCEPsandAtCEPswas constructed and divided into three classes (class I, II and III) of monophyletic clades.Among them, the class III genes contained tenMdCEPs(but noAtCEPs), indicating that the appleCEPgenes have expanded during the evolutionary process (Cannonet al.2004). It is speculated thatCEPgenes may play a more complicated and important role in perennial woody plants(Delayet al. 2013; Robertset al. 2013). In addition, the gene structure ofCEPand the amino acid sequence of CEP domain in different groups were significantly different. For example, the CEP domain in class I contains more proline and glycine. The C-terminal GVG sequence of the CEP domain in class II is much conserved. The N-terminal part of most CEP domain (AYGQTGTGK) in class III showed high amino acid similarity.

    The 27MdCEPgenes were localized to 10 of the 17 apple chromosomes. However, 3MdCEPgenes (MdCEP25,MdCEP26,andMdCEP27) were not conclusively mapped to any chromosome. It is hypothesized that otherMdCEPgenes exist in unknown genomic gaps, because only 81.3%of genomic sequences overlapped in the apple genome (Velascoet al. 2010; Wanget al. 2013). The protein motifs of the appleCEPgene family were analyzed using the MEME Program. As predicted, all MdCEPs contained more than one CEP domain (motif 1) and most MdCEPs contained motif 5 which may be the signal peptide. As shown in Fig. 4-C, the sequences of motifs 1 and 5 were highly conserved. Interestingly, the CEP domains inArabidopsisand apple both contain 15 amino acids and are highly conserved,which means that the genetic evolution of theCEPgenes structures is highly conserved in different plants.

    Gene expression patterns provide important clues of gene function (Kidokoroet al. 2015). AllMdCEPgenes displayed different expression patterns, which represented the distinct roles of individualMdCEPgenes. In particularly,over two-thirds of theMdCEPgenes were highly expressed in the roots. It is hypothesized thatMdCEPgenes might have an indispensable role in regulating the growth and development of apple roots. Previous research found thatCEPgenes inhibit primary root length and lateral root number inArabidopsis(Robertset al. 2013; Robertset al. 2016).InMedicago, ectopic expression ofMtCEP1also reduces the lateral root number, but does not show any significant phenotypic variation after the knockdown ofMtCEP1alone.These data indicate thatCEPsare subject to functional redundancy (Iminet al. 2013; Mohd-Radzmanet al. 2015).

    Plants are constantly under multiple abiotic stresses during their lifespan and have evolved unique strategies to integrate these cues into growth regulation to adapt to adverse environments (Glazebrook 2001; Chenet al. 2002;Zhaoet al. 2012; Suet al. 2013). Research has showed that the small secreted regulatory CEPs play prominent roles in adapting to adverse environmental conditions (Delayet al.2013). Recent studies indicate that CEPs are negative regulators of plant development and that over-expression ofCEPgenes occurs under nutrient or abiotic stress, particularly under nitrogen limitation (Tabataet al. 2014). In this study, we explored the expression ofMdCEPgenes under various growth conditions including nitrogen limitation, NaCl,mannitol, PEG, cold (4°C) and exogenous ABA. The results showed thatMdCEPgenes were induced to varying degrees under different treatments. Notably, nearly two-thirds ofMdCEPgenes were involved in the low nitrogen response,which is consistent with previous reports inArabidopsisandMedicago(Delayet al. 2013; Iminet al. 2013; Tabataet al.2014). Presumably these genes may have similar functions on vegetative growth in plants. MostMdCEPgenes were also significantly up-regulated under NaCl, mannitol, PEG and ABA treatments, indicating thatMdCEPgenes might participate in the drought stress response. Interestingly,the transcript levels ofMdCEP15,MdCEP19andMdCEP20were significantly higher under cold treatment, indicating that these genes may have similar functions in regulating the response to cold stress. These functions are achieved by interacting with membrane-bound receptors on target cells(Etchells and Turner 2010; Hirakawaet al. 2010; Okamotoet al. 2016; Robertset al. 2016). In addition, someMdCEPgenes are associated withAtCEPgenes in the same evolutionary branch and have similar expression patterns. For example, bothMdCEP26andAtCEP5respond significantly to nitrogen limitation.MdCEP25andAtCEP4have a significant response to mannitol treatment (Delayet al. 2013).

    In this study, a bioinformatics approach was applied to perform a genome-wide survey ofCEPgenes in apple. Our data suggest that the 27 CEP proteins of apple form three subgroups according to their CEP domain. In addition, theCEPgene in the same subfamily has similar gene structure and amino acid sequence of CEP domain. The expression patterns ofMdCEPgenes offer a comprehensive framework to analyze the important role of the CEP peptides in plant growth and development. Our results show that the CEP domain sequences may influence gene function. For exam-ple, the expression ofCEPgene in roots was significantly higher in class I. Meanwhile, theMdCEPgenes of class I were induced by PEG treatment specifically. These results constitute the basis for further exploration of the mechanism of theCEPgene family, and indicate that the CEP domain is directly related to the biological function of theCEPgene with higher plants. In addition, peptide signals can be matched to their corresponding receptor through a variety of methods (Murphyet al. 2012). However, only very few peptide ligand-receptors have been shown to pair with corresponding peptides among a large number of possible peptide ligand-receptors (kinases) (Butenkoet al.2009; Songet al. 2016). Additionally, it was shown that one receptor-like protein kinase (RLK) can bind several signaling molecules and that some receptors may functionally replace one another (Murphyet al. 2012). It is unknown why and how these proteins act interchangeably. What we do know regarding peptide signaling potential, including the interaction with receptors and any downstream changes is just the tip of the iceberg. There are still many mysteries that need to be solved.

    5. Conclusion

    In this study, 27 appleCEPgenes were classified. Additionally, we presented a complete analysis of theMdCEPgene family utilizing the apple genome database. Next,qRT-PCR was used to determine the spatial expression ofMdCEPgenes, which exhibited tissue-specific expression patterns. Additionally, the expression level ofMdCEPgenes was increased by nitrogen limitation and a variety of other abiotic stresses. These data may facilitate the functional characterization of theCEPgene family.

    Acknowledgements

    We thank the earmarked fund for the China Agriculture Research System (CARS-28), the National Natural Science Foundation of China (31471854, 31601742), the earmarked fund for the Shandong Province Agriculture Research,China (SDAIT-06-03), the Natural Science Foundation of Shandong Province, China (ZR2011CQ007).

    Appendicesassociated with this paper can be available on http://www.ChinaAgriSci.com/V2/En/appendix.htm

    An J P, Li H H, Song L Q, Su L, Liu X, You C X, Hao Y J. 2016.The molecular cloning and functional characterization of MdMYC2, a bHLH transcription factor in apple.Plant Physiology and Biochemistry, 108, 24-31.

    Bryan A C, Obaidi A, Wierzba M, Tax F E. 2012. XYLEM INTERMIXED WITH PHLOEM1, a leucine-rich repeat receptor-like kinase required for stem growth and vascular development inArabidopsis thaliana.Planta, 235, 111-122.

    Butenko M A, Vie A K, Brembu T, Aalen R B, Bones A M. 2009.Plant peptides in signalling, looking for new partners.Trends in Plant Science, 14, 255-263.

    Cannon S B, Mitra A, Baumgarten A, Young N D, May G. 2004.The roles of segmental and tandem gene duplication in the evolution of large gene families inArabidopsis thaliana.BMC Plant Biology, 4, 10.

    Cao Z H, Zhang S Z, Wang R K, Zhang R F, Hao Y J. 2013.Genome wide analysis of the apple myb transcription factor family allows the identification ofMdoMYB121gene confering abiotic stress tolerance in plants.PLOS ONE,8, e69955.

    Chen W, Provart N J, Glazebrook J, Katagiri F, Chang H S,Eulgem T, Mauch F, Luan S, Zou G, Whitham S A, Budworth P R, Tao Y, Xie Z, Chen X, Lam S, Kreps J A, Harper J F,Ammour A S, Mani B M, Heinlein M, Kobayashi K, Hohn T,Dangl J L, Wang X, Zhua T,et al. 2002. Expression profile matrix ofArabidopsistranscription factor genes suggests their putative functions in response to environmental stresses.The Plant Cell, 14, 559–574.

    Czyzewicz N, Yue K, Beeckman T, De Smet I. 2013. Message in a bottle, small signalling peptide outputs during growth and development.Journal of Experimental Botany, 64,5281–5296.

    Delay C, Imin N, Djordjevic M A. 2013.CEPgenes regulate root and shoot development in response to environmental cues and are specific to seed plants.Journal of Experimental Botany, 64, 5383–5394.

    Dimick P S, Hoskin J C. 1983. Review of apple flavor - state of the art.Critical Reviews in Food Science and Nutrition,18, 387–409.

    Du C, Li X, Chen J, Chen W, Li B, Li C, Wang L, Li J, Zhao X, Lin J. 2016. Receptor kinase complex transmits RALF peptide signal to inhibit root growth inArabidopsis.Proceedings of the National Academy of Sciences of the United States of America, 113, 8326–8334.

    Edgar R C. 2004. MUSCLE, multiple sequence alignment with high accuracy and high throughput.Nucleic Acids Research,32, 1792–1797.

    Etchells J P, Turner S R. 2010. The PXY-CLE41 receptor ligand pair defines a multifunctional pathway that controls the rate and orientation of vascular cell division.Development, 137,767–774.

    Glazebrook J. 2001. Genes controlling expression of defense responses inArabidopsis- 2001 status.Current Opinion in Plant Biology, 4, 301–308.

    Hall B G. 2013. Building phylogenetic trees from molecular data with MEGA.Molecular Biology and Evolution, 30,1229–1235.

    Hirakawa Y, Kondo Y, Fukuda H. 2010. Establishment and maintenance of vascular cell communities through local signaling.Current Opinion in Plant Biology, 14, 17–23.

    Imin N, Mohd-Radzman N A, Ogilvie H A, Djordjevic M A. 2013.The peptide-encodingCEP1gene modulates lateral root and nodule numbers inMedicago truncatula.Journal of Experimental Botany, 64, 5395–5409.

    Katsir L, Davies K A, Bergmann D C, Laux T. 2011. Peptide signaling in plant development.Current Biology, 21,356–364.

    Kidokoro S, Watanabe K, Ohori T, Moriwaki T, Maruyama K,Mizoi J, Myint P S H N, Fujita Y, Sekita S, Shinozaki K. 2015.Soybean DREB1/CBF-type transcription factors function in heat and drought as well as cold stress-responsive gene expression.The Plant Journal, 81, 505–518.

    Kondo Y, Hirakawa Y, Kieber J J, Fukuda H. 2011. CLE peptides can negatively regulate protoxylem vessel formationviacytokinin signaling.Plant & Cell Physiology, 52, 37–48.

    Lee J S, De Smet I. 2016. Fine-tuning development through antagonistic peptides, an emerging theme.Trends in Plant Science, 21, 991–993.

    Lee J S, Hnilova M, Maes M, Lin Y C, Putarjunan A, Han S K,Avila J, Torii K U. 2015. Competitive binding of antagonistic peptides fine-tunes stomatal patterning.Nature, 522,439–443.

    Lee S Y, Chang S S, Shin J H, Kang D H. 2007. Membranefiltration method for enumeration and isolation ofAlicyclobacillusspp. from apple juice.Letters in Applied Microbiology, 45, 540–546.

    Li X, Yin X, Wang H, Li J, Guo C. 2015. Genome-wide identification and analysis of the apple (Malus×domesticaBorkh.)TIFYgene family.Tree Genetics & Genomes, 11,1–13.

    Mohd-Radzman N A, Binos S, Truong T T, Imin N, Mariani M,Djordjevic M A. 2015. Novel MtCEP1 peptides producedin vivodifferentially regulate root development inMedicago truncatula.Journal of Experimental Botany, 66, 5289–5300.

    Mohd-Radzman N A, Laffont C. 2016. Different Pathways act downstream of the CEP peptide receptor CRA2 to regulate lateral root and nodule development.Plant Physiology,171, 2536–2548.

    Murphy E, Smith S, De Smet I. 2012. Small signaling peptides inArabidopsisdevelopment, how cells communicate over a short distance.The Plant Cell, 24, 3198–3217.

    Ogilvie H A, Imin N, Djordjevic M A. 2014. Diversification of the C-TERMINALLY ENCODED PEPTIDE (CEP) gene family in angiosperms, and evolution of plant-family specific CEP genes.BMC Genomics, 15, 870.

    Ohyama K, Ogawa M, Matsubayashi Y. 2008. Identification of a biologically active, small, secreted peptide inArabidopsisbyin silicogene screening, followed by LC-MS-based structure analysis.The Plant Journal, 55, 152–160.

    Okamoto S, Tabata R, Matsubayashi Y. 2016. Long-distance peptide signaling essential for nutrient homeostasis in plants.Current Opinion in Plant Biology, 34, 35–40.

    Roberts I, Smith S, De Rybel B, Van Den Broeke J, Smet W, De Cokere S, Mispelaere M, De Smet I, Beeckman T. 2013. The CEP family in land plants, evolutionary analyses, expression studies, and role inArabidopsisshoot development.Journal of Experimental Botany, 64, 5371–5381.

    Roberts I, Smith S, Stes E, De Rybel B, Staes A, Van De Cotte B, Njo M F, Dedeyne L, Demol H, Lavenus J, Audenaert D,Gevaert K, Beeckman T, Smet L D. 2016. CEP5 and XIP1/CEPR1 regulate lateral root initiation inArabidopsis.Journal of Experimental Botany, 67, 4889–4899.

    Song W, Liu L, Wang J, Wu Z, Zhang H, Tang J, Lin G, Wang Y, Wen X, Li W, Han Z, Guo H, Chai J. 2016. Signature motif-guided identification of receptors for peptide hormones essential for root meristem growth.Cell Research, 26,674–685.

    Su H, Zhang S, Yuan X, Chen C, Wang X F, Hao Y J. 2013.Genome-wide analysis and identification of stressresponsive genes of the NAM-ATAF1, 2-CUC2 transcription factor family in apple.Plant Physiology and Biochemistry,71, 11–21.

    Sui Z, Wang T, Li H, Zhang M, Li Y, Xu R, Xing G, Ni Z, Xin M. 2016. Overexpression of peptide-encodingOsCEP6.1results in pleiotropic effects on growth in rice (O. sativa).Frontiers in Plant Science, 7, 228.

    Tabata R, Sumida K, Yoshii T, Ohyama K, Shinohara H,Matsubayashi Y. 2014. Perception of root-derived peptides by shoot LRR-RKs mediates systemic N-demand signaling.Science, 346, 343–346.

    Velasco R, Zharkikh A, Affourtit J, Dhingra A, Cestaro A,Kalyanaraman A, Fontana P, Bhatnagar S K, Troggio M,Pruss D, Salvi S, Pindo M, Baldi P, Castelletti S, Cavaiuolo M, Coppola G, Costa F, Cova V, Dal Ri A, Goremykin V,et al. 2010. The genome of the domesticated apple(Malus×domesticaBorkh.).Nature Genetics, 42, 833–839.

    Wang J, Guo M Z, Xing L L. 2012. FastJoin, an improved neighbor-joining algorithm.Genetics and Molecular Research, 11, 1909–1922.

    Wang X, Zhang S, Su L, Liu X, Hao Y. 2013. A genome-wide analysis of the LBD (LATERAL ORGAN BOUNDARIES Domain) gene family inMalusdomesticawith a functional characterization ofMdLBD11.PLOS ONE, 8, e57044.

    Xu J, Xing S, Zhang Z, Chen X, Wang X. 2016. Genome-wide identification and expression analysis of the tubby-like protein family in theMalusdomesticagenome.Frontiers in Plant Science, 7, e25198.

    Zhao T, Liang D, Wang P, Liu J, Ma F. 2012. Genome-wide analysis and expression profiling of the DREB transcription factor gene family inMalusunder abiotic stress.Molecular Genetics and Genomics, 287, 423–436.

    黑人猛操日本美女一级片| 日日摸夜夜添夜夜添小说| 国产成人av激情在线播放| 最近最新中文字幕大全电影3 | 欧美在线一区亚洲| 免费观看人在逋| 久久久久久久精品吃奶| 国产黄色免费在线视频| 在线观看免费日韩欧美大片| 久久午夜亚洲精品久久| 欧美日韩一级在线毛片| 两人在一起打扑克的视频| 中出人妻视频一区二区| 高清av免费在线| 伦理电影免费视频| 人人妻人人添人人爽欧美一区卜| 高潮久久久久久久久久久不卡| 免费观看精品视频网站| 亚洲 欧美一区二区三区| 韩国av一区二区三区四区| 手机成人av网站| 妹子高潮喷水视频| 好男人电影高清在线观看| 国产在线观看jvid| 日韩精品免费视频一区二区三区| 欧美午夜高清在线| 国产亚洲欧美在线一区二区| 男女免费视频国产| 男人的好看免费观看在线视频 | 一a级毛片在线观看| 欧美午夜高清在线| 女人久久www免费人成看片| 久久国产精品大桥未久av| 啪啪无遮挡十八禁网站| 一区在线观看完整版| 女人精品久久久久毛片| 高清在线国产一区| 黄片播放在线免费| 在线观看舔阴道视频| 亚洲精品一卡2卡三卡4卡5卡| 婷婷成人精品国产| 精品久久久久久,| 欧美午夜高清在线| 极品教师在线免费播放| 午夜福利在线观看吧| 一进一出抽搐动态| 熟女少妇亚洲综合色aaa.| 人人妻人人添人人爽欧美一区卜| 无限看片的www在线观看| 在线观看66精品国产| 欧美亚洲日本最大视频资源| 亚洲美女黄片视频| 1024香蕉在线观看| 91成年电影在线观看| 丝袜美腿诱惑在线| cao死你这个sao货| 别揉我奶头~嗯~啊~动态视频| 在线av久久热| 丰满迷人的少妇在线观看| 99国产精品99久久久久| 变态另类成人亚洲欧美熟女 | 亚洲国产中文字幕在线视频| 一级毛片女人18水好多| 国产免费男女视频| 国产亚洲一区二区精品| av中文乱码字幕在线| 国产av一区二区精品久久| 国产精品九九99| 国产成人精品久久二区二区91| 亚洲一区二区三区欧美精品| 一级作爱视频免费观看| 色综合欧美亚洲国产小说| 夜夜躁狠狠躁天天躁| 欧美人与性动交α欧美软件| 日韩中文字幕欧美一区二区| 少妇粗大呻吟视频| 亚洲国产看品久久| 高清视频免费观看一区二区| 久久人妻福利社区极品人妻图片| 国产黄色免费在线视频| 国产亚洲欧美98| 亚洲在线自拍视频| 国产免费男女视频| 久久婷婷成人综合色麻豆| 亚洲国产精品sss在线观看 | 国产一区有黄有色的免费视频| 精品乱码久久久久久99久播| 国产伦人伦偷精品视频| 成人黄色视频免费在线看| 中出人妻视频一区二区| 亚洲色图综合在线观看| a级毛片黄视频| 久久久久久久国产电影| 国产av一区二区精品久久| 国产男女超爽视频在线观看| 在线观看午夜福利视频| 国产单亲对白刺激| 黄色丝袜av网址大全| 亚洲国产欧美网| 午夜两性在线视频| 久久国产精品男人的天堂亚洲| 美女高潮到喷水免费观看| 国产成人精品无人区| 欧美人与性动交α欧美精品济南到| 91麻豆av在线| 久久 成人 亚洲| 王馨瑶露胸无遮挡在线观看| 欧美日韩国产mv在线观看视频| 黄色成人免费大全| 老熟妇乱子伦视频在线观看| 三上悠亚av全集在线观看| 99re在线观看精品视频| 美女高潮到喷水免费观看| 久久久久国内视频| 精品久久蜜臀av无| 久久 成人 亚洲| 亚洲熟女毛片儿| 日韩成人在线观看一区二区三区| 亚洲熟妇中文字幕五十中出 | 这个男人来自地球电影免费观看| 夜夜夜夜夜久久久久| av天堂久久9| 在线十欧美十亚洲十日本专区| 亚洲全国av大片| 国产欧美日韩一区二区三| 国产有黄有色有爽视频| 欧美精品一区二区免费开放| 老司机福利观看| 在线国产一区二区在线| 建设人人有责人人尽责人人享有的| 无遮挡黄片免费观看| 亚洲精品美女久久av网站| 国产成人精品无人区| 无人区码免费观看不卡| 久久影院123| 久久久久国产精品人妻aⅴ院 | 老司机午夜十八禁免费视频| 久久久国产一区二区| aaaaa片日本免费| 午夜成年电影在线免费观看| 咕卡用的链子| 老熟妇仑乱视频hdxx| 人妻一区二区av| 高清黄色对白视频在线免费看| 日韩有码中文字幕| 国产精品一区二区在线不卡| www日本在线高清视频| 亚洲专区中文字幕在线| 欧美黄色片欧美黄色片| 大型av网站在线播放| 亚洲人成伊人成综合网2020| 9色porny在线观看| 女性被躁到高潮视频| 日韩免费av在线播放| 欧美黑人欧美精品刺激| 久久亚洲真实| 久久久久视频综合| 一本综合久久免费| 丝袜在线中文字幕| 午夜福利,免费看| 久久亚洲真实| 婷婷丁香在线五月| netflix在线观看网站| 少妇粗大呻吟视频| 高清视频免费观看一区二区| 亚洲aⅴ乱码一区二区在线播放 | 男女下面插进去视频免费观看| 亚洲欧美激情综合另类| av国产精品久久久久影院| 黑人巨大精品欧美一区二区mp4| 精品国产一区二区三区四区第35| 日韩视频一区二区在线观看| 欧美日韩一级在线毛片| 黄色a级毛片大全视频| 亚洲欧美精品综合一区二区三区| 波多野结衣av一区二区av| 精品久久久久久久毛片微露脸| 亚洲色图 男人天堂 中文字幕| 国产高清videossex| av不卡在线播放| 日韩欧美国产一区二区入口| 国产91精品成人一区二区三区| 在线看a的网站| 99热网站在线观看| 亚洲av成人一区二区三| 老熟妇仑乱视频hdxx| 1024香蕉在线观看| 国产黄色免费在线视频| 国产精品 欧美亚洲| 精品福利观看| 美女午夜性视频免费| 欧美激情高清一区二区三区| x7x7x7水蜜桃| netflix在线观看网站| 大型黄色视频在线免费观看| 一二三四社区在线视频社区8| 久久精品亚洲av国产电影网| 最新美女视频免费是黄的| 亚洲专区国产一区二区| 亚洲精品久久成人aⅴ小说| 久久人妻熟女aⅴ| 在线av久久热| 一本综合久久免费| 日韩成人在线观看一区二区三区| 黄色视频,在线免费观看| 免费在线观看亚洲国产| 一级毛片高清免费大全| 国产高清视频在线播放一区| 久久久水蜜桃国产精品网| 一区在线观看完整版| 两个人看的免费小视频| 日韩中文字幕欧美一区二区| 亚洲第一青青草原| 国产精品偷伦视频观看了| 成人影院久久| 在线观看免费视频日本深夜| 国产免费av片在线观看野外av| 1024视频免费在线观看| 日韩欧美在线二视频 | 亚洲五月天丁香| www.精华液| 午夜激情av网站| 精品一品国产午夜福利视频| 大型av网站在线播放| 伦理电影免费视频| 手机成人av网站| 人妻 亚洲 视频| 麻豆av在线久日| 久久久久久人人人人人| a级片在线免费高清观看视频| 男人的好看免费观看在线视频 | 日本wwww免费看| 91在线观看av| 日本五十路高清| 一个人免费在线观看的高清视频| 久久久精品国产亚洲av高清涩受| 看黄色毛片网站| 女人高潮潮喷娇喘18禁视频| 黑丝袜美女国产一区| 夜夜夜夜夜久久久久| 在线十欧美十亚洲十日本专区| 一级a爱片免费观看的视频| 精品免费久久久久久久清纯 | av欧美777| 啪啪无遮挡十八禁网站| www.熟女人妻精品国产| 在线观看66精品国产| 99热只有精品国产| 麻豆国产av国片精品| 亚洲精品av麻豆狂野| 久久久久久久国产电影| 天天添夜夜摸| 成年女人毛片免费观看观看9 | 久热这里只有精品99| 久久天躁狠狠躁夜夜2o2o| 国内久久婷婷六月综合欲色啪| 少妇 在线观看| 久久精品国产清高在天天线| 久久久国产成人免费| 看黄色毛片网站| 建设人人有责人人尽责人人享有的| 天堂动漫精品| 午夜精品在线福利| 欧美日韩国产mv在线观看视频| 新久久久久国产一级毛片| 精品国产乱码久久久久久男人| 成人国语在线视频| 成年版毛片免费区| 黄片小视频在线播放| 99精品久久久久人妻精品| 99久久99久久久精品蜜桃| videos熟女内射| 精品福利永久在线观看| 国产人伦9x9x在线观看| 精品一区二区三区视频在线观看免费 | 欧美乱妇无乱码| 脱女人内裤的视频| 国产精品影院久久| 免费人成视频x8x8入口观看| 熟女少妇亚洲综合色aaa.| 国产亚洲av高清不卡| 国产精品二区激情视频| 欧美日韩一级在线毛片| 久久精品亚洲熟妇少妇任你| 亚洲美女黄片视频| 啪啪无遮挡十八禁网站| 精品第一国产精品| 12—13女人毛片做爰片一| 99精国产麻豆久久婷婷| 老司机福利观看| 一区二区三区激情视频| 日韩三级视频一区二区三区| x7x7x7水蜜桃| 午夜激情av网站| 亚洲精品美女久久久久99蜜臀| x7x7x7水蜜桃| 手机成人av网站| 男女床上黄色一级片免费看| 99久久国产精品久久久| 久久精品国产综合久久久| 国产高清激情床上av| 91国产中文字幕| 美女福利国产在线| 巨乳人妻的诱惑在线观看| 免费看十八禁软件| 亚洲 欧美一区二区三区| 午夜影院日韩av| 亚洲aⅴ乱码一区二区在线播放 | 两人在一起打扑克的视频| 美女国产高潮福利片在线看| 精品乱码久久久久久99久播| 美女 人体艺术 gogo| 欧美国产精品va在线观看不卡| 另类亚洲欧美激情| 人人妻人人澡人人看| 精品卡一卡二卡四卡免费| 国产又色又爽无遮挡免费看| 久久人妻熟女aⅴ| 丰满迷人的少妇在线观看| 日本五十路高清| av超薄肉色丝袜交足视频| 亚洲午夜理论影院| www.999成人在线观看| 久久性视频一级片| 欧美日韩成人在线一区二区| 欧美精品av麻豆av| 9热在线视频观看99| 国产精品一区二区精品视频观看| 女性生殖器流出的白浆| 免费av中文字幕在线| 精品亚洲成国产av| 免费看十八禁软件| 日韩免费高清中文字幕av| 亚洲成人国产一区在线观看| 日韩一卡2卡3卡4卡2021年| 亚洲专区国产一区二区| 久久精品人人爽人人爽视色| 欧美成人免费av一区二区三区 | 国产又色又爽无遮挡免费看| 久9热在线精品视频| 成人永久免费在线观看视频| 亚洲视频免费观看视频| 性少妇av在线| 黄网站色视频无遮挡免费观看| 大陆偷拍与自拍| 国产精品一区二区在线不卡| 国产日韩欧美亚洲二区| 最新的欧美精品一区二区| 亚洲av片天天在线观看| 热99久久久久精品小说推荐| 三级毛片av免费| 亚洲精品成人av观看孕妇| 色94色欧美一区二区| √禁漫天堂资源中文www| 天天躁夜夜躁狠狠躁躁| 人成视频在线观看免费观看| 一级毛片精品| av在线播放免费不卡| 制服诱惑二区| 亚洲国产精品合色在线| 老汉色av国产亚洲站长工具| 91麻豆av在线| 亚洲av成人不卡在线观看播放网| 久久精品亚洲av国产电影网| 国内毛片毛片毛片毛片毛片| 国产精品秋霞免费鲁丝片| 国产精品国产高清国产av | 亚洲国产中文字幕在线视频| 老熟妇乱子伦视频在线观看| 亚洲熟妇中文字幕五十中出 | 欧美日韩成人在线一区二区| 一级毛片精品| 乱人伦中国视频| 国产精品 国内视频| 免费在线观看视频国产中文字幕亚洲| 亚洲av熟女| 99国产精品一区二区三区| 18禁国产床啪视频网站| 男人的好看免费观看在线视频 | 国产一区二区三区视频了| 免费高清在线观看日韩| 韩国av一区二区三区四区| 99国产综合亚洲精品| 一本综合久久免费| 久久久久久久精品吃奶| 热99国产精品久久久久久7| 国产一区二区激情短视频| 一级作爱视频免费观看| 欧美日韩视频精品一区| 亚洲精品av麻豆狂野| 精品视频人人做人人爽| 黄色片一级片一级黄色片| 欧美不卡视频在线免费观看 | 日韩熟女老妇一区二区性免费视频| 人人妻人人澡人人看| 国产精品98久久久久久宅男小说| 人人妻人人澡人人爽人人夜夜| e午夜精品久久久久久久| 亚洲专区中文字幕在线| 视频在线观看一区二区三区| 91在线观看av| 人妻一区二区av| 国产精品美女特级片免费视频播放器 | 国产精品二区激情视频| svipshipincom国产片| 午夜福利欧美成人| 成人18禁在线播放| 99re6热这里在线精品视频| 窝窝影院91人妻| 国产精品久久久人人做人人爽| 久久午夜综合久久蜜桃| 亚洲精品粉嫩美女一区| 精品国产国语对白av| 欧美激情高清一区二区三区| 在线十欧美十亚洲十日本专区| 亚洲av欧美aⅴ国产| 欧美日韩国产mv在线观看视频| 日韩三级视频一区二区三区| 国产高清videossex| 亚洲第一av免费看| 免费在线观看黄色视频的| 少妇 在线观看| 狠狠狠狠99中文字幕| 色综合婷婷激情| www.999成人在线观看| 欧美中文综合在线视频| 精品人妻1区二区| 国产色视频综合| 一级a爱视频在线免费观看| 亚洲性夜色夜夜综合| 看免费av毛片| 91麻豆精品激情在线观看国产 | 黑人猛操日本美女一级片| 精品亚洲成a人片在线观看| 制服诱惑二区| 在线观看66精品国产| 精品免费久久久久久久清纯 | 免费高清在线观看日韩| 精品国产乱子伦一区二区三区| 国产蜜桃级精品一区二区三区 | 国产一区有黄有色的免费视频| 操出白浆在线播放| 亚洲成人免费电影在线观看| 亚洲精品中文字幕一二三四区| 丁香六月欧美| 欧美老熟妇乱子伦牲交| 国产亚洲欧美精品永久| 女人高潮潮喷娇喘18禁视频| 操出白浆在线播放| 不卡一级毛片| 久久精品91无色码中文字幕| 又黄又粗又硬又大视频| 久久久久久久午夜电影 | 午夜福利影视在线免费观看| 国产免费现黄频在线看| 国产1区2区3区精品| 黄色视频不卡| xxxhd国产人妻xxx| 日韩三级视频一区二区三区| 久久久国产一区二区| 久久久精品国产亚洲av高清涩受| 激情视频va一区二区三区| 国产精品久久视频播放| a级毛片黄视频| 亚洲欧美一区二区三区黑人| 国产在视频线精品| 久久精品人人爽人人爽视色| 脱女人内裤的视频| 99精品欧美一区二区三区四区| 黑人猛操日本美女一级片| 啦啦啦免费观看视频1| 欧美黄色淫秽网站| 新久久久久国产一级毛片| 精品国产一区二区久久| 18禁观看日本| 久99久视频精品免费| 欧美亚洲日本最大视频资源| 成人手机av| 丰满饥渴人妻一区二区三| 亚洲精华国产精华精| 国产成+人综合+亚洲专区| 亚洲av电影在线进入| 亚洲成人手机| 色94色欧美一区二区| 美女福利国产在线| 国产精品.久久久| 欧美最黄视频在线播放免费 | 每晚都被弄得嗷嗷叫到高潮| 一区二区三区精品91| 天堂俺去俺来也www色官网| 最近最新免费中文字幕在线| 国产精品99久久99久久久不卡| av网站免费在线观看视频| 色综合婷婷激情| 热99久久久久精品小说推荐| 亚洲avbb在线观看| 在线视频色国产色| 欧美乱妇无乱码| 97人妻天天添夜夜摸| a级毛片在线看网站| av网站在线播放免费| 美女午夜性视频免费| 亚洲综合色网址| 啦啦啦免费观看视频1| 老司机福利观看| 国产淫语在线视频| 老司机福利观看| 亚洲片人在线观看| svipshipincom国产片| 在线永久观看黄色视频| 国内久久婷婷六月综合欲色啪| 久久精品aⅴ一区二区三区四区| av网站在线播放免费| avwww免费| 男女高潮啪啪啪动态图| 两个人免费观看高清视频| 在线播放国产精品三级| 久久精品aⅴ一区二区三区四区| 欧美国产精品一级二级三级| 国产精品国产av在线观看| av一本久久久久| 国产精品久久久人人做人人爽| 一边摸一边抽搐一进一小说 | 国产乱人伦免费视频| 久久99一区二区三区| 一级片'在线观看视频| 色老头精品视频在线观看| 欧美在线一区亚洲| 国产精品免费视频内射| 欧美日韩成人在线一区二区| 婷婷成人精品国产| 国产精品美女特级片免费视频播放器 | 亚洲欧美色中文字幕在线| 色综合欧美亚洲国产小说| 欧美黑人精品巨大| 国产人伦9x9x在线观看| 国产精品二区激情视频| 啦啦啦在线免费观看视频4| 亚洲av日韩精品久久久久久密| 中文字幕高清在线视频| 99精国产麻豆久久婷婷| 中文字幕精品免费在线观看视频| 久久人妻熟女aⅴ| 精品久久久精品久久久| 黑丝袜美女国产一区| 9191精品国产免费久久| 精品国产一区二区三区久久久樱花| 无限看片的www在线观看| 久久久久久久精品吃奶| 亚洲欧美日韩另类电影网站| 91国产中文字幕| 国产成人欧美| 一级,二级,三级黄色视频| 亚洲欧美日韩另类电影网站| 大香蕉久久成人网| 亚洲精品在线美女| 18禁黄网站禁片午夜丰满| www.自偷自拍.com| 国产亚洲欧美98| 又黄又爽又免费观看的视频| 国产主播在线观看一区二区| 高清av免费在线| 免费少妇av软件| 美女 人体艺术 gogo| 亚洲色图综合在线观看| 久久精品人人爽人人爽视色| 无限看片的www在线观看| 人成视频在线观看免费观看| 一级a爱视频在线免费观看| 18禁裸乳无遮挡免费网站照片 | 超碰成人久久| 国产免费现黄频在线看| 免费在线观看日本一区| 国产精品.久久久| 亚洲国产精品合色在线| 人人妻人人澡人人看| 国产一卡二卡三卡精品| 99国产精品一区二区三区| 亚洲精品自拍成人| 日韩欧美国产一区二区入口| 久久热在线av| 亚洲av第一区精品v没综合| 国产精品综合久久久久久久免费 | 午夜福利在线观看吧| 亚洲美女黄片视频| 看免费av毛片| 午夜福利在线观看吧| 免费黄频网站在线观看国产| 亚洲色图综合在线观看| 久久精品成人免费网站| 国产在线观看jvid| 99国产精品99久久久久| 欧美亚洲 丝袜 人妻 在线| 日本欧美视频一区| 啦啦啦在线免费观看视频4| 一边摸一边抽搐一进一小说 | 欧美久久黑人一区二区| 在线看a的网站| 久久精品国产a三级三级三级| 欧美大码av| 叶爱在线成人免费视频播放| 国产精品99久久99久久久不卡| 欧美精品一区二区免费开放| 精品人妻1区二区| 夫妻午夜视频| www.999成人在线观看| 搡老乐熟女国产| 99在线人妻在线中文字幕 | 亚洲国产精品一区二区三区在线| 一边摸一边做爽爽视频免费| 日韩有码中文字幕| 免费少妇av软件| 久久热在线av| 国产男女超爽视频在线观看| 欧美日韩精品网址|