Arnica (Asteraceae) phylogeny revisited using RPB2: complex patterns and multiple d-paralogues

The region coding for the second largest subunit of RNA polymerase II (RPB2) was explored for resolving interspecific relationships in Arnica and lower level taxa in general. The region between exons 17 and 23 was cloned and sequenced for 33 accessions of Arnica and four outgroup taxa. Three paralogues of the RPB2-d copy (RPB2-dA, B and C) were detected in Arnica and outgroup taxa, indicating that the duplications must have occurred before the divergence of Arnica. Parsimony and Bayesian analyses of separate alignments of the three copies reveal complex patterns in Arnica, likely reflecting a history of lineage sorting in combination with apomixis, polyploidization, and possibly hybridization. Cloned sequences of some taxa do not form monophyletic clades within paralogues, but form multiple strongly supported clades with sequences of other taxa. Some well supported groups are present in more than one paralogue and many groups are in line with earlier hypotheses regarding interspecific relationships within the genus. Low levels of homoplasy in combination with relatively high sequence variation indicates that the introns of the RPB2 region could be suitable for phylogenetic studies in low level taxonomy.     

Phylogenetic utility of the nuclear gene arginine decarboxylase: an example from Brassicaceae

Arginine decarboxylase (ADC) is an important enzyme in the production of putrescine and polyamines in plants. It is encoded by a single or low-copy nuclear gene that lacks introns in sequences studied to date. The rate of Adc amino acid sequence evolution is similar to that of ndhF for the angiosperm family studied. Highly conserved regions provide several target sites for PCR priming and sequencing and aid in nucleotide and amino acid sequence alignment across a range of taxonomic levels, while a variable region provides an increased number of potentially informative characters relative to ndhF for the taxa surveyed. The utility of the Adc gene in plant molecular systematic studies is demonstrated by analysis of its partial nucleotide sequences obtained from 13 representatives of Brassicaceae and 3 outgroup taxa, 2 from the mustard oil clade (order Capparales) and 1 from the related order Malvales. Two copies of the Adc gene, Adc1 and Adc2, are found in all members of the Brassicaceae studied to data except the basal genus Aethionema. The resulting Adc gene tree provides robust phylogenetic data regarding relationships within the complex mustard family, as well as independent support for proposed tribal realignments based on other molecular data sets such as those from chloroplast DNA.      

Phylogenetic relationships within the Gentianales based on NDHF and RBCL sequences, with particular reference to the Loganiaceae

Phylogenetic relationships in the Gentianales with focus on Loganiaceae sensu lato are evaluated using parsimony analyses of nucleotide sequence data from the plastid genes rbcL and ndhF. Inter- and intrafamilial relationships in the Gentianales, which consist of the families Apocynaceae (including Asclepiadaceae), Gelsemiaceae, Gentianaceae, Loganiaceae, and Rubiaceae, are studied and receive increased support from the combination of rbcL and ndhF data, which indicate that the family Rubiaceae forms the sister group to the successively nested Gentianaceae, Apocynaceae, and Loganiaceae, all of which are well supported. The family Gelsemiaceae forms a distinct, supported group sister to Apocynaceae. The Loganiaceae sensu stricto form a strongly supported group consisting of 13 genera: Antonia, Bonyunia, Gardneria, Geniostoma, Labordia, Logania, Mitrasacme, Mitreola, Neuburgia, Norrisia, Spigelia, Strychnos, and Usteria. These genera form two well-supported lineages. Several members of Loganiaceae sensu Leeuwenberg and Leenhouts, i.e., Androya, Peltanthera, Plocosperma, Polypremum, and Sanango are clearly not members of the Gentianales. The earlier exclusion of Buddlejaceae (including Buddleja, Emorya, Gomphostigma, and Nicodemia) as well as the reclassification of the genera Nuxia and Retzia to Stilbaceae of the Lamiales are all well supported.     

Migration of endpoints of two genes relative to boundaries between regions of the plastid genome in the grass family (Poaceae)

