Phylogenetic analysis of the nuclear alcohol dehydrogenase (Adh) gene family in Carex section Acrocystis (Cyperaceae) and combined analyses of Adh and nuclear ribosomal ITS and ETS sequences for inferring species relationships

We analyzed sequence variation for the alcohol dehydrogenase (Adh) gene family in Carex section Acrocystis (Cyperaceae) to reconstruct Adh gene trees for Acrocystis species and to characterize the structure of the Adh gene family in Carex. Two Adh loci were included with ITS and ETS sequences in a combined Bayesian inference analysis of Carex section Acrocystis to gain a better understanding of species relationships in the section. In addition, we comment on how the results presented here contribute to our knowledge of the birth-death process of the Adh gene family in angiosperms. It appears that the structure of the Adh gene family in Carex is complex with possibly six loci present in the gene family. Additionally, variation among Acrocystis species within loci is quite low, and there is little phylogenetic resolution in the individual datasets. Bayesian inference analysis of the combined ITS, ETS, Adh1, and Adh2 datasets resulted in a moderately well-supported phylogenetic hypothesis of relationships in the section which is discussed in relation to previous hypotheses of relationships.     

Molecular characterization of pouched amphistome parasites (Trematoda: Gastrothylacidae) using ribosomal ITS2 sequence and secondary structures

Members of the family Gastrothylacidae (Trematoda: Digenea: Paramphistomata) are parasitic in ruminants throughout Africa and Asia. In north-east India, five species of pouched amphistomes, namely Fischoederius cobboldi, F. elongatus, Gastrothylax crumenifer, Carmyerius spatiosus and Velasquezotrema tripurensis, belonging to this family have been reported so far. In the present study, the molecular phylogeny of these five gastrothylacid species is derived using the second internal transcribed spacer (ITS2) sequence and secondary structure analyses. ITS2 sequence analysis was carried out to see the occurrence of interspecific variations among the species. Phylogenetic analyses were performed for primary sequence data alone as well as the combined sequence-structure information using neighbour-joining and Bayesian approaches. The sequence analysis revealed that there exist considerable interspecific variations among the various gastrothylacid fluke species. In contrast, the inferred secondary structures for the five species using minimum free energy modelling showed structural identities, in conformity with the core four-helix domain structure that has been recently identified as common to almost all eukaryotic taxa. The phylogenetic tree reconstructed using combined sequence-structure data showed a better resolution, as compared to the one using sequence data alone, with the gastrothylacid species forming a monophyletic group that is well separated from members of the other family, Paramphistomidae, of the amphistomid flukes group. The study provides the molecular characterization based on primary sequence data of the rDNA ITS2 region of the gastrothylacid amphistome flukes. Results also demonstrate the phylogenetic utility of the ITS2 sequence-secondary structure data for inferences at higher taxonomic levels.     

The ITS2 Database

The internal transcribed spacer 2 (ITS2) has been used as a phylogenetic marker for more than two decades. As ITS2 research mainly focused on the very variable ITS2 sequence, it confined this marker to low-level phylogenetics only. However, the combination of the ITS2 sequence and its highly conserved secondary structure improves the phylogenetic resolution¹ and allows phylogenetic inference at multiple taxonomic ranks, including species delimitation^2-8.The ITS2 Database⁹ presents an exhaustive dataset of internal transcribed spacer 2 sequences from NCBI GenBank¹¹ accurately reannotated¹⁰. Following an annotation by profile Hidden Markov Models (HMMs), the secondary structure of each sequence is predicted. First, it is tested whether a minimum energy based fold¹² (direct fold) results in a correct, four helix conformation. If this is not the case, the structure is predicted by homology modeling¹³. In homology modeling, an already known secondary structure is transferred to another ITS2 sequence, whose secondary structure was not able to fold correctly in a direct fold.The ITS2 Database is not only a database for storage and retrieval of ITS2 sequence-structures. It also provides several tools to process your own ITS2 sequences, including annotation, structural prediction, motif detection and BLAST¹⁴ search on the combined sequence-structure information. Moreover, it integrates trimmed versions of 4SALE^15,16 and ProfDistS¹⁷ for multiple sequence-structure alignment calculation and Neighbor Joining¹⁸ tree reconstruction. Together they form a coherent analysis pipeline from an initial set of sequences to a phylogeny based on sequence and secondary structure.In a nutshell, this workbench simplifies first phylogenetic analyses to only a few mouse-clicks, while additionally providing tools and data for comprehensive large-scale analyses.     

