Extensive ribosomal DNA genic variation in the columnar cactus Lophocereus

Age-Dependent and Tissue-Specific Alterations in the rDNA Clusters of the Panax ginseng C. A. Meyer Cultivated Cell Lines

Long-term cultivation of Panax ginseng cell lines leads to a decreasing synthesis of the biologically active substances used in traditional medicine. To gain insight into the cellular mechanisms which may influence this process, we analyzed variations within the rDNA cluster of the Oriental ginseng cell lines. The cell lines were cultivated for 6 and 24 years; the number of nucleoli and chromosomes was analyzed. The complete 18S rDNA sequences were cloned and sequenced. The nucleotide polymorphism and phylogenetic relations of the sequences were analyzed, and the secondary structures for separate 18S rRNA regions were modeled. The 18S rDNA accumulated mutations during cell cultivation that correlate well with an increase in the number of chromosomes and nucleoli. The patterns of nucleotide diversity are culture-specific and the increasing polymorphism associates with cytosine methylation sites. The secondary structures of some 18S rRNA regions and their interaction can alter during cultivation. The phylogenetic tree topologies are particular for each cell line.The observed alterations in rDNA clusters are associated with a somaclonal variation, leading to changes in the pattern of intracellular synthesis during cell cultivation. The identified divergent rRNAs could provide additional gene expression regulation in P. ginseng cells by forming heterogeneous ribosomes.

RIP mutated ITS genes in populations of Ophiocordyceps sinensis and their implications for molecular systematics

Different hypotheses have been proposed to interpret the observed unusual ITS (internal transcribed spacer) sequences in Ophiocordyceps sinensis. The coexistence of diverged ITS paralogs in a single genome was previously shown by amplifying the ITS region from mono-ascospore isolates using specific primers designed for different ITS paralog groups. Among those paralogs, are AT-biased ITS sequences which were hypothesized to result from repeat-induced point mutation (RIP). This is a process that detects and mutates repetitive DNA and frequently leads to epigenetic silencing, and these mutations have been interpreted as pseudogenes. Here we investigate the occurrence and frequency of ITS pseudogenes in populations of O. sinensis using large-scale sampling, and discusses the implications of ITS pseudogenes for fungal phylogenetic and evolutionary studies. Our results demonstrate a wide distribution of ITS pseudogenes amongst different geographic populations, and indicate how ITS pseudogenes can contribute to the reconstruction of the evolutionary history of the species.

Molecular evolution of the internal transcribed spacers in red oaks (Quercus sect. Lobatae)

Previous studies of the Internal Transcribed Spacers of the nuclear ribosomal DNA (ITS) in sections Quercus (white oaks), Protobalanus (intermediate or golden cup oaks), Cerris (Cerris oaks), and Ilex (Ilex oaks) suggest that ITS regions undergo full concerted evolution in oaks; however, ITS evolution patterns in red oaks (section Lobatae) are unknown due to scant representation in published work. To determine whether full concerted evolution occurs in red oaks, the purpose of this study was to examine ITS sequences from 40 red oak species. The results show incomplete concerted evolution and the presence of three ITS ribotypes of lengths 505, 609, 601 bp, hereafter referred to as ITS-S (small), I ITS-M (medium), and ITS-L (large), respectively. Thirty species had only one ribotype (ITS-M), nine species had two ribotypes (different combinations of ITS-L, ITS-M, and ITS-S), and only one species had all three ribotypes. Furthermore, examination of these three ribotypes showed that only ITS-M is putatively functional and ITS-L and ITS-S are pseudogenes. Bayesian analysis strongly supported (100%) two pseudogenes clades but provided weak support for the monophyly of a putative functional clade (ITS-M); moreover, within the "functional" clade, species relationships were uncertain and, in most cases, sequences from the same species failed to group together. The results of the current study suggest that ITS may not be appropriate for phylogeny reconstruction of red oaks due to low levels of interspecific variation and incomplete concerted evolution.

Disentangling sources of variation in SSU rDNA sequences from single cell analyses of ciliates: impact of copy number variation and experimental error

Small subunit ribosomal DNA (SSU rDNA) is widely used for phylogenetic inference, barcoding and other taxonomy-based analyses. Recent studies indicate that SSU rDNA of ciliates may have a high level of sequence variation within a single cell, which impacts the interpretation of rDNA-based surveys. However, sequence variation can come from a variety of sources including experimental errors, especially the mutations generated by DNA polymerase in PCR. In the present study, we explore the impact of four DNA polymerases on sequence variation and find that low-fidelity polymerases exaggerate the estimates of single-cell sequence variation. Therefore, using a polymerase with high fidelity is essential for surveys of sequence variation. Another source of variation results from errors during amplification of SSU rDNA within the polyploidy somatic macronuclei of ciliates. To investigate further the impact of SSU rDNA copy number variation, we use a high-fidelity polymerase to examine the intra-individual SSU rDNA polymorphism in ciliates with varying levels of macronuclear amplification: Halteria grandinella, Blepharisma americanum and Strombidium stylifer We estimate the rDNA copy numbers of these three species by single-cell quantitative PCR. The results indicate that: (i) sequence variation of SSU rDNA within a single cell is authentic in ciliates, but the level of intra-individual SSU rDNA polymorphism varies greatly among species; (ii) rDNA copy numbers vary greatly among species, even those within the same class; (iii) the average rDNA copy number of Halteria grandinella is about 567 893 (s.d. = 165 481), which is the highest record of rDNA copy number in ciliates to date; and (iv) based on our data and the records from previous studies, it is not always true in ciliates that rDNA copy numbers are positively correlated with cell or genome size.

