Anaeramoebae are a divergent lineage of eukaryotes that shed light on the transition from anaerobic mitochondria to hydrogenosomes

Discoveries of diverse microbial eukaryotes and their inclusion in comprehensive phylogenomic analyses have crucially re-shaped the eukaryotic tree of life in the 21st century.¹ At the deepest level, eukaryotic diversity comprises 9-10 "supergroups." One of these supergroups, the Metamonada, is particularly important to our understanding of the evolutionary dynamics of eukaryotic cells, including the remodeling of mitochondrial function. All metamonads thrive in low-oxygen environments and lack classical aerobic mitochondria, instead possessing mitochondrion-related organelles (MROs) with metabolisms that are adapted to low-oxygen conditions. These MROs lack an organellar genome, do not participate in the Krebs cycle and oxidative phosphorylation,² and often synthesize ATP by substrate-level phosphorylation coupled to hydrogen production.^3,4 The events that occurred during the transition from an oxygen-respiring mitochondrion to a functionally streamlined MRO early in metamonad evolution remain largely unknown. Here, we report transcriptomes of two recently described, enigmatic, anaerobic protists from the genus Anaeramoeba.⁵ Using phylogenomic analysis, we show that these species represent a divergent, phylum-level lineage in the tree of metamonads, emerging as a sister group of the Parabasalia and reordering the deep branching order of the metamonad tree. Metabolic reconstructions of the Anaeramoeba MROs reveal many "classical" mitochondrial features previously not seen in metamonads, including a disulfide relay import system, propionate production, and amino acid metabolism. Our findings suggest that the cenancestor of Metamonada likely had MROs with more classical mitochondrial features than previously anticipated and demonstrate how discoveries of novel lineages of high taxonomic rank continue to transform our understanding of early eukaryote evolution.

The Mastigamoeba balamuthi Genome and the Nature of the Free-Living Ancestor of Entamoeba

The transition of free-living organisms to parasitic organisms is a mysterious process that occurs in all major eukaryotic lineages. Parasites display seemingly unique features associated with their pathogenicity; however, it is important to distinguish ancestral preconditions to parasitism from truly new parasite-specific functions. Here, we sequenced the genome and transcriptome of anaerobic free-living Mastigamoeba balamuthi and performed phylogenomic analysis of four related members of the Archamoebae, including Entamoeba histolytica, an important intestinal pathogen of humans. We aimed to trace gene histories throughout the adaptation of the aerobic ancestor of Archamoebae to anaerobiosis and throughout the transition from a free-living to a parasitic lifestyle. These events were associated with massive gene losses that, in parasitic lineages, resulted in a reduction in structural features, complete losses of some metabolic pathways, and a reduction in metabolic complexity. By reconstructing the features of the common ancestor of Archamoebae, we estimated preconditions for the evolution of parasitism in this lineage. The ancestor could apparently form chitinous cysts, possessed proteolytic enzyme machinery, compartmentalized the sulfate activation pathway in mitochondrion-related organelles, and possessed the components for anaerobic energy metabolism. After the split of Entamoebidae, this lineage gained genes encoding surface membrane proteins that are involved in host–parasite interactions. In contrast, gene gains identified in the M. balamuthi lineage were predominantly associated with polysaccharide catabolic processes. A phylogenetic analysis of acquired genes suggested an essential role of lateral gene transfer in parasite evolution (Entamoeba) and in adaptation to anaerobic aquatic sediments (Mastigamoeba).

Analysis of diverse eukaryotes suggests the existence of an ancestral mitochondrial apparatus derived from the bacterial type II secretion system

