The Complete Mitochondrial Genome of two Tetragnatha Spiders (Araneae: Tetragnathidae): Severe Truncation of tRNAs and Novel Gene Rearrangements in Araneae

Mitogenomes can provide information for phylogenetic analysis and evolutionary biology. The Araneae is one of the largest orders of Arachnida with great economic importance. In order to develop mitogenome data for this significant group, we determined the complete mitogenomes of two long jawed spiders Tetragnatha maxillosa and T. nitens and performed the comparative analysis with previously published spider mitogenomes. The circular mitogenomes are 14578 bp long with A+T content of 74.5% in T. maxillosa and 14639 bp long with A+T content of 74.3% in T. nitens, respectively. Both the mitogenomes contain a standard set of 37 genes and an A+T-rich region with the same gene orientation as the other spider mitogenomes, with the exception of the different gene order by the rearrangement of two tRNAs (trnW and trnG). Most of the tRNAs lose TΨC arm stems and have unpaired amino acid acceptor arms. As interesting features, both trnS^AGN and trnS^UCN lack the dihydrouracil (DHU) arm and long tandem repeat units are presented in the A+T-rich region of both the spider mitogenomes. The phylogenetic relationships of 23 spider mitogenomes based on the concatenated nucleotides sequences of 13 protein-coding genes indicated that the mitogenome sequences could be useful in resolving higher-level relationship of Araneae. The molecular information acquired from the results of this study should be very useful for future researches on mitogenomic evolution and genetic diversities in spiders.     

Complete Mitochondrial Genome of Echinostoma hortense (Digenea: Echinostomatidae)

Echinostoma hortense (Digenea: Echinostomatidae) is one of the intestinal flukes with medical importance in humans. However, the mitochondrial (mt) genome of this fluke has not been known yet. The present study has determined the complete mt genome sequences of E. hortense and assessed the phylogenetic relationships with other digenean species for which the complete mt genome sequences are available in GenBank using concatenated amino acid sequences inferred from 12 protein-coding genes. The mt genome of E. hortense contained 12 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and 1 non-coding region. The length of the mt genome of E. hortense was 14,994 bp, which was somewhat smaller than those of other trematode species. Phylogenetic analyses based on concatenated nucleotide sequence datasets for all 12 protein-coding genes using maximum parsimony (MP) method showed that E. hortense and Hypoderaeum conoideum gathered together, and they were closer to each other than to Fasciolidae and other echinostomatid trematodes. The availability of the complete mt genome sequences of E. hortense provides important genetic markers for diagnostics, population genetics, and evolutionary studies of digeneans.     

Mitochondrial DNA as effective molecular markers for the genetic variation and phylogeny of the family Osteoglossidae

The present study examined the genetic variation of the family Osteoglossidae from different geographical locations based on the mitochondrial NADH dehydrogenase subunit 2 (ND2) and ATPase subunit 6 (ATPase6) genes; we then re-constructed the phylogenetic relationships using the two sequences in combination. The results showed that the partial sequences of mitochondrial ND2 and ATPase6 of the family Osteoglossidae were 813 bp and 669 bp, respectively. A total of 42 species-specific nucleotide positions of the family Osteoglossidae were found to be useful for molecular identification. The sequence variation showed greater differences (8.3% ~ 28.1% for the combined sequences, 8.3% ~ 26.7% for the ND2 gene, and 9.3% ~ 28.7% for the ATPase6 gene) among the different species of Osteoglossidae, and there was a significant association between the genetic difference and geographical location. Phylogenetic analyses using neighbor-joining, Bayesian inference, and maximum parsimony (MP) methods based on the combined sequences of the two genes were able to distinguish the different species and were in agreement with the existing taxonomy based on morphological characters and in association with the geographical distribution among seven species of the family Osteoglossidae.     

Phylogenetic analyses provide the first insights into the evolution of OVATE family proteins in land plants

Background and Aims

The OVATE gene encodes a nuclear-localized regulatory protein belonging to a distinct family of plant-specific proteins known as the OVATE family proteins (OFPs). OVATE was first identified as a key regulator of fruit shape in tomato, with nonsense mutants displaying pear-shaped fruits. However, the role of OFPs in plant development has been poorly characterized.

Methods

