A molecular phylogeny of Caraganeae (Leguminosae, Papilionoideae) reveals insights into new generic and infrageneric delimitations

Abstract Based on sequence data of nuclear ITS and plastid matK, trnL-F and psbA-trnH markers, the phylogeny of the subtribes Caraganinae and Chesneyinae in tribe Caraganeae was inferred. The results support the monophyly of each of the subtribes. Within subtribes Caraganinae, Calophaca and Halimodendron are herein transferred into Caragana to ensure its generic monophyly. The subtribe Chesneyinae is composed of four well-supported genera: Chesneya, Chesniella, Gueldenstaedtia and Tibetia. Based on phylogenetic, morphological, distributional and habitat type evidence, the genus Chesneya was divided into three monophyletic sections: Chesneya sect. Chesneya, Chesneya sect. Pulvinatae and Chesneya sect. Spinosae. Chesneya macrantha is herein transferred into Chesniella. Spongiocarpella is polyphyletic and its generic rank is not maintained. The position of Chesneya was incongruent in the nuclear ITS and the plastid trees. A paternal chloroplast capture event via introgression is hypothesized for the origin of Chesneya, which is postulated to have involved the common ancestor of Chesniella (♂) and that of the Gueldenstaedtia – Tibetia (GUT) clade (♀) as the parents.


Introduction
Caraganeae Ranjbar is a mid-sized tribe in Leguminosae, established by Ranjbar and Karamian (2003) based on fi ve genera: Calophaca Fisch. ex DC., Caragana Fabr., Chesneya Lindl. ex Endl., Gueldenstaedtia Fisch. and Halimodendron Fisch. ex DC., numbers of genera may be altered when treated by diff erent workers (see below). Caraganeae ranges from eastern Europe, central and western Asia to Mongolia, China and the Himalayas, extending northward to Siberia (Lock 2005;Ranjbar et al. 2014). Th is tribe is diagnosed by the asymmetrical axillary peduncles or pedicels attached to the slightly gibbous calyx and dehiscent pods (except for Halimodendron ;Polhill 1981;Ranjbar and Karamian 2003;Ranjbar et al. 2014).
A few recent studies referred to the concept of Caraganeae. Molecular work of Ranjbar et al. (2014) classifi ed Caraganeae into two subtribes: Caraganinae and Chesneyinae Ranjbar, F. Hajmoradi & Waycott. Duan et al. (2015) recognized this tribe based on the genera Calophaca, Caragana and Halimodendron. However, the former was inferred from a limited sampling scheme and few DNA markers, while the latter was subject to the undersampled for Chesneyinae. Hence, the monophyly of this tribe and the division of subtribes need to be further evaluated.
Th e genera Chesneya and Gueldenstaedtia formed a well-supported clade (Sanderson and Wojciechowski 1996), and were treated as the subtribe Chesneyinae (Ranjbar et al. 2014). Within this subtribe, the generic delimitations were controversial, especially concerning the status of Chesniella Boriss. (Borissova 1964), Spongiocarpella Yakovl. et Ulzij. (Yakovlev and Sviazeva 1987), and Tibetia (Ali) H. P. Tsui (Tsui 1979). Th e former two genera were separated from Chesneya, while Tibetia was a segregate of Gueldenstaedtia and has been revised in several studies (Cui 1998;Zhu 2004;Zhu 2005aZhu , 2005bBao and Brach 2010). Zhang et al. (2015b) supported the monophyly of Chesneya and proposed a classifi cation system, but some sections were only weakly supported. Hence, the phylogeny of Chesneyinae and its associated genera needs to be further explored.
We herein employ sequence data from nrDNA ITS and plastid matK, trnL-F and psbA-trnH to a) test the monophyly of Caraganeae and its subtribes; b) estimate the phylogeny of genera in Caraganeae; and c) discuss the taxonomic implications of this phylogeny on the generic and the infrageneric classifi cation of the tribe.
To better resolve the relationships of subtribes Caraganinae and Chesneyinae, 11 Galegeae species (8 genera) and 5 Hedysareae species (4 genera) were also sampled. Cicer microphyllum Royle ex Bentham, Dalbergia hupeana Hance, Lathyrus latifolius L., Robinia pseudoacacia L., Trifolium repens L. and Wisteria sinensis (Sims) Sweet were selected as outgroups based on previous studies (Wojciechowski et al. , 2004Wojciechowski 2003). Sequences of 40 accessions (representing 40 species) were downloaded from Gen-Bank (see Suppl. material 1 for details). Most accessions we sampled were collected from the fi eld or herbarium specimens. Onobrychis arenaria DC. was obtained from seedlings germinated from seeds provided by the Royal Botanic Gardens, Kew.

