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

Generation of Mutants from the 57B Region of Drosophila melanogaster

The 57B region of Drosophila melanogaster includes a cluster of the three homeobox genes orthopedia (otp), Drosophila Retinal homeobox (DRx), and homeobrain (hbn). In an attempt to isolate mutants for these genes, we performed an EMS mutagenesis and isolated lethal mutants from the 57B region, among them mutants for otp, DRx, and hbn. With the help of two newly generated deletions from the 57B region, we mapped additional mutants to specific chromosomal intervals and identified several of these mutants from the 57B region molecularly. In addition, we generated mutants for CG15651 and RIC-3 by gene targeting and mutants for the genes CG9344, CG15649, CG15650, and ND-B14.7 using the CRISPR/Cas9 system. We determined the lethality period during development for most isolated mutants. In total, we analysed alleles from nine different genes from the 57B region of Drosophila, which could now be used to further explore the functions of the corresponding genes in the future.

Structural and Comparative Analyses of Insects Suggest the Presence of an Ultra-Conserved Regulatory Element of the Genes Encoding Vacuolar-Type ATPase Subunits and Assembly Factors

Gene and genome comparison represent an invaluable tool to identify evolutionarily conserved sequences with possible functional significance. In this work, we have analyzed orthologous genes encoding subunits and assembly factors of the V-ATPase complex, an important enzymatic complex of the vacuolar and lysosomal compartments of the eukaryotic cell with storage and recycling functions, respectively, as well as the main pump in the plasma membrane that energizes the epithelial transport in insects. This study involves 70 insect species belonging to eight insect orders. We highlighted the conservation of a short sequence in the genes encoding subunits of the V-ATPase complex and their assembly factors analyzed with respect to their exon-intron organization of those genes. This study offers the possibility to study ultra-conserved regulatory elements under an evolutionary perspective, with the aim of expanding our knowledge on the regulation of complex gene networks at the basis of organellar biogenesis and cellular organization.

Microtubule Organizing Centers Contain Testis-Specific γ-TuRC Proteins in Spermatids of Drosophila

Microtubule nucleation in eukaryotes is primarily promoted by γ-tubulin and the evolutionary conserved protein complex, γ-Tubulin Ring Complex (γ-TuRC). γ-TuRC is part of the centrosome and basal body, which are the best-known microtubule-organizing centers. Centrosomes undergo intensive and dynamic changes during spermatogenesis, as they turn into basal bodies, a prerequisite for axoneme formation during spermatogenesis. Here we describe the existence of a novel, tissue-specific γ-TuRC in Drosophila. We characterize three genes encoding testis-specific components of γ-TuRC (t-γ-TuRC) and find that presence of t-γ-TuRC is essential to male fertility. We show the diverse subcellular distribution of the t-γ-TuRC proteins during post-meiotic development, at first at the centriole adjunct and then also on the anterior tip of the nucleus, and finally, they appear in the tail region, close to the mitochondria. We also prove the physical interactions between the t-γ-TuRC members, γ-tubulin and Mozart1. Our results further indicate heterogeneity in γ-TuRC composition during spermatogenesis and suggest that the different post-meiotic microtubule organizing centers are orchestrated by testis-specific gene products, including t-γ-TuRC.

The mitochondrial aspartate/glutamate carrier (AGC or Aralar1) isoforms in D. melanogaster: biochemical characterization, gene structure, and evolutionary analysis

BACKGROUND

In man two mitochondrial aspartate/glutamate carrier (AGC) isoforms, known as aralar and citrin, are required to accomplish several metabolic pathways. In order to fill the existing gap of knowledge in Drosophila melanogaster, we have studied aralar1 gene, orthologue of human AGC-encoding genes in this organism.

METHODS

The blastp algorithm and the "reciprocal best hit" approach have been used to identify the human orthologue of AGCs in Drosophilidae and non-Drosophilidae. Aralar1 proteins have been overexpressed in Escherichia coli and functionally reconstituted into liposomes for transport assays.

