Mitochondrial glutamate carriers from Drosophila melanogaster: biochemical, evolutionary and modeling studies.

The mitochondrial carriers are members of a family of transport proteins that mediate solute transport across the inner mitochondrial membrane. Two isoforms of the glutamate carriers, GC1 and GC2 (encoded by the SLC25A22 and SLC25A18 genes, respectively), have been identified in humans. Two independent mutations in SLC25A22 are associated with severe epileptic encephalopathy. In the present study we show that two genes (CG18347 and CG12201) phylogenetically related to the human GC encoding genes are present in the D. melanogaster genome. We have functionally characterized the proteins encoded by CG18347 and CG12201, designated as DmGC1p and DmGC2p respectively, by overexpression in Escherichia coli and reconstitution into liposomes. Their transport properties demonstrate that DmGC1p and DmGC2p both catalyze the transport of glutamate across the inner mitochondrial membrane. Computational approaches have been used in order to highlight residues of DmGC1p and DmGC2p involved in substrate binding. Furthermore, gene expression analysis during development and in various adult tissues reveals that CG18347 is ubiquitously expressed in all examined D. melanogaster tissues, while the expression of CG12201 is strongly testis-biased. Finally, we identified mitochondrial glutamate carrier orthologs in 49 eukaryotic species in order to attempt the reconstruction of the evolutionary history of the glutamate carrier function. Comparison of the exon/intron structure and other key features of the analyzed orthologs suggests that eukaryotic glutamate carrier genes descend from an intron-rich ancestral gene already present in the common ancestor of lineages that diverged as early as bilateria and radiata.

Evolution of genes and genomes on the Drosophila phylogeny.

Comparative analysis of multiple genomes in a phylogenetic framework dramatically improves the precision and sensitivity of evolutionary inference, producing more robust results than single-genome analyses can provide. The genomes of 12 Drosophila species, ten of which are presented here for the first time (sechellia, simulans, yakuba, erecta, ananassae, persimilis, willistoni, mojavensis, virilis and grimshawi), illustrate how rates and patterns of sequence divergence across taxa can illuminate evolutionary processes on a genomic scale. These genome sequences augment the formidable genetic tools that have made Drosophila melanogaster a pre-eminent model for animal genetics, and will further catalyse fundamental research on mechanisms of development, cell biology, genetics, disease, neurobiology, behaviour, physiology and evolution. Despite remarkable similarities among these Drosophila species, we identified many putatively non-neutral changes in protein-coding genes, non-coding RNA genes, and cis-regulatory regions. These may prove to underlie differences in the ecology and behaviour of these diverse species.

Exalign: a new method for comparative analysis of exon-intron gene structures.

The evolution of genes is usually studied and reconstructed at the sequence level, that is, by comparing and aligning their genomic, transcript or protein sequences. However, including the exon-intron structure of genes in the analysis can provide further and useful information, for example to draw reliable phylogenetic relationships left unsolved by traditional sequence-based evolutionary studies, or to shed further light on patterns of intron gain and loss. In spite of this, no tool especially devised for this task is currently available. In this work we present Exalign, an algorithm designed to retrieve, compare and search for the exon-intron structure of existing gene annotations, that has been implemented in a software tool freely accessible through a web interface as well as available for download. We present different applications of our method, from the reconstruction of the evolutionary history of homologous gene families to the detection of as of today unknown cases of intron loss in human and rodents, and, remarkably, two never reported intron gain events in human and mouse. The web interface for accessing Exalign is available at http://www.pesolelab.it/exalign/ or http://www.beacon.unimi.it/exalign/

Energy biogenesis: one key for coordinating two genomes.

In metazoan organisms, energy production is the only example of a process that is under dual genetic control: nuclear and mitochondrial. We used a genomic approach to examine how energy genes of both the nuclear and mitochondrial genomes are coordinated, and discovered a novel genetic regulatory circuit in Drosophila melanogaster that is surprisingly simple and parsimonious. This circuit is based on a single DNA regulatory element and can explain both intra- and inter-genomic coordinated expression of genes involved in energy production, including the full complement of mitochondrial and nuclear oxidative phosphorylation genes, and the genes involved in the Krebs cycle.

Characterization of Drosophila ATPsynC mutants as a new model of mitochondrial ATP synthase disorders.