Overlapping genes occur widely in microorganisms and in some plastid genomes, but unique properties are observed when such genes span the boundaries between single-copy and repeat regions. The termini of ndhH and ndhF, situated near opposite ends of the small single-copy region (SSC) in the plastid genomes of grasses (Poaceae), have migrated repeatedly into and out of the adjacent inverted-repeat regions (IR). The two genes are transcribed in the same direction, and the 5' terminus of ndhH extends into the IR in some species, while the 3' terminus of ndhF extends into the IR in others. When both genes extend into the IR, portions of the genes overlap and are encoded by the same nucleotide positions. Fine-scale mapping of the SSC-IR junctions across a sample of 92 grasses and outgroups, integrated into a phylogenetic analysis, indicates that the earliest grasses resembled the related taxa Joinvillea (Joinvilleaceae) and Ecdeiocolea (Ecdeiocoleaceae), with ca. 180 nucleotides of ndhH extending into the IR, and with ndhF confined to the SSC. This structure is maintained in early-diverging grass lineages and in most species of the BEP clade. In the PACMAD clade, ndhH lies completely or nearly completely within the SSC, and ca. 20 nucleotides of ndhF extend into the IR. The nucleotide substitution rate has increased in the PACMAD clade in the portion of ndhH that has migrated into the SSC.     

Presequence Acquisition During Secondary Endocytobiosis and thePossible Role of Introns

Targeting of nucleus-encoded proteins into chloroplasts is mediated by N-terminal presequences. During evolution of plastids from formerly free-living cyanobacteria by endocytobiosis, genes for most plastid proteins have been transferred from the plastid genome to the nucleus and subsequently had to be equipped with such plastid targeting sequences. So far it is unclear how the gene domains coding for presequences and the respective mature proteins may have been assembled. While land plant plastids are supposed to originate from a primary endocytobiosis event (a prokaryotic cyanobacterium was taken up by a eukaryotic cell), organisms with secondary plastids like diatoms experienced a second endocytobiosis step involving a eukaryotic alga taken up by a eukaryotic host cell. In this group of algae, apparently most genes encoding chloroplast proteins have been transferred a second time (from the nucleus of the endosymbiont to the nucleus of the secondary host) and thus must have been equipped with additional targeting signals. We have analyzed cDNAs and the respective genomic DNA fragments of seven plastid preproteins from the diatom Phaeodactylum tricornutum. In all of these genes we found single spliceosomal introns, generally located within the region coding for the N-terminal plastid targeting sequences or shortly downstream of it. The positions of the introns can be related to the putative phylogenetic histories of the respective genes, indicating that the bipartite targeting sequences in these secondary algae might have evolved by recombination events via introns.The nucleotide sequences have been deposited at Genbank under accession numbers AY191862, AY191863, AY191864, AY191865, AY191866, AY191867, and AY191868.     

rRNA二级结构序列用于真菌系统学研究的方法初探

本文首次利用核酸二级结构特征代替核酸碱基作为探讨类群之间亲缘关系的信号,构建了基于结构特征的子囊菌部分类群的系统进化树。该方法以S(规范的碱基对),Q(不规范的碱基对),I(单链),B(侧环),M(多分枝环)和H(发卡结构)为代码将二级结构特征区分为6种不同的亚结构类型,然后将二级结构特征转换为结构序列,并进行结构序列分析。该方法使rRNA不只局限于碱基比较,拓展了其应用范围,为揭示分子的功能与进化的关系提供了线索。结果表明,结构序列分析可用于子囊菌的系统学研究;相对于核酸序列分析,结构分析的结果似乎更加清晰地体现子囊果的演化过程。     

Nucleotide sequence of the B1-hordein gene of barley Hordeum vulgare L

The gene encoding B1 hordein of Hordeum vulgare (cv. Donetsky 4) was cloned and entirely sequenced. It contains no introns and codes for of 293 amino acids long polypeptide with molecular weight 33418. Our clone differs from the previously sequenced B1 hordein genes in some positions within the coding region (there are 4 nucleotide changes and a 12 bp deletion, as compared with the pB11 cDNA clone, and 5 nucleotide changes, as compared with the pBHP184 genomic clone). These changes result in polymorphism of amino acid sequences at 5 positions. 5'-flanking region contains putative regulatory and promoter sequences and differs from that of the pBHP184 clone in 3 positions.     

Elongation factor-1 alpha occurs as two copies in bees: implications for phylogenetic analysis of EF-1 alpha sequences in insects