Coevolution of Retrovarium n. gen. (Digenea: Cryptogonimidae) in Lutjanidae and Haemulidae (Perciformes) in the Indo-West Pacific

We describe 11 new species of cryptogonimids belonging to Retrovarium n. gen., from eight species of Lutjanidae and one species of Haemulidae, from the Great Barrier Reef, French Polynesia and the Maldives. We also transfer Neoparacryptogonimus saccatus (Manter, 1963) and Neoparacryptogonimus sphericus Nahhas, Sey & Nishimoto, 1998 to Retrovarium. The morphologically based taxonomic approach was augmented with DNA sequence data from three nuclear ribosomal DNA regions (28S, ITS1 and ITS2) to explore the species integrity, biogeographic distribution and evolution of the species recognised here. Sequencing included multiple replicates and revealed 11 distinct genotypes which corroborated our morphologically based hypotheses of putative species present in the system. There was no intraspecific variation and all three rDNA regions differed between every combination of species. Two species exhibited wide geographic ranges, having identical rDNA sequences between the Great Barrier Reef and the Maldives, localities separated by over 9600 km. One host species, Symphorus nematophorus, proved to be exceptionally rich, harbouring six species. Minimum evolution analyses were conducted on each of the rDNA datasets independently; minimum evolution, maximum likelihood and Bayesian inference analyses were conducted on a combined (28S, ITS1 and ITS2) dataset for sequence comparison purposes and to explore the evolutionary history of these parasites. To examine the coevolutionary history of this complex, assessment of phylogenetic relationships between the 23 species of Lutjanidae and two species of Haemulidae collected during this survey was performed with data from 16S and cytochrome b mtDNA using Bayesian inference analysis. Despite the high host specificity observed in most of the species, mapping of the parasites on the host phylogeny revealed an absence of strict coevolution or co-descent within this complex. Overall, Retrovarium appears to have had an exceptionally patchy radiation, failing to infect many taxa, infecting species with no readily discernible pattern, and radiating dramatically within one species.     

Molecular phylogenetics and evolution of the endemic Hawaiian genus Adenophorus (Grammitidaceae)

Recent studies of the phylogeny of several groups of native Hawaiian vascular plants have led to significant insights into the origin and evolution of important elements of the Hawaiian flora. No groups of Hawaiian pteridophytes have been subjected previously to rigorous phylogenetic analysis. We conducted a molecular phylogenetic analysis of the endemic Hawaiian fern genus Adenophorus employing DNA sequence variation from three cpDNA fragments: rbcL, atpbeta, and the trnL-trnF intergenic spacer (IGS). In the phylogenetic analyses we employed maximum parsimony and Bayesian inference. Bayesian phylogenetic inference often provided stronger support for hypothetical relationships than did nonparametric bootstrap analyses. Although phylogenetic analyses of individual DNA fragments resulted in different patterns of relationships among species and varying levels of support for various clades, a combined analysis of all three sets of sequences produced one, strongly supported phylogenetic hypothesis. The primary features of that hypothesis are: (1) Adenophorus is monophyletic; (2) subgenus Oligadenus is paraphyletic; (3) the enigmatic endemic Hawaiian species Grammitis tenella is strongly supported as the sister taxon to Adenophorus; (4) highly divided leaf blades are evolutionarily derived in the group and simple leaves are ancestral; and, (5) the biogeographical origin of the common ancestor of the Adenophorus-G. tenella clade remains unresolved, although a neotropical origin seems most likely.     

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.     

A common core of secondary structure of the internal transcribed spacer 2 (ITS2) throughout the Eukaryota