ITS non-concerted evolution and rampant hybridization in the legume genus Lespedeza (Fabaceae)

The internal transcribed spacer (ITS) as one part of nuclear ribosomal DNA is one of the most extensively sequenced molecular markers in plant systematics. The ITS repeats generally exhibit high-level within-individual homogeneity, while relatively small-scale polymorphism of ITS copies within individuals has often been reported in literature. Here, we identified large-scale polymorphism of ITS copies within individuals in the legume genus Lespedeza (Fabaceae). Divergent paralogs of ITS sequences, including putative pseudogenes, recombinants, and multiple functional ITS copies were sometimes detected in the same individual. Thirty-seven ITS pseudogenes could be easily detected according to nucleotide changes in conserved 5.8S motives, the significantly lower GC contents in at least one of three regions, and the lost ability of 5.8S rDNA sequence to fold into a conserved secondary structure. The distribution patterns of the putative functional clones were highly different between the traditionally recognized two subgenera, suggesting different rates of concerted evolution in two subgenera which could be attributable to their different extents/frequencies of hybridization, confirmed by our analysis of the single-copy nuclear gene PGK. These findings have significant implications in using ITS marker for reconstructing phylogeny and studying hybridization.

Evolution of ketosynthase domains of polyketide synthase genes in the Cladonia chlorophaea species complex (Cladoniaceae)

Lichen-forming fungi synthesize a diversity of polyketides, but only a few non-reducing polyketide synthase (PKS) genes from a lichen-forming fungus have been linked with a specific polyketide. While it is a challenge to link the large number of PKS paralogs in fungi with specific products, it might be expected that the PKS paralogs from closely related species would be similar because of recent evolutionary divergence. The objectives of this study were to reconstruct a PKS gene phylogeny of the Cladonia chlorophaea species complex based on the ketosynthase domain, a species phylogeny of the complex, and to explore the presence of PKS gene paralogs among members of the species complex. DNA was isolated from 51 individuals of C. chlorophaea and allies to screen for the presence of 13 PKS paralogs. A 128 sequence PKS gene phylogeny using deduced amino acid sequences estimated from the 13 PKS paralogs and sequences subjected to BLASTx comparisons showed losses of each of two PKS domains (reducing and methylation). This research provided insight into the evolution of PKS genes in the C. chlorophaea group, species evolution in the group, and it identified potential directions for further investigation of polyketide synthesis in the C. chlorophaea species complex.

Multiple ITS copies reveal extensive hybridization within Rheum (Polygonaceae), a genus that has undergone rapid radiation

Background

During adaptive radiation events, characters can arise multiple times due to parallel evolution, but transfer of traits through hybridization provides an alternative explanation for the same character appearing in apparently non-sister lineages. The signature of hybridization can be detected in incongruence between phylogenies derived from different markers, or from the presence of two divergent versions of a nuclear marker such as ITS within one individual.

Methodology/Principal Findings

In this study, we cloned and sequenced ITS regions for 30 species of the genus Rheum, and compared them with a cpDNA phylogeny. Seven species contained two divergent copies of ITS that resolved in different clades from one another in each case, indicating hybridization events too recent for concerted evolution to have homogenised the ITS sequences. Hybridization was also indicated in at least two further species via incongruence in their position between ITS and cpDNA phylogenies. None of the ITS sequences present in these nine species matched those detected in any other species, which provides tentative evidence against recent introgression as an explanation. Rheum globulosum, previously indicated by cpDNA to represent an independent origin of decumbent habit, is indicated by ITS to be part of clade of decumbent species, which acquired cpDNA of another clade via hybridization. However decumbent and glasshouse morphology are confirmed to have arisen three and two times, respectively.

Conclusions

These findings suggested that hybridization among QTP species of Rheum has been extensive, and that a role of hybridization in diversification of Rheum requires investigation.

Non-concerted ITS evolution in fungi, as revealed from the important medicinal fungus Ophiocordyceps sinensis

The internal transcribed spacer (ITS) of nuclear ribosomal DNA (nrDNA) has been widely used as a molecular marker in phylogenetic studies and has been selected as a DNA barcode for fungi. It is generally believed that nrDNA conforms to concerted evolution in most eukaryotes; however, intraindividual-intraspecific polymorphisms of this region were reported in various organisms, suggesting a non-concerted evolutionary process. In Ophiocordyceps sinensis, one of the most valuable medicinal fungi, a remarkable variation of the ITS region has been revealed. Some highly divergent sequences were thought to represent cryptic species, different species or genotypes in previous studies. To clarify the unusual ITS polymorphisms observed in O. sinensis, specific primers were designed to amplify ITS paralogs from pure cultures of both single-ascospore and tissue isolates in this study. All of the available ITS sequences, including those generated by this group and those in GenBank, were analyzed. Several AT-biased ITS paralogs were classified as pseudogenes based on their nucleotide compositions, secondary structures and minimum free energies of their 5.8S rRNAs, substitution rates, phylogenetic positions and gene expression analyses. Furthermore, ITS pseudogenes were amplified with specific primers from 10 of the 28 strains tested, including eight single-ascospore and two tissue isolates. Divergent ITS paralogs were proved to coexist in individual genomes, suggesting a non-concerted mechanism of evolution in the ITS region of O. sinensis. The hypotheses that divergent ITS paralogs represent cryptic or other species or different genotypes were thus rejected.

