The internal transcribed spacer of nuclear ribosomal DNA in the gymnosperm Gnetum

Integration of large and diverse angiosperm DNA fragments into Asian Gnetum mitogenomes

BACKGROUND

Horizontal gene transfer (HGT) events have rarely been reported in gymnosperms. Gnetum is a gymnosperm genus comprising 25‒35 species sympatric with angiosperms in West African, South American, and Southeast Asian rainforests. Only a single acquisition of an angiosperm mitochondrial intron has been documented to date in Asian Gnetum mitogenomes. We wanted to develop a more comprehensive understanding of frequency and fragment length distribution of such events as well as their evolutionary history in this genus.

RESULTS

We sequenced and assembled mitogenomes from five Asian Gnetum species. These genomes vary remarkably in size and foreign DNA content. We identified 15 mitochondrion-derived and five plastid-derived (MTPT) foreign genes. Our phylogenetic analyses strongly indicate that these foreign genes were transferred from diverse eudicots-mostly from the Rubiaceae genus Coptosapelta and ten genera of Malpighiales. This indicates that Asian Gnetum has experienced multiple independent HGT events. Patterns of sequence evolution strongly suggest DNA-mediated transfer between mitochondria as the primary mechanism giving rise to these HGT events. Most Asian Gnetum species are lianas and often entwined with sympatric angiosperms. We therefore propose that close apposition of Gnetum and angiosperm stems presents opportunities for interspecific cell-to-cell contact through friction and wounding, leading to HGT.

CONCLUSIONS

Our study reveals that multiple HGT events have resulted in massive amounts of angiosperm mitochondrial DNA integrated into Asian Gnetum mitogenomes. Gnetum and its neighboring angiosperms are often entwined with each other, possibly accounting for frequent HGT between these two phylogenetically remote lineages.

The high molecular diversity in Noctiluca scintillans is dominated by intra-genomic variations revealed by single cell high-throughput sequencing of 18S rDNA V4

The application of high-throughput sequencing (HTS) technologies has revolutionized research on phytoplankton biodiversity by generating an unprecedented amount of molecular data in marine ecosystem surveys. However, high-level of molecular diversity uncovered in HTS-based metabarcoding analyses may lead to overinterpretation of phytoplankton diversity due to excessive intra-genomic variations (IGVs). The aims in this study are to explore the nature of phytoplankton molecular diversity and to test the hypothesis. We carried out single-cell metabarcoding analysis of 18S rDNA V4 sequences obtained in single Noctiluca scintillans cells isolated from various sites in coastal waters of China. Results showed that each single N. scintillans cell harbored a high level of IGVs with about 100 amplicon sequence variants (ASVs). The large numbers of non-dominant ASVs identified in N. scintillans cells, which might correspond to the larger numbers of ASVs annotated as N. scintillans and showed similar temporal dynamics in metabarcoding analyses, could inflate the inter-species diversity or intra-species genetic diversity. In addition, there were large numbers of additional ASVs that were not annotated as N. scintillans. These non-N. scintillans ASVs might represent diverse preys for N. scintillans, consistent with previous reports that N. scintillans may act as chance predator of a broad-spectrum preys. This single-cell study has unambiguously demonstrated that the existence of high levels of IGVs in N. scintillans and most likely many other phytoplankton species, demonstrating that the majority of the molecular diversity revealed in metabarcoding analysis, which were generally interpreted as the sum of inter-species diversity and intra-species diversity, actually included high levels of IGVs and should be interpreted with caution.

Intragenomic rDNA variation - the product of concerted evolution, mutation, or something in between?

The classical model of concerted evolution states that hundreds to thousands of ribosomal DNA (rDNA) units undergo homogenization, making the multiple copies of the individual units more uniform across the genome than would be expected given mutation frequencies and gene redundancy. While the universality of this over 50-year-old model has been confirmed in a range of organisms, advanced high throughput sequencing techniques have also revealed that rDNA homogenization in many organisms is partial and, in rare cases, even apparently failing. The potential underpinning processes leading to unexpected intragenomic variation have been discussed in a number of studies, but a comprehensive understanding remains to be determined. In this work, we summarize information on variation or polymorphisms in rDNAs across a wide range of taxa amongst animals, fungi, plants, and protists. We discuss the definition and description of concerted evolution and describe whether incomplete concerted evolution of rDNAs predominantly affects coding or non-coding regions of rDNA units and if it leads to the formation of pseudogenes or not. We also discuss the factors contributing to rDNA variation, such as interspecific hybridization, meiotic cycles, rDNA expression status, genome size, and the activity of effector genes involved in genetic recombination, epigenetic modifications, and DNA editing. Finally, we argue that a combination of approaches is needed to target genetic and epigenetic phenomena influencing incomplete concerted evolution, to give a comprehensive understanding of the evolution and functional consequences of intragenomic variation in rDNA.