RESULTS

The transcriptional organization of aralar1 comprises six isoforms, three constitutively expressed (aralar1-RA, RD and RF), and the remaining three distributed during the development or in different tissues (aralar1-RB, RC and RE). Aralar1-PA and Aralar1-PE, representative of all isoforms, have been biochemically characterized. Recombinant Aralar1-PA and Aralar1-PE proteins share similar efficiency to exchange glutamate against aspartate, and same substrate affinities than the human isoforms. Interestingly, although Aralar1-PA and Aralar1-PE diverge only in their EF-hand 8, they greatly differ in their specific activities and substrate specificity.

CONCLUSIONS

The tight regulation of aralar1 transcripts expression and the high request of aspartate and glutamate during early embryogenesis suggest a crucial role of Aralar1 in this Drosophila developmental stage. Furthermore, biochemical characterization and calcium sensitivity have identified Aralar1-PA and Aralar1-PE as the human aralar and citrin counterparts, respectively.

GENERAL SIGNIFICANCE

The functional characterization of the fruit fly mitochondrial AGC transporter represents a crucial step toward a complete understanding of the metabolic events acting during early embryogenesis.

Genomic signatures of globally enhanced gene duplicate accumulation in the megadiverse higher Diptera fueling intralocus sexual conflict resolution

Gene duplication is an important source of evolutionary innovation. To explore the relative impact of gene duplication during the diversification of major insect model system lineages, we performed a comparative analysis of lineage-specific gene duplications in the fruit fly Drosophila melanogaster (Diptera: Brachycera), the mosquito Anopheles gambiae (Diptera: Culicomorpha), the red flour beetle Tribolium castaneum (Coleoptera), and the honeybee Apis mellifera (Hymenoptera). Focusing on close to 6,000 insect core gene families containing maximally six paralogs, we detected a conspicuously higher number of lineage-specific duplications in Drosophila (689) compared to Anopheles (315), Tribolium (386), and Apis (223). Based on analyses of sequence divergence, phylogenetic distribution, and gene ontology information, we present evidence that an increased background rate of gene duplicate accumulation played an exceptional role during the diversification of the higher Diptera (Brachycera), in part by providing enriched opportunities for intralocus sexual conflict resolution, which may have boosted speciation rates during the early radiation of the megadiverse brachyceran subclade Schizophora.

Drosophila melanogaster Mitochondrial Carriers: Similarities and Differences with the Human Carriers

Mitochondrial carriers are a family of structurally related proteins responsible for the exchange of metabolites, cofactors and nucleotides between the cytoplasm and mitochondrial matrix. The in silico analysis of the Drosophila melanogaster genome has highlighted the presence of 48 genes encoding putative mitochondrial carriers, but only 20 have been functionally characterized. Despite most Drosophila mitochondrial carrier genes having human homologs and sharing with them 50% or higher sequence identity, D. melanogaster genes display peculiar differences from their human counterparts: (1) in the fruit fly, many genes encode more transcript isoforms or are duplicated, resulting in the presence of numerous subfamilies in the genome; (2) the expression of the energy-producing genes in D. melanogaster is coordinated from a motif known as Nuclear Respiratory Gene (NRG), a palindromic 8-bp sequence; (3) fruit-fly duplicated genes encoding mitochondrial carriers show a testis-biased expression pattern, probably in order to keep a duplicate copy in the genome. Here, we review the main features, biological activities and role in the metabolism of the D. melanogaster mitochondrial carriers characterized to date, highlighting similarities and differences with their human counterparts. Such knowledge is very important for obtaining an integrated view of mitochondrial function in D. melanogaster metabolism.

Functional insights from the GC-poor genomes of two aphid parasitoids, Aphidius ervi and Lysiphlebus fabarum

BACKGROUND

