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

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

Development of novel low-copy nuclear markers for Hieraciinae (Asteraceae) and their perspective for other tribes

? Premise of the study: The development of three low-copy nuclear markers for low taxonomic level phylogenies in Asteraceae with emphasis on the subtribe Hieraciinae is reported. ? Methods and Results: Marker candidates were selected by comparing a Lactuca complementary DNA (cDNA) library with public DNA sequence databases. Interspecific variation and phylogenetic signal of the selected genes were investigated for diploid taxa from the subtribe Hieraciinae and compared to a reference phylogeny. Their ability to cross-amplify was assessed for other Asteraceae tribes. All three markers had higher variation (2.1-4.5 times) than the internal transcribed spacer (ITS) in Hieraciinae. Cross-amplification was successful in at least seven other tribes of the Asteraceae. Only three cases indicating the presence of paralogs or pseudogenes were detected. ? Conclusions: The results demonstrate the potential of these markers for phylogeny reconstruction in the Hieraciinae as well as in other Asteraceae tribes, especially for very closely related species.     

Senecio changii (Asteraceae: Senecioneae), a New Species from Sichuan,China

Senecio changii (Asteraceae: Senecioneae), a new species from Muli, Sichuan, southwestern China, is described. It is distinguished in Chinese Senecio s.s. by having lyrate-pinnatisect to pinnatisect leaves and a single terminal large discoid capitulum which is somewhat nodding. Evidence from floral micromorphology, karyology and molecular phylogenetic analyses based on the nuclear ITS/ETS sequence data all support its membership within Senecio s.s.     

Complete chloroplast genome sequences of Praxelis (Eupatorium catarium Veldkamp), an important invasive species

Praxelis (Eupatorium catarium Veldkamp) is a new hazardous invasive plant species that has caused serious economic losses and environmental damage in the Northern hemisphere tropical and subtropical regions. Although previous studies focused on detecting the biological characteristics of this plant to prevent its expansion, little effort has been made to understand the impact of Praxelis on the ecosystem in an evolutionary process. The genetic information of Praxelis is required for further phylogenetic identification and evolutionary studies. Here, we report the complete Praxelis chloroplast (cp) genome sequence. The Praxelis chloroplast genome is 151,410 bp in length including a small single-copy region (18,547 bp) and a large single-copy region (85,311 bp) separated by a pair of inverted repeats (IRs; 23,776 bp). The genome contains 85 unique and 18 duplicated genes in the IR region. The gene content and organization are similar to other Asteraceae tribe cp genomes. We also analyzed the whole cp genome sequence, repeat structure, codon usage, contraction of the IR and gene structure/organization features between native and invasive Asteraceae plants, in order to understand the evolution of organelle genomes between native and invasive Asteraceae. Comparative analysis identified the 14 markers containing greater than 2% parsimony-informative characters, indicating that they are potential informative markers for barcoding and phylogenetic analysis. Moreover, a sister relationship between Praxelis and seven other species in Asteraceae was found based on phylogenetic analysis of 28 protein-coding sequences. Complete cp genome information is useful for plant phylogenetic and evolutionary studies within this invasive species and also within the Asteraceae family.     

(Un)Targeted Metabolomics in Asteraceae: Probing the Applicability of Essential‐Oil Profiles of Senecio L. (Senecioneae) Taxa in Chemotaxonomy

The possible applicability of (un)targeted metabolomics (volatile metabolites) for revealing taxonomic/evolutionary relationships among Senecio L. species (Asteraceae; tribe Senecioneae) was explored. Essential‐oil compositional data of selected Senecio/Senecioneae/Asteraceae taxa (93 samples in total) were mutually compared by means of multivariate statistical analysis (MVA), i.e., agglomerative hierarchical clustering and principal component analysis. The MVA input data set included the very first compositional data on the essential oil extracted from the aerial parts of S. viscosus L. as well as on four different Serbian populations of S. vernalis Waldst . & Kit . (oils from aerial parts and roots; eight samples in total). This metabolomic screening of Senecio/Senecioneae/Asteraceae species (herein presented results and data from the literature) pointed to short‐chain alk‐1‐enes (e.g., oct‐1‐ene, non‐1‐ene, and undec‐1‐ene), with up to now restricted general occurrence in Plantae, as characteristic chemotaxonomic markers/targets for future metabolomic studies of Senecio/Senecioneae taxa. The MVA additionally showed that the evolution of the terpene metabolism (volatile mono‐ and sesquiterpenoids) within the Asteraceae tribe Senecioneae was not genera specific. However, the MVA did confirm plant‐organ specific production/accumulation of volatiles within S. vernalis and suggested the existence of at least two volatile chemotypes for this species.     

