Widespread occurrence of organelle genome-encoded 5S rRNAs including permuted molecules

5S Ribosomal RNA (5S rRNA) is a universal component of ribosomes, and the corresponding gene is easily identified in archaeal, bacterial and nuclear genome sequences. However, organelle gene homologs (rrn5) appear to be absent from most mitochondrial and several chloroplast genomes. Here, we re-examine the distribution of organelle rrn5 by building mitochondrion- and plastid-specific covariance models (CMs) with which we screened organelle genome sequences. We not only recover all organelle rrn5 genes annotated in GenBank records, but also identify more than 50 previously unrecognized homologs in mitochondrial genomes of various stramenopiles, red algae, cryptomonads, malawimonads and apusozoans, and surprisingly, in the apicoplast (highly derived plastid) genomes of the coccidian pathogens Toxoplasma gondii and Eimeria tenella. Comparative modeling of RNA secondary structure reveals that mitochondrial 5S rRNAs from brown algae adopt a permuted triskelion shape that has not been seen elsewhere. Expression of the newly predicted rrn5 genes is confirmed experimentally in 10 instances, based on our own and published RNA-Seq data. This study establishes that particularly mitochondrial 5S rRNA has a much broader taxonomic distribution and a much larger structural variability than previously thought. The newly developed CMs will be made available via the Rfam database and the MFannot organelle genome annotator.     

The GC-rich mitochondrial and plastid genomes of the green alga Coccomyxa give insight into the evolution of organelle DNA nucleotide landscape

Most of the available mitochondrial and plastid genome sequences are biased towards adenine and thymine (AT) over guanine and cytosine (GC). Examples of GC-rich organelle DNAs are limited to a small but eclectic list of species, including certain green algae. Here, to gain insight in the evolution of organelle nucleotide landscape, we present the GC-rich mitochondrial and plastid DNAs from the trebouxiophyte green alga Coccomyxa sp. C-169. We compare these sequences with other GC-rich organelle DNAs and argue that the forces biasing them towards G and C are nonadaptive and linked to the metabolic and/or life history features of this species. The Coccomyxa organelle genomes are also used for phylogenetic analyses, which highlight the complexities in trying to resolve the interrelationships among the core chlorophyte green algae, but ultimately favour a sister relationship between the Ulvophyceae and Chlorophyceae, with the Trebouxiophyceae branching at the base of the chlorophyte crown.     

Mitochondrial and Plastid Genomes of the Colonial Green Alga Gonium pectorale Give Insights into the Origins of Organelle DNA Architecture within the Volvocales

Volvocalean green algae have among the most diverse mitochondrial and plastid DNAs (mtDNAs and ptDNAs) from the eukaryotic domain. However, nearly all of the organelle genome data from this group are restricted to unicellular species, like Chlamydomonas reinhardtii, and presently only one multicellular species, the ∼4,000-celled Volvox carteri, has had its organelle DNAs sequenced. The V. carteri organelle genomes are repeat rich, and the ptDNA is the largest plastome ever sequenced. Here, we present the complete mtDNA and ptDNA of the colonial volvocalean Gonium pectorale, which is comprised of ∼16 cells and occupies a phylogenetic position closer to that of V. carteri than C. reinhardtii within the volvocine line. The mtDNA and ptDNA of G. pectorale are circular-mapping AT-rich molecules with respective lengths and coding densities of 16 and 222.6 kilobases and 73 and 44%. They share some features with the organelle DNAs of V. carteri, including palindromic repeats within the plastid compartment, but show more similarities with those of C. reinhardtii, such as a compact mtDNA architecture and relatively low organelle DNA intron contents. Overall, the G. pectorale organelle genomes raise several interesting questions about the origin of linear mitochondrial chromosomes within the Volvocales and the relationship between multicellularity and organelle genome expansion.     

Analysis of mitochondrial and chloroplast genomes in two volvocine algae: Eudorina elegans and Eudorina cylindrica (Volvocaceae,Chlorophyta)