The ongoing characterization of novel species creates the need for a molecular marker which can be used for species- and, simultaneously, for mega-systematics. Recently, the use of the internal transcribed spacer 2 (ITS2) sequence was suggested, as it shows a high divergence in sequence with an assumed conservation in structure. This hypothesis was mainly based on small-scale analyses, comparing a limited number of sequences. Here, we report a large-scale analysis of more than 54,000 currently known ITS2 sequences with the goal to evaluate the hypothesis of a conserved structural core and to assess its use for automated large-scale phylogenetics. Structure prediction revealed that the previously described core structure can be found for more than 5000 sequences in a wide variety of taxa within the eukaryotes, indicating that the core secondary structure is indeed conserved. This conserved structure allowed an automated alignment of extremely divergent sequences as exemplified for the ITS2 sequences of a ctenophorean eumetazoon and a volvocalean green alga. All classified sequences, together with their structures can be accessed at http://www.biozentrum.uni-wuerzburg.de/bioinformatik/projects/ITS2.html. Furthermore, we found that, although sample sequences are known for most major taxa, there exists a profound divergence in coverage, which might become a hindrance for general usage. In summary, our analysis strengthens the potential of ITS2 as a general phylogenetic marker and provides a data source for further ITS2-based analyses.     

Homology modeling revealed more than 20,000 rRNA internal transcribed spacer 2 (ITS2) secondary structures

Structural genomics meets phylogenetics and vice versa: Knowing rRNA secondary structures is a prerequisite for constructing rRNA alignments for inferring phylogenies, and inferring phylogenies is a precondition to understand the evolution of such rRNA secondary structures. Here, both scientific worlds go together. The rRNA internal transcribed spacer 2 (ITS2) region is a widely used phylogenetic marker. Because of its high variability at the sequence level, correct alignments have to take into account structural information. In this study, we examine the extent of the conservation in structure. We present (1) the homology modeled secondary structure of more than 20,000 ITS2 covering about 14,000 species; (2) a computational approach for homology modeling of rRNA structures, which additionally can be applied to other RNA families; and (3) a database providing about 25,000 ITS2 sequences with their associated secondary structures, a refined ITS2 specific general time reversible (GTR) substitution model, and a scoring matrix, available at http://its2.bioapps.biozentrum.uni-wuerzburg.de.     

Systematics of the subfamily Haplorchiinae (Trematoda: Heterphyidae), based on nuclear ribosomal DNA genes and ITS2 region

Phylogenetic relationships of 6 species in the trematode subfamily Haplorchiinae were analyzed using small and large subunit of ribosomal DNA genes (18S rDNA and 28S rDNA) and internal transcribed spacer subunit II (ITS2) region as molecular markers. Maximum Likelihood and Bayesian inference analyses of combined rDNAs and ITS2 indicated a close relationship between the genera Haplorchis and Procerovum, while these two genera were distinct from Stellantchasmus falcatus. These phylogenetic relationships were consistent with the number of testes but not with the characters of the modification of the seminal vesicle or of the ventral sucker. Although three Haplorchis spp. were, together with Procerovum, in the same cluster, their mutual topology was incongruent between rDNA and ITS2 trees. Phylogenetic analyses using other molecular markers with more species are necessary to work out solid phylogenetic relationships among the species in this subfamily.     

Utility of low-copy nuclear markers in phylogenetic reconstruction of Hypericum L. (Hypericaceae)

Primers and sequence variation for two low-copy nuclear genes (LCG) not previously used for phylogenetic inference in the genus Hypericum, PHYC and EMB2765, are presented here in comparison with the fast-evolving nuclear intergenic spacer ITS. Substitution rates in the LCG markers were half those reported in ITS for Hypericum, which might help avoid the problems caused by substitution saturation and difficulties to establish homologies that afflict the latter marker. We included representatives of all major clades within Hypericum and found that levels of phylogenetic resolution, clade support values and internal character consistency were similar to, or even higher than, those of ITS-based phylogenies. The presence of at least two copies in EMB2765 in Hypericum imposed a methodological challenge that was circumvented by the design of an effective clade-specific primer. Both EMB2765 and, especially, PHYC appear to be good alternatives to the ITS marker, confirming the main phylogenetic relationships found in previous studies, but with improved resolution and support values for some basal relationships.     

