Complete sequence of the amphioxus (Branchiostoma lanceolatum) mitochondrial genome: relations to vertebrates

Molecular taxonomy confirms that the northeastern Atlantic and Mediterranean Sea harbor a single lancelet, Branchiostoma lanceolatum (Pallas, 1774) (Cephalochordata: Leptocardii: Branchiostomatidae)

Branchiostomatidae (lancelets or amphioxus) comprises about 30 species, several of which are well-established models in evolutionary development. Our zoological and ecological knowledge of the family is nonetheless limited. Despite evident differences can be found among known populations, the taxonomy of Branchiostoma lanceolatum (type species of the genus Branchiostoma) has never been investigated with modern methods through its range in the northeastern Atlantic and Mediterranean Sea. We address this via a multilocus molecular approach and comparing specimens collected from different European populations. Results obtained here confirm the presence of a single species inhabiting the range between the topotypical localities of B. lanceolatum (Atlantic Ocean) and of its junior synonym B. lubricum (Mediterranean Sea), without evincing geographical structure between populations. This suggests that environment most likely drives the characteristics observed in different geographic areas. The long larval phase and the slow mutation rate in lancelets may have played a key role in the evolutionary history of this iconic species.

Whole-Genome Resequencing of Twenty Branchiostoma belcheri Individuals Provides a Brand-New Variant Dataset for Branchiostoma

As the extant representatives of the basal chordate lineage, amphioxi (including the genera Branchiostoma, Asymmetron and Epigonichthys) play important roles in tracing the state of chordate ancestry. Previous studies have reported that members of the Branchiostoma species have similar morphological phenotypic characteristics, but in contrast, there are high levels of genetic polymorphisms in the populations. Here, we resequenced 20 Branchiostomabelcheri genomes to an average depth of approximately 12.5X using the Illumina HiSeq 2000 platform. In this study, over 52 million variations (~12% of the total genome) were detected in the B. belcheri population, and an average of 12.8 million variations (~3% of the total genome) were detected in each individual, confirming that Branchiostoma is one of the most genetically diverse species sequenced to date. Demographic inference analysis highlighted the role of historical global temperature in the long-term population dynamics of Branchiostoma, and revealed a population expansion at the Greenlandian stage of the current geological epoch. We detected 594 Single nucleotide polymorphism and 148 Indels in the Branchiostoma mitochondrial genome, and further analyzed their genetic mutations. A recent study found that the epithelial cells of the digestive tract in Branchiostoma can directly phagocytize food particles and convert them into absorbable nontoxic nutrients using powerful digestive and immune gene groups. In this study, we predicted all potential mutations in intracellular digestion-associated genes. The results showed that most “probably damaging” mutations were related to rare variants (MAF<0.05) involved in strengthening or weakening the intracellular digestive capacity of Branchiostoma. Due to the extremely high number of polymorphisms in the Branchiostoma genome, our analysis with a depth of approximately 12.5X can only be considered a preliminary analysis. However, the novel variant dataset provided here is a valuable resource for further investigation of phagocytic intracellular digestion in Branchiostoma and determination of the phenotypic and genotypic features of Branchiostoma.

Higher Order Organization of the mtDNA: Beyond Mitochondrial Transcription Factor A

The higher order organization of eukaryotic and prokaryotic genomes is pivotal in the regulation of gene expression. Specifically, chromatin accessibility in eukaryotes and nucleoid accessibility in bacteria are regulated by a cohort of proteins to alter gene expression in response to diverse physiological conditions. By contrast, prior studies have suggested that the mitochondrial genome (mtDNA) is coated solely by mitochondrial transcription factor A (TFAM), whose increased cellular concentration was proposed to be the major determinant of mtDNA packaging in the mitochondrial nucleoid. Nevertheless, recent analysis of DNase-seq and ATAC-seq experiments from multiple human and mouse samples suggest gradual increase in mtDNA occupancy during the course of embryonic development to generate a conserved footprinting pattern which correlate with sites that have low TFAM occupancy in vivo (ChIP-seq) and tend to adopt G-quadruplex structures. These findings, along with recent identification of mtDNA binding by known modulators of chromatin accessibility such as MOF, suggest that mtDNA higher order organization is generated by cross talk with the nuclear regulatory system, may have a role in mtDNA regulation, and is more complex than once thought.

Widespread use of the "ascidian" mitochondrial genetic code in tunicates

Background: Ascidians, a tunicate class, use a mitochondrial genetic code that is distinct from vertebrates and other invertebrates. Though it has been used to translate the coding sequences from other tunicate species on a case-by-case basis, it is has not been investigated whether this can be done systematically. This is an important because a) some tunicate mitochondrial sequences are currently translated with the invertebrate code by repositories such as NCBI GenBank, and b) uncertainties about the genetic code to use can complicate or introduce errors in phylogenetic studies based on translated mitochondrial protein sequences.

Methods: We collected publicly available nucleotide sequences for non-ascidian tunicates including appendicularians such as Oikopleura dioica, translated them using the ascidian mitochondrial code, and built multiple sequence alignments covering all tunicate classes.

Results: All tunicates studied here appear to translate AGR codons to glycine instead of serine (invertebrates) or as a stop codon (vertebrates), as initially described in ascidians. Among Oikopleuridae, we suggest further possible changes in the use of the ATA (Ile → Met) and TGA (Trp → Arg) codons.

Conclusions: We recommend using the ascidian mitochondrial code in automatic translation pipelines of mitochondrial sequences for all tunicates. Further investigation is required for additional species-specific differences.

Mitochondrial genomes and Doubly Uniparental Inheritance: new insights from Musculista senhousia sex-linked mitochondrial DNAs (Bivalvia Mytilidae)

BACKGROUND

