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

Differential gene expression patterns between the head and thorax of Gynaephora aureata are associated with high-altitude adaptation

Grassland caterpillars (Lepidoptera: Erebidae: Gynaephora) are important pests in alpine meadows of the Qinghai-Tibetan Plateau (QTP). These pests have morphological, behavioral, and genetic adaptations for survival in high-altitude environments. However, mechanisms underlying high-altitude adaptation in QTP Gynaephora species remain largely unknown. Here, we performed a comparative analysis of the head and thorax transcriptomes of G. aureata to explore the genetic basis of high-altitude adaptation. We detected 8,736 significantly differentially expressed genes (sDEGs) between the head and thorax, including genes related to carbohydrate metabolism, lipid metabolism, epidermal proteins, and detoxification. These sDEGs were significantly enriched in 312 Gene Ontology terms and 16 KEGG pathways. We identified 73 pigment-associated genes, including 8 rhodopsin-associated genes, 19 ommochrome-associated genes, 1 pteridine-associated gene, 37 melanin-associated genes, and 12 heme-associated genes. These pigment-associated genes were related to the formation of the red head and black thorax of G. aureata. A key gene, yellow-h, in the melanin pathway was significantly upregulated in the thorax, suggesting that it is related to the formation of the black body and contributed to the adaptation of G. aureata to low temperatures and high ultraviolet radiation in the QTP. Another key gene, cardinal, in the ommochrome pathway was significantly upregulated in the head and may be related to red warning color formation. We also identified 107 olfactory-related genes in G. aureata, including genes encoding 29 odorant-binding proteins, 16 chemosensory proteins, 22 odorant receptor proteins, 14 ionotropic receptors, 12 gustatory receptors, 12 odorant degrading enzymes, and 2 sensory neuron membrane proteins. Diversification of olfactory-related genes may be associated with the feeding habits of G. aureata, including larvae dispersal and searching for plant resources available in the QTP. These results provide new insights into high-altitude adaptation of Gynaephora in the QTP and may contribute to the development of new control strategies for these pests.

Temperature and livestock grazing trigger transcriptome responses in bumblebees along an elevational gradient

Climate and land-use changes cause increasing stress to pollinators but the molecular pathways underlying stress responses are poorly understood. Here, we analyzed the transcriptomic response of Bombus lucorum workers to temperature and livestock grazing. Bumblebees sampled along an elevational gradient, and from differently managed grassland sites (livestock grazing vs unmanaged) in the German Alps did not differ in the expression of genes known for thermal stress responses. Instead, metabolic energy production pathways were upregulated in bumblebees sampled in mid- or high elevations or during cool temperatures. Extensive grazing pressure led to an upregulation of genetic pathways involved in immunoregulation and DNA-repair. We conclude that widespread bumblebees are tolerant toward temperature fluctuations in temperate mountain environments. Moderate temperature increases may even release bumblebees from metabolic stress. However, transcriptome responses to even moderate management regimes highlight the completely underestimated complexity of human influence on natural pollinators.

•

Upregulation of energy metabolism pathways in Bombus lucorum with increasing elevation

•

Genes known for thermal stress responses did not change with increased elevation

•

Bombus lucorum are tolerant toward relatively broad temperature fluctuations

•

Grazing lead to an upregulation in genetic information processes in B. lucorum

Evolution of Oxidative Phosphorylation (OXPHOS) Genes Reflecting the Evolutionary and Life Histories of Fig Wasps (Hymenoptera, Chalcidoidea)

Fig wasps are a peculiar group of insects which, for millions of years, have inhabited the enclosed syconia of fig trees. Considering the relatively closed and dark environment of fig syconia, we hypothesize that the fig wasps’ oxidative phosphorylation (OXPHOS) pathway, which is the main oxygen consumption and adenosine triphosphate (ATP) production system, may have adaptively evolved. In this study, we manually annotated the OXPHOS genes of 11 species of fig wasps, and compared the evolutionary patterns of OXPHOS genes for six pollinators and five non-pollinators. Thirteen mitochondrial protein-coding genes and 30 nuclear-coding single-copy orthologous genes were used to analyze the amino acid substitution rate and natural selection. The results showed high amino acid substitution rates of both mitochondrial and nuclear OXPHOS genes in fig wasps, implying the co-evolution of mitochondrial and nuclear genes. Our results further revealed that the OXPHOS-related genes evolved significantly faster in pollinators than in non-pollinators, and five genes had significant positive selection signals in the pollinator lineage, indicating that OXPHOS genes play an important role in the adaptation of pollinators. This study can help us understand the relationship between gene evolution and environmental adaptation.

Comparative transcriptome analysis of three invasive leafminer flies provides insights into interspecific competition

Liriomyza spp. (Diptera: Agromyzidae) represent a group of economically-significant highly polyphagous pests of plants grown in field and greenhouse conditions. Liriomyza spp. share similar biological and morphological characteristics, and complex interspecific interactions have been documented among these species in various geographical regions. Where the displacement of one of these species by the other has been studied, no unique mechanisms have been identified as causing it. The impact of competitive factors (such as, insecticide tolerance, thermotolerance, and adaptability to cropping systems) may be unique to specific geographic regions of Liriomyza spp., but more research is needed to confirm these hypotheses. In this study, RNA-seq was used to determine the transcriptomes of three closely-related leafminers, e.g. L. sativae, L. trifolii, and L. huidobrensis. Over 20 Gb of clean reads were generated and assembled into unique transcriptomes, and 38,747 unigenes were annotated in different databases. In pairwise comparisons, L. trifolii and L. sativae had more up-regulated genes than L. huidobrensis. With respect to common differentially-expressed genes (Co-DEGs), the three leafminers exhibited distinct groups of highly-expressed gene clusters. When genes related to competitive factors were compared, expression patterns in L. trifolii and L. sativae were more closely related to each other than to L. huidobrensis. The data suggest that DEGs involved in competitive factors may play a key role in competition and displacement of leafminers. The divergent genes identified in this study will be valuable in revealing possible mechanisms of invasion, displacement and interspecific competition in Liriomyza spp.

High-resolution transcriptional profiling of Anopheles gambiae spermatogenesis reveals mechanisms of sex chromosome regulation

