The First Draft Genomes of the Ant Formica Exsecta, and Its Wolbachia
Total Page:16
File Type:pdf, Size:1020Kb
Dhaygude et al. BMC Genomics (2019) 20:301 https://doi.org/10.1186/s12864-019-5665-6 RESEARCH ARTICLE Open Access The first draft genomes of the ant Formica exsecta, and its Wolbachia endosymbiont reveal extensive gene transfer from endosymbiont to host Kishor Dhaygude1* , Abhilash Nair1, Helena Johansson1, Yannick Wurm2 and Liselotte Sundström1,3 Abstract Background: Adapting to changes in the environment is the foundation of species survival, and is usually thought to be a gradual process. However, transposable elements (TEs), epigenetic modifications, and/or genetic material acquired from other organisms by means of horizontal gene transfer (HGTs), can also lead to novel adaptive traits. Social insects form dense societies, which attract and maintain extra- and intracellular accessory inhabitants, which may facilitate gene transfer between species. The wood ant Formica exsecta (Formicidae; Hymenoptera), is a common ant species throughout the Palearctic region. The species is a well-established model for studies of ecological characteristics and evolutionary conflict. Results: In this study, we sequenced and assembled draft genomes for F. exsecta and its endosymbiont Wolbachia. The F. exsecta draft genome is 277.7 Mb long; we identify 13,767 protein coding genes, for which we provide gene ontology and protein domain annotations. This is also the first report of a Wolbachia genome from ants, and provides insights into the phylogenetic position of this endosymbiont. We also identified multiple horizontal gene transfer events (HGTs) from Wolbachia to F. exsecta. Some of these HGTs have also occurred in parallel in multiple other insect genomes, highlighting the extent of HGTs in eukaryotes. Conclusion: We present the first draft genome of ant F. exsecta, and its endosymbiont Wolbachia (wFex), and show considerable rates of gene transfer from the symbiont to the host. We expect that especially the F. exsecta genome will be valuable resource in further exploration of the molecular basis of the evolution of social organization. Keywords: Formica exsecta, Genome, Endosymbionts, Transposons, Horizontal gene transfer, Wolbachia Background horizontal gene transfer (HGTs), and this can also lead to Adapting to changes in the environment is the foundation novel adaptive traits [4, 5]. Both mutations and HGTs can of species survival, and is usually thought to be a gradual drive rapid genome evolution [6, 7]. Horizontal gene trans- process. Genomic changes, such as single nucleotide sub- fers have been reported in many taxa, most commonly from stitutions play key roles in adaptive evolution, although bacteria to animals [7], plants [8, 9], fungi [10–12], but the few mutations are beneficial. Besides nucleotide substitu- mechanisms that underpin horizontal gene transfer events, tions, other structural and regulatory units, such as trans- and the mode by which bacterial genetic material is inte- posable elements (TEs) and epigenetic modifications, can grated into the eukaryote genome are not well understood. also act as drivers in adaptation [1–3]. Genetic material Many cases of horizontal gene transfer from bacteria to can also be acquired from other organisms by means of eukaryotes involve intracellular endosymbionts, which are maternally transmitted through oocytes [13, 14]. The most * Correspondence: [email protected] common examples of endosymbiont to host horizontal 1 Organismal and Evolutionary Biology Research Programme, Faculty of gene transfers involve the bacterium Wolbachia, a well de- Biological and environmental sciences, University of Helsinki, P.O. Box 65, FI-00014 Helsinki, Finland scribed intracellular, maternally inherited gram-negative Full list of author information is available at the end of the article © The Author(s). 2019 Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated. Dhaygude et al. BMC Genomics (2019) 20:301 Page 2 of 16 bacterium known to infect over 60% of the investigated in- Results & discussion sect species [15–17]. Wolbachia infection is also prevalent The F. exsecta genome in filarial nematodes, crustaceans, and arachnids [18–20]. The Illumina sequencing libraries from DNA extracted from Wolbachia- host interactions can be mutualistic or patho- testes of males of a F. exsecta colony yielded > 99 gigabases genic [21]. A number of ecdysozoan genomes have been re- of Illumina sequence data. The final genome resulting from ported to contain chromosomal insertions originating from the assembly of these data was 277.7 megabases (Mb) long, Wolbachia, including the mosquito Aedes aegypti [22, 23], encompassing 14,617 scaffolds (Fig. 1)withaN50scaffold the longhorn beetle Monochamus alternatus [24], filarial length of 997.7 kb (Table 1). The number of scaffolds is nematodes of the genera Onchocerca, Brugia,andDirofilaria higher than the number of chromosomes (n =26)reported [20, 25], parasitoid wasps of the genus Nasonia,thefruitfly for F. exsecta [38, 39]. Similarly, the F. exsecta genome as- Drosophila ananassae, the pea aphid Acyrthosiphon pisum sembly is somewhat shorter than genome size estimates ob- [26, 27], and the bean beetle Callosobruchus chinensis [28]. tained by flow cytometry for species in the subfamily Although most of the transferred DNA is probably nonfunc- Formicinae (range: 296–385 Mb) [40]. These discrepancies tional in the host genome [25, 28, 29], some of the trans- are unsurprising given the difficulty of assembling highly re- ferred genes are functional [22]. The functional HGT events petitive gene content from short sequencing reads [41]. In can be categorized into two broad types – one that maintains line with this, the genome assembly length metrics are simi- pre-existing functions in the recipient host, and one that pro- lar to those of the 23 ant genomes that have been published. vides the recipient host with new functionality, including al- The raw data, gapped scaffolds, and annotations underpin- tered host nutrition, protection and adaptation to extreme ning this assembly are deposited on public databases under environments [30]. BioProject PRJNA393850 (accession NPMM00000000). Infection with Wolbachia is widespread in Hymenop- tera. Most hymenopteran Wolbachia infections have the Quantitative assessment of genome assembly cytoplasmic incompatibility phenotype [31], which leads Based on scaffold N50 and N75 statistics, contig size, and to reproductive incompatibility between infected sperm GC content, the F. exsecta genome assembly is comparable and uninfected eggs. The ants that have been investi- in quality and completeness to other sequenced ant genomes gated so far are infected with A-group strains of (Additional file 1: Table S1). All the 248 CEGMA eukaryotic CI-inducing Wolbachia [32–34]. Wenseleers et al. [35] core genes were found, and 241 of these genes were showed that 25 out of 50 species of ants in Java and Su- complete in length. Similarly, 98.5% of 1634 BUSCO Insecta matra screened positive for a single A-group strain of genes were complete in the genome (Table 2). These results Wolbachia. By contrast, a study on a single Swiss popu- held with other BUSCO analysis levels including Eukaryota, lation of the ant Formica exsecta, found that all the ants Arthropoda, and Hymenoptera, with low duplication levels tested were infected with four or five different strains of (2.2 to 5.3%), and a few missing genes (0.6 to 1.27%; Table 2). Wolbachia [32, 36]. Such discrepancies can be due to technical artifacts such as The aims of this study are to produce the first genome sequencing biases or assembly difficulties, as well as to true for the ant genus Formica, to test whether horizontally differences between our F. exsecta sample and the BUSCO transferred genetic elements exist in the genome of the and CEGMA datasets. To further evaluate genome com- ant F. exsecta, and to describe the genomic organization pleteness, we compared the independently generated F. of any such elements. The genus Formica is listed by the exsecta transcriptome [42] to the genome reported here. Global Ant Genome Alliance (GAGA) as one of the More than 98.75% of the 10.999 assembled ESTs mapped high-priority ant taxons to be sequenced [37], owing to − unambiguously to the genome (BLASTn E < 10 50). To- its key taxonomic position, and the ecological and be- gether, these analyses show that the genome assembly has havioral data that are available for the species. We report high completeness. the first whole genome sequencing of this species, and the draft genome sequence of its associated cytoplasmic Wolbachia endosymbiont (wFex). We further report Gene content in the F. exsecta genome the presence of multiple extensive insertions of Wol- We identified 13,637 protein coding genes by com- bachia genetic material in the host genome, and com- bining ab initio, EST-based, and sequence similarity pare the HGTs insertions discovered in the assembled based gene predictions methods. The GC content draft genome to other genomes, to understand the was higher in exons (41.6%)