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.     

Discovery of paralogous nuclear gene sequences coding for the second-largest subunit of RNA polymerase II (RPB2) and their phylogenetic utility in gentianales of the asterids

Paralogous sequences of the RPB2 gene are demonstrated in the angiosperm order Gentianales. Two different copies were found by using different PCR primer pairs targeting a region that corresponds to exons 22-24 in the Arabidopsis RPB2 gene. One of the copies (RPB2-d) lacks introns in this region, whereas the other has introns at locations corresponding to those of green plants previously investigated. When analyzed with other available RPB2 sequences from this region, all 28 RPB2-d sequences obtained from the Gentianales and the four sequences from the Lamiales form a monophyletic group, together with a previously published tomato cDNA sequence. The substitution patterns, relative rates of change, and nucleotide compositions of the two paralogous RPB2 exon regions are similar, and none of them shows any signs of being a pseudogene. Although multiple copies of similar, paralogous sequences can confound phylogenetic interpretations, the lack of introns in RPB2-d make a priori homology assessment easy. The phylogenetic utility of RPB2-d within the Gentianales is evaluated in comparison with the chloroplast genes ndhF and rbcL. The hierarchical information in the RPB2-d region sequenced is more incongruent with that of the plastid genes than the plastid genes are with each other as determined by incongruence length difference tests. In contrast to the plastid genes, parsimony-informative third codon positions of RPB2 have a significantly higher rate of change than first and second positions. Topologically, the trees from the three genes are similar, and the differences are usually only weakly supported. In terms of support, RPB2 gives the highest jackknife support per sequenced nucleotide, whereas ndhF gives the highest Bremer support per sequenced nucleotide. The RPB2-d locus has the potential to be a valuable nuclear marker for determination of phylogenetic relationships within the euasterid I group of plants.     

Angiosperm phylogeny inferred from sequences of four mitochondrial genes

An angiosperm phylogeny was reconstructed in a maximum likelihood analysis of sequences of four mitochondrial genes, atpl, matR, had5, and rps3, from 380 species that represent 376 genera and 296 families of seed plants. It is largely congruent with the phylogeny of angiosperms reconstructed from chloroplast genes atpB, matK, and rbcL, and nuclear 18S rDNA. The basalmost lineage consists of Amborella and Nymphaeales (including Hydatellaceae). Austrobaileyales follow this clade and are sister to the mesangiosperms, which include Chloranthaceae, Ceratophyllum, magnoliids, monocots, and eudicots. With the exception of Chloranthaceae being sister to Ceratophyllum, relationships among these five lineages are not well supported. In eudicots, Ranunculales, Sabiales, Proteales, Trochodendrales, Buxales, Gunnerales, Saxifragales, Vitales, Berberidopsidales, and Dilleniales form a basal grade of lines that diverged before the diversification of rosids and asterids. Within rosids, the COM (Celastrales-Oxalidales-Malpighiales) clade is sister to malvids (or rosid Ⅱ), instead of to the nitrogen-fixing clade as found in all previous large-scale molecular analyses of angiosperms. Santalales and Caryophyllales are members of an expanded asterid clade. This study shows that the mitochondrial genes are informative markers for resolving relationships among genera, families, or higher rank taxa across angiosperms. The low substitution rates and low homoplasy levels of the mitochondrial genes relative to the chloroplast genes, as found in this study, make them particularly useful for reconstructing ancient phylogenetic relationships. A mitochondrial gene-based angiosperm phylogeny provides an independent and essential reference for comparison with hypotheses of angiosperm phylogeny based on chloroplast genes, nuclear genes, and non-molecular data to reconstruct the underlying organismal phylogeny.     

The origin and diversification of angiosperms