Although of high priority for the development of genetic tools to control malaria-transmitting mosquitoes, only a few germline-specific regulatory regions have been characterised to date and the presence of global regulatory mechanisms, such as dosage compensation and meiotic sex chromosome inactivation (MSCI), are mostly assumed from transcriptomic analyses of reproductive tissues or whole gonads. In such studies, samples include a significant portion of somatic tissues inevitably complicating the reconstruction of a defined transcriptional map of gametogenesis. By exploiting recent advances in transgenic technologies and gene editing tools, combined with fluorescence-activated cell sorting and RNA sequencing, we have separated four distinct cell lineages from the Anopheles gambiae male gonads: premeiotic, meiotic (primary and secondary spermatocytes) and postmeiotic. By comparing the overall expression levels of X-linked and autosomal genes across the four populations, we revealed a striking transcriptional repression of the X chromosome coincident with the meiotic phase, classifiable as MSCI, and highlighted genes that may evade silencing. In addition, chromosome-wide median expression ratios of the premeiotic population confirmed the absence of dosage compensation in the male germline. Applying differential expression analysis, we highlighted genes and transcript isoforms enriched at specific timepoints and reconstructed the expression dynamics of the main biological processes regulating the key stages of sperm development and maturation. We generated the first transcriptomic atlas of A. gambiae spermatogenesis that will expand the available toolbox for the genetic engineering of vector control technologies. We also describe an innovative and multidimensional approach to isolate specific cell lineages that can be used for the targeted analysis of other A. gambiae organs or transferred to other medically relevant species and model organisms.

Using Mouse and Drosophila Models to Investigate the Mechanistic Links between Diet, Obesity, Type II Diabetes, and Cancer

Many of the links between diet and cancer are controversial and over simplified. To date, human epidemiological studies consistently reveal that patients who suffer diet-related obesity and/or type II diabetes have an increased risk of cancer, suffer more aggressive cancers, and respond poorly to current therapies. However, the underlying molecular mechanisms that increase cancer risk and decrease the response to cancer therapies in these patients remain largely unknown. Here, we review studies in mouse cancer models in which either dietary or genetic manipulation has been used to model obesity and/or type II diabetes. These studies demonstrate an emerging role for the conserved insulin and insulin-like growth factor signaling pathways as links between diet and cancer progression. However, these models are time consuming to develop and expensive to maintain. As the world faces an epidemic of obesity and type II diabetes we argue that the development of novel animal models is urgently required. We make the case for Drosophila as providing an unparalleled opportunity to combine dietary manipulation with models of human metabolic disease and cancer. Thus, combining diet and cancer models in Drosophila can rapidly and significantly advance our understanding of the conserved molecular mechanisms that link diet and diet-related metabolic disorders to poor cancer patient prognosis.

Variation in Mitochondria-Derived Transcript Levels Associated With DDT Resistance in the 91-R Strain of Drosophila melanogaster (Diptera: Drosophilidae)

The organochloride insecticide dichlorodiphenyltrichloroethane (DDT) and its metabolites can increase cellular levels of reactive oxygen species (ROS), cause mitochondrial dysfunction, and induce apoptosis. The highly DDT-resistant Drosophila melanogaster Meigen 1830 (Drosophila) strain, 91-R, and its susceptible control, 91-C, were used to investigate functional and structural changes among mitochondrial-derived pathways. Resequencing of mitochondrial genomes (mitogenomes) detected no structural differences between 91-R and 91-C, whereas RNA-seq suggested the differential expression of 221 mitochondrial-associated genes. Reverse transcriptase-quantitative PCR validation of 33 candidates confirmed that transcripts for six genes (Cyp12d1-p, Cyp12a4, cyt-c-d, COX5BL, COX7AL, CG17140) were significantly upregulated and two genes (Dif, Rel) were significantly downregulated in 91-R. Among the upregulated genes, four genes are duplicated within the reference genome (cyt-c-d, CG17140, COX5BL, and COX7AL). The predicted functions of the differentially expressed genes, or known functions of closely related genes, suggest that 91-R utilizes existing ROS regulation pathways of the mitochondria to combat increased ROS levels from exposure to DDT. This study represents, to our knowledge, the initial investigation of mitochondrial genome sequence variants and functional adaptations in responses to intense DDT selection and provides insights into potential adaptations of ROS management associated with DDT selection in Drosophila.

Few Nuclear-Encoded Mitochondrial Gene Duplicates Contribute to Male Germline-Specific Functions in Humans

Most of the genes encoding proteins that function in the mitochondria are located in the nucleus and are called nuclear-encoded mitochondrial genes, or N-mt genes. In Drosophila melanogaster, about 23% of N-mt genes fall into gene families, and all duplicates with tissue-biased expression (76%) are testis biased. These genes are enriched for energy-related functions and tend to be older than other duplicated genes in the genome. These patterns reveal strong selection for the retention of new genes for male germline mitochondrial functions. The two main forces that are likely to drive changes in mitochondrial functions are maternal inheritance of mitochondria and male–male competition for fertilization. Both are common among animals, suggesting similar N-mt gene duplication patterns in different species. To test this, we analyzed N-mt genes in the human genome. We find that about 18% of human N-mt genes fall into gene families, but unlike in Drosophila, only 28% of the N-mt duplicates have tissue-biased expression and only 36% of these have testis-biased expression. In addition, human testis-biased duplicated genes are younger than other duplicated genes in the genome and have diverse functions. These contrasting patterns between species might reflect either differences in selective pressures for germline energy-related or other mitochondrial functions during spermatogenesis and fertilization, or differences in the response to similar pressures.

Analysis of Drosophila melanogaster testis transcriptome

Background

The formation of matured and individual sperm involves a series of molecular and spectacular morphological changes of the developing cysts in Drosophila melanogaster testis. Recent advances in RNA Sequencing (RNA-Seq) technology help us to understand the complexity of eukaryotic transcriptomes by dissecting different tissues and developmental stages of organisms. To gain a better understanding of cellular differentiation of spermatogenesis, we applied RNA-Seq to analyse the testis-specific transcriptome, including coding and non-coding genes.

Results

We isolated three different parts of the wild-type testis by dissecting and cutting the different regions: 1.) the apical region, which contains stem cells and developing spermatocytes 2.) the middle region, with enrichment of meiotic cysts 3.) the basal region, which contains elongated post-meiotic cysts with spermatids. Total RNA was isolated from each region and analysed by next-generation sequencing. We collected data from the annotated 17412 Drosophila genes and identified 5381 genes with significant transcript accumulation differences between the regions, representing the main stages of spermatogenesis. We demonstrated for the first time the presence and region specific distribution of 2061 lncRNAs in testis, with 203 significant differences. Using the available modENCODE RNA-Seq data, we determined the tissue specificity indices of Drosophila genes. Combining the indices with our results, we identified genes with region-specific enrichment in testis.

Conclusion