Public databases were searched and a total of 265 putative OVATE protein sequences were identified from 13 sequenced plant genomes that represent the major evolutionary lineages of land plants. A phylogenetic analysis was conducted based on the alignment of the conserved OVATE domain from these 13 selected plant genomes. The expression patterns of tomato SlOFP genes were analysed via quantitative real-time PCR. The pattern of OVATE gene duplication resulting in the expansion of the gene family was determined in arabidopsis, rice and tomato.

Key Results

Genes for OFPs were found to be present in all the sampled land plant genomes, including the early-diverged lineages, mosses and lycophytes. Phylogenetic analysis based on the amino acid sequences of the conserved OVATE domain defined 11 sub-groups of OFPs in angiosperms. Different evolutionary mechanisms are proposed for OVATE family evolution, namely conserved evolution and divergent expansion. Characterization of the AtOFP family in arabidopsis, the OsOFP family in rice and the SlOFP family in tomato provided further details regarding the evolutionary framework and revealed a major contribution of tandem and segmental duplications towards expansion of the OVATE gene family.

Conclusions

This first genome-wide survey on OFPs provides new insights into the evolution of the OVATE protein family and establishes a solid base for future functional genomics studies on this important but poorly characterized regulatory protein family in plants.     

Orientobilharzia turkestanicum is a member of Schistosoma genus based on phylogenetic analysis using ribosomal DNA sequences

In the present study, samples representing Orientobilharzia turkestanicum from cattle, sheep, cashmere goat and goat in Heilongjiang Province, China, were characterized and grouped genetically by sequences of internal transcribed spacer (ITS, including ITS-1 and ITS2) and 28S ribosomal DNA (28S rDNA). The ITS and 28S rDNA were amplified by polymerase chain reaction (PCR) and then sequenced and compared with that of other members of the Schistosomatidae available in GenBank™, and phylogenetic relationships between them were re-constructed using the neighbor-joining and maximum parsimony methods. The lengths of ITS-1, ITS-2 and 28S rDNA sequences for all O. turkestanicum samples from different hosts were 384 bp, 331 bp and 1304 bp, respectively. While the ITS-1 sequences of O. turkestanicum from each of the four different hosts, and ITS-2 of O. turkestanicum from cattle, sheep and cashmere goat were identical, respectively, the ITS-2 of O. turkestanicum from goat differed from that of O. turkestanicum from cattle, sheep and cashmere goat by one nucleotide. The 28S rDNA sequences of O. turkestanicum from sheep and cashmere goat were identical, but differed from that of O. turkestanicum from cattle and goat by two nucleotides, with the latter two also having identical 28S rDNA sequence. Phylogenetic analyses based on the combined sequences of the ITS-1 and ITS-2, or the 28S rDNA sequences placed O. turkestanicum within the genus Schistosoma, and it was phylogenetically closer to the African schistosome group than to the Asian schistosome group. These results should have implications for studying the origin and evolution of O. turkestanicum and other members of the Schistosomatidae.     

Evolution of MHC class I genes in two ancient fish, paddlefish (Polyodon spathula) and Chinese sturgeon (Acipenser sinensis)

Here we present the first isolation of major histocompatibility complex (MHC) class I genes from two ancient fish, paddlefish (Polyodon spathula) and Chinese sturgeon (Acipenser sinensis). Seventeen sequences obtained showed high polymorphism and positive natural selection with d_N/d_S > 1. Evolutionary relationships revealed that sequences from paddlefish and Chinese sturgeon distinguished from other vertebrate class I and had an intermingling of alleles, which indicates that Acipenseriformes have a common ancestral gene of class I and a trans-species polymorphism across Acipenseriformes. We also found clear evidence of recombination among class I genes of paddlefish and Chinese sturgeon.     

Recruitment and diversification of an ecdysozoan family of neuropeptide hormones for black widow spider venom expression