Mitochondrial disorders associated with genetic defects of the ATP synthase are among the most deleterious diseases of the neuromuscular system that primarily manifest in newborns. Nevertheless, the number of established animal models for the elucidation of the molecular mechanisms behind such pathologies is limited. In this paper, we target the Drosophila melanogaster gene encoding for the ATP synthase subunit c, ATPsynC, in order to create a fruit fly model for investigating defects in mitochondrial bioenergetics and to better understand the comprehensive pathological spectrum associated with mitochondrial ATP synthase dysfunctions. Using P-element and EMS mutagenesis, we isolated a set of mutations showing a wide range of effects, from larval lethality to complex pleiotropic phenotypes encompassing developmental delay, early adult lethality, hypoactivity, sterility, hypofertility, aberrant male courtship behavior, locomotor defects and aberrant gonadogenesis. ATPsynC mutations impair ATP synthesis and mitochondrial morphology, and represent a powerful toolkit for the screening of genetic modifiers that can lead to potential therapeutic solutions. Furthermore, the molecular characterization of ATPsynC mutations allowed us to better understand the genetics of the ATPsynC locus and to define three broad pathological consequences of mutations affecting the mitochondrial ATP synthase functionality in Drosophila: i) pre-adult lethality; ii) multi-trait pathology accompanied by early adult lethality; iii) multi-trait adult pathology. We finally predict plausible parallelisms with genetic defects of mitochondrial ATP synthase in humans.

The nuclear OXPHOS genes in insecta: a common evolutionary origin, a common cis-regulatory motif, a common destiny for gene duplicates.

BACKGROUND: When orthologous sequences from species distributed throughout an optimal range of divergence times are available, comparative genomics is a powerful tool to address problems such as the identification of the forces that shape gene structure during evolution, although the functional constraints involved may vary in different genes and lineages. RESULTS: We identified and annotated in the MitoComp2 dataset the orthologs of 68 nuclear genes controlling oxidative phosphorylation in 11 Drosophilidae species and in five non-Drosophilidae insects, and compared them with each other and with their counterparts in three vertebrates (Fugu rubripes, Danio rerio and Homo sapiens) and in the cnidarian Nematostella vectensis, taking into account conservation of gene structure and regulatory motifs, and preservation of gene paralogs in the genome. Comparative analysis indicates that the ancestral insect OXPHOS genes were intron rich and that extensive intron loss and lineage-specific intron gain occurred during evolution. Comparison with vertebrates and cnidarians also shows that many OXPHOS gene introns predate the cnidarian/Bilateria evolutionary split. The nuclear respiratory gene element (NRG) has played a key role in the evolution of the insect OXPHOS genes; it is constantly conserved in the OXPHOS orthologs of all the insect species examined, while their duplicates either completely lack the element or possess only relics of the motif. CONCLUSION: Our observations reinforce the notion that the common ancestor of most animal phyla had intron-rich gene, and suggest that changes in the pattern of expression of the gene facilitate the fixation of duplications in the genome and the development of novel genetic functions.

The MitoDrome database annotates and compares the OXPHOS nuclear genes of Drosophila melanogaster, Drosophila pseudoobscura and Anopheles gambiae.

The oxidative phosphorylation (OXPHOS) is the primary energy-producing process of all aerobic organisms and the only cellular function under the dual control of both the mitochondrial and the nuclear genomes. Functional characterization and evolutionary study of the OXPHOS system is of great importance for the understanding of many as yet unclear aspects of nucleus-mitochondrion genomic co-evolution and co-regulation gene networks. The MitoDrome database is a web-based database which provides genomic annotations about nuclear genes of Drosophila melanogaster encoding for mitochondrial proteins. Recently, MitoDrome has included a new section annotating genomic information about OXPHOS genes in Drosophila pseudoobscura and Anopheles gambiae and their comparative analysis with their Drosophila melanogaster and human counterparts. The introduction of this new comparative annotation section into MitoDrome is expected to be a useful resource for both functional and structural genomics related to the OXPHOS system.

MitoDrome: a database of Drosophila melanogaster nuclear genes encoding proteins targeted to the mitochondrion.