By multiple analyses of our results and integrating existing knowledge about Drosophila melanogaster spermatogenesis to our dataset, we were able to describe transcript composition of different regions of Drosophila testis, including several stage-specific transcripts. We present searchable visualizations that can facilitate the identification of new components that play role in the organisation and composition of different stages of spermatogenesis, including the less known, but complex regulation of post-meiotic stages.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5085-z) contains supplementary material, which is available to authorized users.

Positive and relaxed selection associated with flight evolution and loss in insect transcriptomes

The evolution of powered flight is a major innovation that has facilitated the success of insects. Previously, studies of birds, bats, and insects have detected molecular signatures of differing selection regimes in energy-related genes associated with flight evolution and/or loss. Here, using DNA sequences from more than 1000 nuclear and mitochondrial protein-coding genes obtained from insect transcriptomes, we conduct a broader exploration of which gene categories display positive and relaxed selection at the origin of flight as well as with multiple independent losses of flight. We detected a number of categories of nuclear genes more often under positive selection in the lineage leading to the winged insects (Pterygota), related to catabolic processes such as proteases, as well as splicing-related genes. Flight loss was associated with relaxed selection signatures in splicing genes, mirroring the results for flight evolution. Similar to previous studies of flight loss in various animal taxa, we observed consistently higher nonsynonymous-to-synonymous substitution ratios in mitochondrial genes of flightless lineages, indicative of relaxed selection in energy-related genes. While oxidative phosphorylation genes were not detected as being under selection with the origin of flight specifically, they were most often detected as being under positive selection in holometabolous (complete metamorphosis) insects as compared with other insect lineages. This study supports some convergence in gene-specific selection pressures associated with flight ability, and the exploratory analysis provided some new insights into gene categories potentially associated with the gain and loss of flight in insects.

An in vivo genetic screen in Drosophila identifies the orthologue of human cancer/testis gene SPO11 among a network of targets to inhibit lethal(3)malignant brain tumour growth

Using transgenic RNAi technology, we have screened over 4000 genes to identify targets to inhibit malignant growth caused by the loss of function of lethal(3)malignant brain tumour in Drosophila in vivo. We have identified 131 targets, which belong to a wide range of gene ontologies. Most of these target genes are not significantly overexpressed in mbt tumours hence showing that, rather counterintuitively, tumour-linked overexpression is not a good predictor of functional requirement. Moreover, we have found that most of the genes upregulated in mbt tumours remain overexpressed in tumour-suppressed double-mutant conditions, hence revealing that most of the tumour transcriptome signature is not necessarily correlated with malignant growth. One of the identified target genes is meiotic W68 (mei-W68), the Drosophila orthologue of the human cancer/testis gene Sporulation-specific protein 11 (SPO11), the enzyme that catalyses the formation of meiotic double-strand breaks. We show that Drosophila mei-W68/SPO11 drives oncogenesis by causing DNA damage in a somatic tissue, hence providing the first instance in which a SPO11 orthologue is unequivocally shown to have a pro-tumoural role. Altogether, the results from this screen point to the possibility of investigating the function of human cancer relevant genes in a tractable experimental model organism like Drosophila.

The molecular evolutionary dynamics of oxidative phosphorylation (OXPHOS) genes in Hymenoptera

BACKGROUND

The primary energy-producing pathway in eukaryotic cells, the oxidative phosphorylation (OXPHOS) system, comprises proteins encoded by both mitochondrial and nuclear genes. To maintain the function of the OXPHOS system, the pattern of substitutions in mitochondrial and nuclear genes may not be completely independent. It has been suggested that slightly deleterious substitutions in mitochondrial genes are compensated by substitutions in the interacting nuclear genes due to positive selection. Among the four largest insect orders, Coleoptera (beetles), Hymenoptera (sawflies, wasps, ants, and bees), Diptera (midges, mosquitoes, and flies) and Lepidoptera (moths and butterflies), the mitochondrial genes of Hymenoptera exhibit an exceptionally high amino acid substitution rate while the evolution of nuclear OXPHOS genes is largely unknown. Therefore, Hymenoptera is an excellent model group for testing the hypothesis of positive selection driving the substitution rate of nuclear OXPHOS genes. In this study, we report the evolutionary rate of OXPHOS genes in Hymenoptera and test for evidence of positive selection in nuclear OXPHOS genes of Hymenoptera.

RESULTS

Our analyses revealed that the amino acid substitution rate of mitochondrial and nuclear OXPHOS genes in Hymenoptera is higher than that in other studied insect orders. In contrast, the amino acid substitution rate of non-OXPHOS genes in Hymenoptera is lower than the rate in other insect orders. Overall, we found the dN/dS ratio of the nuclear OXPHOS genes to be higher in Hymenoptera than in other insect orders. However, nuclear OXPHOS genes with high dN/dS ratio did not always exhibit a high amino acid substitution rate. Using branch-site and site model tests, we identified various codon sites that evolved under positive selection in nuclear OXPHOS genes.

CONCLUSIONS

Our results showed that nuclear OXPHOS genes in Hymenoptera are evolving faster than the genes in other three insect orders. The branch test suggested that while some nuclear OXPHOS genes in Hymenoptera show a signature of positive selection, the pattern is not consistent across all nuclear OXPHOS genes. As only few codon sites were under positive selection, we suggested that positive selection might not be the only factor contributing to the rapid evolution of nuclear OXPHOS genes in Hymenoptera.

Genome-Wide Analysis Suggests the Relaxed Purifying Selection Affect the Evolution of WOX Genes in Pyrus bretschneideri, Prunus persica, Prunus mume, and Fragaria vesca

WUSCHEL-related homeobox (WOX) family is one of the largest group of transcription factors (TFs) specifically found in plant kingdom. WOX TFs play an important role in plant development processes and evolutionary novelties. Although the roles of WOXs in Arabidopsis and rice have been well-studied, however, little are known about the relationships among the main clades in the molecular evolution of these genes in Rosaceae. Here, we carried out a genome-wide analysis and identified 14, 10, 10, and 9 of WOX genes from four Rosaceae species (Fragaria vesca, Prunus persica, Prunus mume, and Pyrus bretschneideri, respectively). According to evolutionary analysis, as well as amino acid sequences of their homodomains, these genes were divided into three clades with nine subgroups. Furthermore, due to the conserved structural patterns among these WOX genes, it was proposed that there should exist some highly conserved regions of microsynteny in the four Rosaceae species. Moreover, most of WOX gene pairs were presented with the conserved orientation among syntenic genome regions. In addition, according to substitution models analysis using PMAL software, no significant positive selection was detected, but type I functional divergence was identified among certain amino acids in WOX protein. These results revealed that the relaxed purifying selection might be the main driving force during the evolution of WOX genes in the tested Rosaceae species. Our result will be useful for further precise research on evolution of the WOX genes in family Rosaceae.

