A unique horizontal gene transfer event has provided the octocoral mitochondrial genome with an active mismatch repair gene that has potential for an unusual self-contained function
© Bilewitch and Degnan; licensee BioMed Central Ltd. 2011
Received: 27 April 2011
Accepted: 29 July 2011
Published: 29 July 2011
The mitochondrial genome of the Octocorallia has several characteristics atypical for metazoans, including a novel gene suggested to function in DNA repair. This mtMutS gene is favored for octocoral molecular systematics, due to its high information content. Several hypotheses concerning the origins of mtMutS have been proposed, and remain equivocal, although current weight of support is for a horizontal gene transfer from either an epsilonproteobacterium or a large DNA virus. Here we present new and compelling evidence on the evolutionary origin of mtMutS, and provide the very first data on its activity, functional capacity and stability within the octocoral mitochondrial genome.
The mtMutS gene has the expected conserved amino acids, protein domains and predicted tertiary protein structure. Phylogenetic analysis indicates that mtMutS is not a member of the MSH family and therefore not of eukaryotic origin. MtMutS clusters closely with representatives of the MutS7 lineage; further support for this relationship derives from the sharing of a C-terminal endonuclease domain that confers a self-contained mismatch repair function. Gene expression analyses confirm that mtMutS is actively transcribed in octocorals. Rates of mitochondrial gene evolution in mtMutS-containing octocorals are lower than in their hexacoral sister-group, which lacks the gene, although paradoxically the mtMutS gene itself has higher rates of mutation than other octocoral mitochondrial genes.
The octocoral mtMutS gene is active and codes for a protein with all the necessary components for DNA mismatch repair. A lower rate of mitochondrial evolution, and the presence of a nicking endonuclease domain, both indirectly support a theory of self-sufficient DNA mismatch repair within the octocoral mitochondrion. The ancestral affinity of mtMutS to non-eukaryotic MutS7 provides compelling support for an origin by horizontal gene transfer. The immediate vector of transmission into octocorals can be attributed to either an epsilonproteobacterium in an endosymbiotic association or to a viral infection, although DNA viruses are not currently known to infect both bacteria and eukaryotes, nor mitochondria in particular. In consolidating the first known case of HGT into an animal mitochondrial genome, these findings suggest the need for reconsideration of the means by which metazoan mitochondrial genomes evolve.
Animal mitochondrial genomes are generally conserved in structure, approximate size and gene content. The circular genome of the so-called higher metazoans typically encodes 13 proteins, 22 transfer RNAs and two ribosomal RNAs, while lacking introns [1, 2]. Variation in genome size usually derives only from length variation in the single large non-coding region, so that most mitochondrial genomes range in size from 14 to 17 kb [1, 2]. The phylum Cnidaria, however, provides many interesting and diverse exceptions to this pattern. For example, medusozoans (jellyfish and hydroids; Classes Cubozoa, Scyphozoa, Staurozoa and Hydrozoa) are unique among animals in having a linearized mitochondrial DNA (mtDNA) [3, 4]. Hexacorals (Anthozoa: Hexacorallia) are unique in that their mitochondrial COI and ND5 genes contain introns, and anthozoan mitochondrial genomes on the whole are known to contain a greatly reduced complement of transfer RNA genes (only 1 or 2 tRNAs) [5–9]. Among the octocorals (Anthozoa: Octocorallia), the mitochondrial genome is known to have undergone inversion and reorganization at least three times [10, 11], including twice within a single family . These changes in octocoral gene order appear to have resulted from intramitochondrial recombination , which itself was believed to be unlikely in animals until its discovery in 1997 . As more mitochondrial genome data become available for cnidarians, and for the Octocorallia in particular, it seems that fewer generalizations of mitochondrial structure and evolution seem applicable to these basal metazoans.
Indeed, when the first complete octocoral mitochondrial genomes were sequenced in 1995, another surprise was revealed. The mtDNA of Renilla kolikeri and Sarcophyton glaucum (Anthozoa: Octocorallia) [5, 13] was observed to contain the standard anthozoan complement of oxidative metabolic genes - 13 protein-coding genes, two rRNAs and a single tRNA - but both genomes also contained an unexpected open reading frame (ORF) [5, 14, 15]. Almost 16% of the mtDNA (2949/18453 bp) was occupied by a novel ORF that to this day has not been found in any mitochondrial genome outside of the Octocorallia . Based on 19.7% similarity of the translated sequence to the yeast nuclear Mutation Suppressor Homolog 1 (MSH1) gene, it was suggested that the novel ORF might be active in DNA mismatch repair (MMR) and thus it was first termed 'mtMSH' as a mitochondrial homolog of the eukaryotic MSH family . Alternative names have since been used for this octocoral gene to infer relationship specifically to the MSH1 family ("msh1" [16, 17]) or to encompass broad taxonomic ranges ("mtMutS" ), including MMR genes from both prokaryotes (MutS) and eukaryotes (MSH). We hereafter preferentially use the name 'mtMutS' to refer to the octocoral mitochondrial ORF, since it makes the fewest a priori assumptions of homology to either eukaryotic MSH or prokaryotic MutS MMR gene families, specifically.
Paradoxically, while MutS and MSH function in reducing error rates in DNA replication in eubacteria and eukaryotes, respectively , mtMutS itself possesses the highest level of systematically-informative variation amongst any protein-coding mitochondrial gene examined to date within octocorals . As such, it has played a central role in molecular systematic efforts to revise the taxonomic classification of the Octocorallia. Its utility in phylogenetics was first demonstrated in a species-level analysis  and has subsequently been extended to more comprehensive studies of octocoral genera  and subfamilies . Revisionary systematics and the classification of newly described species have also used sequences of mtMutS, either alone  or in combination with nuclear [23, 24] or other mitochondrial [25, 26] sequences. Despite being rapidly embraced as a phylogenetic tool, however, mtMutS is the only animal mitochondrial gene for which both functionality and origin remain equivocal.
A more recent analysis proposed the closest relative of the octocoral mtMutS to be the MutS sequence from the Acanthamoeba polyphaga nucleocytoplasmic large DNA virus (NCLDV or 'girus'), closely followed by MutS sequences from epsilonproteobacteria . Although no other eukaryotic sequences were included in the analysis, the high level of identity between MutS from the NCLDV and mtMutS (25%), and the presence of a C-terminal HNH endonuclease domain in both, were taken as evidence of common ancestry. This raised the hypothesis that octocorals might have acquired mtMutS via horizontal gene transfer (HGT) from a NCLDV ancestor into the common ancestor of Octocorallia (Figure 1C). An HGT origin of mtMutS is further supported by the most recent and most comprehensive analysis of both MSH and MutS evolution, which included five more NCLDV sequences . This study provided very strong support for the lineage containing octocoral, epsilonproteobacterial and viral MutS sequences, but was unable to discriminate between the bacterial or the viral lineage as the sister group to mtMutS . Thus the most recent common ancestor to mtMutS remains unknown, although it appears almost certain to be either viral or bacterial - not eukaryotic - in origin.
Here we provide further tests of the HGT hypothesis of origin by examining the homology, diversification and evolutionary patterns of the mtMutS gene from multiple octocorals. First, we reconstruct MutS phylogenies using all conserved sections of the gene, rather than only the single domain used previously , thus providing more robust support for the HGT origin. Horizontal gene transfer of a mtMutS precursor from a non-metazoan donor in turn raises multiple questions pertaining to the acquisition, function and evolution of the octocoral gene. In particular, evolutionary theory predicts that the successful integration of HGT- acquired genes into a host genome requires both gene activity and gene longevity . Second, in addressing these predictions, we ask whether the octocoral mtMutS gene is active, and how it has diversified over evolutionary timescales within the foreign host genome . In combination, these analyses significantly expand our understanding of a novel gene that has persisted within octocorals since their most recent common ancestor.
A HGT origin of mtMutSin the common ancestor of Octocorallia
PSI-BLAST returned a set of epsilonproteobacteria and nucleocytoplasmic large double-stranded DNA virus (NCLDV) MutS hits as best matches to complete, translated octocoral mtMutS sequences present in GenBank.
All eukaryotic sequences, with the exception of MSH1, formed gene family-specific monophyletic MSH clades with over 95% support (Figure 2, purple oval). Of all the MSH families, MSH6 shared the most recent common ancestor with octocoral mtMutS in our Bayesian phylogeny. However, none of the eukaryotic MSH families formed a monophyletic lineage with octocorals and, with the exception of MSH2-MSH4-MSH5, formed polyphyletic inter-familial relationships amongst themselves. Additionally, the two included sequences of eukaryotic MSH1 (one each from Nematostella and yeast) formed polyphyletic associations amongst the eubacterial MutS sequences at the distal end of the Bayesian tree to octocoral mtMutS (Figure 2, red oval), indicating they are highly divergent. These arrangements preclude hypotheses of an immediate common ancestry between mtMutS and MSH1 (Figure 1A,B.)
A lack of monophyly between mtMutS and any MSH family, combined with the phylogenetic arrangement of mtMutS within the MutS7 lineage (Figure 2, grey oval) provides further support to a hypothesis of non-eukaryotic origin of mtMutS in the most recent common ancestor of extant octocorals. The sister organism(s) from which the mtMutS precursor was obtained was not identified with confidence in the phylogeny, but the placement of mtMutS within MutS7 strongly supports a hypothesis of horizontal gene transfer.
MtMutSpredicts a functional protein structure
As previously reported for S. glaucum , the HNH endonuclease domain detected by SMART at the C-terminal end of all octocoral mtMutS sequences contains characteristic histidine-histidine-asparagine-histidine (HHNH) amino acids at positions 1676, 1677, 1709 and 1718, respectively (Figure 3). These residues are conserved with HNH-containing MutS sequences found among the epsilon-proteobacteria and NCLDV sequences. HNH domains were not detected by SMART in any other queried MSH or MutS protein sequences, providing further evidence of common ancestry of mtMutS and other MutS7 genes.
MtMutsis expressed and available for translation
MtMutShas evidence of higher rates of mutation than other octocoral mt genes
Alignments of 1597 bp in 169 taxa, 3967 bp in 410 taxa and 3426 bp in 509 taxa were generated for octocoral COI, hexacoral COI and octocoral mtMutS mtDNA sequences, respectively, using all available data present in GenBank. Alignments of 1737 bp in 224 taxa and 1871 bp in 145 taxa were generated for hexacoral and octocoral 18S nuclear DNA sequences, respectively. Mann-Whitney U-tests indicated highly-significant differences (p < 0.0001) between the means of pairwise genetic differences in mtDNA sequence comparisons of octocoral mtMutS vs. octocoral COI and octocoral COI vs. hexacoral COI and between nuclear DNA sequence comparisons of octocoral 18S vs. hexacoral 18S. Tests using only taxa with less than 10% missing data in each comparison yielded identical significant results (data not shown), thus the inclusion of taxa with incomplete gene sequences did not alter the results of the full comparisons.
Among the nuclear 18S sequence comparisons, octocorals display a significantly higher mean pairwise distance than hexacorals (p < 0.0001), opposite to the pattern observed for the mitochondrial genome. Both octocoral and hexacoral 18S distributions were positively skewed and displayed, respectively, the highest and second-highest distances out of all examined datasets. However, both 18S datasets had modal values and the majority of their distributions in a range of distances similar to the distribution of octocoral COI (Figure 7).
These patterns indicate an overall reduced rate of evolution of octocoral COI relative to hexacoral COI since the time that they diverged from their last common ancestor. Given the corroboration of previous observations of low divergence levels in other mitochondrial genes [20, 40, 41], we may extrapolate our observation to the mitochondrial genome as a whole. However, the mtMutS gene itself displays significantly higher pairwise distances than COI, although its modal divergence level is still low compared to that observed between inter-ordinal taxa of scleractinians (Figure 7). Furthermore, the results of the 18S comparison indicate that the reduction in rates of octocoral mitochondrial mutation may not extend to the nuclear genome. Octocorals may, conversely, display higher rates of nuclear mutation than hexacorals, as displayed for our 18S comparison (Figure 7).
Non-eukaryotic origins of mtMutSand mitochondrial HGT
The mitochondrial targeting of MSH1 plus initial observations of sequence similarity between MSH1 and mtMutS formed the basis of early conclusions that the two genes were orthologous . However, this phylogenetic association appears to be artifactual, due to insufficient sequence sampling from prokaryotic sources in these initial studies. The continued expansion of available MutS sequences from an array of prokaryotes allowed a new phylogenetic arrangement for mtMutS - one that placed it in an association with the epsilonproteobacteria and giant double-stranded nucleo-cytoplasmic viruses (NCLDV's or "giruses")  in one of four new MutS lineages, namely MutS7 . Members of the MutS7 lineage are united by Domains I, III and V of MutS plus a C-terminal HNH endonuclease domain. In contrast, MSH1 and other members of its MutS1 lineage contain MutS/MSH domains II and IV and lack an endonuclease domain . Our phylogenetic analysis provides further support for this placement of octocoral mtMutS within MutS7 with high levels of confidence (Figure 2), and provides conclusive evidence that the octocoral mtMutS gene does not share an immediate common ancestor with any eukaryotic MSH family, including MSH1. The exclusion of the HNH domain from our phylogenetic analysis indicates the tree topology is based upon evolution of homologous MutS domains and not on this MutS7 synapomorphy. If mtMutS is active in mitochondrial MMR, this similarity in function with MSH1 would likely be the result of convergent evolution between paralogous genes. For these reasons, the octocoral mitochondrial MutS gene should no longer be referred to as 'msh1', as in the current systematic literature (e.g. [20, 25, 41]), since it carries the connotation of eukaryotic origins. We recommend the adoption of the name mtMutS, after studies which first identified it as a unique mitochondrial MutS gene (this study; [30, 31]). Even more specifically, mtMutS7 could be used to further identify it as a member of the MutS7 HNH endonuclease-containing family .
The evidence for a non-eukaryotic origin of mtMutS in octocorals rules out the possibility of either a mitochondrion-to-nucleus  or nucleus-to-mitochondrion  translocation hypothesis in its relationship to MSH families (see Figure 1B,C). This also necessitates a reinterpretation of the origins of MSH in eukaryotes. It has been proposed that the endosymbiotic incorporation of a mitochondrion-like prokaryotic cell brought MutS into the primordial eukaryotic cell (as MSH1), with translocation to the nucleus and a series of gene duplications producing the remaining MSH families . While MSH1 was proposed to be a basal family among all MSHs, evidence for it once existing within the mitochondrial genome was based on apparently artifactual phylogenetic associations with mtMutS. In addition to our findings (Figure 2), more comprehensive MutS/MSH phylogenies [30, 33] also indicate that MSH1/MutS1 is a highly-derived lineage relative to the other eight families of prokaryotic MutS and six families of eukaryotic MSH. Thus, although a prokaryotic origin of MSH is supported by all current phylogenetic hypotheses [18, 30, 33], evidence for an intermediary role of the mitochondrion in this process is lacking.
Determining the vector of mtMutSHGT
The affinity of a monophyletic octocoral mtMutS lineage to the MutS7 family containing MutS from epsilonproteobacteria and NCLDVs supports the hypothesis of horizontal gene transfer (HGT), as initially suggested  and later evidenced [32, 33]. The presence of the HNH endonuclease domain supports a common ancestry of MutS in these three disparate groups of organisms [32, 33], but the direction and means by which the gene was distributed remain questionable. The MutS7 lineage could have origins either within the basal Eukaryota or the Eubacteria, given that it occupies a clade in-between the MutS1/MSH lineages and the other seven MutS families [: Figure 1]. Within the MutS7 lineage, the branching order of the octocorallian, bacterial and viral MutS lineages can provide clues into the origins of the gene family. Unfortunately, neither the previous comprehensive phylogenetic analysis  nor our current study produced phylogenies with strongly supported topologies for the branching order of MutS from Octocorallia, NCLDVs and epsilonproteobacteria. The octocoral mtMutS and epsilonproteobacterial MutS lineages are both monophyletic with high posterior probability, as is the entire MutS7 lineage, but whether the NCLDV, the epsilonproteobacteria or both are the sister group to octocorallian mtMutS can not yet be determined (Figure 2; [:Figures 1a, S2]). Although less taxonomically-comprehensive for octocorals, a previous study [: Figure 5] provided 89% bootstrap support for A. polyphaga mimiviral MutS as the sister lineage to octocoral mtMutS, with the epsilonproteobacteria forming a basal MutS7 divergence. On this finding, they proposed a scenario whereby mtMutS was transferred from a microbial host into the mitochondrion of the common octocorallian ancestor via a NCLDV vector.
Although viruses infecting the mitochondrial genome ('mitoviruses') have been documented, they have so far been limited to single- or double-stranded RNA viruses and are so far only known to infect fungi . Double-stranded DNA viruses, such as the NCLDVs, are not known to infect host mitochondria, thus a mechanism by which a viral MutS could have entered an octocoral mitochondrion is currently lacking. Interestingly, mitochondria are known to migrate to viral assembly sites in infected primate cells [43, 44]. If similar responses occurred in NCLDV-infected cells, it would place the host mitochondrial genome in proximity to viral DNA. At this time, however, it is unknown whether any MutS7-containing NCLDV is capable of infecting an octocorallian cell. Furthermore, although NCLDVs infect a wide range of hosts , none of them are bacterial and no examples exist of any virus that directly infects both prokaryotic and eukaryotic cells.
An alternative means for octocorals acquiring mtMutS through HGT may not involve NCLDVs at all. Diverse arrays of bacterial genotypes are known to associate with the surfaces and tissues of marine invertebrates, including anthozoans [45, 46]. These microbial communities can include endosymbionts [47, 48], which may confer advantages to host cells such as the fixation of nitrogen . HGT between bacteria is thought to be common in the marine environment  and has been reported to have occurred multiple times between bacteria and nematodes in terrestrial systems . A similar association between a MutS7-containing bacterium and an octocoral may have been responsible for HGT of mtMutS. Epsilonproteobacteria are common constituents of deep-sea environments [50, 51], which also contain diverse octocoral communities (e.g.), therefore associations between the two are possible. The fact that such an occurrence has not yet been documented may simply reflect the lack of study of non-pathogenic bacterial associations in octocorals.
The functionality and function of mtMutS
Evolutionary theory predicts that a successful HGT integration requires that the transferred gene be active in the new host . Our results have demonstrated for the first time that the octocoral mtMutS gene contains all coding necessary for a functional MutS protein. The presence of conserved residues for mismatch recognition in Domain I and ATPase-specific residues in Domain V provide primary evidence for functionality (Figure 3). The presence and conserved order of three MutS protein domains provides secondary evidence (Figure 4) and the conserved conformation model (Figure 5) provides tertiary evidence that mtMutS has the potential to play an active role in the octocoral proteome. Finally, our transcriptional analysis indicates that mtMutS is transcribed and available as mRNA for translation (Figure 6).
Successful integration of a horizontally transferred gene into a host genome also requires longevity over evolutionary time . Phylogenetic analyses clearly reveal that the mtMutS gene has persisted in octocorals since at least the last common ancestor of all extant species (Figure 2; see also ). Given the reduction of gene content in the octocoral mitochondrial genome [2, 5] and the frequency of its reorganization through intramitochondrial recombination events , it can be deduced that there is strong selection for the continued presence of mtMutS. Interestingly, no octocoral has been found to lack mtMutS even though the very existence of the hexacoral sister group demonstrates that anthozoans are capable of persisting over geologic time spans without it. Apparently the presence of mtMutS provides a selective advantage in octocorals, but one that is unnecessary for normal anthozoan mitochondrial function.
The subcellular target of the mtMutS mRNA, and the role of the gene product, remains unknown. The predicted protein structure, including the presence of a mismatch recognition site in Domain I (Figure 3), suggests a role in DNA mismatch repair (MMR) is highly likely. Our findings indicate a reduced rate of mtDNA evolution in at least one gene in the Octocorallia, relative to its sister group, that does not appear to extend to the nuclear genome (Figure 7). The eukaryotic MSH1 gene, for example, is a nuclear-encoded gene with a protein product that is translocated to the mitochondrion, where it performs MMR [27, 53, 54]. The MSH1 protein is known to form homodimers during mismatch binding in the mitochondrion , but its interaction with other proteins, particularly an endonuclease for strand excision, is unknown. The dimerization of mtMutS has not been examined, but it is possible that the presence of an HNH endonuclease domain could confer self-sufficiency to the protein . The domain is of a type II endonuclease, which is capable of nicking a single target strand such as the erroneous daughter strand in a mismatched DNA heteroduplex. Although prokaryotic MutS1 and eukaryotic MSH2, 3 and 6 MMR systems require one or more other proteins for excision (MutL+MutH and MLH, respectively) , the HNH domain on the C-terminal end of mtMutS could replace their actions, obviating the need for additional MMR components. This self-sufficiency could lend itself to mitochondrial MMR activity since MLH proteins associated with nuclear MMR display no active role in mitochondrial MMR [53, 56] and thus may be incapable of translocating across the organellar bilayer.
Implications and future directions
Regardless of the vector of mtMutS HGT, and of specific function, the gene's presence within the mitochondrial genome of octocorals is remarkable in itself. It adds another example of a metazoan subjected to HGT to a list that already includes nematodes  and molluscs . Unlike these other cases, however, mtMutS represents the only recorded instance of HGT into a metazoan mitochondrion. This necessitates a reconsideration of the mechanisms by which mitochondrial genomes evolve since it allows for the possibility of the introduction of foreign genetic material into a system previously thought to be constrained to inheritance through (typically maternal-) descent. It is also noteworthy that the evolution of mtMutS itself seems unconstrained by the reduced rate of mutation seen in other octocoral mitochondrial genes (Figure 7). We have now initiated studies to isolate nuclear MSH genes from octocorals and compare their ancestral patterns of diversification to mtMutS to determine if this pattern is shared by other MMR genes.
We have provided consolidating evidence for a non-eukaryotic origin of mtMutS and have demonstrated the monophyly of this gene among octocorals and other members of the MutS7 lineage. It now seems unequivocal that mtMutS arose once during octocorallian evolution, in the common ancestor of all extant lineages, through a horizontal gene transfer event from either an epsilonproteobacterium or a nucleocytoplasmic large DNA virus. This transfer represents the only instance recorded to date of HGT into a metazoan mitochondrion. The transferred gene appears to have been successfully integrated into the host genome, as evidenced by both its persistence over evolutionary time and its transcriptional activity. The mtMutS protein contains conserved amino acid positions and domains necessary for post-replicative DNA mismatch repair, as well as a unique additional domain that may enable it to function in a self-contained manner. Although mtMutS function remains unknown, we have demonstrated that at least one mitochondrial gene exhibits significantly lower rates of evolution in octocorals compared to hexacorals, providing evidence that mtMutS-mediated MMR may be operating within the octocoral mitochondrial genome.
MtMutS protein structure & gene tree prediction
All complete MutS derived amino acid sequences from octocorals (n = 30) were obtained from GenBank http://www.ncbi.nlm.nih.gov/. These sequences were then used in PSI-BLAST searches to obtain the most similar protein matches from non-octocorallian taxa. In addition, representatives of eubacterial MutS and the six eukaryotic MSH families were also downloaded for comparison. Eubacterial MutS sequences included Escherichia coli and bacterial hits to octocoral BLAST searches. All eukaryotic MSH sequences were obtained from Homo sapiens, Saccharomyces cerevisiae and Nematostella vectensis genomic datasets (Additional file 1).
The presence and identity of protein family domains ('pfam') within octocoral MutS sequences was inferred from queries in the Simple Modular Architecture Resource Tool (SMART) database of EMBL http://smart.embl-heidelberg.de/. The secondary and tertiary structures of protein sequences were inferred through submission to the I-TASSER server http://zhanglab.ccmb.med.umich.edu/I-TASSER/, which predicts structure and function through comparison to databases of known protein function. The results were compared to predicted structures of the most-similar non-octocorallian BLAST results as well as to the known structure of E. coli MutS, as determined by Lamers et al. . In particular, conserved functional residues in each domain were compared for their location within predicted domains and tertiary structures.
All amino acid sequences were aligned with the E-INS-I method in MAFFT , using the BLOSUM45 scoring matrix. The output was manually adjusted in MacClade 4.08 (Sinauer Assoc.) using the SMART pfam results to ensure homologous domains were correctly matched in the final alignment. Conserved functional residues or motifs of residues in each domain were also used to assist in the manual adjustment, which produced an alignment of 101 taxa and 1819 characters (available on request from the authors). Trimming the amino acid alignment of ambiguously aligned regions resulted in a truncated alignment of 592 characters, with most trimmed regions in Domain I and the interdomain region between Domain I and Domain III. Excluded regions typically contained gaps for all but a minority of taxa, and were thus unlikely to be informative in the phylogenetic analysis. Importantly, the C-terminal HNH domain was trimmed from the final alignment and excluded from subsequent phylogenetic analyses since a homologous domain was not found among the included eubacterial MutS and eukaryotic MSH sequences.
MEGA 5.05  was used to calculate maximum parsimony trees and to determine an initial distance correction model for Bayesian phylogenetic analysis. Tree searching used the Close-Neighbor-Interchange method with 1000 initial random trees and a MP search level of three. Clade support within the parsimony tree was calculated using 1000 replicates of nonparametric bootstrapping. The Bayesian analysis was implemented in MrBayes 3.1.2  using four independent MCMC runs, each containing six chains of 107 generations, with two swaps per generation ('nswaps = 2') and trees sampled every 1000 generations. The amino acid distance correction model was adjusted using the 'aamodelpr = mixed' prior in preliminary runs. Burn-in was determined using Tracer 1.5 http://beast.bio.ed.ac.uk/Tracer and the remaining trees were used to construct a 50% majority-rule consensus tree (the 'allcompat' option of 'sumt'). Gene trees obtained by maximum parsimony and Bayesian inference were compared for clade structure and for the determination of the sister taxon to the octocoral mtMutS lineage.
Expression of octocoral mtMutS
Gene expression of mtMutS in octocorals was tested through comparative amplification and sequencing of genomic DNA and reverse-transcribed mRNA. Fresh tissue samples of Isis hippuris, Nephthya sp. and Sarcophyton ehrenbergi were collected from Heron Island Reef, Australia (23.44710°S, 151.91800°E) and were immediately preserved in both 96% ethanol (for DNA extraction) and RNAlater (Sigma-Aldrich Ltd.) and frozen at -80°C (for RNA extraction). Total genomic DNA extraction used a standard 2XCTAB protocol  with a single chloroform extraction. Extraction of total RNA was performed with Trizol (Sigma-Aldrich Ltd.), following the manufacturer's recommended protocol for high-salt precipitation. Total RNA was then treated with 1 μl DNase I (Invitrogen Corp.) prior to reverse-transcription using SuperScriptIII (Invitrogen Corp.) to generate cDNA. Reverse transcriptase-free negative controls were included to check for genomic contamination of total RNA extracts.
Octocoral mtMutS in gDNA and cDNA was amplified using 0.22 mM of dNTPs, 3 units of Taq polymerase, 1X Taq Buffer (50 mM KCl, 10 mM Tris-HCl pH 9.0, 1.5 mM MgCl2), 2.9 mM of MgCl2, 0.16 μM of each primer and either 1 μl (cDNA) or 2 μl (gDNA) of template DNA. The forward primer (mtMutS93_Fwd: AGTTCTATGAACTTTGGCATGAG) was designed by the authors to anneal near the 5'-end (Domain I) of the octocoral mtMutS gene and was paired with an existing reverse primer (Mut3458R: TSGAGCAAAAGCCACTCC) . Although both primers were designed for specificity to mtMutS, BLAST searches were conducted on each to ensure nuclear MSH or other heterologous loci would not be amplified. Reactions were performed under thermocycler conditions of 94°C for 2 min followed by 35 cycles of 94°C for 60 s, 58°C for 60 s and 72°C for 60 s, with a final extension of 5 minutes at 72°C. Template-free negative PCR controls were included in all reactions to ensure observed products were of octocoral origin.
The size of amplified products from gDNA and cDNA were compared on a 1% agarose gel and RT-PCR controls were included to ensure cDNA amplicons were of mRNA origin. Bands of expected size were excised from the gel and dissolved in 4 volumes of 6 M NaI at 55°C for 5 minutes, followed by the addition of 15 μl hydrated silica suspension. The solution was incubated for a further 10 minutes, spun at maximum speed for 30 s and the supernatant was discarded. Pellets were resuspended in 500 μl of buffer (50 mM NaCl, 10 mM Tris-HCl pH 7.5, 2.5 mM EDTA, 50% ethanol), incubated at 55°C for 5 min, spun down and the supernatant was discarded again. The pellets were washed a second time and the pellets were dried prior to resuspension in 20 μl water.
Direct cycle sequencing of these gel-purified PCR products was performed in both forward and reverse directions in 10 μl volumes using 18-30 ng of template DNA, 8 nmol of primer, 1 μl of BigDye Terminator Ready Reaction Mix (Applied Biosystems) and 1.5 μl of buffer. Cycling was performed at 96°C for 2 minutes, followed by 30 cycles of 96°C for 10 s, 50°C for 5 s and 60°C for 4 minutes. Products were cleaned, precipitated and submitted for capillary sequencing at the Australian Genome Research Facility following their guidelines http://www.agrf.org.au/.
Sequences for each colony were assembled and edited in CodonCode Aligner 220.127.116.11 (CodonCode Corp.) and were submitted for BLASTn searches to ensure they were octocoral in origin and not due to contamination. Corresponding mtMutS sequences of gDNA and cDNA for each amplicon were then compared for identity and were deposited into GenBank as accession numbers JN383337-JN383342.
Comparison of Mitochondrial and Nuclear Rates of Evolution Among Anthozoans
All existing octocoral mtMutS and anthozoan cytochrome oxidase subunit I (COI) and small ribosomal subunit (18S) DNA sequences over 200 bp in length were obtained from GenBank. The sets of sequences for mtMutS and COI were aligned according to their amino acid translation frame in MacClade and all three datasets were trimmed to exclude intergenic regions, introns (present in hexacoral COI, including actiniarians [AF000023], scleractinians [FJ345452, FJ345449] and antipatharians [FJ597643]) and any ambiguous sections of alignment. Incomplete 5' and 3' ends were coded as missing data ("?") for sequences not spanning the length of the genes. Redundant sets of sequences were each collapsed into a single, multi-species taxon and the resulting matrix was used to calculate uncorrected pairwise distances among octocoral mtMutS sequences, octocoral COI sequences, hexacoral COI sequences, octocoral 18S sequences and hexacoral 18S sequences in MEGA. Mann-Whitney U-tests were used to test for significant differences in the mean pairwise differences between octocoral mtMutS and octocoral COI, between octocoral COI and hexacoral COI and between octocoral 18S and hexacoral 18S.
We acknowledge the University of Queensland's Heron Island Research Station for facilitating collection of octocorals from Heron Island Reef, and to the Great Barrier Reef Marine Park Authority for providing a permit for collection of octocoral material. We greatly appreciate the helpful discussion and insightful comments provided by members of the joint Degnan Lab, and by two anonymous reviewers. This work was made possible by generous funding support from the Deep Ocean Australia program and the Australian Research Council to SD.
- Wolstenholme DR: Animal Mitochondrial DNA: Structure and Evolution. International Review of Cytology. Edited by: Kwang WJ. 1992, David RW: Academic Press, 141: 173-216.Google Scholar
- Boore JL: Animal mitochondrial genomes. Nucleic Acids Res. 1999, 27: 1767-1780. 10.1093/nar/27.8.1767.View ArticlePubMedPubMed CentralGoogle Scholar
- Warrior R, Gall J: The mitochondrial DNA of Hydra attenuata and Hydra littoralis consists of two linear molecules. Archives des Science, Geneva. 1985, 38: 439-445.Google Scholar
- Bridge D, Cunningham CW, Schierwater B, Desalle R, Buss LW: Class-level relationships in the phylum Cnidaria: Evidence from mitochondrial genome structure. Proc Natl Acad Sci USA. 1992, 89: 8750-8753. 10.1073/pnas.89.18.8750.View ArticlePubMedPubMed CentralGoogle Scholar
- Beagley CT, Macfarlane JL, Pont-Kingdon GA, Okimoto R, Okada NA, Wolstenholme DR: Mitochondrial genomes of anthozoa (Cnidaria). Prog Cell Res. 1995, 5: 149-153.View ArticleGoogle Scholar
- Beagley CT, Okada NA, Wolstenholme DR: Two mitochondrial group I introns in a metazoan, the sea anemone Metridium senile: One intron contains genes for subunits 1 and 3 of NADH dehydrogenase. Proc Natl Acad Sci USA. 1996, 93: 5619-5623. 10.1073/pnas.93.11.5619.View ArticlePubMedPubMed CentralGoogle Scholar
- Beagley CT, Okimoto R, Wolstenholme DR: The mitochondrial genome of the sea anemone Metridium senile (Cnidaria): Introns, a paucity of tRNA genes, and a near-standard genetic code. Genetics. 1998, 148: 1091-1108.PubMedPubMed CentralGoogle Scholar
- Brugler MR, France SC: The complete mitochondrial genome of the black coral Chrysopathes formosa (Cnidaria: Anthozoa: Antipatharia) supports classification of antipatharians within the subclass Hexacorallia. Mol Phylogen Evol. 2007, 42: 776-788. 10.1016/j.ympev.2006.08.016.View ArticleGoogle Scholar
- Sinniger F, Pawlowski J: The partial mitochondrial genome of Leiopathes glaberrima (Hexacorallia: Antipatharia) and the first report of the presence of an intron in COI in black corals. Galaxea, Journal of Coral Reef Studies. 2009, 11: 21-26. 10.3755/galaxea.11.21.View ArticleGoogle Scholar
- Brugler MR, France SC: The mitochondrial genome of a deep-sea bamboo coral (Cnidaria, Anthozoa, Octocorallia, Isididae): genome structure and putative origins of replication are not conserved among octocorals. J Mol Evol. 2008, 67: 125-136. 10.1007/s00239-008-9116-2.View ArticlePubMedGoogle Scholar
- Uda K, Komeda Y, Koyama H, Koga K, Fujita T, Iwasaki N, Suzuki T: Complete mitochondrial genomes of two Japanese precious corals, Paracorallium japonicum and Corallium konojoi (Cnidaria, Octocorallia, Coralliidae): Notable differences in gene arrangement. Gene. 2011, 476: 27-37. 10.1016/j.gene.2011.01.019.View ArticlePubMedGoogle Scholar
- Lunt DH, Hyman BC: Animal mitochondrial DNA recombination. Nature. 1997, 387: 247-247. 10.1038/387247a0.View ArticlePubMedGoogle Scholar
- Beaton MJ, Roger AJ, Cavalier-Smith T: Sequence analysis of the mitochondrial genome of Sarcophyton glaucum: conserved gene order among octocorals. J Mol Evol. 1998, 47: 697-708. 10.1007/PL00006429.View ArticlePubMedGoogle Scholar
- Pont-Kingdon GA, Okada NA, Macfarlane JL, Beagley CT, Wolstenholme DR, Cavalier-Smith T, Clark-Walker GD: A coral mitochondrial mutS gene. Nature. 1995, 375: 109-111.View ArticlePubMedGoogle Scholar
- Pont-Kingdon GA, Okada NA, Macfarlane JL, Beagley CT, Watkins-Sims CD, Cavalier-Smith T, Clark-Walker GD, Wolstenholme DR: Mitochondrial DNA of the coral Sarcophyton glaucum contains a gene for a homologue of bacterial MutS: a possible case of gene transfer from the nucleus to the mitochondrion. J Mol Evol. 1998, 46: 419-431. 10.1007/PL00006321.View ArticlePubMedGoogle Scholar
- Sánchez JA, McFadden CS, France SC, Lasker HR: Molecular phylogenetic analyses of shallow-water Caribbean octocorals. Mar Biol. 2003, 142: 975-987.Google Scholar
- McFadden CS, France SC, Sánchez JA, Alderslade P: A molecular phylogenetic analysis of the Octocorallia (Cnidaria: Anthozoa) based on mitochondrial protein-coding sequences. Mol Phylogen Evol. 2006, 41: 513-527. 10.1016/j.ympev.2006.06.010.View ArticleGoogle Scholar
- Culligan KM, Meyer-Gauen G, Lyons-Weiler J, Hays JB: Evolutionary origin, diversification and specialization of eukaryotic MutS homolog mismatch repair proteins. Nucleic Acids Res. 2000, 28: 463-471. 10.1093/nar/28.2.463.View ArticlePubMedPubMed CentralGoogle Scholar
- Schofield MJ, Hsieh P: DNA mismatch repair: molecular mechanisms and biological function. Ann Rev Microbiol. 2003, 57: 579-608. 10.1146/annurev.micro.57.030502.090847.View ArticleGoogle Scholar
- McFadden CS, Benayahu Y, Pante E, Thoma JN, Nevarez PA, France SC: Limitations of mitochondrial gene barcoding in Octocorallia. Mol Ecol Resources. 2010, 11: 19-31.View ArticleGoogle Scholar
- Wirshing HH, Messing CG, Douady CJ, Reed J, Stanhope MJ, Shivji MS: Molecular evidence for multiple lineages in the gorgonian family Plexauridae (Anthozoa: Octocorallia). Mar Biol. 2005, 147: 497-508. 10.1007/s00227-005-1592-y.View ArticleGoogle Scholar
- McFadden CS, van Ofwegen LP, Beckman EJ, Benayahu Y, Alderslade P: Molecular systematics of the speciose Indo-Pacific soft coral genus, Sinularia (Anthozoa: Ocotocorallia). Invert Biol. 2009, 128: 303-323. 10.1111/j.1744-7410.2009.00179.x.View ArticleGoogle Scholar
- van Ofwegen LP, Groenenberg DSJ: A centuries old problem in nephtheid taxonomy approached using DNA data (Coelenterata: Alcyonacea). Contributions to Zoology. 2007, 76: 153-178.Google Scholar
- Bilewitch JP, Coates KA, Currie DC, Trapido-Rosenthal HG: Molecular and morphological variation supports monotypy of the octocoral Briareum Blainville, 1830 (Octocorallia: Alcyonacea) in the Western Atlantic. Proc Biol Soc Wash. 2010, 123: 93-112. 10.2988/09-22.1.View ArticleGoogle Scholar
- Herrera S, Baco A, Sánchez JA: Molecular systematics of the bubblegum coral genera (Paragorgiidae, Octocorallia) and description of a new deep-sea species. Mol Phylogen Evol. 2010, 55: 123-135. 10.1016/j.ympev.2009.12.007.View ArticleGoogle Scholar
- Sánchez JA, Cairns SD: An unusual new gorgonian coral (Anthozoa: Octocorallia) from the Aleutian Islands, Alaska. Zool Med. 2004, 78: 265-274.Google Scholar
- Reenan RAG, Kolodner RD: Characterization of insertion mutations in the Saccharomyces cerevisiae MSHl and MSH2 genes: evidence for separate mitochondrial and nuclear function. Genetics. 1992, 132: 975-985.PubMedPubMed CentralGoogle Scholar
- Sia EA, Kirkpatrick DT: The yeast MSH1 gene is not involved in DNA repair or recombination during meiosis. DNA Repair. 2005, 4: 253-261. 10.1016/j.dnarep.2004.10.002.View ArticlePubMedGoogle Scholar
- Adams KL, Palmer JD: Evolution of mitochondrial gene content: gene loss and transfer to the nucleus. Mol Phylogen Evol. 2003, 29: 380-395. 10.1016/S1055-7903(03)00194-5.View ArticleGoogle Scholar
- Eisen JA: A phylogenomic study of the MutS family of proteins. Nucleic Acids Res. 1998, 26: 4291-4300. 10.1093/nar/26.18.4291.View ArticlePubMedPubMed CentralGoogle Scholar
- Malik HS, Henikoff S: Dual recognition-incision enzymes might be involved in mismatch repair and meiosis. Trends Biochem Sci. 2000, 25: 414-418. 10.1016/S0968-0004(00)01623-6.View ArticlePubMedGoogle Scholar
- Claverie JM, Grzela R, Lartigue A, Bernadac A, Nitsche S, Vacelet J, Ogata H, Abergel C: Mimivirus and Mimiviridae: giant viruses with an increasing number of potential hosts, including corals and sponges. J Invert Path. 2009, 101: 172-180. 10.1016/j.jip.2009.03.011.View ArticleGoogle Scholar
- Ogata H, Ray J, Toyoda K, Sandaa RA, Nagasaki K, Bratbak G, Claverie JM: Two new subfamilies of DNA mismatch repair proteins (MutS) specifically abundant in the marine environment. Intl Soc Microbial Ecol J. 2011, 2011: 1-9.Google Scholar
- Blaxter M: Symbiont genes in host genomes: fragments with a future?. Cell Host & Microbe. 2007, 2: 211-213. 10.1016/j.chom.2007.09.008.View ArticleGoogle Scholar
- Mayer W, Schuster L, Bartelmes G, Dieterich C, Sommer R: Horizontal gene transfer of microbial cellulases into nematode genomes is associated with functional assimilation and gene turnover. BMC Evol Biol. 2011, 11: 1-10. 10.1186/1471-2148-11-1.View ArticleGoogle Scholar
- Lamers MH, Perrakis A, Enzlin JH, Winterwerp HHK, de Wind N, Sixma TK: The crystal structure of DNA mismatch repair protein MutS binding to a G-T mismatch. Nature. 2000, 407: 711-717. 10.1038/35037523.View ArticlePubMedGoogle Scholar
- Obmolova G, Ban C, Hsieh P, Yang W: Crystal structures of mismatch repair protein MutS and its complex with a substrate DNA. Nature. 2000, 407: 705-710.Google Scholar
- Junop MS, Obmolova G, Rausch K, Hsieh P, Yang W: Composite active site of an ABC ATPase: MutS uses ATP to verify mismatch recognition and authorize DNA repair. Mol Cell. 2001, 7: 1-12. 10.1016/S1097-2765(01)00149-6.View ArticlePubMedGoogle Scholar
- Romano SL, Palumbi SR: Evolution of scleractinian corals inferred from molecular systematics. Science. 1996, 271: 640-642. 10.1126/science.271.5249.640.View ArticleGoogle Scholar
- France SC, Hoover LL: Analysis of variation in mitochondrial DNA sequences (ND3, ND4L, MSH) among Octocorallia (= Alcyonaria) (Cnidaria: Anthozoa). Bull Biol Soc Wash. 2001, 10: 110-118.Google Scholar
- Ham JLvd, Brugler MR, France SC: Exploring the utility of an indel-rich, mitochondrial intergenic region as a molecular barcode for bamboo corals (Octocorallia: Isididae). Marine Genomics. 2009, 2: 183-192. 10.1016/j.margen.2009.10.002.View ArticlePubMedGoogle Scholar
- Rogers HJ, Buck KW, Brasier CM: A mitochondrial target for double-stranded RNA in diseased isolates of the fungus that causes Dutch elm disease. Nature. 1987, 329: 558-560. 10.1038/329558a0.View ArticleGoogle Scholar
- Rojo G, Chamorro M, Salas ML, Vinuela E, Cuezva JM, Salas J: Migration of mitochondria to viral assembly sites in African swine fever virus-infected cells. J Virol. 1998, 72: 7583-7588.PubMedPubMed CentralGoogle Scholar
- Murata T, Goshima F, Daikoku T, Inagaki-Ohara K, Takakuwa H, Kato K, Nishiyama Y: Mitochondrial distribution and function in herpes simplex virus-infected cells. J Gen Virol. 2000, 81: 401-406.View ArticlePubMedGoogle Scholar
- Rohwer F, Seguritan V, Azam F, Knowlton N: Diversity and distribution of coral-associated bacteria. Mar Ecol Prog Ser. 2002, 243: 1-10.View ArticleGoogle Scholar
- Penn K, Wu DY, Eisen JA, Ward N: Characterization of bacterial communities associated with deep-sea corals on Gulf of Alaska seamounts. Appl Environ Microb. 2006, 72: 1680-1683. 10.1128/AEM.72.2.1680-1683.2006.View ArticleGoogle Scholar
- Ainsworth TD, Hoegh-Guldberg O: Bacterial communities closely associated with coral tissues vary under experimental and natural reef conditions and thermal stress. Aquat Biol. 2009, 4: 289-296.View ArticleGoogle Scholar
- Lesser MP, Falcon LI, Rodriguez-Roman A, Enriquez S, Hoegh-Guldberg O, Iglesias-Prieto R: Nitrogen fixation by symbiotic cyanobacteria provides a source of nitrogen for the scleractinian coral Montastraea cavernosa. Mar Ecol Prog Ser. 2007, 346: 143-152.View ArticleGoogle Scholar
- McDaniel LD, Young E, Delaney J, Ruhnau F, Ritchie KB, Paul JH: High frequency of horizontal gene transfer in the oceans. Science. 2010, 330: 50-50. 10.1126/science.1192243.View ArticlePubMedGoogle Scholar
- Li L, Kato C, Horikoshi K: Microbial diversity in sediments collected from the deepest cold-seep area, the Japan Trench. Mar Biotech. 1999, 1: 391-400. 10.1007/PL00011793.View ArticleGoogle Scholar
- Campbell BJ, Smith JL, Hanson TE, Klotz MG, Stein LY, Lee CK, Wu DY, Robinson JM, Khouri HM, Eisen JA, Cary SC: Adaptations to submarine hydrothermal environments exemplified by the genome of Nautilia profundicola. PLoS Genet. 2009, 5: 1-19.View ArticleGoogle Scholar
- Watling L, Auster PJ: Distribution of deep-water Alcyonacea off the northeast coast of the United States. Cold-water Corals and Ecosystems. Edited by: Freiwald A, Roberts JM. 2005, Berlin: Springer-Verlag, 279-296.View ArticleGoogle Scholar
- Chi NW, Kolodner RD: Purification and characterization of MSH1, a yeast mitochondrial protein that binds to DNA mismatches. J Biol Chem. 1994, 269: 29984-29992.PubMedGoogle Scholar
- Reenan RA, Kolodner RD: Isolation and characterization of two Saccharomyces cerevisiae genes encoding homologs of the bacterial HexA and MutS mismatch repair proteins. Genetics. 1992, 132: 963-973.PubMedPubMed CentralGoogle Scholar
- Mookerjee SA, Sia EA: Overlapping contributions of Msh1p and putative recombination proteins Cce1p, Din7p, and Mhr1p in large-scale recombination and genome sorting events in the mitochondrial genome of Saccharomyces cerevisiae. Mutation Res. 2006, 595: 91-106.View ArticlePubMedGoogle Scholar
- Harfe BD, Jinks-Robertson S: DNA mismatch repair and genetic instability. Ann Rev Genet. 2000, 34: 359-399. 10.1146/annurev.genet.34.1.359.View ArticlePubMedGoogle Scholar
- Rumpho ME, Worful JM, Lee J, Kannan K, Tyler MS, Bhattacharya D, Moustafa A, Manhart JR: Horizontal gene transfer of the algal nuclear gene psbO to the photosynthetic sea slug Elysia chlorotica. Proc Natl Acad Sci USA. 2008, 105: 17867-17871. 10.1073/pnas.0804968105.View ArticlePubMedPubMed CentralGoogle Scholar
- Katoh K, Misawa K, Kuma K, Miyata T: MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002, 30: 3059-3066. 10.1093/nar/gkf436.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.View ArticlePubMedGoogle Scholar
- Berntson EA, Bayer FM, McArthur AG, France SC: Phylogenetic relationships within the Octocorallia (Cnidaria: Anthozoa) based on nuclear 18S rRNA sequences. Mar Biol. 2001, 138: 235-246. 10.1007/s002270000457.View ArticleGoogle Scholar
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.