The type 2 secretion system (T2SS) is present in some Gram-negative eubacteria and used to secrete proteins across the outer membrane. Here we report that certain representative heteroloboseans, jakobids, malawimonads and hemimastigotes unexpectedly possess homologues of core T2SS components. We show that at least some of them are present in mitochondria, and their behaviour in biochemical assays is consistent with the presence of a mitochondrial T2SS-derived system (miT2SS). We additionally identified 23 protein families co-occurring with miT2SS in eukaryotes. Seven of these proteins could be directly linked to the core miT2SS by functional data and/or sequence features, whereas others may represent different parts of a broader functional pathway, possibly also involving the peroxisome. Its distribution in eukaryotes and phylogenetic evidence together indicate that the miT2SS-centred pathway is an ancestral eukaryotic trait. Our findings thus have direct implications for the functional properties of the early mitochondrion.

Bacteria use the type 2 secretion system to secrete enzymes and toxins across the outer membrane to the environment. Here the authors analyse the T2SS pathway in three protist lineages and suggest that the early mitochondrion may have been capable of secreting proteins into the cytosol.

The draft nuclear genome sequence and predicted mitochondrial proteome of Andalucia godoyi, a protist with the most gene-rich and bacteria-like mitochondrial genome

Background

Comparative analyses have indicated that the mitochondrion of the last eukaryotic common ancestor likely possessed all the key core structures and functions that are widely conserved throughout the domain Eucarya. To date, such studies have largely focused on animals, fungi, and land plants (primarily multicellular eukaryotes); relatively few mitochondrial proteomes from protists (primarily unicellular eukaryotic microbes) have been examined. To gauge the full extent of mitochondrial structural and functional complexity and to identify potential evolutionary trends in mitochondrial proteomes, more comprehensive explorations of phylogenetically diverse mitochondrial proteomes are required. In this regard, a key group is the jakobids, a clade of protists belonging to the eukaryotic supergroup Discoba, distinguished by having the most gene-rich and most bacteria-like mitochondrial genomes discovered to date.

Results

In this study, we assembled the draft nuclear genome sequence for the jakobid Andalucia godoyi and used a comprehensive in silico approach to infer the nucleus-encoded portion of the mitochondrial proteome of this protist, identifying 864 candidate mitochondrial proteins. The A. godoyi mitochondrial proteome has a complexity that parallels that of other eukaryotes, while exhibiting an unusually large number of ancestral features that have been lost particularly in opisthokont (animal and fungal) mitochondria. Notably, we find no evidence that the A. godoyi nuclear genome has or had a gene encoding a single-subunit, T3/T7 bacteriophage-like RNA polymerase, which functions as the mitochondrial transcriptase in all eukaryotes except the jakobids.

Conclusions

As genome and mitochondrial proteome data have become more widely available, a strikingly punctuate phylogenetic distribution of different mitochondrial components has been revealed, emphasizing that the pathways of mitochondrial proteome evolution are likely complex and lineage-specific. Unraveling this complexity will require comprehensive comparative analyses of mitochondrial proteomes from a phylogenetically broad range of eukaryotes, especially protists. The systematic in silico approach described here offers a valuable adjunct to direct proteomic analysis (e.g., via mass spectrometry), particularly in cases where the latter approach is constrained by sample limitation or other practical considerations.

Microenvironment‑driven resistance to B‑Raf inhibition in a melanoma patient is accompanied by broad changes of gene methylation and expression in distal fibroblasts

