Two RhoGEF isoforms with distinct localisation control furrow position during asymmetric cell division.

Cytokinesis partitions cellular content between daughter cells. It relies on the formation of an acto-myosin contractile ring, whose constriction induces the ingression of the cleavage furrow between the segregated chromatids. Rho1 GTPase and its RhoGEF (Pbl) are essential for this process. However, how Rho1 is regulated to sustain furrow ingression while maintaining correct furrow position remains poorly defined. Here, we show that during asymmetric division of Drosophila neuroblasts, Rho1 is controlled by two Pbl isoforms with distinct localisation. Spindle midzone- and furrow-enriched Pbl-A focuses Rho1 at the furrow to sustain efficient ingression, while Pbl-B pan-plasma membrane localization promotes the broadening of Rho1 activity and the subsequent enrichment of myosin on the entire cortex. This enlarged zone of Rho1 activity is critical to adjust furrow position, thereby preserving correct daughter cell size asymmetry. Our work highlights how the use of isoforms with distinct localisation makes an essential process more robust.

Core genome multilocus sequence typing scheme for Bacillus cereus group bacteria.

Core genome multilocus sequence typing (cgMLST) employs a strategy where the set of orthologous genes common to all members of a group of organisms are used for phylogenetic analysis of the group members. The Bacillus cereus group consists of species with pathogenicity towards insect species as well as warm-blooded animals including humans. While B. cereus is an opportunistic pathogen linked to a range of human disease conditions, including emesis and diarrhoea, Bacillus thuringiensis is an entomopathogenic species with toxicity toward insect larvae, and therefore used as a biological pesticide worldwide. Bacillus anthracis is a classical obligate pathogen causing anthrax, an acute lethal condition in herbivores as well as humans, and which is endemic in many parts of the world. The group also includes a range of additional species, and B. cereus group bacteria have been subject to analysis with a wide variety of phylogenetic typing systems. Here we present, based on analyses of 173 complete genomes from B. cereus group species available in public databases, the identification of a set of 1568 core genes which were used to create a core genome multilocus typing scheme for the group which is implemented in the PubMLST system as an open online database freely available to the community. The new cgMLST system provides unprecedented resolution over existing phylogenetic analysis schemes covering the B. cereus group.

Profiling the intragenic toxicity determinants of toxin-antitoxin systems: revisiting hok/Sok regulation.

Type I toxin-antitoxin systems (T1TAs) are extremely potent bacterial killing systems difficult to characterize using classical approaches. To assess the killing capability of type I toxins and to identify mutations suppressing the toxin expression or activity, we previously developed the FASTBAC-Seq (Functional AnalysiS of Toxin-Antitoxin Systems in BACteria by Deep Sequencing) method in Helicobacter pylori. This method combines a life and death selection with deep sequencing. Here, we adapted and improved our method to investigate T1TAs in the model organism Escherichia coli. As a proof of concept, we revisited the regulation of the plasmidic hok/Sok T1TA system. We revealed the death-inducing phenotype of the Hok toxin when it is expressed from the chromosome in the absence of the antitoxin and recovered previously described intragenic toxicity determinants of this system. We identified nucleotides that are essential for the transcription, translation or activity of Hok. We also discovered single-nucleotide substitutions leading to structural changes affecting either the translation or the stability of the hok mRNA. Overall, we provide the community with an easy-to-use approach to widely characterize TA systems from diverse types and bacteria.

Adaptation of Bacillus thuringiensis to Plant Colonization Affects Differentiation and Toxicity.

