High quality reference genomes are crucial to understanding genome function, structure and evolution. The availability of reference genomes has allowed us to start inferring the role of genetic variation in biology, disease, and biodiversity conservation. However, analyses across organisms demonstrate that a single reference genome is not enough to capture the global genetic diversity present in populations. In this work, we generate 32 high-quality reference genomes for the well-known model species D. melanogaster and focus on the identification and analysis of transposable element variation as they are the most common type of structural variant. We show that integrating the genetic variation across natural populations from five climatic regions increases the number of detected insertions by 58%. Moreover, 26% to 57% of the insertions identified using long-reads were missed by short-reads methods. We also identify hundreds of transposable elements associated with gene expression variation and new TE variants likely to contribute to adaptive evolution in this species. Our results highlight the importance of incorporating the genetic variation present in natural populations to genomic studies, which is essential if we are to understand how genomes function and evolve.

Taeniid cestodes (including the human parasites Echinococcus spp. and Taenia solium) have very few mobile genetic elements (MGEs) in their genome, despite lacking a canonical PIWI pathway. The MGEs of these parasites are virtually unexplored, and nothing is known about their expression and silencing. In this work, we report the discovery of a novel family of small nonautonomous long terminal repeat retrotransposons (also known as terminal-repeat retrotransposons in miniature, TRIMs) which we have named ta-TRIM (taeniid TRIM). ta-TRIMs are only the second family of TRIM elements discovered in animals, and are likely the result of convergent reductive evolution in different taxonomic groups. These elements originated at the base of the taeniid tree and have expanded during taeniid diversification, including after the divergence of closely related species such as Echinococcus multilocularis and Echinococcus granulosus. They are massively expressed in larval stages, from a small proportion of full-length copies and from isolated terminal repeats that show transcriptional read-through into downstream regions, generating novel noncoding RNAs and transcriptional fusions to coding genes. In E. multilocularis, ta-TRIMs are specifically expressed in the germinative cells (the somatic stem cells) during asexual reproduction of metacestode larvae. This would provide a developmental mechanism for insertion of ta-TRIMs into cells that will eventually generate the adult germ line. Future studies of active and inactive ta-TRIM elements could give the first clues on MGE silencing mechanisms in cestodes.

Most signals involved in post-transcriptional regulatory networks are located in the untranslated regions (UTRs) of the mRNAs. Therefore, to deepen our understanding of gene expression regulation, delimitation of these regions with high accuracy is needed. The trypanosomatid lineage includes a variety of parasitic protozoans causing a significant worldwide burden on human health. Given their peculiar mechanisms of gene expression, these organisms depend on post-transcriptional regulation as the main level of gene expression control. In this context, the definition of the UTR regions becomes of key importance. We have developed UTR-mini-exon (UTRme), a graphical user interface (GUI) stand-alone application to identify and annotate 5′ and 3′ UTR regions in a highly accurate way. UTRme implements a multiple scoring system tailored to address the issue of false positive UTR assignment that frequently arise because of the characteristics of the intergenic regions. Even though it was developed for trypanosomatids, the tool can be used to predict 3′ sites in any eukaryote and 5′ UTRs in any organism where trans-splicing occurs (such as the model organism C. elegans). UTRme offers a way for non-bioinformaticians to precisely determine UTRs from transcriptomic data. The tool is freely available via the conda and github repositories.

Differences in social status are often mediated by agonistic encounters between competitors. Robust literature has examined social status-dependent brain gene expression profiles across vertebrates, yet social status and reproductive state are often confounded. it has therefore been challenging to identify the neuromolecular mechanisms underlying social status independent of reproductive state. Weakly electric fish, Gymnotus omarorum, display territorial aggression and social dominance independent of reproductive state. We use wild-derived G. omarorum males to conduct a transcriptomic analysis of non-breeding social dominance relationships. After allowing paired rivals to establish a dominance hierarchy, we profiled the transcriptomes of brain sections containing the preoptic area (region involved in regulating aggressive behaviour) in dominant and subordinate individuals. We identified 16 differentially expressed genes (FDR < 0.05) and numerous genes that co-varied with behavioural traits. We also compared our results with previous reports of differential gene expression in other teleost species. overall, our study establishes G. omarorum as a powerful model system for understanding the neuromolecular bases of social status independent of reproductive state.

