Complete Genome Analysis of Three Live Attenuated Rinderpest Virus Vaccine Strains Derived through Serial Passages in Different Culture Systems

The genomes of three South Korean Rinderpest virus vaccine strains (L72, LA77, and LA96) were analyzed in order to investigate their genetic variability. These three vaccine strains were all derived from the same virus strain origin (Fusan) through repeated passages in different culture systems. The full genome length of the three strains was 15,882 nucleotides, and the sequence similarity between the three South Korean RPV strains at the nucleotide level was 98.1 to 98.9%. The genetic distance between Nakamura III, L72, LA77, LA96, and LATC06 and the Kabete strain was greater than that between the Fusan and Kabete strains for the P, V, and C genes. The difference in pathogenicity among these strains might be due to the V gene, which has a positive (>1) selection ratio based on the analysis of synonymous (dS) and nonsynonymous (dN) substitution rates (dN/dS ratio [ω]).     

Loss of the IR region in conifer plastomes: Changes in the selection pressure and substitution rate of protein‐coding genes

Plastid genomes (plastomes) have a quadripartite structure, but some species have drastically reduced or lost inverted repeat (IR) regions. IR regions are important for genome stability and the evolution rate. In the evolutionary process of gymnosperms, the typical IRs of conifers were lost, possibly affecting the evolutionary rate and selection pressure of genomic protein‐coding genes. In this study, we selected 78 gymnosperm species (51 genera, 13 families) for evolutionary analysis. The selection pressure analysis results showed that negative selection effects were detected in all 50 common genes. Among them, six genes in conifers had higher ω values than non‐conifers, and 12 genes had lower ω values. The evolutionary rate analysis results showed that 9 of 50 common genes differed between conifers and non‐conifers. It is more obvious that in non‐conifers, the rates of psbA (trst, trsv, ratio, dN, dS, and ω) were 2.6‐ to 3.1‐fold of conifers. In conifers, trsv, ratio, dN, dS, and ω of ycf2 were 1.2‐ to 3.6‐fold of non‐conifers. In addition, the evolution rate of ycf2 in the IR was significantly reduced. psbA is undergoing dynamic change, with an abnormally high evolution rate as a small portion of it enters the IR region. Although conifers have lost the typical IR regions, we detected no change in the substitution rate or selection pressure of most protein‐coding genes due to gene function, plant habitat, or newly acquired IRs.     

Phylogenomics and Coalescent Analyses Resolve Extant Seed Plant Relationships

The extant seed plants include more than 260,000 species that belong to five main lineages: angiosperms, conifers, cycads, Ginkgo, and gnetophytes. Despite tremendous effort using molecular data, phylogenetic relationships among these five lineages remain uncertain. Here, we provide the first broad coalescent-based species tree estimation of seed plants using genome-scale nuclear and plastid data By incorporating 305 nuclear genes and 47 plastid genes from 14 species, we identify that i) extant gymnosperms (i.e., conifers, cycads, Ginkgo, and gnetophytes) are monophyletic, ii) gnetophytes exhibit discordant placements within conifers between their nuclear and plastid genomes, and iii) cycads plus Ginkgo form a clade that is sister to all remaining extant gymnosperms. We additionally observe that the placement of Ginkgo inferred from coalescent analyses is congruent across different nucleotide rate partitions. In contrast, the standard concatenation method produces strongly supported, but incongruent placements of Ginkgo between slow- and fast-evolving sites. Specifically, fast-evolving sites yield relationships in conflict with coalescent analyses. We hypothesize that this incongruence may be related to the way in which concatenation methods treat sites with elevated nucleotide substitution rates. More empirical and simulation investigations are needed to understand this potential weakness of concatenation methods.     

Potent antibody-mediated neutralization and evolution of antigenic escape variants of simian immunodeficiency virus strain SIVmac239 in vivo

