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

    Isolation, Expression and Characterization of rbcL Gene from Ulva prolifera J. Agardh (Ulvophyceae, Chlorophyta)

    2015-04-01 01:58:24SHAOZhanruLIWeiGUOHuiandDUANDelin
    Journal of Ocean University of China 2015年6期

    SHAO Zhanru, LI Wei, 2), GUO Hui, 3), and DUAN Delin, *

    ?

    Isolation, Expression and Characterization ofGene fromJ. Agardh (Ulvophyceae, Chlorophyta)

    SHAO Zhanru1), LI Wei1), 2), GUO Hui1), 3), and DUAN Delin1), *

    1),,266071,2),,07743,3),100049,

    is a typical green alga in subtidal areas and can grow tremendously fast. A highly efficient Rubisco enzyme which is encoded by UpRbcL gene may contribute to the rapid growth. In this study, the full-lengthopen reading frame (ORF) was identified, which encoded a protein of 474 amino acids. Phylogenetic analysis ofsequences revealed that Chlorophyta had a closer genetic relationship with higher plants than with Rhodophyta and Phaeophyta. The two distinct residues (aa11 and aa91) were presumed to be unique for Rubisco catalytic activity. The predicted three-dimensional structure showed that one α/β-barrel existed in the C-terminal region, and the sites for Mg2+coordination and CO2fixation were also located in this region. Gene expression profile indicated thatwas expressed at a higher level under light exposure than in darkness. When the culture temperature reached 35℃, the expression level ofwas 2.5-fold lower than at 15℃, and the carboxylase activity exhibited 13.8-fold decrease.was heterologously expressed inand was purified by Ni2+affinity chromatography. The physiological and biochemical characterization of recombinant Rubisco will be explored in the future.

    Ribulose-1, 5-bisphosphate carboxylase/oxygenase large subunit; sequence analysis; real-time PCR;expression;

    1 Introduction

    (Ulvophyceae, Chlorophyta) is one of the typical green algae which niches in subtidal areas along the seashore of China, and commonly attached to solid substratum. In recent years,bloom has been frequently occurring in Yellow Sea (Liu.,2009; Duan.,2012), and caused serious problems on coastal ecosystems (Sun., 2008). Green algal thalli exhibit rapid growth features during their floating state on the sea. It is believed thatpossesses particular physiological and ecological characteristics to maintain long period survival through high efficient photosynthesis (Lin., 2011; Zhang., 2012).

    Ribulose-1, 5-bisphosphate carboxylase/oxygenase (Ru- bisco, EC 4.1.1.39) is a crucial enzyme in Calvin cycle of photosynthesis (Hartman and Harpel, 1994; Houtz and Portis, 2003); it is involved in the first step of carbon fixation in photosynthetic and chemoautotrophic organisms. Rubisco can simultaneously catalyze the carboxylation and oxygenation of ribulose bisphosphate (RuBP). RuBP carboxylation defines gross photosynthesis rate and determines productivity (Jensen and Bahr, 1977). Most reports to date have focused on unicellular algae, including(Dron., 1982),(Ginrich and Hallick, 1985), and, and few molecular studies concerning macroalgal Rubisco are available (Kai., 2002; Ying., 2011). Therefore, we studied the multicellular green algafor better understanding how Rubisco regulates its rapid growth and also the dynamics of real-time expression and enzymatic activity.

    There are four forms of Rubisco (I, II, III and IV), and that of Chlorophyta is believed to be form I which comprises of eight large catalytic subunits (RbcLs) and eight small subunits (RbcSs) (Tabita., 2008). In terrestrial plants and green algae, RbcLs (50–55kD) are encoded by a singlegene present in the chloroplast genome, whereas the rbcSs (12–18kD) are encoded by nuclear genes (Dean., 1989; Spreitzer and Salvucci, 2002). So far, it is reported that RbcLs provide all the catalytically essential residues of Rubisco, while the function of RbcSs is not clear (Kaplan and Reinhold, 1999).

    In this study, we isolated and characterized the gene encoding.RbcL protein. Variation ingene expression was monitored under different conditions, and recombinant protein was purified after it was expressed in.. Our research made contribution to decipher relationship between UpRbcL structure and function, and to explore physiological mechanisms supporting rapid growth of.during floating adaptation.

    2 Materials and Methods

    2.1 Pre-culture and treatments of

    Fresh.were collected from the green algal blooming area on the sea in Qingdao (36?03′N; 120?18′E) in August 2011. The algal filaments were maintained indoor at 25℃, with a continuous illumination of 45μmol m?2s?1under a light exposure of 12:12 L/D for 48h. Then.a cells were collected and subjected to different temperature conditions (10, 15, 20, 25, 30, 35℃) for 1h to test the influence of temperature on gene expression and carboxylase activity. For the light treatment, the alga was incubated at 25℃ under continuous illumination (45μmolm?2s?1) for 48h, and then in dark for 12h. After these treatments, the alga was collected and rinsed with 0.1% DEPC sterile seawater 3 times, frozen in liquid nitrogen and stored at ?80℃ for nucleic acid and crude protein extractions.

    2.2 DNA, RNA Extraction and cDNA Synthesis

    Genomic DNA was extracted with the Plant Genomic DNA Extraction Kit (Tiangen, Beijing, China) following the manufacturer’s instructions. Total RNA extraction was carried out using CTAB (cetyltrimethylammonium bromide) method (Yao., 2009). Extracted RNA was desiccated and dissolved in 30μL RNase-free water, and any possible genomic DNA contamination was removed with the DNase I Kit (Fermentas, Burlington, USA). First strand cDNA was synthesized using the PrimeScript II cDNA Synthesis Kit (Takara, Tokyo, Japan) according to the manufacturer’s instructions, and stored at ?20℃ for subsequentisolation.

    2.3 Cloning of

    One pair of specific primers, UpRbcL-F (5’-ATGG CTCCACAGACTG-3’) and UpRbcL-R (5’-TAATGTA TCAATAGCATCAAATTC-3’), was designed based on thesequence of.(DQ813497) to amplify full-length coding sequence ofcDNA. PCR was performed in a 20μL reaction volume containing 1μL of cDNA, 0.25μLTMDNA polymerase high fidelity (HiFi) (TransGen, Beijing, China), 2μL of 10×HiFi buffer, 1μL of each primer (10μmolL?1), 3.2μL of dNTP (10mmolL?1), and 11.55μL of RNase-free water. The amplification was carried out by predenaturing at 95℃ for 5min, followed by 30 cycles of denaturing at 95℃ for 30s, 47℃ for 30s, 72℃ 2min, and a final extension at 72℃ for 10min. PCR product was sequenced by Sangon (Shanghai, China).

    The same pair of primers, UpRbcL-F and UpRbcL-R, was used to amplify thesequence from genomic DNA, and PCR product was treated through the procedure described above.

    2.4 Analysis of UpRbcL Sequence

    Sequence was analyzed for molecular weight (MW), isoelectric point (pI) with software ORF Finder (Sayers., 2012) and ProtParam (Gasteiger., 2005). The similarity of UpRbcL ORF with others was checked by BLAST analysis (Altschul., 1997) at NCBI database (http://www.ncbi.nlm.nih.gov/blast). Signal peptide of UpRbcL sequence was predicted by SignalP 4.0 Server (Petersen., 2011), and presence of trans-membrane domains was determined with the TMHMM version 2.0 program (Krogh., 2001). Hydrophobicity and hydro- philicity were analyzed by ProtScale program (Gasteiger., 2005), and the secondary structure was predicted by SOPMA (Geourjon and Deleage, 1995). SWISS- MODEL (Peitsch, 1995; Arnold., 2006; Kiefer., 2009) and Pymol Viewer programs were applied to construct the putative 3D-structure of UpRbcL by homology modeling with.

    2.5 Phylogenetic Analysis

    Phylogenetic tree was constructed using the neighbor- joining algorithm by MEGA 5.2 program (Tamura., 2011) and bootstrap values were calculated from 1000 replicates. UpRbcL amino acid sequence was aligned with RbcL sequences from Archaeplastida. Clustal X 2.0 software (Larkin., 2007; Thompson., 1994) was used to perform further alignments.

    2.6 Transcriptional Profile of

    Two specific primers UpRbcL-S (5’-ATGGCTCCACA GACTGAAAC-3’) and UpRbcL-A (5’-CCTGGTTGAG GAGTCATACGG-3’) were applied for amplifying a 139 bp amplicon. Primers UpLhcSR-S (5’-ACTGGTCCATT TGGTTTCTT-3’) and UpLhcSR-A (5’-CCGCTGATGTT ACCCTCCC-3’) were designed according to the coding sequence of light-harvesting chlorophyll a/b-binding protein (LhcSR) (HQ622096.1). Transcriptional levels were normalized by the expression of the gene encoding the small subunit rRNA amplified with the primers 28S-S (5’-AACACGGACCAAGGAGTCTAAC-3’) and 28S-A (5’-GAAACTTCGGAGGGAACCAG-3’).

    RT-qPCR was performed with the SYBRTMII (Takara, Tokyo Japan) on the TP800 Thermal Cycler Dice? (Takara, Tokyo, Japan). The PCR mixture (25μL) contained 12.5μL Premix, 1μL of each primer (10μmolL?1), 2μL of diluted cDNA and 8.5μL of RNase-free water. Thermal cycling protocol was: 95℃ for 30s, followed by 40 cycles of 95℃ for 5s and 58℃ for 30s. Specificity of primers was checked by observing relevant dissociation curve. Three independent biological replicates were carried out for each sample, and relative quantitative values were calculated with 2?△△Ctmethod (Schmittgen., 2000).

    2.7 Measurement of Carboxylase Activity

    Extraction of total protein from.was conducted according to the method described by Rousvoal(2011). Protein content was determined according to the Bradford method (Bradford, 1976) using theProtein Quantitative Kit (TransGen, Beijing, China).

    To determine Rubisco carboxylase activity, we referred to the method described by Lan and Mott (1991), and applied the Plant Rubisco Carboxylation Activity Detection Kit (GenMed Scientifics, China). The reaction was conducted at 25℃ for 5min in a Microplate Reader (BioTek, Vermont, USA). Five μg of algal protein extract was added into the reaction mixture (250μL), and the reaction was initiated by adding the substrate RuBP. Production of NADH was determined by measuring ultraviolet (UV) absorbance at 340nm, and the extinction coefficient of 6.22mmolL?1was used to determine the initial carboxylase activity (Unit: mmolL?1CO2mg?1min?1). Boiled algal protein extract was used as negative control.

    2.8 Heterologous Expression of

    Specific primers UpRbcL-F (5’-ATGGCTCCACAGA CTG-3’) and UpRbcL-R (5’-TAATGTATCAATAGCA TCAAATTC-3’) were used for amplifying. The thermal cycling applied was: 94℃ for 5 min, followed by 30 cycles of 94℃ for 30s, 48℃ for 30s, 72℃ for 2min, and a final extension step of 72℃ for 10min. The positive PCR product was subcloned into the expression vector-E2 (TransGen, Beijing, China). The recombinant plasmid was transformed into.BL21 (DE3) and transformants was spread on antibiotic selection media containing ampicillin (0.1mgmL?1final concentration). Plates were incubated at 37℃ overnight, and positive clones were screened by colony PCR detection.

    To obtain production of the recombinant protein, pre- culture of 2mL of LB medium containing ampicillin was incubated overnight at 37℃, and was then used to start a fresh 200mL LB plus ampicillin culture, until the OD600valuereached to 0.5. Expression ofwas induced by adding 0.5mmolL?1IPTG into the medium at 120 rpm for hours at 37℃.

    In order to optimize the duration of induction and the production of recombinant protein, 500μL of culture were collected every 1h and then centrifuged at 4000rmin?1for 15min. Fifty μL of deionized water and an equal volume of 2×sodium dodecyl-sulfate polyacrylamide gel electrophoresis loading buffer (Tris-HCl 25mmolL?1, DTT 100mmolL?1, SDS 2%, Glycerol 20%, Bromophenol Blue 0.016%) were added to the pellet. After boiling for 10min, mixtures were loaded on a 12% SDS-PAGE (ATTO, Tokyo, Japan). Electrophoresis was carried out at a constant voltage of 80 V for 20min, and then 160V for 60min. Bacterial protein extracts produced from.containing the empty vector were used as negative control.

    2.9 Western Blot Analysis

    After migration on SDS-PAGE, proteins from the gels were electrophoretically transferred to polyvinylidene difluoride (PVDF) membrane (0.45μm) at 30V for 2h with an AE-6675 HorizBlot system (ATTO, Tokyo, Japan). After rinsing with TBST buffer (150mmolL?1NaCl, 20mmolL?1Tris-HCl pH 8.0, 0.05% Tween-20) 3 times, membrane was blocked to avoid non-specific binding by placing it in a freshly prepared 5% solution of non-fat milk (5g dry milk dissolved in up to 100mL TBST buffer), and shaking at 30rmin?1for 2h at room temperature. The membrane was then rinsed twice with TBST buffer for 5min. To detect the target protein, a 1:250 dilution of primary rabbit anti-His antibody (Top- Science, Beijing, China) was incubated with the membrane under gentle agitation for 2h at room temperature. After additional rinsing with TBST, membranes were exposed to a 1:2500 dilution of goat-anti-rabbit IgG antibody conjugated with horseradish peroxidase (Tiandz, Beijing, China) for 1h. DAB sensitive chromogenic reaction (Tiandz, Beijing, China) was used to detect binding of antibodies, and reaction was terminated by transfer of the membrane into deionized water.

    2.10 Purification of Recombinant UpRbcL

    After induction with 0.5mmolL?1IPTG for 5h, LB plus ampicillin culture (200mL) was centrifuged at 4000rmin?1for 15min. Pellet was re-suspended in the buffer which contained 300mmolL?1NaCl, 50mmolL?1NaH2PO4, 10mmolL?1imidazole and 10mmolL?1Tris base (pH 8.0). The suspension was sonicated intermittently and centrifuged at 15000for 20min at 4℃. The supernatant was filtered and loaded on theTMNi-NTA Resin (TransGen, Beijing, China). Elution was conducted by applying successively imidazole solutions (40, 80, 120, 160, 200mmolL?1), and fractions were loaded on 12% SDS-PAGE. All the purification steps were carried out at 4℃. Concentrations of UpRbcL recombinant protein were quantified according to the method described above.

    3 Results

    3.1 Sequence Analysis of UpRbcL

    Full-length open-reading frame of(1425bp), encoding a protein of 474 amino acids, was submitted to GenBank (JQ867403.1). Comparison between sequences obtained from cDNA and genomic DNA showed that there was no intron in. BLAST analysis confirmed thatencoded a protein belonging to Rubisco large superfamily, with a theoretical molecular mass of 52.34kD and an isoelectric point (pI) of 5.82. UpRbcL was considered as a comparatively stable enzyme with an instability index of 38.71.

    Maximal values of the original shearing site (C score), signal peptide (S score), and synthesized shearing site (Y score) were 0.115, 0.105 and 0.123, respectively, meaning that there was no signal peptide in UpRbcL amino acid sequence. With the transmembrane topological struc- ture analysis, all parts of UpRbcL were identified to be outside the membrane.

    In relation to the secondary structure prediction, proportions of α-helix, extended strand, β-turn and random coil of UpRbcL were 42.83%, 15.82%, 9.92%, and 31.43%, respectively, indicating that α-helix and random coil were its major components. The tertiary structure of UpRbcL was constructed with SWISS-MODEL and Pymol Viewer prediction programs based on the structure template of(Taylor, 2001). It was obvious that one α/β-barrel existed in the C-terminal end, and 7 sites (K174, R294, H297, H326, K333, G402 and G403) in this region were responsible for the RuBP activation and catalization (Fig.1). The two adjacent acidic residues D202 and E203 were assumed for Mg2+coordination, and the role of residue K200 was for CO2fixation (Fig.1).

    Fig.1 Predicted 3D structure of UpRbcL. (a) Global tertiary structure of UpRbcL. The white rectangle indicated the α/β-barrel in the C-terminal end. (b) Zoom in on the white rectangle region which contained the barrel and serial important amino acid residues: D202 and E203 were reported to coordinate Mg2+; K200 provided a site for CO2 fixation; amino acids important for catalysis of carboxylation reaction (K174, R294, H297, H326, K333, G402 and G403) were presented in red color.

    3.2 Phylogenetic Analysis of UpRbcLs

    Alignment of amino acid sequences showed that UpRbcL was more closely related to that of other green algae and land plants than to that of,and(Fig.2). The long C-terminal in..and.indicated that brown and red alga sequences might have originated from the same common ancestor. Thr 11 was unique in RbcLs sequences from, while Ala or Val were located at site 11 in other species of Archaeplastida. A deletion of amino acid was observed at site 91. It was deduced that both sites might have function on the catalytic characterization of Rubisco from other species.

    A phylogenetic tree was constructed based on 14sequences from algae and land plants (Fig.3). It shows that UpRbcL was in the subgroup of green algae and shared more than 99% identities with those from.and..alsoexhibited 75.7% identities withand 77.1% with. In reference to the phylogeny of RbcLs,was more closely related to terrestrial plants than to glaucophytes and Rhodophyta.

    3.3 Transcriptional Profile of

    Real-time PCR detection ofandtranscripts was conducted on algal samples harvested after incubation at various temperatures and illuminations. It was shown that both genes were highly expressed after exposure to continuous light illumination for 48h, while their abundance decreased significantly after being kept in dark for 12h (1.5-fold and 7.3-fold change, respectively) (<0.05) (Fig.4a). When the alga was incubated at a range of temperatures, from 10℃ to 35℃, transcript abundance ofwas high at 10℃, 15℃, 20℃ and 25℃, with a maximum at 15℃, and declined significantly (2.5-fold) at 35℃ (<0.05) (Fig.4b).

    3.4 Measurement of Endogenous Rubisco Activity

    After extracting crude protein from.(0.1g each) treated at 15℃, 25℃, 35℃, respectivelywe conducted SDS-PAGE electrophoresis. The bands at 55kD show the potential presence of Rubisco, as it is the major protein within crude extract (Fig.5a). In addition, carboxylase activity of Rubisco was determined using equal amounts of total protein for the three conditions tested. It was observed that the activity reached maximum at 15℃, half at 25℃, and almost completely abolished at 35℃ (Fig.5b).

    3.5 Heterologous Expression and Purification of Recombinant UpRbcL

    After optimizing the protein induction time, heterologous protein was significantly expressed in one hour after inducing with 0.5mmolL?1IPTG, and it became stably expressed within 2 to 5h. After being separated on 12% SDS-PAGE, 55kD protein band was distinctly identified from the recombinant bacterial extract while no other band correspondingly appeared after induction of bacteria containing empty plasmid (Fig.6a). The optimal imidazole concentration for the elution of purified UpRbcL was 80mmolL?1(Fig.6b), and western blot indicated that the 55kD band was positive to the RbcL anti-His antibody (Fig.6c).

    Fig.2 Alignments of Ulva RbcL with those of land plants and algae. The green rectangle indicates a unique characterization at site 11; the red rectangle indicates a Pro deletion distinctively in the Rubisco activase contact region in U. porlifera, respectively; the blue rectangle shows highlights a long C-terminal extension in brown and red algal sequences.

    Fig.3 Phylogenetic tree of RbcL subunits constructed with selected sequences from Archaeplastida. The tree is obtained by the neighbor-joining algorithm using MEGA 5.2 program. Bootstrap values are calculated from 1000 replicates, and the branch lengths are proportional to estimated evolutionary change. The scale bar corresponds to 0.05 estimated amino-acid substitution per site. Green letters show the species having Green-like Rubisco, and red letters indicates Red-like Rubisco. The sequences include Ulva prolifera (JQ867403.1), U. linza (DQ813496.2), U. pertusa (JQ867404.1), Chlorella pyrenoidosa (EU038283.1), Chlamydomonas reinhardtii (ACJ50136.1), Volvox carteri (ACI31253.1), Zea mays (NP_043033.1), Oryza sativa (NP_039391.1), Arabidopsis thaliana (NP_051067.1), Nicotiana tabacum (NP_ 054507.1), Spinacia oleracea (NP_054944.1), Cyanophora biloba (KF631394.1), Cyanophora paradoxa (NP_043240.1), Pyropia yezoensis (DQ227866.1).

    Fig.4Influence of illumination and temperature on UpRbcL and UpLhcSR expression level. (a) Relative transcript abundance of UpRbcL and UpLhcSR exposed to different illuminations. (b) Relative abundance of UpRbcL transcript at different temperatures. All data are expressed as mean±SD (n=3).

    4 Discussion

    Recently, green tides unceasingly occur on large scale in the due seasons each year (Leliaert, 2008). In China, green tide occurs almost every year in Yellow Sea since 2007 (Liu, 2009; Duan, 2012). More and more investigations have been carried out to identify the reason behind this naturally occurring phenomenon. Nevertheless, previous works were concentrated on the morphological features and the origination survey (Zhao., 2011; Duan., 2012). In present study, we isolated and characterized, and carried out the over-ex- pression and purification of this large subunit protein from..

    There was no intron in the sequence of, which was previously observed in.(Ying., 2011). From the evolutionary tree of form I Rubisco, taxa are commonly divided into two groups. One clade included sequences from Chlorophyta and land plants (Green-like Rubisco group), while the other was Rhodophyta (Red- like Rubisco group). Our results showed that although belonging to Archaeplastida, green algae have closer evo- lutionary relationship with land plants than with Glaucophyte algae.

    Being crucial to the photosynthesis, UpRbcL was supposed to be highly efficientand the major contributing factor in sustaining the rapid growth of.. In this study, we found that several crucial residues (K174, R294, H297, H326, K333, G402 and G403) in α/β-barrel were responsible for the catalysis of Rubisco, and these residues were consistent with early reports (Lorimer and Miziorko, 1980; Kellogg.,1997; Cleland., 1998). Previously, genetic screening has identified some residues influencing the specificity and stability of Rubisco: V331/ T342/ G344/ R339/ I393 in loop 6, D473 in C-terminus and A222/ V262/ L290 between large and small subunits (Spreitzer and Salvucci, 2002; Satagopan and Spreitzer, 2004; Karkehabadi, 2007; Andersson and Backlund, 2008). These residues were all found in UpRbcL sequence, which showed high conservation of the active sites under positive selection. However, two amino acid residues (aa11 and aa91) were unique in UpRbcL compared to other RbcL sequences. According to the previous publications (Li., 2005; Portis., 2008), we suspected that the deletion of Pro at site 91 might affect the activation of Rubisco due to the formation of an extra α-helix between the sites 89 and 93. Moreover, the substitution of hydrophobic amino acid A11 (in,,.) or V11 (in land plants) into hydrophilic T11 was presumed to affect the secondary structure of loop 1, and eventually influence the catalytic efficiencies and thermal stability. Site-specific mutagenesis is required to verify the function of these sites in the future.

    When the green tide outbreaks,can float on the sea surface and rapidly grow within a few weeks..may possess unique physiological and biochemical properties that allow it to adapt to persistent illumination on the sea. We applied excessive light exposure, up to 48h, to monitor the transcription ofand. We observed that both genes were expressed at a higher level when alga was exposed to extensive irradiance. This indicated thatandcan be transcribed under continuous illumination, and might be coupled into complicated physiological process for utilizing excessive solar energy for rapid growth.

    Comparative analysis on transcript abundance and enzymatic activity was carried out. The maximumexpression level and the initial activity of Rubisco both emerged in the range of 15 to 20℃, which was consistent with the ecological features of.(Wang, 2007). After 1h at 35℃,mRNA expression and carboxylase activity were highly decreased; indicating limited heat-resisting ability of UpRubisco. Referred to the studies on rice and spinach (Schneider, 1992; Li., 2002), we assumed that these thermal properties might due to the thermal instabilities of Rubisco activase.

    For the prokaryotic expression, the content of recombination enzyme reached maximum whenwas induced with 0.5mmolL?1IPTG for 2h, and the recombinant UpRbcL could be completely eluted with 80mmolL?1imidazole. So far, no attempt to reconstitute active homologous enzyme has succeeded for form I Rubisco to our knowledge (Whitney, 2011), owing to the indispensability of chaperone and small subunit (Tarnawski, 2008; Liu, 2010). Although no activity of UpRubisco could be measured in the present study, we purified UpRbcL and paved the way for the successful overexpression of the whole enzyme UpRubisco.

    In conclusion, this study isolated and characterized Rubisco large subunit gene of. Features detected in its secondary and tertiary structures, such as conservation of important catalytic residues, were consistent with its high photosynthetic efficiencies. Moreover, high transcription levels ofandunder continuous illumination should be part of the mechanisms favoring rapid growth of.. These molecular observations provided clues supporting the fast development of the green algal bloom, and pave the way for further investigation of key physiological processes in..

    Acknowledgements

    This research was financially supported by the Scientific and Technical Supporting Programs of China (2008 BAC49B01) and State Ocean Administration Project (200805069). The authors appreciated Drs. Thierry Tonon and Zimin Hu for revision and suggestions of the manuscript.

    Altschul, S. F., Madden, T. L., Schaffer, A. A., Zhang, J. H., Zhang, Z., Miller, W., and Lipman, D. J., 1997. Gapped BLAST and PSI-BLAST: A new generation of protein database search programs., 25: 3389- 3402.

    Andersson, I., and Backlund, A., 2008. Structure and function of Rubisco., 46: 275-291.

    Arnold, K., Bordoli, L., Kopp, J., and Schwede, T., 2006. The SWISS-MODEL Workspace: A web-based environment for protein structure homology modelling., 22: 195-201.

    Bradford, M. M., 1976. Rapid and sensitive method for quantitation of microgram quantities of protein utilizing principle of protein-dye binding.,72: 248-254.

    Cleland, W. W., Andrews, T. J., Gutteridge, S., Hartman, F. C., and Lorimer, G. H., 1998. Mechanism of Rubisco: The carbamate as general base., 98: 549-561.

    Dean, C., Pichersky, E., and Dunsmuir, P., 1989. Structure, evolution, and regulation of RbcS genes in higher plants., 40: 415-439.

    Dron, M., Rahire, M., and Rochaix, J. D., 1982. Sequence of the chloroplast DNA region ofcon- taining the gene of the large subunit of ribulose bisphosphate carboxylase and parts of its flanking genes.,162: 775-793.

    Duan, W. J., Guo, L. X., Sun, D., Zhu, S. F., Chen, X. F., Zhu, W. R., Xu, T., and Chen, C. F., 2012. Morphological and molecular characterization of free-floating and attached green macroalgaespp. in the Yellow Sea of China.,24: 97-108.

    Gasteiger, E., Hoogland, C., Gattiker, A., Duvaud, S., Wilkins, M. R., Appel, R. D., and Bairoch, A., 2005. Protein iden- tification and analysis tools on the ExPASy Server. In:. Walker, J. M., ed., Humana Press , New Jersey, 571-607.

    Geourjon, C., and Deleage, G., 1995. SOPMA: Significant improvements in protein secondary structure prediction by consensus prediction from multiple alignments., 11 (6): 681-684.

    Ginrich, J. C., and Hallick, R. B., 1985. The Euglena gracilis chloroplast ribulose-1, 5-bisphosphate carboxylase gene. II. The spliced mRNA and its product.,260: 16162-16168.

    Hartman, F. C., and Harpel, M. R., 1994. Structure, function, regulation, and assembly of D-ribulose-1,5-bisphosphate carbo- xylase/oxygenase.,63: 197- 234.

    Houtz, R. L., and Portis, A. R., 2003. The life of ribulose 1,5-bisphosphate carboxylase/oxygenase-posttranslational facts and mysteries., 414: 150-158.

    Jensen, G. R., and Bahr, T. J., 1977. Ribulose 1.5-bisphosphate carboxylase-oxygenase.,28: 379-400.

    Kai, B., Gudrun, K., Christian, W., and Dieter, H., 2002. Solar ultraviolet radiation affects the activity of ribulose-1,5- bisphosphate carboxylase-oxygenase and the composition of photosynthetic and xanthophyll cycle pigments in the inter- tidal green algaL., 215: 502-509.

    Kaplan, A., and Reinhold, L., 1999. CO2concentrating mecha- nisms in photosynthetic microorganisms., 50: 539-570.

    Karkehabadi, S., Satagopan, S., Taylor, T. C., Spreitzer, R. J., and Andersson, I., 2007. Structural analysis of altered large- subunit loop-6/carboxy-terminus interactions that influence catalytic efficiency and CO2/O2specificity of ribulose 1,5- bisphosphate carboxylase/oxygenase.,46: 11080- 11089.

    Kellogg, E. A., and Juliano, N. D., 1997. The structure and function of Rubisco and their implications for systematic studies., 84: 413-428.

    Kiefer, F., Arnold, K., Künzli, M., Bordoli, L., and Schwede, T., 2009. The SWISS-MODEL Repository and associated resources., 37: 387-392.

    Krogh, A., Larsson, B., von Heijne, B., and Sonnhammer, E. L. L., 2001. Predicting transmembrane protein topology with a hidden Markov model: Application to complete genomes.,305 (3): 567-580.

    Lan, Y., and Mott, K. A., 1991. Determination of apparent KM values for ribulose 1,5-bisphosphate carboxylase oxygenase (Rubisco) activase using the spectrophotometric assay of Rubisco activity., 95: 604-609.

    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: 2947-2948.

    Leliaert, F., Malta, E. J., Engelen, A. H., Mineur, F., and De Clerck, O., 2008. Qingdao algal bloom culprit identified., 56: 1516.

    Li, C. H., Salvucci, M. E., and Portis, A. R., 2005. Two residues of Rubisco activase involved in recognition of the Rubisco substrate., 280: 24864-24869.

    Li, G., Mao, H. B., Ruan, X., Xu, Q., Gong, Y. D., Zhang, X. F., and Zhao, N. M., 2002. Association of heat-induced confor- mational change with activity loss of Rubisco.,290: 1128-1132.

    Lin, A. P., Shen, S. D., Wang, G. C., Yi, Q. Q., Qiao, H. J., Niu, J. F., and Pan, G. H., 2011. Comparison of chlorophyll and photosynthesis parameters of floating and attached., 53 (1): 25-34.

    Liu, C. M., Young, A. L., Starling-Windhof, A., Bracher, A., Saschenbrecker, S., Rao, B. V., Rao, K. V., Berninghausen, O., Mielke, T., Hartl, F. U., Beckmann, R., and Hayer-Hartl, M., 2010. Coupled chaperone action in folding and assembly of hexadecameric Rubisco., 463: 197-202.

    Liu, D. Y., Keesing, J. K., Xing, Q. U., and Shi, P., 2009. World’s largest macroalgal bloom caused by expansion of seaweed aquaculture in China.,58: 888-895.

    Lorimer, G. H., and Miziorko, H. M., 1980. Carbamate for- mation on the epsilon-amino group of a lysyl residue as the basis for the activation of ribulosebisphosphate carboxylase by CO2and Mg2+., 19 (23): 5321-5328.

    Peitsch, M. C., 1995. Protein modeling by E-mail.,13: 658-660.

    Petersen, T. N., Brunak, S., von Heijne, G., and Nielsen, H., 2011. SignalP 4.0: Discriminating signal peptides from trans- membrane regions.,8: 785-786.

    Portis, A. R., Li, C. S., Wang, D. F., and Salvucci, M. E., 2008. Regulation of Rubisco activase and its interaction with Rubisco., 59: 1597-1604.

    Rousvoal, S., Groisillier, A., Dittami, S. M., Michel, G., Boyen, C., and Tonon, T., 2011. Mannitol-1-phosphate dehydro- genase activity in, a key role for mannitol synthesis in brown algae.,233: 261-273.

    Satagopan, S., and Spreitzer, R. J., 2004. Substitutions at the Asp-473 latch residue ofribulosebisphosphate carboxylase/oxygenase cause decreases in carboxylation and CO2/O2specificity.,279: 14240-14244.

    Sayers, E. W., Barrett, T., Benson, D. A., Bolton, E., Bryant, S. H., Canese, K., Chetvernin, V., Church, D. M., DiCuccio, M., Federhen, S., Feolo, M., Fingerman, I. M., Geer, L. Y., Helmberg, W., Kapustin, Y., Krasnov, S., Landsman, D., Lipman, D. J., Lu, Z. Y., Madden, T. L., Madej, T., Maglott, D. R., Marchler-Bauer, A., Miller, V., Karsch-Mizrachi, I., Ostell, J., Panchenko, A., Phan, L., Pruitt, K. D., Schuler, G. D., Sequeira, E., Sherry, S. T., Shumway, M., Sirotkin, K., Slotta, D., Souvorov, A., Starchenko, G., Tatusova, T. A., Wagner, L., Wang, Y. L., Wilbur, W. J., Yaschenko, E., and Ye, J., 2012. Database resources of the National Center for Biotechnology Information.,40: D13-D25.

    Schmittgen, T. D., Zakrajsek, B. A., Mills, A. G., Gorn, V., Singer, M. J., and Reed, M. W., 2000. Quantitative reverse transcription-polymerase chain reaction to study mRNA decay: Comparison of endpoint and real-time methods.,285: 194-204.

    Schneider, S. H., 1992. The climatic response to greenhouse gases., 22: 1-32.

    Spreitzer, R. J., and Salvucci, M. E., 2002. Rubisco: Structure, regulatory interactions, and possibilities for a better enzyme., 53: 449-75.

    Sun, S., Wang, F., Li, C. L., Qin, S., Zhou, M. J., Ding, L. P., Pang, S. J., Duan, D. L., Wang, G. C., Yin, B. S., Yu, R. C., Jiang, P., Liu, Z. L., Zhang, G. T., Fei, X. G., and Zhou, M., 2008. Emerging challenges: Massive green algae blooms in the Yellow Sea.,hdl:10101/npre.2008. 2266. 1.

    Tabita, F. R., Satagopan, S., Hanson, T. E., Kreel, N. E., and Scott, S. S., 2008. Distinct form I, II, III, and IV Rubisco proteins from the three kingdoms of life provide clues about Rubisco evolution and structure/function relationships.,59: 1515-1524.

    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: 2731-2739.

    Tarnawski, M., Gubernator, B., Kolesinski, P., and Szczepaniak, A., 2008. Heterologous expression and initial characterization of recombinant RbcX protein fromBP-1 and the role of RbcX in Rubisco assembly., 55: 777-785.

    Taylor, T. C., Backlund, A., Bjorhall, K., Spreitzer, R. J., and Andersson I., 2001 First crystal structure of Rubisco from a green alga,.,276 (51): 48159-48164.

    Thompson, J. D., Higgins, D. G., and Gibson, T. J., 1994. CLUSTAL W: Improving the sensitivity of progressive mul- tiple sequence alignment through sequence weighting, posi- tion-specific gap penalties and weight matrix choice.,22: 4673-4680.

    Wang, J. W., Yan, B. L., Lin, A. P., Hu, J. P., and Shen, S. D., 2007. Ecological factor research on the growth and induction of spores release in(Chlorophyta)., 26 (2): 60-65 (in Chinese with English abstract).

    Whitney, S. M., Houtz, R. L., and Alonso, H., 2011. Advancing our understanding and capacity to engineer nature’s CO2-se- questering enzyme, Rubisco., 155: 27-35.

    Yao, J. T., Fu, W. D., Wang, X. L., and Duan, D. L., 2009. Improved RNA isolation fromAresch (Laminariaceae, Phaeophyta).,21 (2): 233-238.

    Ying, C. Q., Yin, S. J., Shen, Y., Lin, S. J., and He, P. M., 2011. Cloning and analysis of the full-length Rubisco large subunit (rbcL) cDNA from(Chlorophyceae, Chlorophyta)., 54: 303-312.

    Zhang, H. N., Li, W., Li, J. J., Fu, W. D., Yao, J. T., and Duan, D. L., 2012. Characterization and expression analysis of hsp70 gene fromJ. Agardh (Chlorophycophyta, Chlorophyceae)., 5: 53-58.

    Zhao, J., Jiang, P., Liu, Z. Y., Wang, J. F., Cui, Y. L., and Qing, S., 2011. Genetic variation of()(Ulvales, Chlorophyta)-the causative species of the green tides in the Yellow Sea, China., 23: 227-233.

    (Edited by Qiu Yantao)

    DOI 10.1007/s11802-015-2661-6

    ISSN 1672-5182, 2015 14 (6): 1087-1095

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

    (May 9, 2014; revised July 7, 2014; accepted September 10, 2015)

    * Corresponding author. Tel: 0086-532-82898554 E-mail: dlduan@ms.qdio.ac.cn

    99国产综合亚洲精品| 巨乳人妻的诱惑在线观看| 午夜福利,免费看| 日本黄色日本黄色录像| 满18在线观看网站| 首页视频小说图片口味搜索 | 日韩 亚洲 欧美在线| 热re99久久精品国产66热6| 国产淫语在线视频| 亚洲欧美精品自产自拍| 国产精品亚洲av一区麻豆| 99国产精品99久久久久| 丝瓜视频免费看黄片| 国产91精品成人一区二区三区 | 黄色一级大片看看| 免费看十八禁软件| 涩涩av久久男人的天堂| 欧美xxⅹ黑人| 日本a在线网址| 欧美xxⅹ黑人| 母亲3免费完整高清在线观看| 久久女婷五月综合色啪小说| 国产精品国产av在线观看| 一区福利在线观看| 亚洲成人免费av在线播放| 纯流量卡能插随身wifi吗| 一边摸一边做爽爽视频免费| 九色亚洲精品在线播放| 下体分泌物呈黄色| 一边摸一边抽搐一进一出视频| 国产片特级美女逼逼视频| 一区二区三区乱码不卡18| 精品卡一卡二卡四卡免费| 一本大道久久a久久精品| 成人国产一区最新在线观看 | 久久午夜综合久久蜜桃| 99热全是精品| a级片在线免费高清观看视频| 欧美精品一区二区大全| 亚洲国产精品一区三区| 国产一级毛片在线| 国产欧美日韩精品亚洲av| 十八禁人妻一区二区| 午夜久久久在线观看| 黄色视频在线播放观看不卡| av视频免费观看在线观看| 国产福利在线免费观看视频| 美女扒开内裤让男人捅视频| 人妻人人澡人人爽人人| 一个人免费看片子| 亚洲黑人精品在线| 国产一区二区 视频在线| 九色亚洲精品在线播放| 99国产精品免费福利视频| 这个男人来自地球电影免费观看| 亚洲欧洲精品一区二区精品久久久| 777久久人妻少妇嫩草av网站| 日韩 欧美 亚洲 中文字幕| 成年美女黄网站色视频大全免费| 久久久久国产一级毛片高清牌| netflix在线观看网站| 亚洲国产精品国产精品| 肉色欧美久久久久久久蜜桃| 日韩中文字幕视频在线看片| 在现免费观看毛片| 国产xxxxx性猛交| 亚洲一码二码三码区别大吗| 中文字幕最新亚洲高清| 欧美 日韩 精品 国产| 亚洲精品美女久久久久99蜜臀 | 丁香六月天网| 你懂的网址亚洲精品在线观看| 少妇人妻 视频| 少妇猛男粗大的猛烈进出视频| 黄色 视频免费看| 午夜福利,免费看| 欧美日韩视频精品一区| 中文字幕高清在线视频| 国产又爽黄色视频| 免费人妻精品一区二区三区视频| 久久国产精品人妻蜜桃| 亚洲精品国产av蜜桃| 视频区欧美日本亚洲| 亚洲国产av新网站| 不卡av一区二区三区| 日本欧美视频一区| 亚洲精品自拍成人| 免费在线观看日本一区| 精品亚洲成a人片在线观看| 夫妻性生交免费视频一级片| bbb黄色大片| 国产成人精品久久二区二区免费| 在线观看一区二区三区激情| 亚洲国产精品一区二区三区在线| 国产精品欧美亚洲77777| 久久久久网色| 天天影视国产精品| 汤姆久久久久久久影院中文字幕| 我的亚洲天堂| 欧美+亚洲+日韩+国产| 波野结衣二区三区在线| 大片电影免费在线观看免费| 热re99久久国产66热| 亚洲欧美一区二区三区国产| 国产免费又黄又爽又色| 亚洲av欧美aⅴ国产| 国产精品欧美亚洲77777| 午夜两性在线视频| 国产精品一区二区在线观看99| 女人精品久久久久毛片| 欧美精品亚洲一区二区| 欧美人与善性xxx| 亚洲一卡2卡3卡4卡5卡精品中文| 国产成人精品久久二区二区免费| 久久精品熟女亚洲av麻豆精品| 99国产精品一区二区蜜桃av | 亚洲精品美女久久av网站| 日本午夜av视频| 久久99热这里只频精品6学生| 999久久久国产精品视频| 天天添夜夜摸| 一级毛片电影观看| 国产高清视频在线播放一区 | 99香蕉大伊视频| 美女福利国产在线| 久久久国产欧美日韩av| 日韩av在线免费看完整版不卡| 日本猛色少妇xxxxx猛交久久| 日韩中文字幕欧美一区二区 | 国产精品熟女久久久久浪| 欧美精品高潮呻吟av久久| 亚洲七黄色美女视频| 亚洲国产中文字幕在线视频| 性色av一级| 国产欧美亚洲国产| 国产成人欧美在线观看 | 日韩大码丰满熟妇| 亚洲 欧美一区二区三区| 国产爽快片一区二区三区| 色94色欧美一区二区| 99热网站在线观看| 一级黄片播放器| 亚洲精品美女久久久久99蜜臀 | 国产xxxxx性猛交| 久久久久久久久久久久大奶| 午夜福利视频精品| 欧美黑人精品巨大| 视频区图区小说| 亚洲自偷自拍图片 自拍| 考比视频在线观看| 桃花免费在线播放| 亚洲三区欧美一区| 99国产综合亚洲精品| 亚洲精品国产一区二区精华液| 99热全是精品| 你懂的网址亚洲精品在线观看| 国产精品欧美亚洲77777| 国产视频一区二区在线看| 亚洲视频免费观看视频| 性少妇av在线| 天堂俺去俺来也www色官网| 国产精品九九99| 操美女的视频在线观看| 老司机在亚洲福利影院| 亚洲视频免费观看视频| 日韩大片免费观看网站| 超碰97精品在线观看| 午夜福利在线免费观看网站| 亚洲欧洲精品一区二区精品久久久| 国产精品久久久久久精品古装| 亚洲国产精品999| 亚洲精品日韩在线中文字幕| 男女高潮啪啪啪动态图| 午夜免费成人在线视频| 免费女性裸体啪啪无遮挡网站| 日本wwww免费看| 下体分泌物呈黄色| 久久国产精品大桥未久av| 2018国产大陆天天弄谢| 999精品在线视频| 精品国产一区二区三区久久久樱花| 国产男女内射视频| 丝瓜视频免费看黄片| 国产精品久久久久久人妻精品电影 | 丝袜脚勾引网站| 中文字幕人妻丝袜一区二区| 成人三级做爰电影| 婷婷成人精品国产| 交换朋友夫妻互换小说| av国产精品久久久久影院| 丝袜喷水一区| 别揉我奶头~嗯~啊~动态视频 | 无遮挡黄片免费观看| 电影成人av| 看免费av毛片| 国产免费又黄又爽又色| 午夜老司机福利片| 日韩中文字幕欧美一区二区 | 亚洲精品一卡2卡三卡4卡5卡 | 大型av网站在线播放| 国产成人av激情在线播放| 涩涩av久久男人的天堂| 青草久久国产| 精品福利观看| 99精国产麻豆久久婷婷| 亚洲欧美日韩高清在线视频 | 中文字幕制服av| 亚洲国产精品999| 香蕉丝袜av| 麻豆乱淫一区二区| 欧美日韩亚洲高清精品| 国产一区有黄有色的免费视频| 久久人妻熟女aⅴ| 亚洲欧美成人综合另类久久久| 夫妻午夜视频| a级毛片黄视频| 男女午夜视频在线观看| 人人妻人人爽人人添夜夜欢视频| 国产成人av教育| 在线观看免费高清a一片| 自线自在国产av| 欧美激情极品国产一区二区三区| 满18在线观看网站| 成人国语在线视频| 国产99久久九九免费精品| 欧美成人午夜精品| 精品人妻在线不人妻| 大型av网站在线播放| 老司机午夜十八禁免费视频| 亚洲欧洲精品一区二区精品久久久| 丝袜脚勾引网站| 免费在线观看影片大全网站 | 国产成人av激情在线播放| 日日夜夜操网爽| 国产老妇伦熟女老妇高清| 男人舔女人的私密视频| 热99久久久久精品小说推荐| 久久女婷五月综合色啪小说| 国产成人一区二区三区免费视频网站 | 日韩大码丰满熟妇| 亚洲成国产人片在线观看| 亚洲色图综合在线观看| 不卡av一区二区三区| 丝袜在线中文字幕| 亚洲国产精品国产精品| 捣出白浆h1v1| 亚洲色图 男人天堂 中文字幕| 欧美日韩视频高清一区二区三区二| 亚洲欧美清纯卡通| 午夜影院在线不卡| 免费高清在线观看日韩| 亚洲一卡2卡3卡4卡5卡精品中文| 校园人妻丝袜中文字幕| 曰老女人黄片| 国产xxxxx性猛交| 中文字幕人妻丝袜制服| 国产成人av教育| av国产久精品久网站免费入址| 晚上一个人看的免费电影| www.av在线官网国产| 欧美日韩福利视频一区二区| √禁漫天堂资源中文www| av不卡在线播放| 国产日韩欧美在线精品| 人成视频在线观看免费观看| 色视频在线一区二区三区| 一级黄片播放器| 一区二区日韩欧美中文字幕| 国产av精品麻豆| 18禁裸乳无遮挡动漫免费视频| 亚洲精品久久午夜乱码| 老汉色∧v一级毛片| 亚洲精品美女久久久久99蜜臀 | 一二三四在线观看免费中文在| 超碰97精品在线观看| 亚洲精品久久成人aⅴ小说| 天堂中文最新版在线下载| 最新在线观看一区二区三区 | 涩涩av久久男人的天堂| 精品国产一区二区三区久久久樱花| 又大又黄又爽视频免费| 久久久久久久精品精品| 99热全是精品| kizo精华| 久久精品成人免费网站| 啦啦啦在线观看免费高清www| 亚洲欧美中文字幕日韩二区| 91成人精品电影| 免费女性裸体啪啪无遮挡网站| 成人三级做爰电影| 亚洲精品第二区| 五月开心婷婷网| 久久综合国产亚洲精品| 1024香蕉在线观看| 国产人伦9x9x在线观看| 一区二区日韩欧美中文字幕| 国产激情久久老熟女| 97在线人人人人妻| 日本wwww免费看| 每晚都被弄得嗷嗷叫到高潮| 两性夫妻黄色片| 又紧又爽又黄一区二区| 久久免费观看电影| 色婷婷久久久亚洲欧美| 亚洲欧美一区二区三区黑人| 久久精品亚洲熟妇少妇任你| 波野结衣二区三区在线| 色播在线永久视频| 久久久久精品人妻al黑| 久久青草综合色| 国产不卡av网站在线观看| 高清av免费在线| 午夜福利视频在线观看免费| 精品国产一区二区久久| 悠悠久久av| 只有这里有精品99| 久久久精品区二区三区| 大香蕉久久成人网| 亚洲精品在线美女| 久久性视频一级片| 婷婷色av中文字幕| 欧美亚洲 丝袜 人妻 在线| 亚洲午夜精品一区,二区,三区| 少妇裸体淫交视频免费看高清 | 曰老女人黄片| 国产精品久久久久久人妻精品电影 | 精品国产乱码久久久久久小说| 丰满饥渴人妻一区二区三| 美女大奶头黄色视频| 又黄又粗又硬又大视频| 精品国产乱码久久久久久男人| 美女福利国产在线| 国产精品一区二区在线观看99| 亚洲 欧美一区二区三区| 欧美黑人欧美精品刺激| 亚洲情色 制服丝袜| 久久久欧美国产精品| 国产亚洲欧美精品永久| avwww免费| 国产精品人妻久久久影院| 嫩草影视91久久| 免费观看a级毛片全部| 黄网站色视频无遮挡免费观看| 日韩av不卡免费在线播放| 手机成人av网站| 高清黄色对白视频在线免费看| 亚洲精品国产av成人精品| 国产精品.久久久| 国产99久久九九免费精品| 男女床上黄色一级片免费看| 久久青草综合色| 老汉色av国产亚洲站长工具| 日日摸夜夜添夜夜爱| 啦啦啦 在线观看视频| 青青草视频在线视频观看| 国产精品国产av在线观看| 久久久精品94久久精品| 亚洲精品一二三| 伦理电影免费视频| 亚洲男人天堂网一区| 中文字幕最新亚洲高清| 人人妻人人爽人人添夜夜欢视频| 久久久久久久精品精品| 美女高潮到喷水免费观看| 日日爽夜夜爽网站| 性高湖久久久久久久久免费观看| 丁香六月欧美| 国产男女超爽视频在线观看| 91九色精品人成在线观看| 国产高清不卡午夜福利| 女人爽到高潮嗷嗷叫在线视频| 成人亚洲欧美一区二区av| 99久久99久久久精品蜜桃| 成人亚洲精品一区在线观看| 精品国产国语对白av| 色网站视频免费| 一区二区三区乱码不卡18| 国产熟女欧美一区二区| 五月天丁香电影| 欧美av亚洲av综合av国产av| 一区二区av电影网| 午夜免费成人在线视频| 可以免费在线观看a视频的电影网站| 国产精品麻豆人妻色哟哟久久| 首页视频小说图片口味搜索 | 中国国产av一级| 欧美黄色片欧美黄色片| 国产av国产精品国产| 在线天堂中文资源库| 看免费av毛片| 精品久久久久久久毛片微露脸 | 九色亚洲精品在线播放| 一区二区三区精品91| 黄网站色视频无遮挡免费观看| 一区二区三区精品91| 妹子高潮喷水视频| 最黄视频免费看| 婷婷成人精品国产| 久久久精品国产亚洲av高清涩受| 无遮挡黄片免费观看| 久久久久视频综合| 久热这里只有精品99| 伦理电影免费视频| 久久久久精品人妻al黑| 好男人视频免费观看在线| 亚洲人成77777在线视频| 亚洲国产欧美网| 一级毛片黄色毛片免费观看视频| www日本在线高清视频| 久久久久久久久久久久大奶| 一区二区三区激情视频| 亚洲一卡2卡3卡4卡5卡精品中文| av天堂在线播放| 两人在一起打扑克的视频| 在线观看一区二区三区激情| 男女高潮啪啪啪动态图| 久久精品国产亚洲av涩爱| 啦啦啦视频在线资源免费观看| 成人免费观看视频高清| 亚洲欧美清纯卡通| 日韩人妻精品一区2区三区| 97精品久久久久久久久久精品| 大香蕉久久成人网| 成年av动漫网址| 亚洲免费av在线视频| 巨乳人妻的诱惑在线观看| 亚洲国产中文字幕在线视频| 亚洲av美国av| 国产免费视频播放在线视频| 中文字幕人妻丝袜制服| 狂野欧美激情性bbbbbb| 成人黄色视频免费在线看| 中文字幕制服av| 日本vs欧美在线观看视频| 国产在线观看jvid| 欧美另类一区| 久久久久精品人妻al黑| 一级毛片我不卡| 男女下面插进去视频免费观看| 久久国产精品影院| 观看av在线不卡| 亚洲国产精品国产精品| 一边亲一边摸免费视频| 免费在线观看完整版高清| 日日摸夜夜添夜夜爱| 人体艺术视频欧美日本| 美女国产高潮福利片在线看| 久久99精品国语久久久| 欧美+亚洲+日韩+国产| 国产成人欧美在线观看 | 我的亚洲天堂| 国产免费现黄频在线看| 久久精品亚洲熟妇少妇任你| 国产精品国产三级专区第一集| 晚上一个人看的免费电影| 性少妇av在线| 精品人妻在线不人妻| 国产成人av激情在线播放| 色94色欧美一区二区| 亚洲精品在线美女| 只有这里有精品99| 丝瓜视频免费看黄片| 精品久久蜜臀av无| 久久狼人影院| 50天的宝宝边吃奶边哭怎么回事| 亚洲成色77777| 中文字幕精品免费在线观看视频| 日韩av不卡免费在线播放| 国产一级毛片在线| 国产精品亚洲av一区麻豆| 国产精品香港三级国产av潘金莲 | 91精品三级在线观看| 亚洲av片天天在线观看| 欧美精品人与动牲交sv欧美| 热re99久久精品国产66热6| 日本一区二区免费在线视频| 亚洲欧美色中文字幕在线| 成人影院久久| 久久鲁丝午夜福利片| 国语对白做爰xxxⅹ性视频网站| 母亲3免费完整高清在线观看| 超碰97精品在线观看| 操美女的视频在线观看| 久久亚洲精品不卡| 久久精品久久精品一区二区三区| 久久人妻熟女aⅴ| 国产精品免费视频内射| 制服诱惑二区| 午夜福利视频在线观看免费| 男女国产视频网站| 国产一区二区三区综合在线观看| 欧美 亚洲 国产 日韩一| 亚洲精品久久久久久婷婷小说| 男人舔女人的私密视频| 亚洲精品国产一区二区精华液| 亚洲久久久国产精品| 99久久99久久久精品蜜桃| 久久久久久久久久久久大奶| 中文字幕制服av| 久久久久国产一级毛片高清牌| 欧美激情 高清一区二区三区| 精品少妇内射三级| 老鸭窝网址在线观看| 欧美日韩亚洲高清精品| 啦啦啦在线观看免费高清www| 高清视频免费观看一区二区| 婷婷色综合www| 只有这里有精品99| 精品亚洲成国产av| av天堂久久9| 性高湖久久久久久久久免费观看| 曰老女人黄片| 精品国产乱码久久久久久小说| 欧美在线黄色| 免费在线观看日本一区| 男女免费视频国产| 另类精品久久| 国产成人欧美在线观看 | 99九九在线精品视频| 日韩av不卡免费在线播放| 18禁观看日本| 欧美中文综合在线视频| 日韩熟女老妇一区二区性免费视频| 亚洲精品日本国产第一区| 成人手机av| 精品亚洲成国产av| 国产精品一区二区免费欧美 | 爱豆传媒免费全集在线观看| 亚洲av美国av| av不卡在线播放| 午夜av观看不卡| 国产精品麻豆人妻色哟哟久久| 9色porny在线观看| 肉色欧美久久久久久久蜜桃| 国产无遮挡羞羞视频在线观看| 国产成人免费观看mmmm| 午夜福利,免费看| 无限看片的www在线观看| 欧美亚洲 丝袜 人妻 在线| videos熟女内射| 国产精品一区二区在线不卡| 国产高清不卡午夜福利| 麻豆乱淫一区二区| 日韩大片免费观看网站| 亚洲久久久国产精品| 久久精品亚洲熟妇少妇任你| 欧美人与性动交α欧美精品济南到| 国产又爽黄色视频| 亚洲,一卡二卡三卡| 亚洲黑人精品在线| 国产深夜福利视频在线观看| 国产高清国产精品国产三级| 精品国产乱码久久久久久男人| 亚洲精品在线美女| 亚洲综合色网址| 久久鲁丝午夜福利片| 亚洲欧美中文字幕日韩二区| 少妇人妻 视频| 午夜免费男女啪啪视频观看| 国产免费一区二区三区四区乱码| 汤姆久久久久久久影院中文字幕| 亚洲成人免费av在线播放| 丰满迷人的少妇在线观看| 国产亚洲精品久久久久5区| 国产一级毛片在线| 国产激情久久老熟女| 亚洲av日韩精品久久久久久密 | 欧美人与性动交α欧美精品济南到| 日韩 欧美 亚洲 中文字幕| 久久 成人 亚洲| 老鸭窝网址在线观看| 乱人伦中国视频| 在线观看国产h片| 中文字幕人妻丝袜制服| 欧美成人精品欧美一级黄| 国产成人欧美在线观看 | 日本欧美国产在线视频| 别揉我奶头~嗯~啊~动态视频 | 视频在线观看一区二区三区| av在线app专区| 50天的宝宝边吃奶边哭怎么回事| 国产女主播在线喷水免费视频网站| 男人舔女人的私密视频| 国产高清国产精品国产三级| 国产伦人伦偷精品视频| 亚洲国产成人一精品久久久| 精品一区在线观看国产| 中文字幕另类日韩欧美亚洲嫩草| 亚洲成人免费电影在线观看 | 久久久欧美国产精品| av电影中文网址| 日本av免费视频播放| 国产真人三级小视频在线观看| 国产成人a∨麻豆精品| 久久精品成人免费网站| 亚洲国产毛片av蜜桃av| 久久 成人 亚洲| cao死你这个sao货| 亚洲国产精品国产精品| 黄色a级毛片大全视频| 欧美精品人与动牲交sv欧美| 极品少妇高潮喷水抽搐| 黄片播放在线免费| 男的添女的下面高潮视频| 亚洲av成人不卡在线观看播放网 | 另类精品久久| 又粗又硬又长又爽又黄的视频| 90打野战视频偷拍视频| 亚洲三区欧美一区| 欧美日韩亚洲国产一区二区在线观看 | 最黄视频免费看| 国产不卡av网站在线观看| 在现免费观看毛片| 99热国产这里只有精品6| 宅男免费午夜|