Molecular characterization of 'Bhut Jolokia' the hottest chilli

The northeast region of India, considered as 'hot spot' of biodiversity, having unique ecological environment with hot and high-humidity conditions, has given rise to the world's hottest chilli, 'Bhut Jolokia', which is at least two times hotter than Red Savina Habanero in terms of Scoville heat units (SHU). This study was undertaken to determine the distinctiveness of 'Bhut Jolokia' from Capsicum frutescens or Capsicum chinense through sequencing of the ribosomal RNA (rRNA) gene-internal transcribed ((ITS) region along with its phylogenetic analysis. Although a compensatory base change (CBC) in the ITS2 region was not observed between the closely related species of C. frutescens and C. chinense when compared with Bhut Jolokia; phylogenetic analysis using ITS1, 5.8S and ITS2 sequences indicated a distinct clade for all the accessions of 'Bhut Joloikia', while C. frutescens and C. chinense occupied discrete lineages. Further, a unique 13-base deletion was observed in all the representative accessions of 'Bhut Jolokia', making it distinct from all other members within the genus and beyond. The degree of genetic variations along with its extreme pungency might be related to ambient environmental factors of northeastern India.

Intra-genomic ribosomal RNA polymorphism and morphological variation in Elphidium macellum suggests inter-specific hybridization in foraminifera

Elphidium macellum is a benthic foraminifer commonly found in the Patagonian fjords. To test whether its highly variable morphotypes are ecophenotypes or different genotypes, we analysed 70 sequences of the SSU rRNA gene from 25 specimens. Unexpectedly, we identified 11 distinct ribotypes, with up to 5 ribotypes co-occurring within the same specimen. The ribotypes differ by varying blocks of sequence located at the end of stem-loop motifs in the three expansion segments specific to foraminifera. These changes, distinct from typical SNPs and indels, directly affect the structure of the expansion segments. Their mosaic distribution suggests that ribotypes originated by recombination of two or more clusters of ribosomal genes. We propose that this expansion segment polymorphism (ESP) could originate from hybridization of morphologically different populations of Patagonian Elphidium. We speculate that the complex geological history of Patagonia enhanced divergence of coastal foraminiferal species and contributed to increasing genetic and morphological variation.

Molecular phylogeny of Gymnocalycium (Cactaceae): assessment of alternative infrageneric systems, a new subgenus, and trends in the evolution of the genus

PREMISE OF THE STUDY

The South American genus Gymnocalycium (Cactoideae-Trichocereae) demonstrates how the sole use of morphological data in Cactaceae results in conflicts in assessing phylogeny, constructing a taxonomic system, and analyzing trends in the evolution of the genus.

METHODS

Molecular phylogenetic analysis was performed using parsimony and Bayesian methods on a 6195-bp data matrix of plastid DNA sequences (atpI-atpH, petL-psbE, trnK-matK, trnT-trnL-trnF) of 78 samples, including 52 species and infraspecific taxa representing all the subgenera of Gymnocalycium. We assessed morphological character evolution using likelihood methods to optimize characters on a Bayesian tree and to reconstruct possible ancestral states.

KEY RESULTS

The results of the phylogenetic study confirm the monophyly of the genus, while supporting overall the available infrageneric classification based on seed morphology. Analysis showed the subgenera Microsemineum and Macrosemineum to be polyphyletic and paraphyletic. Analysis of morphological characters showed a tendency toward reduction of stem size, reduction in quantity and hardiness of spines, increment of seed size, development of napiform roots, and change from juicy and colorful fruits to dry and green fruits.

CONCLUSIONS

Gymnocalycium saglionis is the only species of Microsemineum and a new name is required to identify the clade including the remaining species of Microsemineum; we propose the name Scabrosemineum in agreement with seed morphology. Identifying morphological trends and environmental features allows for a better understanding of the events that might have influenced the diversification of the genus.

The 5S rDNA family evolves through concerted and birth-and-death evolution in fish genomes: an example from freshwater stingrays

Background

Ribosomal 5S genes are well known for the critical role they play in ribosome folding and functionality. These genes are thought to evolve in a concerted fashion, with high rates of homogenization of gene copies. However, the majority of previous analyses regarding the evolutionary process of rDNA repeats were conducted in invertebrates and plants. Studies have also been conducted on vertebrates, but these analyses were usually restricted to the 18S, 5.8S and 28S rRNA genes. The recent identification of divergent 5S rRNA gene paralogs in the genomes of elasmobranches and teleost fishes indicate that the eukaryotic 5S rRNA gene family has a more complex genomic organization than previously thought. The availability of new sequence data from lower vertebrates such as teleosts and elasmobranches enables an enhanced evolutionary characterization of 5S rDNA among vertebrates.

