New insights into the phylogeny and taxonomy of Chinese Physospermopsis (Apiaceae)

Abstract Physospermopsis (Apiaceae) comprises about 10 species, but its taxonomy and phylogeny are disputed. The genus is mostly distributed in the Himalayas and Hengduan Mountains at high elevation. Earlier molecular studies involving six species of Physospermopsis indicated that this genus is not monophyletic and is nested in the East Asia Clade. Therefore, the aims of this study were to re-assess the phylogenetic position of, and interspecific relationships within, Physospermopsis based on two chloroplast loci (rpl16, rps16) and one nuclear region, the internal transcribed spacers of ribosomal DNA (ITS). Eight species involving 13 populations of Physospermopsis were collected. These were sequenced and analyzed with the sequences of 31 other Apiaceae species obtained from the NCBI to determine phylogenetic relationships using Bayesian inference (BI) and Maximum likelihood (ML). Our study found that Physospermopsis is monophyletic, nested in Pleurospermeae of Apiaceae, sister to Pleurospermum. And we propose that the Physospermopsis clade should be replaced by the East Asia Clade. However, the interspecific relationships within Physospermopsis were not well resolved and the positioning of species was unclear. Diagnostic characteristics to distinguish Physospermopsis species in the field and laboratory are provided for future Physospermopsis phylogenetic studies.