Phylogenetic analysis
Sequences were assembled with Geneious 7.1 (http://www.geneious.com/), and aligned using MUSCLE 3.8.31 (Edgar 2004), followed by manual adjustments in Geneious 7.1. Because the chloroplast markers putatively evolve as a single molecule, sequences of the three plastid markers (matK, trnL-F and psbA-trnH) were directly concatenated. Topological discordance was investigated by comparing the ITS and the concatenated plastid trees (as in García et al. 2014). To further determine the compatibility between these two datasets, an incongruence length diff erence (ILD) test and an approximately unbiased (AU) test were conducted with PAUP* (Swoff ord 2003) and CONSEL (Shimodaira and Hasegawa 2001; using site-wise likelihood values estimated by RA×ML; Stamatakis et al. 2008) programs, respectively. Th e tests retrieved the p values of 0.01 and 0.0001, respectively, suggesting that the incongruence between these two datasets was signifi cant. Th e ITS and the concatenated plastid sequences were thus analyzed separately.
Phylogenetic analyses were carried out using Bayesian inference (BI; Rannala and Yang 1996;Mau et al. 1999) with MrBayes 3.2.5 (Ronquist and Huelsenbeck 2003;Ronquist et al. 2012). Nucleotide substitution model parameters were determined prior to BI using the corrected Akaike information criterion (AIC) in jModeltest 2.1.7. (Posada 2008;Darriba et al. 2012). For the ITS dataset, boundaries of the 5.8S region to the ITS1 and the ITS2 regions were determined by comparison with the published 5.8S sequence of Vicia faba L. (Nazar and Wildeman 1981;Yokota et al. 1989), and the sequence substitution models for the ITS1, 5.8S and ITS2 regions were determined separately. Similarly, the models for each of the three plastid markers were estimated for the best-fi t models, which were used in the BI analysis for concatenated plastid sequences in a partitioned scheme.
In the BI, the Markov chain Monte Carlo (MCMC) search was run by two replicates for 10,000,000 generations, sampling one tree every 1,000 generations. After the fi rst 2,500,000 generations (2,500 trees) were discarded as burn-in, a 50% majorityrule consensus tree and posterior probabilities were obtained among the remaining trees. Results were checked using the program Tracer 1.5 (Rambaut and Drummond 2007) to ensure that plots of the two runs were converging and the value of the eff ective sample size for each replicate was above 200. Maximum likelihood (ML) analyses were conducted using RAxML-MPI v8.2 (Stamatakis 2014) with dataset partition scheme the same as in the BI and the following settings: rapid bootstrap analysis with 1,000 replicates and search for best-scoring ML tree in one program run, starting with a random seed, selecting the GTR model. Bootstrap values (LBS), as well as posterior probabilities (PP) were labeled on the corresponding branches of the Bayesian trees.

Results
Sequence characteristics are shown in Table 1. Our ML results are basically congruent in topology with the corresponding BI trees, the support values of the former were thus labeled on the corresponding branches of the latter (see legend of Figs 1, 2). Th anks to some extra sequences from GenBank (see Suppl. material 1), especially those of Chesneya and Chesniella, the ITS tree ( Fig. 1) was more comprehensively sampled than the plastid tree ( Fig. 2), which was of help to increase the general support of the former.

Nuclear data
In the ITS tree ( Fig. 1

Plastid data
Similar to the ITS results, the plastid tree ( Fig. 2) also showed the monophyly of both subtribes Caraganinae (PP = 1, LBS = 100%) and Chesneyinae (PP = 1, LBS = 100%). Calophaca and Halimodendron were nested in Caragana in diff erent places from the ITS tree, but such placement was weakly supported. Caragana also showed its paraphyly, Table 1. Sequence characteristics with gaps as missing data: alignment length, the number of the constant, variable and potential parsimony-informative (Pi) sites, and the best-fi t nucleotide substitution model determined by AIC.  Figure 1. Bayesian tree of the nrDNA ITS data, showing relationships of genera in subtribes Caraganinae, Chesneyinae and their close relatives. Th e labeled sections of Gueldenstaedtia and Tibetia followed Tsui (1979) and Zhu (2005a), respectively. Bayesian posterior probabilities (PP ≥ 0.95) and maximum likelihood bootstrap (LBS ≥ 70%) are given above and below branches, respectively. Th e asterisk indicates the name of Chesneya macrosperma has not been published, its voucher was storied in LE (details see Zhang et al. 2015b).   Unlike in the ITS tree, Chesneya s.s. and Chesniella were sisters in the plastid tree (PP = 1, LBS = 92%; Fig. 2). As in the ITS tree, the GUT clade (PP = 1, LBS = 100%) contained Gueldenstaedtia (PP = 1, LBS = 100%) and Tibetia (PP = 1, LBS = 100%), with each genus being monophyletic.

Discussion
Caraganeae comprises ca. 100 species distributed in temperate Asia, extending to eastern Europe (Ranjbar and Karamian 2003;Lock 2005). Th e two subtribes (Caraganinae and Chesneyinae) recognized by Ranjbar et al. (2014) are each well-supported in our analyses. However, our results did not recover a monophyletic Caraganeae (Figs 1,  2). Similarly, the previously expanded delimitation of Hedysareae sensu Lock (2005; also see Cardoso et al. 2013), which included the genera of subtribe Caraganinae and tribe Hedysareae sensu Amirahmadi et al. (2014), is not confi rmed herein (Figs 1, 2). Subtribe Caraganinae is composed of Calophaca, Caragana and Halimodendron (Ranjbar et al. 2014 Duan et al. 2015). Caraganinae is also distinct from Hedysareae (as delimited in Amirahmadi et al. 2014 andDuan et al. 2015) based on the following morphological characters: shrubs, rarely small trees; paripinnate, rarely imparipinnate leaves (Calophaca); solitary fl owers, or a few fl owers in fascicles, rarely forming a raceme; pods cylindric, rarely compressed, glabrous or hairy, with dehiscent and twisted valves (except for Halimodendron ;Polhill 1981;Liu et al. 2010b). Caraganinae is also related to the Astragalean clade; yet due to the morphological diversity of the latter, there are few diagnosable features to diff erentiate the Astragalean clade from Caraganinae, except for the twisted valves of Caraganinae (Calophaca and Caragana).
Calophaca morphologically resembles Caragana, and it is only distinguished from the latter by its imparipinnate leaves, rachises without thorns, and relatively denser racemes (Borissova 1945;Liu et al. 2010b). Calophaca contains 5-8 species mainly distributed in mountainous areas of central Asia, with one species extending to eastern Europe, and one endemic to northern China (Borissova 1945;Tutin et al. 1968;Yakovlev et al. 1996;Lock 2005;Liu et al. 2010b;Zhang et al. 2015a). Th e embedded position of Calophaca within Caragana argues that its classifi cation needs to be placed in the broader phylogenetic framework of Caragana, which is supported by our results (Figs 1, 2) and several previous studies (e.g., Zhang et al. 2009Zhang et al. , 2010Zhang et al. , 2015aDuan et al. 2015). We thus merge Calophaca into Caragana and recognize it at the sectional level as Caragana sect. Calophaca (Fisch. ex DC.) L.Duan, J.Wen & Zhao Y.Chang (see Taxonomic Treatment). Th e species-level nomenclatural changes will be made in a follow-up paper.
At the sectional level, our ITS tree (Fig. 1) indicated a strongly supported Car. sect. Calophaca. On the other hand, former Caragana s.s. was divided into fi ve sections mainly based on the combinations of leaf (pinnate or digitate) and petiole/rachis (persistent or caducous) characters (Zhang 1997). Th ree main sections, Car. sect. Bracteolatae, Car. sect. Caragana and Car. sect. Frutescentes, evolved likely accompanying the rapid uplifts of the Qinghai-Tibet Plateau (QTP) at around 8 Ma (Zhang et al. 2009). Th ese three sections also largely correspond to psychrophytic, mesophytic and xerophytic habitats, respectively (Zhang and Fritsch 2010). Our analyses supported the monophyly of the three sections, with Car. sect. Frutescentes only being monophyletic in the plastid tree (also see Zhang et al. 2009;Duan et al. 2015; and see below for an exceptional case in Car. sect. Frutescentes). Our ITS results failed to resolve a monophyletic Car. sect. Frutescentes (Fig. 1), but this may be due to insuffi cient informative sites in the ITS data. Furthermore, we only sampled one series for Car. sect. Spinosae (Car. ser. Spinosae Kom.), thus cannot assess its monophyly (Figs 1, 2). Caragana sect. Jabatae was suggested to have experienced a rapid radiation at 3.4-1.8 Ma (Zhang and Fritsch 2010), which may partly explain its poorly resolved relationships in our trees (Figs 1, 2; also see Zhang et al. 2009;Duan et al. 2015).
At the infra-sectional level, Car. ser. Bracteolatae Kom. and Car. ser. Spinosae are well-supported by our results (not labeled in the trees). Our results are therefore not completely congruent with Zhang et al. (2009), possibly due to diff erences in taxon sampling. Interestingly, a strongly supported psychric group is found within the mainly xeric section Car. sect. Frutescentes (Zhao 2009). Th is group is represented by Car. brevifolia Kom., Car. chinghaiensis Y.X. Liou, Car. densa Kom. and Car. versicolor Benth. (in Fig. 1; but weakly supported in the plastid tree). Most species of Car. sect. Frutescentes range from eastern Europe to northern China, Mongolia and Siberia, however, this abovementioned psychric group is distributed in the southern edge of northern China, extending to Tibet and its neighboring regions. It may represent a vicariant transitional group of Car. sect. Bracteolatae, Car. sect. Jubatae pro parte, Car. sect. Spinosae pro parte (psychrophytic habitat) and Car. sect. Frutescentes. Other cases of vicariant distributions have been noted in Caragana, and vicariance was considered as an important biogeographic pattern for this genus. For example, three closely related species in Car. sect. Caragana, Car. microphylla Lam., Car. intermedia Kuang & H.C.Fu and Car. korshinskii Kom., show non-overlapping to only slightly overlapping distributions in northeast to northwest China (Shue and Hao 1989;Zhang and Wang 1993;Zhang 1998;Chang 2008).
Th is subtribe contains ca. 50 species and diff ers from the Astragalean clade by twisted valves (e.g., in Chesneya), but a few species of Astragalus also have twisted legumes. Taxa of Chesneyinae are distinguished from Hedysareae by their dehiscent pods (Borissova 1945;Yakovlev et al. 1996;Liu et al. 2010a). Th e genera of Chesneyinae are distributed in central and eastern Asia, Tibet, Mongolia and Siberia, extending to eastern Turkey and Armenia ( Fig. 3A; Borissova 1945;Davis 1970;Rechinger 1984;Lock and Schrire 2005;Liu et al. 2010a), which are largely adapted to xerophytic (Chesneya and Chesniella), mesophytic (Gueldenstaedtia) and psychrophytic (Tibetia) habitats, respectively, although some species of Chesneya (see discussion below) and a few of Gueldenstaedtia are psychric taxa. Th e uplift of the QTP and aridifi cation of the former Tethys region might have driven the origination and divergence of genera in the subtribe Chesneyinae (Wen et al. 2014;Meng et al. 2015;Zhang et al. 2015b).

Topological discordance between ITS and plastid trees in subtribe Chesneyinae
Th e ITS and plastid topologies are incongruent within Chesneyinae. Chesneya s.s. formed a clade with the GUT clade in the ITS tree (Fig. 1), whereas it was sister to Chesniella in the plastid tree (Fig. 2). Both relationships were well-supported. Various mechanisms have been proposed to explain discordant topologies between gene trees, such as allopolyploidy, hybridization, horizontal gene transfer, incomplete lineage sorting (ILS), diff erent rate of molecular evolution, and chloroplast capture (Degtjareva et al. 2012;García et al. 2014;Yi et al. 2015).
Allopolyploidy can be ruled out for two reasons. First, taxa within Chesneyinae are diploid (Nie et al. 2002;Yang 2002;Sepet et al. 2014), with no evidence of polyploidy in this subtribe and its allied tribes. Second, deep lineages of Chesneyinae basically display a consistent chromosome number (x = 8; Nie et al. 2002;Sepet et al. 2014), with the only exception of Gueldenstaedtia (x = 7; Yang 2002), which has relatively recently diverged (ca. 15.23 Ma; Zhang et al. 2015b).
ILS and chloroplast capture seem more likely mechanisms for the present case (Tsitrone et al. 2003;Deng et al. 2015;Sun et al. 2015). A time-calibrated phylogeny may facilitate the exploration of the likely mechanism. Incomplete lineage sorting, which rarely occurs in deep lineage , prevails with bifurcation patterns of the shallow lineages of gene trees (especially at the specifi c level; Xu et al. 2012), and usually takes place in groups with relatively recent diversifi cation times (García et al. 2014). Zhang et al. (2015b) estimated that the main clades of subtribe Chesneyinae split at ca. 28 Ma, which is beyond the time frame supporting ILS of ancestral polymorphisms (as suggested by Xu et al. 2012). On the other hand, biogeographic patterns can also be taken into consideration (Goodman et al. 1999). Given peripatry and parapatry may have been involved in the evolution of Chesneyinae, if ILS occurred, the main clades would hardly be resolved with well-supported dichotomy as presented herein. Hence, although ILS could not be completely excluded in this case, we regarded chloroplast capture as the most likely cause for the discordant position of Chesneya s.s.
Compared to the biparental inheritance of the nuclear genome, plastid DNA of angiosperms is usually uniparentally transmitted, especially maternally (Corriveau and Coleman 1988;McCauley et al. 2007;Wicke et al. 2011). Nevertheless, the plastid DNA of the inverted repeat lacking clade (IRLC; see Figs 1, 2; also as in Lavin et al. 1990;Wojciechowski et al. 2000) in Leguminosae was reported to be inherited paternally or biparentally (Zhang et al. 2003), confi rmed by cytoplasmic and phylogenetic studies focusing on Medicago L. (paternal transmission; Schumann and Hancock 1989;Masoud et al. 1990;Havananda et al. 2010) and Wisteria Nutt. (Hu et al. 2005;Trusty et al. 2007). As Chesneya s.s. belongs to IRLC, a paternal inheritance scenario might be the case for the plastid DNA of Chesneya s.s.
We herein hypothesize a chloroplast capture event in the origin of Chesneya s.s. as follows. Th e common ancestor of Chesniella served as the putative paternal parent of Chesneya s.s. (sister to Chesneya s.s. in the plastid tree; Fig. 2). Th e maternal parent most likely was the common ancestor of the GUT clade. Th eir hybrids, with plastid from the paternal parent, may have continuously backcrossed with the maternal parent, and led to Chesneya s.s. inheriting most of the nuclear genome maternally (Fig. 1). Such a chloroplast capture event via introgression likely took place in the Miocene, because the divergence of Chesneya s.s. was dated to be 16.56 Ma and that of Chesniella was estimated as 19.81 Ma (Zhang et al. 2015b).
Analyses of Zhang et al. (2015b) revealed that the divergence of Chesneya and Chesniella most likely occurred around the QTP. Our analysis further indicated the psychric group of Chesneya diverged fi rst in this genus (C. sect. Pulvinatae, see Discussion below). It is probable that the common ancestor of Chesniella adapted to psychro-phytic habitats. However, the extant Chesniella is rarely distributed on the QTP. As for the GUT clade, Gueldenstaedtia possesses a unique chromosome number (x = 7; Yang 2002) within the subtribe. Most species of Gueldenstaedtia are adapted to mesophytic habitats of temperate northern and eastern Asia (Fig. 3A), in contrast to the rest of Chesneyinae, which are psychric or xeric taxa. Such a correlation among the variation of chromosome numbers and adaptation to diff erent habitats has also been recorded in other taxa, such as Hedysarum (Tang 2005;Duan et al. 2015), Passifl ora (Hansen et al. 2006) and Amaryllidaceae (García et al. 2014). But the mechanisms of these types of adaptation need to be further explored with robust phylogenetic, ecological and biogeographic analyses in our future eff orts.

Phylogeny and treatment of Chesneya, Chesniella and Spongiocarpella
Chesneya is the type genus of Chesneyinae, with ca. 35 species (see Fig. 3B-D). Th is genus has its distribution from the Himalayan region to northwestern China and Mongolia, through central and western Asia, westward to Turkey and Armenia (Fig. 3A;Borissova 1945;Davis 1970;Yakovlev et al. 1996;Lock and Schrire 2005;Fig. 3A). Our results suggest that the formerly circumscribed Chesneya, which contains two well-supported but separated parts: the core Chesneya s.s. and the outlier C. macrantha (Fig. 3E) (as in Li 1993 &Larsen 2010), is not monophyletic (Figs 1,  2). Chesneya spinosa P.C.Li (Fig. 3C) of Chesneya s.s. is morphologically similar to C. macrantha (Li 1981). However, C. spinosa is distributed in southern Tibet, while C. macrantha is restricted to the dry lands of Mongolia and northwestern China (Li and Ni 1985;Fu 1989). Th ey occupy psychrophytic and xerophytic habitats, respectively, and are clearly not sister to each other (Figs 1, 2).
Chesneya macrantha is nested within a monophyletic Chesniella according to our ITS tree (Fig. 1), and in the plastid tree, it is sister to the type of Chesniella: Ch. ferganensis (Korsh.) Boriss. (Borissova 1964;see Fig. 2, 3F). Chesneya macrantha shows some distinct morphologies from the other species in Chesniella, including its pulvinate habit and persistent leaf rachis (Li 1993), but this species generally share distribution areas, xerophytic habitats, and some synapomorphies, such as membranous stipules, hairy standard and ovate leafl ets with cuneate apices, with Chesniella (Li and Ni 1985;Fu 1989;Zhu and Larsen 2010). Th erefore, the transfer of Chesneya macrantha to Chesniella is supported by morphological, geographic and phylogenetic evidence (see Taxonomic Treatment). On the other hand, Chesneya was thus re-delimited based on the monophyletic Chesneya s.s.
After its establishment by Lindley (1839), Chesneya was divided into C. sect. Macrocarpon Boriss. and C. sect. Microcarpon Boriss. mainly based on pod morphology (Borissova 1945). Th e latter was segregated as the genus Chesniella by Borissova (1964), and this treatment was followed by Li (1993) and Zhu and Larsen (2010). Zhang et al. (2015b) informally classifi ed Chesneya into fi ve sections without detailed taxonomic treatment. Not all their sections were monophylytic, and the diagnostic characters and distributions of several sections were overlapping to some extent. Th e presently demarcated Chesneya was assigned into three strongly supported sections herein (as in the key of Chesneya proposed by Li 1993; details see Figs 1, 2 and Taxonomic Treatment). Chesneya sect. Macrocarpon possesses non-pulvinate habit, reduced stems, truncate or emarginate leafl et apices and caducous petiole and rachis (Borissova 1945). Th is section is composed of most species of Chesneya, including the type species: C. rytidosperma Jaub. et Spach (see Fig. 2; Borissova 1945;Davis 1970;Rechinger 1984). Chesneya sect. Macrocarpon was thus treated as C. sect. Chesneya (Fig. 3B). Unlike this section, petioles and rachises of C. sect. Pulvinatae M.L.Zhang (Zhang et al., 2015b; see Fig. 3D) are persistent and pubescent. However, most species in C. sect. Pulvinatae have blackened and curved petioles and rachises, while those of one of its species, C. spinosa, are hardened and spiny. Besides, C. spinosa formed a clade separated from C. sect. Pulvinatae. Hence, it is appropriate to segregate this species to form a new monotypic section: C. sect. Spinosae L.Duan, J.Wen & Zhao Y.Chang (see see Fig. 3C and Taxonomic treatment).
Th e infra-sectional relationships within C. sect. Chesneya are basically unresolved in our ITS trees (Fig. 1), and this section is undersampled in the plastid trees (Fig.  2). As for C. sect. Pulvinatae, two accessions of C. nubigena (D.Don) Ali formed a clade, being sister to C. purpurea P.C.Li (Figs 1, 2). Based on such well-supported tree topologies and several morphological diff erences, such as smaller leafl ets and purple corollae, the specifi c status of C. purpurea was retained herein (as in Li 1981Li , 1993. Th e xeric C. sect. Chesneya grows on dry slopes or desert margins of northwestern China, Mongolia and central Asia (see Fig. 3B; Borissova 1945;Rechinger 1984;Lock and Simpson 1991;Yakovlev et al. 1996;Zhu and Larsen 2010). Th is section is morphologically similar to Chesniella (Fig. 3F) and their distributions are more or less overlapping (Borissova 1945;Li, 1993), whereas they are not phylogenetically close to each other (Figs  1, 2). Such a phenomenon may be due to convergent evolution (Degtjareva et al. 2012). Chesneya sect. Spinosae (Fig. 3C) and C. sect. Pulvinatae (Fig. 3D) are restricted to Tibet and adjacent regions, adapting to high-altitude psychrophytic habitats (Ali 1977;Zhu and Larsen 2010). Th e evolutionary history of Chesneya appears complex, whereas the elevation of the QTP and the subsequent aridifi cations may have played an important role (Meng et al. 2015;Zhang et al. 2015b), as in former Calophaca (Zhang et al. 2015a), Caragana (Zhang and Fritsch 2010) and Hedysarum (Shue 1985;Duan et al. 2015).
Most previous workers did not accept the generic status of Chesniella, treating it within Chesneya (Borissova 1945;Li 1981;Rechinger 1984;Zhu and Cao 1986;Fu 1987Fu , 1989Yakovlev 1988;Yakovlev et al. 1991). Nevertheless, Li (1993) and Zhu and Larsen (2010) stated that the former is distinguishable from the latter by non-reduced stems, membranous stipules, obviously smaller calyxes, fl owers and pods. With the inclusion of Ch. macrantha (Fig. 3E), our results justifi ed the monophyly of Chesniella (Figs 1, 2), consistent with Zhang et al. (2015b). Within Chesniella, two well-supported groups were resolved in our ITS tree (Fig. 1). Chesniella macrantha and Ch. mongolica (Maxim.) Boriss. constituted group A, the group B included Ch. ferganensis, Ch. gracilis Boriss. and Ch. tribuloides (Nevski.) Boriss. Th e former confi ned in Mongolia and Inner Mongolia of China, to the contrast, the latter ranged from northwestern China to central Asia, which implied vicariance caused by Altai Moun-tain may drive the divergence of these two groups. However, due to undersampling and distinct morphology of Ch. macrantha in Chesniella, the evolution history and infrageneric taxonomy of this genus needs to be further explored. Yakovlev and Sviazeva (1987) erected Spongiocarpella as a segregate genus from Chesneya in the light of the former's spongiose legumes. Such treatment was followed by Yakovlev (1988), Fu (1989) and Yakovlev et al. (1996), but was rejected by Li (1993), Zhu (1996), Qian (1998) and Zhu and Larsen (2010). Based on fi eld and herbarium studies, we concur with Zhu (1996) that the sponge-like pericarp is an unstable character. Additionally, several species formerly assigned to Spongiocarpella were represented in our study, including Chesneya nubigena (D.Don) Ali, C. Spinosa and Chesniella macrantha. Th ey did not form a monophyletic group (Figs 1, 2). Th us, our data do not support the generic status of Spongiocarpella (as in Zhu 1996;Zhu and Larsen 2010;Ranjbar et al. 2014;Zhang et al. 2015b).

Monophyly of Gueldenstaedtia and Tibetia
Gueldenstaedtia is a small genus comprised of ca. 10 species and is distinguished from Chesneya by its palmately nerved wing petals (vs. pinnately in Chesneya) and non-twisted pod valves (vs. twisted) (see Fig. 3G; Liu et al. 2010a). Th is genus ranges from the Sino-Himalayan region to Mongolia and Siberia (Lock and Schrire 2005;see Fig. 3A). It was established by Fischer (1823) and revised by Fedtschenko (1927), Jacot (1927) and Kitagawa (1936). Ali (1962) divided it into G. subg. Gueldenstaedtia and G. subg. Tibetia Ali, but the latter was elevated to the generic rank by Tsui (1979) based on characters of stems, stipules, styles and seeds (see Fig. 3H). Th e genus Tibetia was generally accepted in subsequent revisions (Shue 1992;Yakovlev et al. 1996;Cui 1998;Wu 1999;Zhu 2004Zhu , 2005aBao and Brach 2010), and it is confi ned to Tibet and the adjacent regions including southern Gansu, southern Qinghai, western Sichuan and northwestern Yunnan of China, northern India, Nepal and Buhtan (Tsui 1979;Grierson and Long 1987;Lock and Schrire 2005;Zhu 2005a;Bao and Brach 2010).
Gueldenstaedtia and Tibetia were each supported to be monophyletic, and the two genera together form the GUT clade (Figs 1, 2). It seems valid to retain the generic status of each genus, which is also supported by karyological studies (Nie et al. 2002;Yang 2002;Zhu 2005b): Gueldenstaedtia (x = 7) vs. Tibetia (x = 8). Within Gueldenstaedtia, three species were sampled (all belonging to G. sect. Gueldenstaedtia according to Tsui 1979), but these species were all treated to be G. verna (Georgi) Boriss. s.l. by some workers (Yakovlev 1988;Zhu 2004;Bao and Brach 2010). Further work is needed to test the delimitation of G. verna s.l.
Within Tibetia, two accessions of T. himalaica (Baker) H.P.Tsui grouped together, which were sister to T. yadongensis H.P.Tsui (Figs 1, 2). Th e tree topology and the morphological characters (e.g., elongate stem and round or retuse leafl ets apex) seem to be consistent with treating T. himalaica as a distinct species (also see Tsui 1979;Cui 1998;Zhu 2005a;Bao and Brach 2010).  Description, distribution and habitat. Th is monotypic section is recognized by its hardened-spiny petiole and rachis. It is restricted in high-altitude psychrophytic rocky slope in southern Tibet.  Chang et al. 2004516 (WUK).

Key to the sections of Chesneya
Distribution and habitat. Dry slopes in Mongolia and Inner Mongolia, Ningxia and Xinjiang of China.