Venoms have attracted enormous attention because of their potent physiological effects and dynamic evolution, including the convergent recruitment of homologous genes for venom expression. Here we provide novel evidence for the recruitment of genes from the Crustacean Hyperglycemic Hormone (CHH) and arthropod Ion Transport Peptide (ITP) superfamily for venom expression in black widow spiders. We characterized latrodectin peptides from venom gland cDNAs from the Western black widow spider (Latrodectus hesperus), the brown widow (Latrodectus geometricus) and cupboard spider (Steatoda grossa). Phylogenetic analyses of these sequences with homologs from other spider, scorpion and wasp venom cDNAs, as well as CHH/ITP neuropeptides, show latrodectins as derived members of the CHH/ITP superfamily. These analyses suggest that CHH/ITP homologs are more widespread in spider venoms, and were recruited for venom expression in two additional arthropod lineages. We also found that the latrodectin 2 gene and nearly all CHH/ITP genes include a phase 2 intron in the same position, supporting latrodectin's placement within the CHH/ITP superfamily. Evolutionary analyses of latrodectins suggest episodes of positive selection along some sequence lineages, and positive and purifying selection on specific codons, supporting its functional importance in widow venom. We consider how this improved understanding of latrodectin evolution informs functional hypotheses regarding its role in black widow venom as well as its potential convergent recruitment for venom expression across arthropods.     

Organellar genomes of the four-toothed moss,Tetraphis pellucida

Background

Mosses are the largest of the three extant clades of gametophyte-dominant land plants and remain poorly studied using comparative genomic methods. Major monophyletic moss lineages are characterised by different types of a spore dehiscence apparatus called the peristome, and the most important unsolved problem in higher-level moss systematics is the branching order of these peristomate clades. Organellar genome sequencing offers the potential to resolve this issue through the provision of both genomic structural characters and a greatly increased quantity of nucleotide substitution characters, as well as to elucidate organellar evolution in mosses. We publish and describe the chloroplast and mitochondrial genomes of Tetraphis pellucida, representative of the most phylogenetically intractable and morphologically isolated peristomate lineage.

Results

Assembly of reads from Illumina SBS and Pacific Biosciences RS sequencing reveals that the Tetraphis chloroplast genome comprises 127,489 bp and the mitochondrial genome 107,730 bp. Although genomic structures are similar to those of the small number of other known moss organellar genomes, the chloroplast lacks the petN gene (in common with Tortula ruralis) and the mitochondrion has only a non-functional pseudogenised remnant of nad7 (uniquely amongst known moss chondromes).

Conclusions

Structural genomic features exist with the potential to be informative for phylogenetic relationships amongst the peristomate moss lineages, and thus organellar genome sequences are urgently required for exemplars from other clades. The unique genomic and morphological features of Tetraphis confirm its importance for resolving one of the major questions in land plant phylogeny and for understanding the evolution of the peristome, a likely key innovation underlying the diversity of mosses. The functional loss of nad7 from the chondrome is now shown to have occurred independently in all three bryophyte clades as well as in the early-diverging tracheophyte Huperzia squarrosa.     

Evolution of the tRNA gene family in mitochondrial genomes of five Meretrix clams (Bivalvia,Veneridae)

In contrast to the extreme conservation of nuclear-encoded tRNAs, organization of the mitochondrial (mt) tRNA gene family in invertebrates is highly dynamic and rapidly evolving. While gene duplication and loss, gene isomerism, recruitment, and rearrangements have occurred sporadically in several invertebrate lineages, little is known regarding the pattern of their evolution. Comparisons of invertebrate mt genomes at a generic level can be extremely helpful in investigating evolutionary patterns of variation, as intermediate stages of the process may be identified. Variation of mitochondrial tRNA organization among Meretrix clams provides good materials to investigate mt tRNA evolution. We characterized the complete mt genome of the lyrate Asiatic hard clam Meretrix lyrata, re-annotated tRNAs of four previously sequenced Meretrix clams, and undertook an intensive comparison of tRNA gene families in these clams. Our results 1) provide evidence that the commonly observed duplication of trnM may have occurred independently in different bivalve lineages and, based on the higher degree of trnM gene similarity, may have occurred more recently than expected; 2) suggest that “horizontal” evolution may have played an important role in tRNA gene family evolution based on frequent gene duplications and gene recruitment events; and 3) reveal the first case of isoacceptor “vertical” tRNA gene recruitment (VTGR) and present the first clear evidence that VTGR allows rapid evolution of tRNAs. We identify the trnS^− UCR gene in Meretrix clams, previously considered missing in this lineage, and speculate that trnS^− UCR lacking the D-arm in both M. lyrata and Meretrix lamarckii may represent the ancestral status. Phylogenetic analysis based on 13 concatenate protein-coding genes provided opportunities to detect rapidly evolved tRNA genes via VTGR and gene isomerism processes. This study suggests that evolution of the mt tRNA gene family in bivalves is more complex than previously thought and that comparison of several congeneric species is a useful strategy in investigating evolutionary patterns and dynamics of mt tRNA genes.     

