Filling in gaps of Drosophila melanogaster urate degradation metabolic pathway using metabolomics approaches: towards the core metabolome of the fruit fly

The primary goal of systems biology is to integrate complex omics data, and data obtained from traditional experimental studies in order to provide a holistic understanding of organismal function. One way of achieving this aim is to generate genome-scale metabolic models (GEMs), which contain information on all metabolites, enzyme-coding genes, and biochemical reactions in a biological system. Drosophila melanogaster GEM has not been reconstructed to date. Constraint-free genome-wide metabolic model of the fruit fly has been reconstructed in our lab, identifying gaps, where no enzyme was identified and metabolites were either only produced or consume. The main focus of the work presented in this thesis was to develop a pipeline for efficient gap filling using metabolomics approaches combined with standard reverse genetics methods, using 5-hydroxyisourate hydrolase (5-HIUH) as an example. 5-HIUH plays a role in urate degradation pathway. Inability to degrade urate can lead to inborn errors of metabolism (IEMs) in humans, including hyperuricemia. Based on sequence analysis Drosophila CG30016 gene was hypothesised to encode 5- HIUH. CG30016 knockout flies were examined to identify Malpighian tubules phenotype, and shortened lifespan might reflect kidney disorders in hyperuricemia in humans. Moreover, LC-MS analysis of mutant tubules revealed that CG30016 is involved in purine metabolism, and specifically urate degradation pathway. However, the exact role of the gene has not been identified, and the complete method for gap filling has not been developed. Nevertheless, thanks to the work presented here, we are a step closer towards the development of a gap-filling pipeline in Drosophila melanogaster GEM. Importantly, the areas that require further optimisation were identified and are the focus of future research. Moreover, LC-MS analysis confirmed that tubules rather than the whole fly were more suitable for metabolomics analysis of purine metabolism. Previously, Dow/Davies lab has generated the most complete tissue-specific transcriptomic atlas for Drosophila – FlyAtlas.org, which provides data on gene expression across multiple tissues of adult fly and larva. FlyAtlas revealed that transcripts of many genes are enriched in specific Drosophila tissues, and that it is possible to deduce the functions of individual tissues within the fly. Based on FlyAtlas data, it has become clear that the fly (like other metazoan species) must be considered as a set of tissues, each 2 with its own distinct transcriptional and functional profile. Moreover, it revealed that for about 30% of the genome, reverse genetic methods (i.e. mutation in an unknown gene followed by observation of phenotype) are only useful if specific tissues are investigated. Based on the FlyAtlas findings, we aimed to build a primary tissue-specific metabolome of the fruit fly, in order to establish whether different Drosophila tissues have different metabolomes and if they correspond to tissue-specific transcriptome of the fruit fly (FlyAtlas.org). Different fly tissues have been dissected and their metabolome elucidated using LC-MS. The results confirmed that tissue metabolomes differ significantly from each other and from the whole fly, and that some of these differences can be correlated to the tissue function. The results illustrate the need to study individual tissues as well as the whole organism. It is clear that some metabolites that play an important role in a given tissue might not be detected in the whole fly sample because their abundance is much lower in comparison to other metabolites present in all tissues, which prevent the detection of the tissue-specific compound.

Deciphering proteome dynamics using cell-type selective metabolic protein labeling in the fruit fly Drosophila melanogaster

Otto-von-Guericke-Universität Magdeburg, Fakultät für Naturwissenschaften, Dissertation, 2016

Transcriptome Profiling of Sexual Maturation and Mating in the Mediterranean Fruit Fly, Ceratitis capitata

Sexual maturation and mating in insects are generally accompanied by major physiological and behavioural changes. Many of these changes are related to the need to locate a mate and subsequently, in the case of females, to switch from mate searching to oviposition behaviour. The prodigious reproductive capacity of the Mediterranean fruit fly, Ceratitis capitata, is one of the factors that has led to its success as an invasive pest species. To identify the molecular changes related to maturation and mating status in male and female medfly, a microarray-based gene expression approach was used to compare the head transcriptomes of sexually immature, mature virgin, and mated individuals. Attention was focused on the changes in abundance of transcripts related to reproduction, behaviour, sensory perception of chemical stimulus, and immune system processes. Broad transcriptional changes were recorded during female maturation, while post-mating transcriptional changes in females were, by contrast, modest. In male medfly, transcriptional changes were consistent both during maturation and as a consequence of mating. Of particular note was the lack of the mating-induced immune responses that have been recorded for Drosophila melanogaster, that may be due to the different reproductive strategies of these species. This study, in addition to increasing our understanding of the molecular machinery behind maturation and mating in the medfly, has identified important gene targets that might be useful in the future management of this pest.

