Inter- and intrasporal nuclear ribosomal gene sequence variation within one isolate of arbuscular mycorrhizal fungus, Diversispora sp.

Most molecular ecological studies of arbuscular mycorrhizal fungi (AMF) have been based on the rRNA gene sequences. However, information about intraspecific nucleotide variation is still limited in these fungi. In this study, we calculated the inter- and intrasporal nucleotide variation of Diversispora sp. EE1 using 78 cloned sequences from four spores within a ca 4960 bp fragment of the nuclear ribosomal operon spanning the near full length small ribosomal subunit (SSU) rRNA gene, the full internal transcribed spacer (ITS: ITS1-5.8S-ITS2) and ca 2740 bp of the large ribosomal subunit (LSU) rRNA gene. Data for each marker region (SSU, ITS and LSU) originated from the very same spores. Sequence variation resulting from point mutations and small indels was recorded in all regions. Highest sequence variation was observed in the ITS region at both the inter- and intrasporal levels. The ITS1 component was more variable than ITS2, whilst the 5.8S gene was the least variable component of the ITS region. Evolutionary divergence of gene copies between spores was intermediate for the LSU and lowest for the SSU. The SSU and the LSU genes had relatively similar evolutionary divergence per spore. Sequence variant richness was not exhaustive for any of the marker regions, indicating that multiple sequences per spore from multiple spores are needed when characterizing a species. This study provides reference sequences for ecological studies, permitting identification of AMF using any of the ribosomal regions or primer systems.     

帘蛤科贝类rDNA内转录间隔区序列的研究

根据18SrDNA、5．8SrDNA和28SrDNA保守序列设计引物,应用聚合酶链式反应（PCR）扩增了文蛤（Meretrix meretrix L．）、青蛤（Cyclina sinensis G）、硬壳蛤（Mercenaria mercenaria L．）和江户布目蛤（Protothaca jedoensis L．）4种帘蛤科贝类的第一内转录间隔区（ITS1）和第二内转录间隔区（ITS2）序列,并进行了测序。结果表明,文蛤、青蛤、硬壳蛤和江户布目蛤的ITS1扩增产物大小分别为978bp、663bp、757bp和942bp,GC含量分别为61．55％、60．78％、62．48％和64．86％～64．97％,其中ITS1序列长度分别为900bp、585bp、679bp和864bp,是迄今已报道双壳贝类中变化范围最大的,GC含量分别为61．67％、61．03％、63．03％和65．51％～65．62％,江户布目蛤种内ITS1序列有个体差异;ITS2扩增产物大小分别为644bp、618～620bp、593bp和513～514bp,GC含量分别为61．18％、61．29％～61．81％、62．73％和61．48％61．60％,其中ITS2序列长度分别为412bp、386～388bp、361bp和281～282bp,GC含量分别为65．29％、65．21％～66．06％、67．87％和67．38％～67．62％,青蛤和江户布目蛤种内ITS2序列有个体差异。4种蛤ITS1和ITS2序列种间差异很大,有明显的长度多态性,ITS2种间序列相似度73．0％～89．1％,与ITS1的种间序列相似度48．7％～81．5％相比略高。此外,在4种蛤ITS1和ITS2序列中各发现2个与rRNA加工有关的保守区。通过对ITS1和ITS2序列的组装获得了4种蛤5．8SrRNA基因完整序列,序列长度都是157bp,GC含量57．96％～58．60％,4种蛤5．8SrRNA基因相对保守,种间序列差异度0-6．0％,共有10个变异位点,其中转换4处,颠换6处,硬壳蛤和江户布目蛤5．8SrRNA基因序列完全相同。以ITS2序列（包含5．8SrRNA和28SrRNA基因部分序列）为标记,调用北极蛤科的Arctica islandica相应序列数据作外群,构建了帘蛤科贝类的系统发育树,其拓扑结构显示江户布目蛤与硬壳蛤亲缘关系最近,青蛤与其他3物种的亲缘关系最远。     

Study on Sequences of Ribosomal DNA Internal Transcribed Spacers of Clams Belonging to the Veneridae Family (Mollusca: Bivalvia)