BackgroundFasciola hepatica is the main agent of fasciolosis, a zoonotic disease affecting livestock worldwide, and an emerging food-borne disease in humans. Even when effective treatments are available, drugs are costly and can result in tolerance, liver damage and normally they do not prevent reinfection. Drug-resistant strains in livestock have been reported in various countries and, more worryingly, drug resistance in human cases has emerged in South America. The present study aims to characterize the transcriptome of two South American resistant isolates, the Cajamarca isolate from Peru, resistant to both triclabendazole and albendazole (TCBZR/ABZR) and the Rubino isolate from Uruguay, resistant to ABZ (TCBZS/ABZR), and compare them to a sensitive strain (Cenapa, Mexico, TCBZS/ABZS) to reveal putative molecular mechanisms leading to drug resistance.ResultsWe observed a major reduction in transcription in the Cajamarca TCBZR/ABZR isolate in comparison to the other isolates. While most of the differentially expressed genes are still unannotated, several trends could be detected. Specific reduction in the expression levels of cytoskeleton proteins was consistent with a role of tubulins as putative targets of triclabendazole (TCBZ). A marked reduction of adenylate cyclase might be underlying pleiotropic effects on diverse metabolic pathways of the parasite. Upregulation of GST mu isoforms suggests this detoxifying mechanism as one of the strategies associated with resistance.ConclusionsOur results stress the value of transcriptomic approaches as a means of providing novel insights to advance the understanding of drug mode of action and drug resistance. The results provide evidence for pleiotropic variations in drug-resistant isolates consistent with early observations of TCBZ and ABZ effects and recent proteomic findings.Electronic supplementary materialThe online version of this article (10.1186/s13071-017-2553-2) contains supplementary material, which is available to authorized users.

High quality reference genomes are crucial to understanding genome function, structure and evolution. The availability of reference genomes has allowed us to start inferring the role of genetic variation in biology, disease, and biodiversity conservation. However, analyses across organisms demonstrate that a single reference genome is not enough to capture the global genetic diversity present in populations. In this work, we generated 32 high-quality reference genomes for the well-known model species D. melanogaster and focused on the identification and analysis of transposable element variation as they are the most common type of structural variant. We showed that integrating the genetic variation across natural populations from five climatic regions increases the number of detected insertions by 58%. Moreover, 26% to 57% of the insertions identified using long-reads were missed by short-reads methods. We also identified hundreds of transposable elements associated with gene expression variation and new TE variants likely to contribute to adaptive evolution in this species. Our results highlight the importance of incorporating the genetic variation present in natural populations to genomic studies, which is essential if we are to understand how genomes function and evolve.

Trypanosomatids regulate gene expression mainly at the post-transcriptional level through processing, exporting and stabilising mRNA and control of translation. In most eukaryotes, protein synthesis is regulated by phosphorylation of eukaryotic initiation factor 2 (eIF2) at serine 51. Phosphorylation halts overall translation by decreasing availability of initiator tRNA met to form translating ribosomes. In trypanosomatids, the N-terminus of eIF2α is extended with threonine 169 the homologous phosphorylated residue. Here, we evaluated whether eIF2α phosphorylation varies during the Trypanosoma cruzi life cycle, the etiological agent of Chagas' disease. Total levels of eIF2α are diminished in infective and non-replicative trypomastigotes compared with proliferative forms from the intestine of the insect vector or amastigotes from mammalian cells, consistent with decreased protein synthesis reported in infective forms. eIF2α phosphorylation increases in proliferative intracellular forms prior to differentiation into trypomastigotes. Parasites overexpressing eIF2α T169A or with an endogenous CRISPR/Cas9-generated eIF2α T169A mutation were created and analysis revealed alterations to the proteome, largely unrelated to the presence of μORF in epimastigotes. eIF2α T169A mutant parasites produced fewer trypomastigotes with lower infectivity than wild type, with increased levels of sialylated mucins and oligomannose glycoproteins, and decreased galactofuranose epitopes and the surface protease GP63 on the cell surface. We conclude that eIF2α expression and phosphorylation levels affect proteins relevant for intracellular progression of T. cruzi.

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.