The incidence of malignant melanoma is rapidly increasing and current medicine is offering only limited options for treatment of the advanced disease. For B‑Raf mutated melanomas, treatment with mutation‑specific drug inhibitors may be used. Unfortunately, tumors frequently acquire resistance to the treatment. Tumor microenvironment, namely cancer‑associated fibroblasts, largely influence this acquired resistance. In the present study, fibroblasts were isolated from a patient suffering from acrolentiginous melanoma (Breslow, 4.0 mm; Clark, IV; B‑Raf V600E mutated). The present study focused on the expression of structural and functional markers of fibroblast activation in melanoma‑associated fibroblasts (MAFs; isolated prior to therapy initiation) as well as in autologous control fibroblasts (ACFs) of the same patient isolated during B‑Raf inhibitor therapy, yet before clinical progression of the disease. Analysis of gene transcription was also performed, as well as DNA methylation status analysis at the genomic scale of both isolates. MAFs were positive for smooth muscle actin (SMA), which is a marker of myofibroblasts and the hallmark of cancer stoma. Surprisingly, ACF isolated from the distant uninvolved skin of the same patient also exhibited strong SMA expression. A similar phenotype was also observed in control dermal fibroblasts (CDFs; from different donors) exclusively following stimulation by transforming growth factor (TGF)‑β1. Immunohistochemistry confirmed that melanoma cells potently produce TGF‑β1. Significant differences were also identified in gene transcription and in DNA methylation status at the genomic scale. Upregulation of SMA was observed in ACF cells at the protein and transcriptional levels. The present results support recent experimental findings that tumor microenvironment is driving resistance to B‑Raf inhibition in patients with melanoma. Such an activated microenvironment may be viable for the growth of circulating melanoma cells.

FMNH2-dependent monooxygenases initiate catabolism of sulfonamides in Microbacterium sp. strain BR1 subsisting on sulfonamide antibiotics

We report a cluster of genes encoding two monooxygenases (SadA and SadB) and one FMN reductase (SadC) that enable Microbacterium sp. strain BR1 and other Actinomycetes to inactivate sulfonamide antibiotics. Our results show that SadA and SadC are responsible for the initial attack of sulfonamide molecules resulting in the release of 4-aminophenol. The latter is further transformed into 1,2,4-trihydroxybenzene by SadB and SadC prior to mineralization and concomitant production of biomass. As the degradation products lack antibiotic activity, the presence of SadA will result in an alleviated bacteriostatic effect of sulfonamides. In addition to the relief from antibiotic stress this bacterium gains access to an additional carbon source when this gene cluster is expressed. As degradation of sulfonamides was also observed when Microbacterium sp. strain BR1 was grown on artificial urine medium, colonization with such strains may impede common sulfonamide treatment during co-infections with pathogens of the urinary tract. This case of biodegradation exemplifies the evolving catabolic capacity of bacteria, given that sulfonamide bacteriostatic are purely of synthetic origin. The wide distribution of this cluster in Actinomycetes and the presence of traA encoding a relaxase in its vicinity suggest that this cluster is mobile and that is rather alarming.

Mitochondrial Genome of Prasinophyte Alga Pyramimonas parkeae

Prasinophytes are a paraphyletic assemblage of nine heterogeneous lineages in the Chlorophyta clade of Archaeplastida. Until now, seven complete mitochondrial genomes have been sequenced from four prasinophyte lineages. Here, we report the mitochondrial genome of Pyramimonas parkeae, the first representative of the prasinophyte clade I. The circular-mapping molecule is 43,294 bp long, AT rich (68.8%), very compact and it comprises two 6,671 bp long inverted repeat regions. The gene content is slightly smaller than the gene-richest prasinophyte mitochondrial genomes. The single identified intron is located in the cytochrome c oxidase subunit 1 gene (cox1). Interestingly, two exons of cox1 are encoded on the same strand of DNA in the reverse order and the mature mRNA is formed by trans-splicing. The phylogenetic analysis using the data set of 6,037 positions assembled from 34 mtDNA-encoded proteins of 48 green algae and plants is not in compliance with the branching order of prasinophyte clades revealed on the basis of 18S rRNA genes and cpDNA-encoded proteins. However, the phylogenetic analyses based on all three genomic elements support the sister position of prasinophyte clades Pyramimonadales and Mamiellales.

A Eukaryote without a Mitochondrial Organelle