Lack of ITS sequence homogenization in Erysimum species (Brassicaceae) with different ploidy levels

The internal transcribed spacers (ITS) exhibit concerted evolution by the fast homogenization of these sequences at the intragenomic level. However, the rate and extension of this process are unclear and might be conditioned by the number and divergence of the different ITS copies. In some cases, such as hybrid species and polyploids, ITS sequence homogenization appears incomplete, resulting in multiple haplotypes within the same organism. Here, we studied the dynamics of concerted evolution in 85 individuals of seven plant species of the genus Erysimum (Brassicaceae) with multiple ploidy levels. We estimated the rate of concerted evolution and the degree of sequence homogenization separately for ITS1 and ITS2 and whether these varied with ploidy. Our results showed incomplete sequence homogenization, especially for polyploid samples, indicating a lack of concerted evolution in these taxa. Homogenization was usually higher in ITS2 than in ITS1, suggesting that concerted evolution operates more efficiently on the former. Furthermore, the hybrid origin of several species appears to contribute to the maintenance of high haplotype diversity, regardless of the level of ploidy. These findings indicate that sequence homogenization of ITS is a dynamic and complex process that might result in varying intra- and inter-genomic diversity levels.

Extensive nrDNA ITS polymorphism in Lycium: Non-concerted evolution and the identification of pseudogenes

The internal transcribed spacer (ITS) is one of the most extensively sequenced molecular markers in plant systematics due to its generally concerted evolution. While non-concerted evolution has been found in some plant taxa, such information is missing in Lycium. Molecular studies of six species and two variants of the genus Lycium revealed high levels of intra- and inter-individual polymorphism in the ITS, indicating non-concerted evolution. All genomic DNA ITS paralogues were identified as putative pseudogenes or functional paralogues through a series of comparisons of sequence features, including length and substitution variation, GC content, secondary structure stability, and the presence of conserved motifs in the 5.8S gene, and the rate of evolution. Approximately, 60% of ITS pseudogenes could be easily detected. Based on phylogenetic analysis, all pseudogenes were highly distinct from their corresponding functional copies, tended to evolve neutrally, and clustered randomly together in the evolutionary tree. The results probably suggest that this ITS non-concerted evolution is related to the recent divergence between tandem repeats within the Lycium genome and hybridization between species. Our study complements those of pseudogenes in plant taxa and provides a theoretical basis for the phylogeny and genetic origin of the genus Lycium while having important implications for the use of ITS molecular markers for phylogenetic reconstruction.

Phylogenetic Analysis and Molecular Diversity of Capsicum Based on rDNA-ITS Region

The genus Capsicum is comprised of 5 domesticated and more than 30 wild species. The region of nuclear ribosomal DNA internal transcribed spacers (rDNA-ITS) has widely been used for species identification, but has rarely been used in Capsicum. In this study, the evaluation of genetic diversity and a phylogenetic analysis were conducted using rDNA-ITS of 28 Capsicum accessions, including five domesticated and two wild species. We surveyed six conventional keys of domesticated species and another five traits in Capsicum accessions. Specific morphological characteristics were found in C. annuum, C. baccatum, and C.pubescens. Three subclones of each accession were sequenced, and rDNA-ITS polymorphisms were detected in all accessions excluding C. annuum, suggesting that incomplete concerted evolution occurred in rDNA-ITS of Capsicum. The genetic diversity was evaluated using nucleotide polymorphism and diversity. C. annuum had the lowest genetic diversity of all species in this study. The phylogenetic tree formed a species-specific clade for C. annuum, C. baccatum, and C. pubescens. The C. chinense clade existed in the C. frutescens clade, implying that it was a cultivated variant of C. frutescens. C. chacoense likely belonged to the C. baccatum complex according to its morphologic and genetic features. This study indicated that the rDNA-ITS region can be used for simple identification of domesticated Capsicum species.

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.

Gnetum chinense, a new species of Gnetaceae from southwestern China

Gnetum chinense sp. nov., a new lianoid species of Gnetaceae, is described from southwestern China. The new species is morphologically similar to G. montanum Markgr. in its oblong elliptic leaves and the ovoid to ellipsoid chlamydosperm, but differs from the latter by its shorter male spikes having fewer involucral collars (7-10 vs. 13-18 in G. montanum). We also did a new molecular analysis using one nuclear marker (i.e. nrITS) and four chloroplast markers (i.e. matK gene, rpoC1 intron, psbB-rps12 IGS, and trnF-trnV IGS). The result suggests that this specific clade is sister to a large clade consisting of all other known Chinese lianoid species of Gnetum except G. parvifolium (Warb.) W.C. Cheng.

Molecular phylogeny of mulberries reconstructed from ITS and two cpDNA sequences