Results

We identified two variant classes of 5S rDNA sequences in the genomes of Potamotrygonidae stingrays, similar to the genomes of other vertebrates. One class of 5S rRNA genes was shared only by elasmobranches. A broad comparative survey among 100 vertebrate species suggests that the 5S rRNA gene variants in fishes originated from rounds of genome duplication. These variants were then maintained or eliminated by birth-and-death mechanisms, under intense purifying selection. Clustered multiple copies of 5S rDNA variants could have arisen due to unequal crossing over mechanisms. Simultaneously, the distinct genome clusters were independently homogenized, resulting in the maintenance of clusters of highly similar repeats through concerted evolution.

Conclusions

We believe that 5S rDNA molecular evolution in fish genomes is driven by a mixed mechanism that integrates birth-and-death and concerted evolution.

DNA barcodes for Mexican Cactaceae, plants under pressure from wild collecting

DNA barcodes could be a useful tool for plant conservation. Of particular importance is the ability to identify unknown plant material, such as from customs seizures of illegally collected specimens. Mexican cacti are an example of a threatened group, under pressure because of wild collection for the xeriscaping trade and private collectors. Mexican cacti also provide a taxonomically and geographically coherent group with which to test DNA barcodes. Here, we sample the matK barcode for 528 species of Cactaceae including approximately 75% of Mexican species and test the utility of the matK region for species-level identification. We find that the matK DNA barcode can be used to identify uniquely 77% of species sampled, and 79-87% of species of particular conservation importance. However, this is far below the desired rate of 95% and there are significant issues for PCR amplification because of the variability of primer sites. Additionally, we test the nuclear ITS regions for the cactus subfamily Opuntioideae and for the genus Ariocarpus (subfamily Cactoideae). We observed higher rates of variation for ITS (86% unique for Opuntioideae sampled) but a much lower PCR success, encountering significant intra-individual polymorphism in Ariocarpus precluding the use of this marker in this taxon. We conclude that the matK region should provide useful information as a DNA barcode for Cactaceae if the problems with primers can be addressed, but matK alone is not sufficiently variable to achieve species-level identification. Additional complementary regions should be investigated as ITS is shown to be unsuitable.

Cyto-nuclear discordance in the phylogeny of Ficus section Galoglychia and host shifts in plant-pollinator associations

BACKGROUND

Hybridization events are relatively common in vascular plants. However, the frequency of these events is unevenly distributed across the plant phylogeny. Plant families in which individual species are pollinated by specific pollinator species are predicted to be less prone to hybridization than other families. However, exceptions may occur within these families, when pollinators shift host-plant species. Indeed, host shifts are expected to increase the rate of hybridization events. Pollinators of Ficus section Galoglychia are suspected to have changed host repeatedly, based on several cases of incongruence between plant phylogeny and taxonomy, and insect phylogeny and taxonomy. We tracked cyto-nuclear discordance across section Galoglychia as evidence for hybridization. To achieve a proper global view, we first clarified the monophyly of section Galoglychia as it had been questioned by recent phylogenetic studies. Moreover, we investigated if fig size could be a factor facilitating host shifts.

RESULTS

Phylogenetic chloroplast and nuclear results demonstrated the monophyly of section Galoglychia. Within section Galoglychia, we detected several cases of statistically significant cyto-nuclear discordance. Discordances concern both terminal nodes of the phylogenetic trees and one deep node defining relationships between subsections. Because nuclear phylogeny is congruent with morphological taxonomy, discordances were caused by the chloroplast phylogeny. Introgressive hybridization was the most likely explanation for these discordances. We also detected that subsections pollinated by several wasp genera had smaller figs and were pollinated by smaller wasps than subsections pollinated by a single wasp genus.

CONCLUSION

As hypothesized, we discovered evidences of past hybridization in Ficus section Galoglychia. Further, introgression was only detected in subsections presenting incongruence between plant and pollinator phylogenies and taxonomy. This supports the hypothesis that host shift is the cause for plant-pollinator incongruence. Moreover, small fig size could facilitate host shifts. Eventually, this study demonstrates that non-coding chloroplast markers are valuable to resolve deep nodes in Ficus phylogeny.

Molecular identification and hidden diversity of novel Daphnia parasites from European lakes

Parasites play important roles in local population dynamics and genetic structure. However, due to insufficient diagnostic tools, detailed host-parasite interactions may remain concealed by hidden parasite diversity in natural systems. Microscopic examination of 19 European lake Daphnia populations revealed the presence of three groups of parasites: fungi, microsporidia, and oomycetes. For most of these parasites no genetic markers have been described so far. Based on sequence similarities of the nuclear small-subunit and internal transcribed spacer (ITS) rRNA gene regions, one fungus, four microsporidian, and nine oomycete taxa were discovered in 147 infected Daphnia (and/or three other zooplankton crustaceans). Additionally, cloning of rRNA gene regions revealed parasite sequence variation within host individuals. This was most pronounced in the ITS region of one microsporidian taxon, where the within-host sequence variation ranged from 1.7% to 5.3% polymorphic sites for parasite isolates from 14 different geographical locations. Interestingly, the parasite isolates from close locations grouped together based on sequence similarities, suggesting that there was parasite dispersal. Taken together, the data obtained in this study revealed hidden diversity of parasite communities in Daphnia lake populations. Moreover, a higher level of resolution for identifying parasite strains makes it possible to test new hypotheses with respect to parasite dispersal, transmission routes, and coinfection.

