Two mitogenomes in Gruiformes (<Emphasis Type="Italic">Amaurornis akool</Emphasis>/<Emphasis Type="Italic">A. phoenicurus</Emphasis>) and the phylogenetic placement of Rallidae

Rallidae, with 34 genera including 142 species, is the largest family in the Gruiformes, the phylogenetic placement of this family was still in debate. The complete mitochondrial genomes (mitogenomes), with many advantageous characters, have become popular markers in phylogenetic analyses. We sequenced the mitogenomes of brown crake (Amaurornis akool) and white-breasted waterhen (Amaurornis phoenicurus), analyzed the genomic characters of mitogenomes in Rallidae, and explored the phylogenetic relationships between Rallidae and other four families in Gruiformes based on mitogenome sequences of 32 species with Bayesian method. The mitogenome of A. akool/A. phoenicurus was 16,950/17,213 bp in length, and contained 37 genes typical to avian mitogenomes and one control region, respectively. The genomic characters of mitogenomes in Rallidae were similar. The phylogenetic results indicated that, among five families, Rallidae had closest relationship with Heliornithidae, which formed a sister taxa to Gruidae, while Rhynochetidae located in the basal lineage. Within Rallidae, Rallina was ancestral clade. Gallirallus & Rallus and Aramides were closely related, Gallicrex & Amaurornis and Fulica & Gallinula had close relationships, and these two taxa formed a sister clade to Porphyrio & Coturnicops. Our phylogenetic analyses provided solid evidence for the phylogenetic placement of Rallidae and the evolutionary relationships among different genus within this family. In addition, the mitogenome data presented here provide useful information for further molecular systematic investigations on Gruiformes as well as conservation biology research of these species.     

Evaluation of the gene encoding translation termination factor eRF3 as a possible phylogenetic marker

The region of the nuclear GSPT2 gene coding for the N and M domains of translation termination factor eRF3b was tested in Rodentia for applicability as a new molecular marker. It cannot be used as a phylogenetic marker at the intrageneric level because of insufficient variability within families and the impossibility of resolving relationships in the family Cricetidae. However, this GSPT2 region allows reliable identification of higher taxa. The phylogenetic relationships among families revealed with the proposed molecular marker is generally in agreement with current concepts. The new marker indicates a close relationship between the genus Acomys and the family Gerbillidae, which is in agreement with other molecular data but contradicts morphological data. Thus, the region of the nuclear GSPT2 gene encoding the N and M domains of eRF3b can serve as an adequate phylogenetic marker in placental mammals at the level of families or higher taxa. It can also be used in solving controversial questions of phylogeny and taxonomy.     

Genome-Wide Identification of the Dicer-Like,Argonaute, and RNA-Dependent RNA Polymerase Gene Families in Cucumber (<Emphasis Type="Italic">Cucumis sativus</Emphasis> L.)

Dicer, Argonaute (AGO), and RNA-dependent RNA polymerase (RDR) comprise the core components of RNA-induced silencing complexes, which trigger RNA silencing. Here, we performed a complete analysis of the cucumber Dicer-like, AGO, and RDR gene families including the gene structure, genomic localization, and phylogenetic relationships among family members. We identified seven CsAGO genes, five CsDCL genes, and eight CsRDR genes in cucumber. Based on phylogenetic analysis, each of these genes families was categorized into three or four clades. The orthologs of CsAGOs, CsDCLs, and CsRDRs were identified in apple, peach, wild strawberry, foxtail millet, and maize, and the evolutionary relationships among the orthologous gene pairs were investigated. We also investigated the expression levels of CsAGOs, CsDCLs, and CsRDRs in various cucumber tissues. All CsAGOs were relatively higher upregulated in leaves and tendrils than in other organs, especially CsAGO1c, CsAGO1d, and CsAGO7. All CsDCL genes were relatively higher upregulated in tendrils, with almost no expression detected for CsDCL1, CsDCL4a, or CsDCL4b in other organs. In addition, CsRDR1a, CsRDR2, CsRDR3, and CsRDR6 had relatively higher upregulation in tendrils, whereas almost all CsRDRs were downregulation in other organs. The results of this study will facilitate further studies of gene silencing pathways in cucumber.     

Genome-wide dissection of the chalcone synthase gene family in <Emphasis Type="Italic">Oryza sativa</Emphasis>