Parasitoid wasps have fascinating life cycles and play an important role in trophic networks, yet little is known about their genome content and function. Parasitoids that infect aphids are an important group with the potential for biological control. Their success depends on adapting to develop inside aphids and overcoming both host aphid defenses and their protective endosymbionts.

RESULTS

We present the de novo genome assemblies, detailed annotation, and comparative analysis of two closely related parasitoid wasps that target pest aphids: Aphidius ervi and Lysiphlebus fabarum (Hymenoptera: Braconidae: Aphidiinae). The genomes are small (139 and 141 Mbp) and the most AT-rich reported thus far for any arthropod (GC content: 25.8 and 23.8%). This nucleotide bias is accompanied by skewed codon usage and is stronger in genes with adult-biased expression. AT-richness may be the consequence of reduced genome size, a near absence of DNA methylation, and energy efficiency. We identify missing desaturase genes, whose absence may underlie mimicry in the cuticular hydrocarbon profile of L. fabarum. We highlight key gene groups including those underlying venom composition, chemosensory perception, and sex determination, as well as potential losses in immune pathway genes.

CONCLUSIONS

These findings are of fundamental interest for insect evolution and biological control applications. They provide a strong foundation for further functional studies into coevolution between parasitoids and their hosts. Both genomes are available at https://bipaa.genouest.org.

Mitochondrial genomic variation drives differential nuclear gene expression in discrete regions of Drosophila gene and protein interaction networks

BACKGROUND

Mitochondria perform many key roles in their eukaryotic hosts, from integrating signaling pathways through to modulating whole organism phenotypes. The > 1 billion years of nuclear and mitochondrial gene co-evolution has necessitated coordinated expression of gene products from both genomes that maintain mitochondrial, and more generally, eukaryotic cellular function. How mitochondrial DNA (mtDNA) variation modifies host fitness has proved a challenging question but has profound implications for evolutionary and medical genetics. In Drosophila, we have previously shown that recently diverged mtDNA haplotypes within-species can have more impact on organismal phenotypes than older, deeply diverged haplotypes from different species. Here, we tested the effects of mtDNA haplotype variation on gene expression in Drosophila under standardized conditions. Using the Drosophila Genetic Reference Panel (DGRP), we constructed a panel of mitonuclear genotypes that consists of factorial variation in nuclear and mtDNA genomes, with mtDNAs originating in D. melanogaster (2x haplotypes) and D. simulans (2x haplotypes).

RESULTS

We show that mtDNA haplotype variation unequivocally alters nuclear gene expression in both females and males, and mitonuclear interactions are pervasive modifying factors for gene expression. There was appreciable overlap between the sexes for mtDNA-sensitive genes, and considerable transcriptional variation attributed to particular mtDNA contrasts. These genes are generally found in low-connectivity gene co-expression networks, occur in gene clusters along chromosomes, are often flanked by non-coding RNA, and are under-represented among housekeeping genes. Finally, we identify the giant (gt) transcription factor motif as a putative regulatory sequence associated with mtDNA-sensitive genes.

CONCLUSIONS

There are predictive conditions for nuclear genes that are influenced by mtDNA variation.

Cecal CircRNAs Are Associated With the Response to Salmonella Enterica Serovar Enteritidis Inoculation in the Chicken

Circular RNAs (circRNAs) are a class of endogenous noncoding RNA, which is different from linear RNA. CircRNA is an RNA molecule with a closed loop structure formed by reverse splicing. CircRNAs have been studied in several organisms, however, the circRNAs associated with the response to Salmonella enterica serovar Enteritidis (SE) inoculation in chickens are still unclear. In the current study, Jining Bairi chickens were inoculated with SE. CircRNAs involved in the response to SE inoculation were identified through next-generation sequencing. Our results showed that there were 5,118 circRNAs identified in the control and treated groups. There were 62 circRNAs significantly differentially expressed following SE inoculation. Functional classification revealed that those significantly differentially expressed circRNAs were associated with immune system process, rhythmic process and signaling following SE inoculation. CircRNAs NC_006091.4: 65510578|65515090, NC_006099.4: 16132825|16236906, and NC_006099.4: 15993284|16006290 play important roles in the response to SE inoculation. The findings in the current study provide evidence that circRNA alterations are involved in the response to SE inoculation in the chicken.

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.

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.