Physical and linkage maps for Drosophila serrata, a model species for studies of clinal adaptation and sexual selection

Drosophila serrata is a member of the montium group, which contains more than 98 species and until recently was considered a subgroup within the melanogaster group. This Drosophila species is an emerging model system for evolutionary quantitative genetics and has been used in studies of species borders, clinal variation and sexual selection. Despite the importance of D. serrata as a model for evolutionary research, our poor understanding of its genome remains a significant limitation. Here, we provide a first-generation gene-based linkage map and a physical map for this species. Consistent with previous studies of other drosophilids we observed strong conservation of genes within chromosome arms homologous with D. melanogaster but major differences in within-arm synteny. These resources will be a useful complement to ongoing genome sequencing efforts and QTL mapping studies in this species

An expressed sequence tag (EST) library for Drosophila serrata, a model system for sexual selection and climatic adaptation studies

The native Australian fly Drosophila serrata belongs to the highly speciose montium subgroup of the melanogaster species group. It has recently emerged as an excellent model system with which to address a number of important questions, including the evolution of traits under sexual selection and traits involved in climatic adaptation along latitudinal gradients. Understanding the molecular genetic basis of such traits has been limited by a lack of genomic resources for this species. Here, we present the first expressed sequence tag (EST) collection for D. serrata that will enable the identification of genes underlying sexually-selected phenotypes and physiological responses to environmental change and may help resolve controversial phylogenetic relationships within the montium subgroup.

Wolbachia infection alters olfactory-cued locomotion in Drosophila spp

Wolbachia pipientis is an endosymbiotic bacterium present in diverse insect species. Although it is well studied for its dramatic effects on host reproductive biology, little is known about its effects on other aspects of host biology, despite its presence in a wide array of host tissues. This study examined the effects of three Wolbachia strains on two different Drosophila species, using a laboratory performance assay for insect locomotion in response to olfactory cues. The results demonstrate that Wolbachia infection can have significant effects on host responsiveness that vary with respect to the Wolbachia strain-host species combination. The wRi strain, native to Drosophila simulans, increases the basal activity level of the host insect as well as its responsiveness to food cues. In contrast, the wMel strain and the virulent wMelPop strain, native to Drosophila melanogaster, cause slight decreases in responsiveness to food cues but do not alter basal activity levels in the host. Surprisingly, the virulent wMelPop strain has very little impact on host responsiveness in D. simulans. This novel strain-host relationship was artificially created previously by transinfection. These findings have implications for understanding the evolution and spread of Wolbachia infections in wild populations and for Wolbachia-based vector-borne disease control strategies currently being developed.

Posttranscriptional gene silencing is not compromised in the Arabidopsis CARPEL FACTORY (DICER-LIKE1) mutant, a homolog of dicer-1 from Drosophila

Posttranscriptional silencing (PTGS) in plants, nematodes, Drosophila, and perhaps all eukaryotes operates by sequence-specific degradation or translational inhibition of the target mRNA. These processes are mediated by duplexed RNA. In Drosophila and nematodes, double-stranded (ds)RNA or self-complementary RNA is processed into fragments of approximately 21 nt by Dicer-1 [1, 2]. These small interfering RNAs (siRNAs) serve as guides to target degradation of homologous single-stranded (ss)RNA [1, 3]. In some cases, the approximately 21 nt guide fragments derived from endogenous, imperfectly self-complementary RNAs cause translational inhibition of their target mRNAs, with which they have substantial, but not perfect sequence complementarity [4-6]. These small temporal RNAs (stRNAs) belong to a class of noncoding microRNAs (miRNAs), 20-24 nt in length, that are found in flies, plants, nematodes, and mammals [4, 6-12]. In nematodes, the Dicer-1 enzyme catalyzes the production of both siRNA and stRNA [2, 13-15]. Mutation of the Arabidopsis Dicer-1 homolog, CARPEL FACTORY (CAF), blocks miRNA production [1, 4, 16-18]. Here, we report that the same caf mutant does not block either PTGS or siRNA production induced by self-complementary hairpin RNA. This suggests either that this mutation only impairs miRNA formation or, more interestingly, that plants have two distinct dicer-like enzymes, one for miRNA and another for siRNAi production.