Intra-individual and intra-species heterogeneity in nuclear rDNA ITS region of Vigna species from subgenus Ceratotropis

The extent of intra-individual and intra-species heterogeneity in the nuclear rDNA internal transcribed spacer (ITS) was investigated among the 'Asiatic Vigna' species (subgenus Ceratotropis). High intra- and inter-individual ITS polymorphism was observed among Vigna radiata accessions, where multiple ITS length variants ranging from approximately 700 to approximately 770 bp were detected on PCR amplification. Subsequent analysis revealed that the variants are 'heteroduplex ITS fragments' generated during the PCR process. Analysis of ITS from wild and cultivated forms of ten Vigna species from subgenus Ceratotropis revealed substantial intra-species divergence in four species: Vigna umbellata, Vigna trilobata, V. radiata and Vigna minima. However, no other species analysed showed intra-individual ITS heterogeneity as observed in V. radiata. The results demonstrate differential evolution of ITS sequence among wild and cultivated forms of V. radiata. Evidence indicates that intra-species hybridization and a slow 'molecular drive' are responsible for this phenomenon. Sequence analysis of 5.8S, ITS1 and ITS2 and secondary-structure analysis of ITS regions indicate that the ITS variants do not belong to pseudogenic rDNA repeat units. Further, reverse transcriptase-PCR (RT-PCR) analysis showed that rDNA repeat units harbouring certain intra-individual ITS variants were transcriptionally inactive, indicating the regulation of these loci by epigenetic gene silencing. The V. radiata ITS variants, when analysed together, did not cause any phylogenetic errors at the species level.

Extensive 5.8S nrDNA polymorphism in Mammillaria (Cactaceae) with special reference to the identification of pseudogenic internal transcribed spacer regions

The internal transcribed spacer (ITS) region (ITS1, 5.8S rDNA, ITS2) represents the most widely applied nuclear marker in eukaryotic phylogenetics. Although this region has been assumed to evolve in concert, the number of investigations revealing high degrees of intra-individual polymorphism connected with the presence of pseudogenes has risen. The 5.8S rDNA is the most important diagnostic marker for functionality of the ITS region. In Mammillaria, intra-individual 5.8S rDNA polymorphisms of up to 36% and up to nine different types have been found. Twenty-eight of 30 cloned genomic Mammillaria sequences were identified as putative pseudogenes. For the identification of pseudogenic ITS regions, in addition to formal tests based on substitution rates, we attempted to focus on functional features of the 5.8S rDNA (5.8S motif, secondary structure). The importance of functional data for the identification of pseudogenes is outlined and discussed. The identification of pseudogenes is essential, because they may cause erroneous phylogenies and taxonomic problems.

Transition-transversion bias is not universal: a counter example from grasshopper pseudogenes

Comparisons of the DNA sequences of metazoa show an excess of transitional over transversional substitutions. Part of this bias is due to the relatively high rate of mutation of methylated cytosines to thymine. Postmutation processes also introduce a bias, particularly selection for codon-usage bias in coding regions. It is generally assumed, however, that there is a universal bias in favour of transitions over transversions, possibly as a result of the underlying chemistry of mutation. Surprisingly, this underlying trend has been evaluated only in two types of metazoan, namely Drosophila and the Mammalia. Here, we investigate a third group, and find no such bias. We characterize the point substitution spectrum in Podisma pedestris, a grasshopper species with a very large genome. The accumulation of mutations was surveyed in two pseudogene families, nuclear mitochondrial and ribosomal DNA sequences. The cytosine-guanine (CpG) dinucleotides exhibit the high transition frequencies expected of methylated sites. The transition rate at other cytosine residues is significantly lower. After accounting for this methylation effect, there is no significant difference between transition and transversion rates. These results contrast with reports from other taxa and lead us to reject the hypothesis of a universal transition/transversion bias. Instead we suggest fundamental interspecific differences in point substitution processes.

Some mutations occur more frequently than others. We need to understand these biases if we are to interpret the differences that have accumulated between species and individuals. Applications include estimating the time since evolutionary lineages diverged and detecting the signature of natural selection in DNA sequences. However, mutational biases have been obscured because, since mutations arose, natural selection has eliminated some whilst allowing others to persist to the present. We therefore study mutations that have accumulated in regions of the genome that are free from selection in a grasshopper with a gigantic genome. All other animal studies using this approach find an excess of mutations between DNA bases having similar biochemical properties (transitions rather than transversions). This bias has been widely interpreted as a consequence of the fundamental biochemical basis of mutation. However, once we exclude mutations associated with DNA methylation, we find no evidence of a transition bias, unlike the few comparable animal studies that make the same correction. We propose that this result indicates previously unanticipated differences between species in the selection on or mutation of their DNA.

A Phylogeographic Split in Buxus balearica (Buxaceae) as Evidenced by Nuclear Ribosomal Markers: When ITS Paralogues Are Welcome