The presence of mitochondria and related organelles in every studied eukaryote supports the view that mitochondria are essential cellular components. Here, we report the genome sequence of a microbial eukaryote, the oxymonad Monocercomonoides sp., which revealed that this organism lacks all hallmark mitochondrial proteins. Crucially, the mitochondrial iron-sulfur cluster assembly pathway, thought to be conserved in virtually all eukaryotic cells, has been replaced by a cytosolic sulfur mobilization system (SUF) acquired by lateral gene transfer from bacteria. In the context of eukaryotic phylogeny, our data suggest that Monocercomonoides is not primitively amitochondrial but has lost the mitochondrion secondarily. This is the first example of a eukaryote lacking any form of a mitochondrion, demonstrating that this organelle is not absolutely essential for the viability of a eukaryotic cell.

A Comparative Analysis of Mitochondrial Genomes in Eustigmatophyte Algae

Eustigmatophyceae (Ochrophyta, Stramenopiles) is a small algal group with species of the genus Nannochloropsis being its best studied representatives. Nuclear and organellar genomes have been recently sequenced for several Nannochloropsis spp., but phylogenetically wider genomic studies are missing for eustigmatophytes. We sequenced mitochondrial genomes (mitogenomes) of three species representing most major eustigmatophyte lineages, Monodopsis sp. MarTras21, Vischeria sp. CAUP Q 202 and Trachydiscus minutus, and carried out their comparative analysis in the context of available data from Nannochloropsis and other stramenopiles, revealing a number of noticeable findings. First, mitogenomes of most eustigmatophytes are highly collinear and similar in the gene content, but extensive rearrangements and loss of three otherwise ubiquitous genes happened in the Vischeria lineage; this correlates with an accelerated evolution of mitochondrial gene sequences in this lineage. Second, eustigmatophytes appear to be the only ochrophyte group with the Atp1 protein encoded by the mitogenome. Third, eustigmatophyte mitogenomes uniquely share a truncated nad11 gene encoding only the C-terminal part of the Nad11 protein, while the N-terminal part is encoded by a separate gene in the nuclear genome. Fourth, UGA as a termination codon and the cognate release factor mRF2 were lost from mitochondria independently by the Nannochloropsis and T. minutus lineages. Finally, the rps3 gene in the mitogenome of Vischeria sp. is interrupted by the UAG codon, but the genome includes a gene for an unusual tRNA with an extended anticodon loop that we speculate may serve as a suppressor tRNA to properly decode the rps3 gene.

First multigene analysis of Archamoebae (Amoebozoa: Conosa) robustly reveals its phylogeny and shows that Entamoebidae represents a deep lineage of the group

Archamoebae is an understudied group of anaerobic free-living or endobiotic protists that constitutes the major anaerobic lineage of the supergroup Amoebozoa. Hitherto, the phylogeny of Archamoebae was based solely on SSU rRNA and actin genes, which did not resolve relationships among the main lineages of the group. Because of this uncertainty, several different scenarios had been proposed for the phylogeny of the Archamoebae. In this study, we present the first multigene phylogenetic analysis that includes members of Pelomyxidae, and Rhizomastixidae. The analysis clearly shows that Mastigamoebidae, Pelomyxidae and Rhizomastixidae form a clade of mostly free-living, amoeboid flagellates, here called Pelobiontida. The predominantly endobiotic and aflagellated Entamoebidae represents a separate, deep-branching lineage, Entamoebida. Therefore, two unique evolutionary events, horizontal transfer of the nitrogen fixation system from bacteria and transfer of the sulfate activation pathway to mitochondrial derivatives, predate the radiation of recent lineages of Archamoebae. The endobiotic lifestyle has arisen at least three times independently during the evolution of the group. We also present new ultrastructural data that clarifies the primary divergence among the family Mastigamoebidae which had previously been inferred from phylogenetic analyses based on SSU rDNA.

Combined Culture-Based and Culture-Independent Approaches Provide Insights into Diversity of Jakobids, an Extremely Plesiomorphic Eukaryotic Lineage