The first and second internal transcribed spacer (ITS1 and ITS2) regions of the ribosomal DNA from four species, Meretrix meretrix L., Cyclina sinensis G., Mercenaria mercenaria L., and Protothaca jedoensis L., belonging to the family Veneridae were amplified by PCR and sequenced. The size of the ITS1 PCR amplification product ranged from 663 bp to 978 bp, with GC contents ranging from 60.78% to 64.97%. The size of the ITS1 sequence ranged from 585 bp to 900 bp, which is the largest range reported thus far in bivalve species, with GC contents ranging from 61.03% to 65.62%. The size of the ITS2 PCR amplification product ranged from 513 bp to 644 bp, with GC contents ranging from 61.29% to 62.73%. The size of the ITS2 sequence ranged from 281 bp to 412 bp, with GC contents ranging from 65.21% to 67.87%. Extensive sequence variation and obvious length polymorphisms were noted for both regions in these species, and sequence similarity of ITS2 was higher than that of ITS1 across species. The complete sequences of 5.8S ribosomal RNA gene were obtained by assembling ITS1 and ITS2 sequences, and the sequence length in all species was 157 bp. The phylogenetic tree of Veneridae clams was reconstructed using ITS2-containing partial sequences of both 5.8S and 28S ribosomal DNA as markers and the corresponding sequence information in Arctica islandica as the outgroup. Tree topologies indicated that P. jedoensis shared a close relationship with M. mercenaria and C. sinensis, a distant relationship with other species.     

Fine-scale geographical structure, intra-individual polymorphism and recombination in nuclear ribosomal internal transcribed spacers in Armeria (Plumbaginaceae)

BACKGROUND AND AIMS: Isolation and drift are the main causes for geographic structure of molecular variation. In contrast, the one found in a previous survey in Armeria (Plumbaginaceae) for nuclear ribosomal ITS multicopy regions was species-independent and has been hypothesized to be due to extensive gene-flow and biased concerted evolution. Since this was inferred from a genus-level phylogenetic analysis, the aim of this study was to check for the occurrence of such structure and the validity of the proposed model at a local scale, in a southern Spanish massif (Sierra Nevada), as well as to examine the evolutionary implications at the organism level. METHODS: In addition to 117 sequences of direct PCR products from genomic DNA, 50 sequences of PCR products from cloned DNA were obtained to analyse cases of intragenomic polymorphisms for the ITS regions. KEY RESULTS: Sequence data confirm the occurrence of a species-independent structure at a local scale and reveal insights through the analysis of contact areas between different ITS copies (ribotypes). A comparison between cloned and direct sequences (a) confirms that, within these contact areas, ITS copies co-occur both in different individuals and within single genomes; and (b) reveals recombination between different copies. CONCLUSIONS: This study supports the utility of direct sequences for detecting intra-individual polymorphism and for partially inferring the ITS copies involved, given previous knowledge of the variability. The main evolutionary implication at the organism level is that gene-flow and concerted evolution shape the geographic structure of ITS variation.     

不同产地瓠瓜品种ITS序列的遗传多样性分析

对国产29个瓠瓜也Lagenariasiceraria（Molina）Standl.页品种的ITS序列进行了扩增及测序,并结合引自GenBank的国产9个瓠瓜品种以及国外6个瓠瓜品种和3个同属种类的ITS序列,对它们的ITS序列长度和GC含量以及变异位点进行比较,在此基础上构建系统发育树并对47个样本间的遗传关系进行研究。结果显示：供试47个样本的ITS序列均由ITS1、5.8SrDNA及ITS2组成,各样本间的ITS序列长度、GC含量以及变异位点差异明显。国产38个瓠瓜品种的ITS序列（包括ITS1、5.8SrDNA及ITS2）长度为619-627bp、GC含量为58.00%-63.32%;国外9个样本的ITS序列长度为591-626bp,GC含量为54.17%-63.26%。序列比对结果显示：国产38个瓠瓜品种的ITS序列同源率为84.6%-100.0%,包含221个变异位点;其中,来源于山东的品种‘砧木2’（‘ZhenmuNo.2’）的ITS序列包含的变异位点最多,与其他品种间的同源率也最低。在系统发育树上,国产38个瓠瓜品种可分为3个分支,来源于山东的品种‘砧木2’和来源于河南的品种‘西瓜砧木1’（‘XiguazhenmuNo.1’）各自聚为第1和第2分支;其余36个品种聚为第3分支。而供试的47个样本则可分为2个分支和5个亚组,第1分支可分为2个亚组,包括国产品种‘砧木2’和产自日本的2个品种;第2分支包含的44个样本则进一步分为3个亚组,国产品种‘西瓜砧木1’和产自法国的品种‘白花瓠瓜’（‘White-floweredgourd’）各自聚为第1和第2亚组,其余的42个样本聚为第3亚组。研究结果表明：供试的不同产地瓠瓜品种间存在丰富的遗传变异和地理分化现象,其ITS序列差异与地理分布有一定关系。     

