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

    A New Problem with Cross-Species Amplification of Microsatellites: Generation of Non-Homologous Products

    2010-12-25 01:03:04YUEGenHuaBalazsKovacsLaszloOrban
    Zoological Research 2010年2期
    關鍵詞:異源微衛(wèi)星同源

    YUE Gen-Hua, Balazs Kovacs, Laszlo Orban

    (1. Reproductive Genomics, Strategic Research Program, Temasek Life Sciences Laboratory, Singapore; 2. Regional University Center of Excellence in Environmental Industry Based on Natural Resources, G?d?ll?, Hungary; 3. Department of Biological Sciences, National University of Singapore, Singapore; 4. Molecular Population Genetics Group, Strategic Research Program, Temasek Life Sciences Laboratory, Singapore)

    A New Problem with Cross-Species Amplification of Microsatellites: Generation of Non-Homologous Products

    YUE Gen-Hua1,4,*, Balazs Kovacs2, Laszlo Orban1,3,*

    (1.Reproductive Genomics, Strategic Research Program, Temasek Life Sciences Laboratory, Singapore; 2.Regional University Center of Excellence in Environmental Industry Based on Natural Resources, G?d?ll?, Hungary; 3.Department of Biological Sciences, National University of Singapore, Singapore; 4.Molecular Population Genetics Group, Strategic Research Program, Temasek Life Sciences Laboratory, Singapore)

    Microsatellites have been widely used in studies on population genetics, ecology and evolutionary biology. However, microsatellites are not always available for the species to be studied and their isolation could be time-consuming. In order to save time and effort researchers often rely on cross-species amplification. We revealed a new problem of microsatellite cross-species amplification in addition to size homoplasy by analyzing the sequences of electromorphs from seven catfish species belonging to three different families (Clariidae, Heteropneustidae and Pimelodidae). A total of 50 different electromorphs were amplified from the seven catfish species by using primers for 4 microsatellite loci isolated from the speciesClarias batrachus. Two hundred and forty PCR-products representing all 50 electromorphs were sequenced and analyzed. Primers for two loci amplified specific products from orthologous loci in all species tested, whereas primers for the other two loci produced specific and polymorphic bands from some non-orthologous loci, even in closely related non-source species. Size homoplasy within the source species was not obvious, whereas extensive size homoplasy across species were detected at three loci, but not at the fourth one. These data suggest that amplification of products from non-orthologous loci and appearance of size homoplasy by cross-amplification are locus dependent, and do not reflect phylogenetic relationship. Amplification of non-orthologous loci and appearance of size homoplasy will lead to obvious complications in phylogenetic interference, population genetic and evolutionary studies. Therefore, we propose that sequence analysis of cross-amplification products should be conducted prior to application of cross-species amplification of microsatellites.

    Microsatellite; Polymorphism; Evolution; Non-orthologous loci

    Microsatellites are short tandem repeat DNA sequences with the unit length of 1 to 6 base pairs (Weber & May, 1989). Because they are highly polymorphic, co-dominant in nature, easy to score by PCR and rather abundant in most organisms studied, they have been widely used for the study of linkage mapping, comparative mapping, demographic structure and phylogenetic history in populations (Goldstein & Schlotterer, 1999; Zhang et al, 2001). However, microsatellites are not always available for the species to be studied and their isolation could be time-consuming (Lin et al, 2008; Wang et al, 2008). In order to save time and effort researchers often rely on cross-species amplification (Chang et al, 2008; Küpper et al, 2008; Kayser et al, 1996; Kijas et al, 1995; Lin et al, 2008). This procedure uses PCR primers complementary to the flanking regions of loci from a extensively studied (source) species to amplify microsatellites from closely (Harr et al, 1998) or sometimes quite distantly related species (Gonzalez-Martinez et al, 2004) for which no such markers are described. One problem related to cross-species amplification is size homoplasy (Anmarkrud et al, 2008; Estoup et al, 1995). PCR products of microsatellite loci with the same fragment length, but different sequence can arise from mutational events (deletion or insertion) in the flanking regions of the repeats or by interruptions in a perfect repeat producing alleles of the same size, which however are not identical by decent. Microsatellite size homoplasy has been reported in a number of papers (Hempel & Peakall, 2003; Makova et al, 2000; van Oppen et al, 2000) and was thought be a major problem of cross-amplification. It seems that size homoplasy increases with time divergence among populations and taxa (Estoup et al, 1995). However, a current study showed that homoplasy at microsatellite electromorphs did not represent a significant problem for many types of population genetics analyses performed by molecular ecologists, as the extensive variability at microsatellite loci often compensated for their homoplasious evolution (Estoup et al, 2002).

    In this paper, we describe a new problem of applying microsatellites for several different taxa. Cross-species amplification of microsatellites generated polymorphic products from non-orthologous loci, which were revealed by sequence analysis of 240 clones representing all 50 electromorphs from four loci in seven species (Clarias batrachus, C. fuscus, C. gariepinus, C. macrocephalus, Heterobranchus longfilis, Heteropneustes fossilisandPhractocephalus hemioliopterus).

    1 Materials and Methods

    1.1 Species and phylogenetic analyses

    Seven species of catfish were used in this study, namely:Clarias batrachus(abbreviation:Cba; the source species),C. fuscus(Cfu),C. gariepinus(Cga),C. macrocephalus(Cma),Heterobranchus longfilis(Hlo),Heteropneustes fossilis(Hfo), andPhractocephalus hemioliopterus(Phe). According to the current taxonomical system, five of the species studied were from the Clariidae family, one (Heteropneustes fossilis) from the Heteropneustidae family, which is closely related to Clariidae and the last (P. hemioliopterus) from the more distant Pimelodidae family. In order to determine the exact evolutionary relationship among the seven catfish species, phylogenetic analyses were conducted on the basis of the partial sequences ofcytbgenes from their mitochondrial genome. The sequences of six speciesC. batrachus[AF235932],C. fuscus[AF416885],C.gariepinus[AF126823],C. macrocephalus[AJ548464],Heterobranchus longfilis[AY995125], andHeteropneustes fossilis[AF126828] were downloaded from Genbank, whereas the one ofP. hemioliopteruswas amplified with PCR and sequenced as described (Agnese & Teugels, 2005). The sequence of the cytb gene of the Asian arowana (Scleropages formosus; DQ023143) was used as an outgroup. All seven sequences were aligned using Clustal_X (Thompson et al, 1997), and a NJ tree was reconstructed using the Kimura-2 parameter model of nucleotide using MEGA 3.0 (Kumar et al, 2001). The partial sequence of the cytb gene ofP. hemioliopteruswas deposited in GenBank under the accession number DQ200272.

    1.2 Sequencing of electromorphs generated by cross-species amplification

    All 50 electromorphs (Tabs. 1-4) generated in an earlier study (Yue et al, 2003) from four microsatellites (Cba01,Cba03,Cba06andCba20) from each of the seven species were used for cloning and sequencing. PCR products (25 μL) were cleaned using a glassmilk-based optimized procedure described earlier (Yue et al, 2007; Yue & Orban, 2001) prior to ligation of the fragments in to the pGEM-T-Easy vector (Promega) and subsequent transformation into XL-10 gold ultracompetent cells (Stratagene). Colonies were subjected to white/blue selection, and the insert of selected white clones was amplified by colony PCR as described (Yue et al, 2000). Un-incorporated PCR primers were removed by treating 5 μL PCR product for each clone with 0.5 unit shrimp alkalic phosphatase (SAP; USB) and 0.2 unit Exonuclease I (ExoI; USB) in 1× SAP buffer at 37℃ for 30 min, followed by a treatment at 80℃ for 15 min to inactivate the enzymes. One μL treated PCR product was directly used as template for sequencing from both directions using a BigDye kit (Applied Biosystems) and either M13 forward or M13 reverse primer in a PTC-100 PCR machine (MJ Research). Electrophoretic separation of the sequencing products was performed by using an ABI3730xl sequencer (Applied Biosystems). In order to exclude the possibility of cloning artifacts, for each electromorph from each species, multiple clones (at least3) were sequenced. Altogether the following number of clones were sequenced for the four microsatellite types:Cba01–107 clones,Cba03–20 clones,Cba06–50 clones andCba20–63 clones. Alignment of sequences was carried out by using Clustal X (Thompson et al, 1997).

    Tab. 1 Electromorphs amplified by the primer pair designed for Cba01 in seven catfish species

    Tab. 2 Electromorphs amplified by the primer pair designed for Cba03 in seven catfish species

    Tab. 3 Electromorphs amplified by the primer pair designed for Cba06 in seven catfish species

    Tab. 4 Electromorphs amplified by the primer pair designed for Cba20 in seven catfish species

    2 Results

    2.1 Phylogenetic relationship of the seven catfish species

    Fig. 1 Phylogenetic relationship among the seven catfish species

    Based on the partial sequences of thecytbgene of the seven species, a NJ tree was constructed (Fig. 1). The three speciesClarias batrachus,C. fuscusandC. macrocephaluswere closely related and clustered into a group. This group was linked to the group ofC. gariepinusandHeterobranchus longifilis. The remaining two species:HeteropneustesfossilisandPhractocephalus hemioliopteruswere distantly related to other five species.

    2.2 Sequence analysis of electromorphs amplified by the Cba01 primer pair

    The primer pair designed to theCba01locus amplified polymorphic products in all seven catfish species tested. Altogether 23 clear bands (eletromorphs) were detected in the seven species (size range: 199-349 bp), their sequencing analyses uncovered the total of 34 different alleles (Tab. 1). InC. fuscus,C. macrocephalusandP. hemioliopterusboth the repeat and the flanking regions exhibited high similarity to source sequences fromC. batrachus(Fig. 2A). On the other hand, the corresponding sequences fromC. gariepinus, andHeteropneustes fossilisspecies were completely different from the source sequences (Fig. 2B), but quite similar among these three species. The length ofHeterobranchus longifilisalleles was similar to those of the source species, but the flanking region and repeats were entirely different (Fig. 2C).

    The 5' and 3' flanking sequences for each allele were nearly identical in different individuals ofC. batrachus,C. fuscus,C. macrocephalusandP. hemioliopterus, respectively. On the other hand, several differences were found between sequences from different species both at the 5' and 3' flanking regions (seven and eight positions, respectively). Most of them seem to have been caused by substitution, whereas the rest by insertion or deletion of a single base pair. A notable feature is, that the repeat structures of this locus were slightly different in these four species: (GC)2(AC)nin the source species, (GC)3GT(GC)5-6(AC)5(GC)0-1(AC)ninC. fuscusandP. hemioliopterus, whereas (GC)2-5(AC)0-1(GC)0-4(AC)0-2GC(AC)ninC. macrocephalus(Fig. 2A). Therefore, the polymorphism at this locus was caused by change in the number of either AC or GC repeat units in different species, resulting in fragments of the same length, but with quite different sequences. Within species, size homoplasy could only be detected inC. macrocephalus, but not in the source species,C. fuscusorP. hemioliopterus.

    Fig. 2 Sequence alignment of some electromorphs (amplified by the primer pair designed for Cba01) from seven different catfish species

    InC. gariepinusandH. fossilis,the sequences of the 7 electromorphs (Tab. 1) were different from those in source species. The flanking sequences were quite similar among different alleles, although the polyA and polyT repeats (located at the 5' and 3' flanking regions, respectively) showed polymorphism both within and among species (Fig. 2B). Moreover, a deletion of 16 bp was detected in the 5' flanking region ofHeteropneustes fossilis(data not shown). InC. gariepinus(but not inH. fossilis) a CAG unit was deleted from the 3' flanking region. A few point mutations, short deletions or insertions have also been detected in the 5' and 3' flanking regions among electromorphs from different species (data not shown). Polymorphism in the repeat at this locus was caused either by a change in the length of polyA stretch in the 5' flanking region, or by the unit number of (GA)n, (GAA)n, (GGA)ncompound repeats or by a deletion of three base pairs CAG and a change in the length of the polyT in the 3' flanking region (Fig. 2B).

    InH. fossilis,the locus appeared to be duplicated, because more than two bands were detected in the PCR product of each individual tested, whereas no such phenomenon was observed in the other two species. The 199 bp allele from all six individuals ofH. fossilistested (Genebank No. AY196549) lacked a 150 bp fragment including the 5' flanking region and even the whole repeat region as compared with the largest allele (Hfo349) (Fig. 2B).

    InHeterobranchus longifilis, the sequences of electromorphs were entirely differently from the alleles of the source species, although the length of the electromorphs was similar to those of the source species (Fig. 2C). The length polymorphism of the electromorphs was caused by the change of number of CT repeats.

    2.3 Sequence analysis of electromorphs amplified by the Cba03 primer pair

    At theCba03locus, a total of three electromorphs (range: 129 - 135 bp) were detected across the seven species (Tab. 2). Sequencing of each electromorph (20 clones) revealed that the sequence of this locus was highly conserved across the catfish species studied (Fig. 3). The polymorphism was caused exclusively by the change in the unit number of the (GGA)nrepeat. At three positions of 3' flanking region, single base pair substitution was also seen in two species (C. macrocephalusandP. hemiolopterus). No size homoplasy was identified among individuals of any species.

    2.4 Sequence analysis of electromorphs amplified by the Cba06 primer pair

    At theCba06locus, a total of 11 electromorphs (range 168 -258 bp) were identified across the seven species (Tab. 3). Their sequence analysis demonstrated that they could be divided into two groups and two individual sequences (Fig.4A - D). Fragments amplified fromC. fuscus(1 allele) andC. macrocephalus(4 alleles) showed an overall high similarity to the source sequence (Fig. 4A). In these two species, an insertion of a 34 bp fragment was detected at the 5' flanking region between the primer and repeats in every allele in comparison to the source sequence. Additional single base pair substitutions, located in the flanking regions were also found. The length polymorphism was caused by the change in the unit number of the (AAC)nrepeat within each species, but among species the length polymorphism could also be caused by change in the extent of polyA in the 3' flanking region or the insertion of a 34 bp fragment into the 5' flanking region. Although no size homoplasy was identified within these two species, its presence was quite obvious among species. For example, the 245 bp electromorph inC. fuscusand that inC. macrocephalusshowed different unit number of CAA-repeats and appearance of a CTA sequence due to an A→T mutation in the latter.

    Fig. 3 Sequence alignment of some electromorphs (amplified by the primer pair designed for Cba03) from seven different catfish species shows no size homoplasy within or among species

    The second group (Fig. 4B) included sequences fromC. gariepinus(2 alleles), andH. longifilis(1). The DNA sequence of the fragments from the two species showed high similarity to each other, but differed from the source sequence both in their flanking regions and repeat motif [(AAC)nvs. (CA)n]. An insertion of five base pairs (CGAAC) was seen in the 5' flanking region of the speciesH. longifilis, as compared the sequences from theC. gariepinus(Fig. 4B). Apart from this insertion, the length polymorphism was caused by the different number of the (AC)nrepeat units in all fragments. Between the two species, single base pair substitution was observed at several positions of the flanking regions. The 168 bp fragment appeared in both species. However comparison of sequences between the two species revealed two different alleles.

    The remaining two sequences (Fig. 4C-D; GenBank Nos. AY196578 and AY196579) originated fromP. hemioliopterusandHeteropneustes fossilis, respectively. They did not show any similarity to the first two groups except the primer binding sites and did not contain repeats.

    2.5 Sequence analysis of electromorphs amplified by the Cba20 primer pair

    Fig. 4 Sequence alignment of some electromorphs (amplified by the primer pair designed for Cba06) from seven different catfish species

    A total of 13 electromorphs (range: 93 - 143 bp) were detected across six species (Tab. 4), but not inP.hemioliopterus. Sequence analysis revealed 10 additional alleles (Fig. 5), without any evidence of homoplasy within the source species. In the 5' flanking region, single base pair substitutions were detected at three positions among species. As compared with the source (4 alleles), sequences fromC. macrocephalus(4) andHeterobranchus longifilis(3) showed an insertion of three basepairs (GTC) in the 3' flanking region. Single basepair substitutions were also detected at two positions of the 3' flanking regions. The repeat region was highly variable within and among species. In the source species, repeat structure for the 95 bp allele was (TC)6GC(TC)2, although longer and shorter alleles showed change in repeat number of longer repeat, the GC(TC)2motif remained constant among all alleles. InC. fuscus(3 alleles), where the (TC)nrepeat was interrupted by a TA unit, the (TC)3upstream from the TA remained unchanged, whereas the downstream (TC)nrepeat showed polymorphism among individuals. InC. gariepinus(2 alleles) the TC repeats were interrupted by GC and TG units at several positions and the polymorphism was caused by the change of the long, upstream TC repeat, whereas the shorter ones remained constant. InC. macrocephalusandHeteropneustesfossilis(3 alleles), the (TC)nrepeat was interrupted by CC, TT and GT motifs, whereas inHeterobranchus longifilisby GC, AG and TG units. The reason for the polymorphism was similar to that described forC. gariepinus.

    Fig. 5 Sequence alignment of some electromorphs (amplified by the primer pair Cba20) from seven different catfish species

    3 Discussion

    Microsatellites are very useful tools for genetic and evolutionary studies. However, their genotyping is based on prior sequence information from the genome to be analyzed. Despite of recent improvements on the procedure (for review see: Zane et al, 2002) the isolation of microsatellites is still cumbersome. One of the possible solutions for this problem is cross-species amplification, which involves the use of primer pairs designed for the flanking region of conserved microsatellites (of a so-called source species) for genotyping in related species amplification (Housley et al, 2006; Kayser et al, 1996; Kijas et al, 1995). Data for several such experiments have been reported in teleosts during the last decade (e.g. Koskinen & Primmer, 1999; Yue et al, 2004; Yue et al, 2003). However all PCR products generated in the non-source species have only been analyzed at the sequence level in a few cases (Kayang et al, 2002; Viard et al, 1998). We have tested the applicability of four conserved microsatellite markers isolated earlier fromC. batrachus(Yue et al, 2003) on six additional catfish species. We found that PCR primer pairs designed for the flanking regions of the fourC. batrachusmicrosatellite loci amplified products in most of the related species. However, sequencing analyses of 240 clones representing 50 electromorphs from seven catfish species revealed a new problem of cross-species amplification of microsatellites: the generation of non-orthologous loci, beside the appearance of size homoplasy. Primer pairs designed for twoC. batrachusloci (Cba03andCba20) amplified highly similar (orthologous) sequence products in all non-source species. On the other hand, those designed for other two loci (Cba01andCba06) yielded polymorphic products with entirely different sequence from some of the distantly related species (e.g.P. hemioliopterusandHeteropneustes fossilis), and even in closely related species (e.g.C. gariepinus) indicating that these bands originated from non-orthologous loci. The amplification of specific products from non-orthologous source was locus-dependent, and did not reflect the phylogenetic relationship. Thus, in the absence of sequence information it would be very difficult to predict whether certain primer pairs will amplify products from orthologous loci in a given non-source species or not. Similar phenomenon was observed earlier in soybean (Peakall et al, 1998) and rice (Chen et al, 2002), but those findings have not been analyzed in detail. Taken together, our data suggest that generation of polymorphic products from non-orthologous loci by cross-species amplification is not a unique feature of certain taxonomic groups in fish, instead it might occur throughout the animal and plant kingdom. Although the mechanisms underlying this phenomenon are not fully understood, they are thought to be related to genome and gene duplication, as well as speciation. Such events are expected be more frequent in fish, since the ancestor of today’s teleosts seems to have experienced an additional round of genome duplication (Meyer & Schartl, 1999; Postlethwait et al, 2000) and chromosome duplications (Chang et al, 2005) after their ancestor has split from that of the other vertebrates. Duplication of microsatellite loci followed by gene conversion can lead to amplification of non-orthologous loci as proposed (Angers et al, 2002).

    Sequencing of all alleles of four microsatellite loci in the source and six non-source species showed that length difference of microsatellites was not restricted to their repeat regions. A longer insertion and several shorter insertions were detected in the flanking region of the loci orthologous toCba06in non-source species. At theCba01,Cba06andCba20loci, a number of alleles from different non-source species showed the same length, but with different sequences. At the same time, atCba03locus electromorphs of the same length represented the same sequences, suggesting that size homoplasy for microsatellite markers produced by cross-species amplification is locus-dependent, it does not reflect the phylogenetic relationship. We also found the tendency of increase in the number of interrupted repeats of orthologous loci in non-source species, as observed by others in different taxonomic groups (e.g. Culver et al, 2001; Di Gaspero et al, 2000; Estoup et al, 1995; Garza et al, 1995; van Oppen et al, 2000). This tendency also seems to be locus-dependent in catfish, since two loci (Cba01andCba20) showed clear interruptions in non-source species, whereas the other two (Cba03andCba06) exhibited no or few interruptions in them.

    Applications of microsatellites to population genetics, ecological and evolutionary studies rely heavily on the models used for explaining the mutational process of these markers. However, all models relay on the assumption that differences between alleles at orthologous loci are due entirely to changes in the number of repeats. In this study, we demonstrated that appearance of size homoplasy and amplification of non-orthologous products by cross-species amplification were locus-dependent, and did not reflect phylogenetic relationships. Therefore, application of cross amplification of microsatellites to population genetics and phylogenetic analyses in distantly or even in closely related species, might make the interpretation of length difference of electromorphs difficult and cause wrong estimation of evolutionary relationship.

    In conclusion, we revealed a new problem of microsatellite cross-species amplification, namely amplification of non-orthologous loci, besides the well-known problem (size homoplasy). The new problem and appearance of size homoplasy will lead to obvious complications for phylogenetic interferences, population genetics, mapping and evolutionary studies. The sequence analysis of products generated by“cross-species primers” should always be performed, as it could reveal previously unrecognized problems and might allow for extracting more information from these loci, thereby increasing their usefulness.

    Agnese, JF, Teugels, GG. 2005. Insight into the phylogeny of African Clariidae (Teleostei, Siluriformes): Implications for their body shape evolution, biogeography, and taxonomy [J].Mol Phylogenet Evol, 36 (3): 546-553.

    Angers B, Gharbi K, Estoup A. 2002. Evidence of gene conversion events between paralogous sequences produced by tetraploidization in Salmoninae fish [J].J Mol Evol, 54 (4): 501-510.

    Anmarkrud JA, Kleven O, Bachmann L, Lifjeld JT. 2008. Microsatellite evolution: Mutations, sequence variation, and homoplasy in the hypervariable avian microsatellite locus HrU10 [J].BMC Evol Biol, 8: 138.

    Chang CH, Hsieh LC, Chen TY, Chen HD, Luo L, Lee HC. 2005. Shannon information in complete genomes [J].J Bioinform Comput Biol, 3 (3): 587-608.

    Chang YM, Kuang YY, Liang LQ, Lu CY, He JG, Su XW. 2008. Searching for protein-coding genes using microsatellites in common carp by comparing to zebrafish EST database [J].Zool Res, 29 (4): 373-378.

    Chen X, Cho YG, McCouch SR. 2002. Sequence divergence of rice microsatellites inOryzaand other plant species [J].Mol GenetGenomics, 268 (3): 331-343.

    Culver M, Menotti-Raymond MA, O'Brien SJ. 2001. Patterns of size homoplasy at 10 microsatellite loci in pumas (Puma concolor) [J].Mol Biol Evol, 18 (6): 1151-1156.

    Di Gaspero G, Peterlunger E, Testolin R, Edwards KJ, Cipriani G. 2000. Conservation of microsatellite loci within the genusVitis[J].Theor Appl Genet, 101 (1-2): 301-308.

    Estoup A, Jarne P, Cornuet JM. 2002. Homoplasy and mutation model at microsatellite loci and their consequences for population genetics analysis [J].Mol Ecol, 11 (9): 1591-1604.

    Estoup A, Tailliez C, Cornuet JM, Solignac M. 1995. Size homoplasy and mutational processes of interrupted microsatellites in two bee species,Apis melliferaandBombus terrestris(Apidae) [J].Mol Biol Evol, 12 (6): 1074-1084.

    Garza JC, Slatkin M, Freimer NB. 1995. Microsatellite allele frequencies in humans and chimpanzees, with implications for constraints on allele size [J].Mol Biol Evol, 12 (4): 594-603.

    Goldstein DB, Schlotterer C. 1999. Microsatellites: Evolution and Applications [M]. Oxford: Oxford University Press.

    Gonzalez-Martinez SC, Robledo-Arnuncio JJ, Collada C, Diaz A, Williams CG, Alia R, Cervera MT. 2004. Cross-amplification and sequence variation of microsatellite loci in Eurasian hard pines [J].Theor Appl Genet, 109 (1): 103-111.

    Harr B, Zangerl B, Brem G, Schlotterer C. 1998. Conservation of locus-specific microsatellite variability across species: A comparison of twoDrosophilasibling species,D. melanogasterandD. simulans[J].Mol Biol Evol, 15 (2): 176-184.

    Hempel K, Peakall R. 2003. Cross-species amplification from crop soybeanGlycine maxprovides informative microsatellite markers for the study of inbreeding wild relatives [J].Genome, 46 (3): 382-393.

    Housley DJ, Zalewski ZA, Beckett SE, Venta PJ. 2006. Design factors that influence PCR amplification success of cross-species primers among 1147 mammalian primer pairs [J].BMC Genomics, 7: 253. Küpper C, Burke T, Székely T, Dawson DA. 2008. Enhanced cross-species utility of conserved microsatellite markers in shorebirds [J].BMC Genomics, 9: 502.

    Kayang BB, Inoue-Murayama M, Hoshi T, Matsuo K, Takahashi H, Minezawa M, Mizutani M, Ito S. 2002. Microsatellite loci in Japanese quail and cross-species amplification in chicken and guinea fowl [J].Genet Sel Evol, 34 (2): 233-253.

    Kayser M, Ritter H, Bercovitch F, Mrug M, Roewer L, Nurnberg P. 1996. Identification of highly polymorphic microsatellites in the rhesus macaqueMacaca mulattaby cross-species amplification [J].Mol Ecol, 5 (1): 157-159.

    Kijas JM, Fowler JC, Thomas MR. 1995. An evaluation of sequence tagged microsatellite site markers for genetic analysis within Citrus and related species [J].Genome, 38 (2): 349-355.

    Koskinen MT, Primmer CR. 1999. Cross-species amplification of salmonid microsatellites which reveal polymorphism in European and Arctic grayling, Salmonidae:Thymallusspp [J].Hereditas, 131 (2): 171-176.

    Kumar S, Tamura K, Jakobsen IB, Nei M. 2001. MEGA2: molecular evolutionary genetics analysis software [J].Bioinformatics, 17: 1244-1245.

    Lin G, Chang A, Yap W, Yue GH. 2008. Characterization and cross-species amplification of microsatellites from the endangered Hawksbill turtle (Eretmochelys imbricate) [J].Conserv Genet, 9: 1071-1073.

    Makova KD, Nekrutenko A, Baker RJ. 2000. Evolution of microsatellite alleles in four species of mice (genusApodemus) [J].J Mol Evol, 51 (2): 166-172.

    Meyer A, Schartl M. 1999. Gene and genome duplications in vertebrates: the one-to-four (- to-eight in fish) rule and the evolution of novel gene functions [J].Curr Opin Cell Biol, 11 (6): 699-704.

    Peakall R, Gilmore S, Keys W, Morgante M, Rafalski A. 1998. Cross-species amplification of soybean (Glycine max) simple sequence repeats (SSRs) within the genus and other legume genera: Implications for the transferability of SSRs in plants [J].Mol Biol Evol, 15 (10): 1275-1287.

    Postlethwait JH, Woods IG, Ngo-Hazelett P, Yan YL, Kelly PD, Chu F, Huang H, Hill-Force A, Talbot WS. 2000. Zebrafish comparative genomics and the origins of vertebrate chromosomes [J].Genome Res, 10 (12): 1890-1902.

    Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. 1997. The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools [J].Nucleic Acids Res, 25 (24): 4876-4882.

    van Oppen MJH, Rico C, Turner GF, Hewitt GM. 2000. Extensive homoplasy, nonstepwise mutations, and shared ancestral polymorphism at a complex microsatellite locus in LakeMalawi cichlids[J].Mol Biol Evol, 17 (4): 489-498.

    Viard F, Franck P, Dubois MP, Estoup A, Jarne P. 1998. Variation of microsatellite size homoplasy across electromorphs, loci, and populations in three invertebrate species [J].J Mol Evol, 47 (1): 42-51.

    Wang HZ, Yin QQ, Feng ZG, Li, DY, Sun XW, Li C. 2008. Construction of fractional genomic libraries and screening microsatellites DNA ofEsox reiehertiDybowski [J].Zool Res, 29 (3): 245-252.

    Weber JL, May PE. 1989. Abundant class of human DNA polymorphisms which can be typed using the polymerase chain-reaction [J].Am J Hum Genet, 44 (3): 388-396.

    Yue GH, Chen F, Orban L. 2000. Rapid isolation and characterization of microsatellites from the genome of Asian arowana (Scleropages formosus, Osteoglossidae, Pisces) [J].Mol Ecol, 9 (7): 1007-1009. Yue GH, David L, Orban L. 2007. Mutation rate and pattern of microsatellites in common carp (Cyprinus carpioL.) [J].Genetica, 129 (3): 329-31.

    Yue GH, Ho MY, Orban L, Komen J. 2004. Microsatellites within genes and ESTs of common carp and their applicability in silver crucian carp [J].Aquaculture, 234 (1-4): 85-98.

    Yue GH, Kovacs B, Orban L. 2003. Microsatellites fromClarias batrachusand their polymorphism in seven additional catfish species [J].Mol Ecol Notes, 3 (3): 465-468.

    Yue GH, Orban L. 2001. Rapid isolation of DNA from fresh and preserved fish scales for polymerase chain reaction [J].Mar Biotechnol, 3 (3): 199-204.

    Zane L, Bargelloni L, Patarnello T. 2002. Strategies for microsatellite isolation: a review [J].Mol Ecol, 11 (1): 1-16.

    Zhang YW, Zhang YP, Aryder O. 2001. Microsatellites and its application [J].Zool Res, 22 (4): 315-320.

    微衛(wèi)星跨物種交叉PCR擴增的一個新問題:擴增非同源產物

    岳根華1,4,*, Balazs Kovacs2, Laszlo Orban1,3,*

    (1.Reproductive Genomics, Strategic Research Program, Temasek Life Sciences Laboratory, Singapore; 2.Regional University Center of Excellence in Environmental Industry Based on Natural Resources, G?d?ll?, Hungary3.Department of Biological Sciences, National University of Singapore, Singapore; 4.Molecular Population Genetics Group, Strategic Research Program, Temasek Life Sciences Laboratory, Singapore)

    微衛(wèi)星已被廣泛應用于群體遺傳學、生態(tài)學和進化生物學研究。然而,一些物種微衛(wèi)星尚未克隆。為了節(jié)省時間和經費,研究人員往往使用一個物種已發(fā)表的微衛(wèi)星引物擴增其近緣物種的微衛(wèi)星。該研究對屬于 3個不同科(Clariidae、Heteropneustidae 和Pimelodidae)的7個鯰魚物種的微衛(wèi)星跨物種PCR擴增產物進行了序列分析,研究發(fā)現(xiàn)擴增非同源(non-orthologous)產物是微衛(wèi)星跨物種PCR擴增的一個新問題。該研究共采用4對胡子鯰微衛(wèi)星座位引物對7個鯰魚物種進行了跨物種PCR擴增。對獲得的204個PCR產物的序列分析結果表明,兩對微衛(wèi)星座位引物擴增了所有7個物種的同源特異產物。而其他兩個座位的引物擴增了特異但非同源的多態(tài)產物,對近緣物種的擴增也獲得類似結果。另外,除胡子鯰等位基因大小異源同型(size homoplasy)的特征不明顯外,其他物種在3個微衛(wèi)星座位都具有這一非常明顯的特征。這些數(shù)據(jù)表明,微衛(wèi)星跨物種間交叉擴增能產生非同源產物;等位基因大小異源同型與微衛(wèi)星座位本身有關,而與物種間的親緣關系無明顯的相關性。微衛(wèi)星跨物種擴增產生的非同源產物和等位基因大小異源同型將使系統(tǒng)發(fā)育、群體遺傳學和進化研究明顯復雜化。因此,在應用微衛(wèi)星跨物種交叉擴增數(shù)據(jù)以前,最好對跨物種交叉擴增產物進行測序驗證。

    2009-09-15;接受日期:2009-12-31

    book=132,ebook=192

    微衛(wèi)星;多態(tài)性;進化;非同源座位

    Q754;Q343.1;Q984.403

    A

    0254-5853-(2010)02-0131-10

    10.3724/SP.J.1141.2010.02131

    date: 2009-09-15; Accepted date: 2009-12-31

    This study was supported financially by the internal research funding from Temasek Life Sciences Laboratory

    *Corresponding authors (通訊作者),F(xiàn)ax: 65-6872-7007,E-mail: genhua@tll.org.sg; laszlo@tll.org.sg

    The authors would like to thank Drs. Graham Mair, Arlo Fast, Bela Urbanyi and Lian Chuan Lim, as well as Ferenc Radics, Judit Raczkevi and Gyula Pasareti for supplying fin clips from various catfish species, and the Strategic Research Program of TLL for financial support. B.K. is grateful for the support of the Temasek Life Sciences Laboratory, the Bolyai Research Fellowship of the Hungarian Academy of Sciences, and the Hungarian Scientific Research Fund (OTKA PD79177)..

    猜你喜歡
    異源微衛(wèi)星同源
    藥食同源
    ——紫 蘇
    兩岸年味連根同源
    華人時刊(2023年1期)2023-03-14 06:43:36
    梅花鹿基因組微衛(wèi)星分布特征研究
    特產研究(2022年6期)2023-01-17 05:05:06
    近岸水體異源遙感反射率產品的融合方法研究
    海洋通報(2022年4期)2022-10-10 07:40:26
    基于EM-PCNN的果園蘋果異源圖像配準方法
    以同源詞看《詩經》的訓釋三則
    臉譜與假面 異源而殊流
    中華戲曲(2020年2期)2020-02-12 05:17:58
    林麝全基因組微衛(wèi)星分布規(guī)律研究
    四川動物(2017年4期)2017-07-31 23:54:19
    虔誠書畫乃同源
    解脂耶氏酵母異源合成番茄紅素的初步研究
    cao死你这个sao货| 日韩欧美国产一区二区入口| 久久久久久久久免费视频了| 丝袜人妻中文字幕| 91av网站免费观看| 一区二区日韩欧美中文字幕| 一区二区三区激情视频| 欧美少妇被猛烈插入视频| tocl精华| 国产亚洲精品久久久久5区| 亚洲专区中文字幕在线| 一级毛片女人18水好多| 午夜福利在线免费观看网站| 久久亚洲精品不卡| 成人av一区二区三区在线看 | 久久精品国产a三级三级三级| 国产精品久久久久成人av| svipshipincom国产片| 国产欧美日韩一区二区精品| 色94色欧美一区二区| 男人添女人高潮全过程视频| 亚洲精品第二区| 国产高清国产精品国产三级| 国产麻豆69| 男女边摸边吃奶| 黄片小视频在线播放| av在线老鸭窝| 夜夜骑夜夜射夜夜干| 多毛熟女@视频| 天天操日日干夜夜撸| av线在线观看网站| 91九色精品人成在线观看| 精品国产一区二区久久| 啦啦啦 在线观看视频| 性色av一级| 国产人伦9x9x在线观看| 丝袜人妻中文字幕| 黄片小视频在线播放| 九色亚洲精品在线播放| 可以免费在线观看a视频的电影网站| 国产成人精品在线电影| 黄色 视频免费看| 91麻豆av在线| 性色av乱码一区二区三区2| 在线观看人妻少妇| 脱女人内裤的视频| 一进一出抽搐动态| 国产成人系列免费观看| 久久国产亚洲av麻豆专区| 午夜精品久久久久久毛片777| 国产精品 国内视频| 岛国毛片在线播放| 丝瓜视频免费看黄片| 高清欧美精品videossex| 国产高清videossex| 菩萨蛮人人尽说江南好唐韦庄| 成在线人永久免费视频| 精品一品国产午夜福利视频| 久久久久久久大尺度免费视频| 黑人操中国人逼视频| 久久久久视频综合| 视频在线观看一区二区三区| 久久久欧美国产精品| 男人舔女人的私密视频| 精品福利观看| 中文字幕av电影在线播放| 黄色片一级片一级黄色片| 亚洲欧美精品自产自拍| 久久亚洲国产成人精品v| 99国产极品粉嫩在线观看| 在线 av 中文字幕| 啦啦啦视频在线资源免费观看| 人妻久久中文字幕网| 精品第一国产精品| 日本一区二区免费在线视频| 亚洲av电影在线进入| 国产欧美日韩一区二区三 | 国产不卡av网站在线观看| 欧美少妇被猛烈插入视频| a级毛片在线看网站| 欧美日本中文国产一区发布| 日韩制服丝袜自拍偷拍| 夜夜骑夜夜射夜夜干| 搡老乐熟女国产| 久久九九热精品免费| 国产精品九九99| avwww免费| 久久人人爽人人片av| 亚洲情色 制服丝袜| 一区二区av电影网| 国产免费福利视频在线观看| 国产主播在线观看一区二区| 热re99久久精品国产66热6| 50天的宝宝边吃奶边哭怎么回事| 亚洲一卡2卡3卡4卡5卡精品中文| 性色av乱码一区二区三区2| 国产精品久久久人人做人人爽| 天天躁夜夜躁狠狠躁躁| 亚洲熟女毛片儿| 在线观看舔阴道视频| 黑人欧美特级aaaaaa片| 少妇人妻久久综合中文| 国产欧美日韩一区二区精品| 1024视频免费在线观看| 国内毛片毛片毛片毛片毛片| 亚洲av日韩精品久久久久久密| 人妻 亚洲 视频| 美女午夜性视频免费| 12—13女人毛片做爰片一| 亚洲五月色婷婷综合| 超碰成人久久| 亚洲av电影在线观看一区二区三区| 成人三级做爰电影| av片东京热男人的天堂| 国产欧美日韩综合在线一区二区| 超色免费av| 国产精品久久久久久人妻精品电影 | 午夜老司机福利片| 啪啪无遮挡十八禁网站| 黑丝袜美女国产一区| 国产一区二区激情短视频 | 1024视频免费在线观看| 久久人妻福利社区极品人妻图片| 亚洲视频免费观看视频| 亚洲全国av大片| 91字幕亚洲| 中文字幕精品免费在线观看视频| 亚洲精品国产一区二区精华液| 丝袜人妻中文字幕| 亚洲视频免费观看视频| 日本欧美视频一区| 国产成人欧美| 亚洲国产精品成人久久小说| 日日摸夜夜添夜夜添小说| 嫁个100分男人电影在线观看| cao死你这个sao货| 中文字幕人妻丝袜一区二区| 精品一品国产午夜福利视频| 亚洲精品成人av观看孕妇| 美女福利国产在线| 伦理电影免费视频| 欧美黄色淫秽网站| 国产精品偷伦视频观看了| 亚洲黑人精品在线| 丝袜美腿诱惑在线| 在线天堂中文资源库| 久久亚洲精品不卡| 亚洲少妇的诱惑av| 久久久久网色| 午夜日韩欧美国产| av在线播放精品| 高清在线国产一区| 国产精品免费大片| 热99re8久久精品国产| 一区二区av电影网| 精品国产乱码久久久久久小说| 欧美成狂野欧美在线观看| 久久久久国内视频| a级毛片在线看网站| 一区福利在线观看| 欧美97在线视频| av线在线观看网站| 人人妻,人人澡人人爽秒播| 久久人人97超碰香蕉20202| 三上悠亚av全集在线观看| 国产一级毛片在线| 天天躁夜夜躁狠狠躁躁| 韩国高清视频一区二区三区| 另类精品久久| 侵犯人妻中文字幕一二三四区| 别揉我奶头~嗯~啊~动态视频 | 后天国语完整版免费观看| 亚洲欧美一区二区三区黑人| 国产一区有黄有色的免费视频| 纯流量卡能插随身wifi吗| 国产无遮挡羞羞视频在线观看| 国产视频一区二区在线看| a 毛片基地| 欧美日韩精品网址| 9热在线视频观看99| 亚洲精品成人av观看孕妇| 侵犯人妻中文字幕一二三四区| 久久久精品区二区三区| 一边摸一边做爽爽视频免费| 老司机靠b影院| 中文字幕制服av| 日韩欧美国产一区二区入口| 国产一区有黄有色的免费视频| 国产精品麻豆人妻色哟哟久久| 极品少妇高潮喷水抽搐| 日韩制服丝袜自拍偷拍| 亚洲熟女毛片儿| 又紧又爽又黄一区二区| 18禁黄网站禁片午夜丰满| 国产淫语在线视频| 国产日韩欧美视频二区| 精品一区二区三区av网在线观看 | 欧美另类一区| 欧美日韩成人在线一区二区| 免费在线观看完整版高清| 伊人久久大香线蕉亚洲五| 一本综合久久免费| 成人三级做爰电影| 日本av免费视频播放| 日本精品一区二区三区蜜桃| 黄网站色视频无遮挡免费观看| 国产高清国产精品国产三级| 国产色视频综合| 女人久久www免费人成看片| 丰满少妇做爰视频| 两个人看的免费小视频| 国产熟女午夜一区二区三区| 国产精品九九99| 欧美av亚洲av综合av国产av| 黄片播放在线免费| 别揉我奶头~嗯~啊~动态视频 | 18禁观看日本| 国产日韩欧美视频二区| 欧美性长视频在线观看| 中文字幕最新亚洲高清| 嫁个100分男人电影在线观看| 国产精品av久久久久免费| 一本一本久久a久久精品综合妖精| 国产精品影院久久| 亚洲 国产 在线| 亚洲少妇的诱惑av| 婷婷色av中文字幕| 国产亚洲精品一区二区www | 夜夜夜夜夜久久久久| cao死你这个sao货| 五月开心婷婷网| kizo精华| 丝瓜视频免费看黄片| 国产精品1区2区在线观看. | 啦啦啦啦在线视频资源| 美女福利国产在线| www.熟女人妻精品国产| 美女脱内裤让男人舔精品视频| 久久精品国产亚洲av香蕉五月 | 国产成人精品久久二区二区91| 日韩,欧美,国产一区二区三区| 考比视频在线观看| 一区二区日韩欧美中文字幕| 精品亚洲成a人片在线观看| 人人妻,人人澡人人爽秒播| 三上悠亚av全集在线观看| 在线观看一区二区三区激情| 久久久久久久大尺度免费视频| 两性午夜刺激爽爽歪歪视频在线观看 | 99国产综合亚洲精品| 老汉色av国产亚洲站长工具| 久久中文看片网| 久久精品aⅴ一区二区三区四区| 在线永久观看黄色视频| 国产av一区二区精品久久| 天堂俺去俺来也www色官网| 久久精品亚洲熟妇少妇任你| 91麻豆av在线| 一区二区av电影网| 欧美成狂野欧美在线观看| 国产男女超爽视频在线观看| 国产精品.久久久| 老汉色∧v一级毛片| 老司机靠b影院| 日韩一区二区三区影片| 国产精品久久久久久精品电影小说| 日日摸夜夜添夜夜添小说| 黄色怎么调成土黄色| 午夜免费鲁丝| 欧美日韩亚洲国产一区二区在线观看 | 国产一级毛片在线| xxxhd国产人妻xxx| 丰满迷人的少妇在线观看| 免费av中文字幕在线| 欧美日韩亚洲国产一区二区在线观看 | 国产日韩一区二区三区精品不卡| 日本撒尿小便嘘嘘汇集6| 久久久精品国产亚洲av高清涩受| 三级毛片av免费| 老司机午夜福利在线观看视频 | 国产av又大| 人成视频在线观看免费观看| 欧美97在线视频| 久久人人97超碰香蕉20202| 欧美亚洲日本最大视频资源| 又紧又爽又黄一区二区| 亚洲精品国产av成人精品| 搡老熟女国产l中国老女人| 日本av手机在线免费观看| 久久午夜综合久久蜜桃| 久久国产精品人妻蜜桃| 国产精品亚洲av一区麻豆| 亚洲色图 男人天堂 中文字幕| 中文字幕另类日韩欧美亚洲嫩草| 最近中文字幕2019免费版| 深夜精品福利| 日韩视频在线欧美| 久久久久久免费高清国产稀缺| 丝袜脚勾引网站| 一本久久精品| 国产精品免费大片| 日韩中文字幕视频在线看片| 三级毛片av免费| 一区二区三区精品91| cao死你这个sao货| 久久人妻福利社区极品人妻图片| 国产在线免费精品| 国产欧美日韩精品亚洲av| 国产av一区二区精品久久| 俄罗斯特黄特色一大片| 99国产综合亚洲精品| 久久精品国产亚洲av香蕉五月 | 91麻豆精品激情在线观看国产 | 黄片大片在线免费观看| 丝袜美腿诱惑在线| 欧美日韩亚洲综合一区二区三区_| 91成年电影在线观看| 成在线人永久免费视频| 99精品久久久久人妻精品| 成人黄色视频免费在线看| 香蕉国产在线看| 免费人妻精品一区二区三区视频| 国产精品久久久久久人妻精品电影 | av视频免费观看在线观看| 亚洲男人天堂网一区| 亚洲成av片中文字幕在线观看| 在线观看免费视频网站a站| www.自偷自拍.com| 女性被躁到高潮视频| 夜夜骑夜夜射夜夜干| 老司机影院成人| 亚洲专区国产一区二区| 欧美成人午夜精品| 777久久人妻少妇嫩草av网站| 一级a爱视频在线免费观看| 亚洲中文字幕日韩| 免费在线观看视频国产中文字幕亚洲 | 亚洲熟女精品中文字幕| 日韩大码丰满熟妇| 中文字幕人妻熟女乱码| 99久久国产精品久久久| 久久精品久久久久久噜噜老黄| 国产97色在线日韩免费| 爱豆传媒免费全集在线观看| 亚洲avbb在线观看| av天堂久久9| 久久久久久久精品精品| 中文字幕高清在线视频| a级毛片黄视频| 国产精品一区二区精品视频观看| 国产精品久久久久久精品电影小说| 日韩一卡2卡3卡4卡2021年| 不卡av一区二区三区| 欧美亚洲日本最大视频资源| 亚洲欧美精品自产自拍| 成年人免费黄色播放视频| a级片在线免费高清观看视频| 欧美午夜高清在线| 国产在线免费精品| 国产精品自产拍在线观看55亚洲 | 男女无遮挡免费网站观看| 热99国产精品久久久久久7| 亚洲精品一卡2卡三卡4卡5卡 | 欧美 亚洲 国产 日韩一| 午夜激情av网站| 久久久久久免费高清国产稀缺| 女人被躁到高潮嗷嗷叫费观| 纯流量卡能插随身wifi吗| 亚洲欧美色中文字幕在线| 国产一卡二卡三卡精品| 国产精品.久久久| 精品久久久久久电影网| 亚洲精品av麻豆狂野| 12—13女人毛片做爰片一| 一本一本久久a久久精品综合妖精| 五月开心婷婷网| 亚洲欧美精品自产自拍| 亚洲av欧美aⅴ国产| 男人操女人黄网站| 亚洲精品久久午夜乱码| 我要看黄色一级片免费的| 国产精品久久久久久精品古装| h视频一区二区三区| 国产精品av久久久久免费| 一区二区三区精品91| 男人操女人黄网站| 在线av久久热| 黄色怎么调成土黄色| 精品第一国产精品| 亚洲国产欧美网| 激情视频va一区二区三区| 国产野战对白在线观看| 少妇精品久久久久久久| 成人av一区二区三区在线看 | 19禁男女啪啪无遮挡网站| 人成视频在线观看免费观看| 后天国语完整版免费观看| 久久久国产一区二区| 国产免费视频播放在线视频| 久久久久久免费高清国产稀缺| 动漫黄色视频在线观看| 日本五十路高清| 十八禁网站网址无遮挡| 在线看a的网站| 99香蕉大伊视频| 亚洲欧洲精品一区二区精品久久久| 国产成人免费观看mmmm| 男女高潮啪啪啪动态图| 国产成人精品久久二区二区免费| 欧美 亚洲 国产 日韩一| 精品少妇一区二区三区视频日本电影| 超碰成人久久| 纯流量卡能插随身wifi吗| 国产av一区二区精品久久| 国产三级黄色录像| 麻豆av在线久日| 最新在线观看一区二区三区| 亚洲色图综合在线观看| 亚洲欧美精品综合一区二区三区| 侵犯人妻中文字幕一二三四区| 少妇精品久久久久久久| 在线天堂中文资源库| 男女下面插进去视频免费观看| 91国产中文字幕| www.熟女人妻精品国产| 9191精品国产免费久久| 动漫黄色视频在线观看| av网站免费在线观看视频| 亚洲精华国产精华精| tocl精华| 黑人猛操日本美女一级片| 一本大道久久a久久精品| 韩国精品一区二区三区| 另类亚洲欧美激情| 永久免费av网站大全| 性少妇av在线| 一级毛片精品| 亚洲五月色婷婷综合| 免费少妇av软件| 成人免费观看视频高清| 精品熟女少妇八av免费久了| 多毛熟女@视频| 欧美 亚洲 国产 日韩一| 欧美性长视频在线观看| 99国产极品粉嫩在线观看| 美女高潮喷水抽搐中文字幕| 嫁个100分男人电影在线观看| 久久久久久亚洲精品国产蜜桃av| 啦啦啦视频在线资源免费观看| 亚洲精品成人av观看孕妇| 亚洲精品国产区一区二| 在线观看免费视频网站a站| av网站免费在线观看视频| 久9热在线精品视频| 国产精品免费大片| 精品一区二区三区av网在线观看 | 国产1区2区3区精品| 国产精品国产三级国产专区5o| 婷婷色av中文字幕| 男人添女人高潮全过程视频| 亚洲中文av在线| 成在线人永久免费视频| 制服诱惑二区| 国产精品秋霞免费鲁丝片| 精品国产乱子伦一区二区三区 | 丝袜脚勾引网站| 美女午夜性视频免费| 大香蕉久久网| 麻豆国产av国片精品| 黄色怎么调成土黄色| 女人被躁到高潮嗷嗷叫费观| 免费少妇av软件| 国产成人精品在线电影| 久久女婷五月综合色啪小说| 欧美精品一区二区免费开放| 国产高清国产精品国产三级| 欧美 亚洲 国产 日韩一| www.自偷自拍.com| 国产精品久久久人人做人人爽| 男女之事视频高清在线观看| 精品熟女少妇八av免费久了| 亚洲情色 制服丝袜| 欧美日韩精品网址| 在线观看人妻少妇| 亚洲av电影在线进入| 在线观看免费午夜福利视频| 色视频在线一区二区三区| 建设人人有责人人尽责人人享有的| 考比视频在线观看| 高清黄色对白视频在线免费看| 国产精品亚洲av一区麻豆| 99九九在线精品视频| 午夜精品国产一区二区电影| 精品少妇黑人巨大在线播放| 国产精品偷伦视频观看了| 亚洲美女黄色视频免费看| 欧美精品av麻豆av| 777久久人妻少妇嫩草av网站| 精品亚洲成a人片在线观看| 欧美激情极品国产一区二区三区| 国产欧美日韩综合在线一区二区| 亚洲成人手机| 一区二区三区激情视频| 亚洲视频免费观看视频| 亚洲精华国产精华精| www.精华液| 久久精品亚洲熟妇少妇任你| 纯流量卡能插随身wifi吗| 久久久欧美国产精品| 亚洲精品自拍成人| 精品第一国产精品| 男男h啪啪无遮挡| 18禁观看日本| 国产成人精品在线电影| 午夜91福利影院| 人人妻,人人澡人人爽秒播| 日韩欧美免费精品| 成年女人毛片免费观看观看9 | 精品一区在线观看国产| 两性午夜刺激爽爽歪歪视频在线观看 | 国产一区二区在线观看av| 男男h啪啪无遮挡| 午夜福利视频在线观看免费| 天堂中文最新版在线下载| 亚洲综合色网址| xxxhd国产人妻xxx| 亚洲少妇的诱惑av| 性高湖久久久久久久久免费观看| 女人精品久久久久毛片| 久久久久久久国产电影| 国产色视频综合| 九色亚洲精品在线播放| 黄频高清免费视频| 精品高清国产在线一区| 深夜精品福利| 欧美少妇被猛烈插入视频| 久久国产精品男人的天堂亚洲| 三级毛片av免费| 老汉色∧v一级毛片| 久久久水蜜桃国产精品网| 一区在线观看完整版| 王馨瑶露胸无遮挡在线观看| 久久精品国产a三级三级三级| 美女高潮喷水抽搐中文字幕| 男人操女人黄网站| 多毛熟女@视频| 日韩 亚洲 欧美在线| 国产淫语在线视频| av网站免费在线观看视频| www.熟女人妻精品国产| 成年人午夜在线观看视频| 日本91视频免费播放| 男女免费视频国产| 成人三级做爰电影| 亚洲专区国产一区二区| 午夜福利,免费看| 少妇人妻久久综合中文| 国产精品久久久久成人av| 亚洲精品久久成人aⅴ小说| 久久久久国内视频| 人人妻人人添人人爽欧美一区卜| 久久人人97超碰香蕉20202| 日韩欧美免费精品| 成人影院久久| 成人av一区二区三区在线看 | 久久久水蜜桃国产精品网| 国产一区二区三区av在线| 我的亚洲天堂| 久久久久久人人人人人| 久久国产精品大桥未久av| 久久久久国产一级毛片高清牌| 黑人巨大精品欧美一区二区mp4| 在线天堂中文资源库| 首页视频小说图片口味搜索| 色老头精品视频在线观看| 亚洲,欧美精品.| 99热网站在线观看| 亚洲精品成人av观看孕妇| 国产精品二区激情视频| 99热网站在线观看| 精品福利观看| 久久精品aⅴ一区二区三区四区| 欧美成狂野欧美在线观看| 欧美精品啪啪一区二区三区 | 18在线观看网站| 日韩中文字幕欧美一区二区| 十八禁人妻一区二区| 九色亚洲精品在线播放| 久热这里只有精品99| 国产欧美日韩一区二区精品| 99久久99久久久精品蜜桃| 国产精品久久久久久精品古装| 日韩熟女老妇一区二区性免费视频| 国产国语露脸激情在线看| 亚洲国产欧美网| 老熟妇仑乱视频hdxx| 永久免费av网站大全| 黄色怎么调成土黄色| 男女午夜视频在线观看| 18在线观看网站| 99久久综合免费| 久久毛片免费看一区二区三区| 精品熟女少妇八av免费久了| 岛国在线观看网站| kizo精华| 中文字幕人妻丝袜制服| 成人影院久久| 精品亚洲乱码少妇综合久久| 一个人免费在线观看的高清视频 | 久久免费观看电影| 高潮久久久久久久久久久不卡| 精品人妻一区二区三区麻豆| 午夜福利,免费看| 成人黄色视频免费在线看| 国产99久久九九免费精品|