Doubly Uniparental Inheritance (DUI) is a fascinating exception to matrilinear inheritance of mitochondrial DNA (mtDNA). Species with DUI are characterized by two distinct mtDNAs that are inherited either through females (F-mtDNA) or through males (M-mtDNA). DUI sex-linked mitochondrial genomes share several unusual features, such as additional protein coding genes and unusual gene duplications/structures, which have been related to the functionality of DUI. Recently, new evidence for DUI was found in the mytilid bivalve Musculista senhousia. This paper describes the complete sex-linked mitochondrial genomes of this species.

RESULTS

Our analysis highlights that both M and F mtDNAs share roughly the same gene content and order, but with some remarkable differences. The Musculista sex-linked mtDNAs have differently organized putative control regions (CR), which include repeats and palindromic motifs, thought to provide sites for DNA-binding proteins involved in the transcriptional machinery. Moreover, in male mtDNA, two cox2 genes were found, one (M-cox2b) 123bp longer.

CONCLUSIONS

The complete mtDNA genome characterization of DUI bivalves is the first step to unravel the complex genetic signals allowing Doubly Uniparental Inheritance, and the evolutionary implications of such an unusual transmission route in mitochondrial genome evolution in Bivalvia. The observed redundancy of the palindromic motifs in Musculista M-mtDNA may have a role on the process by which sperm mtDNA becomes dominant or exclusive of the male germline of DUI species. Moreover, the duplicated M-COX2b gene may have a different, still unknown, function related to DUI, in accordance to what has been already proposed for other DUI species in which a similar cox2 extension has been hypothesized to be a tag for male mitochondria.

The complete mitochondrial genome of the cyclopoid copepod Paracyclopina nana: a highly divergent genome with novel gene order and atypical gene numbers

In this paper, we describe the complete mitogenome of the cyclopoid copepod Paracyclopina nana with emphasis on the highly rearranged gene order and high divergence against published copepod mitogenomes. The P. nana mtDNA is 15,981 bp in length (70.9% AT) and consists of 37 genes (12 protein-coding genes, 2 rRNAs, 23 tRNAs) that are atypical for metazoan mitogenomes. Unusually, it contains an extra tRNA (tRNA-Ala) but it does not contain the ATPase 8 gene. The P. nana mitogenome has a long putative control region with high AT content (1351 bp, 77.0% AT). The Cyt b was considerably short in length, compared to other crustaceans. Compared to typical mitogenomes of arthropods and copepods, the gene order of the P. nana mitogenome is highly rearranged with a novel gene structure. In addition, P. nana has highly divergent mt genes (mostly less than 50%), judged by amino acid substitution. We present the first complete mitogenome sequence from a cyclopoid copepod, thereby increasing our understanding of copepod and crustacean evolution from the mitochondrial point of view.

Evolution of the mitochondrial genome of Metazoa as exemplified by comparison of congeneric species

The mitochondrial genome (mtDNA) of Metazoa is a good model system for evolutionary genomic studies and the availability of more than 1000 sequences provides an almost unique opportunity to decode the mechanisms of genome evolution over a large phylogenetic range. In this paper, we review several structural features of the metazoan mtDNA, such as gene content, genome size, genome architecture and the new parameter of gene strand asymmetry in a phylogenetic framework. The data reviewed here show that: (1) the plasticity of Metazoa mtDNA is higher than previously thought and mainly due to variation in number and location of tRNA genes; (2) an exceptional trend towards stabilization of genomic features occurred in deuterostomes and was exacerbated in vertebrates, where gene content, genome architecture and gene strand asymmetry are almost invariant. Only tunicates exhibit a very high degree of genome variability comparable to that found outside deuterostomes. In order to analyse the genomic evolutionary process at short evolutionary distances, we have also compared mtDNAs of species belonging to the same genus: the variability observed in congeneric species significantly recapitulates the evolutionary dynamics observed at higher taxonomic ranks, especially for taxa showing high levels of genome plasticity and/or fast nucleotide substitution rates. Thus, the analysis of congeneric species promises to be a valuable approach for the assessment of the mtDNA evolutionary trend in poorly or not yet sampled metazoan groups.

The mechanisms of codon reassignments in mitochondrial genetic codes

Many cases of nonstandard genetic codes are known in mitochondrial genomes. We carry out analysis of phylogeny and codon usage of organisms for which the complete mitochondrial genome is available, and we determine the most likely mechanism for codon reassignment in each case. Reassignment events can be classified according to the gain-loss framework. The “gain” represents the appearance of a new tRNA for the reassigned codon or the change of an existing tRNA such that it gains the ability to pair with the codon. The “loss” represents the deletion of a tRNA or the change in a tRNA so that it no longer translates the codon. One possible mechanism is codon disappearance (CD), where the codon disappears from the genome prior to the gain and loss events. In the alternative mechanisms the codon does not disappear. In the unassigned codon mechanism, the loss occurs first, whereas in the ambiguous intermediate mechanism, the gain occurs first. Codon usage analysis gives clear evidence of cases where the codon disappeared at the point of the reassignment and also cases where it did not disappear. CD is the probable explanation for stop to sense reassignments and a small number of reassignments of sense codons. However, the majority of sense-to-sense reassignments cannot be explained by CD. In the latter cases, by analysis of the presence or absence of tRNAs in the genome and of the changes in tRNA sequences, it is sometimes possible to distinguish between the unassigned codon and the ambiguous intermediate mechanisms. We emphasize that not all reassignments follow the same scenario and that it is necessary to consider the details of each case carefully.

A comparative genomics analysis of codon reassignments reveals a link with mitochondrial proteome size and a mechanism of genetic code change via suppressor tRNAs

