Open Access

Polyphyletic origin of the genus Physarum(Physarales, Myxomycetes) revealed by nuclear rDNA mini-chromosome analysis and group I intron synapomorphy

  • Satish CR Nandipati1,
  • Kari Haugli1,
  • Dag H Coucheron1,
  • Edward F Haskins2 and
  • Steinar D Johansen1Email author
BMC Evolutionary Biology201212:166

DOI: 10.1186/1471-2148-12-166

Received: 15 February 2012

Accepted: 15 August 2012

Published: 31 August 2012

Abstract

Background

Physarales represents the largest taxonomic order among the plasmodial slime molds (myxomycetes). Physarales is of particular interest since the two best-studied myxomycete species, Physarum polycephalum and Didymium iridis, belong to this order and are currently subjected to whole genome and transcriptome analyses. Here we report molecular phylogeny based on ribosomal DNA (rDNA) sequences that includes 57 Physarales isolates.

Results

The Physarales nuclear rDNA sequences were found to be loaded with 222 autocatalytic group I introns, which may complicate correct alignments and subsequent phylogenetic tree constructions. Phylogenetic analysis of rDNA sequences depleted of introns confirmed monophyly of the Physarales families Didymiaceae and Physaraceae. Whereas good correlation was noted between phylogeny and taxonomy among the Didymiaceae isolates, significant deviations were seen in Physaraceae. The largest genus, Physarum, was found to be polyphyletic consisting of at least three well supported clades. A synapomorphy, located at the highly conserved G-binding site of L2449 group I intron ribozymes further supported the Physarum clades.

Conclusions

Our results provide molecular relationship of Physarales genera, species, and isolates. This information is important in further interpretations of comparative genomics nd transcriptomics. In addition, the result supports a polyphyletic origin of the genus Physarum and calls for a reevaluation of current taxonomy.

Background

Myxomycetes (plasmodial slime molds) are eukaryotic microorganisms that according to Olive [1] represent one of four main groups of slime molds (Mycetozoans). The myxomycetes consist of about 850 assigned species classified into five orders (Physarales, Stemonitales, Trichiales, Liceales and Echinosteliales) [2]. A typical myxomycete species has a complex sexual life cycle that consists of two vegetative stages; a haploid unicellular stage (amoeba/ flagellate) and a diploid syncytium stage (plasmodium), as well as several dormant and developmental stages [35]. Myxomycetes are commonly found in nature on decaying plant materials where they feed on a variety of bacteria and unicellular eukaryotes as well as dissolved plant nutrients [58]. Identification of species is mainly based on morphological characters including fruiting body structures and sporocarp colors, and lime deposition [6, 9]. More recently, nuclear DNA sequence markers have contributed in resolving relationships among and within taxonomic groups. The ribosomal DNA (rDNA) spacers [10], obligatory group I introns [11], ribosomal RNA (rRNA) genes [1116], and the elongation factor-1α (EF-1α) gene [17] and have all been applied.

The order Physarales is of special interest since it contains the two best-studied myxomycete species at the biochemical and molecular levels (Physarum polycephalum and Didymium iridis) [6, 18]. One unusual feature among the Physarales is the linear multicopy nature of the nuclear rDNA minichromosome, which ranges in size from about 20 kb to 80 kb among species [1922]. Whereas the rDNA minichromosome in D. iridis is only 21 kb and contains a single pre-rRNA transcription unit, a 60-kb palindromic rDNA with two transcription units has been characterized in P. polycephalum. A hallmark of Physarales rDNA minichromosomes is the presence of multiple group I introns within the rRNA coding regions. All isolates investigated to date contain at least two group I introns in the large subunit (LSU) rRNA gene, and high intron loads are exemplified in isolates of Fuligo septica and Diderma niveum which contain 12 and 21 intron insertions, respectively [11, 2225], our unpublished results. Group I introns code for ribozymes (catalytic RNAs) that perform self-splicing by a common molecular mechanism based on a series of transesterification reactions [25]. A group I ribozyme is organized into a well-defined and highly conserved RNA core structure that consists of three helical stacks referred to as the catalytic domain (P3, P7-P9), the substrate domain (P1, P2, P10), and the scaffold domain (P4-P6) [26].

Here we report the presence of 222 rDNA group I introns and molecular phylogeny of 57 Physarales isolates (2 families, 10 genera, 31 defined species) based on combined data sets of nuclear small subunit (SSU) and LSU rRNA gene sequences. The analysis supports a polyphyletic origin of the Physarum genus. Phylogenetic analysis of Physarales is of particular importance since two representative species are currently under whole genome and transcriptome analyses. Whereas P. polycephalum is sequenced at The Genome Institute – Washington University, our laboratory at University of Tromsø investigates D. iridis by deep sequencing technologies. Resolving relationships of genera or species within and between the two Physarales families (Physaraceae and Didymiaceae) are crucial in interpretations of comparative genome and transcriptome data.

Results and discussion

High load of group I introns in myxomycete rDNA