The angiosperms, one of five groups of extant seed plants, are the largest group of land plants. Despite their relatively recent origin, this clade is extremely diverse morphologically and ecologically. However, angiosperms are clearly united by several synapomorphies. During the past 10 years, higher-level relationships of the angiosperms have been resolved. For example, most analyses are consistent in identifying Amborella, Nymphaeaceae, and Austrobaileyales as the basalmost branches of the angiosperm tree. Other basal lineages include Chloranthaceae, magnoliids, and monocots. Approximately three quarters of all angiosperm species belong to the eudicot clade, which is strongly supported by molecular data but united morphologically by a single synapomorphy-triaperturate pollen. Major clades of eudicots include Ranunculales, which are sister to all other eudicots, and a clade of core eudicots, the largest members of which are Saxifragales, Caryophyllales, rosids, and asterids. Despite rapid progress in resolving angiosperm relationships, several significant problems remain: (1) relationships among the monocots, Chloranthaceae, magnoliids, and eudicots, (2) branching order among basal eudicots, (3) relationships among the major clades of core eudicots, (4) relationships within rosids, (5) relationships of the many lineages of parasitic plants, and (6) integration of fossils with extant taxa into a comprehensive tree of angiosperm phylogeny.     

Duplication and paralog sorting of RPB2 and RPB1 genes in core eudicots

RPB1 and RPB2, which encode the largest and second largest subunits of RNA polymerase II, respectively, are essential single copy genes in fungi, animals and most plants. Two paralogs of the RPB2 gene have been found in some groups of angioperms [Oxelman, B., Yoshikawa, N., McConaughy, B.L., Luo, J., Denton, A.L., Hall, B.D., 2004. RPB2 gene phylogeny in flowering plants, with particular emphasis on asterids. Mol. Phylogenet. Evol. 32, 462-479]. Here, we report the results of experiments designed to identify the evolutionary origin of the RPB2 duplicate copies. Through careful sampling and phylogenetic analysis, we were able to construct the RPB2 gene tree in angiosperms and infer the phylogenetic positions of the gene duplication and gene loss events that occurred. Our study shows that an RPB2 gene duplication occurred early in core eudicot evolution, at or near the time of the Buxaceae/Trochodendraceae divergence. Subsequently, multiple gene duplication and paralog sorting events happened independently in different core eudicot taxa. Differential expression of the two RPB2 gene paralogs may explain the preservation of both paralogs in the asterids. One gene (RPB2-i) accounts for most of the RPB2 mRNA made in the flower organs while the other gene (RPB2-d) is predominantly used in the vegetative tissues. We also found two paralogs of the RPB1 gene in some core eudicot species. The RPB1 gene duplication occurred before core eudicot divergence, around the time of RPB2 gene duplication. Several independent RPB1 paralog sorting events happened in different core eudicot taxa; their occurrence was independent of the RPB2 paralog sorting events. Our results suggest that a polyploidization event happened at or near the time of the Buxaceae/Trochodendraceae divergence. We propose that this polyploidization and the partial diploidization processes thereafter may have been the driving force of core eudicot radiation.     

Evolution of a RNA polymerase gene family in Silene (Caryophyllaceae)-incomplete concerted evolution and topological congruence among paralogues

Four low-copy nuclear DNA intron regions from the second largest subunits of the RNA polymerase gene family (RPA2, RPB2, RPD2a, and RPD2b), the internal transcribed spacers (ITSs) from the nuclear ribosomal regions, and the rps16 intron from the chloroplast were sequenced and used in a phylogenetic analysis of 29 species from the tribe Sileneae (Caryophyllaceae). We used a low stringency nested polymerase chain reaction (PCR) approach to overcome the difficulties of constructing specific primers for amplification of the low copy nuclear DNA regions. Maximum parsimony analyses resulted in largely congruent phylogenetic trees for all regions. We tested overall model congruence in a likelihood context using the software PLATO and found that ITSs, RPA2, and RPB2 deviated from the maximum likelihood model for the combined data. The topology parameter was then isolated and topological congruence assessed by nonparametric bootstrapping. No strong topological incongruence was found. The analysis of the combined data sets resolves previously poorly known major relationships within Sileneae. Two paralogues of RPD2 were found, and several independent losses and incomplete concerted evolution were inferred. The among-site rate variation was significantly lower in the RNA polymerase introns than in the rps16 intron and ITSs, a property that is attractive in phylogenetic analyses.     

Exclusive conservation of mitochondrial group II intron nad4i548 among liverworts and its use for phylogenetic studies in this ancient plant clade