A Single-Cell Transcriptome Atlas of the Aging Drosophila Brain

The diversity of cell types and regulatory states in the brain, and how these change during aging, remains largely unknown. We present a single-cell transcriptome atlas of the entire adult Drosophila melanogaster brain sampled across its lifespan. Cell clustering identified 87 initial cell clusters that are further subclustered and validated by targeted cell-sorting. Our data show high granularity and identify a wide range of cell types. Gene network analyses using SCENIC revealed regulatory heterogeneity linked to energy consumption. During aging, RNA content declines exponentially without affecting neuronal identity in old brains. This single-cell brain atlas covers nearly all cells in the normal brain and provides the tools to study cellular diversity alongside other Drosophila and mammalian single-cell datasets in our unique single-cell analysis platform: SCope (http://scope.aertslab.org). These results, together with SCope, allow comprehensive exploration of all transcriptional states of an entire aging brain.

•

A single-cell atlas of the adult fly brain during aging

•

Network inference reveals regulatory states related to oxidative phosphorylation

•

Cell identity is retained during aging despite exponential decline of gene expression

•

SCope: An online tool to explore and compare single-cell datasets across species

The genome of the Hi5 germ cell line from Trichoplusia ni, an agricultural pest and novel model for small RNA biology

We report a draft assembly of the genome of Hi5 cells from the lepidopteran insect pest, Trichoplusia ni, assigning 90.6% of bases to one of 28 chromosomes and predicting 14,037 protein-coding genes. Chemoreception and detoxification gene families reveal T. ni-specific gene expansions that may explain its widespread distribution and rapid adaptation to insecticides. Transcriptome and small RNA data from thorax, ovary, testis, and the germline-derived Hi5 cell line show distinct expression profiles for 295 microRNA- and >393 piRNA-producing loci, as well as 39 genes encoding small RNA pathway proteins. Nearly all of the W chromosome is devoted to piRNA production, and T. ni siRNAs are not 2´-O-methylated. To enable use of Hi5 cells as a model system, we have established genome editing and single-cell cloning protocols. The T. ni genome provides insights into pest control and allows Hi5 cells to become a new tool for studying small RNAs ex vivo.

A common moth called the cabbage looper is becoming increasingly relevant to the scientific community. Its caterpillars are a serious threat to cabbage, broccoli and cauliflower crops, and they have started to resist the pesticides normally used to control them. Moreover, the insect’s germline cells – the ones that will produce sperm and eggs – are used in laboratories as ‘factories’ to artificially produce proteins of interest.

The germline cells also host a group of genetic mechanisms called RNA silencing. One of these processes is known as piRNA, and it protects the genome against ‘jumping genes’. These genetic elements can cause mutations by moving from place to place in the DNA: in germline cells, piRNA suppresses them before the genetic information is transmitted to the next generation. Not all germline cells grow equally well under experimental conditions, or are easy to use to examine piRNA mechanisms in a laboratory. The germline cells from the cabbage looper, on the other hand, have certain characteristics that would make them ideal to study piRNA in insects.

However, the genome of the moth had not yet been fully resolved. This hinders research on new ways of controlling the pest, on how to use the germline cells to produce more useful proteins, or on piRNA.

Decoding a genome requires several steps. First, the entire genetic information is broken in short sections that can then be deciphered. Next, these segments need to be ‘assembled’ – put together, and in the right order, to reconstitute the entire genome. Certain portions of the genome, which are formed of repeats of the same sections, can be difficult to assemble. Finally, the genome must be annotated: the different regions – such as the genes – need to be identified and labeled.

Here, Fu et al. assembled and annotated the genome of the cabbage looper, and in the process developed strategies that could be used for other species with a lot of repeated sequences in their genomes. Having access to the looper’s full genetic information makes it possible to use their germline cells to produce new types of proteins, for example for pharmaceutical purposes. Fu et al. went on to make working with these cells even easier by refining protocols so that modern research techniques, such as the gene-editing technology CRISPR-Cas9, can be used on the looper germline cells.

The mapping of the genome also revealed that the genes involved in removing toxins from the insects’ bodies are rapidly evolving, which may explain why the moths readily become resistant to insecticides. This knowledge could help finding new ways of controlling the pest.

Finally, the genes involved in RNA silencing were labeled: results show that an entire chromosome is the source of piRNAs. Combined with the new protocols developed by Fu et al., this could make cabbage looper germline cells the default option for any research into the piRNA mechanism. How piRNA works in the moth could inform work on human piRNA, as these processes are highly similar across the animal kingdom.

Glial lipid droplets and neurodegeneration in a Drosophila model of complex I deficiency

Mitochondrial defects associated with respiratory chain complex I deficiency lead to heterogeneous fatal syndromes. While the role of NDUFS8, an essential subunit of the core assembly of the complex I, is established in mitochondrial diseases, the mechanisms underlying neuropathology are poorly understood. We developed a Drosophila model of NDUFS8 deficiency by knocking down the expression of its fly homologue in neurons or in glial cells. Downregulating ND23 in neurons resulted in shortened lifespan, and decreased locomotion. Although total brain ATP levels were decreased, histological analysis did not reveal any signs of neurodegeneration except for photoreceptors of the retina. Interestingly, ND23 deficiency-associated phenotypes were rescued by overexpressing the glucose transporter hGluT3 demonstrating that boosting glucose metabolism in neurons was sufficient to bypass altered mitochondrial functions and to confer neuroprotection. We then analyzed the consequences of ND23 knockdown in glial cells. In contrast to neuronal knockdown, loss of ND23 in glia did not lead to significant behavioral defects nor to reduced lifespan, but induced brain degeneration, as visualized by numerous vacuoles found all over the nervous tissue. This phenotype was accompanied by the massive accumulation of lipid droplets at the cortex-neuropile boundaries, suggesting an alteration of lipid metabolism in glia. These results demonstrate that complex I deficiency triggers metabolic alterations both in neurons and glial cells which may contribute to the neuropathology.

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.

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.

Origin and diversification of leucine-rich repeat receptor-like protein kinase (LRR-RLK) genes in plants

Background

Leucine-rich repeat receptor-like protein kinases (LRR-RLKs) are the largest group of receptor-like kinases in plants and play crucial roles in development and stress responses. The evolutionary relationships among LRR-RLK genes have been investigated in flowering plants; however, no comprehensive studies have been performed for these genes in more ancestral groups. The subfamily classification of LRR-RLK genes in plants, the evolutionary history and driving force for the evolution of each LRR-RLK subfamily remain to be understood.

Results

We identified 119 LRR-RLK genes in the Physcomitrella patens moss genome, 67 LRR-RLK genes in the Selaginella moellendorffii lycophyte genome, and no LRR-RLK genes in five green algae genomes. Furthermore, these LRR-RLK sequences, along with previously reported LRR-RLK sequences from Arabidopsis thaliana and Oryza sativa, were subjected to evolutionary analyses. Phylogenetic analyses revealed that plant LRR-RLKs belong to 19 subfamilies, eighteen of which were established in early land plants, and one of which evolved in flowering plants. More importantly, we found that the basic structures of LRR-RLK genes for most subfamilies are established in early land plants and conserved within subfamilies and across different plant lineages, but divergent among subfamilies. In addition, most members of the same subfamily had common protein motif compositions, whereas members of different subfamilies showed variations in protein motif compositions. The unique gene structure and protein motif compositions of each subfamily differentiate the subfamily classifications and, more importantly, provide evidence for functional divergence among LRR-RLK subfamilies. Maximum likelihood analyses showed that some sites within four subfamilies were under positive selection.

Conclusions

Much of the diversity of plant LRR-RLK genes was established in early land plants. Positive selection contributed to the evolution of a few LRR-RLK subfamilies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0891-5) contains supplementary material, which is available to authorized users.

