• <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

    亚洲av免费在线观看| 成人特级av手机在线观看| 日本黄大片高清| 联通29元200g的流量卡| 成人鲁丝片一二三区免费| 午夜亚洲福利在线播放| 国产三级在线视频| 国内精品一区二区在线观看| 在线播放国产精品三级| 国内精品宾馆在线| 丰满乱子伦码专区| 男女做爰动态图高潮gif福利片| 国产白丝娇喘喷水9色精品| 91精品国产九色| 欧美xxxx黑人xx丫x性爽| 国产一区二区三区在线臀色熟女| 黄片wwwwww| 成年女人看的毛片在线观看| 人人妻人人看人人澡| 欧美日韩亚洲国产一区二区在线观看| 俺也久久电影网| 亚洲精品国产成人久久av| 91av网一区二区| 欧美日韩瑟瑟在线播放| 日韩欧美国产在线观看| 一边摸一边抽搐一进一小说| 97碰自拍视频| 国产亚洲精品av在线| 中文字幕人妻熟人妻熟丝袜美| 中亚洲国语对白在线视频| 深夜精品福利| 伊人久久精品亚洲午夜| 国产视频内射| 久久精品综合一区二区三区| av福利片在线观看| 久久久精品大字幕| 麻豆成人午夜福利视频| 色哟哟·www| 成人毛片a级毛片在线播放| 亚洲欧美日韩高清在线视频| 亚洲在线自拍视频| 亚洲av第一区精品v没综合| 欧美又色又爽又黄视频| 久久久久久久久久成人| 特大巨黑吊av在线直播| 亚洲精品影视一区二区三区av| 欧美一区二区亚洲| 成人国产一区最新在线观看| 亚洲人成网站在线播放欧美日韩| 搞女人的毛片| 99精品在免费线老司机午夜| 高清在线国产一区| 色尼玛亚洲综合影院| 欧美成人a在线观看| 亚洲最大成人av| 国产精品国产三级国产av玫瑰| 亚洲成人精品中文字幕电影| 欧美一区二区精品小视频在线| 欧美区成人在线视频| 黄色日韩在线| 欧美3d第一页| 欧美黑人欧美精品刺激| 久久婷婷人人爽人人干人人爱| 五月伊人婷婷丁香| 99久久精品热视频| 长腿黑丝高跟| 久久热精品热| 日韩大尺度精品在线看网址| 国产欧美日韩一区二区精品| 听说在线观看完整版免费高清| 精华霜和精华液先用哪个| 国内精品美女久久久久久| 日本一本二区三区精品| 很黄的视频免费| 啪啪无遮挡十八禁网站| 亚洲av日韩精品久久久久久密| 丰满的人妻完整版| 国产精品久久久久久亚洲av鲁大| 国产高潮美女av| av天堂中文字幕网| 精品人妻偷拍中文字幕| 麻豆久久精品国产亚洲av| 精品不卡国产一区二区三区| 色综合色国产| 99热这里只有精品一区| 久久久久国内视频| 亚洲综合色惰| 亚洲一区二区三区色噜噜| 蜜桃久久精品国产亚洲av| 免费看美女性在线毛片视频| 人妻夜夜爽99麻豆av| 精品无人区乱码1区二区| 男人舔女人下体高潮全视频| 又黄又爽又刺激的免费视频.| 熟女人妻精品中文字幕| 一进一出抽搐gif免费好疼| 国产日本99.免费观看| 国产主播在线观看一区二区| 欧美日本视频| av在线老鸭窝| 亚洲精品456在线播放app | 午夜老司机福利剧场| a级毛片a级免费在线| 美女免费视频网站| ponron亚洲| 又粗又爽又猛毛片免费看| 日本成人三级电影网站| 婷婷精品国产亚洲av| 亚洲一区二区三区色噜噜| 嫩草影视91久久| 亚洲精品成人久久久久久| 日韩中字成人| 搡老岳熟女国产| 久久欧美精品欧美久久欧美| 国产不卡一卡二| 身体一侧抽搐| 动漫黄色视频在线观看| 欧美激情久久久久久爽电影| 久久精品国产亚洲av涩爱 | 亚洲专区中文字幕在线| 91狼人影院| 亚洲图色成人| 乱系列少妇在线播放| 给我免费播放毛片高清在线观看| 可以在线观看毛片的网站| 国产人妻一区二区三区在| 国产视频内射| 又爽又黄无遮挡网站| 99久久精品一区二区三区| 亚洲图色成人| 免费看日本二区| 亚洲人成伊人成综合网2020| 亚洲欧美日韩高清专用| 中国美女看黄片| 欧美zozozo另类| 69av精品久久久久久| 51国产日韩欧美| aaaaa片日本免费| 国产中年淑女户外野战色| 欧美黑人巨大hd| 99久久精品国产国产毛片| 欧美色视频一区免费| 99热这里只有是精品在线观看| 国产v大片淫在线免费观看| 精品久久久噜噜| 免费黄网站久久成人精品| 精品一区二区三区视频在线| 日韩亚洲欧美综合| 久久午夜亚洲精品久久| 成年人黄色毛片网站| 久久久色成人| 夜夜夜夜夜久久久久| 99热6这里只有精品| 一本一本综合久久| 国产亚洲av嫩草精品影院| 久久久久国内视频| 99久久精品国产国产毛片| 夜夜爽天天搞| 最近视频中文字幕2019在线8| avwww免费| 亚洲 国产 在线| .国产精品久久| 日本熟妇午夜| 色5月婷婷丁香| 欧美日本视频| 两性午夜刺激爽爽歪歪视频在线观看| 精品久久久久久久人妻蜜臀av| 99精品在免费线老司机午夜| 成人美女网站在线观看视频| 国产又黄又爽又无遮挡在线| av专区在线播放| 一本精品99久久精品77| 级片在线观看| 联通29元200g的流量卡| 亚洲av中文字字幕乱码综合| 两性午夜刺激爽爽歪歪视频在线观看| 色综合色国产| 精品久久久久久久人妻蜜臀av| 中文字幕人妻熟人妻熟丝袜美| 国产色婷婷99| 国产探花在线观看一区二区| 色综合婷婷激情| 亚洲精华国产精华液的使用体验 | 国产精品亚洲一级av第二区| av天堂中文字幕网| 在线a可以看的网站| 1000部很黄的大片| 小说图片视频综合网站| 欧美高清成人免费视频www| 日韩大尺度精品在线看网址| 日日摸夜夜添夜夜添小说| 日本一本二区三区精品| 国产 一区 欧美 日韩| 99久久成人亚洲精品观看| 国内精品久久久久精免费| 午夜久久久久精精品| 中文字幕久久专区| 亚洲av中文字字幕乱码综合| 在线播放无遮挡| 99久久九九国产精品国产免费| 亚洲经典国产精华液单| 人妻丰满熟妇av一区二区三区| 成人综合一区亚洲| 中文字幕av在线有码专区| 在线观看66精品国产| 亚洲真实伦在线观看| 国产亚洲91精品色在线| 制服丝袜大香蕉在线| 国产男人的电影天堂91| 在线观看美女被高潮喷水网站| 亚洲五月天丁香| 国产亚洲精品久久久com| 老司机福利观看| 午夜免费男女啪啪视频观看 | 国产aⅴ精品一区二区三区波| 久久久久久久精品吃奶| 女人十人毛片免费观看3o分钟| 我的女老师完整版在线观看| 成人特级黄色片久久久久久久| 免费人成在线观看视频色| 午夜久久久久精精品| 又爽又黄无遮挡网站| 99精品久久久久人妻精品| 免费看美女性在线毛片视频| 久久这里只有精品中国| 人妻制服诱惑在线中文字幕| 久久精品国产鲁丝片午夜精品 | 女的被弄到高潮叫床怎么办 | 老司机福利观看| 亚洲av免费高清在线观看| 51国产日韩欧美| 国产精品伦人一区二区| 啦啦啦韩国在线观看视频| 成人一区二区视频在线观看| 国产伦在线观看视频一区| 欧美一区二区国产精品久久精品| 黄片wwwwww| 日韩欧美一区二区三区在线观看| 波多野结衣高清作品| 色综合站精品国产| 黄色视频,在线免费观看| 国产精品久久久久久av不卡| 久久国内精品自在自线图片| 很黄的视频免费| 欧美高清性xxxxhd video| 日韩大尺度精品在线看网址| 搡老熟女国产l中国老女人| 国产精品国产高清国产av| 1000部很黄的大片| 嫩草影视91久久| 亚洲美女黄片视频| 久久久成人免费电影| 亚洲av日韩精品久久久久久密| 国产精品精品国产色婷婷| 亚洲黑人精品在线| 亚洲欧美激情综合另类| 亚洲精品456在线播放app | 中文字幕av成人在线电影| 搡老岳熟女国产| 日韩欧美三级三区| 九九热线精品视视频播放| 久久国内精品自在自线图片| 热99re8久久精品国产| 成人亚洲精品av一区二区| 午夜精品一区二区三区免费看| 欧美人与善性xxx| 搞女人的毛片| 免费观看在线日韩| 久久人人精品亚洲av| 欧美一区二区国产精品久久精品| 天美传媒精品一区二区| 能在线免费观看的黄片| 国产大屁股一区二区在线视频| 欧美高清成人免费视频www| 午夜福利在线观看免费完整高清在 | 国产乱人伦免费视频| 国产精品亚洲一级av第二区| 自拍偷自拍亚洲精品老妇| 精品午夜福利在线看| videossex国产| 国产视频内射| 中文字幕高清在线视频| 日本黄大片高清| 男女下面进入的视频免费午夜| 少妇裸体淫交视频免费看高清| 亚洲四区av| 人妻久久中文字幕网| 色综合婷婷激情| 夜夜看夜夜爽夜夜摸| 91麻豆精品激情在线观看国产| 制服丝袜大香蕉在线| 成人性生交大片免费视频hd| 国产色爽女视频免费观看| 欧美区成人在线视频| 国产精品福利在线免费观看| 淫秽高清视频在线观看| 国产欧美日韩一区二区精品| 色尼玛亚洲综合影院| 国产精品亚洲美女久久久| av在线观看视频网站免费| 91麻豆精品激情在线观看国产| 搡老妇女老女人老熟妇| 日本三级黄在线观看| 可以在线观看毛片的网站| 欧美日韩瑟瑟在线播放| 亚洲成人中文字幕在线播放| 欧美性猛交黑人性爽| 男插女下体视频免费在线播放| 极品教师在线视频| 成人综合一区亚洲| 中文亚洲av片在线观看爽| 国产伦精品一区二区三区视频9| 亚洲熟妇熟女久久| 午夜福利视频1000在线观看| 又爽又黄a免费视频| 网址你懂的国产日韩在线| 精品人妻视频免费看| 欧美在线一区亚洲| 淫妇啪啪啪对白视频| 国产精品一区二区三区四区免费观看 | 日本黄色片子视频| 中文在线观看免费www的网站| 欧美区成人在线视频| 欧美日本亚洲视频在线播放| 极品教师在线免费播放| 久久99热6这里只有精品| 中文资源天堂在线| 麻豆成人av在线观看| 久久婷婷人人爽人人干人人爱| 最近视频中文字幕2019在线8| 别揉我奶头~嗯~啊~动态视频| 级片在线观看| 黄色日韩在线| 美女免费视频网站| 少妇裸体淫交视频免费看高清| 欧美xxxx性猛交bbbb| 成人二区视频| 欧美高清成人免费视频www| 大型黄色视频在线免费观看| 精品久久久久久久久av| 久久午夜亚洲精品久久| 色尼玛亚洲综合影院| 国产伦在线观看视频一区| 日日啪夜夜撸| 有码 亚洲区| 少妇人妻精品综合一区二区 | 日本成人三级电影网站| 男女啪啪激烈高潮av片| 久久久精品欧美日韩精品| 美女高潮的动态| 看片在线看免费视频| 国产极品精品免费视频能看的| 男女做爰动态图高潮gif福利片| 欧美zozozo另类| 国产精品伦人一区二区| 亚洲无线观看免费| 午夜免费成人在线视频| 色综合站精品国产| 久久久国产成人免费| 国产色婷婷99| 麻豆久久精品国产亚洲av| 免费不卡的大黄色大毛片视频在线观看 | 最近最新中文字幕大全电影3| 淫妇啪啪啪对白视频| 成人欧美大片| 成人一区二区视频在线观看| 99精品在免费线老司机午夜| 波多野结衣巨乳人妻| 午夜免费男女啪啪视频观看 | 亚洲欧美日韩卡通动漫| 国产精品久久视频播放| 午夜老司机福利剧场| 99久久无色码亚洲精品果冻| 99riav亚洲国产免费| 亚洲va在线va天堂va国产| 97热精品久久久久久| 久久久国产成人免费| 97超级碰碰碰精品色视频在线观看| 精品久久久噜噜| 国产精品一区二区性色av| 天堂动漫精品| 在线免费十八禁| 亚洲在线自拍视频| 999久久久精品免费观看国产| 亚洲av美国av| 最近在线观看免费完整版| 欧美一级a爱片免费观看看| 亚洲电影在线观看av| 天堂动漫精品| 精品一区二区三区视频在线| 可以在线观看毛片的网站| 亚洲自拍偷在线| 舔av片在线| 人妻久久中文字幕网| 欧美日韩瑟瑟在线播放| 在线免费观看的www视频| 91久久精品国产一区二区成人| 91精品国产九色| 久久香蕉精品热| 长腿黑丝高跟| 中文字幕精品亚洲无线码一区| 不卡视频在线观看欧美| 国产主播在线观看一区二区| 别揉我奶头~嗯~啊~动态视频| 看片在线看免费视频| 免费看a级黄色片| 69av精品久久久久久| 真实男女啪啪啪动态图| 悠悠久久av| 观看美女的网站| 免费黄网站久久成人精品| 欧美三级亚洲精品| 久久精品影院6| 在线免费观看不下载黄p国产 | 草草在线视频免费看| 99riav亚洲国产免费| 直男gayav资源| 一区二区三区激情视频| 九色成人免费人妻av| 亚洲成人久久爱视频| 亚洲成av人片在线播放无| 村上凉子中文字幕在线| 午夜免费激情av| 精品久久久久久久人妻蜜臀av| 2021天堂中文幕一二区在线观| 男人的好看免费观看在线视频| 亚洲欧美日韩高清在线视频| 久久精品人妻少妇| 日韩av在线大香蕉| 久久久精品大字幕| 成人三级黄色视频| 成人无遮挡网站| 天堂网av新在线| 老熟妇仑乱视频hdxx| 国产成人福利小说| 97超级碰碰碰精品色视频在线观看| 女同久久另类99精品国产91| 嫩草影院精品99| 亚洲成人久久爱视频| 午夜视频国产福利| 18禁黄网站禁片免费观看直播| 一本久久中文字幕| 国产伦在线观看视频一区| 日本-黄色视频高清免费观看| 国产一区二区在线av高清观看| 精品一区二区三区视频在线| 一本精品99久久精品77| 亚洲精品亚洲一区二区| 国产成年人精品一区二区| 99热6这里只有精品| 国产精品乱码一区二三区的特点| 精品国产三级普通话版| 午夜精品在线福利| 成人精品一区二区免费| 3wmmmm亚洲av在线观看| 国内精品宾馆在线| 亚洲成人免费电影在线观看| 亚洲av熟女| 国产伦一二天堂av在线观看| 哪里可以看免费的av片| 欧美成人免费av一区二区三区| 国产高潮美女av| 亚洲av日韩精品久久久久久密| 日韩欧美在线二视频| 久久热精品热| 一个人观看的视频www高清免费观看| 直男gayav资源| 日韩大尺度精品在线看网址| av在线观看视频网站免费| 久久天躁狠狠躁夜夜2o2o| 国产欧美日韩精品一区二区| 日韩中字成人| 午夜精品久久久久久毛片777| 老司机深夜福利视频在线观看| 亚洲av美国av| 女人被狂操c到高潮| 99久久九九国产精品国产免费| 长腿黑丝高跟| 亚洲欧美精品综合久久99| 国产熟女欧美一区二区| 亚洲性夜色夜夜综合| 久久人妻av系列| 国产午夜精品久久久久久一区二区三区 | 亚洲av电影不卡..在线观看| 国产免费一级a男人的天堂| 超碰av人人做人人爽久久| 级片在线观看| 2021天堂中文幕一二区在线观| 午夜老司机福利剧场| 亚洲性久久影院| 精品久久久久久久久av| 国产精品伦人一区二区| av在线蜜桃| 亚洲av第一区精品v没综合| 真实男女啪啪啪动态图| 国产一区二区三区av在线 | 黄色配什么色好看| 国产成人福利小说| 亚洲第一电影网av| 亚洲精品亚洲一区二区| 男插女下体视频免费在线播放| 尾随美女入室| 色吧在线观看| 亚洲精品456在线播放app | 男女视频在线观看网站免费| 赤兔流量卡办理| 国产精品免费一区二区三区在线| 桃色一区二区三区在线观看| 欧美日韩黄片免| 久久久久久久精品吃奶| 成人精品一区二区免费| 色吧在线观看| 国产精品久久久久久久久免| 91久久精品电影网| 中亚洲国语对白在线视频| 真人做人爱边吃奶动态| 国产精品人妻久久久影院| 日本在线视频免费播放| 搡老岳熟女国产| 99久国产av精品| 人人妻人人澡欧美一区二区| 亚洲国产色片| 欧美三级亚洲精品| 99久久无色码亚洲精品果冻| 人妻制服诱惑在线中文字幕| 免费观看人在逋| 色5月婷婷丁香| 丰满乱子伦码专区| 九色国产91popny在线| 欧美bdsm另类| 久久精品国产99精品国产亚洲性色| 在线观看av片永久免费下载| 久久久久国内视频| 日韩精品中文字幕看吧| 亚洲性夜色夜夜综合| 欧美潮喷喷水| 亚洲18禁久久av| 老女人水多毛片| 日韩精品有码人妻一区| 成人二区视频| 嫩草影视91久久| 在线免费十八禁| av.在线天堂| 国产成人影院久久av| 欧美不卡视频在线免费观看| 我要看日韩黄色一级片| 99久久精品一区二区三区| 97超视频在线观看视频| 国产精品伦人一区二区| 国产伦一二天堂av在线观看| 国产男人的电影天堂91| 国产精品一区二区免费欧美| 日韩精品有码人妻一区| 国产v大片淫在线免费观看| 欧美区成人在线视频| 日本 欧美在线| 日韩精品中文字幕看吧| 亚洲av一区综合| 国产亚洲精品综合一区在线观看| 综合色av麻豆| 午夜免费男女啪啪视频观看 | 色吧在线观看| 欧美国产日韩亚洲一区| 美女大奶头视频| 69av精品久久久久久| 亚洲图色成人| 无遮挡黄片免费观看| 51国产日韩欧美| 大型黄色视频在线免费观看| 99精品久久久久人妻精品| 永久网站在线| 嫩草影院入口| 亚洲色图av天堂| 久久6这里有精品| 一个人看视频在线观看www免费| 免费看美女性在线毛片视频| 国产免费男女视频| 欧美激情在线99| 成人特级黄色片久久久久久久| 国产欧美日韩一区二区精品| 99久久久亚洲精品蜜臀av| 日本欧美国产在线视频| 999久久久精品免费观看国产| 久久午夜亚洲精品久久| 99久久久亚洲精品蜜臀av| 久久九九热精品免费| 精品一区二区三区人妻视频| 99在线视频只有这里精品首页| 三级国产精品欧美在线观看| 久99久视频精品免费| 又粗又爽又猛毛片免费看| aaaaa片日本免费| 欧美又色又爽又黄视频| 欧美3d第一页| 午夜福利欧美成人| 亚洲美女搞黄在线观看 | 成年免费大片在线观看| 国产精品无大码| 如何舔出高潮| 国产主播在线观看一区二区| 欧美国产日韩亚洲一区| 99热网站在线观看| 国产色爽女视频免费观看| 深夜精品福利| 啦啦啦观看免费观看视频高清| 久久精品综合一区二区三区| 人人妻人人看人人澡| 在线观看美女被高潮喷水网站| 少妇熟女aⅴ在线视频| 欧美zozozo另类| 国产精品1区2区在线观看.| 亚洲黑人精品在线| 变态另类成人亚洲欧美熟女| 毛片一级片免费看久久久久 | 最好的美女福利视频网|