Here, we describe the evolution of antigenic escape variants in a rhesus macaque that developed unusually high neutralizing antibody titers to SIVmac239. By 42 weeks postinfection, 50% neutralization of SIVmac239 was achieved with plasma dilutions of 1:1,000. Testing of purified immunoglobulin confirmed that the neutralizing activity was antibody mediated. Despite the potency of the neutralizing antibody response, the animal displayed a typical viral load profile and progressed to terminal AIDS with a normal time course. Viral envelope sequences from week 16 and week 42 plasma contained an excess of nonsynonymous substitutions, predominantly in V1 and V4, including individual sites with ratios of nonsynonymous to synonymous substitution rates (dN/dS) highly suggestive of strong positive selection. Recombinant viruses encoding envelope sequences isolated from these time points remained resistant to neutralization by all longitudinal plasma samples, revealing the failure of the animal to mount secondary responses to the escaped variants. Substitutions at two sites with significant dN/dS values, one in V1 and one in V4, were independently sufficient to confer nearly complete resistance to neutralization. Substitutions at three additional sites, one in V4 and two in gp41, conferred moderate to high levels of resistance when tested individually. All the amino acid changes leading to escape resulted from single nucleotide substitutions. The observation that antigenic escape resulted from individual, single amino acid replacements at sites well separated in current structural models of Env indicates that the virus can utilize multiple independent pathways to rapidly achieve similar levels of resistance.     

An evolutionary analysis of genome expansion and pathogenicity in Escherichia coli

Background

There are several studies describing loss of genes through reductive evolution in microbes, but how selective forces are associated with genome expansion due to horizontal gene transfer (HGT) has not received similar attention. The aim of this study was therefore to examine how selective pressures influence genome expansion in 53 fully sequenced and assembled Escherichia coli strains. We also explored potential connections between genome expansion and the attainment of virulence factors. This was performed using estimations of several genomic parameters such as AT content, genomic drift (measured using relative entropy), genome size and estimated HGT size, which were subsequently compared to analogous parameters computed from the core genome consisting of 1729 genes common to the 53 E. coli strains. Moreover, we analyzed how selective pressures (quantified using relative entropy and dN/dS), acting on the E. coli core genome, influenced lineage and phylogroup formation.

Results

Hierarchical clustering of dS and dN estimations from the E. coli core genome resulted in phylogenetic trees with topologies in agreement with known E. coli taxonomy and phylogroups. High values of dS, compared to dN, indicate that the E. coli core genome has been subjected to substantial purifying selection over time; significantly more than the non-core part of the genome (p<0.001). This is further supported by a linear association between strain-wise dS and dN values (β = 26.94 ± 0.44, R²~0.98, p<0.001). The non-core part of the genome was also significantly more AT-rich (p<0.001) than the core genome and E. coli genome size correlated with estimated HGT size (p<0.001). In addition, genome size (p<0.001), AT content (p<0.001) as well as estimated HGT size (p<0.005) were all associated with the presence of virulence factors, suggesting that pathogenicity traits in E. coli are largely attained through HGT. No associations were found between selective pressures operating on the E. coli core genome, as estimated using relative entropy, and genome size (p~0.98).

Conclusions

On a larger time frame, genome expansion in E. coli, which is significantly associated with the acquisition of virulence factors, appears to be independent of selective forces operating on the core genome.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-882) contains supplementary material, which is available to authorized users.     

Why Genes Evolve Faster on Secondary Chromosomes in Bacteria

In bacterial genomes composed of more than one chromosome, one replicon is typically larger, harbors more essential genes than the others, and is considered primary. The greater variability of secondary chromosomes among related taxa has led to the theory that they serve as an accessory genome for specific niches or conditions. By this rationale, purifying selection should be weaker on genes on secondary chromosomes because of their reduced necessity or usage. To test this hypothesis we selected bacterial genomes composed of multiple chromosomes from two genera, Burkholderia and Vibrio, and quantified the evolutionary rates (dN and dS) of all orthologs within each genus. Both evolutionary rate parameters were faster among orthologs found on secondary chromosomes than those on the primary chromosome. Further, in every bacterial genome with multiple chromosomes that we studied, genes on secondary chromosomes exhibited significantly weaker codon usage bias than those on primary chromosomes. Faster evolution and reduced codon bias could in turn result from global effects of chromosome position, as genes on secondary chromosomes experience reduced dosage and expression due to their delayed replication, or selection on specific gene attributes. These alternatives were evaluated using orthologs common to genomes with multiple chromosomes and genomes with single chromosomes. Analysis of these ortholog sets suggested that inherently fast-evolving genes tend to be sorted to secondary chromosomes when they arise; however, prolonged evolution on a secondary chromosome further accelerated substitution rates. In summary, secondary chromosomes in bacteria are evolutionary test beds where genes are weakly preserved and evolve more rapidly, likely because they are used less frequently.     

Evidence of the accumulation of allele-specific non-synonymous substitutions in the young region of recombination suppression within the mating-type chromosomes of Neurospora tetrasperma