The potential role of sexual conflict and sexual selection in shaping the genomic distribution of Mito-nuclear genes

Mitochondrial interactions with the nuclear genome represent one of life's most important co-evolved mutualisms. In many organisms, mitochondria are maternally inherited, and in these cases, co-transmission between the mitochondrial and nuclear genes differs across different parts of the nuclear genome, with genes on the X chromosome having two-third probability of co-transmission, compared with one-half for genes on autosomes. These asymmetrical inheritance patterns of mitochondria and different parts of the nuclear genome have the potential to put certain gene combinations in inter-genomic co-adaptation or conflict. Previous work in mammals found strong evidence that the X chromosome has a dearth of genes that interact with the mitochondria (mito-nuclear genes), suggesting that inter-genomic conflict might drive genes off the X onto the autosomes for their male-beneficial effects. Here, we developed this idea to test coadaptation and conflict between mito-nuclear gene combinations across phylogenetically independent sex chromosomes on a far broader scale. We found that, in addition to therian mammals, only Caenorhabditis elegans showed an under-representation of mito-nuclear genes on the sex chromosomes. The remaining species studied showed no overall bias in their distribution of mito-nuclear genes. We discuss possible factors other than inter-genomic conflict that might drive the genomic distribution of mito-nuclear genes.

GLAD: an Online Database of Gene List Annotation for Drosophila