The ITS region of groundwater amphipods: length,secondary structure and phylogenetic information content in Crangonyctoids and Niphargids

The phylogenetic analysis of groundwater amphipods is challenging due to the lack of suitable morphological characters. However, molecular phylogenies based on the 18S and 28S nuclear genes of two Crangonyctoidea species endemic to Iceland, Crymostygius thingvallensis and Crangonyx islandicus, support the taxonomy of these species on the basis of morphological characters. Molecular analyses suggest that the genus Crangonyx is paraphyletic, with the species that is found in Eurasia being highly divergent genetically from the species present in North America and Iceland. Studies of the phylogenetic relationships within the genus Niphargus also warrant further work. The nuclear ITS2 region has recently been proposed as a barcoding marker for plants and animals. In addition, ITS2 has been used to build phylogenies at high taxonomic levels by including its secondary structure. In this study, we want to evaluate the applicability of the ITS region for this group of species and describe its characteristics. The taxonomy of C. thingvallensis, as well as the paraphyly of the genus Crangonyx, is supported herein by phylogenies based on the ITS2 variation. The secondary structure and the length of the ITS2 sequences of the Crangonyctoidea and the Niphargidae species studied are highly variable and are characterized by duplications. The ITS2 sequence of Niphargus plateaui is the longest metazoan sequence deposited in the ITS2 database so far. Although saturation was observed in the nucleotide variation of this marker, the addition of the secondary structure information for the reconstruction of the phylogeny did not add support to the phylogenetic trees. The ITS1 region, which is known to be more variable than ITS2 and bears a large duplication within C. islandicus, was found to be less useful for phylogenetic reconstruction.     

Phylogenetic relationships of New Zealand Asteraceae inferred from ITS sequences

Forty-five sequences from members of all genera of Asteraceae indigenous to New Zealand and 50 published sequences representing the tribal diversity in the family were analyzed to assess the utility of ITS sequences to resolve phylogenetic relationships. Previous studies using chloroplast DNA sequences and morphology provided support for several clades in the Asteraceae, yet the relationships among some of these were uncertain. The results from ITS analysis were largely consistent with these earlier studies. The New Zealand species are included in at least six clades, most of these corresponding to recognized tribes. Our results have also clarified the tribal affinities of a few anomalous genera. Haastia, previously aligned with the Gnaphalieae or the Astereae, is nested in the Senecioneae. Centipeda, previously included in the Astereae or Anthemideae, emerges near the Heliantheae. The relationships of Abrotanella remain unresolved. Received August 8, 2001 Accepted November 6, 2001     

Identification of shared single copy nuclear genes in Arabidopsis, Populus, Vitis and Oryza and their phylogenetic utility across various taxonomic levels

Background

Although the overwhelming majority of genes found in angiosperms are members of gene families, and both gene- and genome-duplication are pervasive forces in plant genomes, some genes are sufficiently distinct from all other genes in a genome that they can be operationally defined as 'single copy'. Using the gene clustering algorithm MCL-tribe, we have identified a set of 959 single copy genes that are shared single copy genes in the genomes of Arabidopsis thaliana, Populus trichocarpa, Vitis vinifera and Oryza sativa. To characterize these genes, we have performed a number of analyses examining GO annotations, coding sequence length, number of exons, number of domains, presence in distant lineages, such as Selaginella and Physcomitrella, and phylogenetic analysis to estimate copy number in other seed plants and to demonstrate their phylogenetic utility. We then provide examples of how these genes may be used in phylogenetic analyses to reconstruct organismal history, both by using extant coverage in EST databases for seed plants and de novo amplification via RT-PCR in the family Brassicaceae.

Results

