Systematics of Mukdenia and Oresitrophe (Saxifragaceae): Insights from genome skimming data

Oresitrophe and Mukdenia (Saxifragaceae) are epilithic sister genera used in traditional Chinese medicine. The taxonomy of Mukdenia, especially of M. acanthifolia, has been controversial. To address this, we produced plastid and mitochondrial data using genome skimming for Mukdenia acanthifolia and Mukdenia rossii, including three individuals of each species. We assembled complete plastomes, mitochondrial CDS and nuclear ribosomal ETS/ITS sequences using these data. Comparative analysis shows that the plastomes of Mukdenia and Oresitrophe are relatively conservative in terms of genome size, structure, gene content, RNA editing sites and codon usage. Five plastid regions that represent hotspots of change (trnH-psbA, psbC-trnS, trnM-atpE, petA-psbJ and ccsA-ndhD) are identified within Mukdenia, and six regions (trnH-psbA, petN-psbM, trnM-atpE, rps16-trnQ, ycf1 and ndhF) contain a higher number of species-specific parsimony-informative sites that may serve as potential DNA barcodes for species identification. To infer phylogenetic relationships between Mukdenia and Oresitrophe, we combined our data with published data based on three different datasets. The monophyly of each species (Oresitrophe rupifraga, M. acanthifolia and M. rossii) and the inferred topology ((M. rossii, M. acanthifolia), O. rupifraga) are well supported in trees reconstructed using the complete plastome sequences, but M. acanthifolia and M. rossii did not form a separate clade in the trees based on ETS + ITS data, while the mitochondrial CDS trees are not well-resolved. We found low recovery of genes in the Angiosperms353 target enrichment panel from our unenriched genome skimming data. Hybridization or incomplete lineage sorting may be the cause of discordance between trees reconstructed from organellar and nuclear data. Considering its morphological distinctiveness and our molecular phylogenetic results, we strongly recommend that M. acanthifolia be treated as a distinct species.     

Phylogenomic conflict analyses in the apple genus Malus s.l. reveal widespread hybridization and allopolyploidy driving diversification,with insights into the complex biogeographic history in the Northern Hemisphere

Phylogenomic evidence from an increasing number of studies has demonstrated that different data sets and analytical approaches often reconstruct strongly supported but conflicting relationships. In this study, 785 single-copy nuclear genes and 75 complete plastomes were used to infer the phylogenetic relationships and estimate the historical biogeography of the apple genus Malus sensu lato, an economically important lineage disjunctly distributed in the Northern Hemisphere and involved in known and suspected hybridization and allopolyploidy events. The nuclear phylogeny recovered the monophyly of Malus s.l. (including Docynia); however, the genus was supported to be biphyletic in the plastid phylogeny. An ancient chloroplast capture event in the Eocene in western North America best explains the cytonuclear discordance. Our conflict analysis demonstrated that ILS, hybridization, and allopolyploidy could explain the widespread nuclear gene tree discordance. One deep hybridization event (Malus doumeri) and one recent event (Malus coronaria) were detected in Malus s.l. Furthermore, our historical biogeographic analysis integrating living and fossil data supported a widespread East Asian-western North American origin of Malus s.l. in the Eocene, followed by several extinction and dispersal events in the Northern Hemisphere. We also propose a general workflow for assessing phylogenomic discordance and biogeographic analysis using deep genome skimming data sets.     

基于叶绿体rps16基因和核基因ITS片段研究肉苁蓉属系统位置

在前人对列当科系统发育研究的基础上,追加了肉苁蓉属(Cistanche)的基因序列数据,运用最大简约法、最大似然法和贝叶斯推断方法探讨了其在列当科中的系统位置及列当科中属间关系。基于rps16基因序列及rps16+ITS联合序列建立了列当科系统发育树,结果显示,肉苁蓉属、列当属(Orobanche)以及草苁蓉属(Boschniakia)聚在同一进化枝内,肉苁蓉属和列当属表现出最近亲缘关系;列当科中的全寄生类群、半寄生类群和非寄生类群分属在3个不同分支中。     

