The complete mitochondrial genome of the grand jackknife clam, Solen grandis (Bivalvia: Solenidae): a novel gene order and unusual non-coding region

Molluscs in general, and bivalves in particular, exhibit an extraordinary degree of mitochondrial gene order variation when compared with other metazoans. The complete mitochondrial genome of Solen grandis (Bivalvia: Solenidae) was determined using long-PCR and genome walking techniques. The entire mitochondrial genome sequence of S. grandis is 16,784 bp in length, and contains 36 genes including 12 protein-coding genes (atp8 is absent), 2 ribosomal RNAs, and 22 tRNAs. All genes are encoded on the same strand. Compared with other species, it bears a novel gene order. Besides these, we find a peculiar non-coding region of 435 bp with a microsatellite-like (TA)₁₂ element, poly-structures and many hairpin structures. In contrast to the available heterodont mitochondrial genomes from GenBank, the complete mtDNA of S. grandis has the shortest cox3 gene, and the longest atp6, nad4, nad5 genes.     

Comparative mitochondrial genome analysis reveals the evolutionary rearrangement mechanism in Brassica

The genus Brassica has many species that are important for oil, vegetable and other food products. Three mitochondrial genome types (mitotype) originated from its common ancestor. In this paper, a B. nigra mitochondrial main circle genome with 232,407 bp was generated through de novo assembly. Synteny analysis showed that the mitochondrial genomes of B. rapa and B. oleracea had a better syntenic relationship than B. nigra. Principal components analysis and development of a phylogenetic tree indicated maternal ancestors of three allotetraploid species in Us triangle of Brassica. Diversified mitotypes were found in allotetraploid B. napus, in which napus‐type B. napus was derived from B. oleracea, while polima‐type B. napus was inherited from B. rapa. In addition, the mitochondrial genome of napus‐type B. napus was closer to botrytis‐type than capitata‐type B. oleracea. The sub‐stoichiometric shifting of several mitochondrial genes suggested that mitochondrial genome rearrangement underwent evolutionary selection during domestication and/or plant breeding. Our findings clarify the role of diploid species in the maternal origin of allotetraploid species in Brassica and suggest the possibility of breeding selection of the mitochondrial genome.     

Interspecies introgressive hybridization in spiny frogs Quasipaa (Family Dicroglossidae) revealed by analyses on multiple mitochondrial and nuclear genes

Introgression may lead to discordant patterns of variation among loci and traits. For example, previous phylogeographic studies on the genus Quasipaa detected signs of genetic introgression from genetically and morphologically divergent Quasipaa shini or Quasipaa spinosa. In this study, we used mitochondrial and nuclear DNA sequence data to verify the widespread introgressive hybridization in the closely related species of the genus Quasipaa, evaluate the level of genetic diversity, and reveal the formation mechanism of introgressive hybridization. In Longsheng, Guangxi Province, signs of asymmetrical nuclear introgression were detected between Quasipaa boulengeri and Q. shini. Unidirectional mitochondrial introgression was revealed from Q. spinosa to Q. shini. By contrast, bidirectional mitochondrial gene introgression was detected between Q. spinosa and Q. shini in Lushan, Jiangxi Province. Our study also detected ancient hybridizations between a female Q. spinosa and a male Q. jiulongensis in Zhejiang Province. Analyses on mitochondrial and nuclear genes verified three candidate cryptic species in Q. spinosa, and a cryptic species may also exist in Q. boulengeri. However, no evidence of introgressive hybridization was found between Q. spinosa and Q. boulengeri. Quasipaa exilispinosa from all the sampling localities appeared to be deeply divergent from other communities. Our results suggest widespread introgressive hybridization in closely related species of Quasipaa and provide a fundamental basis for illumination of the forming mechanism of introgressive hybridization, classification of species, and biodiversity assessment in Quasipaa.     

Complete mitochondrial genomes of prasinophyte algae Pyramimonas parkeae and Cymbomonas tetramitiformis