Mitochondrial-Nuclear Interactions Mediate Sex-Specific Transcriptional Profiles in Drosophila

The assembly and function of mitochondria require coordinated expression from two distinct genomes, the mitochondrial DNA (mtDNA) and nuclear DNA (nDNA). Mutations in either genome can be a source of phenotypic variation, yet their coexpression has been largely overlooked as a source of variation, particularly in the emerging paradigm of mitochondrial replacement therapy. Here we tested how the transcriptome responds to mtDNA and nDNA variation, along with mitonuclear interactions (mtDNA × nDNA) in Drosophila melanogaster. We used two mtDNA haplotypes that differ in a substantial number of single nucleotide polymorphisms, with >100 amino acid differences. We placed each haplotype on each of two D. melanogaster nuclear backgrounds and tested for transcription differences in both sexes. We found that large numbers of transcripts were differentially expressed between nuclear backgrounds, and that mtDNA type altered the expression of nDNA genes, suggesting a retrograde, trans effect of mitochondrial genotype. Females were generally more sensitive to genetic perturbation than males, and males demonstrated an asymmetrical effect of mtDNA in each nuclear background; mtDNA effects were nuclear-background specific. mtDNA-sensitive genes were not enriched in male- or female-limited expression space in either sex. Using a variety of differential expression analyses, we show the responses to mitonuclear covariation to be substantially different between the sexes, yet the mtDNA genes were consistently differentially expressed across nuclear backgrounds and sexes. Our results provide evidence that the main mtDNA effects can be consistent across nuclear backgrounds, but the interactions between mtDNA and nDNA can lead to sex-specific global transcript responses.

A mitochondrial DNA hypomorph of cytochrome oxidase specifically impairs male fertility in Drosophila melanogaster

Due to their strict maternal inheritance in most animals and plants, mitochondrial genomes are predicted to accumulate mutations that are beneficial or neutral in females but harmful in males. Although a few male-harming mtDNA mutations have been identified, consistent with this 'Mother's Curse', their effect on females has been largely unexplored. Here, we identify COII(G177S), a mtDNA hypomorph of cytochrome oxidase II, which specifically impairs male fertility due to defects in sperm development and function without impairing other male or female functions. COII(G177S) represents one of the clearest examples of a 'male-harming' mtDNA mutation in animals and suggest that the hypomorphic mtDNA mutations like COII(G177S) might specifically impair male gametogenesis. Intriguingly, some D. melanogaster nuclear genetic backgrounds can fully rescue COII(G177S) -associated sterility, consistent with previously proposed models that nuclear genomes can regulate the phenotypic manifestation of mtDNA mutations.

Gene flow between Drosophila yakuba and Drosophila santomea in subunit V of cytochrome c oxidase: A potential case of cytonuclear cointrogression

Introgression is the effective exchange of genetic information between species through natural hybridization. Previous genetic analyses of the Drosophila yakuba—D. santomea hybrid zone showed that the mitochondrial genome of D. yakuba had introgressed into D. santomea and completely replaced its native form. Since mitochondrial proteins work intimately with nuclear‐encoded proteins in the oxidative phosphorylation (OXPHOS) pathway, we hypothesized that some nuclear genes in OXPHOS cointrogressed along with the mitochondrial genome. We analyzed nucleotide variation in the 12 nuclear genes that form cytochrome c oxidase (COX) in 33 Drosophila lines. COX is an OXPHOS enzyme composed of both nuclear‐ and mitochondrial‐encoded proteins and shows evidence of cytonuclear coadaptation in some species. Using maximum‐likelihood methods, we detected significant gene flow from D. yakuba to D. santomea for the entire COX complex. Interestingly, the signal of introgression is concentrated in the three nuclear genes composing subunit V, which shows population migration rates significantly greater than the background level of introgression in these species. The detection of introgression in three proteins that work together, interact directly with the mitochondrial‐encoded core, and are critical for early COX assembly suggests this could be a case of cytonuclear cointrogression.

Mitochondrial diseases: Drosophila melanogaster as a model to evaluate potential therapeutics

While often presented as a single entity, mitochondrial diseases comprise a wide range of clinical, biochemical and genetic heterogeneous disorders. Among them, defects in the process of oxidative phosphorylation are the most prevalent. Despite intense research efforts, patients are still without effective treatment. An important part of the development of new therapeutics relies on predictive models of the pathology in order to assess their therapeutic potential. Since mitochondrial diseases are a heterogeneous group of progressive multisystemic disorders that can affect any organ at any time, the development of various in vivo models for the different diseases-associated genes defects will accelerate the search for effective therapeutics. Here, we review existing Drosophila melanogaster models for mitochondrial diseases, with a focus on alterations in oxidative phosphorylation, and discuss the potential of this powerful model organism in the process of drug target discovery. This article is part of a Directed Issue entitled: Energy Metabolism Disorders and Therapies.

Knockdown of ATPsyn-b caused larval growth defect and male infertility in Drosophila

The ATPsyn-b encoding for subunit b of ATP synthase in Drosophila melanogaster is proposed to act in ATP synthesis and phagocytosis, and has been identified as one of the sperm proteins in both Drosophila and mammals. At present, its details of functions in animal growth and spermatogenesis have not been reported. In this study, we knocked down ATPsyn-b using Drosophila lines expressing inducible hairpin RNAi constructs and Gal4 drivers. Ubiquitous knockdown of ATPsyn-b resulted in growth defects in larval stage as the larvae did not grow bigger than the size of normal second-instar larvae. Knockdown in testes did not interrupt the developmental excursion to viable adult flies, however, these male adults were sterile. Analyses of testes revealed disrupted nuclear bundles during spermatogenesis and abnormal shaping in spermatid elongation. There were no mature sperm in the seminal vesicle of ATPsyn-b knockdown male testes. These findings suggest us that ATPsyn-b acts in growth and male fertility of Drosophila.

Transcriptional control of an essential ribozyme in Drosophila reveals an ancient evolutionary divide in animals

Ribonuclease P (RNase P) is an essential enzyme required for 5'-maturation of tRNA. While an RNA-free, protein-based form of RNase P exists in eukaryotes, the ribonucleoprotein (RNP) form is found in all domains of life. The catalytic component of the RNP is an RNA known as RNase P RNA (RPR). Eukaryotic RPR genes are typically transcribed by RNA polymerase III (pol III). Here we showed that the RPR gene in Drosophila, which is annotated in the intron of a pol II-transcribed protein-coding gene, lacks signals for transcription by pol III. Using reporter gene constructs that include the RPR-coding intron from Drosophila, we found that the intron contains all the sequences necessary for production of mature RPR but is dependent on the promoter of the recipient gene for expression. We also demonstrated that the intron-coded RPR copurifies with RNase P and is required for its activity. Analysis of RPR genes in various animal genomes revealed a striking divide in the animal kingdom that separates insects and crustaceans into a single group in which RPR genes lack signals for independent transcription and are embedded in different protein-coding genes. Our findings provide evidence for a genetic event that occurred approximately 500 million years ago in the arthropod lineage, which switched the control of the transcription of RPR from pol III to pol II.