There are 959 single copy nuclear genes shared in Arabidopsis, Populus, Vitis and Oryza ["APVO SSC genes"]. The majority of these genes are also present in the Selaginella and Physcomitrella genomes. Public EST sets for 197 species suggest that most of these genes are present across a diverse collection of seed plants, and appear to exist as single or very low copy genes, though exceptions are seen in recently polyploid taxa and in lineages where there is significant evidence for a shared large-scale duplication event. Genes encoding proteins localized in organelles are more commonly single copy than expected by chance, but the evolutionary forces responsible for this bias are unknown. Regardless of the evolutionary mechanisms responsible for the large number of shared single copy genes in diverse flowering plant lineages, these genes are valuable for phylogenetic and comparative analyses. Eighteen of the APVO SSC single copy genes were amplified in the Brassicaceae using RT-PCR and directly sequenced. Alignments of these sequences provide improved resolution of Brassicaceae phylogeny compared to recent studies using plastid and ITS sequences. An analysis of sequences from 13 APVO SSC genes from 69 species of seed plants, derived mainly from public EST databases, yielded a phylogeny that was largely congruent with prior hypotheses based on multiple plastid sequences. Whereas single gene phylogenies that rely on EST sequences have limited bootstrap support as the result of limited sequence information, concatenated alignments result in phylogenetic trees with strong bootstrap support for already established relationships. Overall, these single copy nuclear genes are promising markers for phylogenetics, and contain a greater proportion of phylogenetically-informative sites than commonly used protein-coding sequences from the plastid or mitochondrial genomes.

Conclusions

Putatively orthologous, shared single copy nuclear genes provide a vast source of new evidence for plant phylogenetics, genome mapping, and other applications, as well as a substantial class of genes for which functional characterization is needed. Preliminary evidence indicates that many of the shared single copy nuclear genes identified in this study may be well suited as markers for addressing phylogenetic hypotheses at a variety of taxonomic levels.     

Phylogenetic congruence between Mollitrichosiphum (Aphididae: Greenideinae) and Buchnera indicates insect–bacteria parallel evolution

We wanted to test whether Mollitrichosiphum, an aphid genus with life cycles on subtropical woody host plants, and Buchnera, the primary endosymbiont of aphids, evolve in parallel. We used three aphid genes (mitochondrial COI, cytochrome oxidase subunit I and Cytb, cytochrome b; nuclear EF1α, translation elongation factor 1 alpha) and two Buchnera genes (16S rDNA; gnd, gluconate‐6‐phosphate dehydrogenase) to reconstruct phylogenies. The congruence between the phylogenetic trees of aphids and Buchnera was then measured. The results present phylogenetic evidence for the parallel evolution of Mollitrichosiphum and Buchnera at the intraspecific as well as the interspecific levels. Our results support the possibility of using endosymbiont genes to study host evolutionary history and biogeographical patterns. We also investigated the usability of the Buchnera gnd gene as a barcoding marker for aphid identification.     

Lines of Evidence for Horizontal Gene Transfer of a Phenazine Producing Operon into Multiple Bacterial Species

Phenazines are secondary metabolites with broad-spectrum antibiotic activity against bacteria, fungi, and eukaryotes. In pseudomonad species, a conserved seven-gene phenazine operon (phzABCDEFG) is required for the conversion of chorismic acid to the broad-spectrum antibiotic phenazine-1-carboxylate. Previous analyses of genes involved in phenazine production from nonpseudomonad species uncovered a high degree of sequence similarity to pseudomonad homologues. The analyses undertaken in this study wished to eluciadate the evolutionary history of genes involved in the production of phenazines. Furthermore, I wanted to determine if the phenazine operon has been transferred through horizontal gene transfer. Analyses of GC content, codon usage patterns, frequency of 3:1 dinucleotides, sequence similarities, and phylogenetic reconstructions were undertaken to map the evolutionary history of phenazine genes from multiple bacterial species. Patchy phyletic distribution, high sequence similarities, and phylogenetic evidence infer that pseudomonad, Streptomyces cinnamonensis, Pantoea agglomerans, Burkholderia cepacia, Pectobacterium atrosepticum, Brevibacterium linens, and Mycobacterium abscessus species all contain a phenazine operon which has most likely been transferred among these species through horizontal gene transfer. The acquisition of an antibiotic-associated operon is significant, as it may increase the relative fitness of the recipient species. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Patterns and causes of incongruence between plastid and nuclear Senecioneae (Asteraceae) phylogenies

One of the longstanding questions in phylogenetic systematics is how to address incongruence among phylogenies obtained from multiple markers and how to determine the causes. This study presents a detailed analysis of incongruent patterns between plastid and ITS/ETS phylogenies of Tribe Senecioneae (Asteraceae). This approach revealed widespread and strongly supported incongruence, which complicates conclusions about evolutionary relationships at all taxonomic levels. The patterns of incongruence that were resolved suggest that incomplete lineage sorting (ILS) and/or ancient hybridization are the most likely explanations. These phenomena are, however, extremely difficult to distinguish because they may result in similar phylogenetic patterns. We present a novel approach to evaluate whether ILS can be excluded as an explanation for incongruent patterns. This coalescence-based method uses molecular dating estimates of the duration of the putative ILS events to determine if invoking ILS as an explanation for incongruence would require unrealistically high effective population sizes. For four of the incongruent patterns identified within the Senecioneae, this approach indicates that ILS cannot be invoked to explain the observed incongruence. Alternatively, these patterns are more realistically explained by ancient hybridization events.     