The colonial volvocine algae span the full range of organizational complexity, from four-celled species to multicellular species, and this group of algae is often used for the study of evolution. In recent years, many organelle genomes have been sequenced using the application of next generation sequencing technology; however, only a few organelle genomes have been reported in colonial volvocine algae. In this study, we determined the organelle genomes of Eudorina elegans and Eudorina cylindrica and analysed the organelle genome size, structure and gene content between these volvocine species. This provided useful information to help us understand the composition of colonial volvocine organelle genomes. Based on the chloroplast genome protein-coding genes, we conducted a phylogenomics analysis of the volvocine algae. The result revealed an unexpected phylogenetic relationship, namely, E. elegans is more closely related to Pleodorina starrii than to E. cylindrica. The substitution rate of volvocine algae was then calculated based on organelle genome protein-coding genes; our analysis suggested the possibility that the two Eudorina species may be under similar evolutionary pressure. Lastly, the synteny analysis of the mitochondrial genome showed that gene arrangements and contents are highly conserved in the family Volvocaceae, and the synteny analysis of the chloroplast genome indicated that the genus Eudorina may have experienced genomic changes.     

Microsatellite analysis in organelle genomes of Chlorophyta

Simple Sequence Repeats (SSRs) or microsatellites constitute a significant portion of genomes however; their significance in organellar genomes has not been completely understood. The availability of organelle genome sequences allows us to understand the organization of SSRs in their genic and intergenic regions. In the present work, SSRs were identified and categorized in 14 mitochondrial and 22 chloroplast genomes of algal species belonging to Chlorophyta. Based on the study, it was observed that number of SSRs in non-coding region were more as compared to coding region and frequency of mononucleotides repeats were highest followed by dinucleotides in both mitochondrial and chloroplast genomes. It was also observed that maximum number of SSRs was found in genes encoding for beta subunit of RNA polymerase in chloroplast genomes and NADH dehydrogenase in mitochondrial genomes. This is the first and original report on whole genomes sequence analysis of organellar genomes of green algae.     

The Complete Chloroplast and Mitochondrial Genomes of the Green Macroalga Ulva sp. UNA00071828 (Ulvophyceae,Chlorophyta)

Sequencing mitochondrial and chloroplast genomes has become an integral part in understanding the genomic machinery and the phylogenetic histories of green algae. Previously, only three chloroplast genomes (Oltmannsiellopsis viridis, Pseudendoclonium akinetum, and Bryopsis hypnoides) and two mitochondrial genomes (O. viridis and P. akinetum) from the class Ulvophyceae have been published. Here, we present the first chloroplast and mitochondrial genomes from the ecologically and economically important marine, green algal genus Ulva. The chloroplast genome of Ulva sp. was 99,983 bp in a circular-mapping molecule that lacked inverted repeats, and thus far, was the smallest ulvophycean plastid genome. This cpDNA was a highly compact, AT-rich genome that contained a total of 102 identified genes (71 protein-coding genes, 28 tRNA genes, and three ribosomal RNA genes). Additionally, five introns were annotated in four genes: atpA (1), petB (1), psbB (2), and rrl (1). The circular-mapping mitochondrial genome of Ulva sp. was 73,493 bp and follows the expanded pattern also seen in other ulvophyceans and trebouxiophyceans. The Ulva sp. mtDNA contained 29 protein-coding genes, 25 tRNA genes, and two rRNA genes for a total of 56 identifiable genes. Ten introns were annotated in this mtDNA: cox1 (4), atp1 (1), nad3 (1), nad5 (1), and rrs (3). Double-cut-and-join (DCJ) values showed that organellar genomes across Chlorophyta are highly rearranged, in contrast to the highly conserved organellar genomes of the red algae (Rhodophyta). A phylogenomic investigation of 51 plastid protein-coding genes showed that Ulvophyceae is not monophyletic, and also placed Oltmannsiellopsis (Oltmannsiellopsidales) and Tetraselmis (Chlorodendrophyceae) closely to Ulva (Ulvales) and Pseudendoclonium (Ulothrichales).     

Cytoplasmic inheritance in green algae: patterns,mechanisms and relation to sex type

Cytological and genetic investigations of two major groups of green algae, chlorophyte and streptophyte green algae, show a predominance of uniparental inheritance of the plastid and mitochondrial genomes in most species. However, in some crosses of isogamous species of Ulva compressa, these genomes are transmitted from mt⁺, mt⁻, and both parents. In species with uniparental organelle inheritance, various mechanisms can eliminate organelles and their DNA during male gametogenesis or after fertilization. Concerning plastid inheritance, two major mechanisms are widespread in green algae: (1) digestion of plastid DNA during male gametogenesis, during fertilization, or after fertilization; and (2) disintegration or fusion of the plastid in the zygote. The first mechanism also eliminates the mitochondrial DNA in anisogamous and oogamous species. These mechanisms would ensure the predominantly uniparental inheritance of organelle genomes in green algae. To trace the evolutionary history of cytoplasmic inheritance in green algae, the relations between uniparental inheritance and sex type were considered in isogamous, anisogamous, and oogamous species using sex-specific features that might be nearly universal among Chlorophyta.     

