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

    石山苣苔屬(苦苣苔科)花形態(tài)演化及分類學(xué)意義

    2017-05-30 10:48:04盧永彬黃俞淞許為斌黃潔劉演向春雷張強(qiáng)
    廣西植物 2017年10期
    關(guān)鍵詞:分類

    盧永彬 黃俞淞 許為斌 黃潔 劉演 向春雷 張強(qiáng)

    摘要: 石山苣苔屬(苦苣苔科)約30種,主要分布于我國南部的石灰?guī)r地區(qū)。目前該屬已知物種數(shù)雖少但花形態(tài)極其多樣,是該科中分類最為困難的類群之一?;诜肿幼C據(jù),其它8個(gè)屬中花形態(tài)迥異的一些物種被并入石山苣苔屬。然而,該屬花形態(tài)的演化趨勢缺乏系統(tǒng)性的研究,傳統(tǒng)分類對(duì)屬的界定與分子系統(tǒng)學(xué)研究結(jié)果相矛盾的原因,以及是否有形態(tài)特征支持新界定的石山苣苔屬還不清楚。該研究中,總共編碼了19種石山苣苔屬植物和9種報(bào)春苣苔屬植物的35個(gè)形態(tài)特征,其中包括26個(gè)花部形態(tài)特征,在分子系統(tǒng)樹上追蹤了它們的演化路徑。結(jié)果表明:無論屬內(nèi)還是屬間,多數(shù)花部形態(tài)特征,尤其以往屬的分類界定特征,在演化過程中變化頻繁且發(fā)生了高度同塑性演化,這是導(dǎo)致傳統(tǒng)形態(tài)分類不自然的關(guān)鍵因素。此外,在觀察研究的所有特征中,花絲和柱頭的差異可能在石山苣苔屬植物共同祖先中經(jīng)歷了演變,或可用于區(qū)分石山苣苔屬與其姐妹報(bào)春苣苔屬的大多數(shù)種類。因此,在苦苣苔科植物的分類學(xué)研究中應(yīng)當(dāng)慎用這些花部性狀作為分類依據(jù),而且應(yīng)對(duì)形態(tài)特征進(jìn)行廣泛地觀察研究,在密集的取樣和分辨率更高、更可靠的系統(tǒng)樹上追蹤它們的演化規(guī)律。更為重要的是,需要進(jìn)一步研究導(dǎo)致復(fù)雜形態(tài)性狀演化的內(nèi)在分子調(diào)控機(jī)理和外在的自然選擇動(dòng)力,最終更加深入地理解石山苣苔屬等典型喀斯特植物的演化過程和機(jī)理。

    關(guān)鍵詞: 分類, 花部特征, 苦苣苔科, 同塑性演化, 喀斯特植物, 石山苣苔屬

    Petrocodon Hance (Gesneriaceae), firstly established in 1883, remained as a monotypic genus (only containing Petrocodon dealbatus) over a century. It is perennial herb and has campanulate corolla with indistinctly twolipped limb, being different from most other Gesneriaceae, which usually possess distinctly twolipped limb with lobes different in size and/or shape. The species is distributed in South and Central China, mostly growing in the crevices of limestone rocks. No more species of the genus were reported in the following century. Till the last decade, several new species of Petrocodon were finally discovered and described from South China (Wei, 2007; Wei et al, 2010).

    The phylogenetic position and the boundary between Petrocodon s.s. and other closelyrelated genera were revealed by the recent molecular phylogenetic studies (Weber et al, 2011b; Xu et al, 2014). The monotypic genera (Calcareoboea C. Y. Wu ex H. W. Li, Paralagarosolen Wei, Dolicholoma D. Fang & W. T. Wang and Tengia Chun) and one to a few species from Wentsaiboea Fang & Qin, Primulina Hance, Lagarosolen Wang and Didymocarpus Wallich were included into Petrocodon and the expanded genus was indicated to be most close to Primulina according to the molecular phylogenetic studies based on the chloroplast trnLtrnF and nuclear ribosome ITS DNA sequences (Weber et al, 2011b; Xu et al, 2014). The redefinition as well as more new species recently published enlarged Petrocodon to contain ca. 30 species (Weber et al, 2011b; Wen et al, 2012; Chen et al, 2014; Hong et al, 2014; Xu et al, 2014; Yu et al, 2015; Guo et al, 2016; Lu et al, 2017). It is more noteworthy is that the expanded genus represent one of the most varied genera in the Old World, comparable to those in the New World (Weber et al, 2011b).

    The morphological traits of the expanded Petrocodon, particularly those of floral parts, are remarkably diverse among the lineages of Old World Gesneriaceae. The floral variation of Petrocodon covers almost all aspects of the floral characters. For example, the symmetry of corolla includes zygomorphy and almost actinomorphy; most species have corolla limb with two adaxial (dorsal) lobes and three abaxial (ventral) lobes but two species of the genus have corolla limb with a very unusual pattern of four adaxial lobes and one abaxial lobe [corresponding to the terminology of fourtoothed upper lip and onetoothed lower lip according to Weber et al (2011b)]; corolla shape is variable from campanulate, tubular, funnelform, urceolate to salverform; corolla coloration is diverse, including pure white, purple, yellow and bright red; stamen numbers per flower are usually two but Pet. scopulorum (previous Tengia scopulorum) and one recently published species Pet. hunanensis (Yu et al, 2015) belonging to the genus possess five and four fertile stamens, respectively. The general appearance of the floral diversity of Petrocodon is as shown in Fig. 1.

    The redefinition of Petrocodon and the inclusion of species from other eight genera, could imply that morphological characters, particularly those diagnostic traits the traditional taxonomy was based on, might experience complicated evolution such as homoplasy (i.e. convergence, parallel or reversal evolution), dramatic morphological divergence of closely related species, and/or plesiomorphy remained in different lineages etc. However, little is known about the evolutionary trends of the diverse morphological characters in the expanded genus and their taxonomic significance, though the previous molecular phylogenetic study sampled the majority of the included species (Weber et al, 2011b; Chen et al, 2014). It remains unknown whether there is any floral or vegetative characters that can differentiate Petrocodon from its sister genus Primulina and other close relatives, and characterize each of the lineages in the molecular phylogenetic tree within the genus. To explore these issues needs carefully examination of the diverse morphological characters of Petrocodon, comparing with other close relatives, and tracing their evolution based on a robust phylogenetic framework suggested by molecular phylogenetic studies.

    In this study, we examined and coded 35 morphological characters (including 26 flora traits) of Petrocodon and some outgroup taxa of Primulina, traced the evolutionary trends of these characters on the phylogenetic tree reconstructed based on DNA sequences downloaded from the public accessible database of National Centre of Biological Information (NCBI). The aims we targeted to explore included: (1) the evolutionary trends how the diverse morphologies evolved; (2) the potential synapomorphy that can characterize and differentiate Petrocodon from its close relatives. We also discuss the implications of the inferred morphological evolution for future classification of Gesneriaceae.

    1Materials and Methods

    1.1 Materials examined

    A total of 35 morphological characters (including 26 floral characters) of 19 Petrocodon species and 9 outgroup species of Primulina were examined and coded numerically. The morphological characters were coded and validated according to the original description of the species, the related monographs such as Flora of China (Wang et al, 1998), Plants of Gesneriaceae in China (Li & Wang, 2005) and Gesneriaceae of South China (Wei et al, 2010), and the observation of transplanted individuals and/or the field investigations.

    We coded different states of the same character with different numerics. For those possessing mosaic character states, we coded them with all contained states. For example, for the corolla coloration, Pet. tiandengensis has corolla with pale blue lobes and white tube, we coded it with these two character states (0 and 2, representing bluish violet and white, respectively.) For the characters lacking of observation and record, we treated them as missing. All the coded traits are listed in Table 1.

    1.2 Selection of DNA markers and sampling strategy strategy

    As the plastid trnLtrnF intergenic spacer and nuclear ribosome internal transcribed spacer (ITS) sequences have been most frequently used for phylogenetic analysis in Gesneriaceae, particularly in the subfamily Didymocarpoideae including Petrocodon from the Old World (Denduangboripant et al, 2007; Wang et al, 2010; Mller et al, 2011a, b; Weber et al, 2011b; Chen et al, 2014), these two DNA fragments were therefore selected and downloaded from GenBank for the present study.

    As shown by some previous studies (Weber et al, 2011b; Chen et al, 2014), the newly redefined genus Primulina is the closest relative to Petrocodon, therefore, all the trnLtrnF and ITS sequences of Petrocodon taxa and a few of Primulina representatives available from GenBank were included. The sequences of some other Chinese Gesneriaceae, namely Allocheilos W. T. Wang, and Oreocharis mileensis which also possess the rare corolla character of four adaxial lobes and one abaxial lobe (Fig. 2) were also included. Some other Old World Gesneriaceae, e.g. Deinocheilos W. T. Wang and Gyrocheilos W. T. Wang with another pattern of special corolla structure (unity of the adaxial lip with normal tripartition of the abaxial lip) as well as several other Oreocharis species were also downloaded and included in the present phylogenetic analyses. Finally, sequences of some representative species of Lysionotus, Cyrtandra and Streptocarpus were included and used as the more distant outgroups according to the previous largerscale phylogenetic analyses (Roalson & Roberts, 2016). The trnLtrnF data included 60 accessions representing 45 species and the ITS data encompassed 54 accessions of 44 species. The GenBank accessions of all the downloaded sequences are as listed in Table S1 in supporting information.

    1.3 Sequence assembling and phylogenetic analyses

    All the sequences were aligned using the software MUSCLE 3.8.31 (Edgar, 2004), and then manually adjusted using Bioedit 5.0.9 (Hall, 1999). The flank regions which were ambiguously aligned or possessed numerous missing nucleotides for the majority of the accessions were deleted.

    The maximum likelihood (ML) in RaxMLVIHPC (Stamatakis, 2006), maximum parsimony (MP) in PAUP 4.10 (Swofford, 2002) and the Bayesian inference (BI) in MrBayes 3.2.6 (Ronquist et al, 2012) were employed to reconstruct the phylogenies. First, we conducted the ML analyses for each of the two regions under the settings of 1 000 rapid bootstrap searches and thereafter a thorough ML search with the substitution model GTR+G. The congruence of the phylogenetic signals of the two loci were assessed by comparing the two ML majorityrule consensus trees using a bootstrap support value of 70% as an arbitrary threshold, meaning that the nodes showing incongruent relationships between trnLtrnF and ITS ML majorityrule consensus trees with bootstrap support value above 70% were assumed to be strong phylogenetic conflict and against combined analysis of the two loci. For the ML analysis of the combined data, the same parameters were settled. For the MP analysis of the combined data, the MP settings were heuristic searches of 1 000 replicates of random sequence addition, tree bisection and reconnection (TBR) swapping and all MP trees saved at each replicate (Multree on); for the bootstrap analysis, 1 000 bootstrap pseudoreplicates were conducted, each with 10 replicates of random sequence addition and a maximum of 5 000 trees saved for each bootstrap pseudoreplicate. For the BI analysis, 100 000 000 generations were run with four chains in two parallel runs and one tree every 5 000 generations were sampled with a burnin of the first 5 000 trees discarded. The convergence of the two parallel runs was guaranteed by the splitting frequency less than 0.005. All other parameters were set as default.

    We reconstructed the ancestral floral morphologies and traced the evolutionary trends based on the molecular phylogenetic tree generated from the combined data with only Petrocodon and Primulina taxa preserved using Fitch maximum parsimony implemented in Mesquite 2.01 (Maddison & Maddison, 2007). The phylogenetic tree to which the morphological characters were mapped was modified to comprise only one representative accession per species, with additional repetitive accessions of the same species pruned. All character states were treated as unordered, with all other settings left as default.

    2Results and Analysis

    2.1 Phylogenetic analysis of Petrocodon and its allies

    The aligned trnLtrnF and ITS datasets were 896 and 744 base pairs including 74 and 286 informative sites, and 72 and 123 variable but parsimony uninformative sites, respectively. The combined matrix thus consisted of 1 640 characters including 555 variable sites with 360 informative sites. The parameters of consistency index (CI), retention index (RI) and homoplasy index (HI) were 0.836, 0.904, and 0.164, and 0.553, 0.736 and 0.447 for the trnLtrnF and ITS data, respectively.

    For the combined data, the ML, BI and MP analyses yielded 50% majorityrule consensus trees with congruent topologies (Fig. 3). The monophyly of Petrocodon was recovered with moderate to high support values (BSML = 65%; PP= 0.98; BSMP = 100%), and the genus was fully supported to be the sister of Primulina (BSML = 100%; PP= 1.00; BSMP = 100%), which are in line with previous phylogenetic analyses (Weber et al, 2011b; Xu et al, 2014; Guo et al, 2016; Lu et al, 2017). The relationships within Petrocodon are as shown in Fig. 3.

    2.2 Floral evolution in Petrocodon

    Mapping morphological characters to the molecular phylogeny clearly indicated most of the floral characters evolved frequently and highly homoplasiously in Petrocodon. The floral symmetry evolved twice in parallel from zygomorphy into almost actinomorphy (Fig. S1 in supporting information); the ancestral corolla shape

    Fig. 1Floral diversity of Petrocodon

    A. Petrocodon tiandengensis; B. Pet. pseudocoriaceifolius; C. Pet. longgangensis; D. Pet. lui; E. Pet. guangxiensis; F. Pet. dealbatus; G. Pet. laxicymosus; H. Pet. coccineus; I. Pet. coriaceifolius; J. Pet. niveolanosus; K. Pet. hancei; L. Pet. hechiensis; M. Pet. jingxiensis; N. Pet. fangianus; O. Pet. jasminiflorus.

    Fig. 2Parallel evolution of four adaxial corolla lobes and one abaxial lobe in separate Gesneriaceae lineages

    A. Petrolodon retroflexus; B. Pet. coccineus; C. Oreocharis mileensis; D. Allocheilos cortusiflorus.

    Fig. 3Maximum likelihood phylogenetic tree based on the combined trnLtrnF and ITS data, showing evolutionary relationships of Petrocodon and its allies

    ML bootstrap support, BI posterior probability and MP bootstrap support values are listed to the corresponding nodes. * stands for the bootstrap support or posterior probability equal to 100% or 1.00; - stands for the bootstrap support or posterior probability below 50% or 0.50.

    Fig. 4Evolutionary trends of four key representative floral characters, showing highly homoplasious evolution of these characters

    Fig. 5Two species pairs which are morphologically similar with each other (especially for the floral characters) but are genetically distant

    A-B. Primulina renifolia; C-D. Petrocodon tiandengensis; E. Primulina tabacum; F. Petrocodon guangxiensis.

    was inferred to be tubular and evolved frequently (a total of 13 times of changes) among variable states, including twice convergently into campanulate and thrice into funnelform; the floral coloration was inferred to be purple ancestrally and experienced 13 times of changes, including six independent transformations in parallel into white; even the very unusual corolla structure of four adaxial lobes and one abaxial lobe evolved twice in parallel from the ancestral and common state of two adaxial lobes and three abaxial lobes in Petrocodon, with additional same two shifts in other Old World Gesneriaceae (Allocheilos and Oreocharis mileensis) also from southern China; the stigma also experienced eight changes, seven from bifid into unity (being solitary) in parallel and one reversal back from unity into bifid (Fig. 4). Most other floral characters were also suggested to experience multiple homoplasious changes (data not shown).

    3Discussion

    3.1 Frequent changes, convergent and parallel evolution of the previously assumed taxonomysignificant floral characters

    Mapping the morphological characters to the phylogenetic tree clearly showed that most of the floral characters of Petrocodon experienced complicated evolution, i.e. frequent changes and multiple convergences and parallel changes, within and beyond the genus (Fig. 4). This seems especially true for those diagnostic characters that were used to delimit the taxonomic units, which could have misled the traditional taxonomic treatments of the genera.

    Floral (corolla) symmetry varied in Gesneriaceae and was assumed to be one of taxonomysignificant characters. Actinomorphic floral corolla symmetry, a relative rare state possessed only by a few Old World Gesneriaceae, was assumed to be a diagnostic character and mainly used to differentiate Primulina tabacum from other close relatives, e.g. the previous genus Chirita. Later, another morphologically similar new species Pri. guangxiensis which has general floral appearance and also almost actinomorphic corolla like Primulina tabacum was described and added to the monotypic genus (Liu et al, 2011). However, the molecular phylogenetic analyses clearly indicated that Pri. tabacum was deeply embedded in the redefined genus Primulina while Pri. guangxiensis was nested within the redefined genus Petrocodon, strongly suggesting parallel evolution of floral symmetry from zygomorphy to actinomorphy. Some other Gesneriaceae in China (e.g. the previous Thamnocharis Wang, Bournea Oliv. and Conandron Sieb. & Zucc.) were also characterized by actinomorphic corolla, whereas, these similarities were suggested to be resulted from morphological convergence (factually parallel evolution) rather than common ancestry according to the molecular phylogenetic analyses (Wang et al, 2010). Even in Petrocodon, the floral symmetry transformation from zygomorphy into almost actinomorphy was also suggested to have happen twice independently, leading to Petrocodon guangxiensis and Pet. scopulorum, respectively. All these suggest that the floral symmetry can evolve rapidly and in parallel and it is not suitable to be used as the diagnostic character to delimit taxonomic units (i.e. genus) in Gesneriaceae.

    The corolla shape, a previously assumed taxonomysignificant character, also evolved in parallel and convergence frequently. The previous Lagarosolen and Paralagarosolen were characterized by narrow and long (i.e. salverform) corolla tube. However, the present ancestral morphological reconstruction indicated that this salverform corolla evolved at least twice in parallel in Petrocodon from ancestral tubular, and it also evolved into urceolate and tubular quickly in different lineages. The campanulate corolla that characterized Petrocodon s.s. and Wentsaiboea was also suggested to evolve twice convergently within the expanded genus Petrocodon, and this convergence also happened between Petrocodon and Primulina. These suggested that the character of corolla shape (salverform and campanulate) that characterized Lagarosolen, Paralagarosolen, Wentsaiboea and Petrocodon s.s. also evolved homoplasiously in Petrocodon and its allies.

    Other characters that were assumed to be characteristic of the previous genera Calcareboea, Dolicholoma and Didymocarpus, whose taxa have been wholly or partially included into Petrocodon, were also suggested to be homoplasious according to the present study. Dolicholoma, a previous monotypic genus, possesses narrow corolla tube like the taxa of the previous Lagarosolen, but differs from the latter by its deep split corolla lobes. However, this character of deep split lobes was also observed in another recently described new species Petrocodon ainsliifolius. According to the previous molecular phylogenetic analysis (Chen et al, 2014), these two species are distantly related to each other, embedding within lineages with quite different characters, respectively, and thus suggesting convergent evolution of the character (discussed also in Weber et al, 2011). Calcareboea, a previous monotypic genus, is characteristic of bright red flower and very special corolla structure of four shallowly divided adaxial lobes and one abaxial lobe. Although the bright red flower is unique to this species in Petrocodon and was indicated to be derived from ancestral bluish violet, the red flower which was presumed to be associated with ornithophily (Weber et al, 2011b) is also possessed by some taxa of Aeschynanthus, another distantlyrelated genus of Gesneriaceae. The other characteristic of four adaxial lobes was also clearly indicated to have evolved in parallel within Petrocodon (at least twice independently) and beyond the genus (additionally twice observed leading to Allocheilos and Oreocharis mileensis). The frequent changes and highly homoplasious evolution of these previously assumed taxonomysignificant floral characters seemed to be the key factors that accounted for the unnaturalness of the previous generic classifications, resulted in taxonomic oversplit of genera.

    3.2 Potential synapomorphy that may characterize and differentiate the redefined Petrocodon from its closest relative Primulina

    Although the redefined genus forms as monophyletic, separated from its closest relative Primulina, little is known about the synapomorphies that can support the new classification. Particularly, these newly redefined sister genera not only encompass remarkable diversity, but also include morphologically similar species reclassified into Petrocodon and Primulina, respectively. For instance, the two morphologically similar taxa, Pri. tabacum and the previous Pri. guangxiensis (namely Petrocodon guangxiensis) were indicated to be embedded in Primulina and Petrocodon, respectively, according to molecular phylogenetic analysis (Xu et al, 2014). Also, Wentsaiboea renifolia, Wen. luochenensis and Wen. tiandenensis were once placed in the same genus by very similar gross vegetative and floral appearance, however, the molecular phylogenetic studies clearly suggested the former two species should be placed in Primulina, while the latter was nested within Petrocodon (Wang et al, 2011; Weber et al, 2011b). These morphologically similar taxa divided into two genera confounded the boundary and brought difficulties to uncover potential diagnostic morphological characters that can distinguish and characterize each of the genera.

    After extensive comparisons of Petrocodon with Primulina taxa, we noticed two possible synapomorphic characters that may be used to differentiate these two genera. Firstly, most Primulina taxa have geniculate filaments of the two fertile stamens, with an exception of Pri. tabacum, which is embedded deeply within the genus and possesses almost straight filaments. This character of geniculate filaments is also shared by some other Gesneriaceae lineages like Hemiboea and Lysionotus taxa, therefore it at least can be inferred to be ancestral of the genus Primulina as a whole, though it remains unknown weather this character shared by different Gesneriaceae lineages (genera) is indicative of plesiomorphy or convergent or parallel evolution. In contrast, most Petrocodon taxa we observed have straight or slightly curved filaments of the fertile stamens, being different from the commonly geniculate filaments of Primulina. However, Petrocodon retroflexus we recently published has special filaments, which are geniculate near the base like most Primulina taxa but with an additional loop near the tip where the anthers attached (Guo et al, 2016). Therefore, this character seemed also to be homoplasious between these two genera though it may differentiate most Petrocodon from Primulina taxa.

    Secondly, Primulina taxa usually have obtrapeziformlike stigma (Chiritoidtype stigma), which is commonly bifid. This type of stigma was assumed to be developed from the lower half of an ancestral biparted stigma, with the upper half reduced (Weber et al, 2011a). This characteristic of stigma is also shared by some other distantlyrelated lineages, particularly those (e.g. Microchirita, Henckelia, Liebigia, etc.) that were unnaturally placed in the previous genus Chirita (Weber et al, 2011a). Petrocodon taxa, however, have different stigma types which usually have undivided or bipartitioned capitate stigma(s) [corresponding to the concept of two stigmas in Wang et al (1998)] and may be a derived state in the ancestor of the lineage leading to Petrocodon. The capitatelike versus obtrapeziformlike stigma may be used to distinguish Petrocodon from Primulina.

    Anyway, these tentatively proposed characteristics of Petrocodon versus Primulina need to be tested by more future extensive examinations of these characters in these two genera and the outgroups and inferences of their evolutionary trends based on more robust phylogenetic tree with denser sampling of the taxa.

    4Conclusions and Prospects

    The present study illustrated complicated morphological evolution in the genus Petrocodon, especially for those floral characters that were assumed to be taxonomysignificant. Frequent changes and highly convergent and parallel evolution of these traits in the genus and its close relatives seemed to be key factors that have misled the traditional taxonomy, and brought difficulties to identify possible synapomorphic characters that can be used to support the new classifications based on molecular phylogenetic studies. More extensive investigations should be carried out on both morphologies and molecular phylogenetic studies. On one side, the available phylogenetic studies in Gesneriaceae were mainly based on limited plastid and the nuclear ribosome ITS DNA sequences (Wang et al, 2011; Weber et al, 2011a, b), which sometimes were likely to yield low resolution or partial or even misleading relationships. So that, denser samplings of DNA sequences from different genomes as well as taxa are required to reconstruct the highly resolved species tree. On the other side, more morphologies including micromorphological and anatomical characters, etc. should be examined and added to the present data set, and their evolutionary trends should be traced on robust phylogenetic framework. More importantly, efforts should also be taken to investigate the underlying causes (e.g. molecular regulatory network from the discipline of evolutionary developmental genetics) and its biological senses (e.g. the fitness adapted to natural selection) of the complicated morphological evolution. All these together would help bridge up the gaps between traditional taxonomy and molecular phylogenetic studies and between phenotypes and genotypes, and advance the understandings of the evolution of this group of typical karst flora and its allies.

    5Supporting Information

    Table S1 GenBank accessions of the sequences used in the present study.

    Fig. S1 Evolutionary trends of four selected flower traits, showing highly homoplasious evolution. The ancestral states of the traits were traced using the maximum parsimony method implemented in Mesquite 2.7.

    AcknowledgmentsWe thank Dong Yang and LI Peng Wei from Institute of Botany, Chinese Academy of Sciences for reading the manuscript and providing helpful suggestions. We also send our gratitude to an anonymous reviewer of the manuscript for insightful suggestions.

    參考文獻(xiàn):

    CHEN WH, MLLER M, SHUI YM, et al, 2014. Three new species of Petrocodon (Gesneriaceae), endemic to the limestone areas of Southwest China, and preliminary insights into the diversification patterns of the genus [J]. Syst Bot, 39(1): 316-330.

    DENDUANGBORIPANT J, CRONK QCB, KOKUBUGATA G, et al, 2007. Variation and inheritance of nuclear ribosomal DNA clusters in Streptocarpus (Gesneriaceae) and their biological and phylogenetic implications [J]. Int J Plant Sci, 168(4): 455-467.

    EDGAR RC, 2004. MUSCLE: multiple sequence alignment with high accuracy and high throughput [J]. Nucl Acid Res, 32(5): 1792-1797.

    GUO J, MENG T, PANG H, et al, 2016. Petrocodon retroflexus sp. nov. (Gesneriaceae) from a karst cave in Guizhou, China [J]. Nord J Bot, 34: 159-164.

    HALL TA, 1999. BioEdit: a userfriendly biological sequence alignment editor and analysis program for Windows 95/98/NT [C]. Nucl Acids Symp Ser, 41: 95-98.

    HONG X, ZHOU SB, WEN F, 2014. Petrocodon villosus (Gesneriaceae), a new species from Guangxi, China [J]. Blumea, 59(1): 33-36.

    LI ZY, WANG YZ, 2005. Plants of Gesneriaceae in China [M]. Zhengzhou: Henan Science and Technology Publishing House.

    LIU Y, XU WB, HUANG YS, 2011. Primulina guangxiensis sp. nov. (Gesneriaceae) from a karst cave in Guangxi, China [J]. Nord J Bot, 29(6): 682-686.

    LU YB, PAN B, HUANG J, et al, 2017. Petrocodon pulchriflorus sp. nov. (Gesneriaceae) from Guangxi, China [J]. Phytotaxa, 291(1): 81-88.

    MLLER M, FORREST A, WEI YG, et al, 2011a. A molecular phylogenetic assessment of the advanced Asiatic and Malesian didymocarpoid Gesneriaceae with focus on nonmonophyletic and monotypic genera [J]. Plant Syst Evol, 292(3-4): 223-248.

    MLLER M, MIDDLETON D, NISHII K, et al, 2011b. A new delineation for Oreocharis incorporating an additional ten genera of Chinese Gesneriaceae [J]. Phytotaxa, 23(1): 1-36.

    MLLER M, PFOSSER M, JANG CG, et al, 2009. A primulinary phylogeny of the ‘Didymocarpoid Gesneriaceae based on three molecular data sets: incongruence with available tribal classification [J]. Am J Bot, 96(5): 989-1010.

    MADDISON WP, MADDISON DR, 2007. Mesquite: A modular system for evolutionary analysis. Version 2.01 [DB/OL]. Available from: http://mesquiteproject.org/mesquite/mesquite.html.

    POSADA D, CRANDALL KA, 1998. Modeltest: testing the model of DNA substitution [J]. Bioinformatics, 14(9): 817-818.

    ROALSON EH, ROBERTS WR, 2016. Distinct processes drive diversification in different clades of Gesneriaceae [J]. Syst Biol, 65(4): 662.

    RONQUIST F, TESLENKO M, VAN DER MARK P, et al, 2012. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space [J]. Syst Biol, 61(3): 539-542.

    STAMATAKIS A, 2006. RAxMLVIHPC: maximum likelihoodbased phylogenetic analyses with thousands of taxa and mixed models [J]. Bioinformatics, 22(21): 2688-2690.

    SWOFFORD DL, 2002. PAUP*. Phylogenetic analysis using parsimony (*and other methods). Version 4 [CP]. Sinauer Associates, Sunderland, Massachusetts.

    WANG WT, PAN KY, LI ZY, et al, 1998. Gesneriaceae [M]//WU ZY, RAVEN PH. Flora of China: Beijing: Science Press; St. Louis: Missouri Botanical Garden Press.

    WANG YZ, LIANG RH, WANG BH, et al, 2010. Origin and phylogenetic relationships of the Old World Gesneriaceae with actinomorphic flowers inferred from ITS and trnLtrnF sequences [J]. Taxon, 59(4): 1044-1052.

    WANG YZ, MAO RB, LIU Y, et al, 2011. Phylogenetic reconstruction of Chirita and allies (Gesneriaceae) with taxonomic treatments [J]. J Syst Evol, 49(1): 50-64.

    WEBER A, MIDDLETON DJ, FORREST A, et al, 2011a. Molecular systematics and remodelling of Chirita and associated genera (Gesneriaceae) [J]. Taxon, 60(3): 767-790.

    WEBER A, WEI YG, PUGLISI C, et al, 2011b. A new definition of the genus Petrocodon (Gesneriaceae) [J]. Phytotaxa, 23(1): 49-67.

    WEI YG, 2007. Petrocodon ferrugineus (Gesneriaceae), a new species from Guangxi, China [J]. Novon, 17(1): 135-137.

    WEI YG, WEN F, MLLER M, et al, 2010. Gesneriaceae of South China [M]. Guangxi Science and Technology Publishing House, Nanning: 1-777.

    WEN F, LIANG GY, WEI YG, 2012. Petrocodon lancifolius (Gesneriaceae), a new species endemic to a central subtropical zone of Guizhou Province, China [J]. Phytotaxa, 49(1): 45-49.

    XU WB, MENG T, ZHANG Q, et al, 2014. Petrocodon (Gesneriaceae) in the limestone karsts of Guangxi, China: three new species and a new combination based on morphological and molecular evidence [J]. Syst Bot, 39(3): 965-974.

    YU XL, LI M, ZHOU JJ, et al, 2015. Petrocodon hunanensis (Gesneriaceae), a new species identified by both morphological and molecular evidence from limestone area in Hunan, China [J]. Phytotaxa, 195(1): 65-72.

    猜你喜歡
    分類
    2021年本刊分類總目錄
    分類算一算
    垃圾分類的困惑你有嗎
    大眾健康(2021年6期)2021-06-08 19:30:06
    星星的分類
    我給資源分分類
    垃圾分類,你準(zhǔn)備好了嗎
    分類討論求坐標(biāo)
    數(shù)據(jù)分析中的分類討論
    按需分類
    教你一招:數(shù)的分類
    观看av在线不卡| 日韩中文字幕视频在线看片| 9色porny在线观看| 国产日韩欧美视频二区| 亚洲精品一二三| 人成视频在线观看免费观看| 狂野欧美激情性xxxx| 啦啦啦中文免费视频观看日本| 亚洲熟女毛片儿| 亚洲精品久久成人aⅴ小说| 少妇裸体淫交视频免费看高清 | 在线看a的网站| 国产精品秋霞免费鲁丝片| 亚洲欧美日韩高清在线视频 | 欧美国产精品va在线观看不卡| 亚洲av在线观看美女高潮| 亚洲成av片中文字幕在线观看| 人成视频在线观看免费观看| 啦啦啦在线观看免费高清www| 午夜日韩欧美国产| 午夜免费成人在线视频| 色网站视频免费| 最黄视频免费看| 国产精品一区二区在线不卡| 国产日韩欧美视频二区| 99国产精品免费福利视频| 国产黄色免费在线视频| 国产日韩一区二区三区精品不卡| 天堂俺去俺来也www色官网| 亚洲视频免费观看视频| 久久精品亚洲av国产电影网| 久久久久久久久免费视频了| 蜜桃国产av成人99| 亚洲成av片中文字幕在线观看| 在线观看免费日韩欧美大片| 欧美日韩视频高清一区二区三区二| 亚洲成人免费电影在线观看 | 亚洲中文日韩欧美视频| 新久久久久国产一级毛片| 性高湖久久久久久久久免费观看| 性少妇av在线| 丝袜美足系列| av在线app专区| 精品久久久久久久毛片微露脸 | 久久久国产一区二区| 老司机深夜福利视频在线观看 | 久久人妻熟女aⅴ| 国产精品亚洲av一区麻豆| 欧美日韩国产mv在线观看视频| www.自偷自拍.com| 亚洲九九香蕉| 亚洲免费av在线视频| 99久久99久久久精品蜜桃| 久久精品国产a三级三级三级| 久久久国产精品麻豆| 国产极品粉嫩免费观看在线| 亚洲欧美一区二区三区久久| 精品亚洲成国产av| 国产免费福利视频在线观看| 日韩av免费高清视频| 色婷婷久久久亚洲欧美| 欧美变态另类bdsm刘玥| 在线观看免费日韩欧美大片| 黄色怎么调成土黄色| 老汉色∧v一级毛片| 精品少妇黑人巨大在线播放| 王馨瑶露胸无遮挡在线观看| 亚洲精品成人av观看孕妇| 丝瓜视频免费看黄片| 十分钟在线观看高清视频www| 日韩一卡2卡3卡4卡2021年| 亚洲成国产人片在线观看| 免费在线观看日本一区| 最近手机中文字幕大全| 一级a爱视频在线免费观看| 久久人人97超碰香蕉20202| 国产一区二区激情短视频 | √禁漫天堂资源中文www| 下体分泌物呈黄色| 熟女av电影| 日本av免费视频播放| 精品少妇久久久久久888优播| 黄色 视频免费看| 久久鲁丝午夜福利片| 中文字幕人妻熟女乱码| av电影中文网址| 搡老岳熟女国产| 欧美变态另类bdsm刘玥| 国产黄色免费在线视频| 国产成人免费无遮挡视频| 伊人久久大香线蕉亚洲五| 午夜福利在线免费观看网站| 国产一区二区在线观看av| 人妻人人澡人人爽人人| 精品视频人人做人人爽| 亚洲一区二区三区欧美精品| 国产99久久九九免费精品| 欧美 日韩 精品 国产| 黄色视频不卡| 日韩大片免费观看网站| 一级黄片播放器| 久久av网站| 亚洲熟女精品中文字幕| 9191精品国产免费久久| 日韩电影二区| 免费看十八禁软件| 一区二区三区四区激情视频| 热99国产精品久久久久久7| 丝袜喷水一区| 国产熟女欧美一区二区| 国产成人一区二区三区免费视频网站 | 777久久人妻少妇嫩草av网站| 久热这里只有精品99| 一区二区三区精品91| 日本av手机在线免费观看| 国产爽快片一区二区三区| 久久久国产精品麻豆| 国产成人欧美在线观看 | 最黄视频免费看| 精品国产国语对白av| 国产精品麻豆人妻色哟哟久久| 大话2 男鬼变身卡| 男女国产视频网站| 999久久久国产精品视频| 大陆偷拍与自拍| 十分钟在线观看高清视频www| 亚洲一区中文字幕在线| 精品少妇黑人巨大在线播放| 一区二区av电影网| 成年人午夜在线观看视频| 免费看不卡的av| 亚洲熟女精品中文字幕| 我要看黄色一级片免费的| 亚洲国产看品久久| 久久九九热精品免费| 久久 成人 亚洲| 欧美av亚洲av综合av国产av| 国产成人精品久久二区二区91| 中文字幕最新亚洲高清| 国产在线一区二区三区精| 成年美女黄网站色视频大全免费| 丰满人妻熟妇乱又伦精品不卡| 国产亚洲一区二区精品| 精品国产一区二区三区久久久樱花| 成人三级做爰电影| 99热全是精品| 免费高清在线观看日韩| 男女高潮啪啪啪动态图| 久久久久久人人人人人| 国产成人精品久久二区二区91| 亚洲av电影在线观看一区二区三区| 欧美 日韩 精品 国产| 欧美日韩亚洲综合一区二区三区_| 大片电影免费在线观看免费| 国产日韩欧美视频二区| 国产高清不卡午夜福利| 国产高清videossex| 每晚都被弄得嗷嗷叫到高潮| 色网站视频免费| 亚洲欧洲日产国产| 亚洲国产欧美日韩在线播放| 你懂的网址亚洲精品在线观看| 男女边吃奶边做爰视频| 亚洲 国产 在线| 国产免费福利视频在线观看| 亚洲国产精品一区三区| 精品亚洲成国产av| 制服诱惑二区| av不卡在线播放| 热99久久久久精品小说推荐| 日本一区二区免费在线视频| 老司机亚洲免费影院| 欧美变态另类bdsm刘玥| 亚洲免费av在线视频| 97精品久久久久久久久久精品| 国产一区二区三区综合在线观看| 国产精品亚洲av一区麻豆| 亚洲精品第二区| 久久精品人人爽人人爽视色| 又粗又硬又长又爽又黄的视频| 日韩中文字幕欧美一区二区 | 国产精品一区二区精品视频观看| 国产欧美日韩精品亚洲av| 老司机午夜十八禁免费视频| 老司机在亚洲福利影院| 亚洲av电影在线进入| 2018国产大陆天天弄谢| 韩国精品一区二区三区| 又大又爽又粗| 欧美 亚洲 国产 日韩一| 18禁国产床啪视频网站| 精品欧美一区二区三区在线| 丰满饥渴人妻一区二区三| 99精国产麻豆久久婷婷| xxxhd国产人妻xxx| 午夜福利一区二区在线看| 日本午夜av视频| 久久精品亚洲熟妇少妇任你| 亚洲久久久国产精品| 久久久久久亚洲精品国产蜜桃av| 99热全是精品| 久久99一区二区三区| 色视频在线一区二区三区| 精品亚洲成国产av| 亚洲人成电影免费在线| 日韩中文字幕欧美一区二区 | 亚洲熟女毛片儿| 人体艺术视频欧美日本| 看免费av毛片| 热99久久久久精品小说推荐| 亚洲伊人久久精品综合| 婷婷色av中文字幕| 亚洲七黄色美女视频| 欧美日韩视频高清一区二区三区二| av在线老鸭窝| 日韩伦理黄色片| 日韩 亚洲 欧美在线| videosex国产| 欧美久久黑人一区二区| 美女中出高潮动态图| 国产欧美日韩综合在线一区二区| 一边摸一边抽搐一进一出视频| www.av在线官网国产| 熟女少妇亚洲综合色aaa.| 男人舔女人的私密视频| 国产色视频综合| 不卡av一区二区三区| 十八禁网站网址无遮挡| 国产高清videossex| 国产野战对白在线观看| 麻豆国产av国片精品| 国产片特级美女逼逼视频| 成年人免费黄色播放视频| 日本猛色少妇xxxxx猛交久久| 成人国产一区最新在线观看 | 欧美黑人精品巨大| 一级毛片 在线播放| 国产亚洲一区二区精品| 一级,二级,三级黄色视频| 一级毛片我不卡| 久久av网站| 又粗又硬又长又爽又黄的视频| 狠狠精品人妻久久久久久综合| 欧美精品一区二区大全| 一本久久精品| 成人亚洲精品一区在线观看| 国产视频首页在线观看| 在线观看人妻少妇| 国产国语露脸激情在线看| 亚洲欧美精品综合一区二区三区| 午夜福利在线免费观看网站| 国产精品一区二区免费欧美 | 日韩人妻精品一区2区三区| 色播在线永久视频| 老汉色av国产亚洲站长工具| 国产精品三级大全| 亚洲少妇的诱惑av| 精品国产乱码久久久久久小说| 一二三四在线观看免费中文在| 久久国产精品大桥未久av| 国产日韩一区二区三区精品不卡| 国产免费又黄又爽又色| 日韩大片免费观看网站| 亚洲少妇的诱惑av| 午夜老司机福利片| 女人被躁到高潮嗷嗷叫费观| 中国美女看黄片| 一本久久精品| 夫妻性生交免费视频一级片| 黑人欧美特级aaaaaa片| 1024香蕉在线观看| av一本久久久久| 天天躁狠狠躁夜夜躁狠狠躁| 日本一区二区免费在线视频| 亚洲中文字幕日韩| 精品福利永久在线观看| 国产老妇伦熟女老妇高清| 亚洲av日韩精品久久久久久密 | 久久精品aⅴ一区二区三区四区| av片东京热男人的天堂| 亚洲图色成人| 女警被强在线播放| 亚洲精品国产av蜜桃| 又大又爽又粗| 一二三四在线观看免费中文在| 肉色欧美久久久久久久蜜桃| 国产1区2区3区精品| 国产欧美日韩一区二区三 | 高潮久久久久久久久久久不卡| 一区福利在线观看| 一本—道久久a久久精品蜜桃钙片| 欧美 日韩 精品 国产| 国产成人影院久久av| 亚洲精品日本国产第一区| 亚洲人成77777在线视频| 脱女人内裤的视频| 一二三四社区在线视频社区8| 久久久久精品国产欧美久久久 | 亚洲久久久国产精品| 男女国产视频网站| 另类亚洲欧美激情| 夫妻午夜视频| 久久精品aⅴ一区二区三区四区| 一本—道久久a久久精品蜜桃钙片| 亚洲专区中文字幕在线| 中文字幕人妻熟女乱码| 中文字幕制服av| 91国产中文字幕| 久久久久国产精品人妻一区二区| 天堂中文最新版在线下载| 王馨瑶露胸无遮挡在线观看| 国产欧美亚洲国产| 久久精品久久精品一区二区三区| 蜜桃在线观看..| 一边亲一边摸免费视频| 国产淫语在线视频| 日韩av免费高清视频| 亚洲欧美色中文字幕在线| 国产一区二区在线观看av| 在线亚洲精品国产二区图片欧美| 91国产中文字幕| 亚洲精品成人av观看孕妇| 电影成人av| 久久久久网色| 人成视频在线观看免费观看| 免费日韩欧美在线观看| 精品国产乱码久久久久久男人| 精品少妇内射三级| videosex国产| 免费一级毛片在线播放高清视频 | 国产91精品成人一区二区三区 | 国产亚洲欧美精品永久| 国产xxxxx性猛交| 欧美黄色淫秽网站| 欧美精品高潮呻吟av久久| 国产成人免费观看mmmm| av在线app专区| 精品一区二区三区四区五区乱码 | 两性夫妻黄色片| 国产日韩欧美视频二区| 欧美精品一区二区免费开放| 久久久精品区二区三区| 日本色播在线视频| 午夜视频精品福利| 国产成人一区二区三区免费视频网站 | 91精品三级在线观看| 国产高清videossex| 精品人妻在线不人妻| 午夜激情久久久久久久| 天天添夜夜摸| 中文字幕色久视频| 亚洲国产中文字幕在线视频| 亚洲三区欧美一区| 亚洲一码二码三码区别大吗| 熟女少妇亚洲综合色aaa.| 99久久人妻综合| 国产成人影院久久av| 欧美 亚洲 国产 日韩一| 在线精品无人区一区二区三| 精品一区二区三区四区五区乱码 | 97在线人人人人妻| 免费人妻精品一区二区三区视频| 大片免费播放器 马上看| 亚洲第一av免费看| 日本色播在线视频| 精品少妇一区二区三区视频日本电影| 岛国毛片在线播放| 国产亚洲欧美精品永久| 妹子高潮喷水视频| videosex国产| www.熟女人妻精品国产| 黑人巨大精品欧美一区二区蜜桃| 亚洲黑人精品在线| 日日摸夜夜添夜夜爱| avwww免费| 精品国产一区二区久久| 亚洲欧美日韩另类电影网站| 亚洲九九香蕉| 久久青草综合色| 两个人免费观看高清视频| 操美女的视频在线观看| 建设人人有责人人尽责人人享有的| 一级毛片黄色毛片免费观看视频| 欧美日韩视频高清一区二区三区二| 久久国产精品男人的天堂亚洲| 99精品久久久久人妻精品| 欧美日韩精品网址| 精品国产乱码久久久久久男人| 午夜免费鲁丝| 1024香蕉在线观看| 亚洲精品久久午夜乱码| 大香蕉久久网| 午夜91福利影院| 国产熟女欧美一区二区| 亚洲国产精品成人久久小说| 电影成人av| 秋霞在线观看毛片| 亚洲天堂av无毛| 老司机深夜福利视频在线观看 | 大片电影免费在线观看免费| h视频一区二区三区| 亚洲欧洲日产国产| 一本一本久久a久久精品综合妖精| 高潮久久久久久久久久久不卡| 九草在线视频观看| 美国免费a级毛片| 黄频高清免费视频| 777久久人妻少妇嫩草av网站| 中国美女看黄片| 亚洲精品在线美女| 人人妻人人澡人人看| 97在线人人人人妻| 可以免费在线观看a视频的电影网站| 国产又爽黄色视频| 国产高清videossex| 人妻一区二区av| 亚洲 欧美一区二区三区| 女人久久www免费人成看片| 色播在线永久视频| 国产成人精品无人区| 777久久人妻少妇嫩草av网站| 这个男人来自地球电影免费观看| 婷婷色av中文字幕| 极品人妻少妇av视频| 丁香六月天网| 一边亲一边摸免费视频| 欧美黄色片欧美黄色片| 亚洲中文日韩欧美视频| 在线看a的网站| 久热爱精品视频在线9| 中文欧美无线码| 每晚都被弄得嗷嗷叫到高潮| 菩萨蛮人人尽说江南好唐韦庄| 一个人免费看片子| 黑人猛操日本美女一级片| 亚洲国产精品999| 日日夜夜操网爽| 波野结衣二区三区在线| 男人爽女人下面视频在线观看| 又大又爽又粗| 天天躁夜夜躁狠狠躁躁| 久久久久国产一级毛片高清牌| 日韩精品免费视频一区二区三区| 国产日韩一区二区三区精品不卡| av有码第一页| 国产精品熟女久久久久浪| 麻豆av在线久日| 蜜桃国产av成人99| 超碰97精品在线观看| 亚洲精品日韩在线中文字幕| 三上悠亚av全集在线观看| 日韩一卡2卡3卡4卡2021年| 天天影视国产精品| 欧美大码av| 黄色一级大片看看| 高清欧美精品videossex| 亚洲国产欧美一区二区综合| 飞空精品影院首页| 一本色道久久久久久精品综合| 国产主播在线观看一区二区 | 精品久久久久久久毛片微露脸 | 老司机靠b影院| 色婷婷av一区二区三区视频| 一区二区日韩欧美中文字幕| 久久99热这里只频精品6学生| 久久久国产精品麻豆| 热99国产精品久久久久久7| 久久精品成人免费网站| 成年av动漫网址| av又黄又爽大尺度在线免费看| kizo精华| 亚洲国产欧美网| 大型av网站在线播放| 麻豆乱淫一区二区| 精品一区二区三区四区五区乱码 | 国产成人免费观看mmmm| 少妇被粗大的猛进出69影院| 丁香六月欧美| 高清黄色对白视频在线免费看| 亚洲中文字幕日韩| 色婷婷av一区二区三区视频| av视频免费观看在线观看| 啦啦啦中文免费视频观看日本| 精品国产乱码久久久久久男人| 久热爱精品视频在线9| 亚洲人成电影免费在线| 久久国产精品男人的天堂亚洲| 亚洲av电影在线观看一区二区三区| videosex国产| 水蜜桃什么品种好| 国产极品粉嫩免费观看在线| 青春草视频在线免费观看| 亚洲情色 制服丝袜| 色94色欧美一区二区| 中文乱码字字幕精品一区二区三区| 91国产中文字幕| 日韩伦理黄色片| 性色av一级| 两性夫妻黄色片| av在线app专区| 久久久国产一区二区| 久久狼人影院| 丰满人妻熟妇乱又伦精品不卡| 欧美日本中文国产一区发布| 中文字幕精品免费在线观看视频| 国产精品偷伦视频观看了| 亚洲国产精品国产精品| 精品一区在线观看国产| 亚洲av电影在线进入| 国产精品99久久99久久久不卡| 69精品国产乱码久久久| 精品一区二区三卡| 一区二区三区四区激情视频| 99国产精品一区二区蜜桃av | 国产精品人妻久久久影院| 亚洲精品一二三| av国产久精品久网站免费入址| 亚洲激情五月婷婷啪啪| 久久精品国产a三级三级三级| www.精华液| 五月开心婷婷网| 一边摸一边抽搐一进一出视频| 69精品国产乱码久久久| 精品视频人人做人人爽| 久久 成人 亚洲| 男的添女的下面高潮视频| 午夜视频精品福利| 亚洲国产中文字幕在线视频| 日韩熟女老妇一区二区性免费视频| 美女国产高潮福利片在线看| 中文字幕精品免费在线观看视频| 欧美精品一区二区免费开放| 搡老岳熟女国产| 久久久久久久精品精品| 在线看a的网站| 免费在线观看影片大全网站 | av在线app专区| 久久午夜综合久久蜜桃| 一级毛片黄色毛片免费观看视频| 久久久久久免费高清国产稀缺| 成人国产av品久久久| 狠狠精品人妻久久久久久综合| 2018国产大陆天天弄谢| 女性生殖器流出的白浆| 亚洲欧美一区二区三区国产| 日本wwww免费看| 91麻豆精品激情在线观看国产 | 欧美精品高潮呻吟av久久| 欧美日韩视频高清一区二区三区二| 国产伦人伦偷精品视频| 亚洲国产欧美日韩在线播放| 亚洲精品国产一区二区精华液| 国产一区二区激情短视频 | 亚洲成人国产一区在线观看 | 亚洲五月婷婷丁香| 日本91视频免费播放| 我的亚洲天堂| 永久免费av网站大全| 亚洲精品国产av蜜桃| svipshipincom国产片| 久久女婷五月综合色啪小说| 亚洲熟女毛片儿| 久久久久视频综合| 国产在线视频一区二区| 精品福利观看| 性少妇av在线| 中文字幕人妻丝袜一区二区| a级毛片黄视频| 午夜老司机福利片| 尾随美女入室| 国产老妇伦熟女老妇高清| 18禁黄网站禁片午夜丰满| 日韩欧美一区视频在线观看| 欧美日韩国产mv在线观看视频| 超色免费av| 欧美av亚洲av综合av国产av| 久久精品久久精品一区二区三区| 电影成人av| 新久久久久国产一级毛片| 我的亚洲天堂| 男的添女的下面高潮视频| 亚洲精品一区蜜桃| 精品一区二区三区av网在线观看 | 少妇被粗大的猛进出69影院| 美国免费a级毛片| 国产精品一国产av| 国产主播在线观看一区二区 | 国产黄色免费在线视频| 波多野结衣av一区二区av| 啦啦啦中文免费视频观看日本| 男的添女的下面高潮视频| 乱人伦中国视频| 亚洲情色 制服丝袜| 99国产精品一区二区蜜桃av | 欧美精品av麻豆av| 欧美日韩亚洲综合一区二区三区_| 老汉色∧v一级毛片| √禁漫天堂资源中文www| 丝袜美腿诱惑在线| 一级毛片电影观看| 成人亚洲欧美一区二区av| 亚洲中文字幕日韩| 高清av免费在线| 国产男女超爽视频在线观看| 999精品在线视频| svipshipincom国产片| 免费观看人在逋| 亚洲专区中文字幕在线| 国产日韩一区二区三区精品不卡| 亚洲一码二码三码区别大吗| 人人妻,人人澡人人爽秒播 | 国产又色又爽无遮挡免| 国产精品熟女久久久久浪| 日韩制服骚丝袜av| 亚洲,欧美,日韩|