Sequences from the ribosomal nuclear internal transcribed spacers (ITS) have been widely used to infer evolutionary hypotheses across a broad range of living organisms. Intraspecific sequence variation is assumed to be absent or negliable in most species, but few detailed studies have been conducted to assess the apportionment of ITS sequence variation within and between plant populations. Buxus balearica was chosen as a model species to assess the levels of infraspecific and intragenomic ITS variation in rare and endangered species occurring in disjunct populations around the Mediterranean basin. Intragenomic polymorphic sites were detected for western and eastern accessions of B. balearica and in two accessions of the sister species B. sempervirens. Overall, 19 different ribotypes were found in B. balearica after sequencing 48 clones, whereas 15 ribotypes were detected in 19 clones of B. sempervirens. The integrity and secondary structure stability of the ribosomal sequences suggest that they are not pseudogenes. The high number of ribotypes recovered through cloning suggested that some sequences could be chimeric or generated in vivo by partial homogenization through gene conversion or unequal crossing-over. Average sequence divergence among B. balearica clones was 0.768%, and the most divergent sequences differed by 1.62%. Available evidence does not suggest that B. balearica paralogues have been obtained from other extant Buxus species through interspecific hybridization. The presence of several ribosomal sequences in box implies that the molecular forces driving the concerted evolution of this multigene family are not fully operational in this genus. Phylogenetic analyses of cloned ITS sequences from B. balearica displayed very poor resolution and only two clades received moderate bootstrap support. Despite the marked intragenomic sequence divergence found, ribosomal data suggest a clear phylogeographic split in B. balearica between western and eastern accessions. The distinct, nonchimeric sequences that are postulated as being present in each biogeographic group suggest that box populations from Anatolia (eastern Mediterranean) are relict.

Ribosomal DNA in the grasshopper Podisma pedestris: escape from concerted evolution

Eukaryote nuclear ribosomal DNA (rDNA) typically exhibits strong concerted evolution: a pattern in which several hundred rDNA sequences within any one species show little or no genetic diversity, whereas the sequences of different species diverge. We report a markedly different pattern in the genome of the grasshopper Podisma pedestris. Single individuals contain several highly divergent ribosomal DNA groups. Analysis of the magnitude of divergence indicates that these groups have coexisted in the Podisma lineage for at least 11 million years. There are two putatively functional groups, each estimated to be at least 4 million years old, and several pseudogene groups, many of which are transcribed. Southern hybridization and real-time PCR experiments show that only one of the putatively functional types occurs at high copy number. However, this group is scarcely amplified under standard PCR conditions, which means that phylogenetic inference on the basis of standard PCR would be severely distorted. The analysis suggests that concerted evolution has been remarkably ineffective in P. pedestris. We propose that this outcome may be related to the species' exceptionally large genome and the associated low rate of deletion per base pair, which may allow pseudogenes to persist.

Toward a Global Phylogeny of the Brassicaceae

The Brassicaceae is a large plant family (338 genera and 3,700 species) of major scientific and economic importance. The taxonomy of this group has been plagued by convergent evolution in nearly every morphological feature used to define tribes and genera. Phylogenetic analysis of 746 nrDNA internal transcribed spacer (ITS) sequences, representing 24 of the 25 currently recognized tribes, 146 genera, and 461 species of Brassicaceae, produced the most comprehensive, single-locus-based phylogenetic analysis of the family published to date. Novel approaches to nrDNA ITS analysis and extensive taxonomic sampling offered a test of monophyly for a large complement of the currently recognized tribes and genera of Brassicaceae. In the most comprehensive analysis, tribes Alysseae, Anchonieae plus Hesperideae, Boechereae, Cardamineae, Eutremeae, Halimolobeae, Iberideae, Noccaeeae, Physarieae, Schizopetaleae, Smelowskieae, and Thlaspideae were all monophyletic. Several broadly defined genera (e.g., Draba and Smelowskia) were supported as monophyletic, whereas others (e.g., Sisymbrium and Alyssum) were clearly polyphyletic. Analyses of ITS data identified several problematic sequences attributable to errors in sample identification or database submission. Results from parsimony ratchet and Bayesian analyses recovered little support for the backbone of the phylogeny, suggesting that many lineages of Brassicaceae have undergone rapid radiations that may ultimately be difficult to resolve with any single locus. However, the development of a preliminary supermatrix including the combination of 10 loci for 65 species provides an initial estimate of intertribal relations and suggests that broad application of such a method will provide greater understanding of relationships in the family.

Non-concerted ITS evolution in Mammillaria (Cactaceae)

Molecular studies of 21 species of the large Cactaceae genus Mammillaria representing a variety of intrageneric taxonomic levels revealed a high degree of intra-individual polymorphism of the internal transcribed spacer region (ITS1, 5.8S rDNA, ITS2). Only a few of these ITS copies belong to apparently functional genes, whereas most are probably non-functional (pseudogenes). As a multiple gene family, the ITS region is subjected to concerted evolution. However, the high degree of intra-individual polymorphism of up to 36% in ITS1 and up to 35% in ITS2 suggests a non-concerted evolution of these loci in Mammillaria. Conserved angiosperm motifs of ITS1 and ITS2 were compared between genomic and cDNA ITS clones of Mammillaria. Some of these motifs (e.g., ITS1 motif 1, 'TGGT' within ITS2) in combination with the determination of GC-content, length comparisons of the spacers and ITS2 secondary structure (helices II and III) are helpful in the identification of pseudogene rDNA regions.