Decoding and analysis of organelle genomes of Indian tea (Camellia assamica) for phylogenetic confirmation

The NCBI database has >15 chloroplast (cp) genome sequences available for different Camellia species but none for C. assamica. There is no report of any mitochondrial (mt) genome in the Camellia genus or Theaceae family. With the strong believes that these organelle genomes can play a great tool for taxonomic and phylogenetic analysis, we successfully assembled and analyzed cp and mt genome of C. assamica. We assembled the complete mt genome of C. assamica in a single circular contig of 707,441 bp length comprising of a total of 66 annotated genes, including 35 protein-coding genes, 29 tRNAs and two rRNAs. The first ever cp genome of C. assamica resulted in a circular contig of 157,353 bp length with a typical quadripartite structure. Phylogenetic analysis based on these organelle genomes showed that C. assamica was closely related to C. sinensis and C. leptophylla. It also supports Caryophyllales as Superasterids.     

The Chlorella variabilis NC64A Genome Reveals Adaptation to Photosymbiosis,Coevolution with Viruses,and Cryptic Sex

Chlorella variabilis NC64A, a unicellular photosynthetic green alga (Trebouxiophyceae), is an intracellular photobiont of Paramecium bursaria and a model system for studying virus/algal interactions. We sequenced its 46-Mb nuclear genome, revealing an expansion of protein families that could have participated in adaptation to symbiosis. NC64A exhibits variations in GC content across its genome that correlate with global expression level, average intron size, and codon usage bias. Although Chlorella species have been assumed to be asexual and nonmotile, the NC64A genome encodes all the known meiosis-specific proteins and a subset of proteins found in flagella. We hypothesize that Chlorella might have retained a flagella-derived structure that could be involved in sexual reproduction. Furthermore, a survey of phytohormone pathways in chlorophyte algae identified algal orthologs of Arabidopsis thaliana genes involved in hormone biosynthesis and signaling, suggesting that these functions were established prior to the evolution of land plants. We show that the ability of Chlorella to produce chitinous cell walls likely resulted from the capture of metabolic genes by horizontal gene transfer from algal viruses, prokaryotes, or fungi. Analysis of the NC64A genome substantially advances our understanding of the green lineage evolution, including the genomic interplay with viruses and symbiosis between eukaryotes.     

Origin and phylogeny of chloroplasts revealed by a simple correlation analysis of complete genomes

The complete sequenced genomes of chloroplast have provided much information on the origin and evolution of this organelle. In this paper we attempt to use these sequences to test a novel approach for phylogenetic analysis of complete genomes based on correlation analysis of compositional vectors. All protein sequences from 21 complete chloroplast genomes are analyzed in comparison with selected archaea, eubacteria, and eukaryotes. The distance-based analysis shows that the chloroplast genomes are most closely related to cyanobacteria, consistent with the endosymbiotic origin of chloroplasts. The chloroplast genomes are separated to two major clades corresponding to chlorophytes (green plants) s.l. and rhodophytes (red algae) s.l. The interrelationships among the chloroplasts are largely in agreement with the current understanding on chloroplast evolution. For instance, the analysis places the chloroplasts of two chromophytes (Guillardia and Odontella) within the rhodophyte lineage, supporting secondary endosymbiosis as the source of these chloroplasts. The relationships among the green algae and land plants in our tree also agree with results from traditional phylogenetic analyses. Thus, this study establishes the value of our simple correlation analysis in elucidating the evolutionary relationships among genomes. It is hoped that this approach will provide insights on comparative genome analysis.     

In silico analysis of SSRs in mitochondrial genomes of plants