Currently, little is known about the origin and early evolution of sex chromosomes. This is largely due to the fact that ancient non-recombining sex chromosomes are highly degenerated, and thus provide little information about the early genomic events in their evolution. The Neurospora tetrasperma mating-type (mat) chromosomes contain a young (<6 Mya) and large region (>6.6 Mb) of suppressed recombination, thereby providing a model system to study early stages of sex chromosome evolution. Here, we examined alleles of 207 genes located on the N. tetrasperma mat a and mat A chromosomes to test for signs of genomic alterations at the protein level in the young region of recombination suppression. We report that the N. tetrasperma mat a and mat A chromosomes have each independently accumulated allele-specific non-synonymous codon substitutions in a time-dependent, and gene-specific manner in the recombinationally suppressed region. In addition, examination of the ratio (ω) of non-synonymous substitutions (dN) to synonymous substitutions (dS) using maximum likelihood analyses, indicates that such changes are associated with relaxed purifying selection, a finding consistent with genomic degeneration. We also reveal that sex specific biases in mutation rates or selection pressures are not necessary for genomic alterations in sex chromosomes, and that recombination suppression in itself is sufficient to explain these results. The present findings extend our current understanding of genomic events associated within the young region of recombination suppression in these fungal sex-regulating chromosomes.     

Mitonuclear Coevolution,but not Nuclear Compensation,Drives Evolution of OXPHOS Complexes in Bivalves

In Metazoa, four out of five complexes involved in oxidative phosphorylation (OXPHOS) are formed by subunits encoded by both the mitochondrial (mtDNA) and nuclear (nuDNA) genomes, leading to the expectation of mitonuclear coevolution. Previous studies have supported coadaptation of mitochondria-encoded (mtOXPHOS) and nuclear-encoded OXPHOS (nuOXPHOS) subunits, often specifically interpreted with regard to the “nuclear compensation hypothesis,” a specific form of mitonuclear coevolution where nuclear genes compensate for deleterious mitochondrial mutations due to less efficient mitochondrial selection. In this study, we analyzed patterns of sequence evolution of 79 OXPHOS subunits in 31 bivalve species, a taxon showing extraordinary mtDNA variability and including species with “doubly uniparental” mtDNA inheritance. Our data showed strong and clear signals of mitonuclear coevolution. NuOXPHOS subunits had concordant topologies with mtOXPHOS subunits, contrary to previous phylogenies based on nuclear genes lacking mt interactions. Evolutionary rates between mt and nuOXPHOS subunits were also highly correlated compared with non-OXPHO-interacting nuclear genes. Nuclear subunits of chimeric OXPHOS complexes (I, III, IV, and V) also had higher dN/dS ratios than Complex II, which is formed exclusively by nuDNA-encoded subunits. However, we did not find evidence of nuclear compensation: mitochondria-encoded subunits showed similar dN/dS ratios compared with nuclear-encoded subunits, contrary to most previously studied bilaterian animals. Moreover, no site-specific signals of compensatory positive selection were detected in nuOXPHOS genes. Our analyses extend the evidence for mitonuclear coevolution to a new taxonomic group, but we propose a reconsideration of the nuclear compensation hypothesis.     

Fast evolution of the retroprocessed mitochondrial rps3 gene in Conifer II and further evidence for the phylogeny of gymnosperms

The popular view that plant mitochondrial genome evolves slowly in sequence has been recently challenged by the extraordinarily high substitution rates of mtDNA documented mainly from several angiosperm genera, but high substitution rate acceleration accompanied with great length variation has been very rarely reported in plant mitochondrial genes. Here, we studied evolution of the mitochondrial rps3 gene that encodes the ribosomal small subunit protein 3 and found a dramatically high variation in both length and sequence of an exon region of it in Conifer II. A sequence comparison between cDNA and genomic DNA showed that there are no RNA editing sites in the Conifer II rps3 gene. Southern blotting analyses of the total DNA and mtDNA, together with the real-time PCR analysis, showed that rps3 exists as a single mitochondrial locus in gymnosperms. It is very likely that the Conifer II rps3 gene has experienced retroprocessing, i.e., the re-integration of its cDNA into the mitochondrial genome, followed by an evolutionary acceleration due to the intron loss. In addition, the phylogenetic analysis of rps3 supports the sister relationship between conifers and Gnetales. In particular, the monophyly of conifer II is strongly supported by the shared loss of two rps3 introns. Our results also indicate that the mitochondrial gene tree would be affected in topology when the “edited” paralogs are analyzed together with their genomic sequences.     