Mitochondria are archetypal eukaryotic organelles that were acquired by endosymbiosis of an ancient species of alpha‐proteobacteria by the last eukaryotic common ancestor. The genetic information contained within the mitochondrial genome has been an important source of information for resolving relationships among eukaryotic taxa. In this study, we utilized mitochondrial and chloroplast genomes to explore relationships among prasinophytes. Prasinophytes are represented by diverse early‐diverging green algae whose physical structures and genomes have the potential to elucidate the traits of the last common ancestor of the Viridiplantae (or Chloroplastida). We constructed de novo mitochondrial genomes for two prasinophyte algal species, Pyramimonas parkeae and Cymbomonas tetramitiformis, representing the prasinophyte clade. Comparisons of genome structure and gene order between these species and to those of other prasinophytes revealed that the mitochondrial genomes of P. parkeae and C. tetramitiformis are more similar to each other than to other prasinophytes, consistent with other molecular inferences of the close relationship between these two species. Phylogenetic analyses using the inferred amino acid sequences of mitochondrial and chloroplast protein‐coding genes resolved a clade consisting of P. parkeae and C. tetramitiformis; and this group (representing the prasinophyte clade I) branched with the clade II, consistent with previous studies based on the use of nuclear gene markers.     

The complete chloroplast and mitochondrial genomes of the diatom Nitzschia palea (Bacillariophyceae) demonstrate high sequence similarity to the endosymbiont organelles of the dinotom Durinskia baltica

Nitzschia palea is a common freshwater diatom used as a bioindicator because of its tolerance of polluted waterways. There is also evidence it may be the tertiary endosymbiont within the “dinotom” dinoflagellate Durinskia baltica. A putative strain of N. palea was collected from a pond on the University of Virginia's College at Wise campus and cultured. For initial identification, three markers were sequenced—nuclear 18S rDNA, the chloroplast 23S rDNA, and rbcL. Morphological characteristics were determined using light and scanning electron microscopy; based on these observations the cells were identified as N. palea and named strain “Wise.” DNA from N. palea was deep sequenced and the chloroplast and mitochondrial genomes assembled. Single gene phylogenies grouped N. palea—Wise within a clearly defined N. palea clade and showed it was most closely related to the strain “SpainA3.” The chloroplast genome of N. palea is 119,447 bp with a quadripartite structure, 135 protein‐coding, 28 tRNA, and 3 rRNA genes. The mitochondrial genome is 37,754 bp with a single repeat region as found in other diatom chondriomes, 37 protein‐coding, 23 tRNA, and 2 rRNA genes. The chloroplast genomes of N. palea and D. baltica have identical gene content, synteny, and a 92.7% pair‐wise sequence similarity with most differences occurring in intergenic regions. The N. palea mitochondrial genome and D. baltica's endosymbiont mitochondrial genome also have identical gene content and order with a sequence similarity of 90.7%. Genome‐based phylogenies demonstrated that D. baltica is more similar to N. palea than any other diatom sequence currently available. These data provide the genome sequences of two organelles for a widespread diatom and show they are very similar to those of Durinskia baltica's endosymbiont.     

The first complete mitochondrial genome sequence of Nanorana parkeri and Nanorana ventripunctata (Amphibia: Anura: Dicroglossidae), with related phylogenetic analyses

Members of the Nanorana genus (family Dicroglossidae) are often referred to as excellent model species with which to study amphibian adaptations to extreme environments and also as excellent keystone taxa for providing insights into the evolution of the Dicroglossidae. However, a complete mitochondrial genome is currently only available for Nanorana pleskei. Thus, we analyzed the complete mitochondrial genomes of Nanorana parkeri and Nanorana ventripunctata to investigate their evolutionary relationships within Nanorana and their phylogenetic position in the family Dicroglossidae. Our results showed that the genomes of N. parkeri (17,837 bp) and N. ventripunctata (18,373 bp) encode 13 protein‐coding genes (PCGs), two ribosomal RNA genes, 23 transfer RNA (tRNA) genes, and a noncoding control region. Overall sequences and genome structure of the two species showed high degree of similarity with N. pleskei, although the motif structures and repeat sequences of the putative control region showed clear differences among these three Nanorana species. In addition, a tandem repeat of the tRNA‐Met gene was found located between the tRNA‐Gln and ND2 genes. On both the 5′ and 3′‐sides, the control region possessed distinct repeat regions; however, the CSB‐2 motif was not found in N. pleskei. Based on the nucleotide sequences of 13 PCGs, our phylogenetic analyses, using Bayesian inference and maximum‐likelihood methods, illustrate the taxonomic status of Nanorana with robust support showing that N. ventripunctata and N. pleskei are more closely related than they are to N. parkeri. In conclusion, our analyses provide a more robust and reliable perspective on the evolutionary history of Dicroglossidae than earlier analyses, which used only a single species (N. pleskei).     