We report the complete sequence of a paralogous copy of elongation factor-1 alpha (EF-1 alpha) in the honeybee, Apis mellifera (Hymenoptera: Apidae). This copy differs from a previously described copy in the positions of five introns and in 25% of the nucleotide sites in the coding regions. The existence of two paralogous copies of EF-1 alpha in Drosophila and Apis suggests that two copies of EF-1 alpha may be widespread in the holometabolous insect orders. To distinguish between a single, ancient gene duplication and parallel, independent fly and bee gene duplications, we performed a phylogenetic analysis of hexapod EF-1 alpha sequences. Unweighted parsimony analysis of nucleotide sequences suggests an ancient gene duplication event, whereas weighted parsimony analysis of nucleotides and unweighted parsimony analysis of amino acids suggests the contrary: that EF-1 alpha underwent parallel gene duplications in the Diptera and the Hymenoptera. The hypothesis of parallel gene duplication is supported both by congruence among nucleotide and amino acid data sets and by topology-dependent permutation tail probability (T-PTP) tests. The resulting tree topologies are also congruent with current views on the relationships among the holometabolous orders included in this study (Diptera, Hymenoptera, and Lepidoptera). More sequences, from diverse orders of holometabolous insects, will be needed to more accurately assess the historical patterns of gene duplication in EF-1 alpha.      

红角辉蝽和紫翅果蝽线粒体完整基因组测序及蝽亚科系统发育分析(半翅目: 蝽科)(英文)

【目的】本研究对红角辉蝽Carbula crassiventris和紫翅果蝽Carpocoris purpureipennis完整线粒体基因组测序,以探究蝽亚科(Pentatominae)线粒体基因组特征并重建其系统发育关系。【方法】使用Illumina MiSeq测序平台测定红角辉蝽和紫翅果蝽线粒体基因组全序列,并进行组装和注释。基于这2个种和其他30个蝽亚科分类单元线粒体基因组的13个蛋白质编码基因的第1和2位密码子以及2个rRNA基因的核苷酸序列,利用贝叶斯和最大似然法重建蝽亚科系统发育树。【结果】红角辉蝽和紫翅果蝽的线粒体基因组全长分别为15 824 和16 575 bp, 包含13个蛋白质编码基因、2个rRNA基因、22个tRNA基因和1个控制区。蝽亚科内线粒体基因组基因排列顺序保守且没有发现基因重排。此外,蝽亚科内的碱基组成、密码子使用和RNA结构均较为保守; 控制区重复序列拥有不同的长度、类型和拷贝数。基于贝叶斯法和最大似然法重建的系统发育树显示二星蝽族(Eysarcorini)、果蝽族(Carpocorini)、稻绿蝽族(Nezarini)和Antestiini构成一个稳定分枝。【结论】系统发育分析支持辉蝽属Carbula应属于二星蝽族,而果蝽属Carpocoris、斑须蝽属Dolycoris和珠蝽属Rubiconia同属于果蝽族。     

Evolutionary relationships of the cup-fungus genus Peziza and Pezizaceae inferred from multiple nuclear genes: RPB2, beta-tubulin, and LSU rDNA

To provide a robust phylogeny of Pezizaceae, partial sequences from two nuclear protein-coding genes, RPB2 (encoding the second largest subunit of RNA polymerase II) and beta-tubulin, were obtained from 69 and 72 specimens, respectively, to analyze with nuclear ribosomal large subunit RNA gene sequences (LSU). The three-gene data set includes 32 species of Peziza, and 27 species from nine additional epigeous and six hypogeous (truffle) pezizaceous genera. Analyses of the combined LSU, RPB2, and beta-tubulin data set using parsimony, maximum likelihood, and Bayesian approaches identify 14 fine-scale lineages within Pezizaceae. Species of Peziza occur in eight of the lineages, spread among other genera of the family, confirming the non-monophyly of the genus. Although parsimony analyses of the three-gene data set produced a nearly completely resolved strict consensus tree, with increased confidence, relationships between the lineages are still resolved with mostly weak bootstrap support. Bayesian analyses of the three-gene data, however, show support for several more inclusive clades, mostly congruent with Bayesian analyses of RPB2. No strongly supported incongruence was found among phylogenies derived from the separate LSU, RPB2, and beta-tubulin data sets. The RPB2 region appeared to be the most informative single gene region based on resolution and clade support, and accounts for the greatest number of potentially parsimony informative characters within the combined data set, followed by the LSU and the beta-tubulin region. The results indicate that third codon positions in beta-tubulin are saturated, especially for sites that provide information about the deeper relationships. Nevertheless, almost all phylogenetic signal in beta-tubulin is due to third positions changes, with almost no signal in first and second codons, and contribute phylogenetic information at the "fine-scale" level within the Pezizaceae. The Pezizaceae is supported as monophyletic in analyses of the three-gene data set, but its sister-group relationships is not resolved with support. The results advocate the use of RPB2 as a marker for ascomycete phylogenetics at the inter-generic level, whereas the beta-tubulin gene appears less useful.     

Horizontal gene transfer of a plastid gene in the non-photosynthetic flowering plants Orobanche and Phelipanche (Orobanchaceae)

