Topological incongruence between nuclear and chloroplast DNA trees suggesting hybridization in the urophyllum group of the genus Fagopyrum (Polygonaceae)

We performed phylogenetic analyses of Fagopyrum species in the urophyllum group based on nucleotide sequences of two nuclear genes, FLORICAULA/LEAFY (FLO/LFY) and AGAMOUS (AG), and three segments of chloroplast DNA (cpDNA), rbcL-accD, trnK intron, and trnC-rpoB spacer. The FLO/LFY and AG sequences turned out to be phylogenetically more informative at the intrageneric level than the cpDNA sequences. Congruence among these gene trees, inferred by a maximum-likelihood (ML) method, demonstrated that topologies were partially incongruent between the nuclear and chloroplast DNA phylogenies. The nuclear DNA sequence data supported a monophyletic relation of F. statice, F. gilesii, and F. jinshaense, whereas the former two species formed another monophyletic relation with the F. capillatum-F. gracilipes-F. gracilipedoides-F. rubifolium clade excluding F. jinshaense in the synthetic cpDNA phylogeny. In addition, two divergent sequences of FLO/LFY were found in F. rubifolium (tetraploid). One of these was sister to F. gracilipedoides and another was sister to F. statice, and a monophyletic relation of these two genes was rejected by a bootstrap analysis. These results suggest that hybridization may have occurred during diversification of Fagopyrum species in the urophyllum group, and that F. rubifolium is possibly allotetraploid species.     

Regulatory elements of the floral homeotic gene AGAMOUS identified by phylogenetic footprinting and shadowing

In Arabidopsis thaliana, cis-regulatory sequences of the floral homeotic gene AGAMOUS (AG) are located in the second intron. This 3-kb intron contains binding sites for two direct activators of AG, LEAFY (LFY) and WUSCHEL (WUS), along with other putative regulatory elements. We have used phylogenetic footprinting and the related technique of phylogenetic shadowing to identify putative cis-regulatory elements in this intron. Among 29 Brassicaceae species, several other motifs, but not the LFY and WUS binding sites identified previously, are largely invariant. Using reporter gene analyses, we tested six of these motifs and found that they are all functionally important for the activity of AG regulatory sequences in A. thaliana. Although there is little obvious sequence similarity outside the Brassicaceae, the intron from cucumber AG has at least partial activity in A. thaliana. Our studies underscore the value of the comparative approach as a tool that complements gene-by-gene promoter dissection but also demonstrate that sequence-based studies alone are insufficient for a complete identification of cis-regulatory sites.     

Molecular phylogeny of Trichomonadidae family inferred from ITS-1, 5.8S rRNA and ITS-2 sequences

The Trichomonads have been the subject of several molecular studies that reported some discrepancies both at the lower and higher taxonomic levels. The purpose of this study was to make an extensive phylogenetic analysis of the Trichomonadidae using ITS-1/5.8S/ITS-2 sequences, to better understand its phylogeny and the usefulness of this marker. ITS-1/5.8S/ITS-2 sequences of 36 strains from 14 species belonging to Trichomonadidae and Monocercomonadidae were analysed, in which 20 were newly determined. Maximum likelihood, maximum parsimony, neighbour joining, and Bayesian phylogenetic methods were employed in order to reconstruct and compare the evolutionary history of this group. Tetratrichomonas gallinarum and four strains of Tetratrichomonas sp. isolated from bull genital organs were found closely related, confirming the classification of the latter, probably as a new species. The monophyly of Tritrichomonadinae and Trichomonadinae subfamilies were corroborated, with the exclusion of Trichomitus batrachorum from the latter since it grouped consistently with Hypotrichomonas acosta. Tritrichomonas foetus, Tritrichomonas suis and potentially also Tritrichomonas mobilensis seemed to correspond to the same species. Monocercomonas sp. and Ditrichomonas honigbergii emerged as independent lineages, with their phylogenetic positions undetermined. Neither Trichomonadidae nor Monocercomonadidae were supported as monophyletic groups. The ITS-1/5.8S/ITS-2 seems to be a reliable locus for phylogenetic studies in the Trichomonadida, mainly at lower taxonomic levels, and at least up to the family level.     