Evolution of Reproductive Organs in Land Plants

LEAFY gene is the positive regulator of the MADS-box genes in flower primordia. The number of MADS-box genes presumably increased by gene duplications before the divergence of ferns and seed plants. Most MADS-box genes in ferns are expressed similarly in both vegetative and reproductive organs, while in gymnosperms, some MADS-box genes are specifically expressed in reproductive organs. This suggests that (1) the increase in the number of MADS-box genes and (2) the subsequent recruitment of some MADS-box genes as homeotic selector genes were important for the evolution of complex reproductive organs. The phylogenetic tree including both angiosperm and gymnosperm MADS-box genes indicates the loss of the A-function genes in the gymnosperm lineage, which is presumably related to the absence of perianths in extant gymnosperms. Comparison of expression patterns of orthologous MADS-box genes in angiosperms, Gnetales, and conifers supports the sister relationship of Gnetales and conifers over that of Gnetales and angiosperms predicted by phylogenetic trees based on amino acid and nucleotide sequences. Received 30 July 1999/ Accepted in revised form 9 September 1999     

Testing the neutral theory of molecular evolution using genomic data: a comparison of the human and bovine transcriptome

Despite growing evidence of rapid evolution in protein coding genes, the contribution of positive selection to intra- and interspecific differences in protein coding regions of the genome is unclear. We attempted to see if genes coding for secreted proteins and genes with narrow expression, specifically those preferentially expressed in the mammary gland, have diverged at a faster rate between domestic cattle (Bos taurus) and humans (Homo sapiens) than other genes and whether positive selection is responsible. Using a large data set, we identified groups of genes based on secretion and expression patterns and compared them for the rate of nonsynonymous (dN) and synonymous (dS) substitutions per site and the number of radical (Dr) and conservative (Dc) amino acid substitutions. We found evidence of rapid evolution in genes with narrow expression, especially for those expressed in the liver and mammary gland and for genes coding for secreted proteins. We compared common human polymorphism data with human-cattle divergence and found that genes with high evolutionary rates in human-cattle divergence also had a large number of common human polymorphisms. This argues against positive selection causing rapid divergence in these groups of genes. In most cases dN/dS ratios were lower in human-cattle divergence than in common human polymorphism presumably due to differences in the effectiveness of purifying selection between long-term divergence and short-term polymorphism.     

A new Permian gnetalean cone as fossil evidence for supporting current molecular phylogeny

• Background and aims The order Gnetales has been the central focus of controversy in seed plant phylogeny. Traditional treatment of morphology supports the anthophyte hypothesis with Gnetales sister to angiosperms but current molecular data reject this hypothesis. A new fossil gnetalean cone, Palaeognetaleana auspicia gen. et sp. nov., is reported from the Upper Permian in North China, and its phylogenic implications are considered.• Methods Samples of cones from the upper part of the Upper Permian redbeds of Baode section, northwestern Shanxi Province, China, were examined.• Key results The cone is characterized by its unusual nature of reproduction that combines features of post-Triassic gnetaleans and some of the Palaeozoic conifers. It is made up of a number of imbricate axillary units, each simply formed by an ovule and a subtending bract, which may be comparable with the axillary seed-scale complex of some of the Palaeozoic conifer cones. The cone exhibits at least a partially bisexual character that appears to have pollen sacs with monosulcate ribbed pollen grains and sessile, asymmetric, and radiospermic ovules. The ovule has an integument of three envelopes: an outer one of pointed scales; a middle sclerified one; and an inner cuticle that extends upward into a micropyle with an oblique tip.• Conclusions The new Permian cone has unequivocal affinity with the Gnetales. The fossil has considerably extended the divergence time of the Gnetales from 140 (210?) back to 270 myr ago and, therefore, provides the first significant fossil evidence to support the current conclusion based on molecular data of seed plants, i.e. monophyletic gymnosperms, comprising the Gnetales are closely related to conifers.Key words: Palaeognetaleana auspicia (gen. et sp. nov.), Gnetales, bisexual fossil cone, phylogeny of seed plants, radiation of gymnosperms, molecular data, Late Permian, North China     

Phylogenetic analysis of the acetyl-CoA carboxylase and 3-phosphoglycerate kinase loci in wheat and other grasses

