Serial replacement of a diatom endosymbiont in the marine dinoflagellate Peridinium quinquecorne (Peridiniales, Dinophyceae)

To infer the phylogeny of both the host and the endosymbiont of Peridinium quinquecorne Abé, the small subunit (SSU) ribosomal DNA (rDNA) from the host and two genes of endosymbiont origin (plastid‐encoded rbcL and nuclear‐encoded SSU rDNA) were determined. The phylogenetic analysis of the host revealed that the marine dinoflagellate P. quinquecorne formed a clade with other diatom‐harbouring dinoflagellates, including Kryptoperidinium foliaceum (Stein) Lindeman, Durinskia baltica (Levander) Carty et Cox and Galeidinium rugatum Tamura et Horiguchi, indicating a single endosymbiotic event for this lineage. Phylogenetic analyses of the endosymbiont in these organisms revealed that the endosymbiont of P. quinquecorne formed a clade with a centric diatom (SSU data indicated it to be closely related to Chaetoceros), whereas the endosymbionts of other three dinoflagellates formed a clade with a pennate diatom. The discrepancy between the host and the endosymbiont phylogenies suggests a secondary replacement of the endosymbiont from a pennate to a centric diatom in P. quinquecorne.     

The dinoflagellates <Emphasis Type="Italic">Durinskia baltica</Emphasis> and <Emphasis Type="Italic">Kryptoperidinium foliaceum</Emphasis> retain functionally overlapping mitochondria from two evolutionarily distinct lineages

Background  

The dinoflagellates Durinskia baltica and Kryptoperidinium foliaceum are distinguished by the presence of a tertiary plastid derived from a diatom endosymbiont. The diatom is fully integrated with the host cell cycle and is so altered in structure as to be difficult to recognize it as a diatom, and yet it retains a number of features normally lost in tertiary and secondary endosymbionts, most notably mitochondria. The dinoflagellate host is also reported to retain mitochondrion-like structures, making these cells unique in retaining two evolutionarily distinct mitochondria. This redundancy raises the question of whether the organelles share any functions in common or have distributed functions between them.     

The complete chloroplast genome of Gracilariopsis lemaneiformis (Rhodophyta) gives new insight into the evolution of family Gracilariaceae

The complete chloroplast genome of Gracilariopsis lemaneiformis was recovered from a Next Generation Sequencing data set. Without quadripartite structure, this chloroplast genome (183,013 bp, 27.40% GC content) contains 202 protein‐coding genes, 34 tRNA genes, 3 rRNA genes, and 1 tmRNA gene. Synteny analysis showed plasmid incorporation regions in chloroplast genomes of three species of family Gracilariaceae and in Grateloupia taiwanensis of family Halymeniaceae. Combined with reported red algal plasmid sequences in nuclear and mitochondrial genomes, we postulated that red algal plasmids may have played an important role in ancient horizontal gene transfer among nuclear, chloroplast, and mitochondrial genomes. Substitution rate analysis showed that purifying selective forces maintaining stability of protein‐coding genes of nine red algal chloroplast genomes over long periods must be strong and that the forces acting on gene groups and single genes of nine red algal chloroplast genomes were similar and consistent. The divergence of Gp. lemaneiformis occurred ~447.98 million years ago (Mya), close to the divergence time of genus Pyropia and Porphyra (443.62 Mya).     

Sterols of the ‘dinotom’ Durinskia baltica (Dinophyceae) are of dinoflagellate origin

‘Dinotoms’ are a relatively small group of dinoflagellates with aberrant tertiary plastids of diatom origin, thus differing from the majority of photosynthetic dinoflagellates which possess the carotenoid pigment peridinin and have secondary plastids of red algal origin. As part of our laboratory's continuing efforts to examine such unusual dinoflagellates in the search for clues to the evolution of their lipid compositions, we have examined the sterol composition of the dinotom Durinskia baltica. As such, we here compared its sterols to those of the previously examined dinotom, Kryptoperidinium foliaceum, more broadly to other photosynthetic, peridinin-containing dinoflagellates, and to the diatom genus Nitzschia, which is the presumed ancestor of the D. baltica dinotom plastid. Sterols are ringed lipids, common to eukaryotes, thought to reinforce phospholipid bilayers. Many peridinin-containing dinoflagellates have sterol compositions which are enriched by the presence of cholesterol (cholest-5-en-3β-ol) and 4α-methyl-substituted sterols such as dinosterol (4α,23,24-trimethyl-5α-cholest-22E-en-3β-ol); this has also been found to be true for K. foliaceum despite its aberrant plastid ancestry. Our objective was to determine if this is also true for D. baltica as only the second dinotom to have its sterols characterized in detail, and to determine if there is any indication of prominent sterols which are uncommon to dinoflagellates, possibly originating from the diatom endosymbiont, as has been demonstrated previously with K. foliaceum and D. baltica chloroplast-associated galactolipids of clear diatom origin. Our results demonstrate that like K. foliaceum, the major sterols of D. baltica are cholesterol, dinosterol, and other 4α-methyl-substituted sterols common to dinoflagellates. Although there were a number of minor sterols, none were found with obvious origin from the diatom endosymbiont, indicating that most originated with the dinoflagellate host itself, most likely before acquisition of the diatom tertiary plastid.     

