Isabel Duarte

The Organellar Genome and Metabolic Potential of the Hydrogen-Producing Mitochondrion of Nyctotherus ovalis

It is generally accepted that hydrogenosomes (hydrogen-producing organelles) evolved from a mitochondrial ancestor. However, until recently, only indirect evidence for this hypothesis was available. Here, we present the almost complete genome of the hydrogen-producing mitochondrion of the anaerobic ciliate Nyctotherus ovalis and show that, except for the notable absence of genes encoding electron transport chain components of Complexes III, IV, and V, it has a gene content similar to the mitochondrial genomes of aerobic ciliates. Analysis of the genome of the hydrogen-producing mitochondrion, in combination with that of more than 9,000 genomic DNA and cDNA sequences, allows a preliminary reconstruction of the organellar metabolism. The sequence data indicate that N. ovalis possesses hydrogen-producing mitochondria that have a truncated, two step (Complex I and II) electron transport chain that uses fumarate as electron acceptor. In addition, components of an extensive protein network for the metabolism of amino acids, defense against oxidative stress, mitochondrial protein synthesis, mitochondrial protein import and processing, and transport of metabolites across the mitochondrial membrane were identified. Genes for MPV17 and ACN9, two hypothetical proteins linked to mitochondrial disease in humans, were also found. The inferred metabolism is remarkably similar to the organellar metabolism of the phylogenetically distant anaerobic Stramenopile Blastocystis. Notably, the Blastocystis organelle and that of the related flagellate Proteromonas lacertae also lack genes encoding components of Complexes III, IV, and V. Thus, our data show that the hydrogenosomes of N. ovalis are highly specialized hydrogen-producing mitochondria.

BOLA1 is an aerobic protein that prevents mitochondrial morphology changes induced by glutathione depletion

AIMS The BolA protein family is widespread among eukaryotes and bacteria. In Escherichia coli, BolA causes a spherical cell shape and is overexpressed during oxidative stress. Here we aim to elucidate the possible role of its human homolog BOLA1 in mitochondrial morphology and thiol redox potential regulation. RESULTS We show that BOLA1 is a mitochondrial protein that counterbalances the effect of L-buthionine-(S,R)-sulfoximine (BSO)-induced glutathione (GSH) depletion on the mitochondrial thiol redox potential. Furthermore, overexpression of BOLA1 nullifies the effect of BSO and S-nitrosocysteine on mitochondrial morphology. Conversely, knockdown of the BOLA1 gene increases the oxidation of mitochondrial thiol groups. Supporting a role of BOLA1 in controlling the mitochondrial thiol redox potential is that BOLA1 orthologs only occur in aerobic eukaryotes. A measured interaction of BOLA1 with the mitochondrial monothiol glutaredoxin GLRX5 provides hints for potential mechanisms behind BOLA1s effect on mitochondrial redox potential. Nevertheless, we have no direct evidence for a role of GLRX5 in BOLA1s function. INNOVATION We implicate a new protein, BOLA1, in the regulation of the mitochondrial thiol redox potential. CONCLUSION BOLA1 is an aerobic, mitochondrial protein that prevents mitochondrial morphology aberrations induced by GSH depletion and reduces the associated oxidative shift of the mitochondrial thiol redox potential.

Evolution and Diversification of the Organellar Release Factor Family

Translation termination is accomplished by proteins of the Class I release factor family (RF) that recognize stop codons and catalyze the ribosomal release of the newly synthesized peptide. Bacteria have two canonical RFs: RF1 recognizes UAA and UAG, RF2 recognizes UAA and UGA. Despite that these two release factor proteins are sufficient for de facto translation termination, the eukaryotic organellar RF protein family, which has evolved from bacterial release factors, has expanded considerably, comprising multiple subfamilies, most of which have not been functionally characterized or formally classified. Here, we integrate multiple sources of information to analyze the remarkable differentiation of the RF family among organelles. We document the origin, phylogenetic distribution and sequence structure features of the mitochondrial and plastidial release factors: mtRF1a, mtRF1, mtRF2a, mtRF2b, mtRF2c, ICT1, C12orf65, pRF1, and pRF2, and review published relevant experimental data. The canonical release factors (mtRF1a, mtRF2a, pRF1, and pRF2) and ICT1 are derived from bacterial ancestors, whereas the others have resulted from gene duplications of another release factor. These new RF family members have all lost one or more specific motifs relevant for bona fide release factor function but are mostly targeted to the same organelle as their ancestor. We also characterize the subset of canonical release factor proteins that bear nonclassical PxT/SPF tripeptide motifs and provide a molecular-model-based rationale for their retained ability to recognize stop codons. Finally, we analyze the coevolution of canonical RFs with the organellar genetic code. Although the RF presence in an organelle and its stop codon usage tend to coevolve, we find three taxa that encode an RF2 without using UGA stop codons, and one reverse scenario, where mamiellales green algae use UGA stop codons in their mitochondria without having a mitochondrial type RF2. For the latter, we put forward a “stop-codon reinvention” hypothesis that involves the retargeting of the plastid release factor to the mitochondrion.