Resolving robust phylogenetic relationships of core Brassicaceae using genome skimming data

The Brassicaceae is an economically and scientifically important family distributed globally, including oilseed rape and the model plant, Arabidopsis thaliana. Although growing molecular data have been used in phylogenetic studies, the relationships among major clades and tribes of Brassicaceae are still controversial. Here, we investigated the core Brassicaceae phylogenetics using 222 plastomes and 235 nrDNA cistrons, including 106 plastomes and 112 nrDNA cistrons assembled from newly sequenced genome skimming data of 112 taxa. The sampling covered 73 genera from 61.5% tribes and four unassigned genera and species. Three well supported lineages LI, LII, and LIII were revealed in our plastomic analyses, with LI sister to LII + LIII. In addition, the monophyly of the newly delimitated LII was strongly supported by three different partition strategies, concatenated methods under Bayesian and Maximum Likelihood analyses. LII comprised 13 tribes, including four tribes previously unassigned to any lineage, that is Biscutelleae as the earliest diverging clade and Cochlearieae as the sister to Megacarpaeeae + Anastaticeae. Within LII, the intertribal relationships were also well resolved, except that a conflicting position of Orychophragmus was detected among different datasets. In LIII, Shehbazia was resolved as a member of Chorisproreae, but Chorisproreae, Dontostemoneae, and Euclidieae were all resolved as paraphyletic, which was also confirmed by nrDNA analyses. Moreover, the loss of the rps16 gene was detected as likely to be a synapomorphy of the tribes Arabideae and Alysseae. Overall, using genome skimming data, we resolved robust phylogenetic relationships of core Brassicaceae and shed new light on the complex evolutionary history of this family.     

Recent advances on phylogenomics of gymnosperms and a new classification

Living gymnosperms comprise four major groups: cycads, Ginkgo, conifers, and gnetophytes. Relationships among/within these lineages have not been fully resolved. Next generation sequencing has made available a large number of sequences, including both plastomes and single-copy nuclear genes, for reconstruction of solid phylogenetic trees. Recent advances in gymnosperm phylogenomic studies have updated our knowledge of gymnosperm systematics. Here, we review major advances of gymnosperm phylogeny over the past 10 years and propose an updated classification of extant gymnosperms. This new classification includes three classes (Cycadopsida, Ginkgoopsida, and Pinopsida), five subclasses (Cycadidae, Ginkgoidae, Cupressidae, Pinidae, and Gnetidae), eight orders (Cycadales, Ginkgoales, Araucariales, Cupressales, Pinales, Ephedrales, Gnetales, and Welwitschiales), 13 families, and 86 genera. We also described six new tribes including Acmopyleae Y. Yang, Austrocedreae Y. Yang, Chamaecyparideae Y. Yang, Microcachrydeae Y. Yang, Papuacedreae Y. Yang, and Prumnopityeae Y. Yang, and made 27 new combinations in the genus Sabina.     

Nuclear and plastid phylogenomic analyses provide insights into the reticulate evolution,species delimitation,and biogeography of the Sino-Japanese disjunctive Diabelia (Caprifoliaceae)

Understanding biological diversity and the mechanisms of the Sino-Japanese disjunctions are major challenges in eastern Asia biogeography. The Sino-Japanese flora has been broadly studied as an ideal model for plant phylogeography. Diabelia Landrein (Caprifoliaceae) is an East Asian genus, with a disjunctive distribution across the Sino-Japanese region. However, relationships within Diabelia remain elusive. In this study, we reconstructed the phylogeny of Diabelia and inferred historical biogeography and evolutionary patterns based on nuclear and plastid sequences from target enrichment and genome skimming approaches, respectively. We found that the main clades within Diabelia were discordant between nuclear and plastid trees. Both nuclear and plastid phylogenetic analyses supported five main clades: Diabelia serrata (Siebold & Zucc.) Landrein, Diabelia tetrasepala (Koidz.) Landrein, Diabelia sanguinea (Makino) Landrein, Diabelia stenophylla (Honda) Landrein, and Diabelia spathulata (Siebold & Zucc.) Landrein. Species network analyses revealed that Diabelia tetrasepala is likely the result of a hybridization event. Divergence time estimation and ancestral area reconstructions showed that Diabelia originated in Japan during the early Miocene, with subsequent vicariance and dispersal events between Japan and Korea, and between Japan and China. Overall, our results support the division of Diabelia into five main clades and the recognition of five species in the genus. This research provides new insights into the species delimitation and speciation processes of taxonomically complex lineages such as Diabelia.     