Fitness Assays Reveal Incomplete Functional Redundancy of the HoxA1 and HoxB1 Paralogs of Mice

Gene targeting techniques have led to the phenotypic characterization of numerous genes; however, many genes show minimal to no phenotypic consequences when disrupted, despite many having highly conserved sequences. The standard explanation for these findings is functional redundancy. A competing hypothesis is that these genes have important ecological functions in natural environments that are not needed under laboratory settings. Here we discriminate between these hypotheses by competing mice (Mus musculus) whose Hoxb1 gene has been replaced by Hoxa1, its highly conserved paralog, against matched wild-type controls in seminatural enclosures. This Hoxb1^A1 swap was reported as a genetic manipulation resulting in no discernible embryonic or physiological phenotype under standard laboratory tests. We observed a transient decline in first litter size for Hoxb1^A1 homozygous mice in breeding cages, but their fitness was consistently and more dramatically reduced when competing against controls within seminatural populations. Specifically, males homozygous for the Hoxb1^A1 swap acquired 10.6% fewer territories and the frequency of the Hoxb1^A1 allele decreased from 0.500 in population founders to 0.419 in their offspring. The decrease in Hoxb1^A1 frequency corresponded with a deficiency of both Hoxb1^A1 homozygous and heterozygous offspring. These data suggest that Hoxb1 and Hoxa1 are more phenotypically divergent than previously reported and support that sub- and/or neofunctionalization has occurred in these paralogous genes leading to a divergence of gene function and incomplete redundancy. Furthermore, this study highlights the importance of obtaining fitness measures of mutants in ecologically relevant conditions to better understand gene function and evolution.     

Proposal to reclassify the Streptomyces albidoflavus clade on the basis of multilocus sequence analysis and DNA–DNA hybridization,and taxonomic elucidation of Streptomyces griseus subsp. solvifaciens

The Streptomyces albidoflavus 16S rRNA gene clade contains 10 species and subspecies with identical 16S rRNA gene sequences and very similar numerical taxonomic data, including Streptomyces griseus subsp. solvifaciens. Type strains of this clade, as well as three CGMCC strains which were received as Streptomyces galilaeus, Streptomyces sioyaensis and Streptomyces vinaceus, respectively, that shared the same 16S rRNA gene sequences with the clade, were subjected to multilocus sequence analysis (MLSA), DNA–DNA hybridization (DDH) and phenotypic characterization for a comprehensive reevaluation. The 13 strains still formed a distinct, albeit loosely related, clade in the phylogenetic tree based on concatenated sequences of aptD, gyrB, recA, rpoB and trpB genes, supported by a high bootstrap value and different tree-making algorithms, with MLSA evolutionary distances ranging from 0 to 0.003. DDH values among these strains were well above the 70% cut-off point for species delineation. Based on the genotypic data of MLSA and DDH, combined with key phenotypic properties in common, it is proposed that the 10 species and subspecies of the S. albidoflavus clade, namely S. albidoflavus, S. canescens, S. champavatii, S. coelicolor, S. felleus, S. globisporus subsp. caucasicus, S. griseus subsp. solvifaciens, S. limosus, S. odorifer and S. sampsonii, should be merged into a single genomic species, for which the name S. albidoflavus is retained, and that the three strains S. galilaeus CGMCC 4.1320, S. sioyaensis CGMCC 4.1306 and S. vinaceus CGMCC 4.1305 should be assigned to S. albidoflavus as well. The results also indicated that MLSA could be the procedure of choice for distinguishing between species within Streptomyces 16S rRNA gene clades.     

Complete mitochondrial genome sequences of three bats species and whole genome mitochondrial analyses reveal patterns of codon bias and lend support to a basal split in Chiroptera