Using a comparative genomics approach we demonstrate a negative correlation between the number of codon reassignments undergone by 222 mitochondrial genomes and the mitochondrial genome size, the number of mitochondrial ORFs, and the sizes of the large and small subunit mitochondrial rRNAs. In addition, we show that the TGA-to-tryptophan codon reassignment, which has occurred 11 times in mitochondrial genomes, is found in mitochondrial genomes smaller than those which have not undergone the reassignment. We therefore propose that mitochondrial codon reassignments occur in a wide range of phyla, particularly in Metazoa, due to a reduced "proteomic constraint" on the mitochondrial genetic code, compared to the nuclear genetic code. The reduced proteomic constraint reflects the small size of the mitochondrial-encoded proteome and allows codon reassignments to occur with less likelihood of lethality. In addition, we demonstrate a striking link between nonsense codon reassignments and the decoding properties of naturally occurring nonsense suppressor tRNAs. This suggests that natural preexisting nonsense suppression facilitated nonsense codon reassignments and constitutes a novel mechanism of genetic code change. These findings explain for the first time the identity of the stop codons and amino acids reassigned in mitochondrial and nuclear genomes. Nonsense suppressor tRNAs provided the raw material for nonsense codon reassignments, implying that the properties of the tRNA anticodon have dictated the identity of nonsense codon reassignments.

The complete mitochondrial genomes of the sea lily Gymnocrinus richeri and the feather star Phanogenia gracilis: Signature nucleotide bias and unique nad4L gene rearrangement within crinoids

Complete DNA sequences have been determined for the mitochondrial genomes of the crinoids Phanogenia gracilis (15892 bp) and Gymnocrinus richeri (15966 bp). The mitochondrial genetic map of the stalkless feather star P. gracilis is identical to that of the comatulid feather star Florometra serratissima (Scouras, A., Smith, M.J., 2001. Mol. Biol. Evol. 18, 61-73). The mitochondrial gene order of the stalked crinoid G. richeri differs from that of F. serratissima and P. gracilis by the transposition of the nad4L protein gene. The G. richeri nad4L mitochondrial map position is unique among metazoa and is likely a derived feature in this stalked crinoid. Nucleotide compositional analyses of protein genes encoded on the major sense strand confirm earlier conclusions regarding a crinoid-distinctive T over C bias. All three crinoids exhibit high T levels in third codon positions, whereas other echinoderm classes favor A or C in the third codon position. The nucleotide bias is reflected in the relative synonymous codon usage patterns of crinoids versus other echinoderms. We suggest that the nucleotide bias of crinoids, in comparison to other echinoderms, indicates that a physical inversion of the origin of replication has occurred in the crinoid lineage. Evolutionary rate tests support the use of the cytochrome b (cob) gene in molecular phylogenetic analyses of echinoderms. A consensus echinoderm tree was generated based on cytochrome b nucleotide alignments that placed the asteroids as a sister group to a clade containing the ophiuroids and the (echinoids+holothuroids) with the crinoids basal to the rest of the echinoderm classes: [Crinoid,(Asteroid,(Ophiuroid,(Echinoid,Holothuroid)))].

The complete mitochondrial genome sequence and characterization of single-nucleotide polymorphisms in the control region of the Asian seabass (Lates calcarifer)

We determined the complete mtDNA nucleotide sequence of Lates calcarifer using the shotgun sequencing method. The mitochondrial DNA (mtDNA) was 16,535 base pairs (bp) in length, and contained 13 protein coding genes, 22 transfer RNAs, 2 ribosomal RNAs, and one major noncoding control region (CR). The CR was unusually short at only 768 bp. A striking feature of the mitochondrial genome was the high G+C content (46.1%), which is among the highest in fish. The gene order was identical to that of a typical vertebrate. Phylogenetic analyses using concatenated amino acid sequences of 12 protein-coding genes of 30 fish species representing 14 suborders clearly showed Lates calcarifer was located in the cluster of fish species from the order Perciformes, supporting the traditional systematic classification. We characterized single-nucleotide polymorphisms (SNPs) in the CR by sequencing the complete CR of 25 individuals obtained from Australia and Singapore. A total of 68 SNPs were detected. Eighteen SNPs were fixed with alternative nucleotides in Australian and Singapore seabass, and these SNPs could be used for differentiating fish from the two countries.

Complete DNA sequence of the mitochondrial genome of the sea-slug, Aplysia californica: conservation of the gene order in Euthyneura

We have sequenced and characterized the complete mitochondrial genome of the sea slug, Aplysia californica, an important model organism in experimental biology and a representative of Anaspidea (Opisthobranchia, Gastropoda). The mitochondrial genome of Aplysia is in the small end of the observed sizes of animal mitochondrial genomes (14,117 bp, NCBI Accession No. NC_005827). The Aplysia genome, like most other mitochondrial genomes, encodes genes for 2 ribosomal subunit RNAs (small and large rRNAs), 22 tRNAs, and 13 protein subunits (cytochrome c oxidase subunits 1-3, cytochrome b apoenzyme, ATP synthase subunits 6 and 8, and NADH dehydrogenase subunits 1-6 and 4L). The gene order is virtually identical between opisthobranchs and pulmonates, with the majority of differences arising from tRNA translocations. In contrast, the gene order from representatives of basal gastropods and other molluscan classes is significantly different from opisthobranchs and pulmonates. The Aplysia genome was compared to all other published molluscan mitochondrial genomes and phylogenetic analyses were carried out using a concatenated protein alignment. Phylogenetic analyses using maximum likelihood based analyses of the well aligned regions of the protein sequences support both monophyly of Euthyneura (a group including both the pulmonates and opisthobranchs) and Opisthobranchia (as a more derived group). The Aplysia mitochondrial genome sequenced here will serve as an important platform in both comparative and neurobiological studies using this model organism.

Identification and characterisation of a homolog of an activation gene for the recombination activating gene 1 (RAG 1) in amphioxus