Generating single-copy nuclear gene data for a recent adaptive radiation

Recent adaptive radiations provide an exceptional opportunity to understand the processes of speciation and adaptation. However, reconstructing the phylogenetic history of recent and rapidly evolving clades often requires the use of multiple, independent gene genealogies. Nuclear introns are an obvious source of the necessary data but their use is often limited because degenerate primers can amplify paralogous loci. To identify PCR primers for a large number of loci in an especially rapid adaptive radiation, that of the flowering plant genus Aquilegia (Ranunculaceae), we developed an efficient method for amplifying multiple single-copy nuclear loci by sequencing a modest number of clones from a cDNA library and designing PCR primers; with one primer anchored in the 3' untranslated region (3'-UTR) and one primer in the coding region of each gene. Variation between paralogous loci evolves more quickly in 3'-UTR regions compared to adjacent exons, and therefore we achieved high specificity for isolating orthologous loci. Furthermore, we were able to identify genes containing large introns by amplifying genes from genomic DNA and comparing the PCR product size to that predicted from their cDNA sequence. In Aquilegia eight out of eleven loci were isolated with this method and six of these loci had introns. Among four genes sequenced for samples spanning the phylogenetic breadth of the genus, we found sequence variation at levels similar to that observed in ITS, further supporting the recent and rapid radiation in Aquilegia. We assessed the orthology of amplification products by phylogenetic congruence among loci, the presence of two well established phylogenetic relationships, and similarity among loci for levels of sequence variation. Higher levels of variation among samples for one locus suggest possible paralogy. Overall, this method provides an efficient means of isolating predominantly single-copy loci from both low and high-copy gene families, providing ample nuclear variation for reconstructing species-level phylogenies in non-model taxa.     

Frequent chloroplast capture among Isodon (Lamiaceae) species in Japan revealed by phylogenies based on variation in chloroplast and nuclear DNA

A recent phylogenetic study based only on chloroplast DNA (cpDNA) variation revealed that populations of an Isodon species are frequently embedded paraphyletically among other Isodon species. This phylogenetic discrepancy between species taxonomy and molecular phylogeny was considered to have resulted from chloroplast DNA captures and/or incomplete lineage sorting. To elucidate which of these factors was mainly responsible for the observed phylogenetic pattern, we performed phylogenetic analyses of multiple populations of Isodon species in Japan using cpDNA variation, three single-copy nuclear genes, and double-digest restriction-site-associated DNA sequencing (ddRAD-seq). Although a species often shared chlorotypes with other species, our phylogenetical analyses based on variation in the three single-copy nuclear genes and the ddRAD-seq data showed that most populations belonging to the same species were monophyletic at the species level, suggesting that chloroplast capture may have frequently occurred between Isodon species. Some populations of an intraspecific taxon were embedded paraphyletically within the species, regardless of the large amount of phylogenetic information in nuclear DNA; this incongruity may have resulted from incomplete lineage sorting.     

Phylogeny of the New World diploid cottons (Gossypium L., Malvaceae) based on sequences of three low-copy nuclear genes

American diploid cottons (Gossypium L., subgenus Houzingenia Fryxell) form a monophyletic group of 13 species distributed mainly in western Mexico, extending into Arizona, Baja California, and with one disjunct species each in the Galapagos Islands and Peru. Prior phylogenetic analyses based on an alcohol dehydrogenase gene (AdhA) and nuclear ribosomal DNA indicated the need for additional data from other molecular markers to resolve phylogenetic relationships within this subgenus. Toward this end, we sequenced three nuclear genes, the anonymous locus A1341, an alcohol dehydrogenase gene (AdhC), and a cellulose synthase gene (CesA1b). Independent and combined analyses resolved clades that are congruent with current taxonomy and previous phylogenies. Our analyses diagnose at least two long distance dispersal events from the Mexican mainland to Baja California, following a rapid radiation of the primary lineages early in the diversification of the subgenus. Molecular data support the proposed recognition of a new species closely related to Gossypium laxum that was recently collected in Mexico.     

Comparative plastome analysis of Blumea,with implications for genome evolution and phylogeny of Asteroideae