Order Chiroptera is a unique group of mammals whose members have attained self-powered flight as their main mode of locomotion. Much speculation persists regarding bat evolution; however, lack of sufficient molecular data hampers evolutionary and conservation studies. Of ~ 1200 species, complete mitochondrial genome sequences are available for only eleven. Additional sequences should be generated if we are to resolve many questions concerning these fascinating mammals. Herein, we describe the complete mitochondrial genomes of three bats: Corynorhinus rafinesquii, Lasiurus borealis and Artibeus lituratus. We also compare the currently available mitochondrial genomes and analyze codon usage in Chiroptera. C. rafinesquii, L. borealis and A. lituratus mitochondrial genomes are 16438 bp, 17048 bp and 16709 bp, respectively. Genome organization and gene arrangements are similar to other bats. Phylogenetic analyses using complete mitochondrial genome sequences support previously established phylogenetic relationships and suggest utility in future studies focusing on the evolutionary aspects of these species. Comprehensive analyses of available bat mitochondrial genomes reveal distinct nucleotide patterns and synonymous codon preferences corresponding to different chiropteran families. These patterns suggest that mutational and selection forces are acting to different extents within Chiroptera and shape their mitochondrial genomes.     

Phylogenetic Relationships Among Selected Heteroderoidea Based on 18S and ITS Ribosomal DNA

In a study of relationships among selected cyst-forming and noncyst-forming species of Heteroderoidea, combined sequences comprised of DNA from part of the conserved 18S ribosomal RNA gene (rDNA) plus the complete ITS rDNA segment were more similar to analyses based on the ITS data alone than to analyses based on the 18S data alone. One of the two noncyst-forming species, Ekphymatodera thomasoni, grouped with cyst-forming species of Heteroderoidea. Bilobodera flexa, also a noncyst-forming species, was separated from all the other taxa by a long branch. Afenestrata koreana, with a weakly sclerotized cyst, grouped closely with H. bifenestra. These observations suggest that phylogenetic analyses using molecular data may aid in our understanding of the evolution of cyst formation in nematodes, including the possibility of secondary loss. The usefulness of molecular phylogenetic analyses in nematodes may depend more on the particular selection of taxa than on mere addition of data from additional genes.     

Molecular evolution of pteridophytes and their relationship to seed plants: Evidence from complete 18S rRNA gene sequences

Complete 18S ribosomal RNA sequence data from representatives of all extant pteridophyte lineages together with RNA sequences from different seed plants were used to infer a molecular phylogeny of vascular plants that included all major land plant lineages. The molecular data indicate that lycopsids are monophyletic and are the earliest diverging group within the vascular land plants, whereasPsilotum nudum is more closely related to the seed plants than to other pteridophyte lineages. The phylogenetic trees based on maximum likelihood, parsimony and distance analyses show substantial agreement with the evolutionary relationships of land plants as interpreted from the fossil record.     

Comparative analyses of the complete mitochondrial genomes of the two ruminant hookworms Bunostomum trigonocephalum and Bunostomum phlebotomum

Bunostomum trigonocephalum and Bunostomum phlebotomum are blood-feeding hookworms of sheep and cattle, causing considerable economic losses to the live stock industries. Studying genetic variability within and among hookworm populations is critical to addressing epidemiological and ecological questions. Mitochondrial (mt) DNA is known to provide useful markers for investigations of population genetics of hookworms, but mt genome sequence data are scant. In the present study, the complete mitochondrial DNA (mtDNA) sequences of the sheep and goat hookworm B. trigonocephalum were determined for the first time, and the mt genome of B. phlebotomum from yak in China was also sequenced for comparative analyses of their gene contents and genome organizations. The lengths of mt DNA sequences of B. trigonocephalum sheep isolate, B.trigonocephalum goat isolate and B. phlebotomum China yak isolate were 13,764 bp, 13,771 bp and 13,803 bp in size, respectively. The identity of the mt genomes was 99.7% between B. trigonocephalum sheep isolate and B. trigonocephalum goat isolate. The identity of B. phlebotomum China yak isolate mt genomes was 85.3% with B. trigonocephalum sheep isolate, and 85.2% with B. trigonocephalum goat isolate. All the mt genes of the two hookworms were transcribed in the same direction and gene arrangements were consistent with those of the GA3 type, including 12 protein-coding genes, 2 rRNA genes and 22 tRNA genes, but lacking ATP synthetase subunit 8 gene. The mt genomes of B. trigonocephalum and B. phlebotomum were similar to prefer bases A and T, the contents of A + T are 76.5% (sheep isolate), 76.4% (goat isolate) and 76.9% (China yak isolate), respectively. Phylogenetic relationships reconstructed using concatenated amino acid sequences of 12 protein-coding genes with three methods (maximum likelihood, Bayesian inference and neighbor joining) revealed that the B. trigonocephalum and B. phlebotomum represent distinct but closely-related species. These data provide novel and useful genetic markers for studying the systematics, and population genetics of the two ruminant hookworms.     