Background

Species in the genus Morus (Moraceae) are deciduous woody plants of great economic importance. The classification and phylogenetic relationships of Morus, especially the abundant mulberry resources in China, is still undetermined. Internal transcribed spacer (ITS) regions are among the most widely used molecular markers in phylogenetic analyses of angiosperms. However, according to the previous phylogenetic analyses of ITS sequences, most of the mulberry accessions collected in China were grouped into the largest clade lacking for phylogenetic resolution. Compared with functional ITS sequences, ITS pseudogenes show higher sequence diversity, so they can provide useful phylogenetic information.

Methods

We sequenced the ITS regions and the chloroplast DNA regions TrnL-TrnF and TrnT-TrnL from 33 mulberry accessions, and performed phylogenetic analyses to explore the evolution of mulberry.

Results

We found ITS pseudogenes in 11 mulberry accessions. In the phylogenetic tree constructed from ITS sequences, clade B was separated into short-type sequence clades (clades 1 and 2), and a long-type sequence clade (clade 3). Pseudogene sequences were separately clustered into two pseudogroups, designated as pseudogroup 1 and pseudogroup 2. The phylogenetic tree generated from cpDNA sequences also separated clade B into two clades.

Conclusions

Two species were separated in clade B. The existence of three connection patterns and incongruent distribution patterns between the phylogenetic trees generated from cpDNA and ITS sequences suggested that the ITS pseudogene sequences connect with genetic information from the female progenitor. Hybridization has played important roles in the evolution of mulberry, resulting in low resolution of the phylogenetic analysis based on ITS sequences. An evolutionary pattern illustrating the evolution history of mulberry is proposed. These findings have significance for the conservation of local mulberry resources. Polyploidy, hybridization, and concerted evolution have all played the roles in the evolution of ITS sequences in mulberry. This study will expand our understanding of mulberry evolution.

The organization of nuclear ribosomal DNA in gnetophytes - physically separate and physically linked arrangements of 35S and 5S genes. A commentary on: 'Remarkable variation of ribosomal DNA organization and copy number in gnetophytes, a distinct lineage of gymnosperms'

This article comments on: Wencai Wang, Tao Wan, Hannes Becher, Alena Kuderova, Ilia J. Leitch, Sonia Garcia, Andrew R. Leitch and Aleš Kovařík. 2019. Remarkable variation of ribosomal DNA organization and copy number in gnetophytes, a distinct lineage of gymnosperms. Annals of Botany 123(5): 767–781.

Remarkable variation of ribosomal DNA organization and copy number in gnetophytes, a distinct lineage of gymnosperms

INTRODUCTION

Gnetophytes, comprising the genera Ephedra, Gnetum and Welwitschia, are an understudied, enigmatic lineage of gymnosperms with a controversial phylogenetic relationship to other seed plants. Here we examined the organization of ribosomal DNA (rDNA) across representative species.

METHODS

We applied high-throughput sequencing approaches to isolate and reconstruct rDNA units and to determine their intragenomic homogeneity. In addition, fluorescent in situ hybridization and Southern blot hybridization techniques were used to reveal the chromosome and genomic organization of rDNA.

KEY RESULTS

The 5S and 35S rRNA genes were separate (S-type) in Gnetum montanum, Gnetum gnemon and Welwitschia mirabilis and linked (L-type) in Ephedra altissima. There was considerable variability in 5S rDNA abundance, ranging from as few as ~4000 (W. mirabilis) to >100 000 (G. montanum) copies. A similar large variation was also observed in 5S rDNA locus numbers (two to 16 sites per diploid cell). 5S rRNA pseudogenes were interspersed between functional genes forming a single unit in E. altissima and G. montanum. Their copy number was comparable or even higher than that of functional 5S rRNA genes. In E. altissima internal transcribed spacers of 35S rDNA were long and intrinsically repetitive while in G. montanum and W. mirabilis they were short without the subrepeats.

CONCLUSIONS

Gnetophytes are distinct from other gymnosperms and angiosperms as they display surprisingly large variability in rDNA organization and rDNA copy and locus numbers between genera, with no relationship between copy numbers and genome sizes apparent. Concerted evolution of 5S rDNA units seems to have led to the amplification of 5S pseudogenes in both G. montanum and E. altissima. Evolutionary patterns of rDNA show both gymnosperm and angiosperm features underlining the diversity of the group.

Sequence-based nomenclature: a reply to Thines et al. and Zamora et al. and provisions for an amended proposal "from the floor" to allow DNA sequences as types of names