Liverworts occupy a pivotal position in land plant (embryophyte) phylogeny as the presumed earliest-branching major clade, sister to all other land plants, including the mosses, hornworts, lycophytes, monilophytes and seed plants. Molecular support for this earliest dichotomy in land plant phylogeny comes from strikingly different occurrences of introns in mitochondrial genes distinguishing liverworts from all other embryophytes. Exceptionally, however, the nad5 gene--the mitochondrial locus hitherto used most widely to elucidate early land plant phylogeny--carries a group I type intron that is shared between liverworts and mosses. We here explored whether a group II intron, the other major type of organellar intron, would similarly be conserved in position across the entire diversity of extant liverworts and could be of use for phylogenetic analyses in this supposedly most ancient embryophyte clade. To this end, we investigated the nad4 gene as a candidate locus possibly featuring different introns in liverworts as opposed to the non-liverwort embryophyte (NLE) lineage. We indeed found group II intron nad4i548 universally conserved in a wide phylogenetic sampling of 55 liverwort taxa, confirming clade specificity and surprising evolutionary stability of plant mitochondrial introns. As expected, intron nad4i548g2 carries phylogenetic information in its variable sequences, which confirms and extends previous cladistic insights on liverwort evolution. We integrate the new nad4 data with those of the previously established mitochondrial nad5 and the chloroplast rbcL and rps4 genes and present a phylogeny based on the fused datasets. Notably, the phylogenetic analyses suggest a reconsideration of previous phylogenetic and taxonomic assignments for the genera Calycularia and Mylia and resolve a sister group relationship of Ptilidiales and Porellales.     

Use of flowering gene <Emphasis Type="Italic">FLOWERING LOCUS T</Emphasis> (<Emphasis Type="Italic">FT</Emphasis>) homologs in the phylogenetic analysis of bambusoid and early diverging grasses

A nuclear gene, FLOWERING LOCUS T (FT) homolog, was cloned from Phyllostachys meyeri as PmFT. Its putative copy number was estimated as four by Southern blot analysis, and the two copies were completely sequenced. Twenty-seven FT homolog sequences of bambusoid and early diverging grasses comprised 172-bp exons, and 357- to 785-bp introns exhibited 0-58.9% pairwise divergence with six modal levels. Parsimony analyses of the FT homologs rooted at Pharus virescens produced six equally parsimonious trees. In the strict consensus tree, five clades were resolved; they were affected by divergence of the intron region rather than exon region. The basal clade was Puelioideae, followed by Olyreae clade including Oryza sativa. Streptogyneae clade combined the Olyreae clade with terminal sister clades of the Bambuseae, i.e., pantropical bamboos and East Asiatic temperate bamboos. The global topology suggested that FT homologs are significant for resolving the tribe level. However, the phylogeny of FT homologs does not resolve monophyly in Bambusoideae because of intercalary positioning by Streptogyneae clade. We discussed the role of FT homologs in controlling the inflorescence architecture and position of Streptogyneae in the bamboo phylogeny.     

Molecular phylogeny of Zygomycota based on EF-1alpha and RPB1 sequences: limitations and utility of alternative markers to rDNA

Earlier molecular phylogenetic analyses based on nuclear small subunit ribosomal DNA (nSSU rDNA) suggest that the Zygomycota are polyphyletic within the Chytridiomycota. However, these analyses failed to resolve almost all interordinal relationships among basal fungi (Chytridiomycota and Zygomycota), due to lack of sufficient characters within the nSSU rDNA. To further elucidate the higher-level phylogeny of Zygomycota, we have sequenced partial RPB1 (DNA dependent RNA polymerase II largest subunit) and EF-1alpha (translation elongation factor 1 alpha) genes from 10 and 3 zygomycete fungi, respectively. Independent molecular phylogenetic analyses were performed based on each sequence by distance and maximum likelihood methods. Although deep phylogenetic relationships among basal fungi still remain poorly resolved using either gene, the RPB1-based phylogeny identified a novel monophyletic clade consisting of the Dimargaritales, Harpellales, and Kickxellales. This result suggests that regularly formed septa (cross walls that divide hyphae into segments) with a lenticular cavity are plesiomorphic for this clade, and indicates the importance of septal pore ultrastructure in zygomycete phylogeny. In addition, a peculiar mucoralean genus Mortierella, which was considered to be distantly related to the other Mucorales based on previous nSSU rDNA analyses, was resolved as the basal most divergence within the Mucorales, consistent with traditional phenotypic-based taxonomy. Although the taxa included in our analysis are restricted, the monophyly of each order suggested by nSSU rDNA phylogeny is supported by the present RPB1-based analysis. These results support the potential use of RPB1 as an alternative marker for fungal phylogenetic studies. Conversely, the overall fungal phylogeny based on EF-1alpha sequence is poorly resolved. A comparison of numbers of observed substitutions versus inferred substitutions within EF-1alpha indicates that this gene is much more saturated than RPB1. This result suggests that the EF-1alpha gene is unsuitable for resolving higher-level phylogenetic relationships within the Fungi.     