Enzymes of the chalcone synthase (CHS) family catalyze the generation of multiple secondary metabolites in fungi, plants, and bacteria. These metabolites have played key roles in antimicrobial activity, UV protection, flower pigmentation, and pollen fertility during the evolutionary process of land plants. We performed a genome-wide investigation about CHS genes in rice (Oryza sativa). The phylogenetic relationships, gene structures, chromosomal locations, and functional predictions of the family members were examined. Twenty-seven CHS family genes (OsCHS01–27) were identified in the rice genome and were found to cluster into six classes according to their phylogenetic relationships. The 27 OsCHS genes were unevenly distributed on six chromosomes, and 17 genes were found in the genome duplication zones with two segmental duplication and five tandem duplication events that may have played key roles in the expansion of the rice CHS gene family. In addition, the OsCHS genes exhibited diverse expression patterns under salicylic acid treatment. Our results revealed that the OsCHS genes exhibit both diversity and conservation in many aspects, which will contribute to further studies of the function of the rice CHS gene family and provide a reference for investigating this family in other plants.     

Two new species of the <Emphasis Type="Italic">Fusarium solani</Emphasis> species complex isolated from compost and hibiscus (<Emphasis Type="Italic">Hibiscus</Emphasis> sp.)

Two new species in the Fusarium solani species complex (FSSC) are described and introduced. The new taxa are represented by German isolates CBS 142481 and CBS 142480 collected from commercial yard waste compost and vascular tissue of a wilting branch of hibiscus, respectively. The phylogenetic relationships of the collected strains to one another and within the FSSC were evaluated based on DNA sequences of 6 gene loci. Due to the limited sequence data available for reference strains in GenBank, however, a multi-gene phylogenetic analysis included partial sequences for the internal transcribed spacer region and intervening 5.8S nrRNA gene (ITS), translation elongation factor 1-alpha (tef1) and the RNA polymerase II second largest subunit (rpb2). Morphological and molecular phylogenetic data independently showed that these strains are distinct populations of the FSSC, nested within Clade 3. Thus, we introduce Fusarium stercicola and Fusarium witzenhausenense as novel species in the complex. In addition, 19 plant species of 7 legume genera were evaluated for their potential to host the newly described taxa. Eighteen plant species were successfully colonized, with 6 and 9 of these being symptomatic hosts for F. stercicola and F. witzenhausenense, respectively. As plants of the family Fabaceae are very distant to the originally sourced material from which the new taxa were recovered, our results suggest that F. stercicola and F. witzenhausenense are not host-specific and are ecologically fit to sustain stable populations in variety of habitats.     

Evolution of blue-flowered species of genus <Emphasis Type="Italic">Linum</Emphasis> based on high-throughput sequencing of ribosomal RNA genes

Background

The species relationships within the genus Linum have already been studied several times by means of different molecular and phylogenetic approaches. Nevertheless, a number of ambiguities in phylogeny of Linum still remain unresolved. In particular, the species relationships within the sections Stellerolinum and Dasylinum need further clarification. Also, the question of independence of the species of the section Adenolinum still remains unanswered. Moreover, the relationships of L. narbonense and other species of the section Linum require further clarification. Additionally, the origin of tetraploid species of the section Linum (2n?=?30) including the cultivated species L. usitatissimum has not been explored. The present study examines the phylogeny of blue-flowered species of Linum by comparisons of 5S rRNA gene sequences as well as ITS1 and ITS2 sequences of 35S rRNA genes.

Results

High-throughput sequencing has been used for analysis of multicopy rRNA gene families. In addition to the molecular phylogenetic analysis, the number and chromosomal localization of 5S and 35S rDNA sites has been determined by FISH.Our findings confirm that L. stelleroides forms a basal branch from the clade of blue-flowered flaxes which is independent of the branch formed by species of the sect. Dasylinum. The current molecular phylogenetic approaches, the cytogenetic analysis as well as different genomic DNA fingerprinting methods applied previously did not discriminate certain species within the sect. Adenolinum. The allotetraploid cultivated species L. usitatissimum and its wild ancestor L. angustifolium (2n?=?30) could originate either as the result of hybridization of two diploid species (2n?=?16) related to the modern L. gandiflorum and L. decumbens, or hybridization of a diploid species (2n?=?16) and a diploid ancestor of modern L. narbonense (2n?=?14).

Conclusions