Unique sequence characteristics account for good DGGE separation of almost full-length 18S rDNAs

A major limiting factor for DGGE-based microbial community studies is that the fragments should not be much longer than 500 bp for successful analysis. However, relatively high-resolution was achieved based on DGGE of the long 18S rDNA fragment (>1500 bp), which might be surprising due to the known decrease in DGGE resolution of DNA molecules with large melted regions. A unique sequence characteristic was found in a specific region (ca. 275 bp, named the NS1-end region) of 18S rDNAs, and fungal communities separated from Hong Qu glutinous rice wine brewing system was used to reveal the relationship between high resolution capacity and the unique sequence characteristics. The results showed that DGGE separation of the long 18S rDNA fragments depended on their NS1-end regions. The region is composed of a sequence-variable and short-length GC-poor region (ca. 160 bp) and a GC-rich region (ca. 110 bp), which contribute to the high resolution capacity achieved for DGGE of the long 18S rDNA fragments. Thus DGGE of the long 18S rDNA fragment is recommended as a target fragment for studies of fungal communities whose 18S rDNAs possess similar sequence characteristics. Good resolution and almost full-length 18S rDNA sequences can thus be obtained to provide more accurate and reliable analysis of fungal communities. Since more sequences are obtained directly from the PCR product through the long rDNA fragment approach, this is a convenient and effective approach for sequence-based analysis without using other complementary methods such as an rDNA clone library method.     

A new insight into the phylogeny of vascular cryptogams with special reference to Selaginella and Isoetes inferred from nuclear ITS/5.8S rDNA sequences

In order to provide a better understanding of the evolutionary history of vascular cryptogams, phylogenetic framework was developed based on ITS1, ITS2 and 5.8S rDNA sequences of 102 extant taxa of vascular cryptogams using Maximum Parsimony (MP) analysis. The analysis revealed high GC content in Isoetaceae (60.5 %) in comparison with Selaginellaceae (54.4 %) that was envisaged to be the result of variation in selection, mutational bias, and biased recombination-associated DNA repair within these two plant lineages during evolution. Transition/transversion ratio was observed to be 0.9 in Isoetaceae, 0.68 in Selaginellaceae and 0.57 among all the 102 taxa belonging to lycophytes and ferns. It is hypothesized that the lycophytes have been separated very early during evolution and therefore acquired independent line of evolution from the other plant lineages. Although Selaginellaceae and Isoetaceae are closely related ancient plant groups, pairwise sequence divergence of sampled taxa on the basis of transition and transversion, and disparity index values per site between sampled sequence pairs pointed towards the differential investment of natural selection process. These lead to high rate of nucleotide substitution within nuclear genome of Selaginellaceae with respect to Isoetaceae. MP phylogenetic tree identified Isoetes subinermis, Isoetes durieui and Salvia microphylla as separate group among the studied taxa due to high sequence variation within these species through the time of evolution. Our result interpreted the polyphyletic origin of ferns and provides valuable information regarding the lycophytes and their fern allies.     

Evidence for cryptic hybridization between different evolutionary lineages of the invasive clam genus Corbicula (Veneroida,Bivalvia)