Macronuclear genome structure of the ciliate Nyctotherus ovalis : Single-gene chromosomes and tiny introns

BackgroundNyctotherus ovalis is a single-celled eukaryote that has hydrogen-producing mitochondria and lives in the hindgut of cockroaches. Like all members of the ciliate taxon, it has two types of nuclei, a micronucleus and a macronucleus. N. ovalis generates its macronuclear chromosomes by forming polytene chromosomes that subsequently develop into macronuclear chromosomes by DNA elimination and rearrangement.ResultsWe examined the structure of these gene-sized macronuclear chromosomes in N. ovalis. We determined the telomeres, subtelomeric regions, UTRs, coding regions and introns by sequencing a large set of macronuclear DNA sequences (4,242) and cDNAs (5,484) and comparing them with each other. The telomeres consist of repeats CCC(AAAACCCC)n, similar to those in spirotrichous ciliates such as Euplotes, Sterkiella (Oxytricha) and Stylonychia. Per sequenced chromosome we found evidence for either a single protein-coding gene, a single tRNA, or the complete ribosomal RNAs cluster. Hence the chromosomes appear to encode single transcripts. In the short subtelomeric regions we identified a few overrepresented motifs that could be involved in gene regulation, but there is no consensus polyadenylation site. The introns are short (21–29 nucleotides), and a significant fraction (1/3) of the tiny introns is conserved in the distantly related ciliate Paramecium tetraurelia. As has been observed in P. tetraurelia, the N. ovalis introns tend to contain in-frame stop codons or have a length that is not dividable by three. This pattern causes premature termination of mRNA translation in the event of intron retention, and potentially degradation of unspliced mRNAs by the nonsense-mediated mRNA decay pathway.ConclusionThe combination of short leaders, tiny introns and single genes leads to very minimal macronuclear chromosomes. The smallest we identified contained only 150 nucleotides.

Structure based hypothesis of a mitochondrial ribosome rescue mechanism.

BackgroundmtRF1 is a vertebrate mitochondrial protein with an unknown function that arose from a duplication of the mitochondrial release factor mtRF1a. To elucidate the function of mtRF1, we determined the positions that are conserved among mtRF1 sequences but that are different in their mtRF1a paralogs. We subsequently modeled the 3D structure of mtRF1a and mtRF1 bound to the ribosome, highlighting the structural implications of these differences to derive a hypothesis for the function of mtRF1.ResultsOur model predicts, in agreement with the experimental data, that the 3D structure of mtRF1a allows it to recognize the stop codons UAA and UAG in the A-site of the ribosome. In contrast, we show that mtRF1 likely can only bind the ribosome when the A-site is devoid of mRNA. Furthermore, while mtRF1a will adopt its catalytic conformation, in which it functions as a peptidyl-tRNA hydrolase in the ribosome, only upon binding of a stop codon in the A-site, mtRF1 appears specifically adapted to assume this extended, peptidyl-tRNA hydrolyzing conformation in the absence of mRNA in the A-site.ConclusionsWe predict that mtRF1 specifically recognizes ribosomes with an empty A-site and is able to function as a peptidyl-tRNA hydrolase in those situations. Stalled ribosomes with empty A-sites that still contain a tRNA bound to a peptide chain can result from the translation of truncated, stop-codon less mRNAs. We hypothesize that mtRF1 recycles such stalled ribosomes, performing a function that is analogous to that of tmRNA in bacteria.ReviewersThis article was reviewed by Dr. Eugene Koonin, Prof. Knud H. Nierhaus (nominated by Dr. Sarah Teichmann) and Dr. Shamil Sunyaev.