GALEIDINIIUM RUGATUM GEN. ET SP. NOV. (DINOPHYCEAE), A NEW COCCOID DINOFLAGELLATE WITH A DIATOM ENDOSYMBIONT1

A new sand‐dwelling dinoflagellate from Palau, Galeidinium rugatum Tamura et Horiguchi gen. et sp. nov., is described. The life cycle of this new alga consists of a dominant nonmotile phase and a brief motile phase. The motile cell transforms itself directly into the nonmotile cell after swimming for a short period, and cell division takes place in the nonmotile phase. The nonmotile cell possesses a dome‐like cell covering, which is wrinkled and equipped with a transverse groove on the surface. The cell has 10–20 chloroplasts and a distinct eyespot. The motile cell is Gymnodinium‐like in shape. The dinoflagellate possesses an endosymbiotic alga to which the chloroplasts belong and which is separated from the host (dinoflagellate) cytoplasm by a unit membrane. The endosymbiont cytoplasm also possesses its own eukaryotic nucleus and mitochondria. The eyespot is surrounded by triple membranes and is located in the host cytoplasm. Photosynthetic pigment analysis, using HPLC, revealed that G. rugatum possesses fucoxanthin as the principal accessory pigment instead of peridinin. The rbcL tree showed that G. rugatum is monophyletic with Durinskia baltica (Levander) Carty et Cox and Kryptoperidinium foliaceum (Stein) Lindemann and that this clade is closely related to the pennate diatom, Cylindrotheca sp. The endosymbiont of G. rugatum is therefore shown to be a diatom. Phylogenetic analysis based on small subunit rDNA sequences demonstrated that G. rugatum, D. baltica, and K. foliaceum, all of which are known to harbor an endosymbiont of diatom origin, are closely related.     

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).     

Cyst‐Theca Relationship and Phylogenetic Position of Impagidinium caspienense Incubated from Caspian Sea Surface Sediments: Relation to Gonyaulax baltica and Evidence for Heterospory within Gonyaulacoid Dinoflagellates

We investigate the cyst‐theca relationship of Impagidinium caspienense. Through an incubation experiment, we succeeded in examining the motile stage. Additional molecular analysis of single‐cyst PCR (LSU and SSU rDNA) reveal that the cyst is related to the species Gonyaulax baltica Ellegaard et al. ( 2002 ). The ability of this species to belong to two types of cyst‐based genera (spiniferate and impagidinioid) suggests that environmental (particularly salinity) and not genetic factors explain the formation of both morphotypes by G. baltica, which provides evidence for heterospory in this species. The affiliation to G. baltica demonstrates that I. caspienense is not endemic to the Caspian Sea. The phylogenetic position of several other gonyaulacoid species is also documented: Impagidinium pallidum, Ataxiodinium choane, Pyxidinopsis psilata, Spiniferites belerius, and Spiniferites ramosus. The LSU and SSU rDNA based phylogenies suggest that the genera Impagidinium and Spiniferites are not monophyletic, and that P. psilata and A. choane are close to Gonyaulax verior and Gonyaulax polygramma, respectively. In addition, this study accentuates the importance of cyst morphology in the classification of the Gonyaulacales.     

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 first complete mitochondrial genome of the Indian Tent Turtle,Pangshura tentoria (Testudines: Geoemydidae): Characterization and comparative analysis

