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).     

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.     

Extensive 5.8S nrDNA polymorphism in <Emphasis Type="Italic">Mammillaria</Emphasis> (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. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

The Evolution of Ribosomal DNA: Divergent Paralogues and Phylogenetic Implications

Although nuclear ribosomal DNA (rDNA) repeats evolve together through concerted evolution, some genomes contain a considerable diversity of paralogous rDNA. This diversity includes not only multiple functional loci but also putative pseudogenes and recombinants. We examined the occurrence of divergent paralogues and recombinants in Gossypium, Nicotiana, Tripsacum, Winteraceae, and Zea ribosomal internal transcribed spacer (ITS) sequences. Some of the divergent paralogues are probably rDNA pseudogenes, since they have low predicted secondary structure stability, high substitution rates, and many deamination-driven substitutions at methylation sites. Under standard PCR conditions, the low stability paralogues amplified well, while many high-stability paralogues amplified poorly. Under highly denaturing PCR conditions (i.e., with dimethylsulfoxide), both low- and high-stability paralogues amplified well. We also found recombination between divergent paralogues. For phylogenetics, divergent ribosomal paralogues can aid in reconstructing ancestral states and thus serve as good outgroups. Divergent paralogues can also provide companion rDNA phylogenies. However, phylogeneticists must discriminate among families of divergent paralogues and recombinants or suffer from muddled and inaccurate organismal phylogenies.     

Nuclear Ribosomal ITS Functional Paralogs Resolve the Phylogenetic Relationships of a Late-Miocene Radiation Cycad Cycas (Cycadaceae)

Cycas is the most widespread and diverse genus among the ancient cycads, but the extant species could be the product of late Miocene rapid radiations. Taxonomic treatments to date for this genus are quite controversial, which makes it difficult to elucidate its evolutionary history. We cloned 161 genomic ITS sequences from 31 species representing all sections of Cycas. The divergent ITS paralogs were examined within each species and identified as putative pseudogenes, recombinants and functional paralogs. Functional paralogs were used to reconstruct phylogenetic relationships with pseudogene sequences as molecular outgroups, since an unambiguous ITS sequence alignment with their closest relatives, the Zamiaceae, is unachievable. A fully resolved and highly supported tree topology was obtained at the section level, with two major clades including six minor clades. The results fully supported the classification scheme proposed by Hill (2004) at the section level, with the minor clades representing his six sections. The two major clades could be recognised as two subgenera. The obtained pattern of phylogenetic relationships, combined with the different seed dispersal capabilities and paleogeography, allowed us to propose a late Miocene rapid radiation of Cycas that might have been promoted by vicariant events associated with the complex topography and orogeny of South China and adjacent regions. In contrast, transoceanic dispersals might have played an important role in the rapid diversification of sect. Cycas, whose members have evolved a spongy layer in their seeds aiding water dispersals.     

四照花亚属的nrDNA ITS致同进化不完全

四照花亚属(Cornus subg.Syncarpea)隶属于山茱萸科山茱萸属(Cornus),我国该亚属共有5种8亚种。为探讨四照花亚属nrDNA ITS序列的致同进化不完全现象及假基因产生的可能原因,分析了该亚属4种(每种1~2个居群)共21个个体的nrDNA ITS序列。结果表明,这些类群的nrDNA ITS存在多态性,通过分析这些nrDNA ITS克隆序列的G+C含量、5.8S保守基序和二级结构最小自由能,推测其可能存在假基因。系统发育研究结果显示所有nrDNA ITS序列分成5个分支,同一个体的不同拷贝被分别置于两个甚至多个分支中,且不同分支显示了不同种间关系。四照花亚属物种个体内部存在nrDNA ITS不完全致同进化,可能归咎于不完全的世系分选(incomplete lineage sorting)、种间杂交或多倍化等进化事件,从而导致基因组内nrITS区序列出现多态性,同时也导致难以通过外部形态来划分亚属内种间界限。     

Phylogeographic congruence between mtDNA and rDNA ITS markers in brown trout

Variation in the internal transcribed spacer (ITS) of rDNA was examined throughout the range of the brown trout (Salmo trutta) to analyze the usefulness of this molecular marker for phylogeographic analysis. The results were compared with those previously obtained with mtDNA, a region exhaustively analyzed along the brown trout distribution. ITS2 was essentially conserved at all populations sampled, no informative characters being detected across the main lineages described in this species. Conversely, ITS1 showed a greater homogenization than other genetic markers at a microgeographic scale, with variation partitioning into several major phylogenetic groups. Phylogeographic patterns were partially congruent between both ITS1 and mtDNA. The main discrepancies were the detection of intra-individual variation and putative recombinant ITS1 sequences in hybridization areas between genetically different, yet historically overlapping, assemblages. Also, the existence of an ancient ITS1 sequence in the Mediterranean-southeastern area (rMEDA), not revealed by mtDNA analysis, was evidenced after rDNA ITS1 analysis.     

Characterization of angiosperm nrDNA polymorphism, paralogy, and pseudogenes

Many early reports of ITS region (ITS 1, 5.8S, and ITS 2) variation in flowering plants indicated that nrDNA arrays within individuals are homogeneous. However, both older and more recent studies have found intra-individual nrDNA polymorphism across a range of plant taxa including presumed non-hybrid diploids. In addition, polymorphic individuals often contain potentially non-functional nrDNA copies (pseudogenes). These findings suggest that complete concerted evolution should not be assumed when embarking on phylogenetic studies using nrDNA sequences. Here we (1). discuss paralogy in relation to species tree reconstruction and conclude that a priori determinations of orthology and paralogy of nrDNA sequences should not be made based on the functionality or lack of functionality of those sequences; (2). discuss why systematists might be particularly interested in identifying and including pseudogene sequences as a test of gene tree sampling; (3). examine the various definitions and characterizations of nrDNA pseudogenes as well as the relative merits and limitations of a subset of pseudogene detection methods and conclude that nucleotide substitution patterns are particularly appropriate for the identification of putative nrDNA pseudogenes; and (4). present and discuss the advantages of a tree-based approach to identifying pseudogenes based on comparisons of sequence substitution patterns from putatively conserved (e.g., 5.8S) and less constrained (e.g., ITS 1 and ITS 2) regions. Application of this approach, through a method employing bootstrap hypothesis testing, and the issues discussed in the paper are illustrated through reanalysis of two previously published matrices. Given the apparent robustness of the test developed and the ease of carrying out percentile bootstrap hypothesis tests, we urge researchers to employ this statistical tool. While our discussion and examples concern the literature on plant systematics, the issues addressed are relevant to studies of nrDNA and other multicopy genes in other taxa.     

云南山茶(Camellia reticulata) nrDNA ITS序列多态性分析

对云南山茶 (Camellia reticulata) 8个品种的核核糖体DNA内转录间隔区(nrDNA ITS)进行克隆测序,将获得的ITS序列进行GC含量、 5.8S区二级结构的稳定性、替代模式、核苷酸多样性及系统发育关系的相关分析.实验结果显示,云南山茶8个品种的ITS序列存在丰富的基因组内多态性,同时包含中性进化的假基因,表明其ITS序列逃离了一致性进化.云南山茶ITS序列多态性的原因可能来自广泛的种间杂交,以及rDNA位点在基因组中有不确定的物理位置.ITS假基因为品种的物种形成研究提供了更全面的遗传证据,同时也提示了利用ITS假基因进行系统发育分析可能会对其真实的系统关系造成影响.     

Recolonization and radiation in Larix (Pinaceae): evidence from nuclear ribosomal DNA paralogues

Gene paralogy frequently causes the conflict between gene tree and species tree, but sometimes the coexistence of a few paralogous copies could provide more markers for tracing the phylogeographical process of some organisms. In the present study, nrDNA ITS paralogues were cloned from all but one species of Larix, an Eocene genus having two sections, Larix and Multiserialis, with a huge circumboreal distribution and an Eastern Asia-Western North America disjunction, respectively. A total of 96 distinct clones, excluding five putative pseudogenes or recombinants, were obtained and used in the gene genealogy analysis. The clones from all Eurasian species of section Larix are mixed together, suggesting that recolonization and recent morphological differentiation could have played important roles in the evolution of this section. In contrast, the species diversification of the Eurasian section Multiserialis may result from radiation in the east Himalayas and its vicinity, considering extensive nrDNA founder effects in this group. Our study also suggests that the distribution pattern analysis of members of multiple gene family would be very useful in tracking the evolutionary history of some taxa with recent origin or rapid radiation that cannot be resolved by other molecular markers.     

Intra-specific heterogeneity of the rDNA internal transcribed spacer in the Simulium damnosum (Diptera: Simuliidae) complex

The internal transcribed spacer (ITS) of the rRNA gene cluster has been used as a model for the study of the action of concerted evolution and molecular drive on repeated sequence families. In contrast to this general finding, preliminary DNA sequence analysis of cloned representatives of the ITS from the West African black fly species complex Simulium damnosum s.1. demonstrated extensive intra-individual and intra-specific polymorphisms. Variability in the ITS was primarily confined to the ITS1 domain. The degree and type of intra-individual and intra-specific variability within the ITS was further characterized using gel electrophoresis, DNA hybridization, and heteroduplex analysis of the PCR products generated from the ITS1 domain. ITS1 copies from individual S. damnosum s.1. differed in length and sequence composition. These results, when taken together, demonstrate that a large degree of intra-individual and intra-specific heterogeneity exists in the ITS of S. damnosum s.1. The intra-individual heterogeneity was greater in the savanna-dwelling than forest-dwelling sibling species of S. damnosum s.1. This heterogeneity may be due in part to inter-breeding among sympatric sibling species, coupled with disturbance of S. damnosum s.1. populations resulting from intensive vector control efforts.      

苏铁nrDNA ITS区的序列多态性:不完全致同进化的证据

本研究对苏铁(Cycas revoluta) nrDNA ITS进行克隆测序, 并以cDNA ITS为参照, 比较分析获得的序列的碱基变异、GC含量、5.8S二级结构的稳定性和5.8S保守基序的有无以及系统发育关系。结果发现苏铁nrDNA ITS存在较高的基因组内多样性, 同时, 这些分化的nrDNA ITS拷贝中包含有假基因的存在, 而且假基因与功能拷贝之间已经形成了较大的遗传分化, 这暗示假基因起源有较长历史。苏铁核仁组织区不仅多达16个, 而且分布在13条染色体上, 这可能是其nrDNA ITS致同进化不完全的主要原因。     

Divergent and reticulate species relationships in Leucaena (Fabaceae) inferred from multiple data sources: insights into polyploid origins and nrDNA polymorphism

Previous analyses of species relationships and polyploid origins in the mimosoid legume genus Leucaena have used chloroplast DNA (cpDNA) restriction site data and morphology. Here we present an analysis of a new DNA sequence data set for the nuclear ribosomal DNA (nrDNA) 5.8S subunit and flanking ITS 1 and ITS 2 spacers, a simultaneous analysis of the morphology, ITS and cpDNA data sets for the diploid species, and a detailed comparison of the cpDNA and ITS gene trees, which include multiple accessions of all five tetraploid species. Significant new insights into species relationships and polyploid origins, including that of the economically important tropical forage tree L. leucocephala, are discussed. Heterogeneous ITS copy types, including 26 putative pseudogene sequences, were found within individuals of four of the five tetraploid and one diploid species. Potential pseudogenes were identified using two pairwise comparison approaches as well as a tree-based method that compares observed and expected proportions of total ITS variation contributed by the 5.8S subunit optimized onto branches of one of the ITS gene trees. Inclusion of putative pseudogene sequences in the analysis provided evidence that some pseudogenes in allopolyploid L. leucocephala are not the result of post-allopolyploidization gene silencing, but were inherited from its putative diploid maternal progenitor L. pulverulenta.     

ITS1 intra-individual variability of Ascaris isolates from Brazil

The zoonotic potential of Ascaris infecting pigs has stimulated studies of molecular epidemiology with internal transcribed spacer 1 (ITS1) as the target. The aim of this study was to determine the value of Ascaris ITS1 as a molecular marker through assessing the intra-individual genetic diversity of Ascaris isolates from two geographical areas of Brazil. DNA was extracted from single isolated eggs, ITS1 PCR was performed, and the PCR products were cloned and sequenced. Clone analysis showed high ITS1 intra-individual variability revealed by 2–4 ITS1 genotypes/haplotypes per sample (egg). Two genotypes, G1 and G6, and 13 new haplotypes were detected and characterized. The most prevalent in humans, G1 and/or the Brazilian G6, were detected in all samples. Except for genotype G1, no relationship was observed between Brazilian ITS1 genotypes/haplotypes and those previously described in China, Bangladesh, Japan, United Kingdom, Australia, and Denmark, with respect to geographic origin or host affiliation. However, an association between the two geographically separated Brazilian ITS1 isolates was observed. The ITS1 intra-individual variability revealed in this study indicated that the use of this genetic region to discriminate human and pig Ascaris genotypes should be reconsidered.     

Intra-individual internal transcribed spacer 1 (ITS1) and ITS2 ribosomal sequence variation linked with multiple rDNA loci: A case of triploid Atractolytocestus huronensis, the monozoic cestode of common carp

Complete sequences of the ribosomal internal transcribed spacers (ITS1 and ITS2) and karyological characters of the monozoic (unsegmented) tapeworm Atractolytocestus huronensis Anthony, 1958 (Cestoda: Caryophyllidea) from Slovakia were analysed, revealing considerable intra-genomic variability and triploidy in all analysed specimens. Analysis of 20 sequences of each ITS1 and ITS2 spacer yielded eight and 10 different sequence types, respectively. In individual tapeworms, two to four ITS1 and three to four ITS2 sequence types were found. Divergent intra-genomic ITS copies were mostly induced by nucleotide substitutions and different numbers of short repetitive motifs within the sequence. In addition, triploidy was found to be a common feature of A. huronensis. The karyotype of Slovakian A. huronensis possesses three sets of chromosomes (3n = 24, n = 4m + 3st + 1 minute chromosome), similar to the previously described triploidy in conspecific tapeworms from North America. Fluorescent in situ hybridisation (FISH) with a ssrDNA probe revealed two distinct rDNA clusters for each homologue of the triplet number 2. To date, A. huronensis is the only cestode species in which intra-individual ITS sequence variants were found in parallel with its triploid nature and multiple rDNA loci. Some of these molecular and genetic features were observed in several other species of basal or nearly basal tapeworms of the orders Caryophyllidea and Diphyllobothriidea, which indicates that the phenomena may be characteristic for evolutionarily lower tapeworms and deserve more attention in future studies.     

Multiple hybrid formation in natural populations: concerted evolution of the internal transcribed spacer of nuclear ribosomal DNA (ITS) in North American Arabis divaricarpa (Brassicaceae)

DNA sequence variation of the internal transcribed spacer (ITS) region of nuclear ribosomal DNA from Arabis holboellii, A. drummondii, and its putative hybrid A. divaricarpa was analyzed to study hybrid speciation in a species system geographically covering nearly the entire North American continent. Based on molecular systematics the investigated species are better combined under the genus Boechera. Multiple intraindividual ITS copies were detected in numerous accessions of A. divaricarpa, and, to a minor extent, in the parental taxa. Comparative phylogenetic analysis demonstrates that reticulate evolution is common. Consequently, concerted evolution of ITS regions resulted in different types of ITS fragments not only in hybrid populations but also in one of the parental taxa, A. holboellii. Hybrid formation often occurred independently at different sites and at different times, which is reflected by ITS copies resampling the original parental sequence variation in different ways. Some biogeographic structuring of genetic diversity is apparent and mirrors postglacial migration routes. Hybridization, reticulation, and apomixis are assumed to be the major forces driving speciation processes in this species complex. Analysis of conserved regions and secondary structures of the ITS region provided no evidence that, in this system, hybrid ITS evolution is predominantly driven in a particular direction. However, two regions in the ITS1 and ITS2, respectively, show higher mutation rates than expected from outgroup comparisons. Strong evidence for the occurrence of apomixis in A. holboellii and A. divaricarpa has come from pollen size measurements and estimations of pollen quality, which favor the hypothesis that A. drummondii served as paternal hybridization partner more frequently than A. holboellii.     

Within- and between-individual sequence variation among ITS1 copies in the meadow grasshopper Chorthippus parallelus indicates frequent intrachromosomal gene conversion

Sequencing multiple copies of the ITS1 region revealed the coexistence of two or more haplotypes within the genome of Chorthippus parallelus. Using a PCR-RFLP approach, the ITS1 numbers and frequencies of haplotypes present in each of 40 individuals were investigated, revealing a consistent lack of homogeneity. For each individual, the level of intra-individual variation was estimated from a sample of 20 ITS1 copies. The level of differentiation in haplotype frequency among individuals was then estimated by maximum likelihood using models based on the Dirichlet distribution. This confirmed the existence of significant levels of variation among individuals within each population studied. The most likely turnover mechanism that could generate this pattern of variation is gene conversion, operating at the intrachromosomal level. Furthermore, the discovery of linkage disequilibrium among the ITS1 haplotypes of C. parallelus suggests that intrachromosomal gene conversion occurs more frequently than interchromosomal recombination. Subspecies of C. parallelus showed significantly different haplotype distributions following about 0.5 Myr of divergence. With respect to the process of concerted evolution, we show that homogenization of repeats is slow relative to speciation, and the standing variation among individuals is sufficient for selection to operate.     

Heterogeneity of three molecular data partition phylogenies of mints related to M. x piperita (Mentha; Lamiaceae)

Phylogenetic reconstructions with molecular tools are now widely used, thanks to advances in PCR and sequencing technologies. The choice of the molecular target still remains a problem because too few comparative data are available. This is particularly true for hybrid taxa, where differential introgression of genome parts leads to incongruity between data sets. We have studied the potential of three data partitions to reconstruct the phylogeny of mints related to M. x piperita. These included nuclear DNA (ITS), chloroplast DNA (non-coding regions trnL intron, intergenic spacers trnL-trnF, and psbA-trnH), and AFLP and ISSR, markers. The taxonomic sampling was composed of hybrids, diploid and polyploid genomes. Since the genealogy of cultivated mint hybrids is known, they represent a model group to compare the usefulness of various molecular markers for phylogeny inference. Incongruities between ITS, chloroplast DNA, and AFLP-ISSR phylogenetic trees were recorded, although DNA fingerprinting data were congruent with morphological classification. Evidence of chloroplast capture events was obtained for M. x piperita. Direct sequencing of ITS led to biased results because of the existence of pseudogenes. Sequencing of cloned ITS further failed to provide evidence of the existence of the two parental copy types for M. x piperita, a sterile hybrid that has had no opportunity for concerted evolution of ITS copies. AFLP-ISSR data clustered M. x piperita with the parent that had the largest genome. This study sheds light on differential of introgression of different genome regions in mint hybrids.     

Variation in the ribosomal internal transcribed spacers and 5.8S rDNA among five species of Acropora (Cnidaria; Scleractinia): patterns of variation consistent with reticulate evolution

The ITS sequences of Acropora spp. are the shortest so far identified in any metazoan and are among the shortest seen in eukaryotes; ITS1 was 70-80 bases, and ITS2 was 100-112 bases. The ITS sequences were also highly variable, but base composition and secondary structure prediction indicate that divergent sequence variants are unlikely to be pseudogenes. The pattern of variation was unusual in several other respects: (1) two distinct ITS2 types were detected in both A. hyacinthus and A. cytherea, species known to hybridize in vitro with high success rates, and a putative intermediate ITS2 form was also detected in A. cytherea; (2) A. valida was found to contain highly (29%) diverged ITS1 variants; and (3) A. longicyathus contained two distinct 5.8S rDNA types. These data are consistent with a reticulate evolutionary history for the genus Acropora.