We studied two Corbicula morphotypes in a syntopic population in the Rhine River in order to reveal their taxonomic, reproductive and phylogenetic relationship, using morphometrics, DAF‐fingerprinting, mitochondrial COI and nuclear ITS1 sequence variation. Morphometric analysis showed that two statistically distinguishable morphotypes with few intermediates were present.Mitochondrial sequence analysis detected two divergent clades. DAF‐fingerprinting revealed three highly distinctive multilocus genotypes. Two of the multilocus genotypes were significantly associated with different morphotypes and mitochondrial lineages. The third genotype B, however, was found in both morphotypes, intermediates and mitochondrial lineages. Conclusive evidence for hybridization came from RFLP analysis of the nuclear ITS1 locus. We interpret the hybrids as F1 hybrids between different evolutionary lineages. Integration of Corbicula sequences from all over the world into Maximum Parsimony analysis suggested a simultaneous radiation resulting in several evolutionary lineages whose species status remained doubtful. An unequivocal taxonomic assignment of the two evolutionary lineages in the Rhine population was therefore not possible.     

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.     

Molecular characterization of hard tick <Emphasis Type="Italic">Haemaphysalis longicornis</Emphasis> from China by sequences of the internal transcribed spacers of ribosomal DNA

In the present study, the entire first and second internal transcribed spacer (ITS-1 and ITS-2) regions of nuclear ribosomal DNA (rDNA) of Haemaphysalis longicornis from China were amplified by polymerase chain reaction. The 45 representative amplicons were sequenced, and sequence variation in the ITS was examined. The ITS sequences of H. longicornis were 3644 bp in size, including the part of 18S rDNA, 28S rDNA sequences and the complete ITS-1, 5.8S rDNA and ITS-2 sequences. Sequence analysis revealed that the ITS-1, 5.8S rDNA and ITS-2 of this hard tick were 1582, 152, and 1610 bp in size, respectively. The intra-specific sequence variations of ITS-1 and ITS-2 within H. longicornis were 0–2 and 0–2.2%; however, the inter-specific sequence differences among members of the genus Haemaphysalis were significantly higher, being 35.1–55.2 and 37–52% for ITS-1 and ITS-2, respectively. The molecular approach employed in this study provides the foundation for further studies of the genetic variation of H. longicornis from different hosts and geographical origins in China.     

Sequence analysis of the ribosomal DNA internal transcribed spacer region in some scallop species (Mollusca: Bivalvia: Pectinidae).

The internal transcribed spacer (ITS) region of the ribosomal DNA from the European scallops Aequipecten opercularis, Mimachlamys varia, Hinnites distortus, and Pecten maximus was PCR amplified and sequenced. For each species, three or five clones were examined. The size ranged between 636 and 713 bp (ITS1, 209-276 bp; 5.8S rRNA gene, 157 bp; ITS2, 270-294 bp) and GC content ranged between 47 and 50% (ITS1, 43-49%; 5.8S rRNA gene, 56-57%; ITS2, 44-49%). Variation within repeats was minimal; only clones from M. varia and P. maximus displayed a few variable sites in ITS2. Among scallops, including Chlamys farreri whose ITS sequence appears in databases, significant variation was observed in both ITS1 and ITS2. Phylogenetic analysis using ITS1, ITS2, or both spacer sequences always yielded trees with similar topology. Aequipecten opercularis and P. maximus grouped in one clade and the other three scallops (C. farreri, M. varia, and H. distortus) in another, where M. varia and H. distortus are the more closely related species. These results provide new insights into the evolutionary relationships of scallop species and corroborate the close evolutionary relationship between the tribes Aequipectinini and Pectinini previously deduced from 18S rDNA sequences.     

Preliminary Analysis of Length and GC Content Variation in the Ribosomal First Internal Transcribed Spacer (ITS1) of Marine Animals

Length and guanine–cytosine (GC) content of the ribosomal first internal transcribed spacer (ITS1) were compared across a wide variety of marine animal species, and its phylogenetic utility was investigated. From a total of 773 individuals representing 599 species, we only failed to amplify the ITS1 sequence from 87 individuals by polymerase chain reaction with universal ITS1 primers. No species was found to have an ITS1 region shorter than 100 bp. In general, the ITS1 sequences of vertebrates were longer (318 to 2,318 bp) and richer in GC content (56.8% to 78%) than those of invertebrates (117 to 1,613 bp and 35.8% to 71.3%, respectively). Specifically, gelatinous animals (Cnidaria and Ctenophora) were observed to have short ITS1 sequences (118 to 422 bp) with lower GC content (35.8% to 61.7%) than the other animal taxa. Mollusca and Crustacea were diverse groups with respect to ITS1 length, ranging from 108 to 1,118 and 182 to 1,613 bp, respectively. No universal relationship between length and GC content was observed. Our data indicated that ITS1 has a limited utility for phylogenetic analysis as obtaining confident sequence alignment was often impossible between different genera of the same family and even between congeneric species.     

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.     