High-throughput sequencing of multicopy rRNA gene families allowed us to make several adjustments to the phylogeny of blue-flowered flax species and also reveal intra- and interspecific divergence of the rRNA gene sequences.

     

Phylogenetic relationships in the genus <Emphasis Type="Italic">Primula</Emphasis> L. (Primulaceae) inferred from the ITS region sequences of nuclear rDNA

The nucleotide sequences of the nuclear rDNA ITS regions were determined for 34 species of the genus Primula L. and one species of the genus Cortusa L., family Primulaceae Vent., and used to infer the phylogenetic relationships among these species. In this analysis species of the Russian flora and the flora of adjacent territories were studied for the first time. The results clarified the taxomic structure of the genus Primula and confirmed the entity of some of its sections; but not the subgenera sensu J.Richards. Our data do not support an independent status of the genus Cortusa, placing it as one of the terminal lineages of the section Cortusoides Balf. f. in the genus Primula.     

Comparative genomic analysis of the <Emphasis Type="Italic">Lipase3</Emphasis> gene family in five plant species reveals distinct evolutionary origins

Lipases are physiologically important and ubiquitous enzymes that share a conserved domain and are classified into eight different families based on their amino acid sequences and fundamental biological properties. The Lipase3 family of lipases was reported to possess a canonical fold typical of α/β hydrolases and a typical catalytic triad, suggesting a distinct evolutionary origin for this family. Genes in the Lipase3 family do not have the same functions, but maintain the conserved Lipase3 domain. There have been extensive studies of Lipase3 structures and functions, but little is known about their evolutionary histories. In this study, all lipases within five plant species were identified, and their phylogenetic relationships and genetic properties were analyzed and used to group them into distinct evolutionary families. Each identified lipase family contained at least one dicot and monocot Lipase3 protein, indicating that the gene family was established before the split of dicots and monocots. Similar intron/exon numbers and predicted protein sequence lengths were found within individual groups. Twenty-four tandem Lipase3 gene duplications were identified, implying that the distinctive function of Lipase3 genes appears to be a consequence of translocation and neofunctionalization after gene duplication. The functional genes EDS1, PAD4, and SAG101 that are reportedly involved in pathogen response were all located in the same group. The nucleotide diversity (Dxy) and the ratio of nonsynonymous to synonymous nucleotide substitutions rates (Ka/Ks) of the three genes were significantly greater than the average across the genomes. We further observed evidence for selection maintaining diversity on three genes in the Toll-Interleukin-1 receptor type of nucleotide binding/leucine-rich repeat immune receptor (TIR-NBS LRR) immunity-response signaling pathway, indicating that they could be vulnerable to pathogen effectors.     

Genotypic characterization of Lactic acid bacteria in gut microbiome of freshwater fish

The present study focused on identification and genotypic characterization of Lactic acid bacteria (LAB) in the intestine of freshwater fish. 76 strains of LAB were isolated and identified by 16S rRNA gene sequences and hsp60 gene sequences as different strains of Lactobacillus plantarum, Lactobacillus pentosus, Lactobacillus fermentum, Lactobacillus delbrueckii subsp. bulgaricus, Lactobacillus brevis, Lactobacillus reuteri, Lactobacillus salivarius, Pediococcus pentosaceus, Pediococcus acidilactici, Weissella paramesenteroides, Weissella cibaria, Enterococcus faecium, and Enterococcus durans. The hsp60 gene showed a higher level of sequence variation among the isolates examined, with lower interspecies sequence similarity providing more resolutions at the species level than the 16S rRNA gene. Phylogenetic tree derived from hsp60 gene sequences with higher bootstrap values at the nodal branches was more consistent as compared to phylogenetic tree constructed from 16S rRNA gene sequences. Closely related species L. plantarum and L. pentosus as well as species L. delbrueckii subsp. bulgaricus and L. fermentum were segregated in different cluster in hsp60 phylogenetic tree whereas such a distribution was not apparent in 16S rRNA phylogenetic tree. In silico restriction analysis revealed a high level of polymorphism within hsp60 gene sequences. Restriction pattern with enzymes AgsI and MseI in hsp60 gene sequences allowed differentiation of all the species including closely related species L. plantarum and L. pentosus, E. faecium and E. durans. In general, hsp60 gene with higher evolutionary divergence proved to be a better phylogenetic marker for the group LAB.     