木根麦冬与林生麦冬5S rRNA基因的进化

为了揭示多拷贝基因的进化方式, 对濒危植物木根麦冬(Ophiopogon xylorrhizus Wang et Dai) 3个居群、6个个体的294个5S rR NA 基因拷贝及其姐妹种林生麦冬(O. sylvicola Wang et Tang)的45个拷贝进行了DNA测序和序列分析,并以这个迄今发表的最大的单个物种的5S rRNA基因数据,以PAUP程序重建了分子系统发育树.结果表明: 1)所得序列呈高度多样性,长度变化在307～548碱基之间,仅13对(3.8%)相同,序列分化指数较高:木根麦冬是0.078,林生麦冬是0.032,两物种间是0.149; 2)100%的统计值支持两物种的5S rRNA 基因分别来自于祖先种的一个拷贝,即"建立者拷贝",这个拷贝在物种形成之后进行了一系列连续的扩增,形成一个直系的基因家族,而祖先种的其他拷贝在物种形成后被丢失; 3)不同拷贝是独立进化的,序列间的一致化过程很弱,这在串联重复的rR N A基因中是罕见的; 4)木根麦冬居群间曾存在频繁的基因交流,使5S rRNA 基因的许多拷贝扩散于不同居群,维持着种内较高的遗传多样性.可以认为是某些近期发生的变化,阻止了居群间的基因交流,导致该物种广泛的自交,发生自交衰退,并最后导致濒危.     

Utility of 17 chloroplast genes for inferring the phylogeny of the basal angiosperms

Sequences from 14 slowly evolving chloroplast genes (including three highly conserved introns) were obtained for representative basal angiosperm and seed-plant taxa, using novel primers described here. These data were combined with published sequences from atpB, rbcL, and newly obtained sequences from ndhF. Combined data from these 17 genes permit sturdy, well-resolved inference of major aspects of basal angiosperm relationships, demonstrating that the new primers are valuable tools for sorting out the deepest events in flowering plant phylogeny. Sequences from the inverted repeat (IR) proved to be particularly reliable (low homoplasy, high retention index). Representatives of Cabomba and Illicium were the first two successive branches of the angiosperms in an initial sampling of 19 exemplar taxa. This result was strongly supported by bootstrap analysis and by two small insertion/deletion events in the slowly evolving introns. Several paleoherb groups (representatives of Piperales) formed a strongly supported clade with taxa representing core woody magnoliids (Laurales, Magnoliales, and Winteraceae). The monophyly of the sampled eudicots and monocots was also well supported. Analyses of three major partitions of the data showed many of the same clades and supported the rooting seen with all the data combined. While Amborella trichopoda was supported as the sister group of the remaining angiosperms when we added Amborella and Nymphaea odorata to the analysis, a strongly conflicting rooting was observed when Amborella alone was added.     

Angiosperm phylogeny inferred from 18S rDNA, vbcL, and atpB sequences