A decade of progress in plant molecular phylogenetics

Over the past decade, botanists have produced several thousand phylogenetic analyses based on molecular data, with particular emphasis on sequencing rbcL, the plastid gene encoding the large subunit of Rubisco (ribulose bisphosphate carboxylase). Because phylogenetic trees retrieved from the three plant genomes (plastid, nuclear and mitochondrial) have been highly congruent, the ‘Angiosperm Phylogeny Group’ has used these DNA-based phylogenetic trees to reclassify all families of flowering plants. However, in addition to taxonomy, these major phylogenetic efforts have also helped to define strategies to reconstruct the ‘tree of life’, and have revealed the size of the ancestral plant genome, uncovered potential candidates for the ancestral flower, identified molecular living fossils, and linked the rate of neutral substitutions with species diversity. With an increased interest in DNA sequencing programmes in non-model organisms, the next decade will hopefully see these phylogenetic findings integrated into new genetic syntheses, from genomes to taxa.     

Out of the Middle East: New phylogenetic insights in the genus Tamarix (Tamaricaceae)

Tamarix is one of the taxonomically most complex genera among the angiosperms, and there is little consensus regarding its infrageneric classification. Here we present the most complete phylogenetic reconstruction of the genus to date. This includes a DNA phylogenetic tree based on nuclear ribosomal ITS, and a plastid DNA phylogeny based on three intergenic spacers (trnS‐trnG, ndhF‐rpl32, and trnQ‐rps16). In total, both nuclear and plastid phylogenetic analyses include more than 70 samples of 39 species from 27 countries, which represent close to 60% of the diversity of the genus. Two complementary trees, based only on one plastid marker, are also included. The first, based on trnS‐trnG, is used to increase the number of species related to T. amplexicaulis. The second, based on ndhF‐rpl32, is used to investigate the separation between T. tetrandra and T. parviflora. The incongruence between the available infrageneric classifications and the molecular results is confirmed. A reticulate evolution is inferred from the trees, showing characters such as vaginate leaves appearing at different stages along the evolutionary history of the genus. The presence of T. canariensis outside the Canary Islands is cast into doubt, and all such records from NW Africa and Europe are here considered to belong to T. gallica. The results also suggest independence of T. karelinii from T. hispida, and T. parviflora from T. tetrandra. Relationships between a number of species are still not resolved, and additional studies will be needed to further refine the complex taxonomy of Tamarix.     

Plastome phylogenomics and phylogenetic diversity of endangered and threatened grassland species (Poaceae) in a North American tallgrass prairie

Native grasslands are one of the most endangered ecosystems in North America. In this study, we examined the ecological and evolutionary roles of endangered and threatened (e/t) grasses by establishing robust evolutionary relationships with other nonthreatened native and introduced grass species of the community. We hypothesized that the phylogenomic distribution of e/t species of grasses in Illinois would be phylogenetically clustered because closely related species would be vulnerable to the same threats and have similar requirements for survival. This study presents the first time a phylogeny based on complete plastome DNA of Poaceae was analyzed by phylogenetic diversity analysis. To avoid the disturbance of e/t populations, DNA was extracted from herbarium specimens. Next‐generation sequencing (NGS) techniques were used to sequence DNA of plastid genomes (plastomes). The resulting phylogenomic tree was analyzed by phylogenetic diversity metrics. The extracted DNA successfully produced complete plastomes demonstrating that herbarium material is a practical source of DNA for genomic studies. The phylogenomic tree was strongly supported and defined Dichanthelium as a separate clade from Panicum. The phylogenetic metrics revealed phylogenetic clustering of e/t species, confirming our hypothesis.     