We present a resource of high quality lists of functionally related Drosophila genes, e.g. based on protein domains (kinases, transcription factors, etc.) or cellular function (e.g. autophagy, signal transduction). To establish these lists, we relied on different inputs, including curation from databases or the literature and mapping from other species. Moreover, as an added curation and quality control step, we asked experts in relevant fields to review many of the lists. The resource is available online for scientists to search and view, and is editable based on community input. Annotation of gene groups is an ongoing effort and scientific need will typically drive decisions regarding which gene lists to pursue. We anticipate that the number of lists will increase over time; that the composition of some lists will grow and/or change over time as new information becomes available; and that the lists will benefit the scientific community, e.g. at experimental design and data analysis stages. Based on this, we present an easily updatable online database, available at www.flyrnai.org/glad, at which gene group lists can be viewed, searched and downloaded.

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.

Drosophila nuclear factor DREF regulates the expression of the mitochondrial DNA helicase and mitochondrial transcription factor B2 but not the mitochondrial translation factor B1

DREF [DRE (DNA replication-related element)-binding factor] controls the transcription of numerous genes in Drosophila, many involved in nuclear DNA (nDNA) replication and cell proliferation, three in mitochondrial DNA (mtDNA) replication and two in mtDNA transcription termination. In this work, we have analysed the involvement of DREF in the expression of the known remaining genes engaged in the minimal mtDNA replication (d-mtDNA helicase) and transcription (the activator d-mtTFB2) machineries and of a gene involved in mitochondrial mRNA translation (d-mtTFB1). We have identified their transcriptional initiation sites and DRE sequences in their promoter regions. Gel-shift and chromatin immunoprecipitation assays demonstrate that DREF interacts in vitro and in vivo with the d-mtDNA helicase and d-mtTFB2, but not with the d-mtTFB1 promoters. Transient transfection assays in Drosophila S2 cells with mutated DRE motifs and truncated promoter regions show that DREF controls the transcription of d-mtDNA helicase and d-mtTFB2, but not that of d-mtTFB1. RNA interference of DREF in S2 cells reinforces these results showing a decrease in the mRNA levels of d-mtDNA helicase and d-mtTFB2 and no changes in those of the d-mtTFB1. These results link the genetic regulation of nuclear DNA replication with the genetic control of mtDNA replication and transcriptional activation in Drosophila.

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.

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.