Simple sequence repeats (SSRs) or microsatellites constitute a countable portion of genomes. However, the significance of SSRs in organelle genomes has not been completely understood. The availability of organelle genome sequences allows us to understand the organization of SSRs in their genic and intergenic regions. In the current study we surveyed the patterns of SSRs in mitochondrial genomes of different taxa of plants. A total of 16 mitochondrial genomes, from algae to angiosperms, have been considered to analyze the pattern of simple sequence repeats present in them. Based on study, the mononucleotide repeats of A/T were found to be more prevalent in mitochondrial genomes over other repeat types. The dinucleotides repeats, TA/AT, were the second most numerous, whereas tri-, tetra-, and pentanucleotide repeats were in less number and present in intronic or intergenic portions only. Mononucleotide repeats prevailed in protein-coding exonic portions of all organisms. These results indicates that microsatellite pattern in mitochondrial genomes is different from nuclear genomes and also focuses on organization and diversity at SSR locuses in mitochondrial genomes. This is the novel report of microsatellite polymorphism in plant mitochondrion on whole genome level.     

Babesia bovis: A comprehensive phylogenetic analysis of plastid-encoded genes supports green algal origin of apicoplasts

Apicomplexan parasites commonly contain a unique, non-photosynthetic plastid-like organelle termed the apicoplast. Previous analyses of other plastid-containing organisms suggest that apicoplasts were derived from a red algal ancestor. In this report, we present an extensive phylogenetic study of apicoplast origins using multiple previously reported apicoplast sequences as well as several sequences recently reported. Phylogenetic analysis of amino acid sequences was used to determine the evolutionary origin of the organelle. A total of nine plastid genes from 37 species were incorporated in our study. The data strongly support a green algal origin for apicoplasts and Euglenozoan plastids. Further, the nearest green algae lineage to the Apicomplexans is the parasite Helicosporidium, suggesting that apicoplasts may have originated by lateral transfer from green algal parasite lineages. The results also substantiate earlier findings that plastids found in Heterokonts such as Odontella and Thalassiosira were derived from a separate secondary endosymbiotic event likely originating from a red algal lineage.     

Genome Analysis of the Janthinobacterium sp. Strain SLB01 from the Diseased Sponge of the Lubomirskia baicalensis

The strain Janthinobacterium sp. SLB01 was isolated from the diseased freshwater sponge Lubomirskia baicalensis (Pallas, 1776) and the draft genome was published previously. The aim of this work is to analyze the genome of the Janthinobacterium sp. SLB01 to search for pathogenicity factors for Baikal sponges. We performed genomic analysis to determine virulence factors, comparing the genome of the strain SLB01 with genomes of other related J. lividum strains from the environment. The strain Janthinobacterium sp. SLB01 contained genes encoding violacein, alpha-amylases, phospholipases, chitinases, collagenases, hemolysin, and a type VI secretion system. In addition, the presence of conservative clusters of genes for the biosynthesis of secondary metabolites of tropodithietic acid and marinocine was found. We present genes for antibiotic resistance, including five genes encoding various lactamases and eight genes for penicillin-binding proteins, which are conserved in all analyzed strains. Major differences were found between the Janthinobacterium sp. SLB01 and J. lividum strains in the spectra of genes for glycosyltransferases and glycoside hydrolases, serine hydrolases, and trypsin-like peptidase, as well as some TonB-dependent siderophore receptors. Thus, the study of the analysis of the genome of the strain SLB01 allows us to conclude that the strain may be one of the pathogens of freshwater sponges.     

Odintifier - A computational method for identifying insertions of organellar origin from modern and ancient high-throughput sequencing data based on haplotype phasing

Background

Cellular organelles with genomes of their own (e.g. plastids and mitochondria) can pass genetic sequences to other organellar genomes within the cell in many species across the eukaryote phylogeny. The extent of the occurrence of these organellar-derived inserted sequences (odins) is still unknown, but if not accounted for in genomic and phylogenetic studies, they can be a source of error. However, if correctly identified, these inserted sequences can be used for evolutionary and comparative genomic studies. Although such insertions can be detected using various laboratory and bioinformatic strategies, there is currently no straightforward way to apply them as a standard organellar genome assembly on next-generation sequencing data. Furthermore, most current methods for identification of such insertions are unsuitable for use on non-model organisms or ancient DNA datasets.

Results

We present a bioinformatic method that uses phasing algorithms to reconstruct both source and inserted organelle sequences. The method was tested in different shotgun and organellar-enriched DNA high-throughput sequencing (HTS) datasets from ancient and modern samples. Specifically, we used datasets from lions (Panthera leo ssp. and Panthera leo leo) to characterize insertions from mitochondrial origin, and from common grapevine (Vitis vinifera) and bugle (Ajuga reptans) to characterize insertions derived from plastid genomes. Comparison of the results against other available organelle genome assembly methods demonstrated that our new method provides an improvement in the sequence assembly.