Mechanisms of Functional and Physical Genome Reduction in Photosynthetic and Nonphotosynthetic Parasitic Plants of the Broomrape Family

Nonphotosynthetic plants possess strongly reconfigured plastomes attributable to convergent losses of photosynthesis and housekeeping genes, making them excellent systems for studying genome evolution under relaxed selective pressures. We report the complete plastomes of 10 photosynthetic and nonphotosynthetic parasites plus their nonparasitic sister from the broomrape family (Orobanchaceae). By reconstructing the history of gene losses and genome reconfigurations, we find that the establishment of obligate parasitism triggers the relaxation of selective constraints. Partly because of independent losses of one inverted repeat region, Orobanchaceae plastomes vary 3.5-fold in size, with 45 kb in American squawroot (Conopholis americana) representing the smallest plastome reported from land plants. Of the 42 to 74 retained unique genes, only 16 protein genes, 15 tRNAs, and four rRNAs are commonly found. Several holoparasites retain ATP synthase genes with intact open reading frames, suggesting a prolonged function in these plants. The loss of photosynthesis alters the chromosomal architecture in that recombinogenic factors accumulate, fostering large-scale chromosomal rearrangements as functional reduction proceeds. The retention of DNA fragments is strongly influenced by both their proximity to genes under selection and the co-occurrence with those in operons, indicating complex constraints beyond gene function that determine the evolutionary survival time of plastid regions in nonphotosynthetic plants.     

Disintegration of the scrophulariaceae

A molecular systematic study of Scrophulariaceae sensu lato using DNA sequences of three plastid genes (rbcL, ndhF, and rps2) revealed at least five distinct monophyletic groups. Thirty-nine genera representing 24 tribes of the Scrophulariaceae s.l. (sensu lato) were analyzed along with representatives of 15 other families of Lamiales. The Scrophulariaceae s.s. (sensu stricto) include part or all of tribes Aptosimeae, Hemimerideae, Leucophylleae, Manuleae, Selagineae, and Verbasceae (= Scrophularieae) and the conventional families Buddlejaceae and Myoporaceae. Veronicaceae includes all or part of tribes Angelonieae, Antirrhineae, Cheloneae, Digitaleae, and Gratioleae and the conventional families Callitrichaceae, Globulariaceae, Hippuridaceae, and Plantaginaceae. The Orobanchaceae include tribes Buchnereae, Rhinantheae, and the conventional Orobanchaceae. All sampled members of Orobanchaceae are parasitic, except Lindenbergia, which is sister to the rest of the family. Family Calceolariaceae Olmstead is newly erected herein to recognize the phylogenetic distinctiveness of tribe Calceolarieae. The Calceolariaceae are close to the base of the Lamiales. The Stilbaceae are expanded by the inclusion of Halleria. Mimulus does not belong in any of these five groups.     

Global versus Chinese perspectives on the phylogeny of the N-fixing clade

There has been increasing interest in integrating a regional tree of life with community assembly rules in the ecological research. This raises questions regarding the impacts of taxon sampling strategies at the regional versus global scales on the topology. To address this concern, we constructed two trees for the nitrogen-fixing clade: (i) a genus-level global tree including 1023 genera; and (ii) a regional tree comprising 303 genera, with taxon sampling limited to China. We used the supermatrix approach and performed maximum likelihood analyses on combined matK, rbcL, and trnL-F plastid sequences. We found that the topology of the global and the regional tree of the N-fixing clade were generally congruent. However, whereas relationships among the four orders obtained with the global tree agreed with the accepted topology obtained in focused analyses with more genes, the regional topology obtained different relationships, albeit weakly supported. At a finer scale, the phylogenetic position of the family Myricaceae was found to be sensitive to sampling density. We expect that internal support throughout the phylogeny could be improved with denser taxon sampling. The taxon sampling approach (global vs. regional) did not have a major impact on fine-level branching patterns of the N-fixing clade. Thus, a well-resolved phylogeny with relatively dense taxon sampling strategy at the regional scale appears, in this case, to be a good representation of the overall phylogenetic pattern and could be used in ecological research. Otherwise, the regional tree should be adjusted according to the correspondingly reliable global tree.     