Filling the gaps in the classification of the Digenea Carus, 1863: systematic position of the Proterodiplostomidae Dubois, 1936 within the superfamily Diplostomoidea Poirier, 1886, inferred from nuclear and mitochondrial DNA sequences

The Diplostomida Olson, Cribb, Tkach, Bray & Littlewood, 2003 is the less diverse order of the two orders within the subclass Digenea Carus, 1863 and is currently classified into three superfamilies, i.e. Brachylaimoidea Joyeux & Foley, 1930, Diplostomoidea Poirier, 1886, and Schistosomatoidea Stiles & Hassall, 1898. Although the suprageneric-level relationships have been elucidated with the use of molecular markers, the lack of representation of some groups obscure the phylogenetic relationships among families, rendering the classification unstable. Here, we tested the phylogenetic position of the family Proterodiplostomidae Dubois, 1936 based on partial 28S rDNA and complete 18S rDNA sequences for Crocodilicola pseudostoma (Willemoes-Suhm, 1870), a crocodile parasite that has been found as a progenetic metacercaria parasitising the pale catfish Rhamdia guatemalensis (Günther) in Mexico and in other siluruforms in the Neotropics. We augmented the representation of the species, genera and families within the Diplostomida, including mostly representatives of the superfamily Diplostomoidea, and assembled a dataset that contains 49 species for the 28S rRNA gene, and 45 species for the 18S rRNA gene. Additionally, we explored the phylogenetic signal of the mitochondrial gene cox1 in reconstructing the phylogenetic relationships of selected members of the superfamily. Our analyses showed that the family Proterodiplostomidae is the sister taxon to the paraphyletic Diplostomidae Poirier, 1886 and Strigeidae Railliet, 1919, with Cyathocotylidae Mühling, 1898 + Brauninidae Wolf, 1903 as their sister group. Analysis of concatenated 18S + 28S sequences revealed the Liolopidae Odhner, 1912 as the basal group of the superfamily Diplostomoidea, although analyses of independent datasets showed that the position of this family remains uncertain. Analysis based on cox1 unequivocally resolved the Proterodiplostomidae as the sister taxon to the Diplostomidae and Strigeidae, although the Cyathocotylidae was nested in a different clade, along with brachylaimoids and schistosomatoids.     

Phylogenetic Inferences and the Evolution of Plastid DNA in Campynemataceae and the Mycoheterotrophic <Emphasis Type="Italic">Corsia dispar</Emphasis> D.L Jones &amp; B. Gray (Corsiaceae)

The phylogenetic positions of the families Campynemataceae and Corsiaceae within the order Liliales remains unclear. To date, molecular data from the plastid genome of Corsiaceae has been obtained exclusively from Arachnitis, for which alignment and phylogenetic inference has proved difficult. The extent of gene conservation among mycoheterotrophic species within Corsiaceae remains unknown. To clarify the phylogenetic position of Campynemataceae and Corsiaceae within Liliales, functional plastid-coding genes of species representing both families have been analyzed. Examination of two phylogenetic data sets of plastid genes employing parsimony, maximum-likelihood, and Bayesian inference methods strongly supported both families forming a basal clade to the remaining taxa of Liliales. The first data set consists of five functional plastid-encoded genes (matK, rps7, rps2, rps19, and rpl2) sequenced from Corsia dispar (Corsiaceae). The data set included 31 species representing all families within Liliales, as well as selected orders that are related closely to Liliales (10 outgroup species from Asparagales, Dioscoreales, and Pandanales). The second phylogenetic analysis was based on 75 plastid genes. This data set included 18 species from Liliales, representing major clades within the order, and 10 outgroup species from Asparagales, Dioscoreales, and Pandanales. In this latter data set, Campynemataceae was represented by 60 plastid-encoded genes sequenced from herbarium material of Campynema lineare. A large proportion of the plastid genome of C. dispar was also sequenced and compared to the plastid genomes of photosynthetic plants within Liliales and mycoheterotrophic plants within Asparagales to explore plastid genome reduction. The plastid genome of C. dispar is in the advanced stages of reduction, which signifies its high dependency on mycorrhizal fungi and is suggestive of a loss in photosynthetic ability. Functional plastid genes found in C. dispar may be applicable to other species in Corsiaceae, which will provide a basis for in-depth molecular analyses of interspecies relationships within the family, once molecular data from other members become available.     