A phylogenetic analysis of a combined data set for 560 angiosperms and seven outgroups based on three genes, 18S rDNA (1855 bp), rbcL (1428 bp), and atpB (1450 bp) representing a total of 4733 bp is presented. Parsimony analysis was expedited by use of a new computer program, the RATCHET. Parsimony jackknifing was performed to assess the support of clades. The combination of three data sets for numerous species has resulted in the most highly resolved and strongly supported topology yet obtained for angiosperms. In contrast to previous analyses based on single genes, much of the spine of the tree and most of the larger clades receive jackknife support 250%. Some of the noneudicots form a grade followed by a strongly supported eudicot clade. The early‐branching angiosperms are Amborellaceae, Nymphaeaceae, and a clade of Austrobaileyaceae, Illiciaceae, and Schi‐sandraceae. The remaining noneudicots, except Ceratophyllaceae, form a weakly supported core eumagnoliid clade comprising six well‐supported subclades: Chloranthaceae, monocots, WinteraceaeICanellaceae, Piperales, Laurales, and Magnoliales. Ceratophyllaceae are sister to the eudicots. Within the well‐supported eudicot clade, the early‐diverging eudicots (e.g. Proteales, Ranunculales, Trochodendraceae, Sabiaceae) form a grade, followed by the core eudicots, the monophyly of which is also strongly supported. The core eudicots comprise six well‐supported subclades: (1) Berberidopsidaceae/Aextoxicaceae; (2) Myrothamnaceae/ Gunneraceae; (3) Saxifragales, which are the sister to Vitaceae (including Leea) plus a strongly supported eurosid clade; (4) Santalales; (5) Caryophyllales, to which Dilleniaceae are sister; and (6) an asterid clade. The relationships among these six subclades of core eudicots do not receive strong support. This large data set has also helped place a number of enigmatic angiosperm families, including Podostemaceae, Aphloiaceae, and Ixerbaceae. This analysis further illustrates the tractability of large data sets and supports a recent, phylogenetically based, ordinal‐level reclassification of the angiosperms based largely, but not exclusively, on molecular (DNA sequence) data.     

A survey of introns in three genes of rotifers

Here we report on the occurrence and position of introns found in three genes of rotifers. A region of the gene for the TATA-box binding protein was examined in three species of Bdelloidea and one of Monogononta. There are two introns in both copies of this gene present in each of the three bdelloids examined – one at a position where introns occur in other eukaryotes and the other at a novel position; the monogonont has no introns in the region examined. A region of the gene encoding the 82 kD heat shock protein was examined in 10 species, with every rotifer class represented. Introns were found in only two species, both bdelloids: one of the species has an intron in all three copies of the gene; the other has an intron in only one of the three copies. Both introns occur at novel positions. The gene for triosephosphate isomerase was examined in one bdelloid. Both copies of the gene in this species contain introns, all at conserved positions: one copy contains five introns, the other copy three. These observations demonstrate the presence of introns in bdelloid rotifers, some in conserved positions, others apparently newly arisen during bdelloid evolution.     

Phylogenetic relationships among arecoid palms (Arecaceae: Arecoideae)

Background and Aims

The Arecoideae is the largest and most diverse of the five subfamilies of palms (Arecaceae/Palmae), containing >50 % of the species in the family. Despite its importance, phylogenetic relationships among Arecoideae are poorly understood. Here the most densely sampled phylogenetic analysis of Arecoideae available to date is presented. The results are used to test the current classification of the subfamily and to identify priority areas for future research.

Methods

DNA sequence data for the low-copy nuclear genes PRK and RPB2 were collected from 190 palm species, covering 103 (96 %) genera of Arecoideae. The data were analysed using the parsimony ratchet, maximum likelihood, and both likelihood and parsimony bootstrapping.

Key Results and Conclusions

Despite the recovery of paralogues and pseudogenes in a small number of taxa, PRK and RPB2 were both highly informative, producing well-resolved phylogenetic trees with many nodes well supported by bootstrap analyses. Simultaneous analyses of the combined data sets provided additional resolution and support. Two areas of incongruence between PRK and RPB2 were strongly supported by the bootstrap relating to the placement of tribes Chamaedoreeae, Iriarteeae and Reinhardtieae; the causes of this incongruence remain uncertain. The current classification within Arecoideae was strongly supported by the present data. Of the 14 tribes and 14 sub-tribes in the classification, only five sub-tribes from tribe Areceae (Basseliniinae, Linospadicinae, Oncospermatinae, Rhopalostylidinae and Verschaffeltiinae) failed to receive support. Three major higher level clades were strongly supported: (1) the RRC clade (Roystoneeae, Reinhardtieae and Cocoseae), (2) the POS clade (Podococceae, Oranieae and Sclerospermeae) and (3) the core arecoid clade (Areceae, Euterpeae, Geonomateae, Leopoldinieae, Manicarieae and Pelagodoxeae). However, new data sources are required to elucidate ambiguities that remain in phylogenetic relationships among and within the major groups of Arecoideae, as well as within the Areceae, the largest tribe in the palm family.     