We have applied a two-gene system based on the sequences of nuclear genes encoding multi-domain plastid acetyl-CoA carboxylase (ACCase) and plastid 3-phosphoglycerate kinase (PGK) to study grass evolution. Our analysis revealed that these genes are single-copy in most of the grass species studied, allowing the establishment of orthologous relationships between them. These relationships are consistent with the known facts of their evolution: the eukaryotic origin of the plastid ACCase, created by duplication of a gene encoding the cytosolic multi-domain ACCase gene early in grass evolution, and the prokaryotic (endosymbiont) origin of the plastid PGK. The major phylogenetic relationships among grasses deduced from the nucleotide sequence comparisons of ACCase and PGK genes are consistent with each other and with the milestones of grass evolution revealed by other methods. Nucleotide substitution rates were calculated based on multiple pairwise sequence comparisons. On a relative basis, with the divergence of the Pooideae and Panicoideae subfamilies set at 60 million years ago (MYA), events leading to the Triticum/Aegilops complex occurred at the following intervals: divergence of Lolium (Lolium rigidum) at 35 MYA, divergence of Hordeum (Hordeum vulgare) at 11 MYA and divergence of Secale (Secale cereale) at 7 MYA. On the same scale, gene duplication leading to the multi-domain plastid ACCase in grasses occurred at 129 MYA, divergence of grass and dicot plastid PGK genes at 137 MYA, and divergence of grass and dicot cytosolic PGK genes at 155 MYA. The ACCase and PGK genes provide a well-understood two-locus system to study grass phylogeny, evolution and systematics.     

C-Glycosylflavones in Gnetum gnemon

The relationship of the Gnetales to other gymnosperms is presently ill understood. The occurrence of C-glycosylflavones as the only detactable flav     

Phylogeny of seed plants based on evidence from eight genes

Relationships among the five groups of extant seed plants (cycads, Ginkgo, conifers, Gnetales, and angiosperms) remain uncertain. To explore relationships among groups of extant seed plants further and to attempt to explain the conflict among molecular data sets, we assembled a data set of four plastid (cpDNA) genes (rbcL, atpB, psaA, and psbB), three mitochondrial (mtDNA) genes (mtSSU, coxI, and atpA), and one nuclear gene (18S rDNA) for 19 exemplars representing the five groups of living seed plants. Analyses of the combined eight-gene data set (15?772 base pairs/taxon) with maximum parsimony (MP), maximum likelihood (ML), and Bayesian approaches reveal a gymnosperm clade that is sister to angiosperms. Within the gymnosperms, a conifer clade includes Gnetales as sister to Pinaceae. Cycads and Ginkgo are either successive sisters to this conifer clade (including Gnetales) or a clade that is sister to conifers and Gnetales. All analyses of the mtDNA partition and ML analyses of the nuclear partition yield very similar topologies. However, MP analyses of the combined cpDNA genes place Gnetales as sister to all other seed plants with strong bootstrap support, whereas ML and Bayesian analyses of the cpDNA data set place Gnetales as sister to Pinaceae. Maximum parsimony and ML analyses of first and second codon positions of the cpDNA partiation also place Gnetales as sister to Pinaceae. In contrast, MP analyses of third codon positions place Gnetales as sister to other seed plants, although ML analyses of third codon positions place Gnetales with Pinaceae. Thus, most of the discrepancies in seed plant topologies involve third codon positions of cpDNA genes. The likelihood ratio (LR) and Shimodaira-Hasegasa (SH) tests were applied to the cpDNA data. The preferred topology based on the LR test is that Gnetales are sister to Pseudotsuga. The SH test based on first and second codon and all three codon positions indicated that there is no significant difference between the best topology (Gnetales sister to Pseudotsuga) and Gnetales sister to a conifer clade. However, there is a significant difference between the best topology and topologies in which Gnetales are sister to the rest of the seed plants or Gnetales sister to angiosperms.     

Whole-gene positive selection, elevated synonymous substitution rates, duplication, and indel evolution of the chloroplast clpP1 gene

Background

Synonymous DNA substitution rates in the plant chloroplast genome are generally relatively slow and lineage dependent. Non-synonymous rates are usually even slower due to purifying selection acting on the genes. Positive selection is expected to speed up non-synonymous substitution rates, whereas synonymous rates are expected to be unaffected. Until recently, positive selection has seldom been observed in chloroplast genes, and large-scale structural rearrangements leading to gene duplications are hitherto supposed to be rare.

Methodology/Principle Findings