Expression of recombination activating genes (RAG) involved in the V (D) J recombination is regulated by the RAG1 gene activator (RGA) in mammals. The sequence of a cDNA clone from an amphioxus cDNA library was found to be homologous to that of RGA from mouse stromal cells with 45% identity. The full-length cDNA sequence comprises 1119 bp and encodes a putative protein of 210 amino acid residues. Characterisation of the amino acid sequence revealed that two MtN3 domains and seven transmembrane spans are present in this protein, indicating a potential role as a plasma membrane protein. This gene is expressed in many tissues and at differential developmental stages. A high expression level of RGA is detected in gonad tissues, and gastrula embryo and adult stages. The presence of the RGA gene in amphioxus suggests that the signal pathway required for the expression of RAG could exist in this primitive protochordate. It also implies that in the related molecules, primitive adaptive immunity may have existed in cephalochordate although the complete machinery of VDJ rearrangement may not be formed.

New Complete Mitochondrial DNA Sequence of the Lancelet Branchiostoma lanceolatum (Cephalochordata) and the Identity of this Species' Sequences

Three mitochondrial (mt) genes were sequenced for two Atlantic lancelet species, Branchiostoma lanceolatum and B. floridae, to examine a serious discrepancy among previously published results of molecular studies: substantial sequence difference in a nuclear gene vs. virtual identity in the mt genome sequence. The results revealed that three mt genes of B. lanceolatum, collected from Helgoland in the North Sea and Naples in the Mediterranean, were quite diverged from those of B. floridae, collected from Tampa Bay, Florida. Therefore, the previously recognized identity in the mt genome between the two species is attributable to misidentification of materials used. To correct this misleading information, the complete mtDNA sequence of B. lanceolatum was determined for an individual from Helgoland.

Evidence that the large noncoding sequence is the main control region of maternally and paternally transmitted mitochondrial genomes of the marine mussel (Mytilus spp.)

Both the maternal (F-type) and paternal (M-type) mitochondrial genomes of the Mytilus species complex M. edulis/galloprovincialis contain a noncoding sequence between the l-rRNA and the tRNA(Tyr) genes, here called the large unassigned region (LUR). The LUR, which is shorter in M genomes, is capable of forming secondary structures and contains motifs of significant sequence similarity with elements known to have specific functions in the sea urchin and the mammalian control region. Such features are not present in other noncoding regions of the F or M Mytilus mtDNA. The LUR can be divided on the basis of indels and nucleotide variation in three domains, which is reminiscent of the tripartite structure of the mammalian control region. These features suggest that the LUR is the main control region of the Mytilus mitochondrial genome. The middle domain has diverged by only 1.5% between F and M genomes, while the average divergence over the whole molecule is approximately 20%. In contrast, the first domain is among the most divergent parts of the genome. This suggests that different parts of the LUR are under different selection constraints that are also different from those acting on the coding parts of the molecule.

Evolution of the Mitochondrial Genome in Cephalochordata as Inferred from Complete Nucleotide Sequences from Two Epigonichthys Species

Complete mitochondrial (mt) DNA sequences of two lancelets, Epigonichthys maldivensis and E. lucayanus, were compared with those of two Branchiostoma lancelets and several deuterostomes previously surveyed. The mt-gene order of E. lucayanus was quite different from that of E. maldivensis, the latter being identical to the two Branchiostoma species. A remarkable genomic change in E. lucayanus mtDNA was an inversion, indicating the possibility of recombination of the mt-genome. Gene rearrangements, probably attributable to tandem genome duplications and subsequent random deletions, were observed in two parts. Short major unassignable sequences of the examined lancelets were regarded as a part of putative regulative elements, judging from some sequence similarity to the conserved sequence block (CSB) in mammalian mtDNA. The considerable mt-genome reorganization in E. lucayanus seemed to have affected the nucleotide substitution pattern, suggested by base composition analyses. The present analysis also suggested that AGR codons in lancelet mtDNA were likely to correspond to serine residue, rather than glycine. Furthermore, the AGG codon, so far reputed to be unassignable in lancelet mtDNA, was found twice in E. maldivensis, indicating the availability of all four AGN codons in some lancelets. This finding lends support to an alternative hypothesis regarding the evolutionary history of AGR-codon assignment in extant chordates, rather than that previously proposed. A molecular phylogenetic tree of the Epigonichthys and Branchiostoma species based on DNA sequences of the 13 mt-protein genes doubted the monophyly of the former genus, unlike the prevailing classification based on their different gonadal arrangements.

The Complete Maternal and Paternal Mitochondrial Genomes of the Mediterranean Mussel Mytilus galloprovincialis: Implications for the Doubly Uniparental Inheritance Mode of mtDNA

The maternal (F) and paternal (M) mitochondrial genomes of the mussel Mytilus galloprovincialis have diverged by about 20% in nucleotide sequence but retained identical gene content and gene arrangement and similar nucleotide composition and codon usage bias. Both lack the ATPase8 subunit gene, have two tRNAs for methionine and a longer open-reading frame for cox3 than seen in other mollusks. Between the F and M genomes, tRNAs are most conserved followed by rRNAs and protein-coding genes, even though the degree of divergence varies considerably among the latter. Divergence at nad3 is exceptionally low most likely because this gene includes the origin of transcription of the lagging strand (O(L)). Noncoding regions are the least conserved with the notable exception of the central domain of the main control region and a segment of another noncoding region immediately following nad3. The amino acid divergence (14%) of the two genomes is smaller than in two other pairs of conspecific genomes that are available in GenBank, that of the clam Venerupis philippinarum (34%) and of the fresh water mussel Inversidens japanensis (50%), suggesting that doubly uniparental inheritance of mtDNA emerged at different times in the three species or that there has been a relatively recent replacement of the male genome by the female in the Mytilus line. The latter hypothesis is supported from phylogenetic and population studies of Mytilidae. That the M genome contains a full complement of genes with no premature termination codons argues against it being a selfish element that rides with the sperm. It is shorter than the F by 118 bp, which apparently cannot account for the postulated replicative advantage of this genome over the F in male gonads. The high similarity of the two genomes explains why the F genome may assume the role of the M genome, but it does not exclude the possibility that for this to happen some M-specific sequences must be transferred on to the F genome by means of recombination. If such sequences exist they would most likely be located in noncoding regions.