A characteristic feature of myxomycete rRNA gene sequences is the high load of autocatalytic group I introns [11, 12, 24, 25, 27, 28]. These introns have to be correctly identified by structural characterizations and removed from coding sequences prior to phylogenetic analysis. Myxomycete group I introns have been found at 23 different, but highly conserved, insertion sites [25], and several of the introns are diverse in sequence and highly complex in organization due to internal protein coding genes, large arrays of direct repeat motifs, and even additional catalytic RNA domains [11, 29].

The rDNA sequences included in this study (near complete SSU rRNA gene and partial LSU rRNA gene; Figure 1) were examined for the presence of group I introns. We observed a massive intron contents (231 group I introns in 81 isolates; Additional file 1: Table S1, Additional file 2: Table S2, and Additional file 3: Table S3). Myxomycetes belonging to the order Physarales (Physaraceae and Didymiaceae families) harbored the majority of intron insertions (222 introns). Here, 55 and 167 introns were recognized in 21 and 36 isolates of Physaraceae and Didymiaceae, respectively (Additional file 3: Table S2 and Additional file 4: Table S3). Interestingly, all Physarales isolates investigated contained two obligatory group I introns at positions L1949 and L2449 in the LSU rRNA gene. Introns at these sites have been found to be strictly vertically inherited [11], have probably gained an essential biological role in benefit for the host cell [25], and possess highly complex structure and organization features at the RNA and DNA levels (see below).
https://static-content.springer.com/image/art%3A10.1186%2F1471-2148-12-166/MediaObjects/12862_2012_Article_2143_Fig1_HTML.jpg
Figure 1

Schematic presentation of the analyzed rDNA regions in myxomycetes. The SSU rRNA and LSU rRNA gene regions used in phylogenetic reconstructions correspond to approximately 1800 nt and 750 nt, respectively. Myxomycete SSU and LSU rRNA genes are frequently interrupted by group I introns in myxomycetes. The position of each intron insertion site in analyzed regions is indicated and numbered according to the E.coli rRNA gene [25].

Relationships among the mycetozoan rDNA