Drosophila 3′ UTRs are more complex than protein-coding sequences

The 3′ UTRs of eukaryotic genes participate in a variety of post-transcriptional (and some transcriptional) regulatory interactions. Some of these interactions are well characterised, but an undetermined number remain to be discovered. While some regulatory sequences in 3′ UTRs may be conserved over long evolutionary time scales, others may have only ephemeral functional significance as regulatory profiles respond to changing selective pressures. Here we propose a sensitive segmentation methodology for investigating patterns of composition and conservation in 3′ UTRs based on comparison of closely related species. We describe encodings of pairwise and three-way alignments integrating information about conservation, GC content and transition/transversion ratios and apply the method to three closely related Drosophila species: D. melanogaster, D. simulans and D. yakuba. Incorporating multiple data types greatly increased the number of segment classes identified compared to similar methods based on conservation or GC content alone. We propose that the number of segments and number of types of segment identified by the method can be used as proxies for functional complexity. Our main finding is that the number of segments and segment classes identified in 3′ UTRs is greater than in the same length of protein-coding sequence, suggesting greater functional complexity in 3′ UTRs. There is thus a need for sustained and extensive efforts by bioinformaticians to delineate functional elements in this important genomic fraction. C code, data and results are available upon request.

Polymorphisms in a desaturase 2 ortholog associate with cuticular hydrocarbon and male mating success variation in a natural population of Drosophila serrata

Elucidating the nature of genetic variation underlying both sexually selected traits and the fitness components of sexual selection is essential to understanding the broader consequences of sexual selection as an evolutionary process. To date, there have been relatively few attempts to connect the genetic variance in sexually selected traits with segregating DNA sequence polymorphisms. We set out to address this in a well-characterized sexual selection system - the cuticular hydrocarbons (CHCs) of Drosophila serrata - using an indirect association study design that allowed simultaneous estimation of the genetic variance in CHCs, sexual fitness and single nucleotide polymorphism (SNP) effects in an outbred population. We cloned and sequenced an ortholog of the D. melanogaster desaturase 2 gene, previously shown to affect CHC biosynthesis in D. melanogaster, and associated 36 SNPs with minor allele frequencies > 0.02 with variance in CHCs and sexual fitness. Three SNPs had significant multivariate associations with CHC phenotype (q-value < 0.05). At these loci, minor alleles had multivariate effects on CHCs that were weakly associated with the multivariate direction of sexual selection operating on these traits. Two of these SNPs had pleiotropic associations with male mating success, suggesting these variants may underlie responses to sexual selection due to this locus. There were 15 significant male mating success associations (q-value < 0.1), and interestingly, we detected a nonrandom pattern in the relationship between allele frequency and direction of effect on male mating success. The minor-frequency allele usually reduced male mating success, suggesting a positive association between male mating success and total fitness at this locus.

Zebrafish as a model of Parkinson's disease

Parkinson’s disease (PD) is the second most common neurodegenerative disease among the elderly. Its etiology is unknown and no disease-modifying drugs are available. Thus, more information concerning its pathogenesis is needed. Among other genes, mutated PTEN-induced kinase 1 (PINK1) has been linked to early-onset and sporadic PD, but its mode of action is poorly understood. Most animal models of PD are based on the use of the neurotoxin 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP). MPTP is metabolized to MPP+ by monoamine oxidase B (MAO B) and causes cell death of dopaminergic neurons in the substantia nigra in mammals. Zebrafish has been a widely used model organism in developmental biology, but is now emerging as a model for human diseases due to its ideal combination of properties. Zebrafish are inexpensive and easy to maintain, develop rapidly, breed in large quantities producing transparent embryos, and are readily manipulated by various methods, particularly genetic ones. In addition, zebrafish are vertebrate animals and results derived from zebrafish may be more applicable to mammals than results from invertebrate genetic models such as Drosophila melanogaster and Caenorhabditis elegans. However, the similarity cannot be taken for granted. The aim of this study was to establish and test a PD model using larval zebrafish. The developing monoaminergic neuronal systems of larval zebrafish were investigated. We identified and classified 17 catecholaminergic and 9 serotonergic neuron populations in the zebrafish brain. A 3-dimensional atlas was created to facilitate future research. Only one gene encoding MAO was found in the zebrafish genome. Zebrafish MAO showed MAO A-type substrate specificity, but non-A-non-B inhibitor specificity. Distribution of MAO in larval and adult zebrafish brains was both diffuse and distinctly cellular. Inhibition of MAO during larval development led to markedly elevated 5-hydroxytryptamine (serotonin, 5-HT) levels, which decreased the locomotion of the fish. MPTP exposure caused a transient loss of cells in specific aminergic cell populations and decreased locomotion. MPTP-induced changes could be rescued by the MAO B inhibitor deprenyl, suggesting a role for MAO in MPTP toxicity. MPP+ affected only one catecholaminergic cell population; thus, the action of MPP+ was more selective than that of MPTP. The zebrafish PINK1 gene was cloned in zebrafish, and morpholino oligonucleotides were used to suppress its expression in larval zebrafish. The functional domains and expression pattern of zebrafish PINK1 resembled those of other vertebrates, suggesting that zebrafish is a feasible model for studying PINK1. Translation inhibition resulted in cell loss of the same catecholaminergic cell populations as MPTP and MPP+. Inactivation of PINK1 sensitized larval zebrafish to subefficacious doses of MPTP, causing a decrease in locomotion and cell loss in one dopaminergic cell population. Zebrafish appears to be a feasible model for studying PD, since its aminergic systems, mode of action of MPTP, and functions of PINK1 resemble those of mammalians. However, the functions of zebrafish MAO differ from the two forms of MAO found in mammals. Future studies using zebrafish PD models should utilize the advantages specific to zebrafish, such as the ability to execute large-scale genetic or drug screens.