Multi‐locus phylogenetics,lineage sorting,and reticulation in Pinus subsection Australes

Premise of the Study

Both incomplete lineage sorting and reticulation have been proposed as causes of phylogenetic incongruence. Disentangling these factors may be most difficult in long‐lived, wind‐pollinated plants with large population sizes and weak reproductive barriers.

Methods

We used solution hybridization for targeted enrichment and massive parallel sequencing to characterize low‐copy‐number nuclear genes and high‐copy‐number plastomes (Hyb‐Seq) in 74 individuals of Pinus subsection Australes, a group of ~30 New World pine species of exceptional ecological and economic importance. We inferred relationships using methods that account for both incomplete lineage sorting and reticulation.

Key Results

Concatenation‐ and coalescent‐based trees inferred from nuclear genes mainly agreed with one another, but they contradicted the plastid DNA tree in recovering the Attenuatae (the California closed‐cone pines) and Oocarpae (the egg‐cone pines of Mexico and Central America) as monophyletic and the Australes sensu stricto (the southern yellow pines) as paraphyletic to the Oocarpae. The plastid tree featured some relationships that were discordant with morphological and geographic evidence and species limits. Incorporating gene flow into the coalescent analyses better fit the data, but evidence supporting the hypothesis that hybridization explains the non‐monophyly of the Attenuatae in the plastid tree was equivocal.

Conclusions

Our analyses document cytonuclear discordance in Pinus subsection Australes. We attribute this discordance to ancient and recent introgression and present a phylogenetic hypothesis in which mostly hierarchical relationships are overlain by gene flow.     

Passiflora plastome sequencing reveals widespread genomic rearrangements

Although past studies have included Passiflora among angiosperm lineages with highly rearranged plastid genomes (plastomes), knowledge about plastome organization in the genus is limited. So far only one draft and one complete plastome have been published. Expanded sampling of Passiflora plastomes is needed to understand the extent of the genomic rearrangement in the genus, which is also unusual in having biparental plastid inheritance and plastome‐genome incompatibility. We sequenced 15 Passiflora plastomes using either Illumina paired‐end or shotgun cloning and Sanger sequencing approaches. Assembled plastomes were annotated using Dual Organellar GenoMe Annotator (DOGMA) and tRNAscan‐SE. The Populus trichocarpa plastome was used as a reference to estimate genomic rearrangements in Passiflora by performing whole genome alignment in progressiveMauve. The phylogenetic distribution of rearrangements was plotted on the maximum likelihood tree generated from 64 plastid encoded protein genes. Inverted repeat (IR) expansion/contraction and loss of the two largest hypothetical open reading frames, ycf1 and ycf2, account for most plastome size variation, which ranges from 139 262 base pairs (bp) in P. biflora to 161 494 bp in P. pittieri. Passiflora plastomes have experienced numerous inversions, gene and intron losses along with multiple independent IR expansions and contractions resulting in a distinct organization in each of the three subgenera examined. Each Passiflora subgenus has a unique plastome structure in terms of gene content, order and size. The phylogenetic distribution of rearrangements shows that Passiflora has experienced widespread genomic changes, suggesting that such events may not be reliable phylogenetic markers.     

基于树木年轮学的古树树龄估算——以敦煌市香水梨为例