Complete mitochondrial genome DNA sequence for two ophiuroids and a holothuroid: the utility of protein gene sequence and gene maps in the analyses of deep deuterostome phylogeny

The complete mitochondrial genome sequences have been determined for the holothuroid Cucumaria miniata and two ophiuroid species Ophiopholis aculeata and Ophiura lütkeni. In addition, the nucleotide sequence of the mitochondrial protein-coding genes for the asteroid Pisaster ochraceus has been completed. Maximum-likelihood and LogDet distance analyses of concatenated protein-coding sequences produced a series of trees that did not conclusively support generally accepted models of echinoderm phylogeny. The ophiuroid data consistently demonstrated accelerated nucleotide divergence rates and lack of stationarity. This confounds the phylogenetic analyses. Molecular investigations using individual protein-coding gene alignments demonstrated that the cytochrome b gene exhibits the least deviation in rate and stationarity and generated some trees consistent with proposed echinoderm phylogenies. Phylogenies based on echinoderm mitochondrial gene rearrangements also proved problematic because of extensive variation in gene order between and within classes. A comparison of the two distinctive ophiuroid mitochondrial gene orders supports the hypothesis that O. lütkeni has a more derived mitochondrial gene order versus O. aculeata. The variation in the echinoderm mitochondrial gene maps reinforces the limitations of the application of mitochondrial gene rearrangements as a global phylogenetic tool.

Complete mtDNA of Ciona intestinalis reveals extensive gene rearrangement and the presence of an atp8 and an extra trnM gene in ascidians

The complete mitochondrial genome (mtDNA) of the model organism Ciona intestinalis (Urochordata, Ascidiacea) has been amplified by long-PCR using specific primers designed on putative mitochondrial transcripts identified from publicly available mitochondrial-like expressed sequence tags. The C. intestinalis mtDNA encodes 39 genes: 2 rRNAs, 13 subunits of the respiratory complexes, including ATPase subunit 8 ( atp8), and 24 tRNAs, including 2 tRNA-Met with anticodons 5'-UAU-3'and 5'-CAU-3', respectively. All genes are transcribed from the same strand. This gene content seems to be a common feature of ascidian mtDNAs, as we have verified the presence of a previously undetected atp8 and of two trnM genes in the two other sequenced ascidian mtDNAs. Extensive gene rearrangement has been found in C. intestinalis with respect not only to the common Vertebrata/Cephalochordata/Hemichordata gene organization but also to other ascidian mtDNAs, including the cogeneric Ciona savignyi. Other features such as the absence of long noncoding regions, the shortness of rRNA genes, the low GC content (21.4%), and the absence of asymmetric base distribution between the two strands suggest that this genome is more similar to those of some protostomes than to deuterostomes.

Ancient phylogenetic separation between Pacific and Atlantic cephalochordates as revealed by mitochondrial genome analysis

The subphylum Cephalochordata (lancelets) is a relatively small taxonomic group in contrast to the subphyla Urochordata and Vertebrata. As an initial step to determine whether lancelets exhibit small genetic divergence in keeping with their conservative body organization or large genetic variation, four Branchiostoma species from the Pacific (B. belcheri and B. malayanum) and Atlantic (B. floridae and B. lanceolatum) Oceans were genetically compared using partial mitochondrial DNA sequences of the cytochrome oxidase c subunit I (COI) and 16S ribosomal RNA (16S rRNA) genes. In both genes, large genetic differences were revealed between the Pacific and Atlantic species, as well as within the former. Two maximum-likelihood trees from the COI and 16S rRNA genes showed that the Pacific and Atlantic lancelets were reciprocally clustered into different clades. Furthermore, both gene trees consistently exhibited deep phylogenetic separation between the two oceans. The estimated divergence time suggested that differentiation may have followed the migration of ancestral lancelets from the Pacific to the Atlantic Oceans via the Tethys Sea.

Phylogeny of hagfish based on the mitochondrial 16S rRNA gene

The phylogenetic relationships among the species belonging to the family Myxinidae are still debatable. The mitochondrial DNA sequences from the large ribosomal RNA gene may be of great value for systematic and phylogenetic studies within families. Partial sequences of the 16S rRNA gene were obtained for comparisons among the following hagfish species, Paramyxine nelsoni, Paramyxine sheni, Paramyxine taiwanae, Paramyxine yangi, Paramyxine cheni, Eptatretus burgeri, Eptatretus stouii, Eptatretus cirrhatus, Myxine glutinosa, Myxine formosana, Myxine circifrons, Myxine sp1, and Myxine sp2. The boundary of four Paramyxine species (P. sheni, P. taiwanae, P. nelsoni, and P. yangi) from 16S rRNA sequences is ambiguous, however, they are valid based on our unpublished isozyme data as well as the gill aperture arrangement pattern. Both NJ and MP trees constructed from the present molecular data indicate that the genus Paramyxine is diphyletic and Eptatretus paraphyletic. The complexity of Eptatretus and Paramyxine in the clade would not be solved until the farther departed P. cheni is included to form a new clade under the genus Eptatretus. The other clade of Myxininae contains but single genus Myxine.

The Ciona intestinalis genome: when the constraints are off

