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.     

Multiple origins and nrDNA internal transcribed spacer homeologue evolution in the Glycine tomentella (Leguminosae) allopolyploid complex

Despite the importance of polyploidy in the evolution of plants, patterns of molecular evolution and genomic interactions following polyploidy are not well understood. Nuclear ribosomal DNA is particularly complex with respect to these genomic interactions. The composition of nrDNA tandem arrays is influenced by intra- and interlocus concerted evolution and their expression is characterized by patterns such as nucleolar dominance. To understand these complex interactions it is important to study them in diverse natural polyploid systems. In this study we use direct sequencing to isolate and characterize nrDNA internal transcribed spacer (ITS) homeologues from multiple accessions of six different races in the Glycine tomentella allopolyploid complex. The results indicate that in most allopolyploid accessions both homeologous nrDNA repeats are present, but that there are significant biases in copy number toward one homeologue, possibly resulting from interlocus concerted evolution. The predominant homeologue often differs between races and between accessions within a race. A phylogenetic analysis of ITS sequences provides evidence for multiple origins in several of the polyploid races. This evidence for diverse patterns of nrDNA molecular evolution and multiple origins of polyploid races will provide a useful system for future studies of natural variation in patterns of nrDNA expression.     

ITS Polymorphisms Shed Light on Hybrid Evolution in Apomictic Plants: A Case Study on the Ranunculus auricomus Complex

The reconstruction of reticulate evolutionary histories in plants is still a major methodological challenge. Sequences of the ITS nrDNA are a popular marker to analyze hybrid relationships, but variation of this multicopy spacer region is affected by concerted evolution, high intraindividual polymorphism, and shifts in mode of reproduction. The relevance of changes in secondary structure is still under dispute. We aim to shed light on the extent of polymorphism within and between sexual species and their putative natural as well as synthetic hybrid derivatives in the Ranunculus auricomus complex to test morphology-based hypotheses of hybrid origin and parentage of taxa. We employed direct sequencing of ITS nrDNA from 68 individuals representing three sexuals, their synthetic hybrids and one sympatric natural apomict, as well as cloning of ITS copies in four representative individuals, RNA secondary structure analysis, and landmark geometric morphometric analysis on leaves. Phylogenetic network analyses indicate additivity of parental ITS variants in both synthetic and natural hybrids. The triploid synthetic hybrids are genetically much closer to their maternal progenitors, probably due to ploidy dosage effects, although exhibiting a paternal-like leaf morphology. The natural hybrids are genetically and morphologically closer to the putative paternal progenitor species. Secondary structures of ITS1-5.8S-ITS2 were rather conserved in all taxa. The observed similarities in ITS polymorphisms suggest that the natural apomict R. variabilis is an ancient hybrid of the diploid sexual species R. notabilis and the sexual species R. cassubicifolius. The additivity pattern shared by R. variabilis and the synthetic hybrids supports an evolutionary and biogeographical scenario that R. variabilis originated from ancient hybridization. Concerted evolution of ITS copies in R. variabilis is incomplete, probably due to a shift to asexual reproduction. Under the condition of comprehensive inter- and intraspecific sampling, ITS polymorphisms are powerful for elucidating reticulate evolutionary histories.     

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

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

ITS假基因对栎属系统学研究的影响及其对分子系统学研究的启示

核糖体rDNA ITS是被子植物系统发育研究中应用最广泛的分子标记之一。以前人们认为同一物种中的ITS序列因致同进化而使不同拷贝高度一致,在分子系统学研究中常以ITS1-5.8S-ITS2序列作为构建系统进化树的基础。近年来,在对一些被子植物的研究中发现这段序列在同一物种中具有多态性,有些拷贝中的5．8S区不具编码功能,人们把含有不具编码功能5．8S区的ITS1-5．8S-ITS2序列定义为ITS假基因序列,它对同源基因致同进化的假设形成了新的挑战。在诸多应用ITS序列重建系统进化关系的研究中,栎属系统学研究因ITS假基因的发现而倍受关注。本文以栎属为例回顾了ITS假基因的发现过程,分析了其对该属系统学研究的影响,为分子生物学在植物系统进化研究中的应用提供一些新的参考。     

Persistent nuclear ribosomal DNA sequence polymorphism in the Amelanchier agamic complex (Rosaceae)

