• <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ù)的分類
    99久久人妻综合| 极品少妇高潮喷水抽搐| 高清欧美精品videossex| 欧美国产精品一级二级三级 | 69av精品久久久久久| 久久精品综合一区二区三区| 视频区图区小说| 成人毛片60女人毛片免费| h日本视频在线播放| 看黄色毛片网站| 亚洲人成网站在线播| 免费在线观看成人毛片| 午夜福利视频1000在线观看| 国产综合精华液| 亚洲av免费在线观看| 国产精品成人在线| 欧美三级亚洲精品| 少妇高潮的动态图| 18+在线观看网站| 国产成人精品一,二区| 蜜桃久久精品国产亚洲av| 麻豆精品久久久久久蜜桃| 午夜福利视频1000在线观看| 一二三四中文在线观看免费高清| 美女xxoo啪啪120秒动态图| 国产成人精品婷婷| 久久久成人免费电影| 一区二区三区精品91| 1000部很黄的大片| 人人妻人人澡人人爽人人夜夜| 欧美+日韩+精品| 另类亚洲欧美激情| 老女人水多毛片| 国产91av在线免费观看| 国产美女午夜福利| 天天躁日日操中文字幕| 少妇猛男粗大的猛烈进出视频 | 精品99又大又爽又粗少妇毛片| 51国产日韩欧美| 美女cb高潮喷水在线观看| av国产精品久久久久影院| 久久久久久国产a免费观看| 免费电影在线观看免费观看| 欧美亚洲 丝袜 人妻 在线| 国产女主播在线喷水免费视频网站| 精品国产一区二区三区久久久樱花 | 成人毛片60女人毛片免费| 成人亚洲欧美一区二区av| 十八禁网站网址无遮挡 | 一区二区av电影网| 国内揄拍国产精品人妻在线| 色网站视频免费| 亚洲aⅴ乱码一区二区在线播放| 亚洲成人一二三区av| 亚洲国产日韩一区二区| 视频区图区小说| 国产有黄有色有爽视频| 一区二区av电影网| 亚洲欧美成人综合另类久久久| 亚洲av福利一区| 亚洲精品国产色婷婷电影| 午夜免费鲁丝| 午夜福利视频1000在线观看| 在线免费观看不下载黄p国产| 亚洲成人一二三区av| 熟妇人妻不卡中文字幕| 乱码一卡2卡4卡精品| 成人欧美大片| 欧美人与善性xxx| av.在线天堂| av播播在线观看一区| 国产av不卡久久| 一区二区三区免费毛片| 欧美一级a爱片免费观看看| 寂寞人妻少妇视频99o| 亚洲一级一片aⅴ在线观看| 插阴视频在线观看视频| a级一级毛片免费在线观看| 亚洲精品自拍成人| 最近中文字幕高清免费大全6| 国产精品伦人一区二区| 国产男女超爽视频在线观看| 国产一区亚洲一区在线观看| 少妇人妻精品综合一区二区| 欧美激情国产日韩精品一区| www.色视频.com| 最后的刺客免费高清国语| 乱码一卡2卡4卡精品| 午夜福利在线观看免费完整高清在| av免费在线看不卡| 简卡轻食公司| 麻豆成人av视频| 久久99蜜桃精品久久| 美女xxoo啪啪120秒动态图| 色视频www国产| 色5月婷婷丁香| 2021天堂中文幕一二区在线观| 久久综合国产亚洲精品| 国产v大片淫在线免费观看| 人妻制服诱惑在线中文字幕| www.色视频.com| 亚洲综合精品二区| 男女那种视频在线观看| 男的添女的下面高潮视频| 国产一区二区亚洲精品在线观看| 麻豆乱淫一区二区| 久久鲁丝午夜福利片| 欧美bdsm另类| 日本一二三区视频观看| 成人午夜精彩视频在线观看| 国产爱豆传媒在线观看| 青青草视频在线视频观看| 建设人人有责人人尽责人人享有的 | 青春草国产在线视频| 久热这里只有精品99| 91久久精品电影网| 18禁在线播放成人免费| 国产精品久久久久久久久免| 伦理电影大哥的女人| 成人国产麻豆网| 久久久久九九精品影院| 一级毛片aaaaaa免费看小| 久久热精品热| 在线观看免费高清a一片| av黄色大香蕉| 男女啪啪激烈高潮av片| 亚洲国产精品成人综合色| 久久久久久久久久久丰满| 亚洲成色77777| 成人国产av品久久久| 人妻 亚洲 视频| 国产爱豆传媒在线观看| av一本久久久久| 亚洲精品第二区| 伦精品一区二区三区| 亚洲成人中文字幕在线播放| 禁无遮挡网站| 国产成人午夜福利电影在线观看| 少妇人妻 视频| 日韩免费高清中文字幕av| 成人一区二区视频在线观看| 哪个播放器可以免费观看大片| 久久韩国三级中文字幕| 99久久精品一区二区三区| 中文欧美无线码| 亚洲高清免费不卡视频| 成人一区二区视频在线观看| 又大又黄又爽视频免费| 日本欧美国产在线视频| 日本熟妇午夜| 又粗又硬又长又爽又黄的视频| 亚洲欧美一区二区三区国产| 免费av不卡在线播放| 国产欧美另类精品又又久久亚洲欧美| 日本黄色片子视频| 久久久久久久久久人人人人人人| 2018国产大陆天天弄谢| 亚洲av不卡在线观看| 美女脱内裤让男人舔精品视频| 男女边吃奶边做爰视频| 亚洲精品乱码久久久v下载方式| 国产精品久久久久久久电影| 午夜免费鲁丝| 美女被艹到高潮喷水动态| 少妇人妻一区二区三区视频| 日本欧美国产在线视频| 色5月婷婷丁香| 国产伦精品一区二区三区视频9| 2021少妇久久久久久久久久久| 91精品一卡2卡3卡4卡| 26uuu在线亚洲综合色| 久久久久九九精品影院| 日韩,欧美,国产一区二区三区| 99热网站在线观看| 婷婷色av中文字幕| av网站免费在线观看视频| 久久久精品94久久精品| 亚洲国产精品国产精品| 国产午夜福利久久久久久| 舔av片在线| 亚洲av在线观看美女高潮| av在线蜜桃| 超碰av人人做人人爽久久| 亚洲国产精品国产精品| 午夜福利在线在线| 亚洲国产最新在线播放| 搞女人的毛片| 麻豆国产97在线/欧美| 国产v大片淫在线免费观看| 夜夜看夜夜爽夜夜摸| 舔av片在线| 亚洲精品一区蜜桃| 亚洲精品久久久久久婷婷小说| 亚洲精品亚洲一区二区| 精品久久久精品久久久| 国产一区二区三区综合在线观看 | 国内精品宾馆在线| 中文资源天堂在线| 亚洲av福利一区| 亚洲精品456在线播放app| 全区人妻精品视频| 久久97久久精品| 夜夜爽夜夜爽视频| www.色视频.com| 99热全是精品| 欧美最新免费一区二区三区| 日韩强制内射视频| 日日啪夜夜爽| 亚洲成人精品中文字幕电影| 久久久久网色| 不卡视频在线观看欧美| 欧美少妇被猛烈插入视频| 久久久久久久国产电影| 欧美精品人与动牲交sv欧美| 国模一区二区三区四区视频| 不卡视频在线观看欧美| 久久ye,这里只有精品| 新久久久久国产一级毛片| 欧美亚洲 丝袜 人妻 在线| 18禁裸乳无遮挡免费网站照片| 在线播放无遮挡| 久久久久久久久大av| 久久99热6这里只有精品| 99视频精品全部免费 在线| 中国国产av一级| 国产 一区 欧美 日韩| 街头女战士在线观看网站| 久久久久九九精品影院| 十八禁网站网址无遮挡 | 国产精品久久久久久精品电影| 国产精品99久久99久久久不卡 | 五月天丁香电影| 日韩国内少妇激情av| 国产av不卡久久| 国产爽快片一区二区三区| 天天躁夜夜躁狠狠久久av| 嫩草影院入口| 亚洲欧美日韩东京热| 大片电影免费在线观看免费| 老师上课跳d突然被开到最大视频| av播播在线观看一区| 热99国产精品久久久久久7| 日韩,欧美,国产一区二区三区| 国产精品三级大全| 真实男女啪啪啪动态图| 久久99热这里只有精品18| 午夜激情福利司机影院| 大话2 男鬼变身卡| 观看美女的网站| 夜夜爽夜夜爽视频| 丰满人妻一区二区三区视频av| 性色av一级| 午夜福利在线观看免费完整高清在| 91在线精品国自产拍蜜月| 秋霞伦理黄片| 777米奇影视久久| 亚洲久久久久久中文字幕| 神马国产精品三级电影在线观看| 五月天丁香电影| 国产黄片美女视频| 女人十人毛片免费观看3o分钟| 少妇的逼水好多| 日日撸夜夜添| 成年av动漫网址| freevideosex欧美| 久久久久精品性色| 三级国产精品欧美在线观看| 春色校园在线视频观看| 黄色日韩在线| 亚洲国产精品成人综合色| 91精品一卡2卡3卡4卡| xxx大片免费视频| 九九爱精品视频在线观看| 日日啪夜夜撸| 97精品久久久久久久久久精品| 国产有黄有色有爽视频| 午夜福利视频精品| 大香蕉久久网| 嫩草影院新地址| 久热这里只有精品99| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 91久久精品电影网| 禁无遮挡网站| av.在线天堂| 一个人观看的视频www高清免费观看| 熟女av电影| 亚洲国产欧美在线一区| 国产高清三级在线| 亚洲精品乱久久久久久| 国产日韩欧美亚洲二区| 九九久久精品国产亚洲av麻豆| 超碰97精品在线观看| 纵有疾风起免费观看全集完整版| 国产精品嫩草影院av在线观看| av福利片在线观看| 嫩草影院入口| 成人亚洲精品一区在线观看 | 精品久久久精品久久久| 三级男女做爰猛烈吃奶摸视频| 最近中文字幕高清免费大全6| 18+在线观看网站| 国产国拍精品亚洲av在线观看| 日本三级黄在线观看| 亚洲天堂av无毛| 99热全是精品| 国产极品天堂在线| 最近的中文字幕免费完整| 永久网站在线| 最近手机中文字幕大全| 美女脱内裤让男人舔精品视频| 丝袜美腿在线中文| 岛国毛片在线播放| 亚洲色图综合在线观看| 九草在线视频观看| 一二三四中文在线观看免费高清| 少妇高潮的动态图| 久久韩国三级中文字幕| 成人亚洲精品一区在线观看 | 大香蕉97超碰在线| 人体艺术视频欧美日本| 能在线免费看毛片的网站| 白带黄色成豆腐渣| 嫩草影院入口| 在线a可以看的网站| 国精品久久久久久国模美| 亚洲国产精品999| 男女国产视频网站| 日韩在线高清观看一区二区三区| 在线观看国产h片| 人妻 亚洲 视频| 又粗又硬又长又爽又黄的视频| 国产亚洲av片在线观看秒播厂| 欧美bdsm另类| 内地一区二区视频在线| 最新中文字幕久久久久| 亚洲图色成人| 国产精品av视频在线免费观看| 亚洲av不卡在线观看| 少妇人妻一区二区三区视频| 又爽又黄无遮挡网站| 午夜福利在线在线| 亚洲精品第二区| 精品人妻偷拍中文字幕| 亚洲人与动物交配视频| 久久久午夜欧美精品| 欧美xxxx性猛交bbbb| 国内少妇人妻偷人精品xxx网站| 国产亚洲精品久久久com| av天堂中文字幕网| 一级毛片我不卡| 亚洲成人av在线免费| 亚洲aⅴ乱码一区二区在线播放| 看免费成人av毛片| 丝袜美腿在线中文| av福利片在线观看| 婷婷色综合www| 人体艺术视频欧美日本| 99久久精品国产国产毛片| 日本免费在线观看一区| 国产亚洲av嫩草精品影院| 亚洲电影在线观看av| 春色校园在线视频观看| 日本免费在线观看一区| 五月天丁香电影| 亚洲av日韩在线播放| 99久久精品国产国产毛片| 欧美日韩一区二区视频在线观看视频在线 | 女的被弄到高潮叫床怎么办| 黄片wwwwww| 三级国产精品片| 国产精品国产三级专区第一集| 一级黄片播放器| 亚洲av中文av极速乱| 97精品久久久久久久久久精品| 日韩av在线免费看完整版不卡| 欧美xxxx性猛交bbbb| 亚洲精品亚洲一区二区| 啦啦啦啦在线视频资源| 亚洲天堂av无毛| 久久久久久久大尺度免费视频| 欧美日韩视频精品一区| 女的被弄到高潮叫床怎么办| 亚洲美女视频黄频| 久久久a久久爽久久v久久| 搡女人真爽免费视频火全软件| 日本三级黄在线观看| 精品一区二区免费观看| 99热这里只有是精品50| 麻豆国产97在线/欧美| 久久久精品免费免费高清| 精品久久国产蜜桃| 免费看不卡的av| 久久久精品欧美日韩精品| 91精品一卡2卡3卡4卡| 国产精品女同一区二区软件| 日本熟妇午夜| 王馨瑶露胸无遮挡在线观看| 一级毛片黄色毛片免费观看视频| 麻豆成人午夜福利视频| 国产黄色免费在线视频| 久久精品久久久久久久性| 国产 一区精品| 香蕉精品网在线| 看非洲黑人一级黄片| 亚洲精品456在线播放app| 黄片wwwwww| 久久精品国产亚洲av天美| 美女国产视频在线观看| 青春草视频在线免费观看| 久久精品综合一区二区三区| 欧美+日韩+精品| 亚洲欧美日韩卡通动漫| 国产熟女欧美一区二区| 国产精品久久久久久精品电影小说 | 精品人妻视频免费看| 亚洲av不卡在线观看| 国产伦精品一区二区三区视频9| 乱码一卡2卡4卡精品| 亚洲国产欧美在线一区| 人人妻人人看人人澡| 成人漫画全彩无遮挡| 日本熟妇午夜| 国语对白做爰xxxⅹ性视频网站| 日本-黄色视频高清免费观看| 国产老妇伦熟女老妇高清| 日韩av不卡免费在线播放| 建设人人有责人人尽责人人享有的 | 国产女主播在线喷水免费视频网站| 国产精品99久久99久久久不卡 | av黄色大香蕉| 2018国产大陆天天弄谢| 国产伦理片在线播放av一区| 国产免费视频播放在线视频| 免费观看av网站的网址| 亚洲精品第二区| 亚洲欧美日韩无卡精品| 久久久久久久久久久丰满| 黑人高潮一二区| 久热这里只有精品99| 欧美一区二区亚洲| 伦理电影大哥的女人| 一本久久精品| 国产永久视频网站| 久热久热在线精品观看| 久久久久久伊人网av| 欧美97在线视频| 狂野欧美激情性xxxx在线观看| 麻豆成人av视频| 在线观看免费高清a一片| 夜夜看夜夜爽夜夜摸| 国产精品一区www在线观看| 亚洲精品,欧美精品| 亚洲精品国产av蜜桃| 色吧在线观看| 国产在线男女| 女人十人毛片免费观看3o分钟| 插逼视频在线观看| 久久精品夜色国产| 国产一区二区在线观看日韩| 老师上课跳d突然被开到最大视频| 久久综合国产亚洲精品| 久久鲁丝午夜福利片| 欧美成人一区二区免费高清观看| 在线看a的网站| 99热全是精品| 搡女人真爽免费视频火全软件| 又黄又爽又刺激的免费视频.| 久久精品综合一区二区三区| 免费大片黄手机在线观看| 一级毛片我不卡| 欧美精品人与动牲交sv欧美| 嘟嘟电影网在线观看| 久久久久久伊人网av| 成人欧美大片| av在线播放精品| 美女脱内裤让男人舔精品视频| 国产精品一区www在线观看| 黄片无遮挡物在线观看| av在线app专区| 又粗又硬又长又爽又黄的视频| 国产精品不卡视频一区二区| 六月丁香七月| av网站免费在线观看视频| 国产精品久久久久久精品古装| 九九久久精品国产亚洲av麻豆| 亚洲三级黄色毛片| 超碰av人人做人人爽久久| 夫妻性生交免费视频一级片| 亚洲成人久久爱视频| 亚洲精品日韩av片在线观看| av在线老鸭窝| 一级二级三级毛片免费看| 久久久久久久久久久免费av| 在线免费十八禁| 欧美日韩在线观看h| av黄色大香蕉| 国产视频首页在线观看| 国产白丝娇喘喷水9色精品| 国产老妇伦熟女老妇高清| 91aial.com中文字幕在线观看| 亚洲精品一区蜜桃| 亚洲精品国产成人久久av| 伦理电影大哥的女人| 久久久久网色| 五月伊人婷婷丁香| 精品一区二区三卡| 免费观看性生交大片5| 最近中文字幕2019免费版| 51国产日韩欧美| 免费播放大片免费观看视频在线观看| 日韩视频在线欧美| 亚洲av一区综合| 99视频精品全部免费 在线| 午夜激情福利司机影院| 精品久久久久久电影网| 国产精品国产av在线观看| 毛片一级片免费看久久久久| 天天躁夜夜躁狠狠久久av| 欧美极品一区二区三区四区| 国产男女内射视频| 九草在线视频观看| 久久久久精品性色| 欧美bdsm另类| av专区在线播放| 一二三四中文在线观看免费高清| 国产高潮美女av| 欧美日韩亚洲高清精品| 国产男女内射视频| 成年版毛片免费区| 国产人妻一区二区三区在| av免费观看日本| 日本一本二区三区精品| 日韩成人av中文字幕在线观看| 一级毛片黄色毛片免费观看视频| 色视频www国产| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产精品一及| 亚洲国产精品999| 久久久久九九精品影院| 国产午夜精品久久久久久一区二区三区| 精品亚洲乱码少妇综合久久| 嫩草影院精品99| 熟女电影av网| 又爽又黄无遮挡网站| 国产一区亚洲一区在线观看| 国产成人freesex在线| 亚洲精品日本国产第一区| eeuss影院久久| 成人欧美大片| 大片免费播放器 马上看| 亚洲精品乱码久久久久久按摩| 亚洲图色成人| 听说在线观看完整版免费高清| 国产有黄有色有爽视频| 免费看av在线观看网站| 日韩三级伦理在线观看| 国产一区二区三区综合在线观看 | 秋霞伦理黄片| 美女内射精品一级片tv| freevideosex欧美| 免费黄频网站在线观看国产| 亚洲av一区综合| 在线看a的网站| 欧美极品一区二区三区四区| 欧美另类一区| 老司机影院毛片| 久久99热6这里只有精品| 日韩三级伦理在线观看| 久久久久国产网址| 美女cb高潮喷水在线观看| 91aial.com中文字幕在线观看| 免费少妇av软件| 嫩草影院新地址| 午夜福利高清视频| 又爽又黄a免费视频| 一级二级三级毛片免费看| 午夜福利在线观看免费完整高清在| 国产v大片淫在线免费观看| 国产成人福利小说| 三级男女做爰猛烈吃奶摸视频| 国产亚洲一区二区精品| 又粗又硬又长又爽又黄的视频| 女人久久www免费人成看片| 午夜激情福利司机影院| 精品国产三级普通话版| 亚洲精品一二三| 99久国产av精品国产电影| 看非洲黑人一级黄片| 午夜免费男女啪啪视频观看| 婷婷色综合www| 五月伊人婷婷丁香| 亚洲人成网站在线观看播放| 亚洲va在线va天堂va国产| 国产真实伦视频高清在线观看| 午夜精品一区二区三区免费看| 久久韩国三级中文字幕| 只有这里有精品99| 精品一区二区三卡| 美女内射精品一级片tv| 日韩欧美精品v在线| 日韩一区二区视频免费看| 久久6这里有精品| 婷婷色麻豆天堂久久| 国产成人精品一,二区| 91久久精品国产一区二区成人| 2022亚洲国产成人精品| 免费少妇av软件| 欧美成人a在线观看| 国产午夜精品久久久久久一区二区三区| 联通29元200g的流量卡| 大片免费播放器 马上看| 欧美日韩在线观看h| 久热这里只有精品99| 美女cb高潮喷水在线观看| 久久精品久久久久久噜噜老黄|