Functional characterization of transcription factor motifs using cross-species comparison across large evolutionary distances

We address the problem of finding statistically significant associations between cis-regulatory motifs and functional gene sets, in order to understand the biological roles of transcription factors. We develop a computational framework for this task, whose features include a new statistical score for motif scanning, the use of different scores for predicting targets of different motifs, and new ways to deal with redundancies among significant motif–function associations. This framework is applied to the recently sequenced genome of the jewel wasp, Nasonia vitripennis, making use of the existing knowledge of motifs and gene annotations in another insect genome, that of the fruitfly. The framework uses cross-species comparison to improve the specificity of its predictions, and does so without relying upon non-coding sequence alignment. It is therefore well suited for comparative genomics across large evolutionary divergences, where existing alignment-based methods are not applicable. We also apply the framework to find motifs associated with socially regulated gene sets in the honeybee, Apis mellifera, using comparisons with Nasonia, a solitary species, to identify honeybee-specific associations.

We develop a computational pipeline for predicting the functions of transcription factor motifs, through DNA sequence analysis. The pipeline is applied to the newly sequenced genome of the jewel wasp, Nasonia vitripennis. It exploits the wealth of molecular data available in another insect species, the fruitfly Drosophila melanogaster, and uses cross-species comparison to its advantage. Our main contribution is to show how this can be done despite the large evolutionary divergence between the two species. The methodology presented here may be applied more generally to other scenarios (genomes) where comparative regulatory genomics must deal with large evolutionary divergences.

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.

Contrasting patterns of selective constraints in nuclear-encoded genes of the oxidative phosphorylation pathway in holometabolous insects and their possible role in hybrid breakdown in Nasonia

The principal energy generating system in animals is the oxidative phosphorylation (OXPHOS) pathway, which depends on the tight interaction of nuclear and mitochondrial encoded genes to function properly. Mitochondrial genes accumulate substitutions more quickly than nuclear genes, yet the impact of selection on mitochondrial genes is significantly reduced relative to nuclear genes due to the non-recombining nature of the mitochondrial genome and its predicted smaller effective population size. It has therefore been hypothesized that the nuclear encoded genes of the OXPHOS pathway are under strong selective pressure to compensate for the accumulation of deleterious nucleotide substitutions in mitochondrial encoded OXPHOS genes; a process known as compensatory co-adaptation. We evaluated this hypothesis by analyzing nuclear encoded OXPHOS genes for signatures of positive selection as well as evolutionary constraints at amino acid sites. We considered OXPHOS genes of six holometabolous insects and their orthologs from three Nasonia parasitoid wasps; hybrids of which suffer from an increased mortality rate caused by cytonuclear genic incompatibilities. Although nuclear OXPHOS genes are typically highly conserved, we found significant evidence for elevated amino acid divergence in four of the 59 studied nuclear encoded OXPHOS genes. We also found that three of these four genes, as well as six other OXPHOS genes, contain amino acid substitutions between Nasonia species at evolutionarily constrained sites. It is possible that these genes account for the reported incompatibility in Nasonia hybrids and their characterization may lead to a better understanding of the role of positive selection in the genetics of speciation.

Gene expression profiling implicates OXPHOS complexes in lifespan extension of flies over-expressing a small mitochondrial chaperone, Hsp22