Mitochondria are organelles present in the cytoplasm of most eukaryotic cells; although they have their own DNA, the majority of the proteins necessary for a functional mitochondrion are coded by the nuclear DNA and only after transcription and translation they are imported in the mitochondrion as proteins. The primary role of the mitochondrion is electron transport and oxidative phosphorylation. Although it has been studied for a long time, the interest of researchers in mitochondria is still alive thanks to the discovery of mitochondrial role in apoptosis, aging and cancer. Aim of the MitoDrome database is to annotate the Drosophila melanogaster nuclear genes coding for mitochondrial proteins in order to contribute to the functional characterization of nuclear genes coding for mitochondrial proteins and to knowledge of gene diseases related to mitochondrial dysfunctions. Indeed D. melanogaster is one of the most studied organisms and a model for the Human genome. Data are derived from the comparison of Human mitochondrial proteins versus the Drosophila genome, ESTs and cDNA sequence data available in the FlyBase database. Links from the MitoDrome entries to the related homologous entries available in MitoNuC will be soon imple-mented. The MitoDrome database is available at http://bighost.area.ba.cnr.it/BIG/MitoDrome. Data are organised in a flat-file format and can be retrieved using the SRS system.

Evolution of nuclearly encoded mitochondrial genes in Metazoa.

All Metazoan nuclear genomes underwent a continuous process of both complete and partial genetic material gain and loss. The forces modulating these events are also subject to the strict interaction between nuclear and mitochondrial (mt) genome. In this context we investigate the evolution of nuclear genes encoding proteins which target the mitochondrion, with a particular attention to genes involved in oxidative phosphorylation (OXPHOS), one of the most ancient and conserved functions. To examine thoroughly the evolutionary strategies that preserve OXPHOS and coordinate the two cellular genomes, a comparative analysis has been carried out for 78 OXPHOS gene families in several Metazoa (insects, tunicates, fishes and mammals). We demonstrate that the duplication rate of OXPHOS genes increases passing from invertebrates to vertebrates consistently with the total increase in genome size, but all species are prone to negatively select OXPHOS duplicates compared to the general trend of nuclear gene families. These results are consistent with the 'balance hypothesis' and, at least in insects, the expression of duplicate genes is low and strongly testis-biased.

A survey of the DNA sequences surrounding the Bari1 repeats in the pericentromeric h39 region of Drosophila melanogaster.

In Drosophila melanogaster, clustered copies of the Bari1 transposon are only present in the pericentromeric h39 region of the second chromosome, where other clusters of repetitive elements, either found organized in large tandem arrays only in the h39 region (Responder, PortoI), or both in the h39 region and in other heterochromatic regions (Hoppel), are also observed. The topological relationship among the repetitive sequences of the h39 region and the nature of the sequences separating its large repeat clusters are at present largely unknown. To get new insights on the sequence composition of the heterochromatin and on the forces governing its origin and maintenance, we have cloned and analyzed part of the DNA sequences flanking the h39 Bari1 repeats. In a region spanning 3 and 9 kb, respectively, from the ends of a Bari1 array we found only single copies of the PortoI and Hoppel transposable elements, and five copies of a variant form of the Responder repeats. No large tandem arrays of any repeated element were present. In addition, a highly conserved 596 bp sequence, that may have a functional role, is present on both sides of the Bari1 repeats. We suggest that the current organization of the h39 heterochromatin implies some topological or functional constraint that prevents the formation of further arrays of repetitive elements in the region.

A comparative study of the porin genes encoding VDAC, a voltage-dependent anion channel protein, in Anopheles gambiae and Drosophila melanogaster.

The protein called voltage-dependent anion-selective channel (VDAC), or mitochondrial porin, forms channels that provide the major pathway for small metabolites across the mitochondrial outer membrane. We have identified and sequenced agporin, a gene of the malaria vector mosquito Anopheles gambiae that conceptually encodes a protein with 73% identity to the VDAC protein encoded by the porin gene in Drosophila melanogaster. By in situ hybridization, we have localized agporin at region 35D on the right arm of A. gambiae chromosome 3, which is homologous to the 2L chromosomal arm of D. melanogaster where the porin gene resides. The comparison of agporin with its putative Drosophila counterpart revealed that both the nucleotide sequence and the structural organization of the two genes are strikingly conserved even though the ancestral lines of A. gambiae and D. melanogaster are thought to have diverged about 250 million years ago. Our results suggest that, while in yeast, plants, and mammals, VDAC isoforms are encoded by small multigene families and are able to compensate for each other at least partially, in A. gambiae a single gene encodes the VDAC protein.

Genetic, functional and evolutionary characterization of scox, the Drosophila melanogaster ortholog of the human SCO1 gene.