The Bacillus cereus group (Bacillus cereus sensu lato) has a diverse ecology, including various species that are vertebrate or invertebrate pathogens. Few isolates from the B. cereus group have however been demonstrated to benefit plant growth. Therefore, it is crucial to explore how bacterial development and pathogenesis evolve during plant colonization. Herein, we investigated Bacillus thuringiensis (Cry-) adaptation to the colonization of Arabidopsis thaliana roots and monitored changes in cellular differentiation in experimentally evolved isolates. Isolates from two populations displayed improved iterative ecesis on roots and increased virulence against insect larvae. Molecular dissection and recreation of a causative mutation revealed the importance of a nonsense mutation in the rho transcription terminator gene. Transcriptome analysis revealed how Rho impacts various B. thuringiensis genes involved in carbohydrate metabolism and virulence. Our work suggests that evolved multicellular aggregates have a fitness advantage over single cells when colonizing plants, creating a trade-off between swimming and multicellularity in evolved lineages, in addition to unrelated alterations in pathogenicity. IMPORTANCE Biologicals-based plant protection relies on the use of safe microbial strains. During application of biologicals to the rhizosphere, microbes adapt to the niche, including genetic mutations shaping the physiology of the cells. Here, the experimental evolution of Bacillus thuringiensis lacking the insecticide crystal toxins was examined on the plant root to reveal how adaptation shapes the differentiation of this bacterium. Interestingly, evolution of certain lineages led to increased hemolysis and insect larva pathogenesis in B. thuringiensis driven by transcriptional rewiring. Further, our detailed study reveals how inactivation of the transcription termination protein Rho promotes aggregation on the plant root in addition to altered differentiation and pathogenesis in B. thuringiensis.

T1TAdb: the database of type I toxin-antitoxin systems.