Intron-spanning primers for the amplification of the nuclear ANT gene in decapod crustaceans

We describe polymerase chain reaction primers that amplify the low-copy nuclear adenine nucleotide transporter gene in decapod crustaceans. These were tested on 35 species from 14 decapod families, and a single polymerase chain reaction product amplified in 32 species. Of 49 sequences generated, only two did not contain an intron, and the longest intron identified was more than 834 nucleotides in length. The amplified fragment is likely to be useful at various taxonomic levels. While the intron is suitable for phylogeographical/population genetic studies and to identify cryptic speciation, the second exon region is sufficiently long to provide signal at both the phylogeographical and phylogenetic levels.     

Tree of life for the genera of Chinese vascular plants

We reconstructed a phylogenetic tree of Chinese vascular plants (Tracheophyta) using sequences of the chloroplast genes atpB, matK, ndhF, and rbcL and mitochondrial matR. We produced a matrix comprising 6098 species and including 13?695 DNA sequences, of which 1803 were newly generated. Our taxonomic sampling spanned 3114 genera representing 323 families of Chinese vascular plants, covering more than 93% of all genera known from China. The comprehensive large phylogeny supports most relationships among and within families recognized by recent molecular phylogenetic studies for lycophytes, ferns (monilophytes), gymnosperms, and angiosperms. For angiosperms, most families in Angiosperm Phylogeny Group IV are supported as monophyletic, except for a paraphyletic Dipterocarpaceae and Santalaceae. The infrafamilial relationships of several large families and monophyly of some large genera are well supported by our dense taxonomic sampling. Our results showed that two species of Eberhardtia are sister to a clade formed by all other taxa of Sapotaceae, except Sarcosperma. We have made our phylogeny of Chinese vascular plants publically available for the creation of subtrees via SoTree (http://www.darwintree.cn/flora/index.shtml), an automated phylogeny assembly tool for ecologists.     

Compound leaf development and evolution in the legumes

Across vascular plants, Class 1 KNOTTED1-like (KNOX1) genes appear to play a critical role in the development of compound leaves. An exception to this trend is found in the Fabaceae, where pea (Pisum sativum) uses UNIFOLIATA, an ortholog of the floral regulators FLORICAULA (FLO) and LEAFY (LFY), in place of KNOX1 genes to regulate compound leaf development. To assess the phylogenetic distribution of KNOX1-independent compound leaf development, a survey of KNOX1 protein expression across the Fabaceae was undertaken. The majority of compound-leafed Fabaceae have expression of KNOX1 proteins associated with developing compound leaves. However, in a large subclade of the Fabaceae, the inverted repeat-lacking clade (IRLC), of which pea is a member, KNOX1 expression is not associated with compound leaves. These data suggest that the FLO/LFY gene may function in place of KNOX1 genes in generating compound leaves throughout the IRLC. The contribution of FLO/LFY to leaf complexity in a member of the Fabaceae outside of the IRLC was examined by reducing expression of FLO/LFY orthologs in transgenic soybean (Glycine max). Transgenic plants with reduced FLO/LFY expression showed only slight reductions in leaflet number. Overexpression of a KNOX1 gene in alfalfa (Medicago sativa), a member of the IRLC, resulted in an increase in leaflet number. This implies that KNOX1 targets, which promote compound leaf development, are present in alfalfa and are still sensitive to KNOX1 regulation. These data suggest that KNOX1 genes and the FLO/LFY gene may have played partially overlapping roles in compound leaf development in ancestral Fabaceae but that the FLO/LFY gene took over this role in the IRLC.     

Multiple data sets, high homoplasy, and the phylogeny of softshell turtles (Testudines: Trionychidae)

We present a phylogenetic hypothesis and novel, rank-free classification for all extant species of softshell turtles (Testudines:Trionychidae). Our data set included DNA sequence data from two mitochondrial protein-coding genes and a approximately 1-kb nuclear intron for 23 of 26 recognized species, and 59 previously published morphological characters for a complimentary set of 24 species. The combined data set provided complete taxonomic coverage for this globally distributed clade of turtles, with incomplete data for a few taxa. Although our taxonomic sampling is complete, most of the modern taxa are representatives of old and very divergent lineages. Thus, due to biological realities, our sampling consists of one or a few representatives of several ancient lineages across a relatively deep phylogenetic tree. Our analyses of the combined data set converge on a set of well-supported relationships, which is in accord with many aspects of traditional softshell systematics including the monophyly of the Cyclanorbinae and Trionychinae. However, our results conflict with other aspects of current taxonomy and indicate that most of the currently recognized tribes are not monophyletic. We use this strong estimate of the phylogeny of softshell turtles for two purposes: (1) as the basis for a novel rank-free classification, and (2) to retrospectively examine strategies for analyzing highly homoplasious mtDNA data in deep phylogenetic problems where increased taxon sampling is not an option. Weeded and weighted parsimony, and model-based techniques, generally improved the phylogenetic performance of highly homoplasious mtDNA sequences, but no single strategy completely mitigated the problems of associated with these highly homoplasious data. Many deep nodes in the softshell turtle phylogeny were confidently recovered only after the addition of largely nonhomoplasious data from the nuclear intron.     

Hierarchical analysis of variation in the mitochondrial 16S rRNA gene among Hymenoptera

Nucleotide sequences from a 434-bp region of the 16S rRNA gene were analyzed for 65 taxa of Hymenoptera (ants, bees, wasps, parasitoid wasps, sawflies) to examine the patterns of variation within the gene fragment and the taxonomic levels for which it shows maximum utility in phylogeny estimation. A hierarchical approach was adopted in the study through comparison of levels of sequence variation among taxa at different taxonomic levels. As previously reported for many holometabolous insects, the 16S data reported here for Hymenoptera are highly AT-rich and exhibit strong site-to-site variation in substitution rate. More precise estimates of the shape parameter (alpha) of the gamma distribution and the proportion of invariant sites were obtained in this study by employing a reference phylogeny and utilizing maximum-likelihood estimation. The effectiveness of this approach to recovering expected phylogenies of selected hymenopteran taxa has been tested against the use of maximum parsimony. This study finds that the 16S gene is most informative for phylogenetic analysis at two different levels: among closely related species or populations, and among tribes, subfamilies, and families. Maximization of the phylogenetic signal extracted from the 16S gene at higher taxonomic levels may require consideration of the base composition bias and the site-to-site rate variation in a maximum-likelihood framework.      

Evolution and phylogenetic utility of alignment gaps within intron sequences of three nuclear genes in bumble bees (Bombus)

To test whether gaps resulting from sequence alignment contain phylogenetic signal concordant with those of base substitutions, we analyzed the occurrence of indel mutations upon a well-resolved, substitution-based tree for three nuclear genes in bumble bees (Bombus, Apidae: Bombini). The regions analyzed were exon and intron sequences of long-wavelength rhodopsin (LW Rh), arginine kinase (ArgK), and elongation factor-1alpha (EF-1alpha) F2 copy genes. LW Rh intron had only a few uninformative gaps, ArgK intron had relatively long gaps that were easily aligned, and EF-1alpha intron had many short gaps, resulting in multiple optimal alignments. The unambiguously aligned gaps within ArgK intron sequences showed no homoplasy upon the substitution-based tree, and phylogenetic signals within ambiguously aligned regions of EF-1alpha intron were highly congruent with those of base substitutions. We further analyzed the contribution of gap characters to phylogenetic reconstruction by incorporating them in parsimony analysis. Inclusion of gap characters consistently improved support for nodes recovered by substitutions, and inclusion of ambiguously aligned regions of EF-1alpha intron resolved several additional nodes, most of which were apical on the phylogeny. We conclude that gaps are an exceptionally reliable source of phylogenetic information that can be used to corroborate and refine phylogenies hypothesized by base substitutions, at least at lower taxonomic levels. At present, full use of gaps in phylogenetic reconstruction is best achieved in parsimony analysis, pending development of well-justified and generally applicable methods for incorporating indels in explicitly model-based methods.     

Potential phylogenetic utility of the low-copy nuclear gene pistillata in dicotyledonous plants: comparison to nrDNA ITS and trnL intron in Sphaerocardamum and other Brassicaceae.

We report the potential phylogenetic utility of DNA sequence data from the last 700 bp of a ca. 1-kb intron of the MADS-box gene pistillata from a sampling of Sphaerocardamum species and other Brassicaceae. These results are compared with nrDNA ITS and the chloroplast trnL intron for the same taxa to demonstrate the potential phylogenetic utility of this pistillata intron and to identify potential historically independent sequences for an ongoing study of relationships within Sphaerocardamum. Analyses of the DNA sequence data for Brassicaceae indicated that pairwise divergences and potentially informative characters were higher in the pistillata intron (0.6-30.8%, 284 characters) and ITS (0-24%, 94 characters) than in the chloroplast trnL intron (0-4.2%, 17 characters). A comparison of Sphaerocardamum sequences identified low divergences and numbers of informative characters for trnL intron (0-2.4%, 1 character) and nrDNA ITS (0-2.5%, 2 characters) and substantially more variation among the pistillata sequences (0.15-3.7%, 19 characters). Phylogenetic analyses of these pistillata sequences fully resolve ingroup relationships without character conflict. Results of pistillata PCR amplifications from a broader dicot sample showed that some primers may be useful in amplifying orthologous pistillata sequences. Ultimately this pistillata intron may be a valuable source of phylogenetic characters at lower taxonomic levels.     

Organization of the chromosome region harboring a FLORICAULA/LEAFY gene in Liriodendron

FLORICAULA/LEAFY (FLO/LFY) plays an important role in the reproductive transition and controls flower spatial patterning by inducing the expression of the ABC floral organ identity genes. In this study, we sequenced two bacterial artificial chromosomes harboring a FLO/LFY and three other genes from yellow-poplar (Liriodendron tulipifera L.) and compared the gene order in this locus between several species. Besides the conserved terminal domains, key residues involved in interactions with DNA bases, backbone, and in dimerization were also conserved in the yellow-poplar FLO/LFY. Phylogenetic analysis of the FLO/LFY amino acid sequences placed yellow-poplar closer to eudicots than to monocotyledonous species. We found that gene content and order in this region of the yellow-poplar genome was more similar to corresponding regions in Vitis vinifera L., Carica papaya L., Populus trichocarpa Torr. & Gray, and Ricinus communis L., regardless of the evolutionary relationship. In addition, evidence for transposition, large insertions, and duplications were found, suggesting multiple and complex mechanisms of basal angiosperm genome evolution.     

Phylogenetic utility of the glycerol-3-phosphate acyltransferase gene: evolution and implications in Paeonia (Paeoniaceae)

The nuclear-encoded chloroplast-expressed glycerol-3-phosphate acyltransferase (GPAT) gene has been found to be single-copy in a number of angiosperm families. In this study we investigated the phylogenetic utility of the GPAT gene at the interspecific level using the genus Paeonia (Paeoniaceae) as an example. An approximately 2.3- to 2.6-kb fragment of the GPAT gene, containing a large intron of more than 2 kb, was amplified, cloned, and sequenced from 19 accessions representing 13 Paeonia species. The GPAT gene phylogeny inferred by parsimony analysis supported interspecific relationships that were previously unresolved, suggesting that large introns of low-copy nuclear genes are particularly informative in the resolution of close relationships at low taxonomic levels. Whereas the GPAT phylogeny is largely congruent with the previous phylogenetic hypothesis of Paeonia, it shows a significant discordance involving the paraphyly of section Paeonia. Given evidence of an ancient duplication and the subsequent silencing of one GPAT locus in P. anomala, this discordance is most likely the result of paralogy. Two distinct genomic clones containing partial GPAT genes were isolated from P. anomala. The GPAT sequence from one clone corresponded to the functional copy of the gene, and the second genomic clone was determined to contain a GPAT pseudogene. The insertion of a retrotransposon in an intron of this pseudogene may have been responsible for the silencing of this GPAT locus in P. anomala. This study suggests that, although it is unlikely that universal nuclear gene markers free from paralogy are usually available, low-copy nuclear genes can be very useful in plant phylogenetic reconstruction, especially at low taxonomic levels, as long as the evolutionary dynamics of the genes are carefully examined.     

Phylogenetic relationships among A-genome species of the genus Oryza revealed by intron sequences of four nuclear genes

The A-genome group in Oryza consists of eight diploid species and is distributed world-wide. Here we reconstructed the phylogeny among the A-genome species based on sequences of nuclear genes and MITE (miniature inverted-repeat transposable elements) insertions. Thirty-seven accessions representing two cultivated and six wild species from the A-genome group were sampled. Introns of four nuclear single-copy genes on different chromosomes were sequenced and analysed by both maximum parsimony (MP) and Bayesian inference methods. All the species except for Oryza rufipogon and Oryza nivara formed a monophyletic group and the Australian endemic Oryza meridionalis was the earliest divergent lineage. Two subspecies of Oryza sativa (ssp. indica and ssp. japonica) formed two separate monophyletic groups, suggestive of their polyphyletic origin. Based on molecular clock approach, we estimated that the divergence of the A-genome group occurred c. 2.0 million years ago (mya) while the two subspecies (indica and japonica) separated c. 0.4 mya. Intron sequences of nuclear genes provide sufficient resolution and are informative for phylogenetic inference at lower taxonomic levels.     

Actin Phylogeny and Intron Distribution in Bangiophyte Red Algae(Rhodoplantae)

The molecular phylogeny of red algal actin genes, with emphasis on the paraphyletic “Bangiophyceae,” was examined and compared to the rhodophyte SSU rDNA phylogeny. Nineteen new genomic actin sequences and seven SSU rDNA sequences were obtained and subjected to diverse phylogenetic analyses (maximum likelihood, distance/neighbor-joining, maximum parsimony, Bayesian analyses, and, with respect to protein sequences, also quartet puzzling). The actin trees confirmed most of the major clades found in the SSU rDNA phylogenies, although with a lower resolution. An actin gene duplication in the florideophycean lineage is reported, presumably related to an increased complexity of sexual reproduction. In addition, the distribution and characteristics of spliceosomal introns found in some of the actin sequences were examined. Introns were found in almost all florideophycean actin genes, whereas only two bangiophyte sequences contained introns. One intron in the florideophycean actin genes was also found in metazoan, and, shifted by one or two nucleotides, in a glaucocystophyte, a cryptophyte, and two fungal actin genes, and thus may be an ancient intron.[Reviewing Editor: Dr. Yves Van de Peer]     

The phylogenetic potential of entire 26S rDNA sequences in plants

18S ribosomal RNA genes are the most widely used nuclear sequences for phylogeny reconstruction at higher taxonomic levels in plants. However, due to a conservative rate of evolution, 18S rDNA alone sometimes provides too few phylogenetically informative characters to resolve relationships adequately. Previous studies using partial sequences have suggested the potential of 26S or large-subunit (LSU) rDNA for phylogeny retrieval at taxonomic levels comparable to those investigated with 18S rDNA. Here we explore the patterns of molecular evolution of entire 26S rDNA sequences and their impact on phylogeny retrieval. We present a protocol for PCR amplification and sequencing of entire (approximately 3.4 kb) 26S rDNA sequences as single amplicons, as well as primers that can be used for amplification and sequencing. These primers proved useful in angiosperms and Gnetales and likely have broader applicability. With these protocols and primers, entire 26S rDNA sequences were generated for a diverse array of 15 seed plants, including basal eudicots, monocots, and higher eudicots, plus two representatives of Gnetales. Comparisons of sequence dissimilarity indicate that expansion segments (or divergence domains) evolve 6.4 to 10.2 times as fast as conserved core regions of 26S rDNA sequences in plants. Additional comparisons indicate that 26S rDNA evolves 1.6 to 2.2 times as fast as and provides 3.3 times as many phylogenetically informative characters as 18S rDNA; compared to the chloroplast gene rbcL, 26S rDNA evolves at 0.44 to 1.0 times its rate and provides 2.0 times as many phylogenetically informative characters. Expansion segment sequences analyzed here evolve 1.2 to 3.0 times faster than rbcL, providing 1.5 times the number of informative characters. Plant expansion segments have a pattern of evolution distinct from that found in animals, exhibiting less cryptic sequence simplicity, a lower frequency of insertion and deletion, and greater phylogenetic potential.