Conclusion

Using datasets from a wide range of species and different levels of complexity we showed that our novel bioinformatic method based on phasing algorithms can be used to achieve the next two goals: i) reference-guided assembly of chloroplast/mitochondrial genomes from HTS data and ii) identification and simultaneous assembly of odins. This method represents the first application of haplotype phasing for automatic detection of odins and reference-based organellar genome assembly.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0682-1) contains supplementary material, which is available to authorized users.     

Chromosome-scale haplotype-phased genome assemblies of the male and female lines of wild asparagus (Asparagus kiusianus), a dioecious plant species

Asparagus kiusianus is a disease-resistant dioecious plant species and a wild relative of garden asparagus (Asparagus officinalis). To enhance A. kiusianus genomic resources, advance plant science, and facilitate asparagus breeding, we determined the genome sequences of the male and female lines of A. kiusianus. Genome sequence reads obtained with a linked-read technology were assembled into four haplotype-phased contig sequences (∼1.6 Gb each) for the male and female lines. The contig sequences were aligned onto the chromosome sequences of garden asparagus to construct pseudomolecule sequences. Approximately 55,000 potential protein-encoding genes were predicted in each genome assembly, and ∼70% of the genome sequence was annotated as repetitive. Comparative analysis of the genomes of the two species revealed structural and sequence variants between the two species as well as between the male and female lines of each species. Genes with high sequence similarity with the male-specific sex determinant gene in A. officinalis, MSE1/AoMYB35/AspTDF1, were presented in the genomes of the male line but absent from the female genome assemblies. Overall, the genome sequence assemblies, gene sequences, and structural and sequence variants determined in this study will reveal the genetic mechanisms underlying sexual differentiation in plants, and will accelerate disease-resistance breeding in garden asparagus.     

Comparative study on mitogenomes of green tide algae

Since 2007, the annual green tide disaster in the Yellow Sea has brought serious economic losses to China. There is no research on the genetic similarities of four constituent species of green tide algae at the genomic level. We previously determined the mitochondrial genomes of Ulva prolifera, Ulva linza and Ulva flexuosa. In the present work, the mitochondrial genome of another green tide (Ulva compressa) was sequenced and analyzed. With the length of 62,311 bp, it contained 29 encoding genes, 26 tRNAs and 10 open reading frames. By comparing these four mitochondrial genomes, we found that U. compressa was quite different from the other three types of Ulva species. However, there were similarities between U. prolifera and U. linza in the number, distribution and homology of open reading frames, evolutionary and codon variation of tRNA, evolutionary relationship and selection pressure of coding genes. Repetitive sequence analysis of simple sequence repeats, tandem repeat and forward repeats further supposed that they have evolved from the same origin. In addition, we directly analyzed gene homologies and translocation of four green tide algae by Mauve alignment. There were gene order rearrangements among them. With fast-evolving genomes, these four green algal mitochondria have both conservatism and variation, thus opening another window for the understanding of origin and evolution of Ulva.     

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.     

Peroxisomal targeting signals in green algae

Peroxisomal enzymatic proteins contain targeting signals (PTS) to enable their import into peroxisomes. These targeting signals have been identified as PTS1 and PTS2 in mammalian, yeast, and higher plant cells; however, no PTS2-like amino acid sequences have been observed in enzymes from the genome database of Cyanidiochyzon merolae (Bangiophyceae), a primitive red algae. In studies on the evolution of PTS, it is important to know when their sequences came to be the peroxisomal targeting signals for all living organisms. To this end, we identified a number of genes in the genome database of the green algae Chlamydomonas reinhardtii, which contains amino acid sequences similar to those found in plant PTS. In order to determine whether these sequences function as PTS in green algae, we expressed modified green fluorescent proteins (GFP) fused to these putative PTS peptides under the cauliflower mosaic virus 35S promoter. To confirm whether granular structures containing GFP–PTS fusion proteins accumulated in the peroxisomes of Closterium ehrenbergii, we observed these cells after the peroxisomes were stained with 3, 3′-diaminobenzidine. Our results confirm that the GFP–PTS fusion proteins indeed accumulated in the peroxisomes of these green algae. These findings suggest that the peroxisomal transport system for PTS1 and PTS2 is conserved in green algal cells and that our fusion proteins can be used to visualize peroxisomes in live cells.