Plastid sequences are among the most widely used in phylogenetic and phylogeographic studies in flowering plants, where they are usually assumed to evolve like non-recombining, uniparentally transmitted, single-copy genes. Among others, this assumption can be violated by intracellular gene transfer (IGT) within cells or by the exchange of genes across mating barriers (horizontal gene transfer, HGT). We report on HGT of a plastid region including rps2, trnL-F, and rbcL in a group of non-photosynthetic flowering plants. Species of the parasitic broomrape genus Phelipanche harbor two copies of rps2, a plastid ribosomal gene, one corresponding to the phylogenetic position of the respective species, the other being horizontally acquired from the related broomrape genus Orobanche. While the vertically transmitted copies probably reside within the plastid genome, the localization of the horizontally acquired copies is not known. With both donor and recipient being parasitic plants, a possible pathway for the exchange of genetic material is via a commonly attacked host.     

Complete plastid genome sequence of the chickpea (Cicer arietinum) and the phylogenetic distribution of rps12 and clpP intron losses among legumes (Leguminosae)

Chickpea (Cicerarietinum, Leguminosae), an important grain legume, is widely used for food and fodder throughout the world. We sequenced the complete plastid genome of chickpea, which is 125,319bp in size, and contains only one copy of the inverted repeat (IR). The genome encodes 108 genes, including 4 rRNAs, 29 tRNAs, and 75 proteins. The genes rps16, infA, and ycf4 are absent in the chickpea plastid genome, and ndhB has an internal stop codon in the 5'exon, similar to other legumes. Two genes have lost their introns, one in the 3'exon of the transpliced gene rps12, and the one between exons 1 and 2 of clpP; this represents the first documented case of the loss of introns from both of these genes in the same plastid genome. An extensive phylogenetic survey of these intron losses was performed on 302 taxa across legumes and the related family Polygalaceae. The clpP intron has been lost exclusively in taxa from the temperate "IR-lacking clade" (IRLC), whereas the rps12 intron has been lost in most members of the IRLC (with the exception of Wisteria, Callerya, Afgekia, and certain species of Millettia, which represent the earliest diverging lineages of this clade), and in the tribe Desmodieae, which is closely related to the tribes Phaseoleae and Psoraleeae. Data provided here suggest that the loss of the rps12 intron occurred after the loss of the IR. The two new genomic changes identified in the present study provide additional support of the monophyly of the IR-loss clade, and resolution of the pattern of the earliest-branching lineages in this clade. The availability of the complete chickpea plastid genome sequence also provides valuable information on intergenic spacer regions among legumes and endogenous regulatory sequences for plastid genetic engineering.     

Comparative evolution of the mitochondrial cytochrome b gene and nuclear beta-fibrinogen intron 7 in woodpeckers

Most molecular phylogenetic studies of vertebrates have been based on DNA sequences of mitochondrial-encoded genes. MtDNA evolves rapidly and is thus particularly useful for resolving relationships among recently evolved groups. However, it has the disadvantage that all of the mitochondrial genes are inherited as a single linkage group so that only one independent gene tree can be inferred regardless of the number of genes sequenced. Introns of nuclear genes are attractive candidates for independent sources of rapidly evolving DNA: they are pervasive, most of their nucleotides appear to be unconstrained by selection, and PCR primers can be designed for sequences in adjacent exons where nucleotide sequences are conserved. We sequenced intron 7 of the beta-fibrinogen gene (beta-fibint7) for a diversity of woodpeckers and compared the phylogenetic signal and nucleotide substitution properties of this DNA sequence with that of mitochondrial-encoded cytochrome b (cyt b) from a previous study. A few indels (insertions and deletions) were found in the beta-fibint7 sequences, but alignment was not difficult, and the indels were phylogentically informative. The beta-fibint7 and cyt b gene trees were nearly identical to each other but differed in significant ways from the traditional woodpecker classification. Cyt b evolves 2.8 times as fast as beta-fibint7 (14. 0 times as fast at third codon positions). Despite its relatively slow substitution rate, the phylogenetic signal in beta-fibint7 is comparable to that in cyt b for woodpeckers, because beta-fibint7 has less base composition bias and more uniform nucleotide substitution probabilities. As a consequence, compared with cyt b, beta-fibint7 nucleotide sites are expected to enter more distinct character states over the course of evolution and have fewer multiple substitutions and lower levels of homoplasy. Moreover, in contrast to cyt b, in which nearly two thirds of nucleotide sites rarely vary among closely related taxa, virtually all beta-fibint7 nucleotide sites appear free of selective constraints, which increases informative sites per unit sequenced. However, the estimated gamma distribution used to model rate variation among sites suggests constraints on some beta-fibint7 sites. This study suggests that introns will be useful for phylogenetic studies of recently evolved groups.     