Testing hypotheses for genealogical discordance in a rainforest lizard

Genealogical discordance, or when different genes tell distinct stories although they evolved under a shared history, often emerges from either coalescent stochasticity or introgression. In this study, we present a strong case of mito‐nuclear genealogical discordance in the Australian rainforest lizard species complex of Saproscincus basiliscus and S. lewisi. One of the lineages that comprises this complex, the Southern S. basiliscus lineage, is deeply divergent at the mitochondrial genome but shows markedly less divergence at the nuclear genome. By placing our results in a comparative context and reconstructing the lineages' demography via multilocus and coalescent‐based approximate Bayesian computation methods, we test hypotheses for how coalescent variance and introgression contribute to this pattern. These analyses suggest that the observed genealogical discordance likely results from introgression. Further, to generate such strong discordance, introgression probably acted in concert with other factors promoting asymmetric gene flow between the mitochondrial and nuclear genomes, such as selection or sex‐biased dispersal. This study offers a framework for testing sources of genealogical discordance and suggests that historical introgression can be an important force shaping the genetic diversity of species and their populations.     

Presence of nuclear copies of mitochondrial origin (NUMTs) in two related species of stingless bee genus Melipona (Hymenoptera: Meliponini)

Transferred copies of mitochondrial DNA (mtDNA) into the nuclear genome (numts) have been reported in several Hymenoptera species, even at a high density in the honey bee nuclear genome. The accidental amplification of numts in phylogenetic studies focused on mtDNA highlights the importance of a correct determination of numts and their related mtDNA sequences. We report here the presence of numts derived from a mitochondrial rDNA 16S gene in the genome of the stingless bee species Melipona colimana and M. fasciata (tribe Meliponini) from Western Mexico. PCR products were cloned in both species obtaining thirty paralogous numts. Numts were identified by the presence of insertions and deletions and the disruption of the 16S secondary structure. Further phylogenetic analyses including alternative mitochondrial cox1 and nuclear ITS1 genes have revealed the presence of another numt (cox1) in the nuclear genome of these two species, and place both as sister lineages within the subgenus Michmelia. This is one of the first studies reporting the presence of numts in Meliponini species, and supports previous studies suggesting frequent transfer of mtDNA to the nuclear genome in Hymenoptera.     

Mitogenomics reveals high synteny and long evolutionary histories of sympatric cryptic nematode species

Species with seemingly identical morphology but with distinct genetic differences are abundant in the marine environment and frequently co‐occur in the same habitat. Such cryptic species are typically delineated using a limited number of mitochondrial and/or nuclear marker genes, which do not yield information on gene order and gene content of the genomes under consideration. We used next‐generation sequencing to study the composition of the mitochondrial genomes of four sympatrically distributed cryptic species of the Litoditis marina species complex (PmI, PmII, PmIII, and PmIV). The ecology, biology, and natural occurrence of these four species are well known, but the evolutionary processes behind this cryptic speciation remain largely unknown. The gene order of the mitochondrial genomes of the four species was conserved, but differences in genome length, gene length, and codon usage were observed. The atp8 gene was lacking in all four species. Phylogenetic analyses confirm that PmI and PmIV are sister species and that PmIII diverged earliest. The most recent common ancestor of the four cryptic species was estimated to have diverged 16 MYA. Synonymous mutations outnumbered nonsynonymous changes in all protein‐encoding genes, with the Complex IV genes (coxI‐III) experiencing the strongest purifying selection. Our mitogenomic results show that morphologically similar species can have long evolutionary histories and that PmIII has several differences in genetic makeup compared to the three other species, which may explain why it is better adapted to higher temperatures than the other species.     

Phylogenetics and comparative analysis of the mitochondrial genomes of three violet‐ringed octopuses