The Mechanisms of Hedgehog Signal Transduction

Studies in both vertebrates and invertebrates have identified proteins of the Hedgehog (Hh) family of secreted signaling molecules as key organizers of tissue patterning. Initially discovered in Drosophila in 1992, Hh family members have been discovered in animals with body plans as diverse as those of mammals, insects and echinoderms. In humans three related Hh genes have been identified: Sonic, Indian and Desert hedgehog (Shh, Ihh and Dhh). Transduction of the Hh signal to the cytoplasm utilizes an unusual mechanism involving consecutive repressive interactions between Hh and its receptor components, Patched (Ptc) and Smoothened (Smo). Several cytoplasmic proteins involved in Hh signal transduction are known in Drosophila, but mammalian homologs are known only for the Cubitus interruptus (Ci) transcription factor (GLI(1-3)) and for the Ci/GLI-associated protein, Suppressor of Fused (Su(fu)). In this study I analyzed the mechanisms of how the Hh receptor Ptc regulates the signal transducer Smo, and how Smo relays the Shh signal from the cell surface to the cytoplasm ultimately leading to the activation of GLI transcription factors. In Drosophila, the kinesin-like protein Costal2 (Cos2) is required for suppression of Hh target gene expression in the absence of ligand, and loss of Cos2 causes embryonic lethality. Cos2 acts by bridging Smo to the Ci. Another protein, Su(Fu) exerts a weak suppressive influence on Ci activity and loss of Su(Fu) causes subtle changes in Drosophila wing pattern. This study revealed that domains in Smo that are critical for Cos2 binding in Drosophila are dispensable for mammalian Smo function. Furthermore, by analyzing the function of Su(Fu) and the closest mouse homologs of Cos2 by protein overexpression and RNA interference I found that inhibition of the Hh response pathway in the absence of ligand does not require Cos2 activity, but instead critically depends on the activity of Su(Fu). These results indicate that a major change in the mechanism of action of a conserved signaling pathway occurred during evolution, probably through phenotypic drift made possible by the existence in some species of two parallel pathways acting between the Hh receptor and the Ci/GLI transcription factors. In a second approach to unravel Hh signaling we cloned > 90% of all human full-length protein kinase cDNAs and constructed the corresponding kinase-activity deficient mutants. Using this kinome resource as a screening tool, two kinases, MAP3K10 and DYRK2 were found to regulate Shh signaling. DYRK2 directly phosphorylated and induced the proteasome dependent degradation of the key Hh-pathway regulated transcription factor, GLI2. MAP3K10, in turn, affected GLI2 indirectly by modulating the activity of DYRK2.

Sex Determination - A Hypothesis Based On Noncoding Dna