Molecular phylogenetic analysis of the ant genus <Emphasis Type="Italic">Formica</Emphasis> L. (Hymenoptera: Formicidae) from Palearctic region

Sixty-five sequences of the mitochondrial DNA cytochrome b gene fragment (759 bp) and 23 sequences of the NADH dehydrogenase subunit 6 gene fragment (224 bp) were compared in ants of the genus Formica L. from different regions of the Palearctic and in Polyergus rufescens Latr. as outgroup. In total, 28 species of the genus Formica were examined. As a result, dated trees with a molecular clock were constructed showing the phylogenetic relationships of Formica ants. The topology of the obtained tree based on the Cyt-b sequences was found to be not consistent with the generally accepted opinion on the Formica rufa and F. rufibarbis groups. New data on the formation history of the present-day fauna of Formica ants of the Palearctic were obtained. It was demonstrated that a considerable fraction of the examined species (about a third) were formed in the Quaternary Period.     

The complete mitochondrial genomes of <Emphasis Type="Italic">Tarsiger cyanurus</Emphasis> and <Emphasis Type="Italic">Phoenicurus auroreus</Emphasis>: a phylogenetic analysis of Passeriformes

Passeriformes is the largest group within aves and the phylogenetic relationships between Passeriformes have caused major disagreement in ornithology. Particularly, the phylogenetic relationships between muscicapoidea and sylvioidea are complex, and their taxonomic boundaries have not been clearly defined. Our aim was to study the status of two bird species: Tarsiger cyanurus and Phoenicurus auroreus. Furthermore, we analyzed the phylogenetic relationships of Passeriformes. Complete mitochondrial DNA (mtDNA) sequences of both species were determined and the lengths were 16,803 (T. cyanurus) and 16,772 bp (P. auroreus), respectively. Thirteen protein-coding genes, 22 tRNA genes, two rRNA genes, and one control region were identified in these mtDNAs. The contents of A and T at the base compositions was significantly higher than the content of G and C, and this AT skew was positive, while the GC skew was negative. The monophyly of Passeriformes is divided into four major clades: Corvoidea, Sylvioidea, Passeroidea, and Musicicapoidea. Paridae should be separated from the superfamily Sylvioidea and placed within the superfamily Muscicapoidea. The family Muscicapidae and Corvida were paraphyly, while Carduelis and Emberiza were grouped as a sister taxon. The relationships between some species of the order passeriformes may remain difficult to resolve despite an effort to collect additional characters for phylogenetic analysis. Current research of avian phylogeny should focus on adding characters and taxa and use both effectively to obtain a better resolution for deeper and shallow nodes.     

Evolution of glutamate dehydrogenase genes: evidence for lateral gene transfer within and between prokaryotes and eukaryotes

Background

Lateral gene transfer can introduce genes with novel functions into genomes or replace genes with functionally similar orthologs or paralogs. Here we present a study of the occurrence of the latter gene replacement phenomenon in the four gene families encoding different classes of glutamate dehydrogenase (GDH), to evaluate and compare the patterns and rates of lateral gene transfer (LGT) in prokaryotes and eukaryotes.

Results

We extend the taxon sampling of gdh genes with nine new eukaryotic sequences and examine the phylogenetic distribution pattern of the various GDH classes in combination with maximum likelihood phylogenetic analyses. The distribution pattern analyses indicate that LGT has played a significant role in the evolution of the four gdh gene families. Indeed, a number of gene transfer events are identified by phylogenetic analyses, including numerous prokaryotic intra-domain transfers, some prokaryotic inter-domain transfers and several inter-domain transfers between prokaryotes and microbial eukaryotes (protists).

Conclusion

LGT has apparently affected eukaryotes and prokaryotes to a similar extent within the gdh gene families. In the absence of indications that the evolution of the gdh gene families is radically different from other families, these results suggest that gene transfer might be an important evolutionary mechanism in microbial eukaryote genome evolution.

     

Genome-wide analysis and characterization of molecular evolution of the <Emphasis Type="Italic">HCT</Emphasis> gene family in pear (<Emphasis Type="Italic">Pyrus bretschneideri</Emphasis>)