Evolution and phylogenetic information content of the ribosomal DNA repeat unit in the Blattodea (Insecta)

The organization, structure, and nucleotide variability of the ribosomal repeat unit was compared among families, genera, and species of cockroaches (Insecta:Blattodea). Sequence comparisons and molecular phylogenetic analyses were used to describe rDNA repeat unit variation at differing taxonomic levels. A reverse similar 1200 bp fragment of the 28S rDNA sequence was assessed for its potential utility in reconstructing higher-level phylogenetic relationships in cockroaches. Parsimony and maximum likelihood analyses of these data strongly support the expected pattern of relationships among cockroach groups. The examined 5' end of the 28S rDNA is shown to be an informative marker for larger studies of cockroach phylogeny. Comparative analysis of the nucleotide sequences of the rDNA internal transcribed spacers (ITS1 and ITS2) among closely related species of Blattella and Periplaneta reveals that ITS sequences can vary widely in primary sequence, length, and folding pattern. Secondary structure estimates for the ITS region of Blattella species indicate that variation in this spacer region can also influence the folding pattern of the 5.8S subunit. These results support the idea that ITS sequences play an important role in the stability and function of the rRNA cluster.     

Phylogeny reconstruction in the gap between too little and too much divergence: the closest relatives of Senecio jacobaea (Asteraceae) according to DNA sequences and AFLPs

The species composition and molecular phylogeny of Senecio sect. Jacobaea (Asteraceae; Senecioneae) were studied to identify the closest relatives of Senecio jacobaea. Maximum parsimony and Bayesian inference analyses of DNA sequence data of the plastid (the trnT-L igs, the trnL intron, two parts of the trnK intron, and the psbA-trnH igs) and nuclear genome (ITS1, 5.8S, and ITS2) showed these markers to be suitable to assess the species composition of sect. Jacobaea, identifying 24 species as members of this section. Of these, nine species were not previously assigned to the section. The selected DNA sequence regions, however, showed too little sequence divergence to be optimal for phylogenetic inference within sect. Jacobaea. In contrast, AFLPs proved to be too variable to be used to study relationships between the basal lineages in sect. Jacobaea. Nonetheless, these markers are very useful to study the phylogeny of S. jacobaea and its closest relatives. The combined use of DNA sequence data and AFLPs allowed us to take a major step towards resolving phylogenetic relationships in sect. Jacobaea, identifying Senecio alpinus, Senecio pancicii and Senecio subalpinus (using DNA sequence data) or Senecio chrysanthemoides (using AFLPs) as the closest relatives of S. jacobaea.     

Phylogenetic analysis of Lappula Moench (Boraginaceae) based on molecular and morphological data

Lappula Moench includes ca. 70 species and exhibits a wide range of nutlet variation. Currently, the evolutionary relationships among species of Lappula have not been examined; therefore, to elucidate phylogenetic relationships and morphological evolution within Lappula and related genera, we conducted phylogenetic analyses with matrices that include 48 species as well as four DNA regions (ITS, trnL-trnF, rpS16 and psbA-trnH) and 18 morphological characters. These matrices were analyzed using maximum parsimony, maximum likelihood, and Bayesian inference methods. Analyses of the combined molecular and morphological data result in a phylogeny that is better resolved than that based solely on molecular sequence data. Phylogenetic results suggest that the current infrageneric classification of Lappula, at least at the subsectional and series level, is artificial. The evolutionary patterns of 18 morphological characters are investigated in a phylogenetic context. In Lappula, nutlet homomorphism and small corollas are resolved as ancestral, while nutlet heteromorphism and larger corollas are derived.     

Application of the secondary structure model of rRNA for phylogeny: D2-D3 expansion segments of the LSU gene of plant-parasitic nematodes from the family Hoplolaimidae Filipjev, 1934