Inference of the oxidative stress network in Anopheles stephensi upon Plasmodium infection

Ookinete invasion of Anopheles midgut is a critical step for malaria transmission; the parasite numbers drop drastically and practically reach a minimum during the parasite's whole life cycle. At this stage, the parasite as well as the vector undergoes immense oxidative stress. Thereafter, the vector undergoes oxidative stress at different time points as the parasite invades its tissues during the parasite development. The present study was undertaken to reconstruct the network of differentially expressed genes involved in oxidative stress in Anopheles stephensi during Plasmodium development and maturation in the midgut. Using high throughput next generation sequencing methods, we generated the transcriptome of the An. stephensi midgut during Plasmodium vinckei petteri oocyst invasion of the midgut epithelium. Further, we utilized large datasets available on public domain on Anopheles during Plasmodium ookinete invasion and Drosophila datasets and arrived upon clusters of genes that may play a role in oxidative stress. Finally, we used support vector machines for the functional prediction of the un-annotated genes of An. stephensi. Integrating the results from all the different data analyses, we identified a total of 516 genes that were involved in oxidative stress in An. stephensi during Plasmodium development. The significantly regulated genes were further extracted from this gene cluster and used to infer an oxidative stress network of An. stephensi. Using system biology approaches, we have been able to ascertain the role of several putative genes in An. stephensi with respect to oxidative stress. Further experimental validations of these genes are underway.

Extensive translation of small Open Reading Frames revealed by Poly-Ribo-Seq

Thousands of small Open Reading Frames (smORFs) with the potential to encode small peptides of fewer than 100 amino acids exist in our genomes. However, the number of smORFs actually translated, and their molecular and functional roles are still unclear. In this study, we present a genome-wide assessment of smORF translation by ribosomal profiling of polysomal fractions in Drosophila. We detect two types of smORFs bound by multiple ribosomes and thus undergoing productive translation. The ‘longer’ smORFs of around 80 amino acids resemble canonical proteins in translational metrics and conservation, and display a propensity to contain transmembrane motifs. The ‘dwarf’ smORFs are in general shorter (around 20 amino-acid long), are mostly found in 5′-UTRs and non-coding RNAs, are less well conserved, and have no bioinformatic indicators of peptide function. Our findings indicate that thousands of smORFs are translated in metazoan genomes, reinforcing the idea that smORFs are an abundant and fundamental genome component.

DOI:http://dx.doi.org/10.7554/eLife.03528.001

To produce a protein, a stretch of DNA must first be transcribed to produce a molecule of messenger RNA (mRNA). The genetic information copied from the DNA is then read three letters at a time, in groups called codons. Each codon either encodes a particular amino acid to be added into a protein or provides further instructions: ‘start codons’ mark the beginning of a protein; ‘stop codons’ mark its end. The DNA between these two points is called an open reading frame (or ORF)—however, not all ORFs produce proteins.

Most proteins are made of several hundred amino acids, but the genomes of animals contain thousands of ORFs that would generate much smaller proteins made of fewer than 100 amino acids, if they were translated. It is, however, unclear how many of these small ORFs are converted into mRNA molecules and functional proteins.

Ribosomes are large molecular machines that translate the code in mRNA molecules and join together the appropriate amino acids in the right order to make a protein. Ribosome profiling is a technique that identifies which mRNA molecules are translated into proteins by determining the sequences of all the mRNA molecules bound to ribosomes at a particular moment. The mRNA sequences can then be compared with the sequence of the whole genome to work out which ORFs they correspond to. Ribosome profiling has been used to detect translated small ORFs, but the method yields a relatively high false positive rate as some mRNAs can bind to ribosomes without being translated.

To better detect small protein-producing ORFs, Aspden et al. developed a technique based on ribosome profiling called Poly-Ribo-Seq. The method takes advantage of the fact that during active translation, clusters of multiple ribosomes, called polysomes, bind mRNAs. Poly-Ribo-Seq isolates these polysomes and determines the sequence bound by each of the ribosomes, thereby reducing the number of false positives.

Applying Poly-Ribo-Seq to cells from the fruit fly Drosophila allowed Aspden et al. to identify two types of short ORF. The first type codes for proteins that are around 80 amino acids long and are translated with the same efficiency as larger ORFs. The sequences of these ORFs are found in other species, match at least in part sequences of known functional ORFs, and the proteins produced are found in specific locations inside cells. These small proteins may contribute to membrane integrity or function. Together, these properties suggest that these mRNAs create functional small proteins.

The second pool consists of very small ORFs (‘dwarf smORFs’) that code for around 20 amino acids, which are not translated so often and do not show many similarities with other species.

As the findings of Aspden et al. suggest that a large fraction of Drosophila small ORFs are translated into proteins, the next challenge will be to determine the roles of these small proteins in cells.

DOI:http://dx.doi.org/10.7554/eLife.03528.002

Mitonuclear interactions: evolutionary consequences over multiple biological scales

Fundamental biological processes hinge on coordinated interactions between genes spanning two obligate genomes--mitochondrial and nuclear. These interactions are key to complex life, and allelic variation that accumulates and persists at the loci embroiled in such intergenomic interactions should therefore be subjected to intense selection to maintain integrity of the mitochondrial electron transport system. Here, we compile evidence that suggests that mitochondrial-nuclear (mitonuclear) allelic interactions are evolutionarily significant modulators of the expression of key health-related and life-history phenotypes, across several biological scales--within species (intra- and interpopulational) and between species. We then introduce a new frontier for the study of mitonuclear interactions--those that occur within individuals, and are fuelled by the mtDNA heteroplasmy and the existence of nuclear-encoded mitochondrial gene duplicates and isoforms. Empirical evidence supports the idea of high-resolution tissue- and environment-specific modulation of intraindividual mitonuclear interactions. Predicting the penetrance, severity and expression patterns of mtDNA-induced mitochondrial diseases remains a conundrum. We contend that a deeper understanding of the dynamics and ramifications of mitonuclear interactions, across all biological levels, will provide key insights that tangibly advance our understanding, not only of core evolutionary processes, but also of the complex genetics underlying human mitochondrial disease.