We used culture-based and culture-independent approaches to discover diversity and ecology of anaerobic jakobids (Excavata: Jakobida), an overlooked, deep-branching lineage of free-living nanoflagellates related to Euglenozoa. Jakobids are among a few lineages of nanoflagellates frequently detected in anoxic habitats by PCR-based studies, however only two strains of a single jakobid species have been isolated from those habitats. We recovered 712 environmental sequences and cultured 21 new isolates of anaerobic jakobids that collectively represent at least ten different species in total, from which four are uncultured. Two cultured species have never been detected by environmental, PCR-based methods. Surprisingly, culture-based and culture-independent approaches were able to reveal a relatively high proportion of overall species diversity of anaerobic jakobids—60 or 80%, respectively. Our phylogenetic analyses based on SSU rDNA and six protein-coding genes showed that anaerobic jakobids constitute a clade of morphologically similar, but genetically and ecologically diverse protists—Stygiellidae fam. nov. Our investigation combines culture-based and environmental molecular-based approaches to capture a wider extent of species diversity and shows Stygiellidae as a group that ordinarily inhabits anoxic, sulfide- and ammonium-rich marine habitats worldwide.

Cubozoan genome illuminates functional diversification of opsins and photoreceptor evolution

Animals sense light primarily by an opsin-based photopigment present in a photoreceptor cell. Cnidaria are arguably the most basal phylum containing a well-developed visual system. The evolutionary history of opsins in the animal kingdom has not yet been resolved. Here, we study the evolution of animal opsins by genome-wide analysis of the cubozoan jellyfish Tripedalia cystophora, a cnidarian possessing complex lens-containing eyes and minor photoreceptors. A large number of opsin genes with distinct tissue- and stage-specific expression were identified. Our phylogenetic analysis unequivocally classifies cubozoan opsins as a sister group to c-opsins and documents lineage-specific expansion of the opsin gene repertoire in the cubozoan genome. Functional analyses provided evidence for the use of the Gs-cAMP signaling pathway in a small set of cubozoan opsins, indicating the possibility that the majority of other cubozoan opsins signal via distinct pathways. Additionally, these tests uncovered subtle differences among individual opsins, suggesting possible fine-tuning for specific photoreceptor tasks. Based on phylogenetic, expression and biochemical analysis we propose that rapid lineage- and species-specific duplications of the intron-less opsin genes and their subsequent functional diversification promoted evolution of a large repertoire of both visual and extraocular photoreceptors in cubozoans.

Updating algal evolutionary relationships through plastid genome sequencing: did alveolate plastids emerge through endosymbiosis of an ochrophyte?

Algae with secondary plastids of a red algal origin, such as ochrophytes (photosynthetic stramenopiles), are diverse and ecologically important, yet their evolutionary history remains controversial. We sequenced plastid genomes of two ochrophytes, Ochromonas sp. CCMP1393 (Chrysophyceae) and Trachydiscus minutus (Eustigmatophyceae). A shared split of the clpC gene as well as phylogenomic analyses of concatenated protein sequences demonstrated that chrysophytes and eustigmatophytes form a clade, the Limnista, exhibiting an unexpectedly elevated rate of plastid gene evolution. Our analyses also indicate that the root of the ochrophyte phylogeny falls between the recently redefined Khakista and Phaeista assemblages. Taking advantage of the expanded sampling of plastid genome sequences, we revisited the phylogenetic position of the plastid of Vitrella brassicaformis, a member of Alveolata with the least derived plastid genome known for the whole group. The results varied depending on the dataset and phylogenetic method employed, but suggested that the Vitrella plastids emerged from a deep ochrophyte lineage rather than being derived vertically from a hypothetical plastid-bearing common ancestor of alveolates and stramenopiles. Thus, we hypothesize that the plastid in Vitrella, and potentially in other alveolates, may have been acquired by an endosymbiosis of an early ochrophyte.

Melanoma cells influence the differentiation pattern of human epidermal keratinocytes

Background