Phylogeny reconstruction in the Caesalpinieae grade (Leguminosae) based on duplicated copies of the sucrose synthase gene and plastid markers

The Caesalpinieae grade (Leguminosae) forms a morphologically and ecologically diverse group of mostly tropical tree species with a complex evolutionary history. This grade comprises several distinct lineages, but the exact delimitation of the group relative to subfamily Mimosoideae and other members of subfamily Caesalpinioideae, as well as phylogenetic relationships among the lineages are uncertain. With the aim of better resolving phylogenetic relationships within the Caesalpinieae grade, we investigated the utility of several nuclear markers developed from genomic studies in the Papilionoideae. We cloned and sequenced the low copy nuclear gene sucrose synthase (SUSY) and combined the data with plastid trnL and matK sequences. SUSY has two paralogs in the Caesalpinieae grade and in the Mimosoideae, but occurs as a single copy in all other legumes tested. Bayesian and maximum likelihood phylogenetic analyses suggest the two nuclear markers are congruent with plastid DNA data. The Caesalpinieae grade is divided into four well-supported clades (Cassia, Caesalpinia, Tachigali and Peltophorum clades), a poorly supported clade of Dimorphandra Group genera, and two paraphyletic groups, one with other Dimorphandra Group genera and the other comprising genera previously recognized as the Umtiza clade. A selection analysis of the paralogs, using selection models from PAML, suggests that SUSY genes are subjected to a purifying selection. One of the SUSY paralogs, under slightly stronger positive selection, may be undergoing subfunctionalization. The low copy SUSY gene is useful for phylogeny reconstruction in the Caesalpinieae despite the presence of duplicate copies. This study confirms that the Caesalpinieae grade is an artificial group, and highlights the need for further analyses of lineages at the base of the Mimosoideae.     

Using RPB1 sequences to improve phylogenetic inference among mushrooms (Inocybe, Agaricales)

An investigation of mushroom phylogeny using the largest subunit of RNA polymerase II gene sequences (RPB1) was conducted in comparison with nuclear ribosomal large subunit RNA gene sequences (nLSU) for the same set of taxa in the genus Inocybe (Agaricales, Basidiomycota). The two data sets, though not significantly incongruent, exhibit conflict among the placement of two taxa that exhibit long branches in the nLSU data set. In contrast, RPB1 terminal branch lengths are rather uniform. Bootstrap support is increased for clades in RPB1. Combined data sets increase the degree of confidence for several relationships. Overall, nLSU data do not yield a robust phylogeny when independently assessed by RPB1 sequences. This multigene study indicates that Inocybe is a monophyletic group composed of at least four distinct lineages-subgenus Mallocybe, section Cervicolores, section Rimosae, and subgenus Inocybe sensu Kühner, Kuyper, non Singer. Within subgenus Inocybe, two additional lineages, one composed of species with smooth basidiospores (clade I) and a second characterized by nodulose-spored species (clade II), are recovered by RPB1 and combined data. The nLSU data recover only clade I. The genera Astrosporina and Inocybella cannot be recognized phylogenetically. "Supersections" Cortinatae and Marginatae are not monophyletic groups.     

Phylogenetic relationships of the dwarf boas and a comparison of Bayesian and bootstrap measures of phylogenetic support

Four New World genera of dwarf boas (Exiliboa, Trachyboa, Tropidophis, and Ungaliophis) have been placed by many systematists in a single group (traditionally called Tropidophiidae). However, the monophyly of this group has been questioned in several studies. Moreover, the overall relationships among basal snake lineages, including the placement of the dwarf boas, are poorly understood. We obtained mtDNA sequence data for 12S, 16S, and intervening tRNA-val genes from 23 species of snakes representing most major snake lineages, including all four genera of New World dwarf boas. We then examined the phylogenetic position of these species by estimating the phylogeny of the basal snakes. Our phylogenetic analysis suggests that New World dwarf boas are not monophyletic. Instead, we find Exiliboa and Ungaliophis to be most closely related to sand boas (Erycinae), boas (Boinae), and advanced snakes (Caenophidea), whereas Tropidophis and Trachyboa form an independent clade that separated relatively early in snake radiation. Our estimate of snake phylogeny differs significantly in other ways from some previous estimates of snake phylogeny. For instance, pythons do not cluster with boas and sand boas, but instead show a strong relationship with Loxocemus and Xenopeltis. Additionally, uropeltids cluster strongly with Cylindrophis, and together are embedded in what has previously been considered the macrostomatan radiation. These relationships are supported by both bootstrapping (parametric and nonparametric approaches) and Bayesian analysis, although Bayesian support values are consistently higher than those obtained from nonparametric bootstrapping. Simulations show that Bayesian support values represent much better estimates of phylogenetic accuracy than do nonparametric bootstrap support values, at least under the conditions of our study.     