Individual plants of several Amelanchier taxa contain many polymorphic nucleotide sites in the internal transcribed spacers (ITS) of nuclear ribosomal DNA (nrDNA). This polymorphism is unusual because it is not recent in origin and thus has resisted homogenization by concerted evolution. Amelanchier ITS sequence polymorphism is hypothesized to be the result of gene flow between two major North American clades resolved by phylogenetic analysis of ITS sequences. Western North American species plus A. humilis and A. sanguinea of eastern North America form one clade (A), and the remaining eastern North American Amelanchier make up clade B. Five eastern North American taxa are polymorphic at many of the nucleotide sites where clades A and B have diverged and are thought to be of hybrid origin, with A. humilis or A. sanguinea as one parent and various members of clade B as the other parent. Morphological evidence suggests that A. humilis is one of the parents of one of the polymorphic taxa, a microspecies that we refer to informally as A. "erecta." Sequences of 21 cloned copies of the ITS1- 5.8S gene-ITS2 region from one A. "erecta" individual are identical to A. humilis sequence or to the clade B consensus sequence, or they are apparent recombinants of A. humilis and clade B ITS repeats. Amelanchier "erecta" and another polymorphic taxon are suspected to be relatively old because both grow several hundred kilometers beyond the range of one of their parents. ITS sequence polymorphisms have apparently persisted in these two taxa perhaps because of polyploidy and/or agamospermy (asexual seed production), which are prevalent in the genus.      

Extensive pyrosequencing reveals frequent intra-genomic variations of internal transcribed spacer regions of nuclear ribosomal DNA

Background

Internal transcribed spacer of nuclear ribosomal DNA (nrDNA) is already one of the most popular phylogenetic and DNA barcoding markers. However, the existence of its multiple copies has complicated such usage and a detailed characterization of intra-genomic variations is critical to address such concerns.

Methodology/Principal Findings

In this study, we used sequence-tagged pyrosequencing and genome-wide analyses to characterize intra-genomic variations of internal transcribed spacer 2 (ITS2) regions from 178 plant species. We discovered that mutation of ITS2 is frequent, with a mean of 35 variants per species. And on average, three of the most abundant variants make up 91% of all ITS2 copies. Moreover, we found different congeneric species share identical variants in 13 genera. Interestingly, different species across different genera also share identical variants. In particular, one minor variant of ITS2 in Eleutherococcus giraldii was found identical to the ITS2 major variant of Panax ginseng, both from Araliaceae family. In addition, DNA barcoding gap analysis showed that the intra-genomic distances were markedly smaller than those of the intra-specific or inter-specific variants. When each of 5543 variants were examined for its species discrimination efficiency, a 97% success rate was obtained at the species level.

Conclusions

Identification of identical ITS2 variants across intra-generic or inter-generic species revealed complex species evolutionary history, possibly, horizontal gene transfer and ancestral hybridization. Although intra-genomic multiple variants are frequently found within each genome, the usage of the major variants alone is sufficient for phylogeny construction and species determination in most cases. Furthermore, the inclusion of minor variants further improves the resolution of species identification.     

Molecular systematics,GISH and the origin of hybrid taxa in Nicotiana (Solanaceae)

Phylogenetic relationships in the genus Nicotiana were investigated using parsimony analyses of the internal transcribed spacer (ITS) regions of nuclear ribosomal DNA (nrDNA). In addition, origins of some amphidiploid taxa in Nicotiana were investigated using the techniques of genomic in situ hybridization (GISH), and the results of both sets of analyses were used to evaluate previous hypotheses about the origins of these taxa. Phylogenetic analyses of the ITS nrDNA data were performed on the entire genus (66 of 77 naturally occurring species, plus three artificial hybrids), comprising both diploid and polyploid taxa, and on the diploid taxa only (35 species) to examine the effects of amphidiploids on estimates of relationships. All taxa, regardless of ploidy, produced clean, single copies of the ITS region, even though some taxa are hybrids. Results are compared with a published plastid (matK) phylogeny using fewer, but many of the same, taxa. The patterns of relationships in Nicotiana, as seen in both analyses, are largely congruent with each other and previous evolutionary ideas based on morphology and cytology, but some important differences are apparent. None of the currently recognized subgenera of Nicotiana is monophyletic and, although most of the currently recognized sections are coherent, others are clearly polyphyletic. Relying solely upon ITS nrDNA analysis to reveal phylogenetic patterns in a complex genus such as Nicotiana is insufficient, and it is clear that conventional analysis of single data sets, such as ITS, is likely to be misleading in at least some respects about evolutionary history. ITS sequences of natural and well-documented amphidiploids are similar or identical to one of their two parents-usually, but not always, the maternal parent-and are not in any sense themselves 'hybrid'. Knowing how ITS evolves in artificial amphidiploids gives insight into what ITS analysis might reveal about naturally occurring amphidiploids of unknown origin, and it is in this perspective that analysis of ITS sequences is highly informative.     