The characterization of a complete mitogenome is widely used in genomics studies for systematics and evolutionary research. However, the sequences and structural motifs contained within the mitogenome of Testudines taxa have rarely been examined. The present study decodes the first complete mitochondrial genome of the Indian Tent Turtle, Pangshura tentoria (16,657 bp) by using next‐generation sequencing. This denovo assembly encodes 37 genes: 13 protein‐coding genes (PCGs), 22 transfer RNA (tRNAs), two ribosomal RNA, and one control region (CR). Most of the genes were encoded on majority strand, except for one PCG (NADH dehydrogenase subunit 6) and eight tRNAs. Most of the PCGs were started with an ATG initiation codon, except for Cytochrome oxidase subunit 1 with “GTG” and NADH dehydrogenase subunit 5 with “ATA.” The termination codons, “TAA” and “AGA” were observed in two subunits of NADH dehydrogenase gene. The relative synonymous codon usage analysis revealed the maximum abundance of alanine, isoleucine, leucine, and threonine. The nonsynonymous/synonymous ratios were <1 in all PCGs, which indicates strong negative selection among all Geoemydid species. The study also found the typical cloverleaf secondary structure in most of the tRNA genes, except for serine with the lack of the conventional DHU arm. The comparative study of Geoemydid mitogenomes revealed the occurrence of tandem repeats was frequent in the 3′ end of CR. Further, two copies of a unique tandem repeat “TTCTCTTT” were identified in P. tentoria. The Bayesian and maximum‐likelihood phylogenetic trees using concatenation of 13 PCGs revealed the close relationships of P. tentoria with Batagur trivittata in the studied dataset. All the Geoemydid species showed distinct clustering with high bootstrap support congruent with previous evolutionary hypotheses. We suggest that the generations of more mitogenomes of Geoemydid species are required, to improve our understanding of their in‐depth phylogenetic and evolutionary relationships.     

Fates of Evolutionarily Distinct,Plastid‐type Glyceraldehyde 3‐phosphate Dehydrogenase Genes in Kareniacean Dinoflagellates

The ancestral kareniacean dinoflagellate has undergone tertiary endosymbiosis, in which the original plastid is replaced by a haptophyte endosymbiont. During this plastid replacement, the endosymbiont genes were most likely flowed into the host dinoflagellate genome (endosymbiotic gene transfer or EGT). Such EGT may have generated the redundancy of functionally homologous genes in the host genome—one has resided in the host genome prior to the haptophyte endosymbiosis, while the other transferred from the endosymbiont genome. However, it remains to be well understood how evolutionarily distinct but functionally homologous genes were dealt in the dinoflagellate genomes bearing haptophyte‐derived plastids. To model the gene evolution after EGT in plastid replacement, we here compared the characteristics of the two evolutionally distinct genes encoding plastid‐type glyceraldehyde 3‐phosphate dehydrogenase (GAPDH) in Karenia brevis and K. mikimotoi bearing haptophyte‐derived tertiary plastids: “gapC1h” acquired from the haptophyte endosymbiont and “gapC1p” inherited from the ancestral dinoflagellate. Our experiments consistently and clearly demonstrated that, in the two species examined, the principal plastid‐type GAPDH is encoded by gapC1h rather than gapC1p. We here propose an evolutionary scheme resolving the EGT‐derived redundancy of genes involved in plastid function and maintenance in the nuclear genomes of dinoflagellates that have undergone plastid replacements. Although K. brevis and K. mikimotoi are closely related to each other, the statuses of the two evolutionarily distinct gapC1 genes in the two Karenia species correspond to different steps in the proposed scheme.     

Flip‐flop organization in the chloroplast genome of Capsosiphon fulvescens (Ulvophyceae,Chlorophyta)

To better understand organelle genome evolution of the ulvophycean green alga Capsosiphon fulvescens, we sequenced and characterized its complete chloroplast genome. The circular chloroplast genome was 111,561 bp in length with 31.3% GC content that contained 108 genes including 77 protein‐coding genes, two copies of rRNA operons, and 27 tRNAs. In this analysis, we found the two types of isoform, called heteroplasmy, were likely caused by a flip‐flop organization. The flip‐flop mechanism may have caused structural variation and gene conversion in the chloroplast genome of C. fulvescens. In a phylogenetic analysis based on all available ulvophycean chloroplast genome data, including a new C. fulvescens genome, we found three major conflicting signals for C. fulvescens and its sister taxon Pseudoneochloris marina within 70 individual genes: (i) monophyly with Ulotrichales, (ii) monophyly with Ulvales, and (iii) monophyly with the clade of Ulotrichales and Ulvales. Although the 70‐gene concatenated phylogeny supported monophyly with Ulvales for both species, these complex phylogenetic signals of individual genes need further investigations using a data‐rich approach (i.e., organelle genome data) from broader taxon sampling.     