The recent genome sequencing of a non-vertebrate deuterostome, the ascidian tunicate Ciona intestinalis, makes a substantial contribution to the fields of evolutionary and developmental biology.1 Tunicates have some of the smallest bilaterian genomes, embryos with relatively few cells, fixed lineages and early determination of cell fates. Initial analyses of the C. intestinalis genome indicate that it has been evolving rapidly. Comparisons with other bilaterians show that C. intestinalis has lost a number of genes, and that many genes linked together in most other bilaterians have become uncoupled. In addition, a number of independent, lineage-specific gene duplications have been detected. These new results, although interesting in themselves, will take on a deeper significance once the genomes of additional invertebrate deuterostomes (e.g. echinoderms, hemichordates and amphioxus) have been sequenced. With such a broadened database, comparative genomics can begin to ask pointed questions about the relationship between the evolution of genomes and the evolution of body plans.

Complete mitochondrial DNA of the hagfish, Eptatretus burgeri: the comparative analysis of mitochondrial DNA sequences strongly supports the cyclostome monophyly

The phylogenetic position of cyclostomes, i.e., the relationships between hagfishes, lampreys, and jawed vertebrates is an unresolved problem. Anatomical data support the paraphyly of cyclostomes, whereas nuclear genes data support monophyly of cyclostomes. Previous results obtained using mitochondrial DNA are ambiguous, presumably due to a lack of informative sequences. By adding the complete mtDNA of a hagfish, Eptatretus burgeri, we have generated a novel data set for sequences of hagfishes and of lampreys. The addition of this mtDNA sequence to the 12 taxa we have already used becomes sufficient to obtain unambiguous results. This data set, which includes sequences of mtDNA of animals closely related to the lamprey/hagfish node, was used in a phylogenetic analysis with two independent statistical approaches and unequivocally supported the monophyly of cyclostomes. Thus molecular data, i.e., our results and those obtained using nuclear genes, conclude that hagfishes and lampreys form a clade.

5S ribosomal RNA

Ribosomes have been visualized in electron micrographs in 1943 but 5S rRNA was discovered 20 years later. The next four decades witnessed big advances in our understanding of the ribosome using biochemical, genetic and low resolution structural approaches. During those times many experimental data accumulates also on 5S rRNA, but its precise function remains unknown. To understand the role of this RNA in ribosome a high-resolution structure is urgently needed. Because the ribosome is a dynamic machine, details on the interaction of 5S rRNA with proteins within entire ribosome are required. Big progress in the structural analysis of ribosome will stimulate further understanding of 5S rRNA.

Gene duplication events producing muscle (M) and brain (B) isoforms of cytoplasmic creatine kinase: cDNA and deduced amino acid sequences from two lower chordates

Creatine kinase (CK) is coded for by at least four loci in higher vertebrates--two cytoplasmic isoforms, muscle (M) and brain (B), and two mitochondrial isoforms, sarcomeric and ubiquitous. M is expressed primarily in skeletal muscle, while B is expressed in a variety of cells, including cardiac and smooth muscle fibers, neurons, transport epithelia, and photoreceptors. M and B subunits form very stable homodimers (MM [M-CK], BB [B-CK]) and heterodimers (MB). M-CK is capable of binding to the M line of the myofibril, thereby creating an energy transfer microcompartment; BB and MB CKs are not. M- and B-like CKs are present in all vertebrates yet examined, including fish. Cytoplasmic, dimeric CKs are widely distributed in the invertebrates. The only available amino acid sequence for an invertebrate dimeric CK, that of the protostome polychaete Chaetopterus variopedatus, is just as similar to the vertebrate M isoform as to the B isoform. Echinoderms lack dimeric, cytoplasmic CKs, which appear to be replaced by a dimeric arginine kinase which evolved secondarily from CK. Thus, it is likely that the gene duplication event producing the M and B isoforms occurred after the divergence of the chordates from echinoderms. To narrow down the timing of this duplication event, we obtained the cDNA and deduced amino acid sequences of dimeric CKs from the tunicate Ciona intestinalis (subphylum Urochordata) and the lancelet Branchiostoma floridae (subphylum Cephalochordata). Our results show that these CKs are strikingly similar to both invertebrate and vertebrate CKs. However, phylogenetic analyses by neighbor-joining and parsimony show that these two enzymes appeared to have diverged before the point of divergence of the M and B isoforms. Thus, the gene duplication event for formation of the muscle and brain isoforms of CK most likely occurred during the radiation of the fish, a time noted for gene duplication events at a variety of other loci.

Euteleost fish genomes are characterized by expansion of gene families

The presence of additional hox clusters in the zebrafish has led to the hypothesis that there was a whole genome duplication at the origin of modern fish. To investigate the generality of this assumption, we analyzed all available actinopterygian fish gene families, and sequenced nuclear receptors from diverse teleost fish. The origin and timing of duplications was systematically determined by phylogenetic analysis. More genes are indeed found in zebrafish than in mouse. This abundance is shared by all major groups of euteleost fish, but not by eels. Phylogenetic analysis shows that it may result from frequent independent duplications, rather than from an ancestral genome duplication. We predict two zebrafish paralogs for most mouse or human genes, thus expressing a note of caution in functional comparison of fish and mammalian genomes. Redundancy appears to be the rule in fish developmental genetics. Finally, our results imply that the outcome of genome projects cannot be extrapolated easily between fish species.

A novel mitochondrial gene order in the crinoid echinoderm Florometra serratissima