SCO proteins are copper-donor chaperones involved in the assembly of mitochondrial cytochrome c oxidase (COX). Mutations in the two human SCO-encoding genes, SCO1 and SCO2, produce tissue-specific COX deficiencies associated with distinct clinical phenotypes. Here, we report the identification and characterization of scox, the single Drosophila melanogaster SCO-encoding gene. Null mutations of the scox gene are associated with larval lethality, while mutations in its 5'UTR are associated with motor dysfunction and female sterile phenotypes. All mutant phenotypes may be rescued by a transgene encompassing wild-type scox. The analysis of the phenotypes associated with the D. melanogaster scox mutations shows that unimpaired COX assembly and activity is required for biological processes that specifically depend on an adequate energy supply. Finally, we identified the SCO1 orthologs in 39 eukaryotic species informative for a tentative reconstruction of the evolutionary history of the SCO function. Comparison of the exon/intron structure and other key features suggest that eukaryotic SCO genes descend from an intron-rich ancestral gene already present in the last common ancestor of lineages that diverged as early as metazoans and flowering plants.

Comparison of the oxidative phosphorylation (OXPHOS) nuclear genes in the genomes of Drosophila melanogaster, Drosophila pseudoobscura and Anopheles gambiae.

BACKGROUND: In eukaryotic cells, oxidative phosphorylation (OXPHOS) uses the products of both nuclear and mitochondrial genes to generate cellular ATP. Interspecies comparative analysis of these genes, which appear to be under strong functional constraints, may shed light on the evolutionary mechanisms that act on a set of genes correlated by function and subcellular localization of their products. RESULTS: We have identified and annotated the Drosophila melanogaster, D. pseudoobscura and Anopheles gambiae orthologs of 78 nuclear genes encoding mitochondrial proteins involved in oxidative phosphorylation by a comparative analysis of their genomic sequences and organization. We have also identified 47 genes in these three dipteran species each of which shares significant sequence homology with one of the above-mentioned OXPHOS orthologs, and which are likely to have originated by duplication during evolution. Gene structure and intron length are essentially conserved in the three species, although gain or loss of introns is common in A. gambiae. In most tissues of D. melanogaster and A. gambiae the expression level of the duplicate gene is much lower than that of the original gene, and in D. melanogaster at least, its expression is almost always strongly testis-biased, in contrast to the soma-biased expression of the parent gene. CONCLUSIONS: Quickly achieving an expression pattern different from the parent genes may be required for new OXPHOS gene duplicates to be maintained in the genome. This may be a general evolutionary mechanism for originating phenotypic changes that could lead to species differentiation.

Molecular clock and gene function.

Molecular phylogenies based on the molecular clock require the comparison of orthologous genes. Orthologous and paralogous genes usually have very different evolutionary fates. In general, orthologs keep the same functions in species, whereas, particularly over a long time span, paralogs diverge functionally and may become pseudogenes or get lost. In eukaryotic genomes, because of the degree of redundancy of genetic information, homologous genes are grouped in gene families, the evolution of which may differ greatly between the various organisms. This implies that each gene in a species does not always have an ortholog in another species and thus, due to multiple duplication events following a speciation, many orthologous clades of paralogs are generated. We are often dealing with a one-to-many or many-to-many relationship between genes. In this paper, we analyze the evolution of two gene families, the p53 gene family and the porin gene family. The evolution of the p53 family shows a one-to-many gene relationship going from invertebrates to vertebrates. In invertebrates only a single gene has been found, while in vertebrates three members of the family, namely p53, p63, and p73, are present. The evolution of porin (VDAC) genes (VDAC1, VDAC2, and VDAC3) is an example of a many-to-many gene relationship going from yeast to mammals. However, the porin gene redundancy found in invertebrates and possibly in some fishes may indicate a tendency to duplicate the genetic material, rather than a real need for function innovation.

Organization and possible origin of the Bari-1 cluster in the heterochromatic h39 region of Drosophila melanogaster.

The molecular organization of the heterochromatic h39 region of the Drosophila melanogaster second chromosome has been investigated by studying two BAC clones identified both by Southern blotting and by FISH experiments as containing tandem arrays of Bari1, a transposable element present only in this region. Such BAC clones appear to contain different portions of the h39 region since they differ in the DNA sequences flanking the Bari1 repeats on both sides. Thus, the 80 Bari1 copies estimated to be present in the h39 region are split into at least two separated subregions. On the basis of the analysis of the flanking sequences a possible mechanism depending on an aberrant activity of the Bari1 transposase is proposed for the genesis of the heterochromatic tandem arrays of the element.