According to Olive [1] the mycetozoans have been divided into four main groups of slime molds; the myxomycetes, dictyostelids, protostelids and acrasids. As an initial approach to investigate the relationships among the mycetozoans groups we performed molecular phylogeny based on unambiguously aligned SSU rRNA genes sequences (1028 nt positions) from 26 defined species representing the four mycetozoan groups as well as various amoebozoans (Table 1). These sequences were either retrieved directly from the sequence database (http://www.ncbi.nlm.nih.gov) or generated in our laboratory by PCR amplification using specific primers followed by Sanger sequence determination (see Methods).
Table 1

Mycetozoan and representative amoebozoan isolates

Species

Isolate(a)

SSU(b)

Accession no(c)

MYCETOZOA (d)

Myxomycetes (plasmodial slime molds)

Physarales (Order)

Didymium iridis

Pan2

+

AJ938153

Physarum polycephalum

Wis1

+

X13160

Stemonitales (Order)

Comatricha nigricapillitia

AMFD114

+

AY643824

Stemonites flavogenita

ATCC24714

+

HE655085

Trichiales (Order)

Arcyria stipata

AMFD257

+

EF513170

Trichia persimilis

--

+

AY643826

Liceales (Order)

Cribraria cancellata

AMFD94

+

EF513177

Echinosteliales (Order)

Echinostelium minutum

ATCC22345

+

HE655087

Protostelids

Soliformovum irregulare

ATCC26826

+

HE655088

Dictyostelids (cellular slime molds)

Acytostelium ellipticum

ATCC22247

+

HE655086

Acytostelium leptosomum

FG12

+

AM168111

Acytostelium subglobosum

LB1

+

AM168110

Dictyostelium discoideum

--

+

K02641

Dictyostelium fasciculatum

SH3

+

AM168087

Dictyostelium medusoides

OH592

+

AM168088

Dictyostelium rhizopodium

AusKY-4

+

AM168063

Acrasids

Acrasis rosea

T-235

+

AF011458

AMOEBOZOA

Amoeba leningradensis

CCAP1503/6

+

AJ314605

Acanthamoeba palestinensis

CCAP1547/1

+

L09599

Entamoeba histolytica

HM1-IMSS

+

X65163

Filamoeba nolandi

ATCC50430

+

AF293896

Gephyramoeba sp.

ATCC50654

+

AF293897

Hartmannella abertawensis

Page180

+

DQ190241

Mastigella commutans

--

+

AF421219

Naegleria gruberi

NEG-M

+

AB298288

Platyamoeba placida

--

+

AY294150

(a) Source of the organism (see Additional file 8: Table S4).

(b) SSU rDNA sequence.

(c) GenBank/EMBL/DDJB accession numbers.

(d) Species are grouped into phyla of Protist.

--, no isolate name given; +, analysed.

Some significant findings are noted from the un-rooted maximum likelihood (ML) tree (Additional file 1: Figure S1). The dictyostelids and myxomycetes were both found to represent monophyletic groups supported by high bootstrap values for ML, maximum parsimony (MP), neighbour joining (NJ), and Bayesian inference (BAY) analyses. Furthermore, each of the five main myxomycete orders (Physarales, Stemonitales, Trichiales, Liceales and Echinosteliales) represented distinct groups within the plasmodial slime molds. These findings corroborate earlier studies on the phylogeny of slime molds [1316, 30, 31]. Finally, Soliformovum irregulare (protostelid) and Acrasis rosea (acrasid) appeared both distantly related from each other and from the dictyostelids and myxomycetes. From this analysis we conclude that the apparent phenotypic similarity between the mycetozoan groups [1] is not reflected in the genetic relationship.

Relationships within the Didymiaceae

Data presented in Additional file 1: Figure S1 and previous studies [1316] supported a monophyletic origin of the order Physarales among the myxomycetes. To gain deeper insights in the relationships among families, genera, species and isolates within the order Physarales we performed molecular phylogeny based on two different data sets of nuclear rRNA gene sequences (Figure 1). The first data set represents an alignment of 2522 nt positions consisting of a near complete SSU rRNA gene in combination with a segment of the LSU rRNA gene (SSU/LSU data set; Additional file 5: Figure S2). The second data set consists of the LSU rRNA gene segment only (765 bp LSU data set; Additional file 6: Figure S3 and Additional file 7: Figure S4).

A representative ML tree based on the SSU/LSU data set is presented in Figure 2. The ML tree includes 17 isolates (4 genera, 8 defined species) of the Didymiaceae family, and supported by NJ/MP/ML > 64% and BAY posterior probability > 0.85 (Table 2). The monophyly of the four genera (Didymium, Lepidoderma, Mucilago and Diderma) was supported by NJ/MP/ML > 93% and BAY posterior probability of 1.0 (Figure 2). Furthermore, phylogenetic analysis was performed on the LSU data set (Table 2) generated from 36 Didymiaceae isolates (4 genera, 17 defined species). A representative ML tree is presented in Figure 3A. This analysis is in general agreement with that of the SSU/LSU data set, but with less statistical support at basal nodes of the Didymium and Diderma genera. However, the Lepidoderma genus was strongly supported in both data sets (Figure 2 and Figure 3A). These findings corroborate more recent studies of Didymiaceae relationships [11, 12], and we conclude that rDNA-based phylogeny is well consistent with the current Didymiaceae taxonomy [2, 32].
https://static-content.springer.com/image/art%3A10.1186%2F1471-2148-12-166/MediaObjects/12862_2012_Article_2143_Fig2_HTML.jpg
Figure 2

ML phylogenetic tree of Physarales isolates based on the combined SSU/LSU data set. The analysis was based on the alignment of 2522 nucleotide positions from 33 isolates representing the Didymiaceae and Physaraceae families as well as Stemonitis flavogenita (order Stemonitales) as out-group. The evolution model is GTR + I + G. Tree reconstruction using NJ, MP, ML, or BAY resulted in essential identical tree topologies. Bootstrap values (> 50%) of 2000 replications are shown above the branches (NJ/MP/ML). Node of high statistical support (bootstrap values of 100%) is indicated by red circle. BAY analysis consists of 1,000,000 generations with trees sampled every 100 generations. BAY posterior probabilities are shown below branches, and branches with posterior probabilities of 1.0 are presented as thick black lines. The Physarales clades (Clades I to III) are indicated.

Table 2

Didymiaceae isolates and rDNA sequences

Species

Isolate

SSU(a)

Accession no(b)

LSU(a)

Accession no(b)

Didymium

D. clavus

It-IG45

na

 

+

HE655047

D. dubium

Fr-K7

+

HE614606

+

AM407422

D. dubium

Fr-K15

+

HE614607

+

AM407423

D. dubium

It-K64

na

 

+

HE655048

D. dubium

Uk-K80

na

 

+

HE655049

D. dubium

Uk-K77

na

 

+

HE655050

D. iridis

Pan1-66

na

 

+

AM407418

D. iridis

Pan2

+

AJ938153

+

AM407414

D. iridis

Pan3-3

na

 

+

AM407419

D. iridis

Hon1-7

+

AJ938152

+

HE655051

D. iridis

CUR1-4

+

AJ938150

+

AM407420

D. iridis

HA4-1

+

AJ938149

+

AM407416

D. iridis

CR19-1

+

AJ938151

+

AM407421

D. iridis

CR8-1

+

AJ938154

+

AM407415

D. squamulosum

Cr10

+

HE614613

+

AM407427

Diderma

D. meyerae

It-K61

+

HE614614

+

HE655059

D. microcarpum

Uk-K93

na

 

+

HE655052

D. niveum

Fr-K10

+

HE614615

+

AM407429

D. niveum

Fr-M26

na

+

 

AM407425

D. niveum

It-K66

+

HE614616

+

HE655060

D. niveum

Uk-K79

+

HE614617

+

HE655061

D. saondersii

Mx-K30

na

 

+

AM407428

D. testaceum

It-IG50

na

 

+

HE655053

Diderma sp.

It-K68

na

 

+

HE655054

Diderma sp.

Fr-K12

na

 

+

AM407426

Diderma sp.

It-K56

na

 

+

HE655057

Diderma sp.

Uk-K78

na

 

+

HE655058

Diderma sp.

It-IG46

na

 

+

HE655055

Diderma sp.(c)

Pr3-1

+

HE614612

+

HE655056

Lepidoderma

L. aggregatum

Uk-K86

na

 

+

HE655062

L. carestianum

Fr-K18

+

HE614609

+

AM407430

L. carestianum

It-K71

+

HE614618

+

HE655063

L. crustaceum

It-K62

+

HE614619

+

HE655064

L. peyerimhoffii

It-K63

na

 

+

HE655065

Lepidoderma sp.

It-K52

na

 

+

HE655066

Mucilago

M. crustacea

No-K94

+

HE614620

+

HE655067

(a) SSU and LSU rDNA sequence.

(b) GenBank/EMBL/DDJB accession numbers.

(c) Reported previous as Didymium anellus in [11, 12].

+, analysed; na, not analysed.

https://static-content.springer.com/image/art%3A10.1186%2F1471-2148-12-166/MediaObjects/12862_2012_Article_2143_Fig3_HTML.jpg
Figure 3

ML phylogenetic tree of Physarales isolates based on the LSU data set. The analysis was based on the alignment of 765 nucleotide positions from 57 isolates representing the Didymiaceae and Physaraceae families. The evolution models are TrNef + I (Didymiaceae) and GTR + I + G (Physaraceae). Tree reconstruction using NJ, MP, ML, or BAY resulted in essential similar tree topologies. Bootstrap values (> 50%) of 2000 replications are shown above the branches (NJ/MP/ML). Node of high statistical support (bootstrap values of 100%) is indicated by red circle. BAY analysis consists of 1,000,000 generations with trees sampled every 100 generations. BAY posterior probabilities are shown below branches, and branches with posterior probabilities of 1.0 are presented as thick black lines. A. ML-tree of Didymiaceae isolates. P. polycephalum was used as an out-group. B. ML-tree of Physaraceae isolates. D. iridis was used as an out-group.

Relationships within Physaraceae

Figure 2 presents molecular phylogeny of 15 isolates (4 genera, 13 defined species) of the family Physaraceae based on the SSU/LSU data set. These isolates were separated from those of Didymiaceae with strong statistical supports represented by NJ/MP/ML bootstraps > 96% and BAY posterior probabilities of 1.0 (Figure 2). Furthermore, an extended analysis including 21 Physaraceae isolates (6 genera; 16 defined species; Table 3) based on the LSU data set is presented in Figure 3B. A different clustering pattern was seen among the Physaraceae genera compared to that of Didymiaceae. The ten investigated Physarum species were found to intersperse among isolates belonging to genera Badhamia, Craterium, Fuligo, Leocarpus and Physarella (Figures 2 and 3B). Three main clades (Clade I to III) with strong statistical support (SSU/LSU data set: NJ/MP/ML bootstraps > 84%, BAY posterior probabilities of 1; LSU data set: NJ/MP/ML bootstraps > 96%, BAY posterior probabilities > 0.94) were found to harbor both Physarum and Badhamia isolates (Figure 3B).
Table 3

Physaraceae isolates and rDNA sequences

Species

Isolate

SSU(a)

Accession no(b)

LSU(a)

Accession no(b)

Badhamia

B. melanospora

Az4-1

+

HE614610

+

HE655068

B. melanospora

Pr-1

+

HE614596

+

HE655069

B. utricularis

--

+

HE614597

+

HE655070

Physaracea-like sp.

Cur1

+

HE614608

+

HE663133

Craterium

C. minutum

It-IG38

na

 

+

HE655081

Fuligo

F. septica

IW-1

na

 

+

HE655082

F. septica

Mx-K28

na

 

+

HE655083

F. septica

NY-1

+

AJ584697

+

AJ584697

Leocarpus

L. fragilis

It-IG39

na

 

+

HE655071

Physarella

P. oblonga

--

+

HE614598

+

HE655072

Physarum

P. albescens

Fr-K2

na

 

+

HE655073

P. bivalve

It-IG42

na

 

+

HE655084

P. cinereum

Idn-2

+

HE614599

+

HE655074

P. compressum

CJ1-1

+

HE614600

+

HE655075

P. compressum

Cr-1

+

HE614601

+

HE655076

P. didermoides

--

+

HE614602

+

HE655077

P. flavicomum

UFF1

+

HE614611

+

X78959

P. polycephalum

Wis1

+

X13160

+

X60211

P. pusillum

CJA3

+

HE614603

+

HE655078

P. rigidum

ATCC22485

+

HE614604

+

HE655079

P. roseum

C1

+

HE614605

+

HE655080

(a) SSU and LSU rDNA sequences.

(b) GenBank/EMBL/DDJB accession numbers.

--, no isolate name given; +, analysed; na, not analysed.

Recent taxonomy updates recognize at least 80 and 20 species of Physarum and Badhamia, respectively [2, 32]. Myxomycete taxonomy has traditionally been based on morphological characteristics of sporocarps and fruiting bodies. However, it has been debated for decades what characters that unambiguously distinguish Physarum and Badhamia genera species [33, 34], and consequently some Physarum species in the current taxonomy has previously been assigned to the Badhamia genus and vice versa [32]. A similar relationship based on rRNA gene sequence phylogeny involving one Badhamia and two Physarum isolates was recently noted, but not further commented [15]. Our findings of a polyphyletic origin of the Physarum genus, with at least three phylogenetic clades interspersed with Badhamia isolates, strongly suggest a reevaluation of the current taxonomy to also include molecular data.

Intron synapomorphy supports Physarumclades

A different approach to validate the Physarum clades is structural analysis of the obligatory group I introns present in the LSU rRNA gene of all investigated Physarales isolates (Additional file 3: Table S2 and Additional file 4: Table S3). Both L1949 and L2449 introns have strict vertical inheritance pattern within the Physarales with potential as genetic markers [11, 24]. We analyzed structural features of L2449 introns at the RNA level in more detail, and representative secondary structure diagrams of the catalytic core are presented in Figure 4A-C. Whereas the approximately 120-nt catalytic core was found conserved in sequence and structure, large extensions of various lengths were observed in all peripheral paired segment regions (P1, P2, P5, P6, P8 and P9). An extreme case was noted in the L2449 intron (2483 nt) in P. didermoides (Figure 4C). Here, direct repeat motifs were found both within the 503 nt P2 extension and the 1268 nt P9 extension (Figure 4D), and exemplifies the complex structural nature of myxomycete group I introns.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2148-12-166/MediaObjects/12862_2012_Article_2143_Fig4_HTML.jpg
Figure 4

RNA structure diagrams of L2449 group I introns in Physarales clades. Paired elements (P1 to P9) are indicated according to the consensus of L2449 group I introns in myxomycetes [11]. The intron core and exon sequences are shown in upper case and lower case letters, respectively. Arrows indicate the 5’ splice site (SS) and the 3’ SS. Sizes of extension sequences at peripheral P-regions are indicated. The synapomorphy at the intron G-binding site in P7 is boxed. A. Clade I (P. pusillum CJA3, B. melanospora Pr-1) L2449 introns that contain the G-binding site motif UCG-CA (green box). B. Clade II (P. polycephalum Wis1, B. utricularis) L2449 introns that contain the G-binding site motif ACG-CU (yellow box). C. Clade III (P. compressum CJ1-1, P. didermoides) L2449 introns that contain the G-binding site motif AUG-CU (blue box). D. Schematic organization of large extension sequences in P2 and P9 of the 2483 nt L2449 intron in P. didermoides (Clade III). The 503 nt P2 extension contains five copies of an approximately 60-nt A motif. The 1268 nt P9 extension has a more complex organization of three different motifs of approximately 40 nt (B motif, 3 copies), 100 nt (C motif, 5 copies) and 55 nt (D motif, 2 copies), respectively. Here, motif C is composed by three sub-motifs (a-b-a-c), and related to motif D (a-c).

The most surprising structural feature, however, was found directly at the guanosine binding site (G-binding site) within the P7 segment. The G-binding site is at the catalytic center of a group I ribozyme, is highly conserved among all group I introns, and contains a universally conserved G:C pair adjacent to a bulged nucleotide (C or A). The G-binding site follows strong sequence co-variation rules, and if the bulged nucleotide is a C the first pair of the P7 segment is always an A:U or a U:A pair [35]. Figure 4A presents structures of Clade I introns (P. pusillum and B. melanospora) which have a U:A pair followed by a bulged C. Clade II introns (P. polycephalum and B. utricularis) have a slightly different G-binding site of an A:U pair followed by a bulged C (Figure 4B). A third G-binding site variant was noted in the Clade III introns (P. compressum and P. didermoides; Figure 4C). Similar to Clade II introns, the Clade III introns have a A:U as the first P7 pair, but followed by a bulged U. The latter feature is highly unusual among group I introns and to our knowledge the Physarum Clade III L2449 introns are the only examples among the more than 22.000 group I introns available in the sequence data base [36, 37]. The fact that G-binding site variants appears to correlate with sequence-based phylogeny suggests that the intron structure character represents a synapomorphic character that further supports a complex origin of the Physarum genus.

Conclusions

Partial and complete rRNA gene sequences have been obtained from a large number of myxomycete isolates and have proven to represent reliable phylogenetic markers. However, myxomycete rRNA gene sequencing is challenged by an extreme high load of group I intron elements, and we report in our data sets 222 introns interrupting coding sequences in Physarales isolates. Several of the introns have complex organisations and were thus difficult to identify at sequence level. Introns are typically located in highly conserved parts of the rRNA genes that directly interfere with the design of universal primers and may result in sampling biases in analysis approaches. We observed congruence between phylogeny and current taxonomy of the Didymiaceae isolates, but a different pattern was observed among the family Physaraceae and includes a polyphyletic origin of the genus Physarum. Three well supported clades containing a mixture of Physarum and Badhamia isolates were noted, suggesting reevaluation of the current taxonomy into new genera. These observations were further supported by an unusual synapomorphic character at the G-binding site of the L2449 group I introns.

Methods

Myxomycetes isolates and culturing

All new isolates reported were collected as sporocarps. Their geographical origin, mode of classification, name of the collector, and if a cell line was obtained, are listed in Additional file 8: Table S4. Cultivation of amoebae was essentially performed as previously described [12, 38].

DNA extraction and sequencing

Total DNA was extracted from growing amoebae as previously described [11] or directly from sporocarps without culturing (summarized in Additional file 8: Table S4). Here, DNA extraction was performed using Berlin Technologies instrument (PreCelly’s 24 tube # VK05). Two to five sporocarp heads were homogenized in 250 μl lysis buffer (4 M Guanidine thiocyanate, 50 mM Tris–HCl (pH 7.5), 10 mM EDTA, 2% SDS, 1% β-mercaptoethanol) at 5000 rpm in 30 seconds for 1 to 4 cycles. The lysed spores were treated twice with phenol-chloroform, and DNA was extracted by ethanol precipitation and re-suspended in TE buffer. The SSU and LSU rRNA gene segments were PCR-amplified from total DNA extracts in several overlapping fragments using multiple sets of specific oligonucleotides (contact SDJ for details). DNA sequences were performed on both strands using automatic sequencing (Big Dye terminator chemistry, Applied Biosystems, Foster City, CA, USA).

Phylogenetic analyses

Three data sets based on nuclear rRNA sequences were used in the molecular phylogeny. The SSU dataset (1028 nucleotide positions), the combined SSU/LSU data set (2522 nucleotide positions), and the LSU data set (765 nucleotide positions) were manually aligned in BioEdit v.7.0.5.3 [39] based on the secondary structure models of P. polycephalum and D. iridis[22, 40]. Phylogenetic trees were built with the methods of NJ applying the Jukes-Cantor nucleotide substitution model based on the number of nucleotide substitution per site [41], and MP using heuristic search with close-neighbor-interchange (CNI) level 3 and generation of 10 random initial trees, in MEGA version 4.0.2 [42]., as well as ML using PhyML interface [43] and nucleotide substitution models selected by jModelTest 0.1 package [44]. The following substitution models were used in ML: TPM2uf + I + G (SSU data set), GTR + I + G (SSU/LSU dataset), TNef + I (LSU dataset; Didymiaceae) and GTR + I + G (LSU dataset; Physaraceae). The reliabilities of tree branching points for NJ, MP and ML trees were evaluated by bootstrap analyses (2000 replications). BAY analyses were performed to reconstruct phylogenetic trees from all datasets using MrBayes, version 3.1.2 [45, 46] with two independent runs and the Metropolis-coupled Markov chain Monte Carlo (MCMCMC) method. Here the evolutionary model GTR + I + G was used for all datasets. A total of 1,000,000 generations were run with sampling every 100 generations. For each dataset the standard deviation of split frequencies was below 0.015 at the end of the run. Twenty-five % of initial trees were discarded and a consensus tree with posterior probabilities was generated from the remaining 15,000 trees.

Declarations

Acknowledgments

We thank Dr. J. Clark (University of Kentucky, USA), Iolanda and Giovanni Manavella (Italy), Dr. G. Shipley (University of Texas, USA), Marianne Meyer (France) and Dr. P. Haugen (University of Tromsø, Norway) for providing various myxomycete sporocarps and cell-lines. Dr. Finn Haugli is acknowledged for comments on the manuscript. This work was supported by grants to SDJ from the Norwegian Research Council and the University of Tromsø.

Authors’ Affiliations

(1)
RNA and Transcriptomics group, Department of Medical Biology, Faculty of Health Sciences, University of Tromsø
(2)
Department of Biology, University of Washington

References

  1. Olive LS, Stoianovitch C: The Mycetozoans. 1975, New York: Academic PressGoogle Scholar
  2. Poulain M, Meyer M, Bozonnet J: Les Myxomycetes. 2011Google Scholar
  3. Dee J: Genes and development in Physarum. Trends Genet. 1987, 3: 208-213.View ArticleGoogle Scholar
  4. Einvik C, Elde M, Johansen S: Group I twintrons: genetic elements in myxomycete and schizopyrenid amoeboflagellate ribosomal DNAs. J Biotechnol. 1998, 64: 63-74. 10.1016/S0168-1656(98)00104-7.PubMedView ArticleGoogle Scholar
  5. Stephenson SL, Fiore-Donno AM, Schnittler M: Myxomycetes in soil. Soil Biol Biochem. 2011, 43: 2237-2242. 10.1016/j.soilbio.2011.07.007.View ArticleGoogle Scholar
  6. Stephenson SL, Stempen H: Myxomycetes: a handbook of slime molds. 1994, Portland, Oregon: TimberGoogle Scholar
  7. Bonner JT: Brainless behavior: a myxomycete chooses a balanced diet. Proc Natl Acad Sci USA. 2010, 107: 5267-5268. 10.1073/pnas.1000861107.PubMedPubMed CentralView ArticleGoogle Scholar
  8. Ko TW, Stephenson SL, Hyde KD, Carlos R, Lumyong S: Patterns of occurrence of myxomycetes on lianas. Fungal Ecol. 2010, 3: 302-310. 10.1016/j.funeco.2009.11.005.View ArticleGoogle Scholar
  9. Takahashi K, Hada Y: Geographical distribution of myxomycetes on coniferous deadwood in relation to air temperature in Japan. Mycoscience. 2010, 51: 281-290. 10.1007/s10267-010-0044-9.View ArticleGoogle Scholar
  10. Martin MP, Lado C, Johansen S: Primers are designed for amplification and direct sequencing of ITS region of rDNA from myxomycetes. Mycologia. 2003, 95: 474-479. 10.2307/3761889.PubMedView ArticleGoogle Scholar
  11. Wikmark OG, Haugen P, Haugli K, Johansen SD: Obligatory group I introns with unusual features at positions 1949 and 2449 in nuclear LSU rDNA of Didymiaceae myxomycetes. Mol Phylogenet Evol. 2007, 43: 596-604. 10.1016/j.ympev.2006.11.004.PubMedView ArticleGoogle Scholar
  12. Wikmark OG, Haugen P, Lundblad EW, Haugli K, Johansen SD: The molecular evolution and structural organization of group I introns at position 1389 in nuclear small subunit rDNA of myxomycetes. J Eukaryot Microbiol. 2007, 54: 49-56. 10.1111/j.1550-7408.2006.00145.x.PubMedView ArticleGoogle Scholar
  13. Fiore-Donno AM, Berney C, Pawlowski J, Baldauf SL: Higher-order phylogeny of plasmodial slime molds (Myxogastria) based on elongation factor 1-A and small subunit rRNA gene sequences. J Eukaryot Microbiol. 2005, 52: 201-210. 10.1111/j.1550-7408.2005.00032.x.PubMedView ArticleGoogle Scholar
  14. Fiore-Donno AM, Meyer M, Baldauf SL, Pawlowski J: Evolution of dark-spored Myxomycetes (slime-molds): molecules versus morphology. Mol Phylogenet Evol. 2008, 46: 878-889. 10.1016/j.ympev.2007.12.011.PubMedView ArticleGoogle Scholar
  15. Fiore-Donno AM, Nikolaev SI, Nelson M, Pawlowski J, Cavalier-Smith T, Baldauf SL: Deep phylogeny and evolution of slime moulds (mycetozoa). Protist. 2010, 161: 55-70. 10.1016/j.protis.2009.05.002.PubMedView ArticleGoogle Scholar
  16. Fiore-Donno AM, Kamono A, Chao EE, Fukui M, Cavalier-Smith T: Invalidation of Hyperamoeba by transferring its species to other genera of Myxogastria. J Eukaryot Microbiol. 2010, 57: 189-196. 10.1111/j.1550-7408.2009.00466.x.PubMedView ArticleGoogle Scholar
  17. Hoppe T, Kutschera U: In the shadow of Darwin: Anton de Bary's origin of myxomycetology and a molecular phylogeny of the plasmodial slime molds. Theory Biosci. 2010, 129: 15-23. 10.1007/s12064-009-0079-7.PubMedView ArticleGoogle Scholar
  18. Aldrich HC, Daniel JW: Cell biology of Physarum and Didymium. 1982, New York: AcademicGoogle Scholar
  19. Vogt VM, Braun R: Structure of Ribosomal DNA in Physarum polycephalum. J Mol Biol. 1976, 106: 567-587. 10.1016/0022-2836(76)90252-7.PubMedView ArticleGoogle Scholar
  20. Ferris PJ: Structure and inheritance of Physarumribosomal DNA.PhD thesis. 1984, Ithaca, NY: Cornell UniversityGoogle Scholar
  21. Silliker ME, Collins OR: Non-mendelian inheritance of mitochondrial DNA and ribosomal DNA in the myxomycete, Didymium iridis. Mol Gen Genet. 1988, 213: 370-378. 10.1007/BF00339605.PubMedView ArticleGoogle Scholar
  22. Johansen S, Johansen T, Haugli F: Extrachromosomal ribosomal DNA of Didymium iridis: sequence analysis of the large subunit ribosomal RNA gene and sub-telomeric region. Curr Genet. 1992, 22: 305-312. 10.1007/BF00317926.PubMedView ArticleGoogle Scholar
  23. Vader A, Naess J, Haugli K, Haugli F, Johansen S: Nucleolar introns from Physarum flavicomum contain insertion elements that may explain how mobile group I introns gained their open reading frames. Nucleic Acids Res. 1994, 22: 4553-4559. 10.1093/nar/22.22.4553.PubMedPubMed CentralView ArticleGoogle Scholar
  24. Lundblad EW, Einvik C, Rønning S, Haugli K, Johansen S: Twelve Group I introns in the same pre-rRNA transcript of the myxomycete Fuligo septica: RNA processing and evolution. Mol Biol Evol. 2004, 21: 1283-1293. 10.1093/molbev/msh126.PubMedView ArticleGoogle Scholar
  25. Nielsen H, Johansen SD: Group I introns: Moving in new directions. RNA Biol. 2009, 6: 375-383. 10.4161/rna.6.4.9334.PubMedView ArticleGoogle Scholar
  26. Vicens Q, Cech TR: Atomic level architecture of group I introns revealed. Trends Biochem Sci. 2006, 31: 41-51. 10.1016/j.tibs.2005.11.008.PubMedView ArticleGoogle Scholar
  27. Johansen S, Johansen T, Haugli F: Structure and evolution of myxomycete nuclear group I introns: a model for horizontal transfer by intron homing. Curr Genet. 1992, 22: 297-304. 10.1007/BF00317925.PubMedView ArticleGoogle Scholar
  28. Haugen P, Reeb V, Lutzoni F, Bhattacharya D: The evolution of homing endonuclease genes and group I introns in nuclear rDNA. Mol Biol Evol. 2004, 21: 129-140.PubMedView ArticleGoogle Scholar
  29. Johansen SD, Haugen P, Nielsen H: Expression of protein-coding genes embedded in ribosomal DNA. Biol Chem. 2007, 388: 679-686.PubMedView ArticleGoogle Scholar
  30. Pawlowski J, Burki F: Untangling the phylogeny of amoeboid protists. J Eukaryot Microbiol. 2009, 56: 16-25. 10.1111/j.1550-7408.2008.00379.x.PubMedView ArticleGoogle Scholar
  31. Lahr DJ, Grant J, Nguyen T, Lin JH, Katz LA: Comprehensive phylogenetic reconstruction of amoebozoa based on concatenated analyses of SSU-rDNA and actin genes. PLoS One. 2011, 6: e22780-10.1371/journal.pone.0022780.PubMedPubMed CentralView ArticleGoogle Scholar
  32. Neubert H, Nowotny W, Baumann K: Die Myxomyceten. 1995, Gomaringen: Karlheinz Baumann VerlagGoogle Scholar
  33. Alexopoulos CJ: Morphology, taxonomy, and phylogeny. Cell biology of Physarum and Didymium. Edited by: Aldrich HC, Daniel JW. 1982, New York: Academic Press, 3-21.View ArticleGoogle Scholar
  34. Clark J, Haskins EF, Stephenson SL: Biosystematics of the myxomycete Badhamia gracilis. Mycologia. 2003, 95: 104-108. 10.2307/3761967.PubMedView ArticleGoogle Scholar
  35. Michel F, Westhof E: Modelling of the three-dimensional architecture of group I catalytic introns based on comparative sequence analysis. J Mol Biol. 1990, 216: 585-610. 10.1016/0022-2836(90)90386-Z.PubMedView ArticleGoogle Scholar
  36. Zhou Y, Lu C, Wu Q-J, Wang Y, Sun Z-T, Deng J-C, Zhang Y: GISSD: group I intron sequence and structure database. Nucleic Acids Res. 2008, 36: D31-D37. 10.1093/nar/gkn052.PubMedPubMed CentralView ArticleGoogle Scholar
  37. Gardner PP, Daub J, Tate JG, Moore BL, Osuch IH, Griffiths-Jones S, Finn RD, Nawrocki EP, Kolbe DL, Eddy SR, Bateman A: Rfam: Wikipedia, clans and the “decimal” release. Nucleic Acids Res. 2011, 39: D131-D145. 10.1093/nar/gkr609.View ArticleGoogle Scholar
  38. Johansen S, Elde M, Vader A, Haugen P, Haugli K, Haugli F: In vivo mobility of a group I twintron in nuclear ribosomal DNA of the myxomycete Didymium iridis. Mol Microbiol. 1997, 24: 737-745. 10.1046/j.1365-2958.1997.3921743.x.PubMedView ArticleGoogle Scholar
  39. Hall TA: BioEdit: a user friendly sequence alignment editor and analysis program for windows 95/98/NT. Nucl Acids Symp Ser. 1999, 41: 95-98.Google Scholar
  40. Johansen T, Johansen S, Haugli FB: Nucleotide sequence of Physarum polycephalum small subunit ribosomal RNA as inferred from the gene sequence: secondary structure and evolutionary implications. Curr Genet. 1988, 14: 265-273. 10.1007/BF00376747.PubMedView ArticleGoogle Scholar
  41. Nei M, Kumar S: Molecular evolution and phylogenetics. 2000, New York: Oxford University PressGoogle Scholar
  42. Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.PubMedView ArticleGoogle Scholar
  43. Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.PubMedView ArticleGoogle Scholar
  44. Posada D: jModelTest: phylogenetic model averaging. Mol Biol Evol. 2008, 25: 1253-1256. 10.1093/molbev/msn083.PubMedView ArticleGoogle Scholar
  45. Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.PubMedView ArticleGoogle Scholar
  46. Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.PubMedView ArticleGoogle Scholar

Copyright

© Nandipati et al.; licensee BioMed Central Ltd. 2012

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Advertisement