A conserved extraordinarily long serine homopolymer in Dictyostelid amoebae

Eukaryotic protein sequences often contain amino-acid homopolymers that consist of a single amino acid repeated from several to dozens of times. Some of these are functional but others may persist largely because of high expansion rates due to DNA slippage. However, very long homopolymers with over a hundred repeats are very rare. We report an extraordinarily long homopolymer consisting of 306 tandem serine repeats from the single-celled eukaryote Dictyostelium discoideum, which also has a multicellular stage. The gene has a paralog with 132 repeats and orthologs, also with high serine repeat numbers, in various other Dictyostelid species. The conserved gene structure and protein sequences suggest that the homopolymer is functional. The high codon diversity and very poor alignment of serine codons in this gene between species similarly indicate functionality. This is because the serine homopolymer is conserved despite much DNA sequence change. A survey of other very long amino-acid homopolymers in eukaryotes shows that high codon diversity is the rule, suggesting that these too may be functional.     

Multilocus sequence analysis of putative Vibrio mediterranei strains and description of Vibrio thalassae sp. nov.

A multilocus sequence analysis based on partial gyrB, mreB, rpoD and pyrH genes was undertaken with 61 putative Vibrio mediterranei/V. shilonii strains from different hosts (mussels, oysters, clams, coral, fish and plankton) or habitat (seawater and sediment) and geographical origins (Mediterranean, Atlantic and Pacific). A consistent grouping was obtained with individual and concatenated gene sequences, and the clade, comprising 54 strains, was split into three subclades by all methods: subclade A (40 strains, including AK1, the former type strain of Vibrio shilonii), subclade B (8 strains) corresponding to the species V. mediterranei, and subclade C (six strains) representing a new species, V. thalassae sp. nov., with strain MD16^T (=CECT 8203^T = KCTC 32373^T) as the proposed type strain.     

Genome Comparison and Phylogenetic Analysis of Orientia tsutsugamushi Strains

Orientia tsutsugamushi (OT) is an obligate intracellular bacterium belonging to the family Rickettsiaceae and is the causative agent of scrub typhus, or Tsutsugamushi disease. The complete genome sequences of two OT strains (Boryong and Ikeda) have recently been determined. In the present study, we performed a fine genome sequence comparison of these strains. Our results indicate that although the core gene set of the family Rickettsiaceae is highly conserved between the two strains, a common set of repetitive sequences have been explosively amplified in both genomes. These amplified repetitive sequences have induced extensive genome shuffling and duplications and deletions of many genes. On the basis of the results of the genome sequence comparison, we selected 11 housekeeping genes and carried out multilocus sequence analysis of OT strains using the nucleotide sequences of these genes. This analysis revealed for the first time the phylogenetic relationships of representative OT strains. Furthermore, the results suggest the presence of an OT lineage with higher potential for virulence, which may explain the clinical and epidemiological differences between ‘classic’ and ‘new’ types of Tsutsugamushi disease in Japan.     

Evaluation of the use of multilocus sequence analysis (MLSA) to resolve taxonomic conflicts within the genus Marichromatium

Four species of marine purple sulfur bacteria of the genus Marichromatium have been validly described. A recent re-analysis of the 16S rRNA-based similarity and genomic DNA–DNA hybridizations (DDH) of the type strains [33] suggested that some of them are so closely related that they can be considered heterotypic synonyms. Here, we report on the evaluation of the multilocus sequence analysis approach (MLSA) for nine Marichromatium strains in order to resolve their intrageneric genealogical relationships. MLSA was based on six protein-coding genes (gyrB, recA, fusA, dnaK, pufM, and soxB), and the results were comparable to DDH. The phylogenetic tree constructed with the concatenated sequences, which also included the 16S rRNA gene and the internal transcriber spacer ITS region (4331 bp), separated the nine strains in four lineages that reflected the four Marichromatium species. The reconstructed phylogenetic tree based on concatenation of six protein-coding genes was also highly congruent with the tree topology based on the 16S rRNA gene.