Employing 454 amplicon pyrosequencing to reveal intragenomic divergence in the internal transcribed spacer rDNA region in fungi

The rDNA internal transcribed spacer (ITS) region has been accepted as a DNA barcoding marker for fungi and is widely used in phylogenetic studies; however, intragenomic ITS variability has been observed in a broad range of taxa, including prokaryotes, plants, animals, and fungi, and this variability has the potential to inflate species richness estimates in molecular investigations of environmental samples. In this study 454 amplicon pyrosequencing of the ITS1 region was applied to 99 phylogenetically diverse axenic single‐spore cultures of fungi (Dikarya: Ascomycota and Basidiomycota) to investigate levels of intragenomic variation. Three species (one Basidiomycota and two Ascomycota), in addition to a positive control species known to contain ITS paralogs, displayed levels of molecular variation indicative of intragenomic variation; taxon inflation due to presumed intragenomic variation was ≈9%. Intragenomic variability in the ITS region appears to be widespread but relatively rare in fungi (≈3–5% of species investigated in this study), suggesting this problem may have minor impacts on species richness estimates relative to PCR and/or pyrosequencing errors. Our results indicate that 454 amplicon pyrosequencing represents a powerful tool for investigating levels of ITS intragenomic variability across taxa, which may be valuable for better understanding the fundamental mechanisms underlying concerted evolution of repetitive DNA regions.     

Molecular characteristics of Camallanus Spp. (Spirurida: Camallanidae) in fishes from China based on its rDNA sequences

In the paper, we explored the intra- and interspecific evolutionary variation among species of Camallanus collected from different fish species in various regions of China. We determined the internal transcribed spacers of ribosomal DNA (ITS rDNA) sequences of these nematodes. The divergence (uncorrected p-distance) of ITS1, ITS2, and ITS rDNA data sets confirmed 2 valid species of Camallanus in China, i.e., C. cotti and C. hypophthalmichthys. The 2 species were distinguished not only by their different morphologies and host ranges but also by a tetranucleotide microsatellite (TTGC)n present in the ITS1 region of C. cotti. Phylogenetic analyses of the nematodes disclosed 2 main clades, corresponding to different individuals of C. cotti and C. hypophthalmichthys from different fish species in various geographical locations, although the interior nodes of each clade received poor support.     

Phylogenetic relationships among the Caribbean members of the Cliona viridis complex (Porifera, Demospongiae, Hadromerida) using nuclear and mitochondrial DNA sequences