Certain recent models of sex determination in mammals, Drosophila melanogaster, Caenorhabditis elegans, and snakes are examined in the light of the hypothesis that the relevant genetic regulatory mechanisms are similar and interrelated. The proposed key element in each of these instances is a noncoding DNA sequence, which serves as a high-affinity binding site for a repressor-like molecule regulating the activity of a major "sex-determining" gene. On this basis it is argued that, in several eukaryotes, (i) certain DNA sequences that are sex-determining are noncoding, in the sense that they are not the structural genes of a sex-determining protein; (ii) in some species these noncoding sequences are present in one sex and absent in the other, while in others their copy number or accessibility to regulatory molecules is significantly unequal between the two sexes; and (iii) this inequality determines whether the embryo develops into a male or a female.

Genome-wide association and genetic functional studies identify autism susceptibility candidate 2 gene (AUTS2) in the regulation of alcohol consumption

Alcohol consumption is a moderately heritable trait, but the genetic basis in humans is largely unknown, despite its clinical and societal importance. We report a genome-wide association study meta-analysis of approximately 2.5 million directly genotyped or imputed SNPs with alcohol consumption (gram per day per kilogram body weight) among 12 population-based samples of European ancestry, comprising 26,316 individuals, with replication genotyping in an additional 21,185 individuals. SNP rs6943555 in autism susceptibility candidate 2 gene (AUTS2) was associated with alcohol consumption at genome-wide significance (P = 4 x 10(-8) to P = 4 x 10(-9)). We found a genotype-specific expression of AUTS2 in 96 human prefrontal cortex samples (P = 0.026) and significant (P < 0.017) differences in expression of AUTS2 in whole-brain extracts of mice selected for differences in voluntary alcohol consumption. Down-regulation of an AUTS2 homolog caused reduced alcohol sensitivity in Drosophila (P < 0.001). Our finding of a regulator of alcohol consumption adds knowledge to our understanding of genetic mechanisms influencing alcohol drinking behavior.

Dosage Compensation And Sex Determination In Drosophila - Mechanism Of Measurement Of The X/A Ratio

We propose a molecular mechanism for the intra-cellular measurement of the ratio of the number of X chromosomes to the number of sets of autosomes, a process central to both sex determination and dosage compensation in Drosophila melanogaster. In addition to the two loci, da and Sxl, which have been shown by Cline (Genetics, 90, 683, 1978)and others to be involved in these processes, we postulate two other loci, one autosomal (ω) and the other, X-linked (π). The product of the autosomal locus da stimulates ω and initiates synthesis of a limited quantity of repressor. Sxl and π ,both of which are X-linked, compete for this repressor as well as for RNA polymerase. It is assumed that Sxl has lower affinity than π for repressor as well as polymerase and that the binding of polymerase to one of these sites modulates the binding affinity of the other site for the enzyme. It can be shown that as a result of these postulated interactions transcription from the Sxl site is proportional to the X/A ratio such that the levels of Sxl+ product are low in males, high in females and intermediate in the intersexes. If, as proposed by Cline, the Sxl- product is an inhibitor of X chromosome activity, this would result in dosage compensation. The model leads to the conclusion that high levels of Sxl+ product promote a female phenotype and low levels, a male phenotype. One interesting consequence of the assumptions on which the model is based is that the level of Sxl+ product in the cell, when examined as a function of increasing repressor concentration, first goes up and then decreases, yielding a bell-shaped curve. This feature of the model provides an explanation for some of the remarkable interactions among mutants at the Sxl, da and mle loci and leads to several predictions. The proposed mechanism may also have relevance to certain other problems, such as size regulation during development, which seem to involve measurement of ratios at the cellular level.

Classification of Nonenzymatic Homologues of Protein Kinases

Protein Kinase-Like Non-kinases (PKLNKs), which are closely related to protein kinases, lack the crucial catalytic aspartate in the catalytic loop, and hence cannot function as protein kinase, have been analysed. Using various sensitive sequence analysis methods, we have recognized 82 PKLNKs from four higher eukaryotic organisms, namely, Homo sapiens, Mus musculus, Rattus norvegicus, and Drosophila melanogaster. On the basis of their domain combination and function, PKLNKs have been classified mainly into four categories: (1) Ligand binding PKLNKs, (2) PKLNKs with extracellular protein-protein interaction domain, (3) PKLNKs involved in dimerization, and (4) PKLNKs with cytoplasmic protein-protein interaction module. While members of the first two classes of PKLNKs have transmembrane domain tethered to the PKLNK domain, members of the other two classes of PKLNKs are cytoplasmic in nature. The current classification scheme hopes to provide a convenient framework to classify the PKLNKs from other eukaryotes which would be helpful in deciphering their roles in cellular processes.