The complete nucleotide sequence of the mitochondrial genome of the crinoid Florometra serratissima has been determined. It is a circular DNA molecule, 16,005 bp in length, containing the genes for 13 proteins, small and large ribosomal RNAs, and 22 transfer RNAs (tRNAs). Three regions of unassigned sequence (UAS) greater than 73 bp have been located. The largest, UAS I, is 432 bp long and exhibits sequence similarity to the putative mitochondrial control regions seen in other animals. UAS II (77 bp) and UAS III (73 bp) are located between the 5' ends of coding sequences and may play roles as bidirectional promoters. Analyses of nucleotide composition revealed that the major peptide-encoding strand is high in T and low in C. This bias is reflected in a specific pattern of codon usage. Molecular phylogenetic analyses based on cytochrome c oxidase (COI, COII, and COIII) amino acid and nucleotide sequences did not resolve all the relationships between echinoderm classes. The overall animal mitochondrial gene content has been maintained in the crinoid, but there is extensive rearrangement with respect to both the echinoid and the asteroid mtDNA gene maps. Florometra serratissima has a novel genome organization in a segment containing most of the tRNA genes, large and small rRNA genes, and the NADH dehydrogenase subunit 1 and 2 genes. Potential pathways and mechanisms for gene rearrangements between mitochondrial gene maps of echinoderm classes and vertebrates are discussed as indicators of early deuterostome phylogeny.

Amphioxus alcohol dehydrogenase is a class 3 form of single type and of structural conservation but with unique developmental expression

The coding region of amphioxus alcohol dehydrogenase class 3 (ADH3) has been characterized from two species, Branchiostoma lanceolatum and Branchiostoma floridae. The species variants have residue differences at positions that result in only marginal functional distinctions. Activity measurements show a class 3 glutathione-dependent formaldehyde dehydrogenase, with kcat/Km values about threefold those of the human class 3 ADH enzyme. Only a single ADH3 form is identified in each of the two amphioxus species, and no ethanol activity ascribed to other classes is detectable, supporting the conclusion that evolution of ethanol-active ADH classes by gene duplications occurred at early vertebrate radiation after the formation of the amphioxus lineage. Similarly, Southern blot analysis indicated that amphioxus ADH3 is encoded by a single gene present in the methylated fraction of the amphioxus genome and northern blots revealed a single 1.4-kb transcript. In situ experiments showed that amphioxus Adh3 expression is restricted to particular cell types in the embryos. Transcripts were first evident at the neurula stage and then located at the larval ventral region, in the intestinal epithelium. This tissue-specific pattern contrasts with the ubiquitous Adh3 expression in mammals.

Mitochondrial evidence on the phylogenetic position of caecilians (Amphibia: Gymnophiona)

The complete nucleotide sequence (17,005 bp) of the mitochondrial genome of the caecilian Typhlonectes natans (Gymnophiona, Amphibia) was determined. This molecule is characterized by two distinctive genomic features: there are seven large 109-bp tandem repeats in the control region, and the sequence for the putative origin of replication of the L strand can potentially fold into two alternative secondary structures (one including part of the tRNA(Cys)). The new sequence data were used to assess the phylogenetic position of caecilians and to gain insights into the origin of living amphibians (frogs, salamanders, and caecilians). Phylogenetic analyses of two data sets-one combining protein-coding genes and the other combining tRNA genes-strongly supported a caecilian + frog clade and, hence, monophyly of modern amphibians. These two data sets could not further resolve relationships among the coelacanth, lungfishes, and tetrapods, but strongly supported diapsid affinities of turtles. Phylogenetic relationships among a larger set of species of frogs, salamanders, and caecilians were estimated with a mitochondrial rRNA data set. Maximum parsimony analysis of this latter data set also recovered monophyly of living amphibians and favored a frog + salamander (Batrachia) relationship. However, bootstrap support was only moderate at these nodes. This is likely due to an extensive among-site rate heterogeneity in the rRNA data set and the narrow window of time in which the three main groups of living amphibians were originated.

The complete nucleotide sequence of the mitochondrial DNA of the agnathan Lampetra fluviatilis: bearings on the phylogeny of cyclostomes

There are two competing theories about the interrelationships of craniates: the cyclostome theory assumes that lampreys and hagfishes are a clade, the cyclostomes, whose sister group is the jawed vertebrates (gnathostomes); the vertebrate theory assumes that lampreys and gnathostomes are a clade, the vertebrates, whose sister group is hagfishes. The vertebrate theory is best supported by a number of unique anatomical and physiological characters. Molecular sequence data from 18S and 28S rRNA genes rather support the cyclostome theory, but mtDNA sequence of Myxine glutinosa rather supports the vertebrate theory. Additional molecular data are thus needed to elucidate this three-taxon problem. We determined the complete nucleotide sequence of the mtDNA of the lamprey Lampetra fluviatilis. The mtDNA of L. fluviatilis possesses the same genomic organization as Petromyzon marinus, which validates this gene order as a synapomorphy of lampreys. The mtDNA sequence of L. fluviatilis was used in combination with relevant mtDNA sequences for an approach to the hagfish/lamprey relationships using the maximum-parsimony, neighbor-joining, and maximum-likelihood methods. Although trees compatible with our present knowledge of the phylogeny of craniates can be reconstructed by using the three methods, the data collected do not support the vertebrate or the cyclostome hypothesis. The present data set does not allow the resolution of this three-taxon problem, and new kinds of data, such as nuclear DNA sequences, need to be collected.

Mitochondrial genomes of Galathealinum, Helobdella, and Platynereis: sequence and gene arrangement comparisons indicate that Pogonophora is not a phylum and Annelida and Arthropoda are not sister taxa