Nodular melanoma is one of the most life threatening tumors with still poor therapeutic outcome. Similarly to other tumors, permissive microenvironment is essential for melanoma progression. Features of this microenvironment are arising from molecular crosstalk between the melanoma cells (MC) and the surrounding cell populations in the context of skin tissue. Here, we study the effect of melanoma cells on human primary keratinocytes (HPK). Presence of MC is as an important modulator of the tumor microenvironment and we compare it to the effect of nonmalignant lowly differentiated cells also originating from neural crest (NCSC).

Methods

Comparative morphometrical and immunohistochemical analysis of epidermis surrounding nodular melanoma (n = 100) was performed. Data were compared to results of transcriptome profiling of in vitro models, in which HPK were co-cultured with MC, normal human melanocytes, and NCSC, respectively. Differentially expressed candidate genes were verified by RT-qPCR. Biological activity of candidate proteins was assessed on cultured HPK.

Results

Epidermis surrounding nodular melanoma exhibits hyperplastic features in 90% of cases. This hyperplastic region exhibits aberrant suprabasal expression of keratin 14 accompanied by loss of keratin 10. We observe that MC and NCSC are able to increase expression of keratins 8, 14, 19, and vimentin in the co-cultured HPK. This in vitro finding partially correlates with pseudoepitheliomatous hyperplasia observed in melanoma biopsies. We provide evidence of FGF-2, CXCL-1, IL-8, and VEGF-A participation in the activity of melanoma cells on keratinocytes.

Conclusion

We conclude that the MC are able to influence locally the differentiation pattern of keratinocytes in vivo as well as in vitro. This interaction further highlights the role of intercellular interactions in melanoma. The reciprocal role of activated keratinocytes on biology of melanoma cells shall be verified in the future.

Electronic supplementary material

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

The mitochondrion-like organelle of Trimastix pyriformis contains the complete glycine cleavage system

All eukaryotic organisms contain mitochondria or organelles that evolved from the same endosymbiotic event like classical mitochondria. Organisms inhabiting low oxygen environments often contain mitochondrial derivates known as hydrogenosomes, mitosomes or neutrally as mitochondrion-like organelles. The detailed investigation has shown unexpected evolutionary plasticity in the biochemistry and protein composition of these organelles in various protists. We investigated the mitochondrion-like organelle in Trimastix pyriformis, a free-living member of one of the three lineages of anaerobic group Metamonada. Using 454 sequencing we have obtained 7 037 contigs from its transcriptome and on the basis of sequence homology and presence of N-terminal extensions we have selected contigs coding for proteins that putatively function in the organelle. Together with the results of a previous transcriptome survey, the list now consists of 23 proteins - mostly enzymes involved in amino acid metabolism, transporters and maturases of proteins and transporters of metabolites. We have no evidence of the production of ATP in the mitochondrion-like organelle of Trimastix but we have obtained experimental evidence for the presence of enzymes of the glycine cleavage system (GCS), which is part of amino acid metabolism. Using homologous antibody we have shown that H-protein of GCS localizes into vesicles in the cell of Trimastix. When overexpressed in yeast, H- and P-protein of GCS and cpn60 were transported into mitochondrion. In case of H-protein we have demonstrated that the first 16 amino acids are necessary for this transport. Glycine cleavage system is at the moment the only experimentally localized pathway in the mitochondrial derivate of Trimastix pyriformis.

Intraspecific sequence comparisons reveal similar rates of non-collinear gene insertion in the B and D genomes of bread wheat

BACKGROUND

Polyploidization is considered one of the main mechanisms of plant genome evolution. The presence of multiple copies of the same gene reduces selection pressure and permits sub-functionalization and neo-functionalization leading to plant diversification, adaptation and speciation. In bread wheat, polyploidization and the prevalence of transposable elements resulted in massive gene duplication and movement. As a result, the number of genes which are non-collinear to genomes of related species seems markedly increased in wheat.

RESULTS