We reply to two recently published, multi-authored opinion papers by opponents of sequence-based nomenclature, namely Zamora et al. (IMA Fungus9: 167-175,2018) and Thines et al. (IMA Fungus9: 177-183, 2018). While we agree with some of the principal arguments brought forward by these authors, we address misconceptions and demonstrate that some of the presumed evidence presented in these papers has been wrongly interpreted. We disagree that allowing sequences as types would fundamentally alter the nature of types, since a similar nature of abstracted features as type is already allowed in the Code (Art. 40.5), namely an illustration. We also disagree that there is a high risk of introducing artifactual taxa, as this risk can be quantified at well below 5 %, considering the various types of high-throughput sequencing errors. Contrary to apparently widespread misconceptions, sequence-based nomenclature cannot be based on similarity-derived OTUs and their consensus sequences, but must be derived from rigorous, multiple alignment-based phylogenetic methods and quantitative, single-marker species recognition algorithms, using original sequence reads; it is therefore identical in its approach to single-marker studies based on physical types, an approach allowed by the Code. We recognize the limitations of the ITS as a single fungal barcoding marker, but point out that these result in a conservative approach, with "false negatives" surpassing "false positives"; a desirable feature of sequence-based nomenclature. Sequence-based nomenclature does not aim at accurately resolving species, but at naming sequences that represent unknown fungal lineages so that these can serve as a means of communication, so ending the untenable situation of an exponentially growing number of unlabeled fungal sequences that fill online repositories. The risks are outweighed by the gains obtained by a reference library of named sequences spanning the full array of fungal diversity. Finally, we elaborate provisions in addition to our original proposal to amend the Code that would take care of the issues brought forward by opponents to this approach. In particular, taking up the idea of the Candidatus status of invalid, provisional names in prokaryote nomenclature, we propose a compromise that would allow valid publication of voucherless, sequence-based names in a consistent manner, but with the obligate designation as "nom. seq." (nomen sequentiae). Such names would not have priority over specimen- or culture-based names unless either epitypified with a physical type or adopted for protection on the recommendation of a committee of the International Commission on the Taxonomy of Fungi following evaluation based on strict quality control of the underlying studies based on established rules or recommendations.

Ten reasons why a sequence-based nomenclature is not useful for fungi anytime soon

The large number of species still to be discovered in fungi, together with an exponentially growing number of environmental sequences that cannot be linked to known taxa, has fuelled the idea that it might be necessary to formally name fungi on the basis of sequence data only. Here we object to this idea due to several shortcomings of the approach, ranging from concerns regarding reproducibility and the violation of general scientific principles to ethical issues. We come to the conclusion that sequence-based nomenclature is potentially harmful for mycology as a discipline. Additionally, a classification based on sequences as types is not within reach anytime soon, because there is a lack of consensus regarding common standards due to the fast pace at which sequencing technologies develop.

A Nuclear Ribosomal DNA Phylogeny of Acer Inferred with Maximum Likelihood, Splits Graphs, and Motif Analysis of 606 Sequences

The multi-copy internal transcribed spacer (ITS) region of nuclear ribosomal DNA is widely used to infer phylogenetic relationships among closely related taxa. Here we use maximum likelihood (ML) and splits graph analyses to extract phylogenetic information from ~ 600 mostly cloned ITS sequences, representing 81 species and subspecies of Acer, and both species of its sister Dipteronia. Additional analyses compared sequence motifs in Acer and several hundred Ana-cardiaceae, Burseraceae, Meliaceae, Rutaceae, and Sapindaceae ITS sequences in GenBank. We also assessed the effects of using smaller data sets of consensus sequences with ambiguity coding (accounting for within-species variation) instead of the full (partly redundant) original sequences. Neighbor-nets and bipartition networks were used to visualize conflict among character state patterns. Species clusters observed in the trees and networks largely agree with morphology-based classifications; of de Jong's (1994) 16 sections, nine are supported in neighbor-net and bipartition networks, and ten by sequence motifs and the ML tree; of his 19 series, 14 are supported in networks, motifs, and the ML tree. Most nodes had higher bootstrap support with matrices of 105 or 40 consensus sequences than with the original matrix. Within-taxon ITS divergence did not differ between diploid and polyploid Acer, and there was little evidence of differentiated parental ITS haplotypes, suggesting that concerted evolution in Acer acts rapidly.

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.

Astonishing 35S rDNA diversity in the gymnosperm species Cycas revoluta Thunb