Knowledge of rRNA structure is increasingly important to assist phylogenetic analysis through reconstructing optimal alignment, utilizing molecule features as an additional source of data and refining appropriate models of evolution of the molecule. We describe a procedure of optimization for alignment and a new coding method for nucleotide sequence data using secondary structure models of the D2 and D3 expansion fragments of the LSU-rRNA gene reconstructed for fifteen nematode species of the agriculturally important and diverse family Hoplolaimidae, order Tylenchida. Using secondary structure information we converted the original sequence data into twenty-eight symbol codes and submitted the transformed data to maximum parsimony analysis. We also applied the original sequence data set for Bayesian inference. This used the doublet model with sixteen states of nucleotide doublets for the stem region and the standard model of DNA substitution with four nucleotide states for loops and bulges. By this approach, we demonstrate that using structural information for phylogenetic analyses led to trees with lower resolved relationships between clades and likely eliminated some artefactual support for misinterpreted relationships, such as paraphyly of Helicotylenchus or Rotylenchus. This study as well as future phylogenetic analyses is herein supported by the development of an on-line database, NEMrRNA, for rRNA molecules in a structural format for nematodes. We also have developed a new computer program, RNAstat, for calculation of nucleotide statistics designed and proposed for phylogenetic studies.     

Phylogenetic analysis of Myobia musculi (Schranck, 1781) by using the 18S small ribosomal subunit sequence

We used high-fidelity PCR to amplify 2 overlapping regions of the ribosomal gene complex from the rodent fur mite Myobia musculi. The amplicons encompassed a large portion of the mite's ribosomal gene complex spanning 3128 nucleotides containing the entire 18S rRNA, internal transcribed spacer (ITS) 1,5.8S rRNA, ITS2, and a portion of the 5'-end of the 28S rRNA. M. musculi's 179-nucleotide 5.8S rRNA nucleotide sequence was not conserved, so this region was identified by conservation of rRNA secondary structure. Maximum likelihood and Bayesian inference phylogenetic analyses were performed by using multiple sequence alignment consisting of 1524 nucleotides of M. musculi 18S rRNA and homologous sequences from 42 prostigmatid mites and the tick Dermacentor andersoni. The phylograms produced by both methods were in agreement regarding terminal, secondary, and some tertiary phylogenetic relationships among mites. Bayesian inference discriminated most infraordinal relationships between Eleutherengona and Parasitengona mites in the suborder Anystina. Basal relationships between suborders Anystina and Eupodina historically determined by comparing differences in anatomic characteristics were less well-supported by our molecular analysis. Our results recapitulated similar 18S rRNA sequence analyses recently reported. Our study supports M. musculi as belonging to the suborder Anystina, infraorder Eleutherenona, and superfamily Cheyletoidea.     

Sequence variation in mitochondrial ATP synthase subunit 6 & 8 and nuclear genes ITS1 and ITS2 in 17 cyprinid species

Sequence variation for the mitochondrial ATP synthase subunit 6 & 8 and the nuclear ITS1 marker was investigated for 17 cyprinid fish species in order to increase the number of possible loci for character sampling. The mitochondrial locus provided appropriate information for further phylogenetic studies, which can be readily applied to other taxa; the nuclear ITS1, however, cannot be recommended for inference of relationships among different taxa. No congruence was obtained between the trees resulting from the mitochondrial and nuclear gene using Bayesian Inference.     

The rpb2 gene represents a viable alternative molecular marker for the analysis of environmental fungal communities