Intragenomic Variation and Evolution of the Internal Transcribed Spacer of the rRNA Operon in Bacteria

Variation in the internal transcribed spacer (ITS) of the rRNA (rrn) operon is increasingly used to infer population-level diversity in bacterial communities. However, intragenomic ITS variation may skew diversity estimates that do not correct for multiple rrn operons within a genome. This study characterizes variation in ITS length, tRNA composition, and intragenomic nucleotide divergence across 155 Bacteria genomes. On average, these genomes encode 4.8 rrn operons (range: 2–15) and contain 2.4 unique ITS length variants (range: 1–12) and 2.8 unique sequence variants (range: 1–12). ITS variation stems primarily from differences in tRNA gene composition, with ITS regions containing tRNA-Ala + tRNA-Ile (48% of sequences), tRNA-Ala or tRNA-Ile (10%), tRNA-Glu (11%), other tRNAs (3%), or no tRNA genes (27%). Intragenomic divergence among paralogous ITS sequences grouped by tRNA composition ranges from 0% to 12.11% (mean: 0.94%). Low divergence values indicate extensive homogenization among ITS copies. In 78% of alignments, divergence is <1%, with 54% showing zero variation and 81% containing at least two identical sequences. ITS homogenization occurs over relatively long sequence tracts, frequently spanning the entire ITS, and is largely independent of the distance (basepairs) between operons. This study underscores the potential contribution of interoperon ITS variation to bacterial microdiversity studies, as well as unequivocally demonstrates the pervasiveness of concerted evolution in the rrn gene family. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Reviewing Editor: Dr. Margaret Riley     

Denaturing gradient gel electrophoresis (DGGE): a rapid and sensitive technique to screen nucleotide sequence variation in populations

We describe a rapid and sensitive method for the detection of nucleotide sequence variation that can be used for large-scale screening of population markers. Denaturing gradient gel electrophoresis (DGGE) detects sequence variants of amplified fragments by the differences in their melting behavior. DGGE detects most single-base substitutions when carried out on products amplified with a primer to which a GC clamp has been added. Although DGGE has been primarily used for the detection of limited numbers of single-base mutations in disease studies, it offers great potential for use in population analysis of genetic markers with greater levels of sequence variation. The methodology described was developed to identify the number and distribution of MHC class I alpha 1 alleles among chinook salmon (Oncorhynchus tshawytscha) populations. DGGE detects 28 of 31 identified alpha 1 sequences, which differ by between 1 and 16 nucleotides and a two-codon indel. By creating a network of control alleles, 22-23 of the MHC alleles can be resolved rapidly and accurately by a single gel run condition, and 27 alleles can be resolved by two gel run conditions. This techniques has been used in surveys scoring alleles from two MHC markers (class I alpha 1 and alpha 2) in 20,000 individuals of chinook and coho (O. kisutch) salmon. A single person in our laboratory now analyzes 160 salmon from one MHC locus per day with DGGE.     

The chloroplast genome sequence from <Emphasis Type="Italic">Eugenia uniflora</Emphasis>, a Myrtaceae from Neotropics

Eugenia uniflora is a plant native to tropical America that holds great ecological and economic importance. The complete chloroplast (cp) genome sequence of Eugenia uniflora, a member of the Neotropical Myrtaceae family, is reported here. The genome is 158,445 bp in length and exhibits a typical quadripartite structure of the large (LSC, 87,459 bp) and small (SSC, 18,318 bp) single-copy regions, separated by a pair of inverted repeats (IRs, 26,334 bp). It contains 111 unique genes, including 77 protein-coding genes, 30 tRNAs and 4 rRNAs. The genome structure, gene order, GC content and codon usage are similar to the typical angiosperm cp genomes. Comparison of the entire cp genomes of E. uniflora L. and three other Myrtaceae revealed an expansion of 43 bp in the intergenic spacer located between the IRA/large single-copy (LSC) border and the first gene of LSC region. Simple sequence repeat (SSR) analysis revealed that most SSRs are AT rich, which contribute to the overall AT richness of the cp genome. Additionally, fewer SSRs are distributed in the protein-coding sequences compared to the noncoding regions. Phylogenetic analysis among 58 species based on 57 cp genes demonstrated a closer relationship between E. uniflora L. and Syzygium cumini (L). Skeels compared to the Eucalyptus clade in the Myrtaceae family. The complete cp genome sequence of E. uniflora reported here has importance for population genetics, as well as phylogenetic and evolutionary studies in this species and other Myrtaceae species from Neotropical regions.     