In all eukaryotes, the highly repeated 35S ribosomal DNA (rDNA) sequences encoding 18S-5.8S-26S ribosomal RNA (rRNA) typically show high levels of intragenomic uniformity due to homogenisation processes, leading to concerted evolution of 35S rDNA repeats. Here, we compared 35S rDNA divergence in several seed plants using next generation sequencing and a range of molecular and cytogenetic approaches. Most species showed similar 35S rDNA homogeneity indicating concerted evolution. However, Cycas revoluta exhibits an extraordinary diversity of rDNA repeats (nucleotide sequence divergence of different copies averaging 12 %), influencing both the coding and non-coding rDNA regions nearly equally. In contrast, its rRNA transcriptome was highly homogeneous suggesting that only a minority of genes (<20 %) encode functional rRNA. The most common SNPs were C > T substitutions located in symmetrical CG and CHG contexts which were also highly methylated. Both functional genes and pseudogenes appear to cluster on chromosomes. The extraordinary high levels of 35S rDNA diversity in C. revoluta, and probably other species of cycads, indicate that the frequency of repeat homogenisation has been much lower in this lineage, compared with all other land plant lineages studied. This has led to the accumulation of methylation-driven mutations and pseudogenisation. Potentially, the reduced homology between paralogs prevented their elimination by homologous recombination, resulting in long-term retention of rDNA pseudogenes in the genome.

Genomic Resources of Three Pulsatilla Species Reveal Evolutionary Hotspots, Species-Specific Sites and Variable Plastid Structure in the Family Ranunculaceae

Background: The European continent is presently colonized by nine species of the genus Pulsatilla, five of which are encountered only in mountainous regions of southwest and south-central Europe. The remaining four species inhabit lowlands in the north-central and eastern parts of the continent. Most plants of the genus Pulsatilla are rare and endangered, which is why most research efforts focused on their biology, ecology and hybridization. The objective of this study was to develop genomic resources, including complete plastid genomes and nuclear rRNA clusters, for three sympatric Pulsatilla species that are most commonly found in Central Europe. The results will supply valuable information about genetic variation, which can be used in the process of designing primers for population studies and conservation genetics research. The complete plastid genomes together with the nuclear rRNA cluster can serve as a useful tool in hybridization studies. Methodology/principal findings: Six complete plastid genomes and nuclear rRNA clusters were sequenced from three species of Pulsatilla using the Illumina sequencing technology. Four junctions between single copy regions and inverted repeats and junctions between the identified locally-collinear blocks (LCB) were confirmed by Sanger sequencing. Pulsatilla genomes of 120 unique genes had a total length of approximately 161–162 kb, and 21 were duplicated in the inverted repeats (IR) region. Comparative plastid genomes of newly-sequenced Pulsatilla and the previously-identified plastomes of Aconitum and Ranunculus species belonging to the family Ranunculaceae revealed several variations in the structure of the genome, but the gene content remained constant. The nuclear rRNA cluster (18S-ITS1-5.8S-ITS2-26S) of studied Pulsatilla species is 5795 bp long. Among five analyzed regions of the rRNA cluster, only Internal Transcribed Spacer 2 (ITS2) enabled the molecular delimitation of closely-related Pulsatilla patens and Pulsatillavernalis. Conclusions/significance: The determination of complete plastid genome and nuclear rRNA cluster sequences in three species of the genus Pulsatilla is an important contribution to our knowledge of the evolution and phylogeography of those endangered taxa. The resulting data can be used to identify regions that are particularly useful for barcoding, phylogenetic and phylogeographic studies. The investigated taxa can be identified at each stage of development based on their species-specific SNPs. The nuclear and plastid genomic resources enable advanced studies on hybridization, including identification of parent species, including their roles in that process. The identified nonsynonymous mutations could play an important role in adaptations to changing environments. The results of the study will also provide valuable information about the evolution of the plastome structure in the family Ranunculaceae.

5.8S rRNA sequence and secondary structure in Monotropa hypopitys and related Ericaceae species

This is the first study to investigate the secondary structure of 5.8S rRNA in M. hypopitys and related Ericaceae species. The identified nucleotide substitutions are localized mostly at the 3' and 5' ends of the gene, in the region of the third hairpin, and do not significantly affect the topology of the secondary structure of the 5.8S rRNA molecule.

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.

Phylogenetic reconstruction of endophytic fungal isolates using internal transcribed spacer 2 (ITS2) region

Endophytic fungi are inhabitants of plants, living most part of their lifecycle asymptomatically which mainly confer protection and ecological advantages to the host plant. In this present study, 48 endophytic fungi were isolated from the leaves of three medicinal plants and characterized based on ITS2 sequence – secondary structure analysis. ITS2 secondary structures were elucidated with minimum free energy method (MFOLD version 3.1) and consensus structure of each genus was generated by 4SALE. ProfDistS was used to generate ITS2 sequence structure based phylogenetic tree respectively. Our elucidated isolates were belonging to Ascomycetes family, representing 5 orders and 6 genera. Colletotrichum/Glomerella spp., Diaporthae/Phomopsis spp., and Alternaria spp., were predominantly observed while Cochliobolus sp., Cladosporium sp., and Emericella sp., were represented by singletons. The constructed phylogenetic tree has well resolved monophyletic groups with >50% bootstrap value support. Secondary structures based fungal systematics improves not only the stability; it also increases the precision of phylogenetic inference. Above ITS2 based phylogenetic analysis was performed for our 48 isolates along with sequences of known ex-types taken from GenBank which confirms the efficiency of the proposed method. Further, we propose it as superlative marker for reconstructing phylogenetic relationships at different taxonomic levels due to their lesser length.

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.

