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

    欧美日韩中文字幕国产精品一区二区三区| 男人的好看免费观看在线视频 | 国内久久婷婷六月综合欲色啪| 熟妇人妻久久中文字幕3abv| 神马国产精品三级电影在线观看 | 18禁裸乳无遮挡免费网站照片 | 久久国产乱子伦精品免费另类| 国产成年人精品一区二区| 99国产精品一区二区蜜桃av| 精品久久久久久成人av| 欧美乱码精品一区二区三区| 成人国产一区最新在线观看| 久久香蕉国产精品| 麻豆国产av国片精品| 亚洲天堂国产精品一区在线| 亚洲人成网站在线播放欧美日韩| 日韩精品免费视频一区二区三区| 久久中文字幕人妻熟女| 首页视频小说图片口味搜索| 精品一区二区三区视频在线观看免费| 国产主播在线观看一区二区| 精品久久久久久久末码| 最新美女视频免费是黄的| 一本大道久久a久久精品| 久久欧美精品欧美久久欧美| 成在线人永久免费视频| 一进一出抽搐动态| 精品久久久久久久人妻蜜臀av| 欧美日韩黄片免| 亚洲第一青青草原| 欧美成人一区二区免费高清观看 | 成年人黄色毛片网站| 日本三级黄在线观看| www日本在线高清视频| 国产av在哪里看| 女性生殖器流出的白浆| 国产极品粉嫩免费观看在线| 麻豆一二三区av精品| 精品少妇一区二区三区视频日本电影| 国产精品免费一区二区三区在线| 好男人在线观看高清免费视频 | 一个人免费在线观看的高清视频| 亚洲成a人片在线一区二区| 精品久久久久久久人妻蜜臀av| 午夜福利高清视频| 亚洲av电影在线进入| 操出白浆在线播放| av电影中文网址| 女人高潮潮喷娇喘18禁视频| 黄频高清免费视频| 亚洲国产中文字幕在线视频| 久久香蕉精品热| 久久99热这里只有精品18| 亚洲熟妇中文字幕五十中出| 黄色视频,在线免费观看| 午夜亚洲福利在线播放| 日韩视频一区二区在线观看| videosex国产| 久久人妻福利社区极品人妻图片| 国产成年人精品一区二区| 久久久久久久午夜电影| 国产黄色小视频在线观看| 日韩高清综合在线| 午夜精品久久久久久毛片777| 三级毛片av免费| 最新在线观看一区二区三区| 熟女少妇亚洲综合色aaa.| 久久午夜综合久久蜜桃| 国产1区2区3区精品| 亚洲午夜精品一区,二区,三区| 久久久久久人人人人人| 国产一区二区三区在线臀色熟女| 欧美黄色淫秽网站| 欧美性长视频在线观看| 香蕉av资源在线| a在线观看视频网站| 中文字幕最新亚洲高清| 成人特级黄色片久久久久久久| 国产一卡二卡三卡精品| 精品熟女少妇八av免费久了| 中文资源天堂在线| 男人舔女人下体高潮全视频| 狠狠狠狠99中文字幕| 欧美成人免费av一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 亚洲美女黄片视频| 一级片免费观看大全| 国产精品 欧美亚洲| 亚洲在线自拍视频| 麻豆av在线久日| 嫩草影院精品99| 男女视频在线观看网站免费 | 欧美丝袜亚洲另类 | 麻豆成人av在线观看| 亚洲五月色婷婷综合| 波多野结衣高清无吗| 巨乳人妻的诱惑在线观看| 中文字幕人妻丝袜一区二区| 久久久国产精品麻豆| 最近最新免费中文字幕在线| 欧美日韩一级在线毛片| 国产成人欧美在线观看| 怎么达到女性高潮| 老熟妇乱子伦视频在线观看| 午夜免费成人在线视频| 欧美在线黄色| 日韩高清综合在线| 中文资源天堂在线| 精品国产美女av久久久久小说| 丝袜人妻中文字幕| 久久热在线av| 成人永久免费在线观看视频| 久久久久国产精品人妻aⅴ院| 无限看片的www在线观看| 一级毛片高清免费大全| 日日夜夜操网爽| netflix在线观看网站| 亚洲av五月六月丁香网| 免费在线观看影片大全网站| av欧美777| 长腿黑丝高跟| 久久香蕉国产精品| 欧美性猛交黑人性爽| 18禁美女被吸乳视频| 成人手机av| 成熟少妇高潮喷水视频| 日日爽夜夜爽网站| 女性被躁到高潮视频| 身体一侧抽搐| 国产精品久久久av美女十八| 精品乱码久久久久久99久播| 中亚洲国语对白在线视频| 又黄又爽又免费观看的视频| 2021天堂中文幕一二区在线观 | 亚洲av五月六月丁香网| 国产av不卡久久| 日韩有码中文字幕| 欧美 亚洲 国产 日韩一| avwww免费| 国产成人影院久久av| 亚洲在线自拍视频| 美女大奶头视频| av超薄肉色丝袜交足视频| 精品久久久久久久久久免费视频| 午夜福利18| 欧美zozozo另类| 国产av一区在线观看免费| 精品一区二区三区视频在线观看免费| 精品免费久久久久久久清纯| 亚洲在线自拍视频| 亚洲av成人不卡在线观看播放网| 中亚洲国语对白在线视频| 亚洲国产中文字幕在线视频| 好男人电影高清在线观看| 免费看美女性在线毛片视频| 91老司机精品| 在线观看一区二区三区| 在线观看免费午夜福利视频| 无限看片的www在线观看| 在线永久观看黄色视频| 这个男人来自地球电影免费观看| 久久 成人 亚洲| 91成人精品电影| 夜夜爽天天搞| 麻豆一二三区av精品| 精品国产一区二区三区四区第35| bbb黄色大片| 国产免费男女视频| 免费av毛片视频| 一进一出抽搐动态| 国产视频一区二区在线看| 女性生殖器流出的白浆| 最新美女视频免费是黄的| 97碰自拍视频| 18美女黄网站色大片免费观看| 中文字幕av电影在线播放| www日本黄色视频网| 俺也久久电影网| 黄色毛片三级朝国网站| 久久精品夜夜夜夜夜久久蜜豆 | xxxwww97欧美| 夜夜看夜夜爽夜夜摸| 听说在线观看完整版免费高清| 欧美另类亚洲清纯唯美| 精品电影一区二区在线| www.www免费av| 真人做人爱边吃奶动态| 色综合亚洲欧美另类图片| 精品熟女少妇八av免费久了| 亚洲五月色婷婷综合| 欧美性猛交黑人性爽| www.精华液| 国产成人欧美| 老鸭窝网址在线观看| 久久青草综合色| 国产在线观看jvid| 三级毛片av免费| 国产精品一区二区精品视频观看| 天天添夜夜摸| 淫妇啪啪啪对白视频| av天堂在线播放| 在线观看午夜福利视频| 国产精品乱码一区二三区的特点| 久久久久久大精品| 久久国产精品人妻蜜桃| 一区福利在线观看| 老司机午夜十八禁免费视频| 国产午夜福利久久久久久| 国产精品自产拍在线观看55亚洲| 两个人看的免费小视频| 99re在线观看精品视频| 欧美色欧美亚洲另类二区| 熟女电影av网| 男女之事视频高清在线观看| 可以在线观看毛片的网站| 丁香欧美五月| aaaaa片日本免费| 在线观看66精品国产| 久久久国产精品麻豆| xxxwww97欧美| 美国免费a级毛片| 国产亚洲av高清不卡| 午夜福利在线在线| 国产精品二区激情视频| a级毛片a级免费在线| 男人操女人黄网站| 亚洲男人的天堂狠狠| 91在线观看av| 久久精品人妻少妇| 久久中文看片网| 久久精品成人免费网站| www.精华液| 亚洲国产精品合色在线| 波多野结衣高清无吗| 91成年电影在线观看| 精品欧美国产一区二区三| 亚洲第一电影网av| 又紧又爽又黄一区二区| 美女午夜性视频免费| 国产成人精品无人区| 精品午夜福利视频在线观看一区| 国产高清视频在线播放一区| 色精品久久人妻99蜜桃| 日韩国内少妇激情av| 久久婷婷人人爽人人干人人爱| 欧美人与性动交α欧美精品济南到| 免费在线观看视频国产中文字幕亚洲| 天堂√8在线中文| 国产v大片淫在线免费观看| 动漫黄色视频在线观看| 国产精品 国内视频| 国产1区2区3区精品| 国产高清激情床上av| 18禁黄网站禁片午夜丰满| 亚洲欧美精品综合一区二区三区| 亚洲av第一区精品v没综合| 成人18禁高潮啪啪吃奶动态图| 人人妻,人人澡人人爽秒播| 久久久久国内视频| 婷婷亚洲欧美| 丝袜美腿诱惑在线| 一区二区三区国产精品乱码| 黄色女人牲交| 精品久久蜜臀av无| 日韩精品青青久久久久久| 亚洲精华国产精华精| 麻豆一二三区av精品| 国产精品免费一区二区三区在线| 午夜久久久在线观看| 亚洲 国产 在线| 日本 欧美在线| 女人爽到高潮嗷嗷叫在线视频| 久久久久久久久久黄片| 黄色丝袜av网址大全| 免费看a级黄色片| 欧美色视频一区免费| 国产伦在线观看视频一区| 免费看a级黄色片| 成人午夜高清在线视频 | 女同久久另类99精品国产91| 亚洲美女黄片视频| 亚洲七黄色美女视频| 久久人妻av系列| 午夜日韩欧美国产| 一进一出抽搐gif免费好疼| 欧美黑人欧美精品刺激| 中亚洲国语对白在线视频| 亚洲精品国产一区二区精华液| 中文亚洲av片在线观看爽| 亚洲国产欧洲综合997久久, | 免费无遮挡裸体视频| 黄色女人牲交| 麻豆成人午夜福利视频| 精品久久久久久久人妻蜜臀av| 91成人精品电影| 久久亚洲真实| 男女那种视频在线观看| 精品电影一区二区在线| 日本a在线网址| 久久午夜综合久久蜜桃| 日韩免费av在线播放| 精品免费久久久久久久清纯| 黄色丝袜av网址大全| 老熟妇仑乱视频hdxx| 中文字幕另类日韩欧美亚洲嫩草| 国产片内射在线| 午夜视频精品福利| 中文字幕另类日韩欧美亚洲嫩草| 女人爽到高潮嗷嗷叫在线视频| 国产精品日韩av在线免费观看| 精品第一国产精品| 国产一卡二卡三卡精品| 中文字幕最新亚洲高清| 国产一卡二卡三卡精品| 九色国产91popny在线| 国产av一区二区精品久久| 特大巨黑吊av在线直播 | 亚洲欧洲精品一区二区精品久久久| 精品久久久久久,| 很黄的视频免费| 国内少妇人妻偷人精品xxx网站 | 这个男人来自地球电影免费观看| 精品不卡国产一区二区三区| 亚洲五月天丁香| 日本熟妇午夜| 国产精品av久久久久免费| 九色国产91popny在线| 国产免费av片在线观看野外av| 国产亚洲欧美精品永久| 亚洲片人在线观看| 99热这里只有精品一区 | 丝袜在线中文字幕| 亚洲av第一区精品v没综合| 热99re8久久精品国产| 亚洲精品中文字幕在线视频| 久久精品成人免费网站| 国产精品98久久久久久宅男小说| 免费人成视频x8x8入口观看| 精品久久久久久成人av| 男人舔女人下体高潮全视频| 欧洲精品卡2卡3卡4卡5卡区| 国产主播在线观看一区二区| 国产片内射在线| 午夜福利一区二区在线看| 91麻豆av在线| 久久国产精品男人的天堂亚洲| 成年免费大片在线观看| 国产精品久久久久久亚洲av鲁大| 日韩一卡2卡3卡4卡2021年| 午夜福利高清视频| 精品无人区乱码1区二区| 久久精品人妻少妇| aaaaa片日本免费| 神马国产精品三级电影在线观看 | 在线观看舔阴道视频| 淫秽高清视频在线观看| 热re99久久国产66热| 色综合站精品国产| 看片在线看免费视频| 天堂√8在线中文| 精品国内亚洲2022精品成人| 国产亚洲欧美在线一区二区| 久久精品国产综合久久久| 亚洲中文av在线| 国产99久久九九免费精品| АⅤ资源中文在线天堂| 波多野结衣巨乳人妻| 黄色视频不卡| 99久久99久久久精品蜜桃| videosex国产| 亚洲一区二区三区不卡视频| 欧美日韩乱码在线| 欧美国产精品va在线观看不卡| 精品国产乱码久久久久久男人| 精品国产国语对白av| 视频区欧美日本亚洲| 波多野结衣巨乳人妻| 午夜福利免费观看在线| 久久精品国产综合久久久| 国产精品二区激情视频| 午夜成年电影在线免费观看| 国产高清激情床上av| 国产激情久久老熟女| 伊人久久大香线蕉亚洲五| 亚洲精品久久国产高清桃花| 中文亚洲av片在线观看爽| 91国产中文字幕| 欧美日韩黄片免| 久久久国产成人免费| 黄色 视频免费看| 男人操女人黄网站| 男人舔奶头视频| 久久这里只有精品19| 国产av在哪里看| 国产成人欧美在线观看| 亚洲片人在线观看| 精品第一国产精品| 亚洲中文av在线| 亚洲av成人一区二区三| 91av网站免费观看| 老熟妇乱子伦视频在线观看| 十八禁人妻一区二区| 久久久久久免费高清国产稀缺| 好男人在线观看高清免费视频 | 精品日产1卡2卡| 国产一区在线观看成人免费| 国产精品 欧美亚洲| 成人午夜高清在线视频 | 别揉我奶头~嗯~啊~动态视频| 欧美三级亚洲精品| 欧美色视频一区免费| 精品电影一区二区在线| 亚洲一区高清亚洲精品| 久久亚洲精品不卡| 国产又黄又爽又无遮挡在线| 国产欧美日韩精品亚洲av| 久久国产亚洲av麻豆专区| 国产精品自产拍在线观看55亚洲| 757午夜福利合集在线观看| 母亲3免费完整高清在线观看| 国产精品1区2区在线观看.| 精品久久久久久久末码| 麻豆久久精品国产亚洲av| 一进一出抽搐动态| 日韩成人在线观看一区二区三区| 日本三级黄在线观看| 亚洲,欧美精品.| 啦啦啦韩国在线观看视频| 国产亚洲欧美在线一区二区| 搞女人的毛片| 视频在线观看一区二区三区| 12—13女人毛片做爰片一| x7x7x7水蜜桃| 久久狼人影院| 人人澡人人妻人| 国产亚洲欧美在线一区二区| 欧美绝顶高潮抽搐喷水| 欧美黑人巨大hd| 亚洲五月色婷婷综合| 欧美色视频一区免费| 欧美又色又爽又黄视频| 校园春色视频在线观看| 18禁国产床啪视频网站| 人人澡人人妻人| 欧美性猛交黑人性爽| 两性午夜刺激爽爽歪歪视频在线观看 | av在线播放免费不卡| 日韩欧美国产在线观看| 久久久久久人人人人人| 无人区码免费观看不卡| 亚洲国产毛片av蜜桃av| svipshipincom国产片| 男女床上黄色一级片免费看| 男男h啪啪无遮挡| 亚洲av第一区精品v没综合| 欧美又色又爽又黄视频| 最新在线观看一区二区三区| 婷婷六月久久综合丁香| 成人一区二区视频在线观看| 久久中文字幕人妻熟女| 黄色成人免费大全| av福利片在线| 欧美日本亚洲视频在线播放| 欧美中文综合在线视频| 亚洲久久久国产精品| 一级毛片高清免费大全| 亚洲一码二码三码区别大吗| 亚洲精品在线美女| 国产麻豆成人av免费视频| 亚洲国产精品999在线| 亚洲国产中文字幕在线视频| 天堂√8在线中文| 美国免费a级毛片| 日本 欧美在线| 老司机靠b影院| 十分钟在线观看高清视频www| 久久精品aⅴ一区二区三区四区| 我的亚洲天堂| 国产精品精品国产色婷婷| 岛国在线观看网站| 后天国语完整版免费观看| 人人妻,人人澡人人爽秒播| 特大巨黑吊av在线直播 | 亚洲中文日韩欧美视频| 日韩av在线大香蕉| 熟女电影av网| 亚洲欧美激情综合另类| 日韩欧美一区二区三区在线观看| 精品久久久久久久久久免费视频| 叶爱在线成人免费视频播放| 欧美 亚洲 国产 日韩一| 18禁裸乳无遮挡免费网站照片 | 成人手机av| av天堂在线播放| 日本三级黄在线观看| 日韩三级视频一区二区三区| 黄色毛片三级朝国网站| cao死你这个sao货| 一进一出抽搐gif免费好疼| 久久草成人影院| 免费在线观看完整版高清| 色综合婷婷激情| 母亲3免费完整高清在线观看| 欧美亚洲日本最大视频资源| 午夜老司机福利片| av欧美777| 成人一区二区视频在线观看| 亚洲av熟女| 午夜两性在线视频| 日本一区二区免费在线视频| 国产精品香港三级国产av潘金莲| 91大片在线观看| videosex国产| av免费在线观看网站| 亚洲成人久久爱视频| av在线播放免费不卡| 精品国内亚洲2022精品成人| 大型av网站在线播放| 老司机靠b影院| 婷婷丁香在线五月| 日韩欧美三级三区| 一本大道久久a久久精品| 男人的好看免费观看在线视频 | 欧美日韩乱码在线| 在线观看www视频免费| 老司机午夜十八禁免费视频| 久久天堂一区二区三区四区| 男男h啪啪无遮挡| 色精品久久人妻99蜜桃| 欧美又色又爽又黄视频| 亚洲av电影不卡..在线观看| 老汉色av国产亚洲站长工具| 又黄又爽又免费观看的视频| 欧美日韩一级在线毛片| 日本 av在线| 人妻丰满熟妇av一区二区三区| 18禁美女被吸乳视频| 999久久久精品免费观看国产| 国产成人欧美在线观看| 国产精品二区激情视频| 欧美日韩瑟瑟在线播放| 操出白浆在线播放| 日韩欧美一区视频在线观看| 欧美一级a爱片免费观看看 | 久久草成人影院| 十八禁网站免费在线| 欧美精品啪啪一区二区三区| 别揉我奶头~嗯~啊~动态视频| 在线观看一区二区三区| 99国产精品99久久久久| 免费看十八禁软件| 久久欧美精品欧美久久欧美| 欧美大码av| 免费观看人在逋| 老汉色∧v一级毛片| 成年女人毛片免费观看观看9| 免费高清视频大片| 国产精品国产高清国产av| 麻豆成人午夜福利视频| 国产高清激情床上av| 亚洲av成人一区二区三| 午夜福利视频1000在线观看| 午夜福利18| 日韩大尺度精品在线看网址| 日韩av在线大香蕉| 亚洲av成人不卡在线观看播放网| 亚洲天堂国产精品一区在线| 国产亚洲精品久久久久5区| 精品一区二区三区视频在线观看免费| 欧美人与性动交α欧美精品济南到| www.自偷自拍.com| e午夜精品久久久久久久| 天天躁狠狠躁夜夜躁狠狠躁| 日韩一卡2卡3卡4卡2021年| 久久久国产欧美日韩av| 欧美一级毛片孕妇| 巨乳人妻的诱惑在线观看| 在线播放国产精品三级| 久久狼人影院| 成人国语在线视频| 久久狼人影院| 高清在线国产一区| 国产精品二区激情视频| 一区二区三区国产精品乱码| 久久中文字幕一级| 美女国产高潮福利片在线看| 黄色视频,在线免费观看| 黄色毛片三级朝国网站| av在线天堂中文字幕| 搞女人的毛片| 成在线人永久免费视频| 一级毛片高清免费大全| 国产亚洲精品久久久久5区| 国产精品免费视频内射| 两性夫妻黄色片| av天堂在线播放| 亚洲男人的天堂狠狠| 欧美黄色片欧美黄色片| 精品高清国产在线一区| 亚洲av熟女| 999久久久精品免费观看国产| 亚洲欧美激情综合另类| 午夜a级毛片| 欧美一级a爱片免费观看看 | 亚洲欧美一区二区三区黑人| 女同久久另类99精品国产91| 国产极品粉嫩免费观看在线| 久久国产乱子伦精品免费另类| 亚洲第一电影网av| 国产一区二区激情短视频| 老司机深夜福利视频在线观看| 亚洲成av人片免费观看| 人人妻人人澡欧美一区二区| 在线播放国产精品三级| 露出奶头的视频|