We used new-generation sequencing (NGS) to generate sequence of a Mb-sized region from wheat chromosome arm 3DS. Sequence assembly of 24 BAC clones resulted in two scaffolds of 1,264,820 and 333,768 bases. The sequence was annotated and compared to the homoeologous region on wheat chromosome 3B and orthologous loci of Brachypodium distachyon and rice. Among 39 coding sequences in the 3DS scaffolds, 32 have a homoeolog on chromosome 3B. In contrast, only fifteen and fourteen orthologs were identified in the corresponding regions in rice and Brachypodium, respectively. Interestingly, five pseudogenes were identified among the non-collinear coding sequences at the 3B locus, while none was found at the 3DS locus.

CONCLUSION

Direct comparison of two Mb-sized regions of the B and D genomes of bread wheat revealed similar rates of non-collinear gene insertion in both genomes with a majority of gene duplications occurring before their divergence. Relatively low proportion of pseudogenes was identified among non-collinear coding sequences. Our data suggest that the pseudogenes did not originate from insertion of non-functional copies, but were formed later during the evolution of hexaploid wheat. Some evidence was found for gene erosion along the B genome locus.

The plastid genome of Eutreptiella provides a window into the process of secondary endosymbiosis of plastid in euglenids

Euglenids are a group of protists that comprises species with diverse feeding modes. One distinct and diversified clade of euglenids is photoautotrophic, and its members bear green secondary plastids. In this paper we present the plastid genome of the euglenid Eutreptiella, which we assembled from 454 sequencing of Eutreptiella gDNA. Comparison of this genome and the only other available plastid genomes of photosynthetic euglenid, Euglena gracilis, revealed that they contain a virtually identical set of 57 protein coding genes, 24 genes fewer than the genome of Pyramimonas parkeae, the closest extant algal relative of the euglenid plastid. Searching within the transcriptomes of Euglena and Eutreptiella showed that 6 of the missing genes were transferred to the nucleus of the euglenid host while 18 have been probably lost completely. Euglena and Eutreptiella represent the deepest bifurcation in the photosynthetic clade, and therefore all these gene transfers and losses must have happened before the last common ancestor of all known photosynthetic euglenids. After the split of Euglena and Eutreptiella only one additional gene loss took place. The conservation of gene content in the two lineages of euglenids is in contrast to the variability of gene order and intron counts, which diversified dramatically. Our results show that the early secondary plastid of euglenids was much more susceptible to gene losses and endosymbiotic gene transfers than the established plastid, which is surprisingly resistant to changes in gene content.

Global transcriptome analysis of the C57BL/6J mouse testis by SAGE: evidence for nonrandom gene order

BACKGROUND

We generated the gene expression profile of the total testis from the adult C57BL/6J male mice using serial analysis of gene expression (SAGE). Two high-quality SAGE libraries containing a total of 76 854 tags were constructed. An extensive bioinformatic analysis and comparison of SAGE transcriptomes of the total testis, testicular somatic cells and other mouse tissues was performed and the theory of male-biased gene accumulation on the X chromosome was tested.

RESULTS

We sorted out 829 genes predominantly expressed from the germinal part and 944 genes from the somatic part of the testis. The genes preferentially and specifically expressed in total testis and testicular somatic cells were identified by comparing the testis SAGE transcriptomes to the available transcriptomes of seven non-testis tissues. We uncovered chromosomal clusters of adjacent genes with preferential expression in total testis and testicular somatic cells by a genome-wide search and found that the clusters encompassed a significantly higher number of genes than expected by chance. We observed a significant 3.2-fold enrichment of the proportion of X-linked genes specific for testicular somatic cells, while the proportions of X-linked genes specific for total testis and for other tissues were comparable. In contrast to the tissue-specific genes, an under-representation of X-linked genes in the total testis transcriptome but not in the transcriptomes of testicular somatic cells and other tissues was detected.

CONCLUSION

Our results provide new evidence in favor of the theory of male-biased genes accumulation on the X chromosome in testicular somatic cells and indicate the opposite action of the meiotic X-inactivation in testicular germ cells.