树龄是古树鉴定的唯一标准。本研究采用树轮年代学方法调查敦煌市莫高镇窦家墩村四组(原杨家园子)梨园中最老的香水梨树的树龄。因老梨树均腐朽中空,调查同时采集不同树龄的香水梨树样芯,建立树木径向生长生理年龄曲线,据此估算缺失段的树轮数。结合准确树轮计数,复原香水梨树树龄。结果表明: 梨园13棵老香水梨树的平均树龄为(280±35)年。由此证明,《敦煌志》中记载的古梨树是真实存在的。按树龄来说,该批梨树至少为三级古树,应妥善保护。本研究为解决敦煌市古梨园存留争议提供了详实的科学依据。     

Phylogenetic relationships in Epidendroideae (Orchidaceae), one of the great flowering plant radiations: progressive specialization and diversification

Background and Aims The largest subfamily of orchids, Epidendroideae, represents one of the most significant diversifications among flowering plants in terms of pollination strategy, vegetative adaptation and number of species. Although many groups in the subfamily have been resolved, significant relationships in the tree remain unclear, limiting conclusions about diversification and creating uncertainty in the classification. This study brings together DNA sequences from nuclear, plastid and mitochrondrial genomes in order to clarify relationships, to test associations of key characters with diversification and to improve the classification.Methods Sequences from seven loci were concatenated in a supermatrix analysis for 312 genera representing most of epidendroid diversity. Maximum-likelihood and parsimony analyses were performed on this matrix and on subsets of the data to generate trees and to investigate the effect of missing values. Statistical character-associated diversification analyses were performed.Key Results Likelihood and parsimony analyses yielded highly resolved trees that are in strong agreement and show significant support for many key clades. Many previously proposed relationships among tribes and subtribes are supported, and some new relationships are revealed. Analyses of subsets of the data suggest that the relatively high number of missing data for the full analysis is not problematic. Diversification analyses show that epiphytism is most strongly associated with diversification among epidendroids, followed by expansion into the New World and anther characters that are involved with pollinator specificity, namely early anther inflexion, cellular pollinium stalks and the superposed pollinium arrangement.Conclusions All tested characters show significant association with speciation in Epidendroideae, suggesting that no single character accounts for the success of this group. Rather, it appears that a succession of key features appeared that have contributed to diversification, sometimes in parallel.     

Molecular markers for phylogenetic applications derived from comparative plastome analysis of Prunus species

Interspecific and intergeneric relationships of Prunus s.l. are still unclear due to low levels of genetic variation among species, and resulting partially unresolved phylogenetic inferences. Here we sequenced and compared six complete plastomes from two subgenera of Prunus in order to choose molecular markers to increase the amount of genetic variation suitable for inference of Prunus phylogeny. The plastomes range between 157 817 and 158 995 bp in length, and we found different levels of inverted repeat (IR) contraction among the three sampled subgenera of Prunus s.l. Most regions in Prunus plastomes considered individually provide low phylogenetic resolution at the subgenus or species level compared to a tree constructed using all 78 coding regions combined. We compared levels of variation among 206 coding regions and noncoding (intergenic and intron) plastid regions and inferred phylogenies from each region considered individually. We then chose using two regions together for future studies of relationships in Prunus, ycf1 and trnT-L, that display high to moderate levels of variation among coding and intergenic regions, respectively, and that individually permit inference of resolved species-level trees in Prunus with moderate to strong branch support. Considered together, these two regions allow inference of the same topology of Prunus inferred using all coding plastid regions combined, with comparable levels of tree support to the full plastome set. These two loci should therefore be useful as a plastid phylogenetic marker set for further inference of relationships within Prunus s.l.     

Prospects on the evolutionary mitogenomics of plants: A case study on the olive family (Oleaceae)