De novo characterization of the gene-rich transcriptomes of two color-polymorphic spiders, Theridion grallator and T. californicum (Araneae: Theridiidae), with special reference to pigment genes

Background

A number of spider species within the family Theridiidae exhibit a dramatic abdominal (opisthosomal) color polymorphism. The polymorphism is inherited in a broadly Mendelian fashion and in some species consists of dozens of discrete morphs that are convergent across taxa and populations. Few genomic resources exist for spiders. Here, as a first necessary step towards identifying the genetic basis for this trait we present the near complete transcriptomes of two species: the Hawaiian happy-face spider Theridion grallator and Theridion californicum. We mined the gene complement for pigment-pathway genes and examined differential expression (DE) between morphs that are unpatterned (plain yellow) and patterned (yellow with superimposed patches of red, white or very dark brown).

Results

By deep sequencing both RNA-seq and normalized cDNA libraries from pooled specimens of each species we were able to assemble a comprehensive gene set for both species that we estimate to be 98-99% complete. It is likely that these species express more than 20,000 protein-coding genes, perhaps 4.5% (ca. 870) of which might be unique to spiders. Mining for pigment-associated Drosophila melanogaster genes indicated the presence of all ommochrome pathway genes and most pteridine pathway genes and DE analyses further indicate a possible role for the pteridine pathway in theridiid color patterning.

Conclusions

Based upon our estimates, T. grallator and T. californicum express a large inventory of protein-coding genes. Our comprehensive assembly illustrates the continuing value of sequencing normalized cDNA libraries in addition to RNA-seq in order to generate a reference transcriptome for non-model species. The identification of pteridine-related genes and their possible involvement in color patterning is a novel finding in spiders and one that suggests a biochemical link between guanine deposits and the pigments exhibited by these species.

An integrative approach to ortholog prediction for disease-focused and other functional studies

Background

Mapping of orthologous genes among species serves an important role in functional genomics by allowing researchers to develop hypotheses about gene function in one species based on what is known about the functions of orthologs in other species. Several tools for predicting orthologous gene relationships are available. However, these tools can give different results and identification of predicted orthologs is not always straightforward.

Results

