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

    Characterization and Expression Analysis of a Complement Component Gene in Sea Cucumber (Apostichopus japonicus)

    2015-04-01 02:11:59CHENZhongZHOUZunchunYANGAifuDONGYingGUANXiaoyanJIANGBeiandWANGBai
    Journal of Ocean University of China 2015年6期

    CHEN Zhong, ZHOU Zunchun, YANG Aifu, DONG Ying, GUAN Xiaoyan,JIANG Bei, and WANG Bai

    ?

    Characterization and Expression Analysis of a Complement Component Gene in Sea Cucumber ()

    CHEN Zhong, ZHOU Zunchun*, YANG Aifu, DONG Ying, GUAN Xiaoyan,JIANG Bei, and WANG Bai

    ,,116023,

    The complement system plays a crucial role in the innate immune system of animals. It can be activated by distinct yet overlapping classical, alternative and lectin pathways.In the alternative pathway, complement factor B (Bf) serves as the catalytic subunit of complement component 3 (C3) convertase, which plays the central role among three activation pathways. In this study, the Bf gene in sea cucumber (), termed, was obtained by rapid amplification of cDNA ends (RACE). The full-length cDNA ofwas 3231bp in length barring the poly (A) tail. It contained an open reading frame (ORF) of 2742bp encoding 913 amino acids, a 105bp 5’-UTR (5’-terminal untranslated region) and a 384bp 3’-UTR. AjBf was a mosaic protein with six CCP (complement control protein) domains, a VWA (von Willebrand factor A) domain, and a serine protease domain. The deduced molecular weight of AjBf protein was 101kDa. Quantitative real time PCR (qRT-PCR) analysis indicated that the expression level ofinwas obviously higher at larval stage than that at embryonic stage. Expression detection in different tissues showed thatexpressed higher in coelomocytes than in other four tissues. In addation,expression in different tissues was induced significantly after LPS or PolyI:C challenge. These results indicated that AjBf plays an important role in immune responses to pathogen infection.

    sea cucumber;; complement factor B; LPS challenge; PolyI:C challenge

    1 Introduction

    The complement system is a highly sophisticated biological reaction system, which plays a crucial role in body defense as a part of both the innate and the adaptive immune systems (Nonaka and Yoshizaki, 2004). This system can be activated by a series of limited proteolytic cascades through three parallel pathways. The classical pathway (CP) is activated by antibody-antigen (Ab-Ag) complexes and is a major effector of the acquired immune system that arose in the jawed vertebrate lineage. The other two, including the alternative pathway (AP) and the lectin pathway (LP), function in innate immune defense, and are considered to be important in invertebrates defending against microbial infections (Nonaka, 2014; Sekine., 2001; Fujita., 2004). These three pathways merge at the proteolytic activation step of the complement component 3 (C3), and the activation fragments of C3 induce several effecter mechanisms against the pathogens and facilitate the immune response (.). In mammalians, complement factor B (Bf) serves as the catalytic subunit of C3 convertase in the alternative pathway (Pangburn, 1983) while in the classical pathway, this function is subjected to C2 (Ishii., 1993; Wei., 2009). In addition, mammalian Bf and C2 have identical modular structures, including three complement control protein (CCP) domains at the N- terminus, one von Willebrand factor A (VWA) domain and one serine protease (SP) domain at the C-terminus (Horiuchi., 1993; Tuckwell., 1997). The VWA and CCP structural modules contain binding sites for C3b and C4b, respectively (Tuckwell., 1997; Xu and Volanakis, 1997).

    The most ancestral pathway to activate C3 is mediated by factor B (Prado-Alvarez., 2009). This mechanism, presenting in sea anemone(.), carpet-shell clam(Prado-Alvarez., 2009), horseshoe craband(Zhu., 2005; Tagawa., 2012), sea urchin(Gross., 1999), ascidian(Yoshizaki., 2005; Nicola., 2014), lamprey(Nonaka, 1994) and amphioxus(He., 2008) resembles the alternative pathway from mammals and constitutes an ‘a(chǎn)rcheo-complement system’. These Bf proteins often have an atypical structure compared with mammalian Bf or C2. For example, sea urchin (Smith., 1998), horseshoe crab (Zhu., 2005) and sea anemone Bf proteins (.) all have five CCP domains, while tunicate Bf possesses four CCP domains with two additional low-density lipoprotein receptor (LDLR) domains. Some marine invertebrates contain more than one isoform of Bf gene. For example, two homologs of Bf gene were found in sea anemone(Kimura., 2009) and horseshoe crab(Tagawa., 2012), while three were found in sea urchin(Hibino., 2006), ascidian(Yoshizaki., 2005) and amphioxus(Huang., 2008), respectively. The presence of multiple Bf genes in marine species may be resulted from gene duplication and exon shuffling (Yoshizaki., 2005). However, the difference in functions among Bf isomers is still unclear. Despite the distinctive characteristics, the described invertebrate Bf proteins have been considered orthologues of the mammalian Bf and C2 (Nonaka and Kimura, 2006).

    Sea cucumber () belongs to echinoderms, which locates between invertebrates and chordates systematically, and occupies an important evolutionary position. It is a reference point for inferring the evolution of deuterostome immune systems. Sea cucumber is also an economically important aquaculture species. In previous studies, a Bf gene of, named-2 (JN- 634069.1) has been cloned and studied on its transcriptional expression after LPS challenge (Zhong., 2012). AjBf-2 shows similarity to the Bf genes of other invertebrates in molecular structure, and is expressed in all determined tissues and has been proved to be involved in the immune response against Gramnegative bacteria. Here, complement Bfgene in, named, was cloned and characterized. The expression patterns ofwere investigated in different tissues at embryonic and larval developmental stages. Furthermore, the gene expression profiles after LPS and PolyI:C challenge were also analyzed in the different tissues of. We also compared the structure and expression patterns ofwith those ofwhich was cloned by Zhong. (2012). Our study aimed to enrich the information of echinoderms complement components that are complementary to comparative and evolutionary aspects of immunity.

    2 Materials and Methods

    2.1 Full-Length cDNA Cloning of

    The partial sequences ofwere obtained from the ESTs ofby cDNA libraries sequencing (Yang., 2009; Zhou., 2014). To obtain the full-length cDNA of, rapid amplication of cDNA ends (RACE) PCR was carried out using a SMART RACE cDNA Amplification Kit (Clontech, Japan). The primer-F and-R were designed according to the known EST sequence of AjBfgene (Table 1). The primer pairs of-F,-R and UPM primer (supplied by Clontech) were used for 3’ RACE and 5’RACE, respectively. The PCR products were purified, cloned into pMD19-T vector (Takara, Japan) and subsequently sequenced. Finally, the full-length cDNA ofwas obtained by assembling the 3’ RACE and 5’RACE fragments.

    Table 1 Primers for RACE and qRT-PCR used in this study

    2.2 Characterization and Phylogenetic Analysis of AjBf

    For sequence analysis, the Bf amino acid sequences were either identified by simple key word searches, or with BLAST searches using sea urchin Bf amino acid sequences at NCBI (http://www.ncbi.nlm.nih.gov/BLAST). Protein sequences retrieved from public database were used for open reading frame (ORF) and domain searches, alignment, and phylogenetic reconstruction. ORF was predicted using Open Reading Frame Finder (http:// www.ncbi.nlm.nih.gov/gorf/gorf.html). The protein domains were predicted by the Simple Modular Architecture Reach Tool (SMART, http://smart.embl-heidelberg.de/). Translation and protein analysis were performed using ExPaSy tools (http://us.expasy.org/tools). Multiple alignment of the Bf was performed using the program ClustalX version 2.0 (Larkin., 2007). Phylogenetic tree was constructed using the neighbor-joining method based on the deduced full-length amino acid sequences within the Molecular Evolutionary Genetics Analysis (MEGA 4.0) package (Tamura., 2007). Data were analyzed using Poisson correction, and gaps were removed by complete deletion. The topological stability of the tree was evaluated by 1000 bootstrap replications.

    2.3 Sampling and Immune Challenge

    Mature sea cucumbers were collected from Guanglu Island (Dalian, China) in early July. Animals were artificially spawned by flowing sea water (20–21℃) stimulation. During the development from fertilized eggs to 1mm-long juvenile, sea cucumbers were cultured in filtered sea water (temperature: 20–21℃, salinity: 32, pH: 7.8). To investigate the expression of AjBfgene at different developmental stages, samples from every developmental stages including unfertilized eggs, fertilized eggs, cellulous stages, blastula prior to hatching, gastrula, early auricularia, auricularia, late auricularia, doliolaria, pentactula and 1-mm long juvenile, were collected by sieving using 60μm filter, pelleted by centrifugation (Labnet Spectrafuge, USA) and then stored in 1.5mL microcentrifuge tubes. All samples were frozen immediately with liquid nitrogen and then stored at ?80℃ prior to RNA isolation.

    To investigate the AjBfgene expression in different tissues, the intestine, respiratory tree, coelomocytes, body wall and tube feet tissues were separated from 15 healthy individuals (3 pools with 5 individuals each pool) (average body weight 10.2g). These tissues from different individuals were mixed respectively. Coelomocytes were collected according to the methods of Santiago-Cardona (Santiago-Cardona,2003). Intestine, respiratory tree, body wall and tube feet were excised, briefly rinsed in filtered seawater, and frozen immediately with liquid nitrogen and then stored at ?80℃.

    To investigate the expression patterns of AjBf gene after immune challenge, Two different immune challenges were conducted by coelomic injection for each sea cucumber with 500μL PolyI:C (100mgmL?1) or LPS (1mgmL?1). PolyI:C (Sigma, P9582) and LPS (Sigma, L2880) were diluted with phosphate buffered saline (PBS, pH= 7.4). The sea cucumbers injected with 500μL PBS were treated as the control group (Ramírez-Gómez, 2009). The sea cucumbers without any injection were treated as the blank group (0h). The different tissues were collected for quantitative real-time PCR (qRT-PCR) analysis after the sea cucumber were injected for 4h, 12h, 24h, 48h and 72h, respectively.At each sampling time, the tissues taken from 15 different individuals (3 pools with 5 individuals each pool) were mixed in three 1.5mL microcentrifuge tubes, and thrown immediately into liquid nitrogen and then stored at ?80℃.

    2.4 Quantitative Real-Time PCR

    Total RNA was isolated using the RNAprep pure Tissue Kit (Tiangen, Beijing, China) according to the manufacturer’s instructions. The quality and quantity of extracted total RNA were measured using the NanoPhotometer (Implen GmbH, Munich, Germany) and agarose gel electrophoresis. First strand cDNA synthesis was performed in a volume of 20μL with 900ng total RNA, 25pmol Oligo dT Primer, 50pmol random 6mers, 4μL 5×PrimeScriptTM buffer, 1μL PrimeScriptTM RT enzyme Mix I (PrimeScriptTM RT reagent Kit, TaKaRa, China). Reactions were conducted by incubating at 37℃ for 15min, and then at 85℃ for 5s to deactivate the enzyme.

    To know the AjBf gene expression profiles ofat different developmental stages, in different tissues and after LPS or PolyI:C challenge, equal amounts of cDNA template were used in qRT-PCR. The cytochrome b (Cytb) gene was used as the reference gene (Yang., 2009), and the relative expression of the AjBf gene was assessed using qRT-PCR in the Mx3000pTMdetection system (Applied Stratagene, USA). The reactions were performed in a total volume of 20μL containing 10μL of 2×SYBR Green Master mix (SYBR PrimeScript_RT-PCR Kit II, TaKaRa), 0.4μL of ROX Reference Dye II, 2μL of cDNA template, and 0.4μmol L?1of each primer (Table 1). The qRT-PCR program was 95℃ for 30s, followed by 40 cycles of 95℃ for 10s, 56℃ for 25s and 72℃ for 25s. Melting curve analysis of amplification products was performed at the end of each PCR reaction to confirm that only one PCR product was amplified and detected. In addition, the amplicons were checked by agarose gel with a 100bp ladder to confirm the correct amplicon sizes.

    2.5 Statistical Analysis

    All the statistical analyses about the expression level ofmRNA were normalized by Cytb gene. We used the Relative Expression Software Tool 384 v.2 (REST) (Pfaffl., 2002) to study themRNA expression profiles. The crossing-point (CP) values were converted to fold differences by relative quanti?cation method. The mathematical model used is based on the PCR efficiencies and the mean crossing point deviation between the control and treatment groups by pairwise fixed reallocation randomization test. The statistically significant difference was defined as<0.05.

    3 Results

    3.1 AjBf Gene Isolation and Sequencing

    Based on our cDNA libraries sequencing of(Yang., 2009), an EST, 1624bp in length, was obtained, which was similar to the partial sequence of sea urchinin Genbank. The full-length cDNA sequence of AjBf gene was obtained by RACE PCR and sequencing. The assembled sequence showed that the full-length cDNA sequence of AjBf gene was 3231bp barring the poly (A) tail. It contained a 105bp 5’-UTR (terminal untranslated region), an open reading frame (ORF) of 2739bp encoding 913 amino acids and a 384bp 3’-UTR with a polyadenylation signal (AATAAA motif) (Fig.1).

    AjBf had a deduced molecular weight of 101kDa. The deduced amino acid sequences of AjBf contained a possible signal peptide, suggesting that it was a secreted protein. Like the other Bf or C2 family members, AjBf was a mosaic protein including CCPs, a VWA domain, and a serine protease domain (Fig.2). It was interesting that AjBf had six CCPs, while both AjBf-2 (Zhong., 2012) and SpBf (sea urchin) had five, and most of vertebrate Bf or C2 proteins had three. In addition, one Arg- Phe (RF) bond, the putative cleavage site for C1s and/or factor D to activate Bf, was located at residues 448–449 before the VWA domain of AjBf (Fig.3). In the Serine Protease domain of AjBf, the critical catalytic triad was His 698, Asp742 and Ser851. There were eight consensus recognition sequences for N-linked glycosylation through- out the AjBf sequence (Fig.1). Four were located in the CCPs, one was found in the region between the CCPs and the VWA domain, one was in the VWA domain, and the other two were in the serine protease domain. A CCP domain was typically about 60 amino acids long with a number of conserved residues. These amino acids in- cluded four cysteines, three glycines, two prolines, two tyrosines (or phenylalanine) and one tryptophan, all of which were located in specific positions within the domain (Reid and Day, 1989). The aligned results of six CCP domains from AjBf were shown in Fig.4. Two disulfide bonds were formed between the four cysteines and maintain the topology of the domain. Although some of the consensus amino acids were missing from AjBf, all four cysteines were present in each CCP, suggesting that these domains fold as expected.

    Fig.1 The deduced protein sequence of AjBf. The full-length cDNA has been deposited in GenBank with accession number HQ993063. The signal peotide is underlined. The nucleotides of initiation codon and the polyadenylation signal are underlined too. The N-linked glycosylation sites are shaded.

    Fig.2 Structural feature of AjBf. Schematic representations of AjBf were predicted by SMART program. Signal peptide determined by the SignalP program is shown in red. Abbreviations indicate as follows: CCP, complement control protein; VWA, von Willebrand factor A; Tryp_SPc, serine protease.

    Fig.3 Alignment of the VWA and serine protease domains from AjBf and the other Bf/C2 family proteins. To simplify the analysis, the CCP domains were deleted, and the remaining domains were aligned with CLUSTAL W program. Domains and amino acids with functional significance are labeled. Accession numbers for sequences used in this alignment can be found in the legend for Fig.5. Aj: A. japonicus; Sp: Strongylocentrotus purpuratus; Hr: Halocynthia roretzi; Hs: Homo sapiens; Mm: Mus musculus; Ol: Oryzias latipes; Xl: Xenopus laevis; Nv: Nematostella vectensis; Lc: Lethenteron camtschaticum; Br: Branchiostoma belcheri tsingtauense.

    3.2 Phylogenetic Analysis

    Pairwise alignments were used to calculate percentage of amino acid similarities and identities between AjBf and all the other Bf porteins without the CCP domains. The AjBf amino acids showed 86.4% similarity with AjBf-2, and 19.2%–31.2% similarities with Bfs in some other species. Consequently, to assess the relationships among the Bf protein sequences, alignments of full length of the AjBf and other Bf/C2 proteins were used for phylogenetic analysis (Fig.5). The evolutionary tree analysis showed that the Bf genes were evolved from a common ancestor in the evolutionary development of two directions, one cluster was the invertebrates and cephalochordates, and another was the Urochordates and vertebrates. Sea cucumber had the closest relationship with sea urchin, and they were involved in one cluster with sea anemone and Branchiostomas.

    Fig.4 Alignment of the AjBf CCPs. The six CCPs from AjBf are aligned to each other with the consensus amino acids shown in bold and at the bottom. Two disulfide bonds maintain the CCP domain structure and are formed between the first and third cysteine, and the second and fourth cysteine. The alignment was done with CLUSTAL W.

    Fig.5 Phylogenetic tree. The tree was constructed using ClustalW and MEGA (3.1). The numbers of branches are bootstrap values for 1000 replicates. The GenBank accession numbers for the sequences are as follows: AjBf-2 (Apostichopus japonicus), AEP68015; SpBf (Strongylocentrotus purpuratus), NP_999700; HrBf (Halocynthia roretzi), AAK00631; HsBf (Homo sapiens), NP_001701; MmBf (Mus musculus), AAC05283; OlBf (Oryzias latipes), NP_001098275; XlBf (Xenopus laevis), BAA06179; NvBf (Nematostella vectensis), BAH22728; LcBf (Lethen- teron camtschaticum), BAA02763; BrBf (Branchiostoma belcheri tsingtauense), ABY28382; HsC2 (Homo sapiens), AAB97607; MmC2 (Mus musculus), AAA63294; BtC2 (Bos taurus), AAI03358; XlC2 (Xenopus laevis), ABB- 85337; OmC2 (Oncorhynchus mykiss), NP001117673.

    3.3 AjBf Gene Expression Analysis

    The qRT-PCR analysis was employed to determine specific expression patterns of AjBf gene at eleven developmental stages of sea cucumber. The expression level of AjBf gene in unfertilized eggs was higher than that in the multi-cell and blastula stages, and AjBfgene expression increased significantly from the gastrulas to the juveniles. The highest expression level was detected in the late auricularias, about 10-fold higher than that in unfertilized eggs (<0.05). The AjBf gene expression levels were obviously higher at larval stages than those at embryonic stages (Fig.6).

    Fig.6 The relative transcript abundance of AjBf gene at different developmental stages of A. japonicus determined through qRT-PCR analysis. 1, Unfertilized egg; 2, Fertilized egg; 3, Cellulous stage; 4, Blastula; 5, Gastrula; 6, Early auricularia; 7, Auricularia; 8, Late auricularia; 9, Doliolaria; 10, Pentactula; 11, Juvenile. Each symbol and vertical bar represents the mean±S.D (n=3). Asterisks represent a significant difference at the level of P<0.05 relative to unfertilized egg.

    The AjBf gene expression in different tissues including intestine, respiratory tree, coelomocytes, tube feet and body wall, were also detected. The results indicated thattranscripts were detected in all tested tissues, the lowest expression level was found in the intestine and thehighest level was found mainly in the coelomocytes, which was about 30-fold higher than that in intestine (Fig.7).

    Fig.7 The relative mRNA levels of AjBf gene in different tissues. I, intestine; BW, body wall; C, coelomocytes; F, tube feet; RT, respiratory tree. Each vertical bar represents the mean±S.D (n=3). Expression levels in all tissues are presented relatively to that in the intestine tissue (1×).

    AjBf gene expression was significantly induced by LPS and PolyI:C in the tested tissues. After LPS challenge, AjBf gene expression in coelomocytes increased at 12h, and peaked at 48h (about 9 fold), then decreased at 72h. It is noted that AjBf gene expression levels were significantly up-regulated at 4h and 72h in body wall. In respiratory tree, AjBf gene expression increased significantly at three time points except of 12h and 48h. While, significant up-regulations of AjBf gene were just observed at 4h and 24h in tube feet. In intestine, the expression levels of AjBf gene were relatively lower at 12h and 72h, and higher at 4h (Fig.8A).

    After PolyI:C challenge, AjBf gene expressed strongly in body wall, 9-fold higher than that of the control group at 4h. In coelomocytes, the expression level increased at 4h after challenge, and reached a peak at 24h (about 3.5 fold), then decreased gradually to the similar level as in the LPS group. AjBf gene expression was significantly up-regulated in intestine at 72h. In tube feet and respiratory tree, it was similar to the LPS challenge group (Fig.8B).

    Fig.8 Fold induction of AjBf gene challenged by LPS (A) and PolyI:C (B) in five tissues. Relative expression of AjBf gene was expressed as fold changes over control samples taken at the same time interval as normalized to change in expression in the Cytb control. Each vertical bar represents the mean±S.D (n=3). Asterisks represent a significant difference at the level of P<0.05 relative to appropriate control.

    4 Discussion

    In the present study, we isolated and characterized the full-length cDNA of the complement factor B-like gene in sea cucumber. Its expression profiles in normal tissues and after LPS and PolyI:C challenges were also analyzed. Complement component C3, a key specific protease in the alternative pathway, had been proved to be present in the sea cucumber (Zhou., 2011). The isolation of Bf and C3 provided a strong evidence for the presence of an alternative activation pathway in sea cucumber.

    Like other members of the Bf/C2 family, AjBf had a mosaic structure including CCPs, a von Willebrand factor A (VWA) domain, and a serine protease domain. The alignment between AjBf and vertebrate Bf/C2 proteins revealed amino acids in conserved positions for serine protease activity, a conserved factor D cleavage site, and Mg2+binding sites.

    We used to believe that five CCP domains is the ancestral structure of Bf/C2 protein family. The recently identified members of Bf/C2 from horseshoe crab and sea urchin have five CCP domains (Zhu., 2005; Smith., 1998), and sea anemone has two Bf homologue with three and five CCP domains, respectively (.), while the vertebrate Bf/C2 members all have three CCP domains, suggesting that vertebrate Bf/C2 may have lost two (Horiuchi., 1993; Hourcade., 1995). In this study, six CCPs were found in AjBf, which was more than AjBf-2 and Bf in other species. Compared with AjBf-2, AjBf has an extra CCP domain and a signal peptide. Six CCPs rather than five may more likely be the ancestral type. Considering the important functions of CCPs in Bf for binding to C3b in vertebrates, the three CCPs close to the VWA domain of AjBf might be more important for binding functions than the three N-terminal CCPs.

    Phylogenetic analysis showed that the Bf gene of sea cucumber has the closest relationship with that of sea urchin, and they are in one cluster with those of sea anemone and Branchiostomas. Based on the sequences alignment and phylogenetic trees, the complement Bf-like homology of sea cucumber is close to Bf rather than C2. In vertebrates, the Bf/C2 families continue to differentiate to Bf and C2. AjBf may be one of the most primitive members of the Bf/C2 protein family.

    The existing of AjBf and AjC3 in sea cucumber indicated that its simple complement pathway is similar to the complement system in amphioxus, which is also composed of C3 and Bf (He., 2008). This simple sea cucumber complement pathway might function like the alternative pathway in vertebrates (Tuckwell., 1997). Thetranscripts were persistent during the embryogenesis. The expression level oftranscripts in unfertilized eggs was significantly higher than those in the multi-cell and blastula stages. It was significantly increased from the gastrulas to the juveniles, and the highest level was detected in the late auricularias (Fig.6). This suggested that complement-mediated defense system may function in unfertilized eggs and early embryos. The increase ofexpression level in the gastrulas period probably meant that the embryos are protected by the fertilization envelope, and do not require a functional defense system. From the gastrulas, because of the organ formation and the change of life habits, thetranscriptional levels increased significantly.

    The gene encoding AjBf expressed specifically in coelomocytes than in other tissues. Coelomocytes have long been considered to be the essential immune effector cells of the immune response in echinoderm, so it is not surprising that AjBf gene showed intensive expression in the coelomocytes (Al-Sharif., 1998). For-2, high level expression was found in both coelomocytes and body wall, while the expression was only in a detectable level in other tissues (Zhong., 2012). The expressions ofand-2 were both up-regulated in coelomocytes after LPS challenge, and followed by a decrease. However, the time to peak value was different.expression reached the peak at 48h, while-2 reached the peak at 12h. In body wall, the expressions ofand-2 both were significantly up-regulated at 4h after LPS challenge. It was inferred that both AjBf and AjBf-2 are involved in the immune response of sea cucumber to LPS challenge, while their functions need to be further studied.

    LPS is a major component of cell walls of gram- negative bacteria and an extremely useful immune-activating substance. Many immune-related genes have been documented to respond to LPS challenge. PolyI:C, the synthetic analog of dsRNA, can be treated as a potent immunoregulatory agent for viral infection. Welsh found that human serum can deplete RNA virus through either the classical complement pathway component C4 or the alternative complement pathway component factor B (Welsh, 1977). Furthermore, it was reported that LPS and dsRNA could effectively increase phagocytes and phagocytic activities in coelomocytes of the sea cucumber(Ramírez-Gómez., 2010). In this study, we found that the AjBf gene expression increased significantly after LPS challenge, which was consistent with gene C3 as we reported early (Zhou., 2011). Thus both Bf and C3 may be involved in the immune response against the bacteria pathogen by alternative pathway. In PolyI:C group,AjBf expression was significantly induced in different tissues, implying that the sea cucumber responds to RNA virus by activating the alternative complement pathway.

    In sea cucumber, coelomocytes are the essential immune effector cells, and body wall is the first barrier against pathogen invasion. Except these two tissues, immune responses of other three tissues, including respiratory tree, intestine and tube feet, were simultaneously detected in our study. Positive results showed that they were all involved in immune responses in sea cucumber.

    Acknowledgements

    This work was supported by the National Natural Science Foundation of China (Nos. 31272687 and 30972 272), Science and Technology Project of Liaoning Province (2014203006), Ocean and Fisheries Project of Liao- ning Province (201301), and Science and Technology Project of Dalian City (2013J21DW025).

    Al-Sharif, W. Z., Sunyer, J. O., Lambris, J. D., and Smith, L. C., 1998. Sea urchin coelomocytes specifically express a homologue of the complement component C3., 160: 2983-2997.

    Franchi, N., and Ballarin, L., 2014. Preliminary characterization of complement in a colonial tunicate: C3, Bf and inhibition of C3 opsonic activity by compstatin., 46 (2): 430-438.

    Fujita, T., Matsushita, M., and Endo, Y., 2004. The lectin-com- plement pathway: Its role in innate immunity and evolution., 198: 185-202.

    Gross, P. S., Al-Sharif, W. Z., Clow, L. A., and Smith, L. C., 1999. Echinoderm immunity and the evolution of the complement system., 23: 429-442.

    He, Y., Tang, B., Zhang, S., Liu, Z., Zhao, B., and Chen, L., 2008. Molecular and immunochemical demonstration of a novel member of Bf/C2 homolog in amphioxus: Implications for involvement of hepatic cecum in acute phase response., 24: 768-778.

    Hibino, T., Loza-Coll, M., Messier, C., Majeske, A. J., Cohen, A. H., Terwilliger, D. P., Buckley, K. M., Brockton, V., Nair, S. V., Berney, K., Fugmann, S. D., Anderson, M. K., Pancer, Z., Cameron, R. A., Smith, L. C., and Rast, J. P., 2006. The immune gene repertoire encoded in the purple sea urchin genome., 300 (1): 349-365.

    Horiuchi, T., Kim, S., Matsumoto, M., Watanabe, I., Fujita, S., and Volanakis, J. E., 1993. Human complement factor B: cDNA cloning, nucleotide sequencing, phenotypic conversion by site-directed mutagenesis and expression., 30: 1587-1592.

    Hourcade, D. E., Wagner, L. M., and Oglesby, T. J., 1995. Analysis of the short consensus repeats of human complement factor B by site-directed mutagenesis.,270: 19716-19722.

    Huang, S., Yuan, S., Guo, L., Yu, Y., Li, J., Wu, T., Liu, T., Yang, M., Wu, K., Liu, H., Ge, J., Yu, Y., Huang, H., Dong, M., Yu, C., Chen S., and Xu, A., 2008. Genomic analysis of the immune gene repertoire of amphioxus reveals extraordinary innate complexity and diversity.,18 (7): 1112- 1126.

    Ishii, Y., Zhu, Z. B., Macon, K. J., and Volanakis, J. E., 1993. Structure of the human C2 gene., 151 (1): 170-174.

    andMulti-com- ponent complement system of Cnidaria: C3, Bf, and MASP genes expressed in the endodermal tissues of a sea anemone, . , 214: 165-

    Larkin, M. A., Blackshields, G., Brown, N. P., Chenna, R., Mc- Gettigan, 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: 2947-2948.

    Nonaka, M., 1994. Molecular analysis of the lamprey complement system., 4: 437-446.

    Nonaka, M., 2014. Evolution of the complement system., 80: 31-43.

    Nonaka, M., and Kimura, A., 2006. Genomic view of the evolution of the complement system., 58: 701-713.

    Nonaka, M., and Yoshizaki, F., 2004. Primitive complement system of invertebrates.,198: 203-215.

    Pangburn, M. K., 1983. Activation of complement via the alternative pathway.,42: 139-143.

    Pfaffl, M. W., Horgan, G. W., and Dempfle, L., 2002. Relative expression software tool (REST(c)) for group-wise comparison and statistical analysis of relative expression results in real-time PCR., 30 (9): e36.

    Prado-Alvarez, M., Rotllant, J., Gestal, C., Novoa, B., and Figueras, A., 2009. Characterization of a C3 and a factor B-like in the carpet-shell clam,., 26: 305-315.

    Ramírez-Gómez, F., Aponte-Rivera, F., Méndez-Castaner, L., and García-Arrarás J. E., 2010. Changes inpopulations following immune stimulation with different molecular patterns., 29: 175-185.

    Ramírez-Gómez, F., Ortíz-Pineda, P. A., Rivera-Cardona, G., andGarcía-Arrarás, J. E., 2009. LPS-induced genes in intestinal tissue of the sea cucumber., 4: 61-78.

    Reid, K. B., and Day, A. J., 1989. Structure-function relationships of the complement components., 10: 177-180.

    Santiago-Cardona, P. G., Berríos, C. A., Ramírez, F., andGar- cía-Arrarás, J. E., 2003. Lipopolysaccharides induce intestinal serum amyloid A expression in the sea cucumber.,27: 105-110.

    Sekine, H., Kenjo, A., Azumi, K., Ohi, G., Takahashi, M., Kasukawa, R., Ichikawa, N., Nakata, M., Mizuochi, T., Matsushita, M., Endo, Y., and Fujita, T., 2001. An ancient lectin-dependent complement system in an ascidian: Novel lectin isolated from the plasma of the solitary ascidian, Halocynthia roretzi.,167: 4504-4510.

    Smith, L. C., Shih, C. S., and Dachenhausen, S. G., 1998. Coe- lomocytes express SpBf, a homologue of factor B, the second component in the sea urchin complement system., 161: 6784-6793.

    Tagawa, K., Yoshihara, T., Shibata, T., Kitazaki, K., Endo, Y., Fujita, T., Koshiba, T., and Kawabata, S., 2012. Microbe- specific C3b deposition in the horseshoe crab complement system in a C2/factor B-dependent or -independent manner., 7 (5): e36783.

    Tamura, K., Dudley, J., Nei, M., and Kumar, S., 2007. MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0., 24: 1596-1599.

    Tuckwell, D. S., Xu, Y., Newham, P., Humphries, M. J., and Volanakis, J. E., 1997. Surface loops adjacent to the cation- binding site of the complement factor B von Willebrand factor type A module determine C3b binding specificity., 36: 6605-6613.

    Wei, W., Wu, H., Xu, H., Xu, T., Zhang, X., Chang, K., and Zhang, Y., 2009. Cloning and molecular characterization of two complement Bf/C2 genes in large yellow croaker ()., 27: 285- 295.

    Welsh, R. M., 1977. Host cell modification of lymphocytic choriomeningitis virus and Newcastle disease virus altering viral inactivation by human complement.,118: 348-354.

    Xu, Y., and Volanakis, J. E., 1997. Contribution of the complement control protein modules of C2 in C4b binding assessed by analysis of C2/factor B chimeras., 158: 5958-5965.

    Yang, A. F., Zhou, Z. C., He, C. B., Hu, J. J., Chen, Z., Gao, X. G., Dong, Y., Jiang, B., Liu. W. D., Guan, X. Y., and Wang, X. Y., 2009. Analysis of expressed sequence tags from body wall, intestine and respiratory tree of sea cucumber ()., 296: 193-199.

    Yoshizaki, F. Y., Ikawa, S., Satake, M., Satoh, N., and Nonaka, M., 2005. Structure and the evolutionary implication of the triplicated complement factor B genes of a urochordate ascidian,., 56: 930-942.

    Zhong, L., Zhang, F., and Chang, Y., 2012. Gene cloning and function analysis of complement B factor-2 of., 33: 504-513.

    Zhou, Z. C., Dong, Y., Sun, H. J., Yang, A. F., Chen, Z., Gao, S., Jiang, J. W., Guan, X. Y., Jiang, B., and Wang, B., 2014. Transcriptome sequencing of sea cucumber (Apostichopus japonicus) and the identification of gene-associated markers., 14 (1): 127-138.

    Zhou, Z. C., Sun, D. P., Yang, A. F., Dong, Y., Chen, Z., Wang, X. Y., Guan, X. Y., Jiang, B., and Wang, B., 2011. Molecular characterization and expression analysis of a complement component 3 in the sea cucumber ()., 31 (4): 540-547.

    Zhu, Y., Thangamani, S., Ho, B., and Ding, J. L., 2005. The ancient origin of the complement system.,24: 382-394.

    (Edited by Qiu Yantao)

    DOI 10.1007/s11802-015-2696-8

    ISSN 1672-5182, 2015 14 (6): 1096-1104

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

    (July 11, 2014; revised November 5, 2014; accepted September 1, 2015)

    * Corresponding author. Tel: 0086-411-84678047 E-mail: zunchunz@hotmail.com

    精品一区二区三区av网在线观看| 日本黄色视频三级网站网址| 精品国产美女av久久久久小说| 黄色视频,在线免费观看| 99久久无色码亚洲精品果冻| 97超级碰碰碰精品色视频在线观看| 看免费av毛片| 麻豆成人av在线观看| avwww免费| 一区福利在线观看| 亚洲av第一区精品v没综合| 99热6这里只有精品| 国产色婷婷99| 999久久久精品免费观看国产| 久久精品国产亚洲av香蕉五月| 51午夜福利影视在线观看| 国产麻豆成人av免费视频| 夜夜躁狠狠躁天天躁| 欧美日韩综合久久久久久 | 免费观看精品视频网站| 国产亚洲精品一区二区www| 久久久久久人人人人人| 成年女人永久免费观看视频| av女优亚洲男人天堂| 男女那种视频在线观看| 国产一区二区三区视频了| 国产91精品成人一区二区三区| 欧美黑人巨大hd| 老司机福利观看| 欧美日韩国产亚洲二区| 成人三级黄色视频| 18禁美女被吸乳视频| 最新在线观看一区二区三区| 亚洲午夜理论影院| 亚洲成人精品中文字幕电影| 两个人视频免费观看高清| 美女高潮的动态| 天堂av国产一区二区熟女人妻| 国内精品美女久久久久久| 久久中文看片网| 中国美女看黄片| 久久亚洲精品不卡| 69人妻影院| 国产精品亚洲av一区麻豆| netflix在线观看网站| 嫩草影视91久久| 久久香蕉国产精品| 麻豆成人午夜福利视频| 在线国产一区二区在线| 国产精品98久久久久久宅男小说| 国产精品影院久久| 黄色日韩在线| www日本在线高清视频| 国产精品久久久久久亚洲av鲁大| 91字幕亚洲| 国产私拍福利视频在线观看| 我要搜黄色片| 精品欧美国产一区二区三| 亚洲人成网站在线播| 日日夜夜操网爽| 无限看片的www在线观看| 夜夜看夜夜爽夜夜摸| 亚洲美女黄片视频| 国产成年人精品一区二区| 天美传媒精品一区二区| 男人和女人高潮做爰伦理| 人妻久久中文字幕网| 啦啦啦韩国在线观看视频| 岛国在线观看网站| 一区二区三区高清视频在线| 日韩 欧美 亚洲 中文字幕| 国产精品一区二区三区四区久久| 淫秽高清视频在线观看| 午夜福利在线观看吧| 日韩有码中文字幕| 欧美日韩亚洲国产一区二区在线观看| 小蜜桃在线观看免费完整版高清| 少妇高潮的动态图| 国产高潮美女av| 午夜福利欧美成人| 99在线人妻在线中文字幕| 亚洲av熟女| 午夜激情福利司机影院| 亚洲国产中文字幕在线视频| 人人妻人人看人人澡| h日本视频在线播放| 亚洲 国产 在线| 久久这里只有精品中国| 床上黄色一级片| 亚洲五月天丁香| av女优亚洲男人天堂| 欧美中文综合在线视频| 99热精品在线国产| 久久久久国产精品人妻aⅴ院| 我的老师免费观看完整版| 一级a爱片免费观看的视频| www.999成人在线观看| 国内毛片毛片毛片毛片毛片| 亚洲最大成人手机在线| 欧美色欧美亚洲另类二区| 黄色丝袜av网址大全| 香蕉久久夜色| av天堂中文字幕网| 久久欧美精品欧美久久欧美| 99久久精品热视频| 欧美日韩中文字幕国产精品一区二区三区| 国产精品久久久久久人妻精品电影| 一进一出抽搐动态| 日韩欧美在线二视频| 99精品在免费线老司机午夜| 亚洲国产精品合色在线| 亚洲一区二区三区色噜噜| 亚洲成人中文字幕在线播放| 免费电影在线观看免费观看| 午夜福利免费观看在线| 最近最新免费中文字幕在线| 老司机午夜十八禁免费视频| 男女那种视频在线观看| 欧美日韩瑟瑟在线播放| 91字幕亚洲| 免费av毛片视频| 性色av乱码一区二区三区2| 在线观看日韩欧美| 亚洲人成网站在线播| 9191精品国产免费久久| 麻豆成人av在线观看| 午夜福利18| 国产精品一区二区三区四区久久| 欧美激情在线99| aaaaa片日本免费| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 99久久无色码亚洲精品果冻| aaaaa片日本免费| 免费在线观看日本一区| 小说图片视频综合网站| 搡老岳熟女国产| 中文字幕熟女人妻在线| 99在线视频只有这里精品首页| 国产av一区在线观看免费| 欧美日本视频| 国内少妇人妻偷人精品xxx网站| 免费电影在线观看免费观看| 亚洲黑人精品在线| 亚洲成a人片在线一区二区| 99精品在免费线老司机午夜| www日本在线高清视频| 色在线成人网| 哪里可以看免费的av片| 久久久久国内视频| 亚洲人成网站在线播| 黄色成人免费大全| 国内揄拍国产精品人妻在线| 香蕉丝袜av| 狂野欧美白嫩少妇大欣赏| 最新美女视频免费是黄的| 成人国产综合亚洲| 婷婷精品国产亚洲av在线| 亚洲av美国av| 亚洲七黄色美女视频| 99热6这里只有精品| 亚洲男人的天堂狠狠| 美女被艹到高潮喷水动态| 精品一区二区三区人妻视频| 亚洲激情在线av| 成人午夜高清在线视频| 久久草成人影院| 9191精品国产免费久久| 国产高清视频在线播放一区| 久久久久久久久大av| 极品教师在线免费播放| 18禁黄网站禁片午夜丰满| 精品久久久久久久久久免费视频| 婷婷精品国产亚洲av| 精品久久久久久,| 女人高潮潮喷娇喘18禁视频| 久久精品国产自在天天线| 午夜影院日韩av| 日韩免费av在线播放| 精品一区二区三区视频在线观看免费| 亚洲中文字幕一区二区三区有码在线看| 我的老师免费观看完整版| 成人国产一区最新在线观看| 97碰自拍视频| 日本五十路高清| 男人舔女人下体高潮全视频| 国产精品,欧美在线| 亚洲熟妇熟女久久| 黄色视频,在线免费观看| 日韩高清综合在线| 国产aⅴ精品一区二区三区波| 亚洲精品一卡2卡三卡4卡5卡| 在线观看一区二区三区| 99国产综合亚洲精品| 无人区码免费观看不卡| 搡老熟女国产l中国老女人| 色在线成人网| 丝袜美腿在线中文| svipshipincom国产片| 亚洲色图av天堂| 在线观看一区二区三区| 老鸭窝网址在线观看| 99热这里只有精品一区| 国内精品一区二区在线观看| 久久国产乱子伦精品免费另类| 国产精品女同一区二区软件 | 国产不卡一卡二| 两个人视频免费观看高清| 成熟少妇高潮喷水视频| 欧美精品啪啪一区二区三区| 五月伊人婷婷丁香| 偷拍熟女少妇极品色| 国产69精品久久久久777片| 欧美精品啪啪一区二区三区| 欧美在线黄色| 啪啪无遮挡十八禁网站| 精品久久久久久成人av| 精品福利观看| 久久久成人免费电影| 母亲3免费完整高清在线观看| 国产精品永久免费网站| 黄色成人免费大全| 国产精品香港三级国产av潘金莲| 18美女黄网站色大片免费观看| 国产综合懂色| 久久香蕉精品热| 国产精品,欧美在线| 别揉我奶头~嗯~啊~动态视频| 窝窝影院91人妻| 精品人妻偷拍中文字幕| 丰满乱子伦码专区| 国产成人av教育| 亚洲一区二区三区不卡视频| 亚洲乱码一区二区免费版| 少妇裸体淫交视频免费看高清| 99热这里只有精品一区| 精品不卡国产一区二区三区| 男女做爰动态图高潮gif福利片| 老鸭窝网址在线观看| 国产精品久久电影中文字幕| 成人高潮视频无遮挡免费网站| 国产一区在线观看成人免费| 69av精品久久久久久| 午夜福利成人在线免费观看| av天堂中文字幕网| 热99在线观看视频| 精品久久久久久久末码| 级片在线观看| 亚洲 欧美 日韩 在线 免费| 免费电影在线观看免费观看| 亚洲精华国产精华精| 亚洲av电影不卡..在线观看| 精品福利观看| 亚洲片人在线观看| 色精品久久人妻99蜜桃| 国产久久久一区二区三区| 又紧又爽又黄一区二区| 亚洲18禁久久av| 国产欧美日韩一区二区精品| 可以在线观看毛片的网站| 午夜福利高清视频| 综合色av麻豆| 一级黄片播放器| 国产午夜精品论理片| 丁香欧美五月| 免费观看人在逋| 欧美另类亚洲清纯唯美| 网址你懂的国产日韩在线| 国产免费av片在线观看野外av| 国产精品久久视频播放| 真实男女啪啪啪动态图| 久久久久久人人人人人| 久久精品影院6| bbb黄色大片| 午夜精品一区二区三区免费看| 精品电影一区二区在线| 我要搜黄色片| 亚洲男人的天堂狠狠| 熟女电影av网| 国产av不卡久久| 51国产日韩欧美| 好看av亚洲va欧美ⅴa在| 亚洲成人精品中文字幕电影| 成人午夜高清在线视频| 午夜久久久久精精品| 亚洲成人久久性| 丁香欧美五月| 在线观看美女被高潮喷水网站 | 欧美成人性av电影在线观看| 国内久久婷婷六月综合欲色啪| 99久久九九国产精品国产免费| 久久婷婷人人爽人人干人人爱| 免费观看的影片在线观看| 成人一区二区视频在线观看| 大型黄色视频在线免费观看| 听说在线观看完整版免费高清| 婷婷丁香在线五月| 国产亚洲精品久久久久久毛片| 色吧在线观看| 黄色片一级片一级黄色片| 亚洲最大成人中文| 美女黄网站色视频| 麻豆成人午夜福利视频| 一本综合久久免费| 国产成人系列免费观看| 欧美成人a在线观看| 精品久久久久久久久久免费视频| 国产乱人视频| 欧美成人性av电影在线观看| 国产综合懂色| 大型黄色视频在线免费观看| 免费大片18禁| 欧美中文综合在线视频| 国产99白浆流出| 最近视频中文字幕2019在线8| 成熟少妇高潮喷水视频| 午夜两性在线视频| 一个人免费在线观看的高清视频| www国产在线视频色| 黄色日韩在线| 欧美日韩中文字幕国产精品一区二区三区| 亚洲精品456在线播放app | 国产真实乱freesex| 99久久99久久久精品蜜桃| 国内揄拍国产精品人妻在线| 高清毛片免费观看视频网站| 午夜日韩欧美国产| 久久精品夜夜夜夜夜久久蜜豆| 99精品久久久久人妻精品| 日韩国内少妇激情av| 国产精品免费一区二区三区在线| 精品一区二区三区视频在线 | 亚洲国产欧美人成| 性欧美人与动物交配| 亚洲av二区三区四区| 久久精品影院6| 国内精品久久久久久久电影| 激情在线观看视频在线高清| 日韩欧美三级三区| 一进一出抽搐动态| 又黄又爽又免费观看的视频| 99热只有精品国产| 日本撒尿小便嘘嘘汇集6| 国产精品一区二区三区四区免费观看 | 波多野结衣巨乳人妻| 一级a爱片免费观看的视频| 日韩欧美精品免费久久 | 香蕉久久夜色| 精品国产三级普通话版| 精品一区二区三区视频在线 | 偷拍熟女少妇极品色| 国产高清videossex| 亚洲精品亚洲一区二区| 国产精品女同一区二区软件 | 九色成人免费人妻av| 首页视频小说图片口味搜索| 欧美日韩综合久久久久久 | 人妻丰满熟妇av一区二区三区| 午夜免费成人在线视频| 午夜老司机福利剧场| 中文字幕av在线有码专区| 久久亚洲精品不卡| 日韩有码中文字幕| 国产精品一区二区三区四区免费观看 | 午夜免费成人在线视频| 性色avwww在线观看| 亚洲午夜理论影院| 久久久成人免费电影| 国产日本99.免费观看| 小说图片视频综合网站| 乱人视频在线观看| 亚洲天堂国产精品一区在线| 最好的美女福利视频网| 久久九九热精品免费| 国产精品国产高清国产av| 97超视频在线观看视频| 久久精品国产99精品国产亚洲性色| 1000部很黄的大片| 51午夜福利影视在线观看| 夜夜躁狠狠躁天天躁| 18禁黄网站禁片免费观看直播| 精品欧美国产一区二区三| 成人国产综合亚洲| 亚洲18禁久久av| 日本三级黄在线观看| 国产乱人视频| 18禁裸乳无遮挡免费网站照片| 免费看光身美女| 天堂√8在线中文| 99热这里只有精品一区| 欧美+亚洲+日韩+国产| 中国美女看黄片| 亚洲国产日韩欧美精品在线观看 | 亚洲第一电影网av| 国产野战对白在线观看| 国产亚洲精品久久久com| 老司机深夜福利视频在线观看| 欧美黑人欧美精品刺激| 欧美黄色淫秽网站| 亚洲成人久久性| 午夜福利在线观看吧| 天天一区二区日本电影三级| 国产精品三级大全| 丰满人妻熟妇乱又伦精品不卡| 人人妻人人澡欧美一区二区| 1024手机看黄色片| 在线a可以看的网站| 香蕉丝袜av| 天堂网av新在线| 日韩欧美一区二区三区在线观看| 97超级碰碰碰精品色视频在线观看| 亚洲中文字幕一区二区三区有码在线看| 欧美在线一区亚洲| netflix在线观看网站| 亚洲熟妇中文字幕五十中出| 法律面前人人平等表现在哪些方面| 亚洲精品一区av在线观看| 欧美日韩福利视频一区二区| 三级男女做爰猛烈吃奶摸视频| 日本免费a在线| 中文字幕av在线有码专区| 精品乱码久久久久久99久播| 欧美高清成人免费视频www| 日韩成人在线观看一区二区三区| 日韩欧美精品免费久久 | 精品福利观看| 性色av乱码一区二区三区2| 精品久久久久久久人妻蜜臀av| 久久精品夜夜夜夜夜久久蜜豆| 亚洲成a人片在线一区二区| 每晚都被弄得嗷嗷叫到高潮| 日韩高清综合在线| 国产亚洲精品综合一区在线观看| 亚洲国产欧洲综合997久久,| svipshipincom国产片| 91av网一区二区| 老司机午夜十八禁免费视频| av在线蜜桃| 高清毛片免费观看视频网站| 欧美激情久久久久久爽电影| 亚洲国产欧美人成| 国产色爽女视频免费观看| 热99在线观看视频| 久久久国产成人精品二区| 男人舔奶头视频| 欧美日韩福利视频一区二区| 亚洲 欧美 日韩 在线 免费| 一区二区三区免费毛片| 欧美激情久久久久久爽电影| 日本a在线网址| 成人永久免费在线观看视频| 丰满乱子伦码专区| 丁香六月欧美| www.999成人在线观看| 国产在线精品亚洲第一网站| 久久99热这里只有精品18| 中文字幕人妻熟人妻熟丝袜美 | 一区二区三区免费毛片| 老司机深夜福利视频在线观看| 国产一区二区三区在线臀色熟女| 日日夜夜操网爽| 国产欧美日韩精品亚洲av| 制服人妻中文乱码| 丰满的人妻完整版| 亚洲精品乱码久久久v下载方式 | 免费看日本二区| 午夜福利在线在线| 51国产日韩欧美| 国产爱豆传媒在线观看| 久久久国产成人免费| 欧美日韩黄片免| 国产亚洲精品久久久com| av天堂中文字幕网| 啪啪无遮挡十八禁网站| 在线十欧美十亚洲十日本专区| 叶爱在线成人免费视频播放| 亚洲成av人片在线播放无| 九九热线精品视视频播放| 国产伦精品一区二区三区四那| 亚洲18禁久久av| 黑人欧美特级aaaaaa片| 免费av毛片视频| 丰满人妻熟妇乱又伦精品不卡| 男女下面进入的视频免费午夜| 老司机午夜福利在线观看视频| 欧美日本亚洲视频在线播放| 亚洲欧美激情综合另类| 欧美+亚洲+日韩+国产| 热99re8久久精品国产| 悠悠久久av| 好男人在线观看高清免费视频| 欧美日韩黄片免| 精品欧美国产一区二区三| 亚洲一区二区三区色噜噜| 真人一进一出gif抽搐免费| av天堂中文字幕网| 国产精品 国内视频| 老司机在亚洲福利影院| 亚洲国产精品成人综合色| 国产精品国产高清国产av| 岛国在线免费视频观看| 好看av亚洲va欧美ⅴa在| 欧美性感艳星| 欧美日韩国产亚洲二区| av在线天堂中文字幕| netflix在线观看网站| 天天添夜夜摸| 国产亚洲精品av在线| 国产精品精品国产色婷婷| 国产亚洲精品综合一区在线观看| 午夜福利免费观看在线| 天堂av国产一区二区熟女人妻| 欧美性猛交黑人性爽| 欧美高清成人免费视频www| 最近视频中文字幕2019在线8| 91在线观看av| 乱人视频在线观看| 中文资源天堂在线| 亚洲无线在线观看| 中文字幕高清在线视频| 九九热线精品视视频播放| 人人妻人人看人人澡| 露出奶头的视频| 欧美乱码精品一区二区三区| 亚洲av一区综合| 18+在线观看网站| 精品久久久久久久毛片微露脸| 在线观看舔阴道视频| 制服丝袜大香蕉在线| netflix在线观看网站| 国产一区二区在线av高清观看| 啦啦啦观看免费观看视频高清| 国产精品一区二区三区四区久久| 免费观看精品视频网站| 床上黄色一级片| 少妇的逼水好多| 乱人视频在线观看| 国产免费男女视频| 欧美zozozo另类| 成年版毛片免费区| 18禁美女被吸乳视频| 亚洲av二区三区四区| 国产成人影院久久av| 身体一侧抽搐| 国产精品美女特级片免费视频播放器| 好男人电影高清在线观看| av片东京热男人的天堂| 亚洲自拍偷在线| 成熟少妇高潮喷水视频| 久久精品国产亚洲av香蕉五月| 级片在线观看| 美女黄网站色视频| 天堂网av新在线| 久久久国产精品麻豆| 91在线观看av| 亚洲国产精品成人综合色| 国产精品亚洲一级av第二区| 最近最新免费中文字幕在线| 亚洲一区二区三区不卡视频| 国产黄a三级三级三级人| 看片在线看免费视频| 日韩欧美三级三区| 法律面前人人平等表现在哪些方面| 国产一区二区在线观看日韩 | 亚洲中文字幕日韩| 亚洲国产色片| 欧美日本亚洲视频在线播放| 老司机福利观看| 国产精品一及| 国产激情偷乱视频一区二区| 国产精品精品国产色婷婷| 国产精品乱码一区二三区的特点| 久久亚洲真实| 久久香蕉精品热| 别揉我奶头~嗯~啊~动态视频| av中文乱码字幕在线| 成人亚洲精品av一区二区| 午夜免费观看网址| 日韩大尺度精品在线看网址| 色综合站精品国产| 国产亚洲精品av在线| 狠狠狠狠99中文字幕| 日本成人三级电影网站| 欧美色视频一区免费| av天堂在线播放| 亚洲avbb在线观看| 亚洲精品色激情综合| 亚洲成av人片在线播放无| 亚洲国产精品成人综合色| 日韩中文字幕欧美一区二区| 亚洲成人久久爱视频| 亚洲成人精品中文字幕电影| 亚洲男人的天堂狠狠| 日本精品一区二区三区蜜桃| 色综合婷婷激情| 欧美激情在线99| 亚洲 欧美 日韩 在线 免费| 三级毛片av免费| 精华霜和精华液先用哪个| 中文字幕人妻熟人妻熟丝袜美 | 久久久久久久亚洲中文字幕 | 国产毛片a区久久久久| 18禁美女被吸乳视频| 丰满人妻熟妇乱又伦精品不卡| 美女黄网站色视频| 国产精品野战在线观看| 在线观看美女被高潮喷水网站 | 国产精品野战在线观看| 观看美女的网站| 国产精品亚洲一级av第二区| 可以在线观看的亚洲视频| 日韩欧美精品v在线| 级片在线观看| 欧美性猛交黑人性爽| 少妇高潮的动态图| 美女cb高潮喷水在线观看|