Nuclear ribosomal DNA internal transcribed spacer 1 (ITS1) in Picea (Pinaceae): sequence divergence and structure

The nrDNA ITS1 of Picea is 2747-3271 bp, the longest known of all plants. We obtained 24 cloned ITS1 sequences from six individuals of Picea glehnii, Picea mariana, Picea orientalis, and Picea rubens. Mean sequence divergence within these individuals (0.018+/-0.009) is more than half that between the species (0.031+/-0.011) and may be maintained against concerted evolution by separation of Picea 18S-26S rDNA repeats on multiple chromosomes. Picea ITS1 contains three subrepeats with a motif (5'-GGCCACCCTAGTC) that is conserved across Pinaceae. Two subrepeats are tandem, remote from the third, and more closely related and significantly more similar to one another than either is to the third subrepeat. This correlation between similarity and proximity may be the result of subrepeat duplication or concerted evolution within rDNA repeats. In inferred secondary structures, subrepeats generally form long hairpins, with a portion of the Pinaceae conserved motif in the terminal loop, and tandem subrepeats pair with one another over most of their length. Coalescence of ITS1 sequences occurs in P. orientalis but not in the other species.

Recent origin and phylogenetic utility of divergent ITS putative pseudogenes: a case study from Naucleeae (Rubiaceae)

The internal transcribed spacer (ITS) of nuclear ribosomal DNA has been widely used by systematists for reconstructing phylogenies of closely related taxa. Although the occurrence of ITS putative pseudogenes is well documented for many groups of animals and plants, the potential utility of these pseudogenes in phylogenetic analyses has often been underestimated or even ignored in part because of deletions that make unambiguous alignment difficult. In addition, long branches often can lead to spurious relationships, particularly in parsimony analyses. We have discovered unusually high levels of ITS polymorphism (up to 30%, 40%, and 14%, respectively) in three tropical tree species of the coffee family (Rubiaceae), Adinauclea fagifolia, Haldina cordifolia, and Mitragyna rubrostipulata. Both secondary structure stability and patterns of nucleotide substitutions in a highly conserved region (5.8S gene) were used for distinguishing presumed functional sequences from putative pseudogenes. The combination of both criteria was the most powerful approach. The sequences from A. fagifolia appear to be a mix of functional genes and highly distinct putative pseudogenes, whereas those from H. cordifolia and M. rubrostipulata were identified as putative pseudogenes. We explored the potential utility of the identified putative pseudogenes in the phylogenetic analyses of Naucleeae sensu lato. Both Bayesian and parsimony trees identified the same monophyletic groups and indicated that the polymorphisms do not transcend species boundaries, implying that they do not predate the divergence of these three species. The resulting trees are similar to those produced by previous analyses of chloroplast genes. In contrast to results of previous studies therefore, divergent putative pseudogenes can be useful for phylogenetic analyses, especially when no sequences of their functional counterparts are available. Our studies clearly show that ITS polymorphism may not necessarily mislead phylogenetic inference. Despite using many different PCR conditions (different primers, higher denaturing temperatures, and absence or presence of DMSO and BSA-TMACl), we recovered only a few functional ITS copies from A. fagifolia and none from H. cordifolia and M. rubrostipulata, which suggests that PCR selection is occurring and/or the presumed functional alleles are located at minor loci (with few ribosomal DNA copies).

The rDNA ITS region in the lessepsian marine angiosperm Halophila stipulacea (Forssk.) Aschers. (Hydrocharitaceae): intragenomic variability and putative pseudogenic sequences

Halophila stipulacea is a dioecious marine angiosperm, widely distributed along the western coasts of the Indian Ocean and the Red Sea. This species is thought to be a Lessepsian immigrant that entered the Mediterranean Sea from the Red Sea after the opening of the Suez Canal (1869). Previous studies have revealed both high phenotypic and genetic variability in Halophila stipulacea populations from the western Mediterranean basin. In order to test the hypothesis of a Lessepsian introduction, we compare genetic polymorphism between putative native (Red Sea) and introduced (Mediterranean) populations through rDNA ITS region (ITS1-5.8S-ITS2) sequence analysis. A high degree of intraindividual variability of ITS sequences was found. Most of the intragenomic polymorphism was due to pseudogenic sequences, present in almost all individuals. Features of ITS functional sequences and pseudogenes are described. Possible causes for the lack of homogenization of ITS paralogues within individuals are discussed.

Ribosomal ITS sequences and plant phylogenetic inference

One of the most popular sequences for phylogenetic inference at the generic and infrageneric levels in plants is the internal transcribed spacer (ITS) region of the 18S-5.8S-26S nuclear ribosomal cistron. The prominence of this source of nuclear DNA sequence data is underscored by a survey of phylogenetic publications involving comparisons at the genus level or below, which reveals that of 244 papers published over the last five years, 66% included ITS sequence data. Perhaps even more striking is the fact that 34% of all published phylogenetic hypothesis have been based exclusively on ITS sequences. Notwithstanding the many important contributions of ITS sequence data to phylogenetic understanding and knowledge of genome relationships, a number of molecular genetic processes impact ITS sequences in ways that may mislead phylogenetic inference. These molecular genetic processes are reviewed here, drawing attention to both underlying mechanism and phylogenetic implications. Among the most prevalent complications for phylogenetic inference is the existence in many plant genomes of extensive sequence variation, arising from ancient or recent array duplication events, genomic harboring of pseudogenes in various states of decay, and/or incomplete intra- or inter-array homogenization. These phenomena separately and collectively create a network of paralogous sequence relationships potentially confounding accurate phylogenetic reconstruction. Homoplasy is shown to be higher in ITS than in other DNA sequence data sets, most likely because of orthology/paralogy conflation, compensatory base changes, problems in alignment due to indel accumulation, sequencing errors, or some combination of these phenomena. Despite the near-universal usage of ITS sequence data in plant phylogenetic studies, its complex and unpredictable evolutionary behavior reduce its utility for phylogenetic analysis. It is suggested that more robust insights are likely to emerge from the use of single-copy or low-copy nuclear genes.

Historical vicariance and postglacial colonization effects on the evolution of genetic structure in Lophocereus, a Sonoran Desert columnar cactus

Distinguishing the historical effects of gene migration and vicariance on contemporary genetic structure is problematic without testable biogeographic hypotheses based on preexisting geological and environmental evidence. The availability of such hypotheses for North America's Sonoran Desert has contributed to our understanding of the effect of historical vicariance and dispersal events on the diversification of this region's vertebrate biota but have not yet been applied to its flora. In this paper we describe a detailed allozyme analysis of the population genetic structure and phylogeography of the Sonoran Desert columnar cactus, Lophocereus schottii (senita). Inferred phylogroup distributions reflect two historical vicariance events: (1) a middle Pliocene northward transgression of the Sea of Cortéz that is reflected in well-supported Baja California peninsular and continental phylogroups but not in current taxonomic treatments of the species; and (2) a late Pliocene transpeninsular seaway across southern Baja that is reflected in tentative support for peninsular and southern Cape Region phylogroups corresponding to taxonomic varieties L. schottii var. schottii and L. schottii var. australis, respectively. A middle Pleistocene midpeninsular seaway hypothesized to explain congruent phylogroup distributions in several vertebrate taxa is not reflected in L. s. var. schottii, nor is the distinction of a third variety, L. s. var. tenuis, from continental populations of L. s. var. schottii. Linear regression of pairwise estimates of interpopulation differentiation (M and F(ST)/[1 - F(ST)]) on interpopulation geographic distance revealed significant evidence of isolation by distance within peninsular and continental phylogroups but not between them, consistent with historical vicariance between but not within these regions. We also found significant evidence of isolation by distance between putative L. s. var. schottii and L. s. var. australis phylogroups, suggesting that reproductive isolation between peninsular and Cape Region forms is incomplete. Within peninsular, but not continental, phylogroups, northward range expansion from southern Pleistocene refugia is reflected in significant declines in genetic variation with increasing latitude and in an area phenogram in which populations are progressively nested from south (ancestral) to north (descendant) along the Baja peninsula. Although the geographic concordance of phylogenetic topologies suggests that ancient vicariance events, and not dispersal, have primarily influenced the biogeographic distributions of Baja's vertebrate biota, the phylogeographic structure of L. schottii suggests that Sonoran Desert plant species may exhibit genetic signatures of postglacial range expansion and gene flow as well as vicariance.

Phylogenetic origins of Lophocereus (Cactaceae) and the senita cactus-senita moth pollination mutualism

Recent ecological research has revealed that the Sonoran Desert columnar cactus Lophocereus and the pyralid moth Upiga virescens form an obligate pollination mutualism, a rare but important case of coevolution. To investigate the phylogenetic origins of this unusual pollination system, we used molecular sequence data to reconstruct the phylogeny of the four taxa within the genus Lophocereus and to determine the phylogenetic position of Lophocereus within the North American columnar cacti (tribe Pachycereeae). Our analysis included Lophocereus, six Pachycereus species, Carnegiea gigantea, and Neobuxbaumia tetetzo within the subtribe Pachycereinae, and Stenocereus thurberi as an outgroup within the Stenocereinae. Extensive screening of chloroplast and mitochondrial genomes failed to reveal sequence variation within Lophocereus. At a deeper phylogenetic level, however, we found strong support for the placement of Lophocereus within Pachycereus as sister group to the hummingbird-pollinated P. marginatus. We discuss possible hypotheses that may explain the transition from bat pollination (ancestral) to moth and hummingbird pollination in Lophocereus and P. marginatus, respectively. Additional phylogenetic analyses suggest that the genus Pachycereus should be expanded to include Lophocereus, Carnegiea, Neobuxbaumia, and perhaps other species, whereas P. hollianus may need to be excluded from this clade. Future study will be needed to test taxonomic distinctions within Lophocereus, to test for parallel cladogenesis between phylogroups within Lophocereus and Upiga, and to fully delineate the genus Pachycereus and relationships among genera in the Pachycereinae.