Species complexes - groups of closely related species in which intraspecific and interspecific variability overlap - have generated considerable interest and study. Frequently, members of a species complex do not have complete reproductive isolation; therefore, the complex may go through extensive gene flow. In the Caribbean Sea, some encrusting and excavating sponges of the genus Cliona (Porifera, Hadromerida, Clionaidae) are grouped within the great "Cliona viridis" complex because of their morphological similarities. This study examined the evolutionary relationships of the Caribbean members of this complex (C. caribbaea, C. tenuis, C. aprica and C. varians) and related taxa based on nuclear (ITS1 and ITS2) and mitochondrial (3' end of ND6) DNA sequences. The intragenomic ITS variation and its secondary structures were evaluated using a mixed approach of Denaturing Gradient Gel Electrophoresis (DGGE), DNA sequencing and secondary structure prediction. Considerable intragenomic variation was found in all the species, with apparently functional ITS1 and ITS2 secondary structures. Despite the subtle but clear morphological differentiation in these excavating sponges, the intragenomic copies of C. caribbaea, C. tenuis and C. aprica had a polyphyletic placement in the ITS1 and ITS2 genealogies and very low divergence. Therefore, it is clear that these species constitute a species complex (herein called Ct-complex). Genetic distances within the Ct-complex revealed that an important part of the interspecific variation overlapped with intraspecific variation, suggesting either incomplete lineage sorting or extensive gene flow. In contrast, C. varians and an unidentified "Pione" species emerged as monophyletic clades, being the closest sister groups to the Ct-complex. Additionally, our results support that C. laticavicola and C. delitrix conform a monophyletic group, but absence of reciprocal monophyly in these species suggests they may be life stages or ecophenotypes of a single species or they have diverged recently. Our work showed that the 3' end of the ND6 mitochondrial gene was highly conserved and not suitable for phylogenetic analysis at the interspecific level.     

Compensatory Base Changes in ITS2 Secondary Structures Correlate with the Biological Species Concept Despite Intragenomic Variability in ITS2 Sequences – A Proof of Concept

Compensatory base changes (CBCs) in internal transcribed spacer 2 (ITS2) rDNA secondary structures correlate with Ernst Mayr’s biological species concept. This hypothesis also referred to as the CBC species concept recently was subjected to large-scale testing, indicating two distinct probabilities. (1) If there is a CBC then there are two different species with a probability of ∼0.93. (2) If there is no CBC then there is the same species with a probability of ∼0.76. In ITS2 research, however, the main problem is the multicopy nature of ITS2 sequences. Most recently, 454 pyrosequencing data have been used to characterize more than 5000 intragenomic variations of ITS2 regions from 178 plant species, demonstrating that mutation of ITS2 is frequent, with a mean of 35 variants per species, respectively per individual organism. In this study, using those 454 data, the CBC criterion is reconsidered in the light of intragenomic variability, a proof of concept, a necessary criterion, expecting no intragenomic CBCs in variant ITS2 copies. In accordance with the CBC species concept, we could demonstrate that the probability that there is no intragenomic CBC is ∼0.99.     

Low variability of internal transcribed spacer rDNA and trnL (UAA) intron sequences of several taxa in the Festuca ovina aggregate (POACEAE)

Identification and classification of numerous Festuca species is still a difficult problem due to the close morphological resemblance. The most difficult fine fescues to identify belong to the Festuca ovina aggregate, which is the largest group in the genus Festuca. Many taxons are considered to be separate species based on quantitative taxonomic characters, differences in ploidy level or the structure of sclerenchyma cells. In order to evaluate the taxonomic value of DNA-based markers, sequence analysis of the internal transcribed spacer (ITS1-5.8S-ITS2) region and the chloroplast trnL (UAA) intron was performed in the ten most problematic fine fescues belonging to the Festuca ovina aggregate. Intraspecific ITS variants were found in a single case while in other cases only intragenomic ITS polymorphisms were detected with 1-2 ambiguous positions. Among the sequences of the trnL (UAA) intron even intragenomic polymorphisms were not detected in any of the Festuca species studied. Thus, the results do not support the species status of these ten taxa.     

Assessing Symbiodinium diversity in scleractinian corals via next‐generation sequencing‐based genotyping of the ITS2 rDNA region

The persistence of coral reef ecosystems relies on the symbiotic relationship between scleractinian corals and intracellular, photosynthetic dinoflagellates in the genus Symbiodinium. Genetic evidence indicates that these symbionts are biologically diverse and exhibit discrete patterns of environmental and host distribution. This makes the assessment of Symbiodinium diversity critical to understanding the symbiosis ecology of corals. Here, we applied pyrosequencing to the elucidation of Symbiodinium diversity via analysis of the internal transcribed spacer 2 (ITS2) region, a multicopy genetic marker commonly used to analyse Symbiodinium diversity. Replicated data generated from isoclonal Symbiodinium cultures showed that all genomes contained numerous, yet mostly rare, ITS2 sequence variants. Pyrosequencing data were consistent with more traditional denaturing gradient gel electrophoresis (DGGE) approaches to the screening of ITS2 PCR amplifications, where the most common sequences appeared as the most intense bands. Further, we developed an operational taxonomic unit (OTU)‐based pipeline for Symbiodinium ITS2 diversity typing to provisionally resolve ecologically discrete entities from intragenomic variation. A genetic distance cut‐off of 0.03 collapsed intragenomic ITS2 variants of isoclonal cultures into single OTUs. When applied to the analysis of field‐collected coral samples, our analyses confirm that much of the commonly observed Symbiodinium ITS2 diversity can be attributed to intragenomic variation. We conclude that by analysing Symbiodinium populations in an OTU‐based framework, we can improve objectivity, comparability and simplicity when assessing ITS2 diversity in field‐based studies.     

Genomic structure of the autotetraploid oat species Avena macrostachya inferred from comparative analysis of the ITS1 and ITS2 sequences: on the oat karyotype evolution during the early stages of the Avena species divergence

To examine the genomic structure of Avena macrostachya, internal transcribed spacers, ITS1 and ITS2, as well as nuclear 5.8S tRNA genes from three oat species with AsAs karyotype (A. wiestii, A. hirtula, and A. atlantica), and those from A. longiglumis (AlAl), A. canariensis (AcAc), A. ventricosa (CvCv), A. pilosa, and A. clauda (CpCp) were sequenced. All species of the genus Avena examined represented a monophyletic group (bootstrap index = 98), within which two branches, i.e., species with A- and C-genomes, were distinguished (bootstrap indices = 100). The subject of our study, A. macrostachya, albeit belonging to the phylogenetic branch of C-genome oat species (karyotype with submetacentic and subacrocentric chromosomes), has preserved an isobrachyal karyotype, (i.e., that containing metacentric chromosomes), probably typical of the common Avena ancestor. It was suggested to classify the A. macrostachya genome as a specific form of C-genome, Cm-genome. Among the species from other genera studied, Arrhenatherum elatius was found to be the closest to Avena in ITS1 and ITS structure. Phylogenetic relationships between Avena and Helictotrichon remain intriguingly uncertain. The HPR389153 sequence from H. pratense genome was closest to the ITS1 sequences specific to the Avena A-genomes (p-distance = 0.0237), while the differences of this sequence from the ITS1 of A. macrostachya reached 0.1221. On the other hand, HAD389117 from H. adsurgens was close to the ITS1 specific to Avena C-genomes (p-distance = 0.0189), while its differences from the A-genome specific ITS1 sequences reached 0.1221. It seems likely that the appearance of highly polyploid (2n = 12-21x) species of H. pratense and H. adsurgens could be associated with interspecific hybridization involving Mediterranean oat species carrying A- and C-genomes. A hypothesis on the pathways of Avena chromosomes evolution during the early stages the oat species divergence is proposed.     

Improved resolution of reef-coral endosymbiont (Symbiodinium) species diversity, ecology, and evolution through psbA non-coding region genotyping

Ribosomal DNA sequence data abounds from numerous studies on the dinoflagellate endosymbionts of corals, and yet the multi-copy nature and intragenomic variability of rRNA genes and spacers confound interpretations of symbiont diversity and ecology. Making consistent sense of extensive sequence variation in a meaningful ecological and evolutionary context would benefit from the application of additional genetic markers. Sequences of the non-coding region of the plastid psbA minicircle (psbA(ncr)) were used to independently examine symbiont genotypic and species diversity found within and between colonies of Hawaiian reef corals in the genus Montipora. A single psbA(ncr) haplotype was recovered in most samples through direct sequencing (~80-90%) and members of the same internal transcribed spacer region 2 (ITS2) type were phylogenetically differentiated from other ITS2 types by substantial psbA(ncr) sequence divergence. The repeated sequencing of bacterially-cloned fragments of psbA(ncr) from samples and clonal cultures often recovered a single numerically common haplotype accompanied by rare, highly-similar, sequence variants. When sequence artifacts of cloning and intragenomic variation are factored out, these data indicate that most colonies harbored one dominant Symbiodinium genotype. The cloning and sequencing of ITS2 DNA amplified from these same samples recovered numerically abundant variants (that are diagnostic of distinct Symbiodinium lineages), but also generated a large amount of sequences comprising PCR/cloning artifacts combined with ancestral and/or rare variants that, if incorporated into phylogenetic reconstructions, confound how small sequence differences are interpreted. Finally, psbA(ncr) sequence data from a broad sampling of Symbiodinium diversity obtained from various corals throughout the Indo-Pacific were concordant with ITS lineage membership (defined by denaturing gradient gel electrophoresis screening), yet exhibited substantially greater sequence divergence and revealed strong phylogeographic structure corresponding to major biogeographic provinces. The detailed genetic resolution provided by psbA(ncr) data brings further clarity to the ecology, evolution, and systematics of symbiotic dinoflagellates.     

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.     

Molecular identification of species in Juglandaceae:A tiered method

DNA barcoding is a method of species identification and recognition using DNA sequence data. A tiered or multilocus method has been recommended for barcoding plant species. In this study, we sampled 196 individuals representing 9 genera and 54 species of Juglandaceae to investigate the utility of the four potential barcoding loci (rbcL, matK, trnH-psbA, and internal transcribed spacer (ITS)). Our results show that all four DNA regions are easy to amplify and sequence. In the four tested DNA regions, ITS has the most variable information, and rbcL has the least. At generic level, seven of nine genera can be efficiently identified by matK. At species level, ITS has higher interspecific p-distance than the trnH-psbA region. Difficult to align in the whole family, ITS showed heterogeneous variability among different genera. Except for the monotypic genera (Cyclocarya, Annamocarya, Platycarya), ITS appeared to have limited power for species identification within the Carya and Engelhardia complex, and have no power for Juglans or Pterocarya. Overall, our results confirmed that a multilocus tiered method for plant barcoding was applicable and practicable. With higher priority, matK is proposed as the first-tier DNA region for genus discrimination, and the second locus at species level should have enough stable variable characters.     

Molecular identification of species in Juglandaceae: A tiered method

DNA barcoding is a method of species identification and recognition using DNA sequence data. A tiered or multilocus method has been recommended for barcoding plant species. In this study, we sampled 196 individuals representing 9 genera and 54 species of Juglandaceae to investigate the utility of the four potential barcoding loci (rbcL, matK, trnH-psbA, and internal transcribed spacer (ITS)). Our results show that all four DNA regions are easy to amplify and sequence. In the four tested DNA regions, ITS has the most variable information, and rbcL has the least. At generic level, seven of nine genera can be efficiently identified by matK. At species level, ITS has higher interspecific p-distance than the trnH-psbA region. Difficult to align in the whole family, ITS showed heterogeneous variability among different genera. Except for the monotypic genera (Cyclocarya, Annamocarya, Platycarya), ITS appeared to have limited power for species identification within the Carya and Engelhardia complex, and have no power for Juglans or Pterocarya. Overall, our results confirmed that a multilocus tiered method for plant barcoding was applicable and practicable. With higher priority, matK is proposed as the first-tier DNA region for genus discrimination, and the second locus at species level should have enough stable variable characters.     

SymPortal: A novel analytical framework and platform for coral algal symbiont next‐generation sequencing ITS2 profiling

We present SymPortal (SymPortal.org), a novel analytical framework and platform for genetically resolving the algal symbionts of reef corals using next‐generation sequencing (NGS) data of the ITS2 rDNA. Although the ITS2 marker is widely used to genetically characterize taxa within the family Symbiodiniaceae (formerly the genus Symbiodinium), the multicopy nature of the marker complicates its use. Commonly, the intragenomic diversity resultant from this multicopy nature is collapsed by analytical approaches, thereby focusing on only the most abundant sequences. In contrast, SymPortal employs logic to identify within‐sample informative intragenomic sequences, which we have termed ‘defining intragenomic variants' (DIVs), to identify ITS2‐type profiles representative of putative Symbiodiniaceae taxa. By making use of this intragenomic ITS2 diversity, SymPortal is able to resolve genetic delineations using the ITS2 marker at a level that was previously only possible by using additional genetic markers. We demonstrate this by comparing this novel approach to the most commonly used alternative approach for NGS ITS2 data, the 97% similarity clustering to operational taxonomic units (OTUs). The SymPortal platform accepts NGS raw sequencing data as input to provide an easy‐to‐use, standardization‐enforced, and community‐driven framework that integrates with a database to gain resolving power with increased use. We consider that SymPortal, in conjunction with ongoing large‐scale sampling and sequencing efforts, should play an instrumental role in making future sampling efforts more comparable and in maximizing their efficacy in working towards the classification of the global Symbiodiniaceae diversity.     

INTRAGENOMIC NRDNA ITS POLYMORPHISM IN THE STEPHANODISCUS NIAGARAE COMPLEX (BACILLARIOPHYCEAE)

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.     

Intragenomic variation in nuclear ribosomal markers and its implication in species delimitation,identification and barcoding in fungi

Intragenomic variation is the molecular variation within the genome among repetitive DNA. As a multigene family, nuclear ribosomal DNA (rDNA) has been widely used in fungal taxonomy for their ease in amplification and suitable variability to attain various levels of taxonomic resolution. At the intraspecific level, rDNA is believed to be under concerted evolution and the internal transcribed spacers (ITS) region is actually accepted as a universal barcoding marker for fungi. However, documentation of intragenomic variation of rDNA indicated that it can be problematic in species delimitation and identification. Fungal taxonomic studies have not generally taken into account the intragenomic variation of rDNA in a systematic manner. In this review, our objective is to address the definition, the origin and the mechanisms for maintenance of intragenomic variation, as well as its implication in the domain of fungal molecular taxonomy, particularly for species delimitation, identification and DNA barcoding. With advanced sequencing technologies (second and third generations), we also addressed how these technologies can be used to study the intragenomic variation of rDNA and also how the intragenomic variation will impact on DNA barcoding via high-throughput sequencing.     

A molecular phylogeny of the genus Ichthyocotylurus (Digenea, Strigeidae)

Three nucleotide data sets, two nuclear (ribosomal internal transcribed spacer regions 1 and 2, ITS1 and ITS2) and one mitochondrial (cytochrome c oxidase subunit 1, CO1), were analysed using distance matrix and maximum likelihood methods to determine the inter-relationships amongst the four species attributed to the genus Ichthyocotylurus Odening, 1969. Sequence data obtained from all gene loci investigated supported the position of Ichthyocotylurus variegatus as a species discrete from Ichthyocotylurus platycephalus. Phylogenetic analyses yielded congruent trees, with I. variegatus isolates comprising a common clade to which I. platycephalus constitutes a sister taxon. Ichthyocotylurus erraticus and Ichthyocotylurus pileatus were found to demonstrate a similarly close inter-specific relationship. The greatest intra-generic divergence occurred in the CO1 region (16% variability), with resultant disparities in three to eight encoded amino acids. PCR amplification yielded multiple ITS1 products for all Ichthyocotylurus spp. Analyses of equivalent-sized amplicons showed 5.4% intra-generic variation and several point mutations between I. variegatus isolates from different geographical localities and from different piscine hosts. The ITS2 locus was extremely conserved, with less than 1% variation between species. No intra-specific variation was recorded for any CO1 or ITS2 sequences.     

樟属植物ITS序列多态性分析

对樟科樟属(Cinnamomum Schaeffer) 17个代表样本的核糖体DNA内转录间隔区(nrDNA ITS)进行克隆测序。对获得的87条不同ITS序列的长度变异、GC含量、5.8S区二级结构的稳定性、遗传距离、进化模式以及系统发育关系进行了相关分析。研究结果显示, ITS序列在樟属植物内存在明显的多态性, 87条序列中的22条序列被鉴定为假基因序列, 其余65条序列为功能基因序列; 假基因序列采用中性进化模式, 变异明显大于功能序列。ITS序列在樟属植物中出现一致性进化不完全和假基因现象也可能发生在樟科其它类群中, 这可能是导致樟科植物ITS序列直接测序方式成功率低的重要原因。     

Intragenomic variation in the ITS rDNA region obscures phylogenetic relationships and inflates estimates of operational taxonomic units in genus Laetiporus

Regions of rDNA are commonly used to infer phylogenetic relationships among fungal species and as DNA barcodes for identification. These regions occur in large tandem arrays, and concerted evolution is believed to reduce intragenomic variation among copies within these arrays, although some variation still might exist. Phylogenetic studies typically use consensus sequencing, which effectively conceals most intragenomic variation, but cloned sequences containing intragenomic variation are becoming prevalent in DNA databases. To understand effects of using cloned rDNA sequences in phylogenetic analyses we amplified and cloned the ITS region from pure cultures of six Laetiporus species and one Wolfiporia species (Basidiomycota, Polyporales). An average of 66 clones were selected randomly and sequenced from 21 cultures, producing a total of 1399 interpretable sequences. Significant variation (≥ 5% variation in sequence similarity) was observed among ITS copies within six cultures from three species clades (L. cincinnatus, L. sp. clade J, and Wolfiporia dilatohypha) and phylogenetic analyses with the cloned sequences produced different trees relative to analyses with consensus sequences. Cloned sequences from L. cincinnatus fell into more than one species clade and numerous cloned L. cincinnatus sequences fell into entirely new clades, which if analyzed on their own most likely would be recognized as "undescribed" or "novel" taxa. The use of a 95% cut off for defining operational taxonomic units (OTUs) produced seven Laetiporus OTUs with consensus ITS sequences and 20 OTUs with cloned ITS sequences. The use of cloned rDNA sequences might be problematic in fungal phylogenetic analyses, as well as in fungal bar-coding initiatives and efforts to detect fungal pathogens in environmental samples.     

南方灵芝ITS、EF1-α和RPB2序列多态性研究

内转录间隔区（ITS）、延伸因子1-α（EF1-α）与RNA聚合酶II第二大亚基（RPB2）是真菌系统学研究中常用的基因片段。已有研究发现,在同一个体内ITS也可能存在差异,但EF1-α和RPB2是否也存在同样的现象却鲜有报道。本研究通过克隆测序,分析了南方灵芝Ganoderma australe的ITS、EF1-α和RPB2在同一个体中的序列差异,并探讨这种序列差异是否会影响分子鉴定的准确性。结果表明,在南方灵芝供试样品中,ITS、EF1-α和RPB2序列都有可能存在个体内变异,但个体内变异程度不尽相同。基于供试样品和克隆测序结果发现,ITS序列在同一个体内的变异最高可达2.3%,而RPB2序列在同一个体内的差异仅0.5%,EF1-α序列在同一个体内序列差异也可高达1.8%。ITS和EF1-α序列在同一个体内的差异可能超过了灵芝属部分物种间的差异,对于部分灵芝属物种来说直接利用克隆序列进行分子鉴定或系统发育分析可能存在一些问题。     

The rDNA ITS Region in the Lessepsian Marine Angiosperm <Emphasis Type="BoldItalic">Halophila stipulacea</Emphasis> (Forssk.) Aschers. (Hydrocharitaceae): Intragenomic Variability and Putative Pseudogenic Sequences

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