TGFβ-induced switch from adipogenic to osteogenic differentiation of human mesenchymal stem cells: identification of drug targets for prevention of fat cell differentiation

BackgroundPatients suffering from osteoporosis show an increased number of adipocytes in their bone marrow, concomitant with a reduction in the pool of human mesenchymal stem cells (hMSCs) that are able to differentiate into osteoblasts, thus leading to suppressed osteogenesis.MethodsIn order to be able to interfere with this process, we have investigated in-vitro culture conditions whereby adipogenic differentiation of hMSCs is impaired and osteogenic differentiation is promoted. By means of gene expression microarray analysis, we have investigated genes which are potential targets for prevention of fat cell differentiation.ResultsOur data show that BMP2 promotes both adipogenic and osteogenic differentiation of hMSCs, while transforming growth factor beta (TGFβ) inhibits differentiation into both lineages. However, when cells are cultured under adipogenic differentiation conditions, which contain cAMP-enhancing agents such as IBMX of PGE2, TGFβ promotes osteogenic differentiation, while at the same time inhibiting adipogenic differentiation. Gene expression and immunoblot analysis indicated that IBMX-induced suppression of HDAC5 levels plays an important role in the inhibitory effect of TGFβ on osteogenic differentiation. By means of gene expression microarray analysis, we have investigated genes which are downregulated by TGFβ under adipogenic differentiation conditions and may therefore be potential targets for prevention of fat cell differentiation. We thus identified nine genes for which FDA-approved drugs are available. Our results show that drugs directed against the nuclear hormone receptor PPARG, the metalloproteinase ADAMTS5, and the aldo-keto reductase AKR1B10 inhibit adipogenic differentiation in a dose-dependent manner, although in contrast to TGFβ they do not appear to promote osteogenic differentiation.ConclusionsThe approach chosen in this study has resulted in the identification of new targets for inhibition of fat cell differentiation, which may not only be relevant for prevention of osteoporosis, but also of obesity.

StemChecker: a web-based tool to discover and explore stemness signatures in gene sets

Stem cells present unique regenerative abilities, offering great potential for treatment of prevalent pathologies such as diabetes, neurodegenerative and heart diseases. Various research groups dedicated significant effort to identify sets of genes—so-called stemness signatures—considered essential to define stem cells. However, their usage has been hindered by the lack of comprehensive resources and easy-to-use tools. For this we developed StemChecker, a novel stemness analysis tool, based on the curation of nearly fifty published stemness signatures defined by gene expression, RNAi screens, Transcription Factor (TF) binding sites, literature reviews and computational approaches. StemChecker allows researchers to explore the presence of stemness signatures in user-defined gene sets, without carrying-out lengthy literature curation or data processing. To assist in exploring underlying regulatory mechanisms, we collected over 80 target gene sets of TFs associated with pluri- or multipotency. StemChecker presents an intuitive graphical display, as well as detailed statistical results in table format, which helps revealing transcriptionally regulatory programs, indicating the putative involvement of stemness-associated processes in diseases like cancer. Overall, StemChecker substantially expands the available repertoire of online tools, designed to assist the stem cell biology, developmental biology, regenerative medicine and human disease research community. StemChecker is freely accessible at http://stemchecker.sysbiolab.eu.

A three-dimensional topology of complex I inferred from evolutionary correlations