Complete Sequence of the Duckweed (<Emphasis Type="Italic">Lemna minor</Emphasis>) Chloroplast Genome: Structural Organization and Phylogenetic Relationships to Other Angiosperms

The complete nucleotide sequence of the duckweed (Lemna minor) chloroplast genome (cpDNA) was determined. The cpDNA is a circular molecule of 165,955 bp containing a pair of 31,223-bp inverted repeat regions (IRs), which are separated by small and large single-copy regions of 89,906 and 13,603 bp, respectively. The entire gene pool and relative positions of 112 genes (78 protein-encoding genes, 30 tRNA genes, and 4 rRNA genes) are almost identical to those of Amborella trichopoda cpDNA; the minor difference is the absence of infA and ycf15 genes in the duckweed cpDNA. The inverted repeat is expanded to include ycf1 and rps15 genes; this pattern is unique and does not occur in any other sequenced cpDNA of land plants. As in basal angiosperms and eudicots, but not in other monocots, the borders between IRs and a large single-copy region are located upstream of rps19 and downstream of trnH, so that trnH is not included in IRs. The model of rearrangements of the chloroplast genome during the evolution of monocots is proposed as the result of the comparison of cpDNA structures in duckweed and other monocots. The phylogenetic analyses of 61 protein-coding genes from 38 plastid genome sequences provided strong support for the monophyly of monocots and position of Lemna as the next diverging lineage of monocots after Acorales. Our analyses also provided support for Amborella as a sister to all other angiosperms, but in the bayesian phylogeny inference based on the first two codon positions Amborella united with Nymphaeales.     

The phylogeny of the Asteridae sensu lato based on chloroplast ndhF gene sequences

A phylogenetic study of Asteridae sensu lato was conducted based on chloroplast ndhF gene sequences for 116 ingroup and 13 outgroup species. Prior molecular studies based on rbcL sequences identified terminal groups corresponding to families, but were unable to resolve relationships among them. These results are largely consistent with earlier rbcL studies, but provide much greater resolution and stronger bootstrap support throughout the tree. The parsimony analysis found eight equally parsimonious trees, all of which recognize four major clades with the following relationship: (Cornales (Ericales (Euasterids I, Euasterids II))). Euasterids I includes (Garryales ((Solanales, Boraginaceae) (Gentianales, Lamiales))), although with weak support for relationships among the named clades. Euasterids II includes (Aquifoliales (Asterales (Apiales, Dipsacales))) with strong support for these relationships. Relationships within Ericales are weakly supported and merit further attention.     

A supermatrix approach provides a comprehensive genus-level phylogeny for Gentianales

Gentianales consist of Apocynaceae, Gelsemiaceae, Gentianaceae, Loganiaceae, and Rubiaceae, of which the majority are woody plants in tropical and subtropical areas. Despite extensive efforts in reconstructing the phylogeny of Gentianales based on molecular data, some interfamily and intrafamily relationships remain uncertain. We reconstructed the genus-level phylogeny of Gentianales based on the supermatrix of eight plastid markers (rbcL, matK, atpB, ndhF, rpl16, rps16, thetrnL-trnF region, and atpB-rbcL spacer) and one mitochondrial gene (matR) using maximum likelihood. The major clades and their relationships retrieved in the present study concur with those of previous studies. All of the five families of Gentianales are monophyletic with strong support. We resolved Rubiaceae as sister to the remaining families in Gentianales and showed support for the sister relationship between Loganiaceae and Apocynaceae. Our results provide new insights into relationships among intrafamilial clades. For example, within Rubiaceae we found that Craterispermeae were sister to Morindeae + (Palicoureeae + Psychotrieae) and that Theligoneae were sister to Putorieae. Within Gentianaceae, our phylogeny revealed that Gentianeae were sister to Helieae and Potalieae, and subtribe Lisianthiinae were sister to Potaliinae and Faroinae. Within Loganiaceae, we found Neuburgia as sister to Spigelieae. Within Apocynaceae, our results supported Amsonieae as sister to Melodineae, and Hunterieae as sister to a clade comprising Plumerieae + (Carisseae + APSA). We also confirmed the monophyly of Perplocoideae and the relationships among Baisseeae + (Secamonoideae + Asclepiadoideae).     

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.