Although the commonly used internal transcribed spacer region of rDNA (ITS) is well suited for taxonomic identification of fungi, the information on the relative abundance of taxa and diversity is negatively affected by the multicopy nature of rDNA and the existence of ITS paralogues. Moreover, due to high variability, ITS sequences cannot be used for phylogenetic analyses of unrelated taxa. The part of single‐copy gene encoding the second largest subunit of RNA polymerase II (rpb2) was thus compared with first spacer of ITS as an alternative marker for the analysis of fungal communities in spruce forest topsoil, and their applicability was tested on a comprehensive mock community. In soil, rpb2 exhibited broad taxonomic coverage of the entire fungal tree of life including basal fungal lineages. The gene exhibited sufficient variation for the use in phylogenetic analyses and taxonomic assignments, although it amplifies also paralogues. The fungal taxon spectra obtained with rbp2 region and ITS1 corresponded, but sequence abundance differed widely, especially in the basal lineages. The proportions of OTU counts and read counts of major fungal groups were close to the reality when rpb2 was used as a molecular marker while they were strongly biased towards the Basidiomycota when using the ITS primers ITS1/ITS4. Although the taxonomic placement of rbp2 sequences is currently more difficult than that of the ITS sequences, its discriminative power, quantitative representation of community composition and suitability for phylogenetic analyses represent significant advantages.     

Toward a better understanding of the infrageneric relationships in Cortinarius (Agaricales, Basidiomycota)

Research on the molecular systematics of Cortinarius, a species-rich mushroom genus with nearly global distribution, is just beginning. The present study explores infrageneric relationships using rDNA ITS and LSU sequence data. One large dataset of 132 rDNA ITS sequences and one combined da-taset with 54 rDNA ITS and LSU sequences were generated. Hebeloma was used as outgroup. Bayesian analyses and maximum-likelihood (ML) analyses were carried out. Bayesian phylogenetic inference performed equally well or better than ML, especially in large datasets. The phylogenetic analysis of the combined dataset with species representing all currently recognized subgenera recovered seven well-supported clades (Bayesian posterior probabilities BPP > 90%). These major clades are: /Myxacium s.l., /subg. Cortinarius, the /phlegmacioid clade (including the subclades /Phlegmacium and /Delibuti), the /calochroid clade (/Calochroi, /Ochroleuci and /Allutus), the /telamonioid clade (/Telamonia, /Orellani, /Anomali), /Dermocybe s.l. and /Myxotelamonia. Our results show that Cortinarius consists of many lineages, but the relationships among these clades could not be elucidated. On one hand, the low divergence in rDNA sequences can be held responsible for this; on the other hand, taxon sampling is problematic in Cortinarius phylogeny. Because of the incredibly high diversity (~2000 Cortinarius species), our sampling included <5% of the known species. By choosing type species of subgenera and sections, our sampling is strongly biased toward Northern Hemisphere taxa. More extensive taxon sampling, especially of species from the Southern Hemisphere, is essential to resolve the phylogeny of this important genus of ectomycorrhizal fungi.     

Molecular identification and phylogenetic analysis of nuclear rDNA sequences among three opisthorchid liver fluke species (Opisthorchiidae: Trematoda)

In this study, we describe the development of a fast and accurate molecular identification system for human-associated liver fluke species (Opisthorchis viverrini, Opisthorchis felineus, and Clonorchis sinensis) using the PCR-RFLP analysis of the 18S-ITS1-5.8S nuclear ribosomal DNA region. Based on sequence variation in the target rDNA region, we selected three species-specific restriction enzymes within the ITS1 regions, generating different restriction profiles among the species: MunI for O. viverrini, NheI for O. felineus, and XhoI for C. sinensis, respectively. Each restriction enzyme generated different-sized fragments specific to the species examined, but no intraspecific polymorphism or cross-reaction between the species was detected in their restriction pattern. These results indicate that PCR-linked restriction analysis of the ITS1 region allows for the rapid and reliable molecular identification among these opisthorchid taxa. In addition, phylogenetic analysis of rDNA sequences using different methods (MP, ML, NJ, and Bayesian inference) displayed O. viverrini and O. felineus as a sister group, but this relationship was not strongly supported. The failure of recovering a robust phylogeny may be due to the relatively small number of synapomorphic characters shared among the species, yielding weak phylogenetic signal. Alternatively, rapid speciation within a very short period time could be another explanation for the relatively poorly resolved relationships among these species. Our data are insufficient for discriminating between sudden cladogenesis and other potential causes of poor resolution. Further information from independent loci might help resolve this phylogeny.