Evolution of rDNA ITS1 and ITS2 sequences and RNA secondary structures within members of the fungal genera Grosmannia and Leptographium

The two internal transcribed spacers (ITS) of the nuclear ribosomal (r) DNA tandem repeat were examined in ophiostomatoid fungi belonging to the genera Grosmannia and Leptographium and closely-related taxa. Although the DNA sequence of the ITS region evolves rapidly, core features of the RNA secondary structure of the ITS1 and ITS2 segments are conserved. The results demonstrate that structural conservation of GC-rich helical regions is facilitated primarily through compensatory base changes (CBCs), hemi-CBCs, and compensating insertions/deletions (indels), although slippage of the RNA strand is potentially an additional mechanism for maintaining basepairing interactions. The major conclusion of the structural analysis of both ITS segments is that two factors appear to be involved in limiting the type of changes observed: a high GC bias for both ITS1 and ITS2 and structural constraints at the RNA level.

Reconsidering the generation time hypothesis based on nuclear ribosomal ITS sequence comparisons in annual and perennial angiosperms

Background

Differences in plant annual/perennial habit are hypothesized to cause a generation time effect on divergence rates. Previous studies that compared rates of divergence for internal transcribed spacer (ITS1 and ITS2) sequences of nuclear ribosomal DNA (nrDNA) in angiosperms have reached contradictory conclusions about whether differences in generation times (or other life history features) are associated with divergence rate heterogeneity. We compared annual/perennial ITS divergence rates using published sequence data, employing sampling criteria to control for possible artifacts that might obscure any actual rate variation caused by annual/perennial differences.

Results

Relative rate tests employing ITS sequences from 16 phylogenetically-independent annual/perennial species pairs rejected rate homogeneity in only a few comparisons, with annuals more frequently exhibiting faster substitution rates. Treating branch length differences categorically (annual faster or perennial faster regardless of magnitude) with a sign test often indicated an excess of annuals with faster substitution rates. Annuals showed an approximately 1.6-fold rate acceleration in nucleotide substitution models for ITS. Relative rates of three nuclear loci and two chloroplast regions for the annual Arabidopsis thaliana compared with two closely related Arabidopsis perennials indicated that divergence was faster for the annual. In contrast, A. thaliana ITS divergence rates were sometimes faster and sometimes slower than the perennial. In simulations, divergence rate differences of at least 3.5-fold were required to reject rate constancy in > 80 % of replicates using a nucleotide substitution model observed for the combination of ITS1 and ITS2. Simulations also showed that categorical treatment of branch length differences detected rate heterogeneity > 80% of the time with a 1.5-fold or greater rate difference.

Conclusion

Although rate homogeneity was not rejected in many comparisons, in cases of significant rate heterogeneity annuals frequently exhibited faster substitution rates. Our results suggest that annual taxa may exhibit a less than 2-fold rate acceleration at ITS. Since the rate difference is small and ITS lacks statistical power to reject rate homogeneity, further studies with greater power will be required to adequately test the hypothesis that annual and perennial plants have heterogeneous substitution rates. Arabidopsis sequence data suggest that relative rate tests based on multiple loci may be able to distinguish a weak acceleration in annual plants. The failure to detect rate heterogeneity with ITS in past studies may be largely a product of low statistical power.

Molecular evolution and phylogenetic utility of the internal transcribed spacer 2 (ITS2) in Calyptratae (Diptera: Brachycera)

The resolution potential of internal transcribed spacer 2 (ITS2) at deeper levels remains controversial. In this study, 105 ITS2 sequences of 55 species in Calyptratae were analyzed to examine the phylogenetic utility of the spacer above the subfamily level and to further understand its evolutionary characteristics. We predicted the secondary structure of each sequence using the minimum-energy algorithm and constructed two data matrixes for phylogenetic analysis. The ITS2 regions of Calyptratae display strong A-T bias and slight variation in length. The tandem and dispersed repeats embedded in the spacers possibly resulted from replication slippage or transposition. Most foldings conformed to the four-domain model. Sequence comparison in combination with the secondary structures revealed six conserved motifs. Covariation analysis from the conserved motifs indicated that the secondary structure restrains the sequence evolution of the spacer. The deep-level phylogeny derived from the ITS2 data largely agreed with the phylogenetic hypotheses from morphologic and other molecular evidence. Our analyses suggest that the accordant resolutions generated from different analyses can be used to infer deep-level phylogenetic relations.

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.

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.

General functions to transform associate data to host data, and their use in phylogenetic inference from sequences with intra-individual variability

