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

    久久久久久久久久久丰满| 18禁在线播放成人免费| 国产 一区 欧美 日韩| 人人妻人人添人人爽欧美一区卜 | 夜夜爽夜夜爽视频| 久久久久精品性色| 人妻少妇偷人精品九色| 亚洲精品国产av成人精品| 久久国产亚洲av麻豆专区| 校园人妻丝袜中文字幕| 国产高清有码在线观看视频| 天美传媒精品一区二区| 在线观看美女被高潮喷水网站| 亚洲精品国产成人久久av| 少妇人妻久久综合中文| 男人舔奶头视频| 这个男人来自地球电影免费观看 | 成人亚洲精品一区在线观看 | av女优亚洲男人天堂| 国产高清有码在线观看视频| 菩萨蛮人人尽说江南好唐韦庄| 22中文网久久字幕| 日本免费在线观看一区| 国产精品国产av在线观看| 国产爽快片一区二区三区| 日韩亚洲欧美综合| 亚洲婷婷狠狠爱综合网| 欧美一区二区亚洲| 亚洲精品色激情综合| 最黄视频免费看| 成人毛片a级毛片在线播放| 人人妻人人澡人人爽人人夜夜| 联通29元200g的流量卡| 日本黄色片子视频| 国产精品嫩草影院av在线观看| 少妇裸体淫交视频免费看高清| 自拍欧美九色日韩亚洲蝌蚪91 | 美女视频免费永久观看网站| 最近手机中文字幕大全| 中文字幕av成人在线电影| 黄色配什么色好看| 精品人妻熟女av久视频| 午夜激情福利司机影院| 能在线免费看毛片的网站| 中国国产av一级| 国产精品国产三级国产av玫瑰| 亚洲成人中文字幕在线播放| 久久精品人妻少妇| 精品熟女少妇av免费看| 久久国产精品大桥未久av | 国产在视频线精品| 女人十人毛片免费观看3o分钟| 国产视频首页在线观看| 18禁裸乳无遮挡动漫免费视频| 成人国产麻豆网| 国产熟女欧美一区二区| 一个人看的www免费观看视频| 日韩成人伦理影院| av在线老鸭窝| 丝瓜视频免费看黄片| 天堂俺去俺来也www色官网| 久久国内精品自在自线图片| 大片免费播放器 马上看| 观看美女的网站| 日韩强制内射视频| 欧美精品国产亚洲| 亚洲av男天堂| 黄片无遮挡物在线观看| 亚洲国产高清在线一区二区三| 一个人看视频在线观看www免费| 91精品一卡2卡3卡4卡| 黄色欧美视频在线观看| 乱码一卡2卡4卡精品| 欧美极品一区二区三区四区| 观看免费一级毛片| 中文字幕亚洲精品专区| 干丝袜人妻中文字幕| 久久久久精品性色| 精品一区二区三区视频在线| 99久国产av精品国产电影| 高清av免费在线| 日韩大片免费观看网站| 少妇人妻 视频| 午夜福利高清视频| 肉色欧美久久久久久久蜜桃| av福利片在线观看| 联通29元200g的流量卡| 国产色爽女视频免费观看| 极品教师在线视频| 美女视频免费永久观看网站| 亚洲av综合色区一区| 国产精品99久久99久久久不卡 | 国产国拍精品亚洲av在线观看| 欧美最新免费一区二区三区| 国产一区二区三区综合在线观看 | 久久99热6这里只有精品| 国产精品久久久久久精品电影小说 | 一个人免费看片子| 亚洲美女搞黄在线观看| 国产精品99久久99久久久不卡 | 国产高清有码在线观看视频| 日日摸夜夜添夜夜爱| 卡戴珊不雅视频在线播放| 天美传媒精品一区二区| 成人18禁高潮啪啪吃奶动态图 | 爱豆传媒免费全集在线观看| 少妇丰满av| 国产男人的电影天堂91| 久久人人爽av亚洲精品天堂 | 国产精品人妻久久久久久| 久久婷婷青草| 国产淫语在线视频| 精品一区二区三卡| 成人综合一区亚洲| av专区在线播放| 日日啪夜夜爽| 一级av片app| 男人添女人高潮全过程视频| 夫妻性生交免费视频一级片| 五月伊人婷婷丁香| 免费不卡的大黄色大毛片视频在线观看| 亚洲婷婷狠狠爱综合网| 亚洲精品中文字幕在线视频 | 男女下面进入的视频免费午夜| 国产亚洲午夜精品一区二区久久| 国产精品久久久久成人av| 色视频www国产| 精品视频人人做人人爽| 国产高清不卡午夜福利| 国产精品偷伦视频观看了| 深夜a级毛片| 久久毛片免费看一区二区三区| 日韩视频在线欧美| 老熟女久久久| 精品99又大又爽又粗少妇毛片| 欧美成人午夜免费资源| videos熟女内射| 乱码一卡2卡4卡精品| 联通29元200g的流量卡| 亚洲欧美成人精品一区二区| 国产伦理片在线播放av一区| 在线播放无遮挡| 亚洲欧美成人精品一区二区| 久久av网站| 免费看光身美女| 岛国毛片在线播放| .国产精品久久| 黑丝袜美女国产一区| 国产有黄有色有爽视频| 永久免费av网站大全| 午夜福利在线在线| 哪个播放器可以免费观看大片| 久久精品国产自在天天线| 视频区图区小说| 七月丁香在线播放| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲精品乱码久久久久久按摩| 日韩人妻高清精品专区| 天天躁夜夜躁狠狠久久av| 欧美亚洲 丝袜 人妻 在线| 黄片wwwwww| 99热网站在线观看| 国产男女内射视频| 九九久久精品国产亚洲av麻豆| 搡女人真爽免费视频火全软件| 色婷婷av一区二区三区视频| 最近最新中文字幕免费大全7| av又黄又爽大尺度在线免费看| 深夜a级毛片| 成人国产麻豆网| 欧美精品一区二区免费开放| 日韩伦理黄色片| 日本wwww免费看| 色吧在线观看| 黑人高潮一二区| 午夜激情久久久久久久| 久久久久视频综合| 麻豆国产97在线/欧美| 亚洲精品国产色婷婷电影| 亚洲国产精品一区三区| av又黄又爽大尺度在线免费看| 中文精品一卡2卡3卡4更新| av免费观看日本| 国产黄频视频在线观看| 色网站视频免费| www.av在线官网国产| www.av在线官网国产| 夜夜看夜夜爽夜夜摸| av网站免费在线观看视频| 高清黄色对白视频在线免费看 | 久久精品国产a三级三级三级| 涩涩av久久男人的天堂| 伦理电影大哥的女人| 大陆偷拍与自拍| 亚洲精品自拍成人| 国产精品一区二区性色av| 久久亚洲国产成人精品v| 黄色配什么色好看| 大陆偷拍与自拍| 欧美精品一区二区大全| av免费观看日本| 国产精品国产三级国产av玫瑰| 激情五月婷婷亚洲| 久久精品久久精品一区二区三区| 一级毛片我不卡| 国产美女午夜福利| 免费在线观看成人毛片| 久久久久国产精品人妻一区二区| 性色av一级| 亚洲av在线观看美女高潮| 亚洲精品中文字幕在线视频 | 中文字幕制服av| 国产成人精品一,二区| 国产av精品麻豆| 国产精品秋霞免费鲁丝片| 女人十人毛片免费观看3o分钟| 国产一级毛片在线| 日日摸夜夜添夜夜添av毛片| 日韩av不卡免费在线播放| 国产国拍精品亚洲av在线观看| 99久国产av精品国产电影| 精品亚洲成国产av| 插逼视频在线观看| 国产真实伦视频高清在线观看| 高清av免费在线| 亚洲欧美中文字幕日韩二区| 亚洲av成人精品一区久久| 老师上课跳d突然被开到最大视频| 亚洲成色77777| 国产免费一级a男人的天堂| 最近最新中文字幕免费大全7| 国产在视频线精品| av在线老鸭窝| 伊人久久精品亚洲午夜| av在线蜜桃| 午夜激情久久久久久久| av国产久精品久网站免费入址| 国产精品蜜桃在线观看| 亚洲第一av免费看| 人妻制服诱惑在线中文字幕| 中文天堂在线官网| 成人黄色视频免费在线看| 99九九线精品视频在线观看视频| 午夜福利影视在线免费观看| 日本黄色片子视频| 成人美女网站在线观看视频| 日本黄色日本黄色录像| 男女免费视频国产| 99热国产这里只有精品6| 九色成人免费人妻av| 天美传媒精品一区二区| 日本欧美国产在线视频| 国产精品伦人一区二区| 乱码一卡2卡4卡精品| 在线精品无人区一区二区三 | 亚洲av成人精品一二三区| 啦啦啦在线观看免费高清www| 91在线精品国自产拍蜜月| 亚洲精品国产色婷婷电影| 亚洲成人手机| 色视频www国产| 国产免费又黄又爽又色| 亚洲国产日韩一区二区| 边亲边吃奶的免费视频| 在线观看一区二区三区| 高清av免费在线| 人人妻人人添人人爽欧美一区卜 | 精品亚洲成国产av| 色5月婷婷丁香| 国产成人精品久久久久久| 成人黄色视频免费在线看| 亚洲,一卡二卡三卡| 亚洲欧美中文字幕日韩二区| 丝瓜视频免费看黄片| 日韩成人伦理影院| 人妻夜夜爽99麻豆av| 精品国产乱码久久久久久小说| 啦啦啦在线观看免费高清www| 日韩成人av中文字幕在线观看| 在线 av 中文字幕| 18禁动态无遮挡网站| 妹子高潮喷水视频| 内地一区二区视频在线| a级一级毛片免费在线观看| 人妻制服诱惑在线中文字幕| 亚洲精品自拍成人| av又黄又爽大尺度在线免费看| 一二三四中文在线观看免费高清| 亚洲人与动物交配视频| av黄色大香蕉| 日本黄色日本黄色录像| 国产精品秋霞免费鲁丝片| 视频中文字幕在线观看| 国产人妻一区二区三区在| 精品亚洲成a人片在线观看 | 国产精品99久久久久久久久| 只有这里有精品99| 精品少妇久久久久久888优播| 蜜桃亚洲精品一区二区三区| av视频免费观看在线观看| h视频一区二区三区| 在线观看一区二区三区激情| 国产在线免费精品| 国产精品福利在线免费观看| 一级a做视频免费观看| 最近手机中文字幕大全| 精品99又大又爽又粗少妇毛片| 尤物成人国产欧美一区二区三区| 小蜜桃在线观看免费完整版高清| 欧美三级亚洲精品| 少妇人妻久久综合中文| 亚洲欧美日韩卡通动漫| 国产亚洲一区二区精品| 大又大粗又爽又黄少妇毛片口| 亚洲人成网站高清观看| 老司机影院毛片| 日韩不卡一区二区三区视频在线| 日韩免费高清中文字幕av| 国产欧美日韩精品一区二区| 亚洲精品,欧美精品| 天天躁日日操中文字幕| 久久久久久久亚洲中文字幕| 男女边摸边吃奶| 亚洲精品久久午夜乱码| 免费黄频网站在线观看国产| 亚洲精品自拍成人| 国产成人免费无遮挡视频| 亚洲av成人精品一区久久| 久久久午夜欧美精品| 91精品伊人久久大香线蕉| 亚洲av成人精品一区久久| 国产亚洲5aaaaa淫片| 日日摸夜夜添夜夜爱| 欧美97在线视频| 久久国产乱子免费精品| 亚洲国产最新在线播放| 少妇的逼水好多| 精品久久久久久久久亚洲| 最后的刺客免费高清国语| 免费看日本二区| 日韩成人伦理影院| 亚洲激情五月婷婷啪啪| 黑人高潮一二区| 国产精品久久久久久av不卡| 久久精品夜色国产| 日韩电影二区| 国产精品一区二区性色av| 亚洲精品乱久久久久久| 王馨瑶露胸无遮挡在线观看| 亚洲av综合色区一区| 18+在线观看网站| 亚洲精品aⅴ在线观看| 国产精品秋霞免费鲁丝片| av国产精品久久久久影院| 免费观看的影片在线观看| 国产大屁股一区二区在线视频| 国产在线免费精品| 精品人妻视频免费看| 激情五月婷婷亚洲| 下体分泌物呈黄色| 黄色欧美视频在线观看| 亚洲欧美一区二区三区黑人 | 纯流量卡能插随身wifi吗| 亚洲久久久国产精品| 国产有黄有色有爽视频| 草草在线视频免费看| 97在线视频观看| av不卡在线播放| 亚洲国产毛片av蜜桃av| 人体艺术视频欧美日本| 少妇丰满av| 黑人高潮一二区| 18禁裸乳无遮挡免费网站照片| 国产精品麻豆人妻色哟哟久久| 免费看光身美女| 国产一区二区三区av在线| 久久久久人妻精品一区果冻| 熟女人妻精品中文字幕| 亚洲三级黄色毛片| av网站免费在线观看视频| 久久久国产一区二区| 欧美日韩国产mv在线观看视频 | 欧美97在线视频| 肉色欧美久久久久久久蜜桃| 精品久久久久久电影网| 国产成人免费无遮挡视频| 久久久久久久国产电影| 人人妻人人爽人人添夜夜欢视频 | 五月开心婷婷网| 免费人妻精品一区二区三区视频| 777米奇影视久久| 一本一本综合久久| 国产大屁股一区二区在线视频| 人人妻人人澡人人爽人人夜夜| 精品国产乱码久久久久久小说| 午夜免费鲁丝| 性色av一级| 新久久久久国产一级毛片| 一个人看的www免费观看视频| 观看av在线不卡| 美女内射精品一级片tv| 性色av一级| 在线精品无人区一区二区三 | 边亲边吃奶的免费视频| 亚洲自偷自拍三级| 亚洲国产色片| 精品一品国产午夜福利视频| 欧美三级亚洲精品| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲欧美日韩另类电影网站 | 成年av动漫网址| 内射极品少妇av片p| 国产精品国产三级专区第一集| 少妇人妻 视频| 精品久久久久久久末码| 久久精品国产鲁丝片午夜精品| 精品视频人人做人人爽| 看十八女毛片水多多多| 久久久久久久久大av| 激情 狠狠 欧美| 51国产日韩欧美| av视频免费观看在线观看| 一区二区三区精品91| 天堂8中文在线网| 日韩一区二区三区影片| 直男gayav资源| av福利片在线观看| 亚洲伊人久久精品综合| 91精品伊人久久大香线蕉| 99热这里只有精品一区| 九草在线视频观看| 美女主播在线视频| 91久久精品国产一区二区成人| 中文字幕亚洲精品专区| 人人妻人人澡人人爽人人夜夜| 精品国产乱码久久久久久小说| av在线播放精品| 晚上一个人看的免费电影| h日本视频在线播放| 亚洲av成人精品一二三区| 18禁裸乳无遮挡免费网站照片| 亚洲欧美日韩另类电影网站 | 狂野欧美白嫩少妇大欣赏| 亚洲av电影在线观看一区二区三区| 51国产日韩欧美| 大话2 男鬼变身卡| 一级毛片我不卡| 伦理电影大哥的女人| 观看美女的网站| 亚洲人成网站在线观看播放| 777米奇影视久久| 日日啪夜夜爽| 视频区图区小说| 91精品国产九色| 蜜臀久久99精品久久宅男| 人妻少妇偷人精品九色| 久久人妻熟女aⅴ| 毛片一级片免费看久久久久| 国产伦在线观看视频一区| 丰满乱子伦码专区| 男女无遮挡免费网站观看| 各种免费的搞黄视频| 亚洲图色成人| 99热国产这里只有精品6| 免费看日本二区| 国产午夜精品久久久久久一区二区三区| 少妇人妻精品综合一区二区| 22中文网久久字幕| 丝袜脚勾引网站| 天堂中文最新版在线下载| 大片免费播放器 马上看| 黄色配什么色好看| 午夜免费鲁丝| 亚洲国产精品999| 国产亚洲最大av| 日韩视频在线欧美| 欧美精品国产亚洲| 新久久久久国产一级毛片| 亚洲成人中文字幕在线播放| 韩国高清视频一区二区三区| 亚洲一区二区三区欧美精品| 亚洲国产精品成人久久小说| 插逼视频在线观看| 久久午夜福利片| 成人亚洲精品一区在线观看 | 国产 精品1| 三级国产精品欧美在线观看| 国内精品宾馆在线| 午夜福利在线观看免费完整高清在| 春色校园在线视频观看| 欧美xxⅹ黑人| 亚洲av免费高清在线观看| 夜夜看夜夜爽夜夜摸| 精品久久久精品久久久| 黑丝袜美女国产一区| 亚洲高清免费不卡视频| 性高湖久久久久久久久免费观看| 91在线精品国自产拍蜜月| 99国产精品免费福利视频| 国产综合精华液| 久久亚洲国产成人精品v| 我要看日韩黄色一级片| 日韩在线高清观看一区二区三区| 国产精品三级大全| 日韩电影二区| 免费久久久久久久精品成人欧美视频 | 美女内射精品一级片tv| 五月开心婷婷网| 女人十人毛片免费观看3o分钟| 看免费成人av毛片| 久久人人爽av亚洲精品天堂 | 人人妻人人添人人爽欧美一区卜 | 日韩精品有码人妻一区| 香蕉精品网在线| 日本av免费视频播放| 啦啦啦啦在线视频资源| 亚洲美女黄色视频免费看| 婷婷色av中文字幕| 亚洲精品自拍成人| 99热全是精品| 一区二区三区精品91| 国产精品一二三区在线看| 国产精品偷伦视频观看了| 老司机影院成人| 三级国产精品欧美在线观看| 18禁在线无遮挡免费观看视频| 国产男女超爽视频在线观看| 六月丁香七月| 欧美区成人在线视频| 亚洲av电影在线观看一区二区三区| 国产av一区二区精品久久 | 人妻系列 视频| 这个男人来自地球电影免费观看 | 男人舔奶头视频| 少妇 在线观看| 久久99精品国语久久久| 麻豆精品久久久久久蜜桃| 亚洲av中文av极速乱| 亚洲欧美日韩另类电影网站 | av在线播放精品| 看十八女毛片水多多多| 欧美日韩国产mv在线观看视频 | 亚洲四区av| 免费观看在线日韩| 成年人午夜在线观看视频| 久久热精品热| 超碰97精品在线观看| 天天躁夜夜躁狠狠久久av| 国产精品国产三级国产专区5o| 亚洲熟女精品中文字幕| 精品99又大又爽又粗少妇毛片| 亚洲综合色惰| 欧美一级a爱片免费观看看| 丰满乱子伦码专区| 内射极品少妇av片p| 日日摸夜夜添夜夜爱| 欧美区成人在线视频| 男人爽女人下面视频在线观看| 亚洲精品乱码久久久久久按摩| 性色av一级| 日本wwww免费看| 亚洲人成网站高清观看| 亚洲欧洲日产国产| 人妻一区二区av| 国产免费视频播放在线视频| 久久久久久久久久人人人人人人| 国语对白做爰xxxⅹ性视频网站| 夜夜看夜夜爽夜夜摸| 在线免费十八禁| 日韩av在线免费看完整版不卡| 国产亚洲一区二区精品| 26uuu在线亚洲综合色| videos熟女内射| 草草在线视频免费看| 亚洲国产毛片av蜜桃av| 夫妻午夜视频| 99re6热这里在线精品视频| 日日摸夜夜添夜夜爱| 全区人妻精品视频| 精品人妻熟女av久视频| 亚洲精品,欧美精品| 成年免费大片在线观看| av国产免费在线观看| 97在线人人人人妻| 国产亚洲欧美精品永久| 少妇的逼水好多| 日本爱情动作片www.在线观看| 国产精品一区二区在线观看99| 亚洲四区av| 亚洲精品日韩在线中文字幕| 蜜桃亚洲精品一区二区三区| 精品久久久久久电影网| 男人添女人高潮全过程视频| 成人漫画全彩无遮挡| 黄色配什么色好看| 日韩免费高清中文字幕av| 22中文网久久字幕| 美女视频免费永久观看网站| 国产伦精品一区二区三区视频9| 99热这里只有是精品50| 日韩强制内射视频| 不卡视频在线观看欧美| 性色avwww在线观看| 久久热精品热| 香蕉精品网在线| 97超视频在线观看视频| 国产精品国产av在线观看| 黄色怎么调成土黄色| 人人妻人人澡人人爽人人夜夜| 国产精品嫩草影院av在线观看| 2022亚洲国产成人精品| 免费黄网站久久成人精品| 在线观看三级黄色| 国产精品国产av在线观看| 日日摸夜夜添夜夜爱|