We report a simple but effective tool, the Drosophila RNAi Screening Center Integrative Ortholog Prediction Tool (DIOPT; http://www.flyrnai.org/diopt), for rapid identification of orthologs. DIOPT integrates existing approaches, facilitating rapid identification of orthologs among human, mouse, zebrafish, C. elegans, Drosophila, and S. cerevisiae. As compared to individual tools, DIOPT shows increased sensitivity with only a modest decrease in specificity. Moreover, the flexibility built into the DIOPT graphical user interface allows researchers with different goals to appropriately 'cast a wide net' or limit results to highest confidence predictions. DIOPT also displays protein and domain alignments, including percent amino acid identity, for predicted ortholog pairs. This helps users identify the most appropriate matches among multiple possible orthologs. To facilitate using model organisms for functional analysis of human disease-associated genes, we used DIOPT to predict high-confidence orthologs of disease genes in Online Mendelian Inheritance in Man (OMIM) and genes in genome-wide association study (GWAS) data sets. The results are accessible through the DIOPT diseases and traits query tool (DIOPT-DIST; http://www.flyrnai.org/diopt-dist).

Conclusions

DIOPT and DIOPT-DIST are useful resources for researchers working with model organisms, especially those who are interested in exploiting model organisms such as Drosophila to study the functions of human disease genes.

Mitochondrial genome deletions and minicircles are common in lice (Insecta: Phthiraptera)

Background

The gene composition, gene order and structure of the mitochondrial genome are remarkably stable across bilaterian animals. Lice (Insecta: Phthiraptera) are a major exception to this genomic stability in that the canonical single chromosome with 37 genes found in almost all other bilaterians has been lost in multiple lineages in favour of multiple, minicircular chromosomes with less than 37 genes on each chromosome.

Results

Minicircular mt genomes are found in six of the ten louse species examined to date and three types of minicircles were identified: heteroplasmic minicircles which coexist with full sized mt genomes (type 1); multigene chromosomes with short, simple control regions, we infer that the genome consists of several such chromosomes (type 2); and multiple, single to three gene chromosomes with large, complex control regions (type 3). Mapping minicircle types onto a phylogenetic tree of lice fails to show a pattern of their occurrence consistent with an evolutionary series of minicircle types. Analysis of the nuclear-encoded, mitochondrially-targetted genes inferred from the body louse, Pediculus, suggests that the loss of mitochondrial single-stranded binding protein (mtSSB) may be responsible for the presence of minicircles in at least species with the most derived type 3 minicircles (Pediculus, Damalinia).

Conclusions

Minicircular mt genomes are common in lice and appear to have arisen multiple times within the group. Life history adaptive explanations which attribute minicircular mt genomes in lice to the adoption of blood-feeding in the Anoplura are not supported by this expanded data set as minicircles are found in multiple non-blood feeding louse groups but are not found in the blood-feeding genus Heterodoxus. In contrast, a mechanist explanation based on the loss of mtSSB suggests that minicircles may be selectively favoured due to the incapacity of the mt replisome to synthesize long replicative products without mtSSB and thus the loss of this gene lead to the formation of minicircles in lice.

Analyses of nuclearly encoded mitochondrial genes suggest gene duplication as a mechanism for resolving intralocus sexually antagonistic conflict in Drosophila

Gene duplication is probably the most important mechanism for generating new gene functions. However, gene duplication has been overlooked as a potentially effective way to resolve genetic conflicts. Here, we analyze the entire set of Drosophila melanogaster nuclearly encoded mitochondrial duplicate genes and show that both RNA- and DNA-mediated mitochondrial gene duplications exhibit an unexpectedly high rate of relocation (change in location between parental and duplicated gene) as well as an extreme tendency to avoid the X chromosome. These trends are likely related to our observation that relocated genes tend to have testis-specific expression. We also infer that these trends hold across the entire Drosophila genus. Importantly, analyses of gene ontology and functional interaction networks show that there is an overrepresentation of energy production-related functions in these mitochondrial duplicates. We discuss different hypotheses to explain our results and conclude that our findings substantiate the hypothesis that gene duplication for male germline function is likely a mechanism to resolve intralocus sexually antagonistic conflicts that we propose are common in testis. In the case of nuclearly encoded mitochondrial duplicates, our hypothesis is that past sexually antagonistic conflict related to mitochondrial energy function in Drosophila was resolved by gene duplication.

Drcd-1 related: a positively selected spermatogenesis retrogene in Drosophila

Gene duplication is a major force driving genome evolution, and examples of this mode of evolution and of the functions of duplicated genes are needed to reveal general patterns. Here, our study focuses on a particular retrogene (i.e., CG9573) that originated about 5-13 million years ago that we have named Drcd-1 related. It originated in Drosophila through retroposition of the parental gene Required for cell differentiation 1 of Drosophila (Drcd-1; CG14213), which is a known transcription cofactor. Drcd-1r is only present in D. melanogaster, D. simulans, D. sechellia, and D. mauritiana. Drcd-1r is an X to autosome retroposition event. Many retrogenes are X to autosome copies and it has been shown that positive selection underlies this bias. We sought to understand Drcd-1r mode of evolution and function to contribute to the understanding of the selective pressures acting on X to autosome retrogenes. Drcd-1r overlaps with another gene, it is within the 3' UTR of the gene CG13102 and is encoded in the opposite orientation. We have studied the characteristics of the transcripts and quantified expression of CG13102 and Drcd-1r in wild-type flies. We found that Drcd-1r is transcribed specifically in testes. We also studied the molecular evolution of Drcd-1r and Drcd-1 and found that the parental gene has evolved under very strong purifying selection but the retrogene has evolved very rapidly (Ka/Ks ~1) under both positive and purifying selection, as revealed using divergence and polymorphism data. These results indicate that Drcd-1r has a novel function in the Drosophila testes. To further explore Drcd-1r function we used a strain containing a P element inserted in the region where CG13102 and Drcd-1r are located that shows recessive male sterility. Analysis of this strain reveals the difficulties that can be encountered in studying the functions of genes with overlapping transcripts. Avenues for studying of the function of this gene are proposed.

The biochemical properties of the mitochondrial thiamine pyrophosphate carrier from Drosophila melanogaster

The mitochondrial carriers are a family of transport proteins that shuttle metabolites, nucleotides and cofactors across the inner mitochondrial membrane. The genome of Drosophila melanogaster encodes at least 46 members of this family. Only five of these have been characterized, whereas the transport functions of the remainder cannot be assessed with certainty. In the present study, we report the functional identification of two D. melanogaster genes distantly related to the human and yeast thiamine pyrophosphate carrier (TPC) genes as well as the corresponding expression pattern throughout development. Furthermore, the functional characterization of the D. melanogaster mitochondrial thiamine pyrophosphate carrier protein (DmTpc1p) is described. DmTpc1p was over-expressed in bacteria, the purified protein was reconstituted into liposomes, and its transport properties and kinetic parameters were characterized. Reconstituted DmTpc1p transports thiamine pyrophosphate and, to a lesser extent, pyrophosphate, ADP, ATP and other nucleotides. The expression of DmTpc1p in Saccharomyces cerevisiaeTPC1 null mutant abolishes the growth defect on fermentable carbon sources. The main role of DmTpc1p is to import thiamine pyrophosphate into mitochondria by exchange with intramitochondrial ATP and/or ADP.

The telomere-linked helicase (TLH) gene family in Magnaporthe oryzae: revised gene structure reveals a novel TLH-specific protein motif

Telomere-linked RecQ helicase (TLH) genes have been identified in several fungi, where they occur as small gene families with each member copy residing within ~10 kb of a telomere. Here we describe the characterization of all 11 TLH gene copies in the reference strain of the fungus Magnaporthe oryzae. A consensus gene prediction revealed that the previously reported TLH1 gene is actually a mutated copy, and the full-length gene is almost two times longer. Only four full-length TLH genes were present in the strain that was analyzed, with the remaining copies containing premature stops caused by point mutations, indels, transposon insertions, and a telomere truncation. Interestingly, all of the TLH gene copies possessed numerous mutations indicative of the action of the repeat-induced point mutation process. However, there was evidence of purifying selection indicating maintenance of gene function. Alignment of full-length proteins from M. oryzae, Schizosaccharomyces pombe and M. anisopliae revealed the presence of a novel, highly conserved protein motif which suggests that the telomere-linked helicases have different functions and/or substrates to the RecQ helicases encoded by "internal" genes.

Duplicated proteasome subunit genes in Drosophila and their roles in spermatogenesis

The proteasome is a large, multisubunit complex that acts as the cell's 'protein-degrading machine' in the ubiquitin-mediated proteolytic pathway for regulated protein turnover. Although proteasomes are usually thought of as being homogeneous structures, recent studies have revealed their more dynamic and heterogeneous nature. For example, in a number of plant and animal species, multiple isoforms of several proteasome subunits, encoded by paralogous genes, have been discovered, and in some cases, these alternative isoforms have been shown to be functionally distinct from their conventional counterparts. A particularly striking example of this phenomenon is seen in Drosophila melanogaster, where 12 of the 33 subunits that make up the 26S proteasome holoenzyme are represented in the genome by multiple paralogous genes. Remarkably, in every case, the 'extra' genes are expressed in a testis-specific manner. Here, we describe the extent and nature of these testis-specific gene duplications and discuss their functional significance, and speculate on why this situation might have evolved.

Nuclear respiratory factor 1 controls myocyte enhancer factor 2A transcription to provide a mechanism for coordinate expression of respiratory chain subunits

Nuclear respiratory factors NRF1 and NRF2 regulate the expression of nuclear genes encoding heme biosynthetic enzymes, proteins required for mitochondrial genome transcription and protein import, and numerous respiratory chain subunits. NRFs thereby coordinate the expression of nuclear and mitochondrial genes relevant to mitochondrial biogenesis and respiration. Only two of the nuclear-encoded respiratory chain subunits have evolutionarily conserved tissue-specific forms: the cytochrome c oxidase (COX) subunits VIa and VIIa heart/muscle (H) and ubiquitous (L) isoforms. We used genome comparisons to conclude that the promoter regions of COX6A(H) and COX7A(H) lack NRF sites but have conserved myocyte enhancer factor 2 (MEF2) elements. We show that MEF2A mRNA is induced with forced expression of NRF1 and that the MEF2A 5'-regulatory region contains an evolutionarily conserved canonical element that binds endogenous NRF1 in chromatin immunoprecipitation (ChIP) assays. NRF1 regulates MEF2A promoter-reporters according to overexpression, RNA interference underexpression, and promoter element mutation studies. As there are four mammalian MEF2 isotypes, we used an isoform-specific antibody in ChIP to confirm MEF2A binding to the COX6A(H) promoter. These findings support a role for MEF2A as an intermediary in coordinating respiratory chain subunit expression in heart and muscle through a NRF1 --> MEF2A --> COX(H) transcriptional cascade. MEF2A also bound the MEF2A and PPARGC1A promoters in ChIP, placing it within a feedback loop with PGC1alpha in controlling NRF1 activity. Interruption of this cascade and loop may account for striated muscle mitochondrial defects in mef2a null mice. Our findings also account for the previously described indirect regulation by NRF1 of other MEF2 targets in muscle such as GLUT4.

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.

Drosophila melanogaster as a model system to study mitochondrial biology

Mitochondria play an essential role in cellular homeostasis. Although in the last few decades our knowledge of mitochondria has increased substantially, the mechanisms involved in the control of mitochondrial biogenesis remain largely unknown. The powerful genetics of Drosophila combined with a wealth of available cell and molecular biology techniques, make this organism an excellent system to study mitochondria. In this chapter we will review briefly the opportunities that Drosophila offers as a model system and describe in detail how to purify mitochondria from Drosophila and to perform the analysis of developmental gene expression using in situ hybridization.

Metazoan OXPHOS gene families: Evolutionary forces at the level of mitochondrial and nuclear genomes

Mitochondrial and nuclear DNAs contribute to encode the whole mitochondrial protein complement. The two genomes possess highly divergent features and properties, but the forces influencing their evolution, even if different, require strong coordination. The gene content of mitochondrial genome in all Metazoa is in a frozen state with only few exceptions and thus mitochondrial genome plasticity especially concerns some molecular features, i.e. base composition, codon usage, evolutionary rates. In contrast the high plasticity of nuclear genomes is particularly evident at the macroscopic level, since its redundancy represents the main feature able to introduce genetic material for evolutionary innovations. In this context, genes involved in oxidative phosphorylation (OXPHOS) represent a classical example of the different evolutionary behaviour of mitochondrial and nuclear genomes. The simple DNA sequence of Cytochrome c oxidase I (encoded by the mitochondrial genome) seems to be able to distinguish intra- and inter-species relations between organisms (DNA Barcode). Some OXPHOS subunits (cytochrome c, subunit c of ATP synthase and MLRQ) are encoded by several nuclear duplicated genes which still represent the trace of an ancient segmental/genome duplication event at the origin of vertebrates.

Expression and localization in spermatozoa of the mitochondrial porin isoform 2 in Drosophila melanogaster

Mitochondrial porins or VDACs (voltage-dependent anion-selective channels) are transmembrane pore-forming proteins. In eukaryotic genomes multiple genes coding for VDAC homologues have been discovered, but their function remains unknown. In Drosophila melanogaster three additional genes homologous to the gene porin have been found. In a previous report we have expressed in vitro Porin 2 (gene GC17137) and we have found that the reconstituted protein shows pore-forming activity but it is cation-selective and poorly dependent from voltage. In this work we have characterized the expression pattern of Porin 2. Amplification upon germinal and somatic or stage specific mRNA showed that the highest transcription level of Porin 2 is in testis. Western blot analysis performed with antibodies raised against the recombinant Porin 2 confirmed a high level of expression in the fly spermatozoa. Immuno-histochemical studies indicate that Porin 2 is selectively present in spermatozoa tail, where the mitochondria are located, but not in spermatocytes. A lethal mutant of D. melanogaster carrying a P-element in the first intron of the porin (Porin 1) gene hinders the expression of both Porin 1 and 2. Our results suggest that Porin 2 is truly expressed and that it is required for functional germinal tissues.

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.

Evolution of ATP synthase subunit c and cytochrome c gene families in selected Metazoan classes

To investigate the integrated evolution of mitochondrial and nuclear genomes in the eukaryotic cell, we have focused our attention on OXPHOS (oxidative phosphorylation) gene families which encode proteins involved in the main mitochondrial function. The present study reports the phylogenetic analysis of two OXPHOS gene families: ATP synthase subunit c (or lipid binding protein, LBP) and Cytochrome c (Cytc). Both gene families possess a higher expansion trend than the typically low duplication rate of OXPHOS genes in Metazoa, but follow a completely different evolutionary history, especially in mammals. LBP is represented by three well conserved isoforms in all mammals (P1, P2, P3): only P3 possesses a clearly conserved isoform in all Vertebrates, P1 and P2 were already present before the bird-mammal divergence and there are preliminary evidence from the in silico analysis that P1, the most evolutionary divergent isoform, is poorly expressed and not regulated by NRF1. In contrast, Cytc family presents at least two duplicated genes in all the analysed Vertebrates, is subject to a high expansion trend, especially of processed pseudogenes in mammals, and some events of gain and loss of function can be supposed.

The role of mitochondrial respiration in physiological and evolutionary adaptation

Aerobic mitochondria serve as the power sources of eukaryotes by producing ATP through oxidative phosphorylation (OXPHOS). The enzymes involved in OXPHOS are multisubunit complexes encoded by both nuclear and mitochondrial DNA. Thus, regulation of respiration is necessarily a highly coordinated process that must organize production, assembly and function of mitochondria to meet an organism's energetic needs. Here I review the role of OXPHOS in metabolic adaptation and diversification of higher animals. On a physiological timescale, endocrine-initiated signaling pathways allow organisms to modulate respiratory enzyme concentration and function under changing environmental conditions. On an evolutionary timescale, mitochondrial enzymes are targets of natural selection, balancing cytonuclear coevolutionary constraints against physiological innovation. By synthesizing our knowledge of biochemistry, physiology and evolution of respiratory regulation, I propose that we can now explore questions at the interface of these fields, from molecular translation of environmental cues to selection on mitochondrial haplotype variation.

MiGenes: a searchable interspecies database of mitochondrial proteins curated using gene ontology annotation

MOTIVATION

There has been an explosion of interest in the role of mitochondria in programmed cell death and other fundamental pathological processes underlying the development of human diseases. Nevertheless, the inventory of mitochondrial proteins encoded in the nuclear genome remains incomplete, providing an impediment to mitochondrial research at the interface with systems biology. We created the MiGenes database to further define the scope of the mitochondrial proteome in humans and model organisms including mice, rats, flies and worms as well as budding and fission yeasts. MiGenes is intended to stimulate mitochondrial research using model organisms.

SUMMARY

MiGenes is a large-scale relational database that is automatically updated to keep pace with advances in mitochondrial proteomics and is curated to assure that the designation of proteins as mitochondrial reflects gene ontology (GO) annotations supported by high-quality evidence codes. A set of postulates is proposed to help define which proteins are authentic components of mitochondria. MiGenes incorporates >1160 new GO annotations to human, mouse and rat protein records, 370 of which represent the first GO annotation reflecting a mitochondrial localization. MiGenes employs a flexible search interface that permits batchwise accession number searches to support high-throughput proteomic studies. A web interface is provided to permit members of the mitochondrial research community to suggest modifications in protein annotations or mitochondrial status.

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.