Background

Amongst the most commonly used molecular markers for plant phylogenetic studies are the nuclear ribosomal internal transcribed spacers (ITS). Intra-individual variability of these multicopy regions is a very common phenomenon in plants, the causes of which are debated in literature. Phylogenetic reconstruction under these conditions is inherently difficult. Our approach is to consider this problem as a special case of the general biological question of how to infer the characteristics of hosts (represented here by plant individuals) from features of their associates (represented by cloned sequences here).

Results

Six general transformation functions are introduced, covering the transformation of associate characters to discrete and continuous host characters, and the transformation of associate distances to host distances. A pure distance-based framework is established in which these transformation functions are applied to ITS sequences collected from the angiosperm genera Acer, Fagus and Zelkova. The formulae are also applied to allelic data of three different loci obtained from Rosa spp. The functions are validated by (1) phylogeny-independent measures of treelikeness; (2) correlation with independent host characters; (3) visualization using splits graphs and comparison with published data on the test organisms. The results agree well with these three measures and the datasets examined as well as with the theoretical predictions and previous results in the literature. High-quality distance matrices are obtained with four of the six transformation formulae. We demonstrate that one of them represents a generalization of the Sørensen coefficient, which is widely applied in ecology.

Conclusion

Because of their generality, the transformation functions may be applied to a wide range of biological problems that are interpretable in terms of hosts and associates. Regarding cloned sequences, the formulae have a high potential to accurately reflect evolutionary relationships within angiosperm genera, and to identify hybrids and ancestral taxa. These results corroborate earlier ones which showed that treelikeness measures are a valuable tool in comparative studies of biological distance functions.

Domain II hairpin structure in ITS1 sequences as an aid in differentiating recently evolved animal and plant pathogenic fungi

The hypothesis that ITS structural features can be used to define fungal groups, where sequence analysis is unsatisfactory, was examined in plant and animal pathogenic fungi. Structural models of ITS1 regions were predicted for presumed closely related species in Colletotrichum and Trichophyton anamorphs of Arthroderma species. Structural alignment of models and comparison with ITS sequence analysis identified a variable region in a conserved hairpin formed from a common inverted repeat. Thirteen different hairpin structure models were obtained for Colletotrichum species and five different models were obtained for Trichophyton species. The different structure types could be matched to individual species and species complexes as defined by ITS sequence analysis.

Its evolution in Platanus (Platanaceae): homoeologues, pseudogenes and ancient hybridization

BACKGROUND AND AIMS

Platanaceae is an old family of angiosperms extending back to the Early Cretaceous but consisting of a single extant genus, Platanus. Species of Platanus have long been known to hybridize, and the London plane, Platanus x hispanica, is a well-known example of a hybrid species that formed in historical times. In addition, morphological studies have suggested past interspecific or interlineage hybridization and reticulation as possibly important factors in the evolution of the genus. This study aims at unravelling the complex evolutionary information contained in internal transcribed spacer (ITS) sequences.

METHODS

The ITS of the 35S nuclear ribosomal DNA are biparentally inherited, multi-copy markers with a high potential for resolving intrageneric relationships especially when ancient hybridization (reticulation) is involved. Phylogenetic trees, splits graphs and motif analysis are used to extract phylogenetic information from 223 cloned ITS sequences, representing ten species and varieties of Platanus. Non-pseudogenous and pseudogenous sequence motives are assessed to explain how different evolutionary modes contribute to possibly conflicting character state patterns in the ITS.

KEY RESULTS

It was found that putative non-functional ITS copies ('pseudogenes') form distinct groups in phylograms and splits graphs, and that pseudogenous lineages reflect ancient hybridization events conserved in the ITS. Specifically, pseudogenous clones of an 'Atlantic' North American clade share sequence motives with non-pseudogenous clones of the western ('Pacific') North American P. racemosa species aggregate. In addition, evidence was found for recent lateral gene flow as a possible factor in the evolution of the central Mexican P. rzedowskii.

CONCLUSIONS

Broad ITS data sets that cover intra- and interindividual variability reveal past and ongoing speciation processes in Platanus. Evolutionary pathways can be visualized with splits graphs, but not with bifurcating trees.

Species identification of some common necrophagous flies in Guangdong province, southern China based on the rDNA internal transcribed spacer 2 (ITS2)

Recent studies suggest that sequence analysis technique displays a tempting foreground in identifying unknown specimens of necrophagous flies. In this study, we analyzed 63 complete ITS2 sequences concerning 29 fly species to evaluate the identification potential of the ITS2 region, among of which 41 sequence entries were obtained by sequencing and 22 sequence entries were available on the line. Additionally, phenetic method was recommended to substitute for phylogenetic method because it is very difficult to align the ITS2 sequences. The neighbor-joining tree generated by clustalx1.81 allowed us to differentiate each species. Meanwhile the tree topology also suggested that the ITS2 region showed no resolution for the distinction of geographical populations of some species. The overlapping between intra- and interspecific variation revealed by sequence analyses did not affect species identification. High sequence homology between some congeneric species required further sequencing for forensic practice.

