RESUMEN
Transcriptomes are dynamic, with cells, tissues, and body parts expressing particular sets of transcripts. Transposable elements (TEs) are a known source of transcriptome diversity; however, studies often focus on a particular type of chimeric transcript, analyze single body parts or cell types, or are based on incomplete TE annotations from a single reference genome. In this work, we have implemented a method based on de novo transcriptome assembly that minimizes the potential sources of errors while identifying a comprehensive set of gene-TE chimeras. We applied this method to the head, gut, and ovary dissected from five Drosophila melanogaster natural strains, with individual reference genomes available. We found that â¼19% of body part-specific transcripts are gene-TE chimeras. Overall, chimeric transcripts contribute a mean of 43% to the total gene expression, and they provide protein domains for DNA binding, catalytic activity, and DNA polymerase activity. Our comprehensive data set is a rich resource for follow-up analysis. Moreover, because TEs are present in virtually all species sequenced to date, their role in spatially restricted transcript expression is likely not exclusive to the species analyzed in this work.
Asunto(s)
Drosophila melanogaster , Drosophila , Animales , Femenino , Drosophila/genética , Drosophila melanogaster/genética , Transcriptoma , Ovario , Elementos Transponibles de ADN/genéticaRESUMEN
The advent of long-read sequencing technologies has allowed the generation of multiple high-quality de novo genome assemblies for multiple species, including well-known model species such as Drosophila melanogaster. Genome assemblies for multiple individuals of the same species are key to discover the genetic diversity present in natural populations, especially the one generated by transposable elements, the most common type of structural variant. Despite the availability of multiple genomic data sets for D. melanogaster populations, we lack an efficient visual tool to display different genome assemblies simultaneously. In this work, we present DrosOmics, a population genomic-oriented browser currently containing 52 high-quality reference genomes of D. melanogaster, including annotations from a highly reliable set of transposable elements, and functional transcriptomics and epigenomics data for 26 genomes. DrosOmics is based on JBrowse 2, a highly scalable platform, which allows the visualization of multiple assemblies at once, key to unraveling structural and functional features of D. melanogaster natural populations. DrosOmics is an open access browser and is freely available at http://gonzalezlab.eu/drosomics.
Asunto(s)
Elementos Transponibles de ADN , Drosophila melanogaster , Animales , Drosophila melanogaster/genética , GenómicaRESUMEN
RNA sequencing (RNAseq) methodology has experienced a burst of technological developments in the last decade, which has opened up opportunities for studying the mechanisms of adaptation to environmental factors at both the organismal and cellular level. Selecting the most suitable experimental approach for specific research questions and model systems can, however, be a challenge and researchers in ecology and evolution are commonly faced with the choice of whether to study gene expression variation in whole bodies, specific tissues, and/or single cells. A wide range of sometimes polarised opinions exists over which approach is best. Here, we highlight the advantages and disadvantages of each of these approaches to provide a guide to help researchers make informed decisions and maximise the power of their study. Using illustrative examples of various ecological and evolutionary research questions, we guide the readers through the different RNAseq approaches and help them identify the most suitable design for their own projects.
RESUMEN
Adaptive challenges that humans faced as they expanded across the globe left specific molecular footprints that can be decoded in our today's genomes. Different sets of metrics are used to identify genomic regions that have undergone selection. However, there are fewer methods capable of pinpointing the allele ultimately responsible for this selection. Here, we present PopHumanVar, an interactive online application that is designed to facilitate the exploration and thorough analysis of candidate genomic regions by integrating both functional and population genomics data currently available. PopHumanVar generates useful summary reports of prioritized variants that are putatively causal of recent selective sweeps. It compiles data and graphically represents different layers of information, including natural selection statistics, as well as functional annotations and genealogical estimations of variant age, for biallelic single nucleotide variants (SNVs) of the 1000 Genomes Project phase 3. Specifically, PopHumanVar amasses SNV-based information from GEVA, SnpEFF, GWAS Catalog, ClinVar, RegulomeDB and DisGeNET databases, as well as accurate estimations of iHS, nSL and iSAFE statistics. Notably, PopHumanVar can successfully identify known causal variants of frequently reported candidate selection regions, including EDAR in East-Asians, ACKR1 (DARC) in Africans and LCT/MCM6 in Europeans. PopHumanVar is open and freely available at https://pophumanvar.uab.cat.
Asunto(s)
Bases de Datos Genéticas , Genoma Humano/genética , Selección Genética/genética , Programas Informáticos , Adaptación Fisiológica/genética , Biología Computacional , Genómica , Humanos , Polimorfismo de Nucleótido Simple/genéticaRESUMEN
BACKGROUND: Escalation in industrialization and anthropogenic activity have resulted in an increase of pollutants released into the environment. Of these pollutants, heavy metals such as copper are particularly concerning due to their bio-accumulative nature. Due to its highly heterogeneous distribution and its dual nature as an essential micronutrient and toxic element, the genetic basis of copper tolerance is likely shaped by a complex interplay of genetic and environmental factors. RESULTS: In this study, we utilized the natural variation present in multiple populations of Drosophila melanogaster collected across Europe to screen for variation in copper tolerance. We found that latitude and the degree of urbanization at the collection sites, rather than any other combination of environmental factors, were linked to copper tolerance. While previously identified copper-related genes were not differentially expressed in tolerant vs. sensitive strains, genes involved in metabolism, reproduction, and protease induction contributed to the differential stress response. Additionally, the greatest transcriptomic and physiological responses to copper toxicity were seen in the midgut, where we found that preservation of gut acidity is strongly linked to greater tolerance. Finally, we identified transposable element insertions likely to play a role in copper stress response. CONCLUSIONS: Overall, by combining genome-wide approaches with environmental association analysis, and functional analysis of candidate genes, our study provides a unique perspective on the genetic and environmental factors that shape copper tolerance in natural D. melanogaster populations and identifies new genes, transposable elements, and physiological traits involved in this complex phenotype.
Asunto(s)
Cobre , Drosophila , Animales , Cobre/toxicidad , Drosophila melanogaster/genética , Genómica , Europa (Continente)RESUMEN
Drosophila melanogaster is a leading model in population genetics and genomics, and a growing number of whole-genome data sets from natural populations of this species have been published over the last years. A major challenge is the integration of disparate data sets, often generated using different sequencing technologies and bioinformatic pipelines, which hampers our ability to address questions about the evolution of this species. Here we address these issues by developing a bioinformatics pipeline that maps pooled sequencing (Pool-Seq) reads from D. melanogaster to a hologenome consisting of fly and symbiont genomes and estimates allele frequencies using either a heuristic (PoolSNP) or a probabilistic variant caller (SNAPE-pooled). We use this pipeline to generate the largest data repository of genomic data available for D. melanogaster to date, encompassing 271 previously published and unpublished population samples from over 100 locations in >20 countries on four continents. Several of these locations have been sampled at different seasons across multiple years. This data set, which we call Drosophila Evolution over Space and Time (DEST), is coupled with sampling and environmental metadata. A web-based genome browser and web portal provide easy access to the SNP data set. We further provide guidelines on how to use Pool-Seq data for model-based demographic inference. Our aim is to provide this scalable platform as a community resource which can be easily extended via future efforts for an even more extensive cosmopolitan data set. Our resource will enable population geneticists to analyze spatiotemporal genetic patterns and evolutionary dynamics of D. melanogaster populations in unprecedented detail.
Asunto(s)
Drosophila melanogaster , Metagenómica , Animales , Drosophila melanogaster/genética , Frecuencia de los Genes , Genética de Población , GenómicaRESUMEN
The McDonald and Kreitman test (MKT) is one of the most powerful and widely used methods to detect and quantify recurrent natural selection using DNA sequence data. Here we present iMKT (acronym for integrative McDonald and Kreitman test), a novel web-based service performing four distinct MKT types. It allows the detection and estimation of four different selection regimes -adaptive, neutral, strongly deleterious and weakly deleterious- acting on any genomic sequence. iMKT can analyze both user's own population genomic data and pre-loaded Drosophila melanogaster and human sequences of protein-coding genes obtained from the largest population genomic datasets to date. Advanced options in the website allow testing complex hypotheses such as the application example showed here: do genes located in high recombination regions undergo higher rates of adaptation? We aim that iMKT will become a reference site tool for the study of evolutionary adaptation in massive population genomics datasets, especially in Drosophila and humans. iMKT is a free resource online at https://imkt.uab.cat.
Asunto(s)
Adaptación Fisiológica/genética , Drosophila melanogaster/genética , Genoma , Recombinación Genética , Selección Genética , Análisis de Secuencia de ADN/estadística & datos numéricos , Alelos , Animales , Evolución Biológica , Conjuntos de Datos como Asunto , Frecuencia de los Genes , Humanos , Metagenómica , Polimorfismo GenéticoRESUMEN
Since the migrations that led humans to colonize Earth, our species has faced frequent adaptive challenges that have left signatures in the landscape of genetic variation and that we can identify in our today's genomes. Here, we (i) perform an outlier approach on eight different population genetic statistics for 22 non-admixed human populations of the Phase III of the 1000 Genomes Project to detect selective sweeps at different historical ages, as well as events of recurrent positive selection in the human lineage; and (ii) create PopHumanScan, an online catalog that compiles and annotates all candidate regions under selection to facilitate their validation and thoroughly analysis. Well-known examples of human genetic adaptation published elsewhere are included in the catalog, as well as hundreds of other attractive candidates that will require further investigation. Designed as a collaborative database, PopHumanScan aims to become a central repository to share information, guide future studies and help advance our understanding of how selection has modeled our genomes as a response to changes in the environment or lifestyle of human populations. PopHumanScan is open and freely available at https://pophumanscan.uab.cat.
Asunto(s)
Adaptación Fisiológica/genética , Biología Computacional/métodos , Genética de Población/métodos , Genoma Humano/genética , Selección Genética , Bases de Datos Genéticas , Evolución Molecular , Genómica/métodos , Humanos , Internet , Desequilibrio de Ligamiento , Modelos GenéticosRESUMEN
We present a survey of selection across Drosophila melanogaster embryonic anatomy. Our approach integrates genomic variation, spatial gene expression patterns, and development with the aim of mapping adaptation over the entire embryo's anatomy. Our adaptation map is based on analyzing spatial gene expression information for 5,969 genes (from text-based annotations of in situ hybridization data directly from the BDGP database, Tomancak et al. 2007) and the polymorphism and divergence in these genes (from the project DGRP, Mackay et al. 2012).The proportion of nonsynonymous substitutions that are adaptive, neutral, or slightly deleterious are estimated for the set of genes expressed in each embryonic anatomical structure using the distribution of fitness effects-alpha method (Eyre-Walker and Keightley 2009). This method is a robust derivative of the McDonald and Kreitman test (McDonald and Kreitman 1991). We also explore whether different anatomical structures differ in the phylogenetic age, codon usage, or expression bias of the genes they express and whether genes expressed in many anatomical structures show more adaptive substitutions than other genes.We found that: 1) most of the digestive system and ectoderm-derived structures are under selective constraint, 2) the germ line and some specific mesoderm-derived structures show high rates of adaptive substitution, and 3) the genes that are expressed in a small number of anatomical structures show higher expression bias, lower phylogenetic ages, and less constraint.
Asunto(s)
Drosophila melanogaster/anatomía & histología , Drosophila melanogaster/embriología , Drosophila melanogaster/genética , Sustitución de Aminoácidos/genética , Animales , Proteínas de Drosophila/genética , Evolución Molecular , Regulación del Desarrollo de la Expresión Génica/genética , Estudios de Asociación Genética/métodos , Variación Estructural del Genoma/genética , Genómica/métodos , Modelos Genéticos , Filogenia , Polimorfismo Genético/genética , Recombinación Genética/genética , Selección Genética/genética , Análisis Espacio-TemporalRESUMEN
Hill-Robertson interference (HRi) is expected to reduce the efficiency of natural selection when two or more linked selected sites do not segregate freely, but no attempt has been done so far to quantify the overall impact of HRi on the rate of adaptive evolution for any given genome. In this work, we estimate how much HRi impedes the rate of adaptive evolution in the coding genome of Drosophila melanogaster. We compiled a data set of 6,141 autosomal protein-coding genes from Drosophila, from which polymorphism levels in D. melanogaster and divergence out to D. yakuba were estimated. The rate of adaptive evolution was calculated using a derivative of the McDonald-Kreitman test that controls for slightly deleterious mutations. We find that the rate of adaptive amino acid substitution at a given position of the genome is positively correlated to both the rate of recombination and the mutation rate, and negatively correlated to the gene density of the region. These correlations are robust to controlling for each other, for synonymous codon bias and for gene functions related to immune response and testes. We show that HRi diminishes the rate of adaptive evolution by approximately 27%. Interestingly, genes with low mutation rates embedded in gene poor regions lose approximately 17% of their adaptive substitutions whereas genes with high mutation rates embedded in gene rich regions lose approximately 60%. We conclude that HRi hampers the rate of adaptive evolution in Drosophila and that the variation in recombination, mutation, and gene density along the genome affects the HRi effect.
Asunto(s)
Adaptación Biológica , Evolución Biológica , Drosophila/genética , Selección Genética , Animales , Drosophila melanogaster/genética , Genoma de los Insectos , Modelos Genéticos , Mutación , Sistemas de Lectura Abierta , Polimorfismo Genético , Recombinación GenéticaRESUMEN
The epigenome is the suite of interacting chemical marks and molecules that helps to shape patterns of development, phenotypic plasticity and gene regulation, in part due to its responsiveness to environmental stimuli. There is increasing interest in understanding the functional and evolutionary importance of this sensitivity under ecologically realistic conditions. Observations that epigenetic variation abounds in natural populations have prompted speculation that it may facilitate evolutionary responses to rapid environmental perturbations, such as those occurring under climate change. A frequent point of contention is whether epigenetic variants reflect genetic variation or are independent of it. The genome and epigenome often appear tightly linked and interdependent. While many epigenetic changes are genetically determined, the converse is also true, with DNA sequence changes influenced by the presence of epigenetic marks. Understanding how the epigenome, genome and environment interact with one another is therefore an essential step in explaining the broader evolutionary consequences of epigenomic variation. Drawing on results from experimental and comparative studies carried out in diverse plant and animal species, we synthesize our current understanding of how these factors interact to shape phenotypic variation in natural populations, with a focus on identifying similarities and differences between taxonomic groups. We describe the main components of the epigenome and how they vary within and between taxa. We review how variation in the epigenome interacts with genetic features and environmental determinants, with a focus on the role of transposable elements (TEs) in integrating the epigenome, genome and environment. And we look at recent studies investigating the functional and evolutionary consequences of these interactions. Although epigenetic differentiation in nature is likely often a result of drift or selection on stochastic epimutations, there is growing evidence that a significant fraction of it can be stably inherited and could therefore contribute to evolution independently of genetic change.
RESUMEN
Transposable elements (TEs), also known as transposons, are repetitive DNA sequences, present in virtually all organisms, that can move from one genomic position to another. TEs can be a source of mutations with important consequences for organisms. Despite their interest, its repetitive nature has made their study very challenging. However, the emergence of new sequencing technologies that allow obtaining long-read sequences, has improved the in silico de novo detection and annotation of TEs. The de novo annotation of TEs has already been performed in several organisms including the fruit fly Drosophila melanogaster. Yet, experimental validation can be used to confirm the presence of TEs in specific D. melanogaster natural populations. Here, we present a step-by-step protocol to experimentally validate by polymerase chain reaction (PCR) the presence and/or absence of TEs in natural populations of D. melanogaster. This detailed protocol has been implemented in the participant high schools of the Citizen Fly Lab activity that is part of the international citizen science project Melanogaster: Catch the Fly! ( https://melanogaster.eu ). Specifically, the students collaborate with the scientists of the European Drosophila Population Genomics Consortium (DrosEU) in the experimental validation of new genetic variants, previously identified using bioinformatic techniques.
Asunto(s)
Elementos Transponibles de ADN , Drosophila melanogaster , Humanos , Animales , Elementos Transponibles de ADN/genética , Drosophila melanogaster/genética , Reacción en Cadena de la Polimerasa , Drosophila , GenómicaRESUMEN
INTRODUCTION: Stroke, the second leading cause of death worldwide, is a complex disease influenced by many risk factors among which we can find reactive oxygen species (ROS). Since mitochondria are the main producers of cellular ROS, nowadays studies are trying to elucidate the role of these organelles and its DNA (mtDNA) variation in stroke risk. The aim of the present study was to perform a comprehensive evaluation of the association between mtDNA mutations and mtDNA content and stroke risk. MATERIAL AND METHODS: Homoplasmic and heteroplasmic mutations of the mtDNA were analysed in a case-controls study using 110 S cases and their corresponding control individuals. Mitochondrial DNA copy number (mtDNA-CN) was analysed in 73 of those case-control pairs. RESULTS: Our results suggest that haplogroup V, specifically variants m.72C > T, m.4580G > A, m.15904C > T and m.16298 T > C have a protective role in relation to stroke risk. On the contrary, variants m.73A > G, m.11719G > A and m.14766C > T appear to be genetic risk factors for stroke. In this study, we found no statistically significant association between stroke risk and mitochondrial DNA copy number. CONCLUSIONS: These results demonstrate the possible role of mtDNA genetics on the pathogenesis of stroke, probably through alterations in mitochondrial ROS production.
Asunto(s)
ADN Mitocondrial , Accidente Cerebrovascular , Humanos , ADN Mitocondrial/genética , Especies Reactivas de Oxígeno , Estudios de Casos y Controles , Predisposición Genética a la Enfermedad , Haplotipos , Mitocondrias/genética , Accidente Cerebrovascular/genéticaRESUMEN
The McDonald and Kreitman test is one of the most powerful and widely used methods to detect and quantify recurrent natural selection in DNA sequence data. One of its main limitations is the underestimation of positive selection due to the presence of slightly deleterious variants segregating at low frequencies. Although several approaches have been developed to overcome this limitation, most of them work on gene pooled analyses. Here, we present the imputed McDonald and Kreitman test (impMKT), a new straightforward approach for the detection of positive selection and other selection components of the distribution of fitness effects at the gene level. We compare imputed McDonald and Kreitman test with other widely used McDonald and Kreitman test approaches considering both simulated and empirical data. By applying imputed McDonald and Kreitman test to humans and Drosophila data at the gene level, we substantially increase the statistical evidence of positive selection with respect to previous approaches (e.g. by 50% and 157% compared with the McDonald and Kreitman test in Drosophila and humans, respectively). Finally, we review the minimum number of genes required to obtain a reliable estimation of the proportion of adaptive substitution (α) in gene pooled analyses by using the imputed McDonald and Kreitman test compared with other McDonald and Kreitman test implementations. Because of its simplicity and increased power to detect recurrent positive selection on genes, we propose the imputed McDonald and Kreitman test as the first straightforward approach for testing specific evolutionary hypotheses at the gene level. The software implementation and population genomics data are available at the web-server imkt.uab.cat.
Asunto(s)
Evolución Biológica , Selección Genética , Animales , Drosophila/genética , Evolución Molecular , Humanos , Metagenómica , Programas InformáticosRESUMEN
Previous studies of the evolution of genes expressed at different life-cycle stages of Drosophila melanogaster have not been able to disentangle adaptive from nonadaptive substitutions when using nonsynonymous sites. Here, we overcome this limitation by combining whole-genome polymorphism data from D. melanogaster and divergence data between D. melanogaster and Drosophila yakuba. For the set of genes expressed at different life-cycle stages of D. melanogaster, as reported in modENCODE, we estimate the ratio of substitutions relative to polymorphism between nonsynonymous and synonymous sites (α) and then α is discomposed into the ratio of adaptive (ωa) and nonadaptive (ωna) substitutions to synonymous substitutions. We find that the genes expressed in mid- and late-embryonic development are the most conserved, whereas those expressed in early development and postembryonic stages are the least conserved. Importantly, we found that low conservation in early development is due to high rates of nonadaptive substitutions (high ωna), whereas in postembryonic stages it is due, instead, to high rates of adaptive substitutions (high ωa). By using estimates of different genomic features (codon bias, average intron length, exon number, recombination rate, among others), we also find that genes expressed in mid- and late-embryonic development show the most complex architecture: they are larger, have more exons, more transcripts, and longer introns. In addition, these genes are broadly expressed among all stages. We suggest that all these genomic features are related to the conservation of mid- and late-embryonic development. Globally, our study supports the hourglass pattern of conservation and adaptation over the life-cycle.