We found high substitution rates in the exons of the plastid clpP1 gene in Oenothera (the Evening Primrose family) and three separate lineages in the tribe Sileneae (Caryophyllaceae, the Carnation family). Introns have been lost in some of the lineages, but where present, the intron sequences have substitution rates similar to those found in other introns of their genomes. The elevated substitution rates of clpP1 are associated with statistically significant whole-gene positive selection in three branches of the phylogeny. In two of the lineages we found multiple copies of the gene. Neighboring genes present in the duplicated fragments do not show signs of elevated substitution rates or positive selection. Although non-synonymous substitutions account for most of the increase in substitution rates, synonymous rates are also markedly elevated in some lineages. Whereas plant clpP1 genes experiencing negative (purifying) selection are characterized by having very conserved lengths, genes under positive selection often have large insertions of more or less repetitive amino acid sequence motifs.

Conclusions/Significance

We found positive selection of the clpP1 gene in various plant lineages to correlated with repeated duplication of the clpP1 gene and surrounding regions, repetitive amino acid sequences, and increase in synonymous substitution rates. The present study sheds light on the controversial issue of whether negative or positive selection is to be expected after gene duplications by providing evidence for the latter alternative. The observed increase in synonymous substitution rates in some of the lineages indicates that the detection of positive selection may be obscured under such circumstances. Future studies are required to explore the functional significance of the large inserted repeated amino acid motifs, as well as the possibility that synonymous substitution rates may be affected by positive selection.     

Phylogeny and Divergence Times of Gymnosperms Inferred from Single-Copy Nuclear Genes

Phylogenetic reconstruction is fundamental to study evolutionary biology and historical biogeography. However, there was not a molecular phylogeny of gymnosperms represented by extensive sampling at the genus level, and most published phylogenies of this group were constructed based on cytoplasmic DNA markers and/or the multi-copy nuclear ribosomal DNA. In this study, we use LFY and NLY, two single-copy nuclear genes that originated from an ancient gene duplication in the ancestor of seed plants, to reconstruct the phylogeny and estimate divergence times of gymnosperms based on a complete sampling of extant genera. The results indicate that the combined LFY and NLY coding sequences can resolve interfamilial relationships of gymnosperms and intergeneric relationships of most families. Moreover, the addition of intron sequences can improve the resolution in Podocarpaceae but not in cycads, although divergence times of the cycad genera are similar to or longer than those of the Podocarpaceae genera. Our study strongly supports cycads as the basal-most lineage of gymnosperms rather than sister to Ginkgoaceae, and a sister relationship between Podocarpaceae and Araucariaceae and between Cephalotaxaceae-Taxaceae and Cupressaceae. In addition, intergeneric relationships of some families that were controversial, and the relationships between Taxaceae and Cephalotaxaceae and between conifers and Gnetales are discussed based on the nuclear gene evidence. The molecular dating analysis suggests that drastic extinctions occurred in the early evolution of gymnosperms, and extant coniferous genera in the Northern Hemisphere are older than those in the Southern Hemisphere on average. This study provides an evolutionary framework for future studies on gymnosperms.     

Complete Plastid Genome of the Brown Alga Costaria costata (Laminariales,Phaeophyceae)

Costaria costata is a commercially and industrially important brown alga. In this study, we used next-generation sequencing to determine the complete plastid genome of C. costata. The genome consists of a 129,947 bp circular DNA molecule with an A+T content of 69.13% encoding a standard set of six ribosomal RNA genes, 27 transfer RNA genes, and 137 protein-coding genes with two conserved open reading frames (ORFs). The overall genome structure of C. costata is nearly the same as those of Saccharina japonica and Undaria pinnatifida. The plastid genomes of these three algal species retain a strong conservation of the GTG start codon while infrequently using TGA as a stop codon. In this regard, they differ substantially from the plastid genomes of Ectocarpus siliculosus and Fucus vesiculosus. Analysis of the nucleic acid substitution rates of the Laminariales plastid genes revealed that the petF gene has the highest substitution rate and the petN gene contains no substitution over its complete length. The variation in plastid genes between C. costata and S. japonica is lower than that between C. costata and U. pinnatifida as well as that between U. pinnatifida and S. japonica. Phylogenetic analyses demonstrated that C. costata and U. pinnatifida have a closer genetic relationship. We also identified two gene length mutations caused by the insertion or deletion of repeated sequences, which suggest a mechanism of gene length mutation that may be one of the key explanations for the genetic variation in plastid genomes.