The genus Blumea (Asteroideae, Asteraceae) comprises about 100 species, including herbs, shrubs, and small trees. Previous studies have been unable to resolve taxonomic issues and the phylogeny of the genus Blumea due to the low polymorphism of molecular markers. Therefore, suitable polymorphic regions need to be identified. Here, we de novo assembled plastomes of the three Blumea species B. oxyodonta, B. tenella, and B. balsamifera and compared them with 26 other species of Asteroideae after correction of annotations. These species have quadripartite plastomes with similar gene content, genome organization, and inverted repeat contraction and expansion comprising 113 genes, including 80 protein‐coding, 29 transfer RNA, and 4 ribosomal RNA genes. The comparative analysis of codon usage, amino acid frequency, microsatellite repeats, oligonucleotide repeats, and transition and transversion substitutions has revealed high resemblance among the newly assembled species of Blumea. We identified 10 highly polymorphic regions with nucleotide diversity above 0.02, including rps16‐trnQ, ycf1, ndhF‐rpl32, petN‐psbM, and rpl32‐trnL, and they may be suitable for the development of robust, authentic, and cost‐effective markers for barcoding and inference of the phylogeny of the genus Blumea. Among these highly polymorphic regions, five regions also co‐occurred with oligonucleotide repeats and support use of repeats as a proxy for the identification of polymorphic loci. The phylogenetic analysis revealed a close relationship between Blumea and Pluchea within the tribe Inuleae. At tribe level, our phylogeny supports a sister relationship between Astereae and Anthemideae rooted as Gnaphalieae, Calenduleae, and Senecioneae. These results are contradictory to recent studies which reported a sister relationship between “Senecioneae and Anthemideae” and “Astereae and Gnaphalieae” or a sister relationship between Astereae and Gnaphalieae rooted as Calenduleae, Anthemideae, and then Senecioneae using nuclear genome sequences. The conflicting phylogenetic signals observed at the tribal level between plastidt and nuclear genome data require further investigation.     

Discordant mitochondrial and nuclear gene phylogenies in emydid turtles: implications for speciation and conservation

Do phylogenies and branch lengths based on mitochondrial DNA (mtDNA) provide a reasonable approximation to those based on multiple nuclear loci? In the present study, we show widespread discordance between phylogenies based on mtDNA (two genes) and nuclear DNA (nucDNA; six loci) in a phylogenetic analysis of the turtle family Emydidae. We also find an unusual type of discordance involving the unexpected homogeneity of mtDNA sequences across species within genera. Of the 36 clades in the combined nucDNA phylogeny, 24 are contradicted by the mtDNA phylogeny, and six are strongly contested by each data set. Two genera (Graptemys, Pseudemys) show remarkably low mtDNA divergence among species, whereas the combined nuclear data show deep divergences and (for Pseudemys) strongly supported clades. These latter results suggest that the mitochondrial data alone are highly misleading about the rate of speciation in these genera and also about the species status of endangered Graptemys and Pseudemys species. In addition, despite a strongly supported phylogeny from the combined nuclear genes, we find extensive discordance between this tree and individual nuclear gene trees. Overall, the results obtained illustrate the potential dangers of making inferences about phylogeny, speciation, divergence times, and conservation from mtDNA data alone (or even from single nuclear genes), and suggest the benefits of using large numbers of unlinked nuclear loci. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 99 , 445–461.     

Optimization of Multilocus Sequence Analysis for Identification of Species in the Genus Vibrio

Multilocus sequence analysis (MLSA) is an important method for identification of taxa that are not well differentiated by 16S rRNA gene sequences alone. In this procedure, concatenated sequences of selected genes are constructed and then analyzed. The effects that the number and the order of genes used in MLSA have on reconstruction of phylogenetic relationships were examined. The recA, rpoA, gapA, 16S rRNA gene, gyrB, and ftsZ sequences from 56 species of the genus Vibrio were used to construct molecular phylogenies, and these were evaluated individually and using various gene combinations. Phylogenies from two-gene sequences employing recA and rpoA in both possible gene orders were different. The addition of the gapA gene sequence, producing all six possible concatenated sequences, reduced the differences in phylogenies to degrees of statistical (bootstrap) support for some nodes. The overall statistical support for the phylogenetic tree, assayed on the basis of a reliability score (calculated from the number of nodes having bootstrap values of ≥80 divided by the total number of nodes) increased with increasing numbers of genes used, up to a maximum of four. No further improvement was observed from addition of the fifth gene sequence (ftsZ), and addition of the sixth gene (gyrB) resulted in lower proportions of strongly supported nodes. Reductions in the numbers of strongly supported nodes were also observed when maximum parsimony was employed for tree construction. Use of a small number of gene sequences in MLSA resulted in accurate identification of Vibrio species.