The mitogenome is rarely used to reconstruct the evolutionary history of plants, contrary to nuclear and plastid markers. Here, we evaluate the usefulness of mitochondrial DNA for molecular evolutionary studies in Oleaceae, in which cases of cytoplasmic male sterility (CMS) and of potentially contrasted organelle inheritance are known. We compare the diversity and the evolution of mitochondrial and chloroplast genomes by focusing on the olive complex and related genera. Using high‐throughput techniques, we reconstructed complete mitogenomes (ca. 0.7 Mb) and plastomes (ca. 156 kb) for six olive accessions and one Chionanthus. A highly variable organization of mitogenomes was observed at the species level. In olive, two specific chimeric genes were identified in the mitogenome of lineage E3 and may be involved in CMS. Plastid‐derived regions (mtpt) were observed in all reconstructed mitogenomes. Through phylogenetic reconstruction, we demonstrate that multiple integrations of mtpt regions have occurred in Oleaceae, but mtpt regions shared by all members of the olive complex derive from a common ancestor. We then assembled 52 conserved mitochondrial gene regions and complete plastomes of ten additional accessions belonging to tribes Oleeae, Fontanesieae and Forsythieae. Phylogenetic congruence between topologies based on mitochondrial regions and plastomes suggests a strong disequilibrium linkage between both organellar genomes. Finally, while phylogenetic reconstruction based on plastomes fails to resolve the evolutionary history of maternal olive lineages in the Mediterranean area, their phylogenetic relationships were successfully resolved with complete mitogenomes. Overall, our study demonstrates the great potential of using mitochondrial DNA in plant phylogeographic and metagenomic studies.     

利用低拷贝核基因重建菊科紫菀亚科族间系统发育关系

紫菀亚科(Asteroideae)是菊科最大的一个亚科, 包含的种数多于被子植物的绝大多数科。目前, 紫菀亚科族间的系统发育关系主要依赖于叶绿体基因信息, 但是叶绿体基因为单亲遗传, 并不能完整反映进化历史。鉴于杂交现象在菊科普遍存在, 故利用核基因可以反映更完整的紫菀亚科进化历史。该研究首次使用从转录组数据(20个新测+11个从NCBI数据库下载)中筛选出的47个直系同源低拷贝核基因来研究紫菀亚科的系统发育关系, 共选取了29个物种, 代表了紫菀亚科20个族中的13个族。用超矩阵分析方法和溯祖推测分析方法各获得了1个稳定的紫菀亚科系统树, 每个树上绝大多数分支都得到了高度支持, 且2个树之间没有明显的冲突。新的紫菀亚科族间系统发育关系揭示了千里光超族应并入紫菀超族, 春黄菊族可能是千里光族与紫菀族杂交起源的, 金鸡菊族很可能也是杂交起源的。该研究结果显示低拷贝核基因可以更好地解决科以下分类阶元的系统发育关系, 对菊科乃至被子植物其它科的系统发育研究具有重要的借鉴意义。     

The first complete plastid genome of Burmannia disticha L. from the mycoheterotrophic monocot family Burmanniaceae

Burmanniaceae is one major group within the monocot order Dioscoreales that has not had its plastome sequenced. Members of Burmanniaceae are mostly achlorophyllous, although the genus Burmannia also includes autotrophs. Here, we report sequencing and analysis of the first Burmanniaceae plastid genome from Burmannia disticha L.. This plastome is 157,480 bp and was assembled as a circular sequence with the typical quadripartite structure of plant plastid genomes. This plastome has a regular number of potentially functional genes with a total of 111, including 78 protein coding genes, 4 ribosomal RNA (rRNA) genes, and 29 tRNA genes. The ratio of the total length of genic:intergenic DNA is 1.58:1, and the mean length of intergenic regions is 398 bp, the longest being 1918 bp. The overall GC content of the B. disticha plastome is 34.90%, and the IR regions in B. disticha are more GC rich (39.50%) than the LSC (32.30%) and SSC (28.80%) regions. Phylogenetic analysis of protein-coding sequences from plastomes of related species in the order Dioscoreales support a clade comprising Burmanniaceae and Dioscoreaceae. This phylogenetic placement is congruent with previous findings based on nuclear and mitochondrial evidence.