木根麦冬与林生麦冬5S rRNA基因的进化

为了揭示多拷贝基因的进化方式,对濒危植物木根麦冬（Ophiopogon xylorrhizus Wang et Dai)3个居群、6个个体的294个5S rRNA基因拷贝及其姐妹种林生麦冬（O.sylvicola Wang et Tang) 的45个拷贝进行了DNA测序和序列分析,并以这个迄今发表的最大的单个物种的5S rDNA基因数据,以PAUP程序重建了分子系统发育树。结果表明：1）所得序列呈高度多样性,长度变化在307－548碱基之间,仅13对（3．8％）相同,序列分化指数较高：木根麦冬是0．078,林生麦冬是0．032,两物种间是0．149;2）100％的统计值支持两物种的5S rRNA基因分别来自于祖先种的一个拷贝,即“建立者拷贝”,这个拷贝在物种形成之后进行了一系列连续的扩增,形成一个直系的基因家族,而祖先种的其他拷贝在物种形成后被丢失;3）不同拷贝是独立进化的,序列间的一致化过程很弱,这在串联重复的rRNA基因中是罕见的;4）木根麦冬居群间曾存在频繁的基因交流,使5S rRNA基因的许多拷贝扩散于不同居群,维持着种内较高的遗传多样性。可以认为是某些近期发生的变化,阻止了居群间的基因交流,导致该物种广泛的自交,发生自交衰退,并最后导致濒危。     

Intron evolution in a phylogenetic perspective: Divergent trends in the two copies of the duplicated def gene in Impatiens L. (Balsaminaceae)

The history of MADS box genes is well-known in angiosperms. While duplication events and gene losses occur frequently, gene structure and intron positions are very conserved. We investigated all six introns in a duplicated MADS box gene (deficiens, def) in selected Impatiens taxa, thereby assessing intron features. For the first time, our study provides a comparison of molecular changes in all introns of a gene from a phylogenetic perspective. Interestingly, a uniform pattern of molecular evolution in the introns of each copy was not observed, but intron length increases, decreases, and size retention can be found in each copy. A tendency to accumulate long autapomorphic indels is also present, thus, a longer intron length does not reflect a higher number of parsimony-informative characters. Substitution rates vary between introns of each gene copy. While four of the six introns of def1 exhibit a change in their substitution rate, five of the six def2 introns maintain their rates throughout the genus albeit at different levels. In MADS box genes several regulatory sequences are found residing in introns. Thus, presence of putative regulatory motifs was investigated. Most of them are not conserved in position and usually present in only one of the gene copies. In addition, the potential for phylogenetic reconstruction of introns in both def copies is shortly discussed.     

Complete sequence of Euglena gracilis chloroplast DNA.

We report the complete DNA sequence of the Euglena gracilis, Pringsheim strain Z chloroplast genome. This circular DNA is 143,170 bp, counting only one copy of a 54 bp tandem repeat sequence that is present in variable copy number within a single culture. The overall organization of the genome involves a tandem array of three complete and one partial ribosomal RNA operons, and a large single copy region. There are genes for the 16S, 5S, and 23S rRNAs of the 70S chloroplast ribosomes, 27 different tRNA species, 21 ribosomal proteins plus the gene for elongation factor EF-Tu, three RNA polymerase subunits, and 27 known photosynthesis-related polypeptides. Several putative genes of unknown function have also been identified, including five within large introns, and five with amino acid sequence similarity to genes in other organisms. This genome contains at least 149 introns. There are 72 individual group II introns, 46 individual group III introns, 10 group II introns and 18 group III introns that are components of twintrons (introns-within-introns), and three additional introns suspected to be twintrons composed of multiple group II and/or group III introns, but not yet characterized. At least 54,804 bp, or 38.3% of the total DNA content is represented by introns.