Aging is a complex process accompanied by a decreased capacity to tolerate and respond to various stresses. Heat shock proteins as part of cell defense mechanisms are up-regulated following stress. In Drosophila, the mitochondrial Hsp22 is preferentially up-regulated in aged flies. Its over-expression results in an extension of lifespan and an increased resistance to stress. Hsp22 has chaperone-like activity in vitro, but the mechanism(s) by which it increases lifespan in flies are unknown. Genome-wide analysis was performed on long-lived Hsp22+ and control flies to unveil transcriptional changes brought by Hsp22. Transcriptomes obtained at 45days, 90% and 50% survival were then compared between them to focus more on genes up- or down-regulated in presence of higher levels of hsp22 mRNA. Hsp22+ flies display an up-regulation of genes mainly related to mitochondrial energy production and protein biosynthesis, two functions normally down-regulated during aging. Interestingly, among the 26 genes up-regulated in Hsp22+ flies, 7 genes encode for mitochondrial proteins, 5 of which being involved in OXPHOS complexes. Other genes that could influence aging such as CG5002, dGCC185 and GstS1 also displayed a regulation linked to Hsp22 expression. The up-regulation of genes of the OXPHOS system in Hsp22+ flies suggest that mitochondrial homeostasis is at the center of Hsp22 beneficial effects on lifespan.

Specific retention of the protostome-specific PsGEF may parallel with the evolution of mushroom bodies in insect and lophotrochozoan brains

Background

Gene gain and subsequent retention or loss during evolution may be one of the underlying mechanisms involved in generating the diversity of metazoan nervous systems. However, the causal relationships acting therein have not been studied extensively.

Results

We identified the gene PsGEF (protostome-specific GEF), which is present in all the sequenced genomes of insects and limpet but absent in those of sea anemones, deuterostomes, and nematodes. In Drosophila melanogaster, PsGEF encodes a short version of a protein with the C2 and PDZ domains, as well as a long version with the C2, PDZ, and RhoGEF domains through alternative splicing. Intriguingly, the exons encoding the RhoGEF domain are specifically deleted in the Daphnia pulex genome, suggesting that Daphnia PsGEF contains only the C2 and PDZ domains. Thus, the distribution of PsGEF containing the C2, PDZ, and RhoGEF domains among metazoans appears to coincide with the presence of mushroom bodies. Mushroom bodies are prominent neuropils involved in the processing of multiple sensory inputs as well as associative learning in the insect, platyhelminth, and annelid brains. In the adult Drosophila brain, PsGEF is expressed in mushroom bodies, antennal lobe, and optic lobe, where it is necessary for the correct axon branch formation of alpha/beta neurons in mushroom bodies. PsGEF genetically interacts with Rac1 but not other Rho family members, and the RhoGEF domain of PsGEF induces actin polymerization in the membrane, thus resulting in the membrane ruffling that is observed in cultured cells with activated forms of Rac.

Conclusion

The specific acquisition of PsGEF by the last common ancestor of protostomes followed by its retention or loss in specific animal species during evolution demonstrates that there are some structural and/or functional features common between insect and lophotrochozoan nervous systems (for example, mushroom bodies), which are absent in all deuterostomes and cnidarians. PsGEF is therefore one of genes associated with the diversity of metazoan nervous systems.

Comparative genomics allows the discovery of cis-regulatory elements in mosquitoes

The discovery and mapping of cis-regulatory elements is important for understanding regulation of gene transcription in mosquito vectors of human diseases. Genome sequence data are available for 3 species, Aedes aegypti, Anopheles gambiae, and Culex quinquefasciatus (Diptera: Culicidae), representing 2 subfamilies (Culicinae and Anophelinae) that are estimated to have diverged 145 to 200 million years ago. Comparative genomics tools were used to screen genomic DNA fragments located in the 5'-end flanking regions of orthologous genes. These analyses resulted in the identification of 137 sequences, designated "mosquito motifs," 7 to 9 nucleotides in length, representing 18 families of putative cis-regulatory elements conserved significantly among the 3 species when compared to the fruit fly, Drosophila melanogaster. Forty-one of the motifs were implicated previously in experiments as sites for binding transcription factors or functioning in the regulation of mosquito gene expression. Further analyses revealed associations between specific motifs and expression profiles, particularly in those genes that show increased or decreased mRNA abundance in females following a blood meal, and those accumulating transcription products exclusively or preferentially in the midgut, fat bodies, or ovaries. These results validate the methodology and support a relationship between the discovered motifs and the conservation of hematophagy in mosquitoes.