Ribosomal DNA variation,recombination and inheritance in the basidiomycete Trichaptum abietinum: implications for reticulate evolution

Two divergent nuclear ribosomal DNA (nrDNA) types, designated alpha and beta, were found distributed in 11 North European populations of the basidiomycete Trichaptum abietinum. These types differed by a 220 bp indel in the internal transcribed spacer 1 (ITS1) sequence and a number of linked substitutions and small indel motives in the internal transcribed and intergenic spacers (ITS1, ITS2, IGS1 and IGS2). The alpha and beta haplotypes co-occurred in heterozygous somatic individuals (dikaryons) and segregated in a Mendelian fashion in monokaryotic single spore progenies. This result suggests that the haplotypes are encoded in different nuclei of field-collected dikaryons and inherited as a single locus. No meiotic recombinants were observed among the sequenced monokaryons. Population genetic analyses by PCR-RFLP revealed that a low frequency of evolutionary intermediate nrDNA types also existed in natural populations, presumably as a result of meiotic recombination of alpha and beta nrDNA. The existence of divergent nrDNA types in T. abietinum could be a result of a former independent evolution followed by a hybridization event. Phylogenetic analyses of ITS sequences suggest that the sister taxon T. fusco-violaceum has been involved in the evolutionary history of T. abietinum. Sequence polymorphisms observed in the translation elongation factor 1alpha (efa) and glyceraldehyde-3-phosphate dehydrogenase (gpd) genes, did not reveal two well-defined types of these genes. The results are discussed in the light of other evolutionary mechanisms as well.     

Ribosomal DNA polymorphisms in the yeast Geotrichum candidum

The dimorphic yeast Geotrichum candidum (teleomorph: Galactomyces candidus) is commonly used to inoculate washed-rind and bloomy-rind cheeses. However, little is known about the phylogenetic lineage of this microorganism. We have sequenced the complete 18S, 5.8S, 26S ribosomal RNA genes and their internal transcribed spacers (ITS1) and ITS2 regions (5126 nucleotides) from 18 G. candidum strains from various environmental niches, with a focus on dairy strains. Multiple sequence alignments revealed the presence of 60 polymorphic sites, which is generally unusual for ribosomal DNA (rDNA) within a given species because of the concerted evolution mechanism. This mechanism drives genetic homogenization to prevent the divergent evolution of rDNA copies within individuals. While the polymorphisms observed were mainly substitutions, one insertion/deletion (indel) polymorphism was detected in ITS1. No polymorphic sites were detected downstream from this indel site, that is, in 5.8S and ITS2. More surprisingly, many sequence electrophoregrams generated during the sequencing of the rDNA had dual peaks, suggesting that many individuals exhibited intragenomic rDNA variability. The ITS1-5.8S-ITS2 regions of four strains were cloned. The sequence analysis of 68 clones revealed 32 different ITS1-5.8S-ITS2 variants within these four strains. Depending on the strain, from four to twelve variants were detected, indicating that multiple rDNA copies were present in the genomes of these G. candidum strains. These results contribute to the debate concerning the use of the ITS region for barcoding fungi and suggest that community profiling techniques based on rDNA should be used with caution.     

Non-concerted ITS evolution, early origin and phylogenetic utility of ITS pseudogenes in Pyrus

Molecular studies of 19 species of the genus Pyrus L. revealed different degrees of intra-individual polymorphism of the internal transcribed spacer (ITS 1, 5.8S rDNA and ITS 2) region due to the existence of putative non-functional copies (pseudogenes), putative recombinants and non-concerted evolution among functional copies. Different types of ITS pseudogenes displaying lower GC content and unstable secondary structure were preferentially amplified under normal PCR conditions. Functional ITS copies were successfully obtained in all investigated accessions under the modified PCR conditions. All pseudogenes were highly divergent from their corresponding functional copies and formed a monophyletic group in the phylogenetic tree based on all paralogs, indicating they were of relatively early origin. Functional ITS copies led to confused and poorly resolved phylogeny as a result of low sequence divergence, existence of unidentified ancient recombinants and a high degree of intra-individual functional ITS polymorphism, while certain types of pseudogenes and some relict pseudogenes offered more credible clues for the evolutionary history of Pyrus species.     

Phylogenetic relationships and biogeographic patterns in North American members of Carex section Acrocystis (Cyperaceae) using nrDNA ITS and ETS sequence data

