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

    Molecular Cloning, Expression Pattern, and 3D Structural Prediction of the Cold Inducible RNA - Binding Protein (CIRP) in Japanese Flounder (Paralichthys olivaceus)

    2015-03-31 02:55:37YANGXiaoGAOJinningMALimanLIZanWANGWenjiWANGZhongkaiYUHaiyangQIJieWANGXuboWANGZhigangandZHANGQuanqi
    Journal of Ocean University of China 2015年1期

    YANG Xiao, GAO Jinning, MA Liman, LI Zan, WANG Wenji, WANG Zhongkai,YU Haiyang, QI Jie, WANG Xubo, WANG Zhigang, and ZHANG Quanqi

    ?

    Molecular Cloning, Expression Pattern, and 3D Structural Prediction of the Cold Inducible RNA-Binding Protein (CIRP) in Japanese Flounder ()

    YANG Xiao, GAO Jinning, MA Liman, LI Zan, WANG Wenji, WANG Zhongkai,YU Haiyang, QI Jie, WANG Xubo, WANG Zhigang, and ZHANG Quanqi*

    ,,,266003,..

    Cold-inducible RNA-binding protein (CIRP) is a kind of RNA binding proteins that plays important roles in many physiological processes. Thehas been widely studied in mammals and amphibians since it was first cloned from mammals. On the contrary, there are little reports in teleosts. In this study, the Pogene of the Japanese flounder was cloned and sequenced. The genomic sequence consists of seven exons and six introns. The putative PoCIRP protein of flounder was 198 amino acid residues long containing the RNA recognition motif (RRM). Phylogenetic analysis showed that the flounder PoCIRP is highly conserved with other teleost CIRPs. The 5’ flanking sequence was cloned by genome walking and many transcription factor binding sites were identified. There is a CpGs region located in promoter and exon I region and the methylation state is low. Quantitative real-time PCR analysis uncovered that Pogene was widely expressed in adult tissues with the highest expression level in the ovary. The mRNA of the Powas maternally deposited and the expression level of the gene was regulated up during the gastrula and neurula stages. In order to gain the information how the protein interacts with mRNA, we performed the modeling of the 3D structure of the flounder PoCIRP. The results showed a cleft existing the surface of the molecular. Taken together, the results indicate that the CIRP is a multifunctional molecular in teleosts and the findings about the structure provide valuable information for understanding the basis of this protein’s function.

    cold-inducible RNA-binding protein; CIRP;; expression pattern; 3D modeling; CpGs

    1 Introduction

    Cold-inducible RNA-binding protein (CIRP, also named CIRBP) was first cloned from mammals as a cold-shock protein by Nishiyama. ( 1997). The CIRP is one of the highly conserved glycine-rich RNA-binding proteins which contain an amino-terminal domain named consensus sequence of RNA-binding domain (CS-RBD) or RNA recognition motif (RRM) and a carboxyl-terminal glycine-rich domain (GRD) (Bergeron, 1993).

    CIRP is not only induced in response to mild hypothermia(Nishiyama, 1998a; Tong, 2013) but also involved in many other biological processes (Liu, 2013; Nishiyama., 1997). It could be induced by DNA damage through UV irradiation (Yang & Carrier, 2001), osmotic stress(Pan, 2004), low oxygen level (Wellmann, 2004), light (Sugimoto and Jiang, 2008), diseases (cancers and ischemia) (Hamid, 2003, Zhou, 2009) and damages (Xue, 1999){Yang, 2001 #14730}. In ad-dition to these stress induction, CIRP is expressed in brain associated with hibernation (Saito, 2000) and may play important roles in pronephrogenesis and neurogenesis during the embryonic period of(Uochi and Asashima, 1998). Recent study also found that it may modulate circadian gene expression., 2013, Morf, 2012) and is associated with circadian rhythm (Morf., 2012; Nishiyama, 1998b;Ripperger, 2000). In the mouse testis the expression ofduring spermatogenesis and after heat stress has special effects in the primary spermatocytes through regulating target mRNA stabilization (Nishiyama., 1998a; Xia, 2012). In, however, CIRP was found to be the major cytoplasmic RNA-bind- ing protein (Matsumoto, 2000). This may suggests that the CIRP may have different roles between the males and females.

    Although CIRPs have been explored in many species including mammals and amphibians, teleost CIRP and the methylation level of this gene have not been reported in detail. The CIRP has been reported to be induced by osmotic stress in salmon (Pan., 2004). Screening of EST-SSRs markers revealed that CIRP is related to growth traits of common carp (Ying, 2009). A cDNA-AFLP analysis identified the CIRP gene as one of female specific expressed sex-related genes in flatfish (Taboada, 2012). As is known that female flatfish grows faster and is bigger than the males, the related reports implied that the CIRP may be involved in this phenomenon. In order to better understand its biology functions in fishes and its relationship with sex, we cloned the full-length cDNA, the genomic sequence and the 5’ promoter region of CIRP gene from Japanese flounder(Po), a commercially important marine fish for both fisheries and aquaculture in East Asian countries. We quantified its expression in different adult tissues and at several early developmental stages to get knowledge about its possible functions during embryogenesis as well as adult growth. We also analyzed the potential regulatory motif and the methylation level of CpGs located in the promoter region. Phylogenetic reconstruction and 3D structural prediction were also applied to provide more information.

    2 Materials and Methods

    2.1 Animals and Sampling

    The fishes were sampled from Haiyang Yellow Sea Aquatic Product Co., Ltd, Yantai, Shandong, China. Six adults (three females and three males) were randomly selected from the tank. Eight tissues (heart, liver, spleen, kidney, whole bran, intestine, ovary, testis and skeletal muscle) were dissected and snapped frozen in liquid nitrogen and then kept at ?80℃ before using. Artificially fertilized eggs were incubated at 15℃ in tank with continuous aeration. Embryos of different developmental stages (including unfertilized eggs, one cell period, eight cells period, multi cells stages, high blastula, gastrula, neurula, tail-bud stage, heart-beating stage and hatching stage) were collected and directly preserved in 1mL RNAwait liquid (Solarbio, SR0020, Shanghai) and then kept at ?80℃ before use. A small piece of muscle tissue was saved at ?80℃ for genomic DNA extraction.

    2.2 Extraction of Genomic DNA and RNA

    Genomic DNA was extracted from muscle by phenol chloroform extracting method. Briefly, the tissues grinded in liquid nitrogen were added to 600μL TNES (125mmolL?1NaCl, 10mmolL?1Tris-HCl, 10mmolL?1EDTA, SDS 0.5%, Urea 4molL?1, pH7.4) buffer. Then the samples were digested overnight at 56℃ by 6μL proteinase K (20mg mL?1) at the final concentration of 200μgmL?1. The DNA was then isolated with phenol and chloroform/isoamyl alcohol, precipitated with isopropanol and washed with 70% ethanol. The isolated DNA were then dissolved in ddH2O, degraded with RNaseA at 37℃ for 1h, and then preserved at ?20℃ until use.

    Total RNAs were extracted from adult tissues and embryos of different developmental stages using Trizol reagent (Invitrogen) according to the manufacturer’s instruction and kept at ?80℃before use.

    2.3 cDNA Cloning and Sequence Analysis ofPo

    Total RNA (1μg) from the ovary tissue was reverse transcribed to synthesize the first strand cDNA using MMLV-Reverse-Transcriptase (RNase H?) (TaKaRa, Dalian, China). Primers for cloning the initial fragment of PomRNA were designed according to conserved sequences predicted in other teleosts (Table 1). The PCR program was as follows: 94℃ pre-denaturation for 10min; 30cycles of 94℃ denaturation for 30s, 51℃ annealing for 30s, 72℃ extension for 30s; and then final elongation at 72 for 10min.

    In order to get the full length of the cDNA sequence,gene specific primers for 5’ and 3’ RACE were designed according to the obtained sequence fragment (Table 1). We conducted the experiments to obtain the 5’ and 3’ unknown cDNA regions by using the SMART RACE PCR Amplification Kit (Clontech, CA , USA) following the manufacture’s protocol. The RACE PCR conditions were: 94℃ for 5min; 24 cycles of 94℃ for 30s, 58℃ for 30s, 72℃ for 3min; and a final elongation at 72℃ for 10min. The full-length of the coding sequence was amplified with two specifically designed primers (Table 1) using the following program: 94℃ for 5min; 30 cycles of 94℃ for 30s, 55℃ for 30s, 72℃ for 3min; and a final elongation at 72℃ for 10min.

    Table 1 Primers used in this study

    All the PCR reactions were performed using a GeneAmp PCR system 9700 (Applied Biosystems, Forster City, CA, USA), and the amplification products were separated by 1% agarose gel electrophoresis, purified by using the Gel DNA recovery Kit (Zymoclean, CA, USA), cloned to pMD18-T vector (TaKaRa, Dalian, China), and transformed to competentTrans5α strain (TransGenBiotech, Beijing, China). Positive clones containing the plasmid with inserts of the appropriate sizes were identified by colony PCR using M13Fw and M13Rv primers and sequenced.

    The ORF Finder (http://www.ncbi.nlm.nih.gov/gorf/ gorf.html) was used to predict the open reading frame (ORF).

    2.4 Analysis of Gene Sequence and Promoter of Po

    In order to know the genomic structure of Powe amplified the genomic DNA fragment to get the introns with primers as showed in Table 1. The PCR cycles were described as follows:94℃ for 5min; 30 cycles of 94℃ for 30s, 55℃ for 30s, 72℃ for 3min; and a final elongation at 72℃ for 10min. The promoter region was amplified by using Genome Walking Kit (TaKaRa, Dalian, China) with gene specific primers listed in Table 1. The PCR program was designed following the instructions of the Kit.

    The gene structure of the exon-intron was estimated by the Spidey program in NCBI homepage (http://ncbi.nlm. nih.gov/IEB/Research/Ostell/Spidey/). The putative transcription start site (TSS), TATA box and CAAT box were predicted by using the tools of Promoter 2.0 Prediction Server (http://www.cbs.dtu.dk/services/Promoter/) and SIGSCAN (http://www-bimas.cit.nih.gov/molbio/signal/). Alibaba v2.1 (http://www.gene-regulation.com/pub/programs/alibaba2/index.html) was used to predict the transcription factor (TF) binding sites. The CpGs island was predicted by EMBOSS Cpgplot (http://www.ebi.ac.uk/ Tools/seqstats/emboss_cpgplot/).

    2.5 Methylation Level of the CpGs in Promoter and Exon 1 Region

    The DNA samples of six tissues were extracted from 3 females and 3 males respectively using the method mentioned above. The samples were treated with the EZ DNA Methylation-GoldTMKit following the manual instruction, and used for PCR with the primers listed in Table 1. The amplicons were purified, cloned and sequenced as described above. For each sample, 10 positive clones were sequenced and the results were analyzed with BDPC (Bisulfite sequencing Data Presentation and Compilation) (http://services.ibc.uni-stuttgart.de/BDPC/).

    2.6 Quantitative Real-Time PCR (qPCR) of Po

    Total RNA extracted from different tissues of six healthy adults (three females and three males) and different developmental stages was used to analyze the expression pattern using qPCR. The embryos or unfertilized eggs at different development stages were pooled for qPCR with each pool containing 20 eggs or embryos. Three pools were prepared for each stage as triplicates. The cDNA were generated from RNA samples through the method described above. The primers (Table 1) used for qPCR were designed with Primer Express v3.0 (Applied Biosystems). The amplification conditions were performed according to the manual book of SYBR Premix?EX TaqTMII (TaKaRa, Bio Inc., Shiga, Japan). 18S rRNA was selected as the reference gene (Zhong, 2008) and used to normalize the mRNA levels of target genes in all samples. Negative control was performed at the same time but no template was added. The amplicfication was performed in an ABI Prism 7500 Sequence Detection System (Applied Biosystems, Forster, CA, USA). The program was as follows: 1 cycle of 95℃ for 2min; 40 cycles of 95℃ for 15s, 60℃ for 45s; 1cycle of 95℃ 15s, 60℃ 1min, 95℃ 15s, 60℃ 15s.

    We constructed standard curve with a serial dilution of solution which contained certain amount of plasmids inserted with either 18Sor the Pogene fragments. Experiment with each sample was conducted in triplicate and the means were used for further analyses using the relative standard curve.

    2.7 Sequence Alignment and Phylogenetic Reconstruction

    We collected homologous sequence of other species’ CIRPs, and made the multiple alignments of the CIRPs by ClustalX 2.1software (Larkin, 2007). The phylogenetic reconstruction and molecular evolutionary analyses were conducted by using the MEGA version 5 (Tamura, 2011). And the tree was constructed using the neighbor joining method with 1000 bootstrap replicates.

    2.8 Analysis of the PoProtein Sequence and Modeling of Its 3D Structure

    The protein sequence was analyzed for its molecular mass and theoretical isoelectric point using the Compute PI/Mw tool (http://web.expasy.org/compute_pi/). The P- fam (http://pfam.sanger.ac.uk/) databases were used to scan the domains in the protein sequence.

    Three-dimension structure of the Japanese flounder Poprotein was predicted using SWISS-MODEL server (http://www.expasy.org/swissmod/SWISS-MODEL. html) and the PDBsum Generate (http://www.ebi.ac.uk/ thornton-srv/ databases/ pdbsum/Generate.html) was used to predict the molecular surface and cleft in the molecular.

    2.9 Statistical Analysis

    The expression results were analyzed with SPSS19.0 software. Data reported were logarithmically transformed before statistical analysis. Significant differences between samples were determined by one-way or two-way ANO- VA.multiple comparisons were carried with LSD multiple range tests. Data obtained from between samples were considered significantly different when< 0.05. The data are displayed using the mean±standard error of the mean.

    3 Results and Discussions

    3.1 The Cloning of the mRNA Sequence of the Japanese Flounder(Po)

    A 597bp long cDNA fragment of the Pogene was cloned from ovary with the primers cirp-m-fw/rv. The unknown 5’ and 3’ regions of the cDNA fragments were obtained through the RACE reaction with the primers RACE-5-rv and RACE-3-fw, respectively. The full length of the cDNA sequence was obtained by assembling all the sequences according to the overlaps between the fragments and verified using the gene specific primers cirp-w-fw/rv. The mRNA sequence of the Pois of 1596bp with an 85bp 5’ UTR (untranslated region), a 914bp 3’UTR and a 597bp coding sequence (CDS) predicted by ORF Finder (Fig.1). A non-canonical polyade-nylation signal site (ATTAAA) was found in the 3’ UTR near the polyA site.

    Fig.1 The mRNA sequence of PoCIRP. CDS is shown in capital letters, whereas the 5’UTR and 3’UTR are shown in small letters. The deduced amino acid sequence is shown by single letter code of amino acids below the CDS. The non-canonical polyadenylation signal site (ATTAAA) is boxed.

    The predicted amino acid sequence of the PoCIRP consists of 198 amino acid residues with the calculated molecular mass of 20.6kDa. The PoCIRP was predicted to contain the RNA recognition motif (RRM_1, PF00076) from the 7th to 77th residues and several Arg-Gly-Gly (RGG) repeats in the glycine-rich region, which is the same as reported in salmon, human and mice (De Leeuw, 2007; Jin, 2009). BLASTP analysis showed that the PoCIRP was similar to other CIRPs of mammals, amphibians and teleosts. The sequence identities were 47%, 75%, 73%, 69%, compared to,,, respectively. The genomic, mRNA sequence and the predicted protein sequence were submitted to NCBI and the accession NO. is KJ739525.

    3.2 Phylogenetic Analysis of Poin Vertebrates

    In order to understand the relationships of PoCIRP among the vertebrate CIRPs, we constructed the phylogenetic tree with CIRPs of 16 sequences using human RBM3, a member of highly conserved family of RNA binding proteins, as outgroup (Fig.2). The CIRPs were divided into four clusters which agree with taxonomic status of these species, and the PoCIRP was located with- in the teleosts group with 76% bootstrap support, indicating that this protein is one of the fish CIRP orthologs and these proteins are highly conserved proteins.

    Fig.2 The phylogenetic tree based on amino sequence alignment for CIRPs in vertebrates. Bootstrap values (%) indicate 1000 replicates. Homo sapiens RBM3 was used as outgroup. The number of each internal branch is the local bootstrap probability. The GenBank accession number of the proteins are as follows: Homo sapiens (BAA11212.1), Pongo abelii (NP_001124692.1), Sus scrofa (NP_001233126.1), Rattus norvegicus (BAA19092.1), Mus musculus (BAA11213.1), Cricetulus griseus (NP_001231002.1), Gallus gallus (AGS43093.1), Xenopus laevis (A: NP_001079794.1, B: NP_ 001080069.1), Xenopus (Silurana) tropicalis (NP_001017228.1), Oncorhynchus mykiss (ACO08387.1), Salmo salar (ACI 33459.1), Danio rerio (AAI15178.1), Paralichthys olivaceus (KJ739525), Anoplopoma fimbria (ACQ58801.1), Osmerus mordax (ACO09027.1), RBM3 human (P98179.1).

    3.3 Identification and Analysis of the Genomic Structure of the Po

    Whole genome sequence of Powas obtained with specific primers (Table 1). The introns were identified by comparing the mRNA sequence with the gene fragment (Fig.3). The Pogene had six introns which were 661 bp, 110bp, 161bp, 95bp, 94bp and 516bp, respectively. All the introns were found to be spliced conforming to the GT-AT rule. The 5’ UTR was encoded by the exon I and a part of the exon II, the 3’ UTR was encoded by most of the exon VII and the CDS was encoded by the rest parts of the exons.

    Nearly about 4000bp of the upstream region of Powas sequenced by genome walking method. The putative transcription start site (TSS), TATA box and CAAT box predicted by Promoter 2.0 Prediction Server and SIGSCAN were shown in Fig.4. There was a high likely prediction of the TSS site upstream, the sequenced start site of the mRNA. As there was the possibility that the mRNA molecular was not complete when it was extracted, the predicted TSS may be right. Many transcription factor binding sites were predicted, such as Oct-1, AP-1, SP1, C/EBP,. (Fig.4). The regulation regions of the Pogene may be conserved as it has been reported that there are some AP-1 binding sites in enhancer of the genes (Sumitomo, 2012). Many HSF and several SRY binding sites could be found in the promoter region. The existence of these HSF binding sites implied that the gene express level of Pocould be regulated by the temperature shift. In our previous study we also found many HSF and SRY binding sites in the promoter region of heat shock protein (HSP) 70KDa (Qi., 2014). Although the genetic sex determination of the Japanese flounder was reported to be XX/XY type, it could easily be influenced by high temperature (Yama- moto, 1999). Temperature shift might change the expression levels of CIRP and HSP through the HSF. The existence of several SRY binding sites in the promoter region might indicate that these genes may also influence/or be associated with gonad development and sex differentiation.

    Fig.3 The genomic structure of the PoCIRP gene. The exons are shown as boxes, the introns are represented with black lines and the upstream region is represented with gray line. The CDS regions are shown with black boxes and the UTR regions are shown with white boxes.

    Fig.4 The predicted TSS, TATA box, CpGs island and transcription factor binding sites. The binding sites of different transcription factors are marked with different symbols. The vertical bars on the second line are used to represent the CpGs sites.

    3.4 The mRNA Expression Pattern of Poin Adult Tissues

    We performed the experiments to study the expression level of PomRNA in different flounder tissues, and the RT-qPCR was designed using the same amounts of RNAs from different tissues. The result revealed that PomRNA was expressed in all the tissues tested. The highest expression was detected in gonads, especially high expression was found in the ovary. The difference between testis and ovary was significant (<0.05). In all other tissues the Poexpression levels were not very high (Fig.5). Similar result has been reported in turbot, in which the predictedEST is specifically expressed in females (Taboada., 2012). Inwas reported to be found specifically in oocytes (Matsumoto., 2000). Our results, together with the previously re- ported ones, indicate that CIRP may be important for the physiological functions of ovary development or may have special functions in oogenesis and oocyte maturation.

    Fig.5 The quantitative analysis of PoCIRP expression in adult tissues. The relative expression variance is shown as ration (the amounts of PoCIRP were normalized to the corresponding 18S rRNA values).The data are shown as mean ± SD (n=3). Two way ANOVA analysis was performed to detect the significant differences between means. Columns with different letters show significant difference (P<0.05).

    The existence ofmRNA in many tissues of flounder such as brain and heart suggests that this gene might play many other roles in this fish. CIRP may play roles of protection for testis, brain and heart as reported in earlier papers (Liu, 2010; Xia., 2012).

    3.5 Methylation Analysis of the CpGs in the Pro-moter and Exon 1 Region

    The 229nt putative CpGs islands located in the region of exon I were predicted using EMBOSS Cpgplot. The relative positions between the predicted CpGs and the TF binding sites are showed in Fig.4. It happened that some of the CpG sites were adjacent to the TF binding site, such as C/EBP at 3567bp, EmBP-1 at 3584bp, myogenin at 3596bp and Sp1 at 3639bp (Fig.4). This indicates that the methylation of these sites may have influence on the expression of the CIRP gene. The promoter and exon I region’s GpGs methylation levels range from 0.0%-1.2% in different tissues of male and female (Fig.6).The detected methylation CpGs sites are located in the sites of 7th, 8th, 9th, 10th, 11st, 12nd, 14th, 15th CpG site. The detected methylation sites were related to the C/EBP, REB1, Sp-1 and EmBP-1. We compared the methylation level between the two factors: sexes and various tissues based on two way ANOVA analysis. The-value was 0.425 and the-value was 0.827 between male and female, which indicates that there is no significant difference between the sexes in methylation level. Meanwhile, the methylation level among different tissues also showed no significant difference as the-value was 0.048 and the-value was 0.828. It is known that the DNA methylation status can influence the gene expression level (Tate and Bird, 1993). So the results indicate that the Pomight be widely expressed in males’ and females’ different tissues and this has been proved in 3.4. And the expression level of the Pois not mainly regulated by the methylation status of the CpGs site.

    Fig.6 The methylation levels of the CpGs sites in the promoter region of both sex and different tissues. Ten clones were analyzed per fish. Result is showed as mean ± SE. The Y-axis value is the percentage of the methylated sites. Groups with different letters are significantly different (P<0.05). The methylation levels of all the groups show no significant difference.

    3.6 The Expression Pattern of PoDuring Developmental Stages

    The qPCR analysis showed a low level of maternal deposit ofmRNA in the fertilized eggs of Japanese flounder. The mRNA level was kept steady until multicell stages, began to increase at blastula stage and increased sharply at gastrula stage (Fig.7). The Poexpression reached the highest level at the neurula stage, and decreased afterward reaching relatively lower levels from tail-bud stage to hatching (Fig.7). One-way ANOVA analysis showed significant variations between different stages (<0.05).

    It has been reported that CIRP participates in many processes during the developmental stages (Peng, 2006; Tanaka, 2000; van Venrooy, 2008).was found to be expressed transiently in developing pronephros and neural tissues (Uochi and Asashima, 1998). However, there has been no direct experimental evidence showing the developmental regulation ofexpression in teleost fish. In the present study, the highest PomRNA level was observed at neurula stages. This is in agreement with the result obtained in, which may suggest its possible function in the development of the neural system in amphibians and teleosts. It is also noticeable that the embryonic expression, or the increase of expression, of Pobegan at blastula stage. This timing of increase is earlier than the differentiation of germ layers, but is the same as the expression of some genes specific for embryonic stem cells (ESC) and/or primordial germ cells (PGC) in Japanese flounder (Gao, 2013; Wu, 2014). Therefore, whether or not the start of Poexpression is associated with the formation of primitive cells for different germ layers, or with the appearance of ESC and PGC, will be of great interest.

    Fig.7 The quantitative analyses of PoCIRP expression during the developmental stages. The relative expression variance was shown as ration (the amounts of PoCIRP were normalized to the corresponding 18S rRNA values). The data was shown as mean±SD (n=3). One way ANOVA analysis was performed to detect the significant differences between means. Columns with different letters show significant difference (P<0.05).

    3.7 3D Modeling of CIRP

    The 3D model of PoCIRP was predicted with the Swiss-Model Workspace (Fig.8A). The solution stucture of RRM domain (1×5sA) was selected as the template model (Fig.8B) and the sequence identity was 69.51% between PoCIRP and the model sequence. The modelled residue sites ranged from 2 to 83, containing the RNA recognition motif. For assessment part, QMEANscore4 was 0.78, showing high quality of the model, QMEAN Z-score was ?0.17. All the assessment indicated that the 3D structure of PoCIRP model was of a good quality. When superposed to the template structure, the model showed very good overlapping (Fig.8C). By comparing the two models, we found that the target protein had similar composition to the second structures, such as helices and beta-strands, and that the spatial structure and location of these compositions were nearly the same. This also indicates that CIRPs are highly conserved proteins. With PDBsum Generate tools, the potential RNA binding site and the molecular surface cleft of PoCIRP were successfully displayed. There was a cleft between the binding domains (Fig.8D). This predicted structure may explain that the PoCIRP (RNA binding protein) interacts with mRNAs though the cleft.

    Fig.8 A, Model of the PoCIRP predicted by swiss program; B, The 3D structure of the template sequence(1×5sA) in the RCSB database; C, A structure superposition of the PoCIRP and the template. The yellow band stands for the template sequence and the green band stands for the target sequence; D, The molecular surface and cleft of the PoCIRP predicted by PDBsum generate.

    4 Conclusion

    In this study, we identified the mRNA and genomic sequences of Pofromwhich is a highly conserved protein in teleosts. Many TF binding sites were predicted, especially the HSF binding sites and SRY binding sites in the enhncer and promoter region of the Pogene. This could explain two phenomena. The first is that the expression level of thegoes down when temperature rises (Nishiyama., 1998a) and the second is that the temperature has influence on the sex ratio of Japanese flounder. In addition, a putative CpGs island was predicted and low methylation level was identified. The qPCR results verified a wide expression in all tissues with the highest expression level existing in ovary. Expression profiling during development stages indicated that the CIRP might be involved in the development of the neural system. The 3D structure prediction also gave information about the mechanism of the interaction between CIRP and mRNAs. The present study provides the basic information about the Poof the Japanese flounder and will benefit further studies on the function of this gene.

    Acknowledgements

    The work was supported by the National High Technology R&D Program of China (2012AA10A402) and the National Natural Science Foundation of China (31172385).

    Bergeron, D., Beauseigle, D., and Bellemare, G., 1993. Sequence and expression of a gene encoding a protein with RNA-binding and glycine-rich domains in., 1216 (1): 123-125.

    De Leeuw, F., Zhang, T., Wauquier, C., Huez, G., Kruys, V., and Gueydan, C., 2007. The cold-inducible RNA-binding protein migrates from the nucleus to cytoplasmic stress granules by a methylation-dependent mechanism and acts as a translational repressor., 313 (20): 4130-4144, DOI: 10.1016/j.yexcr.2007.09.017.

    Gao, J., Wang, J., Jiang, J., Fan, L., Wang, W., Liu, J., Zhang, Q., and Wang, X., 2013. Identification and characterization of a nanog homolog in Japanese flounder ()., 531 (2): 411-421, DOI: 10.1016/j.gene.2013.08.030.

    Gu, Y., Cao, D. C., Zhang, Y., Chang, Y., Li, O., Zhang, X., Lu, C., and Sun, X., 2009. Screening of EST-SSRs markers related to growth traits in common carp (L.)., (1): 15- 22 (in Chinese with English abstract).

    Hamid, A. A., Mandai, M., Fujita, J., Nanbu, K., Kariya, M., Kusakari, T., Fukuhara, K., and Fujii, S., 2003. Expression of cold-inducible RNA-binding protein in the normal endo- metrium, endometrial hyperplasia, and endometrial carcino- ma., 22 (3): 240-247, DOI: 10.1097/01.pgp.0000070851.25718.ec.

    Jin, F., Pan, Y., Li, S., Yang, H., Ji, H., Zhao, Q., and Yin, W., 2009. Cloning and sequence analysis of cold inducible RNA-binding protein cDNA from testis tissue in BALB/C Mice., 15 (1): 87-90, DOI: 10.3724/SP.J.1145.2009.00087.

    Larkin, M. A., Blackshields, G., Brown, N. P., Chenna, R., McGettigan, P. A., McWilliam, H., Valentin, F., Wallace, I. M., Wilm, A., Lopez, R., Thompson, J. D., Gibson, T. J., and Higgins, D. G., 2007. Clustal W and clustal X version 2.0., 23 (21): 2947-2948, DOI: 10.1093/bioinfor- matics/btm404.

    Li, S. Z., Jin, F. H., and Pang, Y., 2009. Cloning and Sequence Analysis of Cold Inducible RNA-binding Protein cDNA from Testis Tissue in BALB/C Mice., 15 (1): 87-90.

    Liu, A. J., Zhang, Z. W., Li, A. M., and Xue, J. H., 2010. Effects of hypothermia and cerebral ischemia on cold- inducible RNA-binding protein mRNA expression in rat brain., 13 (47): 104-110, DOI: 10.1016/j.brainres. 2010.05.029.

    Liu, Y., Hu, W., Murakawa, Y., Yin, J., Wang, G., Landthaler, M., and Yan, J., 2013. Cold-induced RNA-binding proteins regulate circadian gene expression by controlling alternative polyadenylation., (3): 2054, DOI: 10.1038/ srep02054.

    Matsumoto, K., Aoki, K., Dohmae, N., Takio, K., and Tsujimoto, M., 2000. CIRP2, a major cytoplasmic RNA-binding protein inoocytes., 28 (23): 4689- 4697, DOI: 10.1093/nar/28.23.4689.

    Morf, J., Rey, G., Schneider, K., Stratmann, M., Fujita, J., Naef, F., and Schibler, U., 2012. Cold-inducible RNA-binding protein modulates circadian gene expression posttran- scrip- tionally., 338 (6105): 379-383, DOI: 10.1126/ science.1217726.

    Nishiyama, H., Danno, S., Kaneko, Y., Itoh, K., Yokoi, H., Fukumoto, M., Okuno, H., Millan, J. L., Matsuda, T., Yoshida, O., and Fujita, J., 1998a. Decreased expression of cold-inducible RNA-binding protein (CIRP) in male germ cells at elevated temperature., 152 (1): 289-296.

    Nishiyama, H., Itoh, K., Kaneko, Y., Kishishita, M., Yoshida, O., and Fujita, J., 1997. A glycine-rich RNA-binding protein mediating cold-inducible suppression of mammalian cell growth., 137 (4): 899-908.

    Nishiyama, H., Xue, J. H., Sato, T., Fukuyama, H., Mizuno, N., Houtani, T., Sugimoto, T., and Fujita, J., 1998b. Diurnal change of the cold-inducible RNA-binding protein (CIRP) expression in mouse brain., 245 (2): 534-538, DOI: 10.1006/ bbrc.1998.8482.

    Pan, F., Zarate, J., Choudhury, A., Rupprecht, R., and Bradley, T. M., 2004. Osmotic stress of salmon stimulates upregulation of a cold inducible RNA binding protein (CIRP) similar to that of mammals and amphibians., 86 (7): 451-461, DOI: 10.1016/j.biochi.2004.06.006.

    Peng, Y., Yang, P. H., Tanner, J. A., Huang, J. D., Li, M., Lee, H. F., Xu, R. H., Kung, H. F., and Lin, M. C. M., 2006. Cold-inducible RNA binding protein is required for the expression of adhesion molecules and embryonic cell movement in., 344 (1): 416-424, DOI: 10.1016/ j.bbrc.2006.03.086.

    Qi, J., Liu, X., Liu, J., Yu, H., Wang, W., Wang, Z., Zhang, Q., 2014. Molecular characterization of heat shock protein 70 (HSP 70) promoter in Japanese flounder (), and the association of Pohsp70 SNPs with heat-resistant trait., 39 (2): 503-511, DOI: 10.1016/j.fsi.2014.05.038.

    Ripperger, J. A., Shearman, L. P., Reppert, S. M., and Schibler, U., 2000. CLOCK, an essential pacemaker component, controls expression of the circadian transcription factor DBP., 14 (6): 679-689, DOI: 10.1101/gad. 14.6.679.

    Saito, T., Sugimoto, K., Adachi, Y., Wu, Q., and Mori, K. J., 2000. Cloning and characterization of amphibian cold inducible RNA-binding protein.-, 125 (2): 237-245, DOI: 10.1016/s0305-0491(99)00174-1.

    Sugimoto, K., and Jiang, H., 2008. Cold stress and light signals induce the expression of cold-inducible RNA binding protein (CIRP) in the brain and eye of the Japanese treefrog ()..,, 151 (4): 628-636, DOI: 10.1016/j.cbpa.2008.07.027.

    Sumitomo, Y., Higashitsuji, H., Liu, Y., Fujita, T., Sakurai, T., Candeias, M. M., Itoh, K., Chiba, T., and Fujita, J., 2012. Identification of a novel enhancer that binds Sp1 and contributes to induction of cold-inducible RNA-binding protein (CIRP) expression in mammalian cells., 12 (1): 72, DOI: 10.1186/1472-6750-12-72.

    Taboada, X., Robledo, D., Del Palacio, L., Rodeiro, A., Felip, A., Martinez, P., and Vinas, A., 2012. Comparative expression analysis in mature gonads, liver and brain of turbot () by cDNA-AFLPS., 492 (1): 250-261, DOI: 10.1016/j.gene.2011.10.020.

    Tamura, K., Peterson, D., Peterson, N., Stecher, G., Nei, M., and Kumar, S., 2011. MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods., 28 (10): 2731-2739, DOI: 10.1093/molbev/msr 121.

    Tanaka, K. J., Kawamura, H., Matsugu, H., and Nishikata, T., 2000. An ascidian glycine-rich RNA binding protein is not induced by temperature stress but is expressed under a genetic program during embryogenesis., 243 (1-2): 207- 214, DOI: 10.1016/s0378-1119(99)00567-3.

    Tate, P. H., and Bird, A. P., 1993. Effects of DNA methylation on DNA-binding proteins and gene expression., 3 (2): 226-231.

    Tong, G., Endersfelder, S., Rosenthal, L. M., Wollersheim, S., Sauer, I. M., Buhrer, C., Berger, F., and Schmitt, K. R. L., 2013. Effects of moderate and deep hypothermia on RNA- binding proteins RBM3 and CIRP expressions in murine hippocampal brain slices., 15 (04): 74-84, DOI: 10.1016/j.brainres.2013.01.041.

    Uochi, T., and Asashima, M., 1998.(homolog of cold-inducible RNA-binding protein) is expressed tran- siently in developing pronephros and neural tissue., 211 (2): 245-250, DOI: 10.1016/s0378-1119(98)00102-4.

    van Venrooy, S., Fichtner, D., Kunz, M., Wedlich, D., and Gradl, D., 2008. Cold-inducible RNA binding protein (CIRP), a novel XTcf-3 specific target gene regulates neural development in., (8): 77, DOI: 10.1186/1471-213x-8-77.

    Wellmann, S., Buhrer, C., Moderegger, E., Zelmer, A., Kirschner, R., Koehne, P., Fujita, J., and Seeger, K., 2004. Oxygen-regulated expression of the RNA-binding proteins RBM3 and CIRP by a HIF-1-independent mechanism., 117 (Pt 9): 1785-1794, DOI: 10.1242/jcs.01026.

    Wu, X., Wang, Z., Jiang, J., Gao, J., Wang, J., Zhou, X., and Zhang, Q., 2014. Cloning, expression promoter analysis of vasa gene in Japanese flounder ()..,, (167): 41-50, DOI: 10.1016/j. cbpb.2013.06.004.

    Xia, Z., Zheng, X., Zheng, H., Liu, X., Yang, Z., and Wang, X., 2012. Cold-inducible RNA-binding protein (CIRP) regulates target mRNA stabilization in the mouse testis., 586 (19): 3299-3308, DOI: 10.1016/j.febslet.2012.07.004.

    Xue, J. H., Nonoguchi, K., Fukumoto, M., Sato, T., Nishiyama, H., Higashitsuji, H., Itoh, K., and Fujita, J., 1999. Effects of ischemia and H2O2 on the cold stress protein CIRP expression in rat neuronal cells., 27 (11-12): 1238-1244, DOI: 10.1016/s0891-5849 (99)00158-6.

    Yamamoto, E., 1999. Studies on sex-manipulation and production of cloned populations in hirame,(Temminck et Schlegel)., 173 (1-4): 235-246, DOI: http://dx.doi.org/10.1016/S0044-8486(98) 00448-7.

    Yang, C. L., and Carrier, F., 2001. The UV-inducible RNA- binding protein A18 (A18 hnRNP) plays a protective role in the genotoxic stress response., 276 (50): 47277-47284, DOI: 10.1074/jbc.M 105396200.

    Zhong, Q., Zhang, Q., Wang, Z., Qi, J., Chen, Y., Li, S., Sun, Y., Li, C., and Lan, X., 2008. Expression profiling and validation of potential reference genes duringembryogenesis.(), 10 (3): 310-318, DOI: 10.1007/s10126-007-9064-7.

    Zhou, K. W., Zheng, X. M., Yang, Z. W., Zhang, L., and Chen, H. D., 2009. Overexpression of CIRP may reduce testicular damage induced by cryptorchidism., 32 (2): E103-E111.

    (Edited by Ji Dechun)

    DOI 10.1007/s11802-015-2622-0

    ISSN 1672-5182, 2015 14 (1): 161-170

    ? Ocean University of China, Science Press and Springer-Verlag Berlin Heidelberg 2015

    (March 11, 2014; revised May 19, 2014; accepted June 12, 2014)

    * Corresponding author. Tel: 0086-532-82031806 E-mail: qzhang@ouc.edu.cn

    国产国拍精品亚洲av在线观看| 日本一二三区视频观看| 男人舔奶头视频| 亚洲av免费高清在线观看| 一区二区三区四区激情视频 | 黑人高潮一二区| av中文乱码字幕在线| 欧美激情国产日韩精品一区| 黄色欧美视频在线观看| 一级黄片播放器| 午夜爱爱视频在线播放| 国产精品爽爽va在线观看网站| 久久久久久伊人网av| 日韩欧美精品v在线| av天堂在线播放| 国内少妇人妻偷人精品xxx网站| 人妻夜夜爽99麻豆av| 亚洲人成网站在线播放欧美日韩| 99国产极品粉嫩在线观看| 国产一区二区在线观看日韩| 亚洲成人久久性| 两个人视频免费观看高清| 老师上课跳d突然被开到最大视频| 日韩欧美在线乱码| 免费观看在线日韩| av国产免费在线观看| 久久99热这里只有精品18| 亚洲av免费在线观看| 欧美精品国产亚洲| 久久中文看片网| 有码 亚洲区| 热99在线观看视频| eeuss影院久久| 久久精品夜夜夜夜夜久久蜜豆| 久久99热这里只有精品18| 日韩中字成人| 亚洲中文字幕一区二区三区有码在线看| 日本黄色视频三级网站网址| 国产精品1区2区在线观看.| 欧美不卡视频在线免费观看| 美女被艹到高潮喷水动态| 亚洲欧美成人精品一区二区| 小蜜桃在线观看免费完整版高清| 最近最新中文字幕大全电影3| 免费黄网站久久成人精品| 欧美成人精品欧美一级黄| 老司机影院成人| 99riav亚洲国产免费| 春色校园在线视频观看| 一本久久中文字幕| 久久午夜亚洲精品久久| 俺也久久电影网| 美女cb高潮喷水在线观看| 亚洲欧美清纯卡通| 在线免费十八禁| 人妻少妇偷人精品九色| 别揉我奶头 嗯啊视频| 久久久久免费精品人妻一区二区| 国产精品亚洲一级av第二区| 精品人妻视频免费看| 精品久久久久久久久av| 亚洲无线在线观看| 高清毛片免费观看视频网站| 国产综合懂色| 亚洲国产精品合色在线| 搡老岳熟女国产| 国产黄色小视频在线观看| 极品教师在线视频| 成人高潮视频无遮挡免费网站| 成人三级黄色视频| 国产精品三级大全| 97热精品久久久久久| 国产日本99.免费观看| 天堂网av新在线| 国内精品一区二区在线观看| 中文字幕av成人在线电影| 色吧在线观看| 变态另类丝袜制服| 男插女下体视频免费在线播放| 精品午夜福利视频在线观看一区| 中文字幕av在线有码专区| 国产午夜福利久久久久久| 成人国产麻豆网| 欧美一区二区亚洲| 熟女人妻精品中文字幕| 又黄又爽又刺激的免费视频.| 亚洲av五月六月丁香网| 村上凉子中文字幕在线| 露出奶头的视频| 老司机影院成人| 国产综合懂色| 国产精品1区2区在线观看.| 97人妻精品一区二区三区麻豆| 变态另类丝袜制服| 99久久精品国产国产毛片| 亚洲久久久久久中文字幕| 欧美绝顶高潮抽搐喷水| 久久久午夜欧美精品| 国产男人的电影天堂91| 18禁裸乳无遮挡免费网站照片| 欧美成人免费av一区二区三区| 国产精品不卡视频一区二区| 嫩草影视91久久| 国产欧美日韩精品亚洲av| videossex国产| 色av中文字幕| 欧美+日韩+精品| 久久欧美精品欧美久久欧美| 亚洲熟妇熟女久久| 国产亚洲av嫩草精品影院| 91精品国产九色| 亚洲三级黄色毛片| 老司机午夜福利在线观看视频| 日韩精品青青久久久久久| 少妇高潮的动态图| 不卡视频在线观看欧美| 看十八女毛片水多多多| 搡老熟女国产l中国老女人| 两性午夜刺激爽爽歪歪视频在线观看| 黄色日韩在线| 少妇熟女aⅴ在线视频| 小说图片视频综合网站| av在线蜜桃| 久久久午夜欧美精品| 亚洲av中文av极速乱| 午夜福利在线在线| 久久精品国产亚洲av香蕉五月| 亚洲av美国av| 三级毛片av免费| 亚洲av成人av| 久久久久久大精品| 在线看三级毛片| 欧美色欧美亚洲另类二区| 直男gayav资源| 观看免费一级毛片| 午夜视频国产福利| 亚洲,欧美,日韩| h日本视频在线播放| av.在线天堂| 国内少妇人妻偷人精品xxx网站| 日本与韩国留学比较| 丰满人妻一区二区三区视频av| 激情 狠狠 欧美| 一级黄片播放器| 在线国产一区二区在线| 欧美一区二区精品小视频在线| 国产国拍精品亚洲av在线观看| 别揉我奶头~嗯~啊~动态视频| 黄色欧美视频在线观看| 日韩,欧美,国产一区二区三区 | 深爱激情五月婷婷| 在线看三级毛片| 国产精品人妻久久久久久| 男女啪啪激烈高潮av片| 久久精品综合一区二区三区| 最近的中文字幕免费完整| 精品久久久久久久人妻蜜臀av| 成人亚洲欧美一区二区av| 精品国产三级普通话版| 嫩草影院入口| 亚洲美女搞黄在线观看 | 俄罗斯特黄特色一大片| 国产麻豆成人av免费视频| 青春草视频在线免费观看| 国产精品1区2区在线观看.| 热99在线观看视频| 欧美成人一区二区免费高清观看| 高清毛片免费看| 人妻制服诱惑在线中文字幕| 精华霜和精华液先用哪个| 成人漫画全彩无遮挡| 亚洲欧美日韩高清专用| 亚洲成a人片在线一区二区| 国产精品av视频在线免费观看| 国产精品乱码一区二三区的特点| 国产精品1区2区在线观看.| 成人av一区二区三区在线看| 少妇人妻一区二区三区视频| 无遮挡黄片免费观看| 国产精品一及| 免费看美女性在线毛片视频| 成人无遮挡网站| 国产精品美女特级片免费视频播放器| 久久草成人影院| 大又大粗又爽又黄少妇毛片口| 国产探花极品一区二区| 插逼视频在线观看| 日韩人妻高清精品专区| 欧美日韩精品成人综合77777| av在线天堂中文字幕| 欧美高清性xxxxhd video| 亚洲一区二区三区色噜噜| 成人漫画全彩无遮挡| av黄色大香蕉| 男人舔奶头视频| 性插视频无遮挡在线免费观看| 国产色婷婷99| 亚洲成人中文字幕在线播放| 日本欧美国产在线视频| 国产成人freesex在线 | 国产久久久一区二区三区| 悠悠久久av| 亚洲av第一区精品v没综合| 国产aⅴ精品一区二区三区波| 欧美+亚洲+日韩+国产| 男女视频在线观看网站免费| 免费一级毛片在线播放高清视频| 国产不卡一卡二| 精品欧美国产一区二区三| 久久久久久久久久成人| 亚洲av美国av| 老司机午夜福利在线观看视频| 亚洲中文字幕一区二区三区有码在线看| 综合色丁香网| 久久久a久久爽久久v久久| 亚洲自拍偷在线| 午夜精品在线福利| 天天一区二区日本电影三级| 最近手机中文字幕大全| 亚洲人成网站在线观看播放| 婷婷六月久久综合丁香| 韩国av在线不卡| 在线观看免费视频日本深夜| 亚洲最大成人av| 亚洲天堂国产精品一区在线| 2021天堂中文幕一二区在线观| 国产成人影院久久av| 成人一区二区视频在线观看| 大型黄色视频在线免费观看| 热99在线观看视频| 成人高潮视频无遮挡免费网站| 久久热精品热| 搡女人真爽免费视频火全软件 | 精品一区二区三区视频在线| 可以在线观看毛片的网站| 国产高潮美女av| 国产精品日韩av在线免费观看| 在现免费观看毛片| 国产色婷婷99| 亚洲中文字幕一区二区三区有码在线看| 欧美精品国产亚洲| 日韩 亚洲 欧美在线| 国产亚洲精品综合一区在线观看| 国产三级在线视频| 丝袜美腿在线中文| 99久国产av精品| 亚洲中文日韩欧美视频| 色在线成人网| 一个人观看的视频www高清免费观看| 国产成人a∨麻豆精品| 久久久久久伊人网av| 午夜福利在线观看吧| 真实男女啪啪啪动态图| 深夜a级毛片| 国产亚洲91精品色在线| 午夜激情福利司机影院| 久久久精品大字幕| 免费在线观看影片大全网站| 国产成年人精品一区二区| 久久午夜福利片| 国产在视频线在精品| 高清午夜精品一区二区三区 | 亚洲美女视频黄频| 在线观看午夜福利视频| 国产aⅴ精品一区二区三区波| 精品一区二区三区视频在线观看免费| 久久人人爽人人爽人人片va| 亚洲成人久久爱视频| 啦啦啦韩国在线观看视频| 男女啪啪激烈高潮av片| 国产精品女同一区二区软件| 欧美成人免费av一区二区三区| 亚洲经典国产精华液单| 免费av毛片视频| 夜夜夜夜夜久久久久| 成人av一区二区三区在线看| 最后的刺客免费高清国语| 亚洲欧美日韩东京热| 丰满乱子伦码专区| 久久久精品大字幕| 免费av毛片视频| 欧美成人免费av一区二区三区| 波多野结衣高清无吗| 男人和女人高潮做爰伦理| 亚洲色图av天堂| 偷拍熟女少妇极品色| 日日摸夜夜添夜夜添小说| 亚洲一级一片aⅴ在线观看| 女的被弄到高潮叫床怎么办| 午夜福利高清视频| 免费大片18禁| 精品午夜福利在线看| 五月玫瑰六月丁香| 久久久精品欧美日韩精品| 久久久国产成人精品二区| 成人综合一区亚洲| 国产亚洲精品av在线| 男女做爰动态图高潮gif福利片| 在线免费十八禁| 国产亚洲精品av在线| 人人妻人人澡欧美一区二区| 亚洲最大成人av| 国产亚洲精品久久久com| 国产精品永久免费网站| 寂寞人妻少妇视频99o| 亚洲精品色激情综合| 日韩高清综合在线| 亚洲av中文字字幕乱码综合| 成人三级黄色视频| 中国国产av一级| 精品人妻一区二区三区麻豆 | 在线观看免费视频日本深夜| 久久久精品94久久精品| 欧美成人免费av一区二区三区| 日本熟妇午夜| av国产免费在线观看| 99久国产av精品国产电影| 亚洲性久久影院| 少妇被粗大猛烈的视频| 亚洲国产色片| 搡老熟女国产l中国老女人| 亚洲成人精品中文字幕电影| 午夜福利在线观看吧| 免费观看精品视频网站| 国产午夜福利久久久久久| 亚洲精品国产成人久久av| 欧美成人精品欧美一级黄| 99视频精品全部免费 在线| 国产成人福利小说| 特大巨黑吊av在线直播| 寂寞人妻少妇视频99o| 少妇裸体淫交视频免费看高清| 欧美另类亚洲清纯唯美| 久久精品91蜜桃| 国产欧美日韩精品亚洲av| 国产69精品久久久久777片| av天堂中文字幕网| АⅤ资源中文在线天堂| 国产单亲对白刺激| 欧美不卡视频在线免费观看| 亚洲av美国av| a级毛色黄片| 国产老妇女一区| 深爱激情五月婷婷| 成人综合一区亚洲| 国产精品一区二区三区四区久久| 又黄又爽又刺激的免费视频.| 高清日韩中文字幕在线| 国产在视频线在精品| 菩萨蛮人人尽说江南好唐韦庄 | 一级毛片久久久久久久久女| 午夜日韩欧美国产| 女生性感内裤真人,穿戴方法视频| 午夜精品一区二区三区免费看| 在线免费观看的www视频| 国产综合懂色| 亚洲av美国av| 精品久久久噜噜| 天堂动漫精品| 亚洲在线观看片| 亚洲国产欧洲综合997久久,| 成人综合一区亚洲| 黑人高潮一二区| 亚洲av成人精品一区久久| 国产精品一及| 十八禁国产超污无遮挡网站| 国产精品一及| 男女那种视频在线观看| 五月玫瑰六月丁香| 啦啦啦观看免费观看视频高清| 国产精品人妻久久久影院| 精品久久久噜噜| 久久精品影院6| 亚洲丝袜综合中文字幕| 亚洲精品粉嫩美女一区| 免费人成在线观看视频色| 哪里可以看免费的av片| 国产精品电影一区二区三区| 精品久久国产蜜桃| 午夜免费男女啪啪视频观看 | av天堂中文字幕网| 日韩高清综合在线| 久久精品91蜜桃| 男人和女人高潮做爰伦理| 免费在线观看影片大全网站| 在线免费十八禁| 日本免费一区二区三区高清不卡| 亚洲成人久久爱视频| 日韩人妻高清精品专区| 韩国av在线不卡| 国产三级在线视频| 国产精品久久久久久久久免| 婷婷色综合大香蕉| 老司机午夜福利在线观看视频| 搡老岳熟女国产| 老司机福利观看| 久久亚洲精品不卡| 国产大屁股一区二区在线视频| 国产精品免费一区二区三区在线| 亚洲在线自拍视频| 亚洲av成人av| 看黄色毛片网站| 97在线视频观看| 日韩欧美在线乱码| 天美传媒精品一区二区| 国产毛片a区久久久久| 亚洲高清免费不卡视频| 亚洲国产高清在线一区二区三| 男人和女人高潮做爰伦理| 日本成人三级电影网站| 亚洲最大成人手机在线| 99久久久亚洲精品蜜臀av| 国产麻豆成人av免费视频| 亚洲成人久久性| 插阴视频在线观看视频| 久久久精品大字幕| 一卡2卡三卡四卡精品乱码亚洲| 日韩欧美 国产精品| 性插视频无遮挡在线免费观看| 亚洲图色成人| 午夜福利18| 国产片特级美女逼逼视频| 天天一区二区日本电影三级| 成人特级av手机在线观看| 内地一区二区视频在线| 波多野结衣高清无吗| 亚洲av免费高清在线观看| 波野结衣二区三区在线| 亚洲av一区综合| 成人精品一区二区免费| 国产成年人精品一区二区| 在线播放国产精品三级| eeuss影院久久| 久久99热6这里只有精品| 国产高清视频在线播放一区| 久久久久九九精品影院| 免费无遮挡裸体视频| 18禁在线播放成人免费| 日日摸夜夜添夜夜添小说| 深夜精品福利| 国产av一区在线观看免费| 97人妻精品一区二区三区麻豆| 午夜激情欧美在线| 国产高清视频在线播放一区| 午夜免费激情av| 欧美日韩一区二区视频在线观看视频在线 | 2021天堂中文幕一二区在线观| 人妻久久中文字幕网| www.色视频.com| 欧美一级a爱片免费观看看| 久久人人爽人人爽人人片va| 久久亚洲国产成人精品v| 日本a在线网址| 亚洲在线自拍视频| av女优亚洲男人天堂| 成人亚洲欧美一区二区av| 亚洲av熟女| 国产伦一二天堂av在线观看| 亚洲熟妇熟女久久| 一本精品99久久精品77| 久久亚洲精品不卡| 老司机福利观看| 五月伊人婷婷丁香| 国产成年人精品一区二区| 国内精品宾馆在线| 99视频精品全部免费 在线| 91久久精品国产一区二区成人| 亚洲精华国产精华液的使用体验 | av国产免费在线观看| 国产精品一二三区在线看| 日韩欧美 国产精品| 日韩高清综合在线| 国产精品国产高清国产av| 综合色丁香网| 日日摸夜夜添夜夜添av毛片| 麻豆国产97在线/欧美| 国产精品一区二区免费欧美| 99热网站在线观看| 内地一区二区视频在线| 男人舔女人下体高潮全视频| 欧美成人精品欧美一级黄| 村上凉子中文字幕在线| 99热这里只有是精品在线观看| 97碰自拍视频| 欧美xxxx性猛交bbbb| 成人毛片a级毛片在线播放| 久久九九热精品免费| 高清午夜精品一区二区三区 | 熟女人妻精品中文字幕| 我的老师免费观看完整版| 久久午夜亚洲精品久久| 日本-黄色视频高清免费观看| 国模一区二区三区四区视频| 欧美最新免费一区二区三区| 18禁在线无遮挡免费观看视频 | 晚上一个人看的免费电影| 在线免费观看的www视频| 国产亚洲精品综合一区在线观看| 精品久久久久久久末码| 久久热精品热| 九九爱精品视频在线观看| 天堂√8在线中文| 特级一级黄色大片| 国产亚洲精品久久久com| 精品一区二区三区视频在线观看免费| 99久久精品国产国产毛片| 在现免费观看毛片| 在线观看免费视频日本深夜| 国产亚洲av嫩草精品影院| 女同久久另类99精品国产91| 国产黄色小视频在线观看| 久久久久免费精品人妻一区二区| 亚洲国产欧洲综合997久久,| 可以在线观看毛片的网站| 欧美性猛交╳xxx乱大交人| 91在线观看av| 九九热线精品视视频播放| 午夜视频国产福利| 成人特级av手机在线观看| 中文字幕av在线有码专区| 亚洲精品国产成人久久av| 国产亚洲精品综合一区在线观看| 狂野欧美激情性xxxx在线观看| 亚洲精华国产精华液的使用体验 | 99在线人妻在线中文字幕| 欧美潮喷喷水| 国国产精品蜜臀av免费| 中文字幕av在线有码专区| 最近视频中文字幕2019在线8| av在线亚洲专区| 69人妻影院| 国产白丝娇喘喷水9色精品| 亚洲高清免费不卡视频| 在线天堂最新版资源| 久久九九热精品免费| 成人av一区二区三区在线看| 国产美女午夜福利| av天堂在线播放| 久久精品国产99精品国产亚洲性色| 色噜噜av男人的天堂激情| 免费人成在线观看视频色| 69av精品久久久久久| 国产精品一区二区性色av| 人妻久久中文字幕网| 国产黄a三级三级三级人| 亚洲欧美日韩东京热| 亚洲国产精品成人久久小说 | 嫩草影院入口| 99国产极品粉嫩在线观看| 乱码一卡2卡4卡精品| 我要看日韩黄色一级片| av.在线天堂| 老司机午夜福利在线观看视频| 国产成年人精品一区二区| 免费看光身美女| 别揉我奶头~嗯~啊~动态视频| 最新中文字幕久久久久| 人人妻人人看人人澡| 中国美女看黄片| 欧美又色又爽又黄视频| 国产高清三级在线| 国产高清有码在线观看视频| 天天躁夜夜躁狠狠久久av| 在线播放国产精品三级| 成人漫画全彩无遮挡| 永久网站在线| 精品久久久久久成人av| 狂野欧美激情性xxxx在线观看| 熟女电影av网| 成年女人毛片免费观看观看9| 波多野结衣高清作品| 国产精品国产高清国产av| 51国产日韩欧美| 一进一出好大好爽视频| 中文资源天堂在线| 毛片女人毛片| 少妇丰满av| 麻豆一二三区av精品| 一级毛片我不卡| 91av网一区二区| 精品国内亚洲2022精品成人| 亚洲aⅴ乱码一区二区在线播放| 精品人妻偷拍中文字幕| av黄色大香蕉| 搡老妇女老女人老熟妇| 中文字幕av成人在线电影| 直男gayav资源| 免费大片18禁| 久久天躁狠狠躁夜夜2o2o| 婷婷精品国产亚洲av在线| 99热这里只有是精品50| 日韩av在线大香蕉| 日产精品乱码卡一卡2卡三| 国产高清三级在线| 丝袜美腿在线中文| av中文乱码字幕在线| 人妻丰满熟妇av一区二区三区| 国产精品久久久久久av不卡| 免费看a级黄色片| 大又大粗又爽又黄少妇毛片口| 亚洲高清免费不卡视频| 亚洲美女黄片视频| 成人欧美大片| 国产成人影院久久av| 五月伊人婷婷丁香| 午夜日韩欧美国产| 男人和女人高潮做爰伦理| 草草在线视频免费看| 国产伦在线观看视频一区| 午夜爱爱视频在线播放| 我的女老师完整版在线观看| 内射极品少妇av片p| 久久精品夜色国产| 深夜精品福利| 亚洲精品一卡2卡三卡4卡5卡|