Similar morphological characters and little molecular data of Amphioctopus rex, A. neglectus and A. cf. ovulum resulted in their unknown phylogenetic statuses and equivocal relationships. In this study, the complete mitochondrial genomes of these three species collected in Chinese waters were sequenced and compared with each other to clarify the relationships among them. The lengths of the mitochondrial genomes varied from 15,646 bp to 15,814 bp, and the A + T content and GC skew for protein‐coding genes showed little variation. In contrast, both a dendrogram based on codon usage and the gene arrangements of the three octopuses showed that A. rex was more closely related to A. neglectus than to A. cf. ovulum. Five data sets and two methods (maximum likelihood and Bayesian inference) were utilized for the first time to explore the phylogenetic relationships among these three species in Octopodidae. The results indicated that a data set combining protein‐coding genes and RNA genes (PR) was optimal for analysing the relationships among 43 cephalopods. All of the phylogenetic trees divided the cephalopods into 10 taxa and supported the monophyly of Oegopsida, Myopsida, Sepiidae and Octopodidae. In this study, Idiosepiidae was classified as sister to Sepiolidae. Trees constructed using all data sets robustly supported the monophyly of the genus Amphioctopus. Notably, A. rex was more closely related to A. neglectus than to A. cf. ovulum, although these three species share the characteristic of violet rings on dark ocelli.     

Phylogeography and lineage‐specific patterns of genetic diversity and molecular evolution in a group of North American skinks

Geography influences the evolutionary trajectory of species by mediating opportunities for hybridization, gene flow, demographic shifts and adaptation. We sought to understand how geography and introgression can generate species‐specific patterns of genetic diversity by examining phylogeographical relationships in the North American skink species Plestiodon multivirgatus and P. tetragrammus (Squamata: Scincidae). Using a multilocus dataset (three mitochondrial genes, four nuclear genes; a total of 3455 bp) we discovered mito‐nuclear discordance, consistent with mtDNA introgression. We further tested for evidence of species‐wide mtDNA introgression by using comparisons of genetic diversity, selection tests and extended Bayesian skyline analyses. Our findings suggest that P. multivirgatus acquired its mitochondrial genome from P. tetragrammus after their initial divergence. This putative species‐wide mitochondrial capture was further evidenced by statistically indistinguishable substitution rates between mtDNA and nDNA in P. multivirgatus. This rate discrepancy was observed in P. multivirgatus but not P. tetragrammus, which has important implications for studies that combine mtDNA and nDNA sequences when inferring time since divergence between taxa. Our findings suggest that by facilitating opportunities for interspecific introgression, geography can alter the course of molecular evolution between recently diverged lineages.     

The Mitochondrial Genomes of Two Scallops, Argopecten irradians and Chlamys farreri (Mollusca: Bivalvia): The Most Highly Rearranged Gene Order in the Family Pectinidae

Molluscs in general and bivalves in particular, exhibit an extraordinary degree of mitochondrial gene order variation when compared with other metazoans. Here, we determined the mitochondrial genomes of two scallops Argopecten irradians and Chlamys farreri. The complete mitochondrial genome of A. irradians is 16,211 nts in length and the nearly complete mitochondrial genome of C. farreri is 20,789 nts in length. Both of the genomes contain 35 genes including 12 protein-coding genes, 2 ribosomal RNAs, and 21 transfer RNAs. In contrast to the typical animal mitochondrial genome, both of them lack one protein-coding gene atp8 and two trnSs, but show an additional copy of trnF in A. irradians and of trnM in C. farreri, respectively. Gene order and genome content were compared among the four sequenced scallops. Gene arrangement of C. farreri closely resembles that of Mizuhopecten yessoensis. However, two genomes of C. farreri and A. irradians show only three small identical gene blocks and two genomes of A. irradians and Placopecten magellanicus share only one gene block. Comparison of the gene arrangement demonstrated that the four scallops share few identical gene blocks although they belong to the same family. This feature is seldom observed in Metazoa, even in other molluscan classes. The dramatic gene rearrangement often occurs in bivalves, especially in marine bivalves. In addition, comparisons of genomic character among bivalves are also presented.     

Resequencing and comparison of whole mitochondrial genome to gain insight into the evolutionary status of the Shennongjia golden snub‐nosed monkey (SNJ R. roxellana)