Carex section Acrocystis currently includes 27 taxa in North America. Recent phylogenetic studies have suggested that the North American and some but not all of the Eurasian species form a clade. Relationships and biogeographic patterns among species in this core-Acrocystis group are explored here using nuclear ribosomal (nrDNA) internal transcribed spacer region (ITS) and nrDNA external transcribed spacer region (ETS) sequence data. While maximum parsimony analysis of the ITS and ETS data provides only a moderately resolved branching structure for species relationships within the core-Acrocystis clade, maximum likelihood analysis provides a more resolved hypothesis of relationships in the section. The core-Acrocystis clade consists of a grade of Eurasian and primarily western North American species, with a well-supported clade of only eastern North American species nested within this grade. ITS and ETS types do not coalesce within many species or species complexes. Possible explanations for the non-coalescent nature of ITS and ETS copies in Acrocystis are explored, including lineage sorting, hybridization, and cryptic species.     

BIOLUMINESCENCE IN THE MARINE DINOFLAGELLATE PYROCYSTIS LUNULA

Intragenomic nrDNA ITS variability was examined in representative species of the Stephanodiscus niagarae complex. PCR was performed on DNA extracted from monoclonal cultures, and sequence variability determined using standard cloning techniques. Preliminary data from all species reveal intragenomic polymorphism. ITS sequence data is being used to delineate closely related species of algae and polymorphisms are potentially problematic for phylogenetic reconstruction and identity. We are exploring the extent of polymorphism in other freshwater Thalassiorsiraceae, including Stephanodiscus hantzschii fo. tenuis and Cyclostephanos invisitatus.     

Analysis of nrDNA polymorphism in closely related diploid sexual, tetraploid sexual and polyploid agamospermous species

Nuclear sequences of ITS1-5.8S-ITS2 region of rDNA may be an important source of phylogenetically informative data provided that nrDNA is cloned and the character of sequence variation of clones is properly analyzed. nrDNA of selected Taraxacum sections was studied to show sequence variation differences among diploid sexual, tetraploid sexual and polyploid agamospermous species. We examined nucleotide characteristics, substitution pattern, secondary structure, and the phylogenetic utility of ITS1-5.8S-ITS2 from 301 clones of 32 species representing 11 sections. The most divergent sequences of ITS1&2 differed by 17.1% and in 5.8S only by 3.7%. The ITS1-5.8S-ITS2 characteristics, integrity and also stability of secondary structures confirmed that pseudogenes are not responsible for the above variation. The within-individual polymorphism of clones implies that the concerted evolution of ITS cistron of agamospermous polyploid Taraxacum is remarkably suppressed. Sequences of ITS clones proved to be a useful tool for mapping pathways of complex reticulation (polyploid hybridity) in agamospermous Taraxacum.     

ITS sequence variation and concerted evolution in the natural hybrid species Malus toringoides

The internal transcribed spacer (ITS) region of nuclear ribosomal DNA (nrDNA) is one of the most used molecular characters in plant systematics. Our previous studies based on morphological analysis and ITS sequence variation suggested that Malus toringoides (Rehd.) Hughes is derived from hybridization between M. transitoria (Batal.) Schneid. and M. kansuensis (Batal.) Schneid. To further understand the variation pattern of ITS sequences in M. toringoides, and to elucidate the evolutionary processes that affect ITS sequence variation after hybridization, we sampled 99 accessions from multiple populations of the hybrid and parental species, and then obtained totally 254 ITS sequences by cloning and sequencing. Our ITS variation data demonstrates three outcomes of ITS repeats after hybrid speciation. ～ 27–41% of M. toringoides have only M. transitoria type ITS sequence, ～ 40–70% have M. transitoria type ITS sequence plus one or two chimeric ITS sequences generated by recombination between parental ITS sequences, and six accessions retain both parental type ITS sequences. The plausible evolutionary processes that created the observed ITS variations were inferred to be the joint actions of recombination, concerted evolution, pseudogenization and backcrossing. Our study provides further understandings of the variation model of ITS repeats after hybridization as well as the evolution of M. toringoides after its hybrid speciation.     

Comparative ITS and AFLP analysis of diploid Cardamine (Brassicaceae) taxa from closely related polyploid complexes