The genome sequence of Samia ricini,a new model species of lepidopteran insect

Samia ricini, a gigantic saturniid moth, has the potential to be a novel lepidopteran model species. Samia ricini is far more resistant to diseases than the current model species Bombyx mori, and therefore can be more easily reared. In addition, genetic resources available for S. ricini rival those for B. mori: at least 26 ecoraces of S. ricini are reported and S. ricini can hybridize with wild Samia species, which are distributed throughout Asian countries, and produce fertile progenies. Physiological traits such as food preference, integument colour and larval spot pattern differ among S. ricini strains and wild Samia species so that those traits can be targeted in forward genetic analyses. To facilitate genetic research in S. ricini, we determined its whole genome sequence. The assembled genome of S. ricini was 458 Mb with 155 scaffolds, and the scaffold N50 length of the assembly was ~ 21 Mb. In total, 16,702 protein coding genes were predicted. While the S. ricini genome was mostly collinear with that of B. mori with some rearrangements and few S. ricini‐specific genes were discovered, chorion genes and fibroin genes seemed to have expanded in the S. ricini lineage. As the first step of genetic analyses, causal genes for “Blue,” “Yellow,” “Spot,” and “Red cocoon” phenotypes were mapped to chromosomes.     

Are glacial refugia hotspots of speciation and cytonuclear discordances? Answers from the genomic phylogeography of Spanish common frogs

Subdivided Pleistocene glacial refugia, best known as “refugia within refugia”, provided opportunities for diverging populations to evolve into incipient species and/or to hybridize and merge following range shifts tracking the climatic fluctuations, potentially promoting extensive cytonuclear discordances and “ghost” mtDNA lineages. Here, we tested which of these opposing evolutionary outcomes prevails in northern Iberian areas hosting multiple historical refugia of common frogs (Rana cf. temporaria), based on a genomic phylogeography approach (mtDNA barcoding and RAD‐sequencing). We found evidence for both incipient speciation events and massive cytonuclear discordances. On the one hand, populations from northwestern Spain (Galicia and Asturias, assigned to the regional endemic R. parvipalmata), are deeply‐diverged at mitochondrial and nuclear genomes (~4 My of independent evolution), and barely admix with northeastern populations (assigned to R. temporaria sensu stricto) across a narrow hybrid zone (~25 km) located in the Cantabrian Mountains, suggesting that they represent distinct species. On the other hand, the most divergent mtDNA clade, widespread in Cantabria and the Basque country, shares its nuclear genome with other R. temporaria s. s. lineages. Patterns of population expansions and isolation‐by‐distance among these populations are consistent with past mitochondrial capture and/or drift in generating and maintaining this ghost mitochondrial lineage. This remarkable case study emphasizes the complex evolutionary history that shaped the present genetic diversity of refugial populations, and stresses the need to revisit their phylogeography by genomic approaches, in order to make informed taxonomic inferences.     

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.     

Development of genomic resources for Nothofagus species using next‐generation sequencing data

Using next‐generation sequencing, we developed the first whole‐genome resources for two hybridizing Nothofagus species of the Patagonian forests that crucially lack genomic data, despite their ecological and industrial value. A de novo assembly strategy combining base quality control and optimization of the putative chloroplast gene map yielded ~32 000 contigs from 43% of the reads produced. With 12.5% of assembled reads, we covered ~96% of the chloroplast genome and ~70% of the mitochondrial gene content, providing functional and structural annotations for 112 and 52 genes, respectively. Functional annotation was possible on 15% of the contigs, with ~1750 potentially novel nuclear genes identified for Nothofagus species. We estimated that the new resources (13.41 Mb in total) included ~4000 gene regions representing ~6.5% of the expected genic partition of the genome, the remaining contigs potentially being nongenic DNA. A high‐quality single nucleotide polymorphisms resource was developed by comparing various filtering methods, and preliminary results indicate a strong conservation of cpDNA genomes in contrast to numerous exclusive nuclear polymorphisms in both species. Finally, we characterized 2274 potential simple sequence repeat (SSR) loci, designed primers for 769 of them and validated nine of 29 loci in 42 individuals per species. Nothofagus obliqua had more alleles (4.89) on average than N. nervosa (2.89), 8 SSRs were efficient to discriminate species, and three were successfully transferred in three other Nothofagus species. These resources will greatly help for future inferences of demographic, adaptive and hybridizing events in Nothofagus species, and for conserving and managing natural populations.