Shennongjia Rhinopithecus roxellana (SNJ R. roxellana) is the smallest geographical population of R. roxellana. The phylogenetic relationships among its genera and species and the biogeographic processes leading to their current distribution are largely unclear. To address these issues, we resequenced and obtained a new, complete mitochondrial genome of SNJ R. roxellana by next‐generation sequencing and standard Sanger sequencing. We analyzed the gene composition, constructed a phylogenetic tree, inferred the divergence ages based on complete mitochondrial genome sequences, and analyzed the genetic divergence of 13 functional mtDNA genes. The phylogenetic tree and divergence ages showed that R. avunculus (the Tonkin snub‐nosed monkey) was the first to diverge from the Rhinopithecus genus ca. 2.47 million years ago (Ma). Rhinopithecus bieti and Rhinopithecus strykeri formed sister groups, and the second divergence from the Rhinopithecus genus occurred ca. 1.90 Ma. R. roxellana and R. brelichi diverged from the Rhinopithecus genus third, ca. 1.57 Ma. SNJ R. roxellana was the last to diverge within R. roxellana species in 0.08 Ma, and the most recent common ancestor of R. roxellana is 0.10 Ma. The analyses on gene composition showed SNJ R. roxellana was the newest geographic population of R. roxellana. The work will help to develop a more accurate protection policy for SNJ R. roxellana and facilitate further research on selection and adaptation of R. roxellana.     

Phylogenetic positioning of the Antarctic alga Prasiola crispa (Trebouxiophyceae) using organellar genomes and their structural analysis

Antarctica is one of the most difficult habitats for sustaining life on earth; organisms that live there have developed different strategies for survival. Among these organisms is the green alga Prasiola crispa, belonging to the class Trebouxiophyceae. The literature on P. crispa taxonomy is scarce, and many gaps in the evolutionary relationship with its closest relatives remain. The goal of this study was to analyze the evolutionary relationships between P. crispa and other green algae using plastid and mitochondrial genomes. In addition, we analyzed the synteny conservation of these genomes of P. crispa with those of closely related species. Based on the plastid genome, P. crispa grouped with Prasiolopsis sp. SAG 84.81, another Trebouxiophyceaen species from the Prasiola clade. Based on the mitochondrial genome analysis, P. crispa grouped with other Trebouxiophyceaen species but had a basal position. The structure of the P. crispa chloroplast genome had low synteny with Prasiolopsis sp. SAG 84.81, despite some conserved gene blocks. The same was observed in the mitochondrial genome compared with Coccomyxa subellipsoidea C‐169. We were able to establish the phylogenetic position of P. crispa with other species of Trebouxiophyceae using its genomes. In addition, we described the plasticity of these genomes using a structural analysis. The plastid and mitochondrial genomes of P. crispa will be useful for further genetic studies, phylogenetic analysis and resource protection of P. crispa as well as for further phylogenetic analysis of Trebouxiophyceaen green algae.     

Mitogenomic phylogeny of Nassarius (Gastropoda: Neogastropoda)

Nassariids (Family Nassariidae) are a group of marine snails that are distributed worldwide, with their maximum species diversity in tropical regions, particularly the Indo‐Pacific. However, the traditional taxonomy of Nassariidae defined by shell or radula characters is usually inconsistent with little phylogenetic signal. In the present study, the complete mitochondrial (mt) genomes of nine Nassarius species were sequenced and compared with other eight nassariid species previously reported. All nassariid mt genomes showed the same gene order as in most caenogastropods and shared a very similar pattern with respect to genome size, nucleotide composition and AT contents. A deletion of three nucleotides in nad6 gene was detected in Nassarius jacksonianus and Nassarius acuticostus, and this feature also provided implications for nassariid phylogeny. The genetic distance analysis and reconstructed phylogeny revealed a distant relationship between N. jacksonianus or N. acuticostus and other members in Nassarius. The mitogenomic phylogeny recovered the evolutionary relationships within Nassarius with high statistical support. In addition, a chronogram was reconstructed under an uncorrelated relaxed molecular clock, which dated the divergence among main lineages of Nassarius during ~31 MYA.     

Complete mitochondrial genome sequence of the Himalayan Griffon,Gyps himalayensis (Accipitriformes: Accipitridae): Sequence,structure, and phylogenetic analyses