Dating dispersal and radiation in the gymnosperm Gnetum (Gnetales)--clock calibration when outgroup relationships are uncertain

Most implementations of molecular clocks require resolved topologies. However, one of the Bayesian relaxed clock approaches accepts input topologies that include polytomies. We explored the effects of resolved and polytomous input topologies in a rate-heterogeneous sequence data set for Gnetum, a member of the seed plant lineage Gnetales. Gnetum has 10 species in South America, 1 in tropical West Africa, and 20 to 25 in tropical Asia, and explanations for the ages of these disjunctions involve long-distance dispersal and/or the breakup of Gondwana. To resolve relationships within Gnetum, we sequenced most of its species for six loci from the chloroplast (rbcL, matK, and the trnT-trnF region), the nucleus (rITS/5.8S and the LEAFY gene second intron), and the mitochondrion (nad1 gene second intron). Because Gnetum has no fossil record, we relied on fossils from other Gnetales and from the seed plant lineages conifers, Ginkgo, cycads, and angiosperms to constrain a molecular clock and obtain absolute times for within-Gnetum divergence events. Relationships among Gnetales and the other seed plant lineages are still unresolved, and we therefore used differently resolved topologies, including one that contained a basal polytomy among gymnosperms. For a small set of Gnetales exemplars (n = 13) in which rbcL and matK satisfied the clock assumption, we also obtained time estimates from a strict clock, calibrated with one outgroup fossil. The changing hierarchical relationships among seed plants (and accordingly changing placements of distant fossils) resulted in small changes of within-Gnetum estimates because topologically closest constraints overrode more distant constraints. Regardless of the seed plant topology assumed, relaxed clock estimates suggest that the extant clades of Gnetum began diverging from each other during the Upper Oligocene. Strict clock estimates imply a mid-Miocene divergence. These estimates, together with the phylogeny for Gnetum from the six combined data sets, imply that the single African species of Gnetum is not a remnant of a once Gondwanan distribution. Miocene and Pliocene range expansions are inferred for the Asian subclades of Gnetum, which stem from an ancestor that arrived from Africa. These findings fit with seed dispersal by water in several species of Gnetum, morphological similarities among apparently young species, and incomplete concerted evolution in the nuclear ITS region.

Wax plants disentangled: a phylogeny of Hoya (Marsdenieae, Apocynaceae) inferred from nuclear and chloroplast DNA sequences

Hoya (Marsdenieae, Apocynaceae) includes at least 200 species distributed from India to the Pacific Islands. We here infer major species groups in the genus based on combined sequences from the chloroplast atpB-rbcL spacer, the trnL region, and nuclear ribosomal DNA ITS region for 42 taxa of Hoya and close relatives. To assess levels of ITS polymorphism, ITS sequences for a third of the accessions were obtained by cloning. Most ITS clones grouped by species, indicating that speciation in Hoya usually predates ITS duplication. One ITS sequence of H. carnosa, however, grouped with a sequence of the morphologically similar H. pubicalyx, pointing to recent hybridization or the persistence of paralogous copies through a speciation event. The topology resulting from the combined chloroplast and nuclear data recovers some morphology-based sections, such as Acanthostemma and Eriostemma, as well as a well-supported Australian/New Guinean clade. The combined data also suggest that morphological adaptations for ant-symbiosis evolved at least three times within Hoya.

The chloroplast trnT-trnF region in the seed plant lineage Gnetales

The trnT-trnF region is located in the large single-copy region of the chloroplast genome. It consists of the trnL intron, a group I intron, and the trnT-trnL and trnL-trnF intergenic spacers. We analyzed the evolution of the region in the three genera of the gymnosperm lineage Gnetales (Gnetum, Welwitschia, and Ephedra), with especially dense sampling in Gnetum for which we sequenced 41 accessions, representing most of the 25-35 species. The trnL intron has a conserved secondary structure and contains elements that are homologous across land plants, while the spacers are so variable in length and composition that homology cannot be found even among the three genera. Palindromic sequences that form hairpin structures were detected in the trnL-trnF spacer, but neither spacer contained promoter elements for the tRNA genes. The absence of promoters, presence of hairpin structures in the trnL-trnF spacer, and high sequence variation in both spacers together suggest that trnT and trnF are independently transcribed. Our model for the expression and processing of the genes tRNA(Thr)(UGU), tRNA(Leu)(UAA), and tRNA(Phe) (GAA) therefore attributes the seemingly neutral evolution of the two spacers to their escape from functional constraints.