We report a contiguous region of more than half (> 7,500 nt) of the mitochondrial genomes for Platynereis dumerii (Annelida: Polychaeta), Helobdella robusta (Annelida: Hirudinida), and Galathealinum brachiosum (Pogonophora: Perviata). The relative arrangements of all 22 genes identified for Helobdella and Galathealinum are identical to one another and to their arrangements in the mtDNA of the previously studied oligochaete annelid Lumbricus. In contrast, Platynereis differs from these taxa in the positions of several tRNA genes and in having two additional tRNA genes (trnC and trnM) and a large noncoding sequence in this region. Comparisons of relative gene arrangements and of the nucleotide and inferred amino acid sequences among these and other published taxa provide strong support for an annelid-mollusk clade that excludes arthropods, and for the inclusion of pogonophorans within Annelida, rather than giving them separate phylum status. Gene arrangement comparisons include the first use of a recently described method on previously unpublished data. Although a variety of alternative initiation codons are typically used by mitochondrial protein-encoding genes, ATG appears to be the initiator for all but one reported here. The large noncoding region (1,091 nt) identified in Platynereis has no significant sequence similarity to the noncoding region of Lumbricus, although each contains runs of TA dinucleotides and of homopolymers, which could potentially serve as signaling elements. There is strong bias for synonymous codon usage in Helobdella and especially in Galathealinum. In this latter taxon, 5 codons are completely unused, 13 are used three or fewer times, and G appears at third codon positions in only 26 of the 2,236 codons. Nucleotide composition bias appears to influence amino acid composition of the proteins.

5S ribosomal RNA database Y2K

This paper presents the updated version (Y2K) of the database of ribosomal 5S ribonucleic acids (5S rRNA) and their genes (5S rDNA), http://rose.man/poznan.pl/5SData/index.html. This edition of the database contains 1985primary structures of 5S rRNA and 5S rDNA. They include 60 archaebacterial, 470 eubacterial, 63 plastid, nine mitochondrial and 1383 eukaryotic sequences. The nucleotide sequences of the 5S rRNAs or 5S rDNAs are divided according to the taxonomic position of the source organisms.

Analysis of the complete mitochondrial DNA sequence of the brachiopod terebratulina retusa places Brachiopoda within the protostomes

Brachiopod phylogeny is still a controversial subject. Analyses using nuclear 18SrRNA and mitochondrial 12SrDNA sequences place them within the protostomes but some recent interpretations of morphological data support a relationship with deuterostomes. In order to investigate brachiopod affinities within the metazoa further, we compared the gene arrangement on the brachiopod mitochondrial genome with several metazoan taxa. The complete (15 451 bp) mitochondrial DNA (mtDNA) sequence of the articulate brachiopod Terebratulina retusa was determined from two overlapping long polymerase chain reaction products. All the genes are encoded on the same strand and gene order comparisons showed that.only one major rearrangement is required to interconvert the T. retusa and Katharina tunicata (Mollusca: Polvplacophora) mitochondrial genomes. The partial mtDNA sequence of the prosobranch mollusc Littorina saxatilis shows complete congruence with the T. rehtusa gene arrangement with regard to the ribosomal and protein coding genes. This high similarity in gene arrangement is the first to be reported within the protostomes. Sequence analyses of mitochondrial protein coding genes also support a close relationship of the brachiopod with molluscs and annelids, thus supporting the clade Lophotrochozoa. Though being highly informative, sequence analyses of the mitochondrial protein coding genes failed to resolve the branching order within the lophotrochozoa.

Evolutionary genomics in Metazoa: the mitochondrial DNA as a model system

One of the most important aspects of mitochondrial (mt) genome evolution in Metazoa is constancy of size and gene content of mtDNA, whose plasticity is maintained through a great variety of gene rearrangements probably mediated by tRNA genes. The trend of mtDNA to maintain the same genetic structure within a phylum (e.g., Chordata) is generally accepted, although more recent reports show that a considerable number of transpositions are observed also between closely related organisms. Base composition of mtDNA is extremely variable. Genome GC content is often low and, when it increases, the two complementary bases distribute asymmetrically, creating, particularly in vertebrates, a negative GC-skew. In mammals, we have found coding strand base composition and average degree of gene conservation to be related to the asymmetric replication mechanism of mtDNA. A quantitative measurement of mtDNA evolutionary rate has revealed that each of the various components has a different evolutionary rate. Non-synonymous rates are gene specific and fall in a range comparable to that of nuclear genes, whereas synonymous rates are about 22-fold higher in mt than in nuclear genes. tRNA genes are among the most conserved but, when compared to their nuclear counterparts, they evolve 100 times faster. Finally, we describe some molecular phylogenetic reconstructions which have produced unexpected outcomes, and might change our vision of the classification of living organisms.

Animal mitochondrial genomes

Animal mitochondrial DNA is a small, extrachromosomal genome, typically approximately 16 kb in size. With few exceptions, all animal mitochondrial genomes contain the same 37 genes: two for rRNAs, 13 for proteins and 22 for tRNAs. The products of these genes, along with RNAs and proteins imported from the cytoplasm, endow mitochondria with their own systems for DNA replication, transcription, mRNA processing and translation of proteins. The study of these genomes as they function in mitochondrial systems-'mitochondrial genomics'-serves as a model for genome evolution. Furthermore, the comparison of animal mitochondrial gene arrangements has become a very powerful means for inferring ancient evolutionary relationships, since rearrangements appear to be unique, generally rare events that are unlikely to arise independently in separate evolutionary lineages. Complete mitochondrial gene arrangements have been published for 58 chordate species and 29 non-chordate species, and partial arrangements for hundreds of other taxa. This review compares and summarizes these gene arrangements and points out some of the questions that may be addressed by comparing mitochondrial systems.

5S Ribosomal RNA Data Bank

This paper presents the updated version of the data base of ribosomal 5S ribonucleic acids (5S rRNA) and their genes (5S rDNA). This edition of the data bank contains 1889 primary structures of 5S rRNA and 5S rDNA. These include 60 archaebacterial, 439 eubacterial, 63 plastid, 9 mitochondrial and 1318 eukaryotic sequences. The nucleotide sequences of 5S rRNAs or 5S rDNAs are divided according to the taxonomic position of organisms. The sequences stored in the database can be viewed and retrieved using the taxonomic browser at the URL: http://rose.man.poznan.pl/5SData/5SRNA.html++ +