This is the first study to describe the mitochondrial genome of the Himalayan Griffon, Gyps himalayensis, which is an Old World vulture belonging to the family Accipitridae and occurring along the Himalayas and the adjoining Tibetan Plateau. Its mitogenome is a closed circular molecule 17,381 bp in size containing 13 protein‐coding genes, 22 tRNA coding genes, two rRNA‐coding genes, a control region (CR), and an extra pseudo‐control region (CCR) that are conserved in most Accipitridae mitogenomes. The overall base composition of the G. himalayensis mitogenome is 24.55% A, 29.49% T, 31.59% C, and 14.37% G, which is typical for bird mitochondrial genomes. The alignment of the Accipitridae species control regions showed high levels of genetic variation and abundant AT content. At the 5′ end of the domain I region, a long continuous poly‐C sequence was found. Two tandem repeats were found in the pseudo‐control regions. Phylogenetic analysis with Bayesian inference and maximum likelihood based on 13 protein‐coding genes indicated that the relationships at the family level were (Falconidae + (Cathartidae + (Sagittariidae + (Accipitridae + Pandionidae))). In the Accipitridae clade, G. himalayensis is more closely related to Aegypius monachus than to Spilornis cheela. The complete mitogenome of G. himalayensis provides a potentially useful resource for further exploration of the taxonomic status and phylogenetic history of Gyps species.     

The Chimonanthus salicifolius genome provides insight into magnoliid evolution and flavonoid biosynthesis

Chimonanthus salicifolius, a member of the Calycanthaceae of magnoliids, is one of the most famous medicinal plants in Eastern China. Here, we report a chromosome‐level genome assembly of C. salicifolius, comprising 820.1 Mb of genomic sequence with a contig N50 of 2.3 Mb and containing 36 651 annotated protein‐coding genes. Phylogenetic analyses revealed that magnoliids were sister to the eudicots. Two rounds of ancient whole‐genome duplication were inferred in the C. salicifolious genome. One is shared by Calycanthaceae after its divergence with Lauraceae, and the other is in the ancestry of Magnoliales and Laurales. Notably, long genes with > 20 kb in length were much more prevalent in the magnoliid genomes compared with other angiosperms, which could be caused by the length expansion of introns inserted by transposon elements. Homologous genes within the flavonoid pathway for C. salicifolius were identified, and correlation of the gene expression and the contents of flavonoid metabolites revealed potential critical genes involved in flavonoids biosynthesis. This study not only provides an additional whole‐genome sequence from the magnoliids, but also opens the door to functional genomic research and molecular breeding of C. salicifolius.     

The complete nucleotide sequence of the mitochondrial genome of <Emphasis Type="Italic">Phthonandria atrilineata</Emphasis> (Lepidoptera: Geometridae)

Using long-polymerase chain reaction (Long-PCR) method, we determined the complete nucleotide sequence of the mitochondrial genome (mitogenome) of Phthonandria atrilineata. The complete mtDNA from P. atrilineata was 15,499 base pairs in length and contained 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes, and a control region. The P. atrilineata genes were in the same order and orientation as the completely sequenced mitogenomes of other lepidopteran species. The nucleotide composition of P. atrilineata mitogenome was biased toward A + T nucleotides (81.02%), and the 13 PCGs show different A + T contents that range from 73.25% (cox1) to 92.12% (atp8). Phthonandria had the canonical set of 22 tRNA genes, that fold in the typical cloverleaf structure described for metazoan mt tRNAs, with the unique exception of trnS(AGN). The phylogenetic relationships were reconstructed with the concatenated sequences of the 13 PCGs of the mitochondrial genome, which confirmed that P. atrilineata is most closely related to the superfamily Bombycoidea.     

Phylogeographic characteristics of hypervariable regions in the mitochondrial genome of a cosmopolitan,bloom‐forming raphidophyte,Heterosigma akashiwo

The global distribution of phytoplankton is defined by many events, including long‐term evolutionary processes and shorter time span processes (e.g., global climate change). Furthermore, human‐assisted, unintentional dispersion, including the transport of live fish and spat for aquaculture, and transfer of aquatic microorganisms contained in ship ballast water, may aid the spread of phytoplankton. To understand the phylogeographic history of a species, the development of useful molecular markers is crucial. We previously reported a hypervariable mitochondrial gene in the cosmopolitan bloom‐forming alga, Heterosigma akashiwo. In this study, we identified two additional hypervariable segments in the H. akashiwo mitochondrial genome, one a protein coding sequence, and the other an intergenic region, by comparing the whole mitogenomes of strains obtained from various geographic origins. Interestingly, the newly identified hypervariable protein coding sequence was a paralogue of the previously identified gene, and both sequences showed tendencies to correlate with latitude of geographic origin. However, the hypervariable intergenic sequence did not show a clear correlation with origin. Our results demonstrated that the protein coding sequences may serve as useful tools for understanding the phylogeographic history of H. akashiwo, and they may crucially function in adaptation of the species to the environment.