Type I toxin-antitoxin (T1TA) systems constitute a large class of genetic modules with antisense RNA (asRNA)-mediated regulation of gene expression. They are widespread in bacteria and consist of an mRNA coding for a toxic protein and a noncoding asRNA that acts as an antitoxin preventing the synthesis of the toxin by directly base-pairing to its cognate mRNA. The co- and post-transcriptional regulation of T1TA systems is intimately linked to RNA sequence and structure, therefore it is essential to have an accurate annotation of the mRNA and asRNA molecules to understand this regulation. However, most T1TA systems have been identified by means of bioinformatic analyses solely based on the toxin protein sequences, and there is no central repository of information on their specific RNA features. Here we present the first database dedicated to type I TA systems, named T1TAdb. It is an open-access web database (https://d-lab.arna.cnrs.fr/t1tadb) with a collection of ∼1900 loci in ∼500 bacterial strains in which a toxin-coding sequence has been previously identified. RNA molecules were annotated with a bioinformatic procedure based on key determinants of the mRNA structure and the genetic organization of the T1TA loci. Besides RNA and protein secondary structure predictions, T1TAdb also identifies promoter, ribosome-binding, and mRNA-asRNA interaction sites. It also includes tools for comparative analysis, such as sequence similarity search and computation of structural multiple alignments, which are annotated with covariation information. To our knowledge, T1TAdb represents the largest collection of features, sequences, and structural annotations on this class of genetic modules.

Phylogenetic and functional diversity of aldehyde-alcohol dehydrogenases in microalgae.

KEY MESSAGE: The study shows the biochemical and enzymatic divergence between the two aldehyde-alcohol dehydrogenases of the alga Polytomella sp., shedding light on novel aspects of the enzyme evolution amid unicellular eukaryotes. Aldehyde-alcohol dehydrogenases (ADHEs) are large metalloenzymes that typically perform the two-step reduction of acetyl-CoA into ethanol. These enzymes consist of an N-terminal acetylating aldehyde dehydrogenase domain (ALDH) and a C-terminal alcohol dehydrogenase (ADH) domain. ADHEs are present in various bacterial phyla as well as in some unicellular eukaryotes. Here we focus on ADHEs in microalgae, a diverse and polyphyletic group of plastid-bearing unicellular eukaryotes. Genome survey shows the uneven distribution of the ADHE gene among free-living algae, and the presence of two distinct genes in various species. We show that the non-photosynthetic Chlorophyte alga Polytomella sp. SAG 198.80 harbors two genes for ADHE-like enzymes with divergent C-terminal ADH domains. Immunoblots indicate that both ADHEs accumulate in Polytomella cells growing aerobically on acetate or ethanol. ADHE1 of ~ 105-kDa is found in particulate fractions, whereas ADHE2 of ~ 95-kDa is mostly soluble. The study of the recombinant enzymes revealed that ADHE1 has both the ALDH and ADH activities, while ADHE2 has only the ALDH activity. Phylogeny shows that the divergence occurred close to the root of the Polytomella genus within a clade formed by the majority of the Chlorophyte ADHE sequences, next to the cyanobacterial clade. The potential diversification of function in Polytomella spp. unveiled here likely took place after the loss of photosynthesis. Overall, our study provides a glimpse at the complex evolutionary history of the ADHE in microalgae which includes (i) acquisition via different gene donors, (ii) gene duplication and (iii) independent evolution of one of the two enzymatic domains.

MogR Is a Ubiquitous Transcriptional Repressor Affecting Motility, Biofilm Formation and Virulence in Bacillus thuringiensis.

Flagellar motility is considered an important virulence factor in different pathogenic bacteria. In Listeria monocytogenes the transcriptional repressor MogR regulates motility in a temperature-dependent manner, directly repressing flagellar- and chemotaxis genes. The only other bacteria known to carry a mogR homolog are members of the Bacillus cereus group, which includes motile species such as B. cereus and Bacillus thuringiensis as well as the non-motile species Bacillus anthracis, Bacillus mycoides and Bacillus pseudomycoides. Furthermore, the main motility locus in B. cereus group bacteria, carrying the genes for flagellar synthesis, appears to be more closely related to L. monocytogenes than to Bacillus subtilis, which belongs to a separate phylogenetic group of Bacilli and does not carry a mogR ortholog. Here, we show that in B. thuringiensis, MogR overexpression results in non-motile cells devoid of flagella. Global gene expression profiling showed that 110 genes were differentially regulated by MogR overexpression, including flagellar motility genes, but also genes associated with virulence, stress response and biofilm lifestyle. Accordingly, phenotypic assays showed that MogR also affects cytotoxicity and biofilm formation in B. thuringiensis. Overexpression of a MogR variant mutated in two amino acids within the putative DNA binding domain restored phenotypes to those of an empty vector control. In accordance, introduction of these mutations resulted in complete loss in MogR binding to its candidate flagellar locus target site in vitro. In contrast to L. monocytogenes, MogR appears to be regulated in a growth-phase dependent and temperature-independent manner in B. thuringiensis 407. Interestingly, mogR was found to be conserved also in non-motile B. cereus group species such as B. mycoides and B. pseudomycoides, which both carry major gene deletions in the flagellar motility locus and where in B. pseudomycoides mogR is the only gene retained. Furthermore, mogR is expressed in non-motile B. anthracis. Altogether this provides indications of an expanded set of functions for MogR in B. cereus group species, beyond motility regulation. In conclusion, MogR constitutes a novel B. thuringiensis pleiotropic transcriptional regulator, acting as a repressor of motility genes, and affecting the expression of a variety of additional genes involved in biofilm formation and virulence.

Structural Alignment and Covariation Analysis of RNA Sequences.

RNA molecules adopt defined structural conformations that are essential to exert their function. During the course of evolution, the structure of a given RNA can be maintained via compensatory base-pair changes that occur among covarying nucleotides in paired regions. Therefore, for comparative, structural, and evolutionary studies of RNA molecules, numerous computational tools have been developed to incorporate structural information into sequence alignments and a number of tools have been developed to study covariation. The bioinformatic protocol presented here explains how to use some of these tools to generate a secondary-structure-aware multiple alignment of RNA sequences and to annotate the alignment to examine the conservation and covariation of structural elements among the sequences.

The BF4 and p71 antenna mutants from Chlamydomonas reinhardtii.

Two pale green mutants of the green alga Chlamydomonas reinhardtii, which have been used over the years in many photosynthesis studies, the BF4 and p71 mutants, were characterized and their mutated gene identified in the nuclear genome. The BF4 mutant is defective in the insertase Alb3.1 whereas p71 is defective in cpSRP43. The two mutants showed strikingly similar deficiencies in most of the peripheral antenna proteins associated with either photosystem I or photosystem 2. As a result the two photosystems have a reduced antenna size with photosystem 2 being the most affected. Still up to 20% of the antenna proteins remain in these strains, with the heterodimer Lhca5/Lhca6 showing a lower sensitivity to these mutations. We discuss these phenotypes in light of those of other allelic mutants that have been described in the literature and suggest that eventhough the cpSRP route serves as the main biogenesis pathway for antenna proteins, there should be an escape pathway which remains to be genetically identified.

A genetic selection reveals functional metastable structures embedded in a toxin-encoding mRNA.

Post-transcriptional regulation plays important roles to fine-tune gene expression in bacteria. In particular, regulation of type I toxin-antitoxin (TA) systems is achieved through sophisticated mechanisms involving toxin mRNA folding. Here, we set up a genetic approach to decipher the molecular underpinnings behind the regulation of a type I TA in Helicobacter pylori. We used the lethality induced by chromosomal inactivation of the antitoxin to select mutations that suppress toxicity. We found that single point mutations are sufficient to allow cell survival. Mutations located either in the 5' untranslated region or within the open reading frame of the toxin hamper its translation by stabilizing stem-loop structures that sequester the Shine-Dalgarno sequence. We propose that these short hairpins correspond to metastable structures that are transiently formed during transcription to avoid premature toxin expression. This work uncovers the co-transcriptional inhibition of translation as an additional layer of TA regulation in bacteria.

Mapping and characterization of G-quadruplexes in the genome of the social amoeba Dictyostelium discoideum.

G-quadruplexes (G4) are non-canonical DNA and/or RNA secondary structures formed in guanine-rich regions. Given their over-representation in specific regions in the genome such as promoters and telomeres, they are likely to play important roles in key processes such as transcription, replication or RNA maturation. Putative G4-forming sequences (G4FS) have been reported in humans, yeast, bacteria, viruses and many organisms. Here we present the first mapping of G-quadruplex sequences in Dictyostelium discoideum, the social amoeba. 'Dicty' is an ameboid protozoan with a small (34 Mb) and extremely AT rich genome (78%). As a consequence, very few G4-prone motifs are expected. An in silico analysis of the Dictyostelium genome with the G4Hunter software detected 249-1055 G4-prone motifs, depending on G4Hunter chosen threshold. Interestingly, despite an even lower GC content (as compared to the whole Dicty genome), the density of G4 motifs in Dictyostelium promoters and introns is significantly higher than in the rest of the genome. Fourteen selected sequences located in important genes were characterized by a combination of biophysical and biochemical techniques. Our data show that these sequences form highly stable G4 structures under physiological conditions. Five Dictyostelium genes containing G4-prone motifs in their promoters were studied for the effect of a new G4-binding porphyrin derivative on their expression. Our results demonstrated that the new ligand significantly decreased their expression. Overall, our results constitute the first step to adopt Dictyostelium discoideum as a 'G4-poor' model for studies on G-quadruplexes.

FASTBAC-Seq: Functional Analysis of Toxin-Antitoxin Systems in Bacteria by Deep Sequencing.

As the number of bacterial genomes and transcriptomes increases, so does the number of newly identified toxin-antitoxin (TA) systems. However, their functional characterization remains challenging, often requiring the use of overexpression vectors that can lead to misinterpretations of in vivo results. To fill this gap, we developed a systematic approach called FASTBAC-Seq (Functional AnalysiS of Toxin-Antitoxin Systems in BACteria by Deep Sequencing). Combining life/death phenotypic selection with next-generation sequencing, FASTBAC-Seq allows the rapid identification of loss-of-function (toxicity) mutations in toxin-encoding genes belonging to TA loci with nucleotide resolution. Here, we present the setup used on the first-time application of FASBACT-Seq to characterize a member of the aapA/IsoA family of type I TA systems hosted on the chromosome of the major human gastric pathogen Helicobacter pylori. We propose FASBACT-Seq as a powerful tool for the functional characterization of TA systems that can in addition uncover key elements for the understanding of gene expression regulation in bacteria.

Modulation of oncogenic miRNA biogenesis using functionalized polyamines.

MicroRNAs are key factors in the regulation of gene expression and their deregulation has been directly linked to various pathologies such as cancer. The use of small molecules to tackle the overexpression of oncogenic miRNAs has proved its efficacy and holds the promise for therapeutic applications. Here we describe the screening of a 640-compound library and the identification of polyamine derivatives interfering with in vitro Dicer-mediated processing of the oncogenic miR-372 precursor (pre-miR-372). The most active inhibitor is a spermine-amidine conjugate that binds to the pre-miR-372 with a KD of 0.15 µM, and inhibits its in vitro processing with a IC50 of 1.06 µM. The inhibition of miR-372 biogenesis was confirmed in gastric cancer cells overexpressing miR-372 and a specific inhibition of proliferation through de-repression of the tumor suppressor LATS2 protein, a miR-372 target, was observed. This compound modifies the expression of a small set of miRNAs and its selective biological activity has been confirmed in patient-derived ex vivo cultures of gastric carcinoma. Polyamine derivatives are promising starting materials for future studies about the inhibition of oncogenic miRNAs and, to the best of our knowledge, this is the first report about the application of functionalized polyamines as miRNAs interfering agents.

Regulation of RNA polymerase III transcription during transformation of human IMR90 fibroblasts with defined genetic elements.

RNA polymerase (Pol) III transcribes small untranslated RNAs that are essential for cellular homeostasis and growth. Its activity is regulated by inactivation of tumor suppressor proteins and overexpression of the oncogene c-MYC, but the concerted action of these tumor-promoting factors on Pol III transcription has not yet been assessed. In order to comprehensively analyse the regulation of Pol III transcription during tumorigenesis we employ a model system that relies on the expression of five genetic elements to achieve cellular transformation. Expression of these elements in six distinct transformation intermediate cell lines leads to the inactivation of TP53, RB1, and protein phosphatase 2A, as well as the activation of RAS and the protection of telomeres by TERT, thereby conducting to full tumoral transformation of IMR90 fibroblasts. Transformation is accompanied by moderately enhanced levels of a subset of Pol III-transcribed RNAs (7SK; MRP; H1). In addition, mRNA and/or protein levels of several Pol III subunits and transcription factors are upregulated, including increased protein levels of TFIIIB and TFIIIC subunits, of SNAPC1 and of Pol III subunits. Strikingly, the expression of POLR3G and of SNAPC1 is strongly enhanced during transformation in this cellular transformation model. Collectively, our data indicate that increased expression of several components of the Pol III transcription system accompanied by a 2-fold increase in steady state levels of a subset of Pol III RNAs is sufficient for sustaining tumor formation.

Quantitative RNA-seq meta-analysis of alternative exon usage in C. elegans.

Almost 20 years after the completion of the C. elegans genome sequence, gene structure annotation is still an ongoing process with new evidence for gene variants still being regularly uncovered by additional in-depth transcriptome studies. While alternative splice forms can allow a single gene to encode several functional isoforms, the question of how much spurious splicing is tolerated is still heavily debated. Here we gathered a compendium of 1682 publicly available C. elegans RNA-seq data sets to increase the dynamic range of detection of RNA isoforms, and obtained robust measurements of the relative abundance of each splicing event. While most of the splicing reads come from reproducibly detected splicing events, a large fraction of purported junctions is only supported by a very low number of reads. We devised an automated curation method that takes into account the expression level of each gene to discriminate robust splicing events from potential biological noise. We found that rarely used splice sites disproportionately come from highly expressed genes and are significantly less conserved in other nematode genomes than splice sites with a higher usage frequency. Our increased detection power confirmed trans-splicing for at least 84% of C. elegans protein coding genes. The genes for which trans-splicing was not observed are overwhelmingly low expression genes, suggesting that the mechanism is pervasive but not fully captured by organism-wide RNA-seq. We generated annotated gene models including quantitative exon usage information for the entire C. elegans genome. This allows users to visualize at a glance the relative expression of each isoform for their gene of interest.

The Deep Thioredoxome in Chlamydomonas reinhardtii: New Insights into Redox Regulation.

Thiol-based redox post-translational modifications have emerged as important mechanisms of signaling and regulation in all organisms, and thioredoxin plays a key role by controlling the thiol-disulfide status of target proteins. Recent redox proteomic studies revealed hundreds of proteins regulated by glutathionylation and nitrosylation in the unicellular green alga Chlamydomonas reinhardtii, while much less is known about the thioredoxin interactome in this organism. By combining qualitative and quantitative proteomic analyses, we have comprehensively investigated the Chlamydomonas thioredoxome and 1188 targets have been identified. They participate in a wide range of metabolic pathways and cellular processes. This study broadens not only the redox regulation to new enzymes involved in well-known thioredoxin-regulated metabolic pathways but also sheds light on cellular processes for which data supporting redox regulation are scarce (aromatic amino acid biosynthesis, nuclear transport, etc). Moreover, we characterized 1052 thioredoxin-dependent regulatory sites and showed that these data constitute a valuable resource for future functional studies in Chlamydomonas. By comparing this thioredoxome with proteomic data for glutathionylation and nitrosylation at the protein and cysteine levels, this work confirms the existence of a complex redox regulation network in Chlamydomonas and provides evidence of a tremendous selectivity of redox post-translational modifications for specific cysteine residues.

Analysis of the lipid body proteome of the oleaginous alga Lobosphaera incisa.

BACKGROUND: Lobosphaera incisa (L. incisa) is an oleaginous microalga that stores triacylglycerol (TAG) rich in arachidonic acid in lipid bodies (LBs). This organelle is gaining attention in algal research, since evidence is accumulating that proteins attached to its surface fulfill important functions in TAG storage and metabolism. RESULTS: Here, the composition of the LB proteome in L incisa was investigated by comparing different cell fractions in a semiquantitative proteomics approach. After applying stringent filters to the proteomics data in order to remove contaminating proteins from the list of possible LB proteins (LBPs), heterologous expression of candidate proteins in tobacco pollen tubes, allowed us to confirm 3 true LBPs: A member of the algal Major Lipid Droplet Protein family, a small protein of unknown function and a putative lipase. In addition, a TAG lipase that belongs to the SUGAR DEPENDENT 1 family of TAG lipases known from oilseed plants was identified. Its activity was verified by functional complementation of an Arabidopsis thaliana mutant lacking the major seed TAG lipases. CONCLUSIONS: Here we describe 3 LBPs as well as a TAG lipase from the oleaginous microalga L. incisa and discuss their possible involvement in LB metabolism. This study highlights the importance of filtering LB proteome datasets and verifying the subcellular localization one by one, so that contaminating proteins can be recognized as such. Our dataset can serve as a valuable resource in the identification of additional LBPs, shedding more light on the intriguing roles of LBs in microalgae.

The putative drug efflux systems of the Bacillus cereus group.

The Bacillus cereus group of bacteria includes seven closely related species, three of which, B. anthracis, B. cereus and B. thuringiensis, are pathogens of humans, animals and/or insects. Preliminary investigations into the transport capabilities of different bacterial lineages suggested that genes encoding putative efflux systems were unusually abundant in the B. cereus group compared to other bacteria. To explore the drug efflux potential of the B. cereus group all putative efflux systems were identified in the genomes of prototypical strains of B. cereus, B. anthracis and B. thuringiensis using our Transporter Automated Annotation Pipeline. More than 90 putative drug efflux systems were found within each of these strains, accounting for up to 2.7% of their protein coding potential. Comparative analyses demonstrated that the efflux systems are highly conserved between these species; 70-80% of the putative efflux pumps were shared between all three strains studied. Furthermore, 82% of the putative efflux system proteins encoded by the prototypical B. cereus strain ATCC 14579 (type strain) were found to be conserved in at least 80% of 169 B. cereus group strains that have high quality genome sequences available. However, only a handful of these efflux pumps have been functionally characterized. Deletion of individual efflux pump genes from B. cereus typically had little impact to drug resistance phenotypes or the general fitness of the strains, possibly because of the large numbers of alternative efflux systems that may have overlapping substrate specificities. Therefore, to gain insight into the possible transport functions of efflux systems in B. cereus, we undertook large-scale qRT-PCR analyses of efflux pump gene expression following drug shocks and other stress treatments. Clustering of gene expression changes identified several groups of similarly regulated systems that may have overlapping drug resistance functions. In this article we review current knowledge of the small molecule efflux pumps encoded by the B. cereus group and suggest the likely functions of numerous uncharacterised pumps.

The complete sequence of the chloroplast genome of the green microalga Lobosphaera (Parietochloris) incisa.

We hereby report the complete chloroplast genome sequence of the green unicellular alga Lobosphaera (Parietochloris) incisa (strain SAG 2468). The genome consists of a circular chromosome of 156,028 bp, which is 72% A-T rich and does not contain a large rRNA-encoding inverted repeat. It is predicted to encode a total of 111 genes including 78 protein-coding, three rRNA, and 30 tRNA genes. The genome sequence also carries a self-splicing group I intron and a group II intron remnant. Overall, the gene and intron content of the L. incisa chloroplast genome is highly similar to that of other species of Trebouxiophyceae. In contrast, the L. incisa chloroplast genome harbors 88 copies of various intergenic dispersed DNA repeat sequences that are all unique to L. incisa.

The complete mitochondrial genome sequence of the green microalga Lobosphaera (Parietochloris) incisa reveals a new type of palindromic repetitive repeat.

BACKGROUND: Lobosphaera incisa, formerly known as Myrmecia incisa and then Parietochloris incisa, is an oleaginous unicellular green alga belonging to the class Trebouxiophyceae (Chlorophyta). It is the richest known plant source of arachidonic acid, an ω-6 poly-unsaturated fatty acid valued by the pharmaceutical and baby-food industries. It is therefore an organism of high biotechnological interest, and we recently reported the sequence of its chloroplast genome. RESULTS: We now report the complete sequence of the mitochondrial genome of L. incisa from high-throughput Illumina short-read sequencing. The circular chromosome of 69,997 bp is predicted to encode a total of 64 genes, some harboring specific self-splicing group I and group II introns. Overall, the gene content is highly similar to that of the mitochondrial genomes of other Trebouxiophyceae, with 34 protein-coding, 3 rRNA, and 27 tRNA genes. Genes are distributed in two clusters located on different DNA strands, a bipartite arrangement that suggests expression from two divergent promoters yielding polycistronic primary transcripts. The L. incisa mitochondrial genome contains families of intergenic dispersed DNA repeat sequences that are not shared with other known mitochondrial genomes of Trebouxiophyceae. The most peculiar feature of the genome is a repetitive palindromic repeat, the LIMP (L. Incisa Mitochondrial Palindrome), found 19 times in the genome. It is formed by repetitions of an AACCA pentanucleotide, followed by an invariant 7-nt loop and a complementary repeat of the TGGTT motif. Analysis of the genome sequencing reads indicates that the LIMP can be a substrate for large-scale genomic rearrangements. We speculate that LIMPs can act as origins of replication. Deep sequencing of the L. incisa transcriptome also suggests that the LIMPs with long stems are sites of transcript processing. The genome also contains five copies of a related palindromic repeat, the HyLIMP, with a 10-nt motif related to that of the LIMP. CONCLUSIONS: The mitochondrial genome of L. incisa encodes a unique type of repetitive palindromic repeat sequence, the LIMP, which can mediate genome rearrangements and play a role in mitochondrial gene expression. Experimental studies are needed to confirm and further characterize the functional role(s) of the LIMP.