Multiple ITS Haplotypes in the Genome of the Lichenized Basidiomycete Cora inversa (Hygrophoraceae): Fact or Artifact?

The internal transcribed spacer region (ITS) of the nuclear rDNA cistron represents the barcoding locus for Fungi. Intragenomic variation of this multicopy gene can interfere with accurate phylogenetic reconstruction of biological entities. We investigated the amount and nature of this variation for the lichenized fungus Cora inversa in the Hygrophoraceae (Basidiomycota: Agaricales), analyzing base call and length variation in ITS1 454 pyrosequencing data of three samples of the target mycobiont, for a total of 16,665 reads obtained from three separate repeats of the same samples under different conditions. Using multiple fixed alignment methods (PaPaRa) and maximum likelihood phylogenetic analysis (RAxML), we assessed phylogenetic relationships of the obtained reads, together with Sanger ITS sequences from the same samples. Phylogenetic analysis showed that all ITS1 reads belonged to a single species, C. inversa. Pyrosequencing data showed 266 insertion sites in addition to the 325 sites expected from Sanger sequences, for a total of 15,654 insertions (0.94 insertions per read). An additional 3,279 substitutions relative to the Sanger sequences were detected in the dataset, out of 5,461,125 bases to be called. Up to 99.3 % of the observed indels in the dataset could be interpreted as 454 pyrosequencing errors, approximately 65 % corresponding to incorrectly recovered homopolymer segments, and 35 % to carry-forward-incomplete-extension errors. Comparison of automated clustering and alignment-based phylogenetic analysis demonstrated that clustering of these reads produced a 35-fold overestimation of biological diversity in the dataset at the 95 % similarity threshold level, whereas phylogenetic analysis using a maximum likelihood approach accurately recovered a single biological entity. We conclude that variation detected in 454 pyrosequencing data must be interpreted with great care and that a combination of a sufficiently large number of reads per taxon, a set of Sanger references for the same taxon, and at least two runs under different emulsion PCR and sequencing conditions, are necessary to reliably separate biological variation from 454 sequencing errors. Our study shows that clustering methods are highly sensitive to artifactual sequence variation and inadequate to properly recover biological diversity in a dataset, if sequencing errors are substantial and not removed prior to clustering analysis.     

Assessing the potential of regions of the nuclear and mitochondrial genome to develop a "molecular tool box" for the detection and characterization of Phytophthora species

Four different intergenic regions of mitochondrial DNA (mt-IGS), a fragment of the intergenic spacer (IGS) region of the rDNA (rDNA-IGS), and a fragment of the ras-related protein (Ypt1) gene were amplified and sequenced from a panel of 31 Phytophthora species representing the most significant forest pathogens and the breadth of diversity in the genus. Over 80 kbp of novel sequences were generated and alignments showed very variable (introns and non-coding regions) as well as conserved coding regions. The mitochondrial DNA regions had an AT/GC ratio ranging from 67.2 to 89.0% and were appropriate for diagnostic development and phylogeographic analysis. The IGS fragment was less variable but still appropriate to discriminate amongst some important forest pathogens. The introns of the Ypt1 gene were sufficiently polymorphic for the development of molecular markers for almost all Phytophthora species, with more conserved flanking coding regions appropriate for the design of Phytophthora genus-specific primers. In general, phylogenetic analysis of the sequence alignments grouped species in clades that matched those based on the ITS regions of the rDNA. In many cases the resolution was improved over ITS but in other cases sequences were too variable to align accurately and yielded phylograms inconsistent with other data. Key studies on the intraspecific variation and primer specificity remain. However the research has already yielded an enormous dataset for the identification, detection and study of the molecular evolution of Phytophthora species.