Introduction
Physospermopsis H. Wolff (1925: 276) has been reported to contain about 10 species, with eight species distributed in China (Pan and Watson 2005). However, 11 species were reported by Pimenov and Kljuykov (Pimenov and Kljuykov 2000a, b, c;Pimenov 2017). There has been difficulty in interpretation of diverse morphology to diagnose species and even limits of the genus. Most species of this genus occur in the Himalayas and Hengduan Mountains, and of these, four are endemic to the Hengduan Mountains (Wang and Pu 1992). In China, most species of Physospermopsis grow in open forests, scrubs, grasslands and alpine meadows at elevations of 2250-4800 m (Wolff 1929;Mukherjee 1982;Farille and Malla 1985;Pan and Watson 2005). Physospermopsis is characterized by having a long, conic taproot, ribbed stem, pinnate, rarely entire leaf blade, prominent, leaf-like bracts, variable bracteoles, minute calyx teeth, emerald young fruits, ovoid to broadly ovoid mature fruit with slightly cordate base (Wolff 1925;Pan and Watson 2005). Based on an analysis of previous research (Wang and Pu 1992;Liu 2005, 2006), Physospermopsis species usually possess a pericarp with wavy stria or reticulate ornamentation, prominent or inconspicuous fruit ribs, diverse carpoderms and endosperms, and pollen morphology showing a trend from rhomboidal type to rectangular type.
Previous studies on Physospermopsis have been extensive, including on micromorphology, anatomy and pollen morphology (Wang and Pu 1992;Liu 2005, 2006). However, previous molecular phylogenetic analyses of Physospermopsis have only involved a small number of taxa mostly using internal transcribed spacer (ITS) sequences (Downie et al. 2000;Calviño et al. 2006;Zhou et al. 2008Zhou et al. , 2009Downie et al. 2010;Valiejo-Roman et al. 2012). Phylogeny of Physospermopsis has been disputed with Downie et al. (2000) placing P. kingdon-wardii (H.Wolff) C. Norman (1938: 231) and P. rubrinervis (Franchet) C. Norman (1938: 231) in the Komarovia clade based on the materials collected from Yunnan, China. Then, Calviño et al. (2006) provisionally positioned Physospermopsis in the Physospermopsis clade based on more comprehensive maximum parsimony (MP) analyses of ITS sequences, which arose as a weakly supported sister group to the Komarovia clade. Later, Zhou et al. (2008) studied five species of Physospermopsis and concluded that Physospermopsis was not a monophyletic group. Zhou et al. (2008) placed P. kingdon-wardii and P. rubrinervis in the East Asia clade and referred them to Trachydium J. Lindley (1835: 232). Zhou et al. (2008) also concluded that P. cuneata H. Wolff (1929: 126) was nested in Pimpinelleae and should be close to Pimpinella C.Linnaeus (1753: 263), while P. muliensis R.H. Shan & S.L.Liou (1979: 105) and P. shaniana C.Y. Wu & F.T.Pu (1993: 1285 were allied within Pleurospermeae. Additionally, the East Asia clade was proposed as the synonym of Physospermopsis clade due to almost all of its species being primarily distributed in East Asia (Zhou et al. 2008). The following year, Zhou et al. (2009) added P. delavayi (the nomenclatural type of Physospermopsis) to their previous analyses and placed it in Pleurospermeae. Downie et al. (2010) decided that Physospermopsis should be placed in the Physospermopsis clade (East Asia clade) and was not monophyletic, but did not include the type species in their analyses potentially influencing their conclusions. Valiejo-Roman et al. (2012) conducted a molecular phylogenetic analysis of the genus Pleurospermum G.F.Hoffmann (Hoffmann 1814) and its related genera, including three Physospermopsis species.
Physospermopsis is a taxonomically complex genus whose generic limits with Pleurospermum, Tongoloa H. Wolff (1925: 279), and Trachydium are problematic (Pan and Watson 2005). Therefore, misidentification was common due to the absence of convincing morphological evidence, limitations of collected materials and examinations of type specimens. Additionally, until now there has been no comprehensive analysis using molecular phylogenetics and morphology within the one study. Therefore, we aimed to determine an accurate phylogeny of Physospermopsis and infrageneric relationships within Physospermopsis based on molecular, morphology data and combined analysis linking phylogeny and morphology. We acquired accurate data by collecting field specimens of eight Physospermopsis species involving 13 populations from their type localities and adjacent areas. Species were identified by field observations, validation with herbarium specimens and primary literature.
In the field investigations, we sampled three populations of P. delavayi, two populations of P. rubrinervis, two populations of P. shaniana, one population of P. obtusiuscula (1938: 231) and one population of P. nana (2000: 538) in Yunnan Province. We sampled one population of P. kingdon-wardii and one population of P. obtusiuscula in Tibet. One population of P. alepidioides (H. Wolff & Hand.-Mazz.) R.H. Shan (1941: 187) and one population of P. muliensis were sampled in Sichuan Province. All populations were collected from the type locality and adjacent regions, and the features of every species were closely matched with the types and original descriptions (de Candolle 1830;Franchet 1894;Wolff 1929;Wolff et Handel-Mazzetti 1933;Shan et Liou 1979). The specific collection information are listed in Appendix 1.
Fruits, leaf segments and specimens from these eight species of Physospermopsis were collected in the field for morphological study. Morphological analyses of leaves and fruits based on herbarium specimens or formaldehyde-acetic acid-alcohol (FAA) preserved material were photographed by a stereomicroscope Nikon SMZ25 (Japan). The morphological data were measured using KaryoType (Altnordu et al. 2016).

Taxon sampling
We sampled 13 populations, representing eight species of Physospermopsis in our phylogenetic analysis, and obtained 31 sequences of other Apiaceae species from the NCBI (Appendix 1). Based on previous research (Zhou et al. 2009), Bupleurum krylovianum B.K. Schischkin (1935Schischkin ( : 2010 and Bupleurum rockii H. Wolff (1929: 187) were selected as the outgroup for studying the phylogenetic position of Physospermopsis. We chose Pl. franchetianum W.B. Hemsley (1892: 307) and Pl. wrightianum H. Boissieu (1903: 847) as the outgroup for studying interspecific relationships within Physospermopsis. The DNA sequences of two chloroplast loci (rpl16, rps16) and one nuclear region, the internal transcribed spacers of ribosomal DNA (ITS), were used for phylogenetic analyses. According to the research to date (Zhou et al. 2008(Zhou et al. , 2009Downie et al. 2010;Guo et al. 2018;Panahi et al. 2018), these three markers should be sufficient to obtain the general information about relationships within the genus and its phylogenetic position within the family Apiaceae.

DNA extraction and sequencing
The fresh leaves were collected from field specimens in Yunnan, Sichuan and Tibet, China. Voucher specimens were deposited in the Herbarium of Sichuan University (SZ). Total genomic DNA was extracted from silica-dried leaves with plant genomic DNA kit (Cwbio Biosciences, Beijing, China). The universal primers ITS4 (5'-TCC TCCGCT TAT TGA TAT GC-3') and ITS5(5'-GGA AGT AAA AGT CGT AAC AAG G-3'; White et al. 1990) were used to amplify the entire internal transcribed sequences. The rpl16 intron region was amplified using primers F71(5'-GCT ATG CTT AGT GTG TGA CTC GTT G-3') and R1516 (5'-CCC TTC ATT CTT CTA TGT TG-3') (Jordan et al. 1996;Kelchner and Clark 1997). The rps16 sequences were amplified with primers rps16 3'exon (5'-CCT GTA GGY TGN GCN CCY TT-3') and rps16 5'exon (5'-AAA CGA TGT GGN AGN AAR CA-3') (Downie and Katz-Downie 1999). PCR amplification was implemented in a 30 μL volume reaction, including 3 μL total DNA, 1.5 μL forward primer, 1.5 μL reverse primer, 15 μL 2×Taq MasterMix (Cwbio, Beijing, China), and 9 μL ddH2O. The amplification of the ITS region was obtained by initial denaturation for 3 min at 94 °C, followed by 30 cycles of 45 s at 94 °C, 70 s at 54 °C, and 90 s at 72 °C, then final extension of 10 min at 72 °C. Amplification of cpDNA intron regions was obtained by initial denaturation for 3 min at 94 °C, followed by 36 cycles of 45 s at 94 °C,70 s at 58.5 °C, and 90 s at 72 °C, then final extension of 10 min at 72 °C. All PCR products were separated using a 1.5% (w/v) agarose TAE gel and sent to Sangon (Shanghai, China) for sequencing. New sequences obtained for this study have been deposited in GenBank. GenBank accession numbers are provided in the Appendix 1.

Data analysis
We used SegMan7 (Burland 2000) to assemble ITS and cpDNA sequences. ClustalX (Jeanmougin et al. 1998) was used to align DNA sequences with manual adjustment. We then used MEGA7 (Kumar et al. 2016) to manually adjust and obtain ITS and cpDNA datasets. Gaps were positioned to minimize nucleotide mismatches. Bayesian inference (BI) and Maximum likelihood (ML) methods were used for phylogenetic analyses, using MrBayes v3.2 (Ronquist et al. 2012) and RAxML v8.2.4 (Stamatakis 2014), respectively. Before undertaking BI analyses, MrModeltest version 2.2 (Nylander 2004) was used to determine the best model of nucleotide substitution and the GTR+G model under the Akaike Information Criterion (Akaike 1974) was selected. Bayesian analyses were performed over 20 million generations with a variant of Markov Chain Monte Carlo (MCMC) method and the trees were saved to a file every 1,000 generations. The first 20% trees were discarded as "burn-in" and the remaining 80% trees were used to build a majority-rule consensus tree based on analysis of the program Tracer v1.4 (Drummond and Rambaut 2007). ML analyses were performed using RAxML v8.2.4 with the GTR+G model and 1,000 bootstrap replicates. We constructed the BI tree with ITS data from all 44 taxa to test the systematic position of Physospermopsis. And we mapped some valuable morphological characteristics of Physospermopsis on phylogenetic tree, including leaves, bracts and bracteoles, ribs of fruits. The BI and ML trees were constructed for analysis of interspecific relationships within Physospermopsis using ITS and plastid datasets from the 13 Physospermopsis populations we sampled, one Physospermopsis species and the two Pleurospermum species downloaded from NCBI. Detailed information on the investigated taxa can be found in the Appendix 1.

Morphological characteristics of Physospermopsis
Through observations in the field, the most important characteristic to identify Physospermopsis species was prominent bracts and bracteoles. Physospermopsis shaniana, P. nana, P. muliensis, P. rubrinervis, P. obtusiuscula and P. kingdon-wardii usually have leaf-like bracts and bracteoles ( Fig. 1A3-F3). While P. alepidioides and P. delavayi possess lanceolate or oblong bracts and bracteoles with a 2-3-lobed apex and dark purple margin (Fig. 1G3, H3). Furthermore, leaf shape varies with species and can be obovate-lanceolate (e.g. P. alepidioides), triangular (e.g. P. rubrinervis), obovate-orbicular (e.g. P. delavayi) or linear-lanceolate (e.g. P. nana) segments ( Fig. 1B2, D2, G2, H2). Besides, the leaves of P. kingdon-wardii and P. obtusiuscula are 2-pinnate and ovateoblong, and have 2-6 pairs of ovate pinnae with pinnatisect margin. Physospermopsis muliensis and P. shaniana possess 3-5 pairs pinnae with pinnatifid margin, narrowly winged petioles and narrow and purple-red sheaths. Fruit morphology was recorded prior to alcohol preservation because the alcohol altered the color slightly (as is seen in photographs). The fruits of Physospermopsis were emerald green or chartreuse, ovoid to broadly ovoid, and typically had a slightly cordate base, a gradually narrowed and laterally flattened apex, with filiform or prominent ribs. Fruit shape and size of all Physospermopsis species were similar except that P. kingdon-wardii had fruit half the size of other species and very prominent and sinuate ribs. Physospermopsis nana and Morphological characters of Physospermopsis A1-H1 habit A2-H2 basal leaf A3-H3 umbel A4-H4 mericarps A1 habit of P. shaniana B1 habit of P. nana C1 habit of P. muliensis D1 habit of P. rubrinervis E1 habit of P. obtusiuscula F1 habit of P. kingdon-wardii G1 habit of P. delavayi H1 habit of P. alepidioides A2 basal leaf of P. shaniana B2 basal leaf of P. nana C2 basal leaf of P. muliensis D2 basal leaf of P. rubrinervis E2 basal leaf of P. obtusiuscula F2 basal leaf of P. kingdon-wardii G2 basal leaf of P. delavayi H2 basal leaf of P. alepidioides A3 umbel of P. shaniana B3 umbel of P. nana C3 umbel of P. muliensis D3 umbel of P. rubrinervis E3 umbel of P. obtusiuscula F3 umbel of P. kingdon-wardii G3 umbel of P. delavayi H3 umbel of P. alepidioides A4 mericarps of P. shaniana B4 mericarps of P. nana C4 mericarps of P. muliensis D4 mericarps of P. rubrinervis E4 mericarps of P. obtusiuscula F4 mericarps of P. kingdon-wardii G4 mericarps of P. delavayi H4 mericarps of P. alepidioides. P. muliensis fruits had relatively prominent and filiform ribs, but P. muliensis fruits had scattered warts especially on the ribs and P. nana had smaller fruit. Physospermopsis obtusiuscula fruits were ovoid with narrowly winged and sinuolate ribs. The fruit of P. delavayi had an obvious cordate base, and filiform and less prominent ribs. Physospermopsis alepidioides, P. rubrinervis and P. shaniana had ovoid, verucose fruits with prominent ribs, but P. alepidioides did not have a cordate base, while the other two species had a slightly cordate base. The fruit of P. shaniana had many small warts distinguishing it from P. rubrinervis. For easy reading and comparison, the main morphological characteristics were listed in Table 1.
The ITS dataset tree topologies generated from BI and ML analyses were consistent. Therefore, only the BI tree with posterior probabilities (PP, 0-1) and bootstrap support values (BS, 0-100%) is illustrated in Fig. 3A. The first to differentiate from Physospermopsis was P. muktinathensis, which is distributed in Nepal. Three populations of P. delavayi and one of P. alepidioides united as a strongly supported (BI-PP = 1; ML-BS = 100%) group. Physospermopsis obtusiuscula was supported as a sister group to P. kingdon-wardii (BI-PP = 1; ML-BS = 100%). Physospermopsis rubrinervis, P. muliensis, P. nana and P. shaniana were allied in all trees (BI-PP = 1; ML-BS = 100%). However, clear interspecific relationships between P. rubrinervis, P. muliensis, P. nana and P. shaniana were not strongly supported by ML or BI analyses.  The cpDNA dataset tree topologies inferred by BI and ML analyses were consistent (Fig. 3B). However, results of a partition homogeneity test for the ITS and cpDNA datasets indicated that these genomes provide significantly different phylogenetic estimates. The taxa involved in this conflict are highlighted in Fig. 3. There was no chloroplast data for P. muktinathensis. The first to differentiate were P. obtusiuscula and P. kingdon-wardii (BI-PP = 1; ML-BS = 99%). The relationships of the three P. delavayi populations differed from the ITS dataset tree topology, although this cpDNA dataset relationship was not strongly supported. The cpDNA dataset tree topologies indicated that LYG population was closer to the HB population (BI-PP = 0.33; ML-BS = 47%), while LYG was closer to LGH in the ITS dataset tree topologies (BI-PP = 0.34; ML-BS = 70%). Additionally, the relationships between P. rubrinervis, P. muliensis, P. nana were not consistent with the ITS tree, where P. nana allied with P. rubrinervis in the cpDNA tree (BI-PP = 1; ML-BS = 72%), whereas P. nana allied with P. muliensis in the ITS tree (BI-PP = 1; ML-BS = 94%).

The Phylogenetic position of Physospermopsis and relationship between Physospermopsis and Pleurospermum
Physospermopsis is monophyletic. The reasons for previous designations as a polyphyletic genus were likely attributable to the misidentification of several species (e.g. P. rubrinervis, P. kingdon-wardii, P. cuneata). Besides, P. cuneata is a poorly known species and unusual within the genus for its lack of conspicuous bracts and bracteoles, and therefore the phylogenetic placement of it is highly controversial. However, the most recent consensus is that P. cuneata should not be placed in Physospermopsis (Zhou et al. 2008;Zhou et al. 2009;Pimenov 2017). So previous molecular studies only involved five physospermopsis species which were widely accepted; we added another three physospermopsis species in this study, including P. alepidioides, P. obtusiuscula, and P. nana. Evidence obtained through more precise checking of generic type, infrageneric types and extensive herbarium specimens, literature and field investigations, analyzing morphological characters, and ITS and cpDNA evidence. This comprehensiveness allows us to be con- fident that Physospermopsis is monophyletic and nested in Pleurospermeae. In addition, we propose that the Physospermopsis clade should be replaced by the East Asia Clade.
The molecular results indicated that Physospermopsis is closest to Pleurospermum. Morphologically, Pleurospermum usually possess numerous bracts and bracteoles with white scarious margins, conspicuous or obsolete calyx teeth, white or purple-red petals with clawed base and narrow apex, prominent, acute ribs (Pan and Watson 2005). However, we found that Physospermopsis differed from Pleurospermum by less prominent and even inconspicuous fruit ribs, and the bracts and bracteoles did not have white scarious margins, resulting in an obvious, diagnostic boundary between Pleurospermum and Physospermopsis. The closeness of the two genera is also evidenced in pollen morphology. Wang and Pu (1992) found P. alepidioides and P. muliensis pollen to be rhomboidal and similar to several Pleurospermum species whereas other Physospermopsis species (P. rubrinervis and P. delavayi) have more advanced rectangular types. In addition, Pan and Watson (2005) identified several Physospermopsis species (e.g. P. obtusiuscula) with morphological similarities to Pleurospermum species, including having long fruit ribs and bulgy fruit walls, while other species had flattening of fruit and reduced wall thickness. Consequently, Physospermopsis is sister to Pleurospermum.

Interspecific relationships within Physospermopsis
The morphological characteristics mapped on the phylogenetic tree indicated that most closely related species have similar morphological characteristics. For example, P. rubrinervis, P. muliensis and P. shaniana are highly consistent on leaves, bracts and bracteoles, ribs on fruits (Fig. 2). Similarly, these species are the geographically sympatric species (Fig. 4). Resolution of the relationships between these species will only be achieved through continued studies, which may be difficult due to their geographic and morphological similarities. However, we can learn that P. nana, P. rubrinervis, P. muliensis and P. shaniana are the more advanced species in Physospermopsis. The morphological characters of P. nana are the most particular; these might be caused by hybridization with Pleurospermum species.
The interspecific relationships between certain species within Physospermopsis are evident based on the consistencies between ITS and cpDNA trees. For instance, P. alepidioides showed a close affinity to P. delavayi in phylogenetic tree, and they have similar bracts and bracteoles (entire or 2-3-lobed at apex, with dark purple margin) (Figs 1, 2). However, differing leaf shapes can be used to easily distinguish these two species because P. alepidioides has an undivided leaf with sparsely serrated margin (Fig. 1H2) and P. delavayi has a pinnate leaf (Fig. 1G2). Physospermopsis kingdon-wardii is sister to P. obtusiuscula, which is consistent with their geographic closeness. Physospermopsis kingdon-wardii appears more morphologically similar to P. obtusiuscula (including the leaves, bracts and bracteoles), but differs in its reduced stem, small stature and small fruits with prominent and sinuate ribs (Fig. 1).
The topologies of the ITS and cpDNA trees differed in the positioning of P. delavayi, P. muliensis and P. nana (Fig. 3). This inconsistency between nrDNA ITS and cpDNA data has been reported in some studies of Apiaceae (Lee and Downie 2006;Zhou et al. 2008Zhou et al. , 2009Spalik et al. 2009;Bone et al. 2011;Yi et al. 2015;Panahi et al. 2018). This difference generally has been caused by incomplete lineage sorting, hybridization, homoplastic substitutions and introgression. Since we did not sample by lineage and execute gene flow analysis, what caused the inconsistency cannot be determined. Previous studies have indicated that Pleurospermeae occupies a relative position in the base of the Apioideae (Zhou et al. 2008(Zhou et al. , 2009Downie et al. 2010), the differentiation time should be earlier. Thus, for Physospermopsis, we infer the more effective reason for the inconsistency between nrDNA ITS and cpDNA data is hybridization. A further study based on widely sampling and deeper analysisis needed. However, several diagnostic characteristics can be utilized in the field and laboratory to separate them. P. rubrinervis can easily be recognized by purple-red nerves on the leaves, bracts and bracteoles with purple-red margin (Fig. 1). Physospermopsis muliensis possesses a slender, branched stem and narrowly winged basal petioles with narrow sheaths (Fig. 1). Physospermopsis nana has bracts and bracteoles with white scarious margins and linear-lanceolate segments with membranous-margined sheaths (Fig. 1). The stem of P. shaniana was reduced and branched at the base, and had prominent bracts 1-2-pinnate with developed, broad sheaths (Fig. 1).  Diagnostic characters. Physospermopsis alepidioides usually possesses an entire leaf blade with a sparsely serrated margin. The shape of the entire leaf segment is an obvious diagnostic characteristic to distinguish it from other Physospermopsis species. The stem of it is velutinous.

Chinese Physospermopsis species
Distribution. Endemic to China, Sichuan (Fig. 4) Diagnostic characters. Physospermopsis rubrinervis usually possesses dark purple, sparsely branching stems. The basal blade is ovate to broadly ovate in outline, having almost purple-red nerves.
Habitat. This species grows in the forest edge or rhododendron shrubs at an elevation of 2800-4800 m. Additional
Diagnostic characters. P. kingdon-wardii is similar to P. obtusiuscula in shape of basal leaves, but the stem of P. kingdon-wardii is reduced. The fruits are smaller than other species, and the immature fruits sometimes have sparse minute warts. Additionally, the ribs are prominent, often sinuate.

Conclusion
Physospermopsis is monophyletic and nested in Pleurospermeae, sister to Pleurospermum. Although the interspecific relationships within Physospermopsis were not well resolved and the positioning of species was unclear, the relationships of P. alepidioides and P. delavayi, P. kingdon-wardii and P. obtusiuscula are close. Diagnostic characteristics for distinguishing the species in the field and laboratory are provided for necessary morphological and molecular research in future Physospermopsis phylogenetic studies.