Lignin is a major component of stone cells in pear fruit, which significantly affects fruit quality. Hydroxycinnamoyl CoA: shikimate hydroxycinnamoyl transferase (HCT), a recently discovered enzyme in plants, is an important gene that participates in the formation of lignin. Although HCT gene cloning and expression patterns have been studied in several species, including pear, there is still no extensive genome-wide bioinformatics analysis on the whole gene family, and the evolutionary history of HCT gene family is still unknown. A total of 82 HCT genes were identified in pear, most of which have one or two exons, and all with the conserved HXXXD motif and transferase domains. Based on the structural characteristics and phylogenetic analysis of these sequences, the HCT gene family genes could be classified into four main groups. Structural analysis also revealed that 25 % of HCT genes share a MYB binding site. Expansion of the HCT gene family mostly occurred before the divergence between Arabidopsis and Rosaceae, with whole-genome duplication or segmental duplication events playing the most important role in the expansion of the HCT gene family in pear. At the same time, purifying selection also played a critical role in the evolution of HCT genes. Five of the 82 HCT genes were verified by qRT-PCR to correspond to the pattern of stone cell formation during pear fruit development. The genome-wide identification, chromosome localization, gene structures, synteny, and expression analyses of pear HCT genes provide an overall insight into HCT gene family and their potential involvement in growth and development of stone cells.     

The threonine-tRNA ligase gene region is applicable in classification,typing, and phylogenetic analysis of bifidobacteria

In the modern era, molecular genetic techniques are crucial in ecological studies, as well as in the classification, typing, and phylogenetic analysis of prokaryotes. These techniques are mainly aimed at whole genome comparisons and PCR-derived experiments, including amplifying the 16S rRNA and other various housekeeping genes used in taxonomy, as well as MLST (multilocus sequence typing) and MLSA (multilocus sequence analysis) of different taxonomic bacterial groups. The gene encoding threonine-tRNA ligase (thrS) is a gene potentially applicable as an identification and phylogenetic marker in bacteria. It is widely distributed in bacterial genomes and is subject to evolutionary selection pressure due to its important function in protein synthesis. In this study, specific primers were used to amplify a thrS gene fragment (~740 bp) in 36 type and 30 wild strains classified under family Bifidobacteriaceae. The full-length gene has not yet been considered as a possible identification, classification, and phylogenetic marker in bifidobacteria. The thrS sequences revealed higher sequence variability (82.7% of pairwise identities) among members of the family than that shown by 16S rRNA gene sequences (96.0%). Although discrepancies were found between the thrS-derived and previously reported whole genome phylogenetic analyses, the main phylogenetic groups of bifidobacteria were properly assigned. Most wild strains of bifidobacteria were better differentiated based on their thrS sequences than on their 16S rRNA gene identities. Phylogenetic confidence of the evaluated gene with respect to other alternative genetic markers widely used in taxonomy of bifidobacteria (fusA, GroELhsp60, pyrG, and rplB genes) was confirmed using the localized incongruence difference - Templeton analysis.     

Phylogeny of anaerobic fungi (phylum <Emphasis Type="Italic">Neocallimastigomycota</Emphasis>), with contributions from yak in China

The phylum Neocallimastigomycota contains eight genera (about 20 species) of strictly anaerobic fungi. The evolutionary relationships of these genera are uncertain due to insufficient sequence data to infer their phylogenies. Based on morphology and molecular phylogeny, thirteen isolates obtained from yak faeces and rumen digesta in China were assigned to Neocallimastix frontalis (nine isolates), Orpinomyces joyonii (two isolates) and Caecomyces sp. (two isolates), respectively. The phylogenetic relationships of the eight genera were evaluated using complete ITS and partial LSU sequences, compared to the ITS1 region which has been widely used in this phylum in the past. Five monophyletic lineages corresponding to six of the eight genera were statistically supported. Isolates of Caecomyces and Cyllamyces were present in a single lineage and could not be separated properly. Members of Neocallimastigomycota with uniflagellate zoospores represented by Piromyces were polyphyletic. The Piromyces-like genus Oontomyces was consistently closely related to the traditional Anaeromyces, and separated the latter genus into two clades. The phylogenetic position of the Piromyces-like genus Buwchfawromyces remained unresolved. Orpinomyces and Neocallimastix, sharing polyflagellate zoospores, were supported as sister genera in the LSU phylogeny. Apparently ITS, specifically ITS1 alone, is not a good marker to resolve the generic affinities of the studied fungi. The LSU sequences are easier to align and appear to work well to resolve generic relationships. This study provides a comparative phylogenetic revision of Neocallimastigomycota isolates known from culture and sequence data.