BackgroundThe quaternary structure of eukaryotic NADH:ubiquinone oxidoreductase (complex I), the largest complex of the oxidative phosphorylation, is still mostly unresolved. Furthermore, it is unknown where transiently bound assembly factors interact with complex I. We therefore asked whether the evolution of complex I contains information about its 3D topology and the binding positions of its assembly factors. We approached these questions by correlating the evolutionary rates of eukaryotic complex I subunits using the mirror-tree method and mapping the results into a 3D representation by multidimensional scaling.ResultsMore than 60% of the evolutionary correlation among the conserved seven subunits of the complex I matrix arm can be explained by the physical distance between the subunits. The three-dimensional evolutionary model of the eukaryotic conserved matrix arm has a striking similarity to the matrix arm quaternary structure in the bacterium Thermus thermophilus (rmsd=19 Å) and supports the previous finding that in eukaryotes the N-module is turned relative to the Q-module when compared to bacteria. By contrast, the evolutionary rates contained little information about the structure of the membrane arm. A large evolutionary model of 45 subunits and assembly factors allows to predict subunit positions and interactions (rmsd = 52.6 Å). The model supports an interaction of NDUFAF3, C8orf38 and C2orf56 during the assembly of the proximal matrix arm and the membrane arm. The model further suggests a tight relationship between the assembly factor NUBPL and NDUFA2, which both have been linked to iron-sulfur cluster assembly, as well as between NDUFA12 and its paralog, the assembly factor NDUFAF2.ConclusionsThe physical distance between subunits of complex I is a major correlate of the rate of protein evolution in the complex I matrix arm and is sufficient to infer parts of the complex’s structure with high accuracy. The resulting evolutionary model predicts the positions of a number of subunits and assembly factors.

The intestinal proteome of diabetic and control children is enriched with different microbial and host proteins.

In this study the intestinal microbial proteome of children with established type 1 diabetes (T1D) was compared with the proteome of healthy children (Control) with the aim to identify differences in the activity of the intestinal microbiota that will not only contribute to a deeper knowledge of the functionality of the gut in these children but may also provide new approaches to improve the Control of the disease. Fecal protein extracts collected from three T1D children (aged of 9.3 ± 0.6) and three Control children (aged of 9.3 ± 1.5) were analyzed using a combination of two-dimensional gel electrophoresis (2-DE) and spectral counting. The results evidenced markedly differences between the intestinal proteome of T1D children and the Control. The T1D microbial intestinal proteome was enriched with proteins of clostridial cluster XVa and cluster IV and Bacteroides. In contrast, the Control proteome was enriched with bifidobacterial proteins. In both groups proteins with moonlight function were observed. Human proteins also distinguished the two groups with T1D children depleted in exocrine pancreatic enzymes.In this study, the intestinal microbial proteome of children with established type 1 diabetes (T1D) was compared with the proteome of healthy children (Control) with the aim to identify differences in the activity of the intestinal microbiota that not only will contribute to a deeper knowledge of the functionality of the gut in these children but also may provide new approaches to improve the control of the disease. Faecal protein extracts collected from three T1D children (aged 9.3±0.6 years) and three Control children (aged 9.3±1.5 years) were analysed using a combination of 2D gel electrophoresis and spectral counting. The results evidenced markedly differences between the intestinal proteome of T1D children and the Control. The T1D microbial intestinal proteome was enriched with proteins of clostridial cluster XVa and cluster IV and Bacteroides. In contrast, the Control proteome was enriched with bifidobacterial proteins. In both groups, proteins with moonlight function were observed. Human proteins also distinguished the two groups with T1D children depleted in exocrine pancreatic enzymes.

StemMapper: a curated gene expression database for stem cell lineage analysis.

Transcriptomic data have become a fundamental resource for stem cell (SC) biologists as well as for a wider research audience studying SC-related processes such as aging, embryonic development and prevalent diseases including cancer, diabetes and neurodegenerative diseases. Access and analysis of the growing amount of freely available transcriptomics datasets for SCs, however, are not trivial tasks. Here, we present StemMapper, a manually curated gene expression database and comprehensive resource for SC research, built on integrated data for different lineages of human and mouse SCs. It is based on careful selection, standardized processing and stringent quality control of relevant transcriptomics datasets to minimize artefacts, and includes currently over 960 transcriptomes covering a broad range of SC types. Each of the integrated datasets was individually inspected and manually curated. StemMappers user-friendly interface enables fast querying, comparison, and interactive visualization of quality-controlled SC gene expression data in a comprehensive manner. A proof-of-principle analysis discovering novel putative astrocyte/neural SC lineage markers exemplifies the utility of the integrated data resource. We believe that StemMapper can open the way for new insights and advances in SC research by greatly simplifying the access and analysis of SC transcriptomic data. StemMapper is freely accessible at http://stemmapper.sysbiolab.eu.