BACKGROUND AND AIMS: Diploid representatives from the related polyploid complexes of Cardamine amara, C. pratensis and C. raphanifolia (Brassicaceae), were studied to elucidate phylogenetic relationships among the complexes and among the individual taxa included. METHODS: Two independent molecular data sets were used: nucleotide sequences from the internal transcribed spacers (ITS) of nrDNA, and amplified fragment length polymorphism (AFLP) markers. Seventeen diploid taxa from the studied groups were sampled. KEY RESULTS: Both ITS and AFLP analyses provided congruent results in inferred relationships, and revealed two main lineages. While the C. amara group, consisting of C. wiedemanniana and four subspecies of C. amara, was resolved as a well-supported monophyletic group, taxa from the C. pratensis and C. tenera groups (the latter representing diploid taxa of the complex of C. raphanifolia) all appeared together in a single clade/cluster with no support for the recognition of either of the groups. Intra-individual polymorphisms and patterns of nucleotide variation in the ITS region in C. uliginosa and C. tenera, together with the distribution of AFLP bands, indicate ancient hybridization and introgression among these Caucasian diploids. CONCLUSIONS: The lack of supported hierarchical structure suggests that extensive reticulate evolution between these groups, even at the diploid level, has occurred (although an alternative explanation, namely ancestral polymorphism in ITS data, cannot be completely excluded). Several implications for the investigation of the polyploid complexes of concern are drawn. When tracing origins of polyploid taxa, a much more complex scenario should be expected, taking into account all relatives as potential parents, irrespective of the group in which they are classified.     

粉拟青霉种内nrDNA ITS分析

通过对20株粉拟青霉(Paecilomyces farinosus)ITS1-5．8S-ITS2(rDNA)区域序列测定,确定了粉拟青霉ITS序列,而韩国学者测定的粉拟青霉ITS序列应为细脚拟青霉(P．tenuipes)的序列。序列比较发现,韩国2株未定名的拟青霉(Paecilomyces spp．)菌株(KACC40219、KACC40221)应为粉拟青霉。基于本研究构建的邻接树推断,粉拟青霉的有性型可能是一种虫草。粉拟青霉的起源应为单源的。不同的粉拟青霉菌株的ITS序列具有多态性,源于同一地区的菌株的ITS变异也较大。ITS序列的证据表明,粉拟青霉菌株间的差异与地理来源及寄主均无相关性。     

Evolution of the Euphrasia transmorrisonensis complex (Orobanchaceae) in alpine areas of Taiwan

Aims To unravel isolation and differentiation of the genetic structure of the Euphrasia transmorrisonensis complex, a showy herb, among alpine regions of mountain peaks in subtropical Taiwan and to infer its evolutionary history. Location Alpine ecosystems of high‐montane regions of Taiwan. Methods Phylogenetic analyses of the trnL intron and the trnL–trnF intergenic spacer of chloroplast (cp) DNA, and the intertranscribed spacer (ITS) of nuclear ribosomal (nr) DNA between 18S and 26S were carried out on 18 populations of the E. transmorrisonensis complex in Taiwan. Results In total, 10 haplotypes for cpDNA and 14 haplotypes for nrDNA were detected. Three population groups located in the northern, north‐eastern, and south‐central regions of the Central Mountain Range (CMR) were revealed according to the frequencies of haplotypes and haplotype lineages of nrDNA. Balancing selection might have played a role in the evolution of Euphrasia in Taiwan. Main conclusions By integrating the spatial‐genetic patterns of cpDNA and nrDNA, two possible evolutionary histories of Euphrasia in Taiwan were inferred. The favourable hypotheses for interpreting the data suggest at least three origins of the E. transmorrisonensis complex in Taiwan, corresponding to each nuclear lineage in the northern (II), northern/north‐eastern (I), and central/southern regions (III) with subsequent hybridization between lineages I and II and lineages II and III. These lineage boundaries are strengthened by the finding that haplotypes of C derived from cpDNA were found in the geographical region of lineage II of nrDNA, while haplotypes of A derived from cpDNA were found in the region of lineage III of nrDNA. Thus, the origin of chloroplasts exclusive to lineages II and III supports their long‐term isolation from one another.     

Molecular authentication of ginseng cultivars by comparison of internal transcribed spacer and 5.8S rDNA sequences

Molecular authentication among three Panax species and within cultivars and accessions of P. ginseng was investigated using the DNA sequence in the ribosomal ITS1–5.8S–ITS2 region. Four single-nucleotide polymorphisms were identified between P. ginseng and other Panax species. In the electrophoresis profile, obtained after digestion with the enzyme TaqI, three fingerprinting patterns were obtained from cultivars and accessions of Panax species. Consequently, this authentication procedure based upon the restriction fragment length polymorphism in the ribosomal ITS1–5.8S–ITS2 region can now be utilized to differentiate these Panax species as well as major Korean cultivars such as Gopoong and Kumpoong from other cultivars and accessions in Panax species at the DNA level. O. T. Kim and K. H. Bang contributed equally to this paper.