Metatranscriptomics and pyrosequencing facilitate discovery of potential viral natural enemies of the invasive Caribbean crazy ant, Nylanderia pubens

Antiviral immune response reveals host-specific virus infections in natural ant populations

Hosts can carry many viruses in their bodies, but not all of them cause disease. We studied ants as a social host to determine both their overall viral repertoire and the subset of actively infecting viruses across natural populations of three subfamilies: the Argentine ant (Linepithema humile, Dolichoderinae), the invasive garden ant (Lasius neglectus, Formicinae) and the red ant (Myrmica rubra, Myrmicinae). We used a dual sequencing strategy to reconstruct complete virus genomes by RNA-seq and to simultaneously determine the small interfering RNAs (siRNAs) by small RNA sequencing (sRNA-seq), which constitute the host antiviral RNAi immune response. This approach led to the discovery of 41 novel viruses in ants and revealed a host ant-specific RNAi response (21 vs. 22 nt siRNAs) in the different ant species. The efficiency of the RNAi response (sRNA/RNA read count ratio) depended on the virus and the respective ant species, but not its population. Overall, we found the highest virus abundance and diversity per population in Li. humile, followed by La. neglectus and M. rubra. Argentine ants also shared a high proportion of viruses between populations, whilst overlap was nearly absent in M. rubra. Only one of the 59 viruses was found to infect two of the ant species as hosts, revealing high host-specificity in active infections. In contrast, six viruses actively infected one ant species, but were found as contaminants only in the others. Disentangling spillover of disease-causing infection from non-infecting contamination across species is providing relevant information for disease ecology and ecosystem management.

The nuclear and mitochondrial genomes of Frieseomelitta varia - a highly eusocial stingless bee (Meliponini) with a permanently sterile worker caste

BACKGROUND

Most of our understanding on the social behavior and genomics of bees and other social insects is centered on the Western honey bee, Apis mellifera. The genus Apis, however, is a highly derived branch comprising less than a dozen species, four of which genomically characterized. In contrast, for the equally highly eusocial, yet taxonomically and biologically more diverse Meliponini, a full genome sequence was so far available for a single Melipona species only. We present here the genome sequence of Frieseomelitta varia, a stingless bee that has, as a peculiarity, a completely sterile worker caste.

RESULTS

The assembly of 243,974,526 high quality Illumina reads resulted in a predicted assembled genome size of 275 Mb composed of 2173 scaffolds. A BUSCO analysis for the 10,526 predicted genes showed that these represent 96.6% of the expected hymenopteran orthologs. We also predicted 169,371 repetitive genomic components, 2083 putative transposable elements, and 1946 genes for non-coding RNAs, largely long non-coding RNAs. The mitochondrial genome comprises 15,144 bp, encoding 13 proteins, 22 tRNAs and 2 rRNAs. We observed considerable rearrangement in the mitochondrial gene order compared to other bees. For an in-depth analysis of genes related to social biology, we manually checked the annotations for 533 automatically predicted gene models, including 127 genes related to reproductive processes, 104 to development, and 174 immunity-related genes. We also performed specific searches for genes containing transcription factor domains and genes related to neurogenesis and chemosensory communication.

CONCLUSIONS

The total genome size for F. varia is similar to the sequenced genomes of other bees. Using specific prediction methods, we identified a large number of repetitive genome components and long non-coding RNAs, which could provide the molecular basis for gene regulatory plasticity, including worker reproduction. The remarkable reshuffling in gene order in the mitochondrial genome suggests that stingless bees may be a hotspot for mtDNA evolution. Hence, while being just the second stingless bee genome sequenced, we expect that subsequent targeting of a selected set of species from this diverse clade of highly eusocial bees will reveal relevant evolutionary signals and trends related to eusociality in these important pollinators.

Near-Complete Genome Sequences of New Strains of Nylanderia Fulva Virus 1 from Solenopsis invicta

Nylanderia fulva virus 1 (NfV-1) is a single-stranded positive-sense RNA virus that infects the tawny crazy ant. Three near-complete genomes of NfV-1SI (92% to 94% nucleotide identity to reference strain NfV-1) found infecting the red imported fire ant were determined. The genomes have 10,904 to 10,918 nucleotides and include most of the coding region for the polyprotein.

Nylanderia fulva virus 1 (NfV-1) is a single-stranded positive-sense RNA virus that infects the tawny crazy ant. Three near-complete genomes of NfV-1SI (92% to 94% nucleotide identity to reference strain NfV-1) found infecting the red imported fire ant were determined. The genomes have 10,904 to 10,918 nucleotides and include most of the coding region for the polyprotein.

Season- and caste-specific variation in RNA viruses in the invasive Argentine ant European supercolony

The Argentine ant (Linepithema humile, Mayr) is a highly invasive species. Recently, several RNA viruses have been identified in samples from invasive Argentine ant colonies. Using quantitative PCR, we investigated variation in the levels of these viruses in the main European supercolony over the course of a year. We discovered that virus prevalence and amounts of viral RNA were affected by season and caste: ants had more virus types during warm versus cold months, and queens had more virus types and higher virus prevalence than did workers or males. This seasonal variation was largely due to the appearance of positive-strand RNA viruses in the summer and their subsequent disappearance in the winter. The prevalences of positive-strand RNA viruses were positively correlated with worker foraging activity. We hypothesise that during warmer months, ants are more active and more numerous and, as a result, they have more conspecific and heterospecific interactions that promote virus transmission.

Supercolonial structure of invasive populations of the tawny crazy ant Nylanderia fulva in the US

BACKGROUND

Social insects are among the most serious invasive pests in the world, particularly successful at monopolizing environmental resources to outcompete native species and achieve ecological dominance. The invasive success of some social insects is enhanced by their unicolonial structure, under which the presence of numerous queens and the lack of aggression against non-nestmates allow high worker densities, colony growth, and survival while eliminating intra-specific competition. In this study, we investigated the population genetics, colony structure and levels of aggression in the tawny crazy ant, Nylanderia fulva, which was recently introduced into the United States from South America.

RESULTS

We found that this species experienced a genetic bottleneck during its invasion lowering its genetic diversity by 60%. Our results show that the introduction of N. fulva is associated with a shift in colony structure. This species exhibits a multicolonial organization in its native range, with colonies clearly separated from one another, whereas it displays a unicolonial system with no clear boundaries among nests in its invasive range. We uncovered an absence of genetic differentiation among populations across the entire invasive range, and a lack of aggressive behaviors towards conspecifics from different nests, even ones separated by several hundreds of kilometers.

CONCLUSIONS

Overall, these results suggest that across its entire invasive range in the U.S.A., this species forms a single supercolony spreading more than 2000 km. In each invasive nest, we found several, up to hundreds, of reproductive queens, each being mated with a single male. The many reproductive queens per nests, together with the free movement of individuals between nests, leads to a relatedness coefficient among nestmate workers close to zero in introduced populations, calling into question the stability of this unicolonial system in which indirect fitness benefits to workers is apparently absent.

Meta-Omics Tools in the World of Insect-Microorganism Interactions

Microorganisms are able to influence several aspects of insects’ life, and this statement is gaining increasing strength, as research demonstrates it daily. At the same time, new sequencing technologies are now available at a lower cost per base, and bioinformatic procedures are becoming more user-friendly. This is triggering a huge effort in studying the microbial diversity associated to insects, and especially to economically important insect pests. The importance of the microbiome has been widely acknowledged for a wide range of animals, and also for insects this topic is gaining considerable importance. In addition to bacterial-associates, the insect-associated fungal communities are also gaining attention, especially those including plant pathogens. The use of meta-omics tools is not restricted to the description of the microbial world, but it can be also used in bio-surveillance, food safety assessment, or even to bring novelties to the industry. This mini-review aims to give a wide overview of how meta-omics tools are fostering advances in research on insect-microorganism interactions.

Indirect evidence of pathogen-associated altered oocyte production in queens of the invasive yellow crazy ant, Anoplolepis gracilipes, in Arnhem Land, Australia

Anoplolepis gracilipes is one of the six most widespread and pestiferous invasive ant species. Populations of this invader in Arnhem Land, Australia have been observed to decline, but the reasons behind these declines are not known. We investigated if there is evidence of a pathogen that could be responsible for killing ant queens or affecting their reproductive output. We measured queen number per nest, fecundity and fat content of queens from A. gracilipes populations in various stages of decline or expansion. We found no significant difference in any of these variables among populations. However, 23% of queens were found to have melanized nodules, a cellular immune response, in their ovaries and fat bodies. The melanized nodules found in dissected queens are highly likely to indicate the presence of pathogens or parasites capable of infecting A. gracilipes. Queens with nodules had significantly fewer oocytes in their ovaries, but nodule presence was not associated with low ant population abundances. Although the microorganism responsible for the nodules is as yet unidentified, this is the first evidence of the presence of a pathogenic microorganism in the invasive ant A. gracilipes that may be affecting reproduction.

Prospecting for viral natural enemies of the fire ant Solenopsis invicta in Argentina

Metagenomics and next generation sequencing were employed to discover new virus natural enemies of the fire ant, Solenopsis invicta Buren in its native range (i.e., Formosa, Argentina) with the ultimate goal of testing and releasing new viral pathogens into U.S. S. invicta populations to provide natural, sustainable control of this ant. RNA was purified from worker ants from 182 S. invicta colonies, which was pooled into 4 groups according to location. A library was created from each group and sequenced using Illumina Miseq technology. After a series of winnowing methods to remove S. invicta genes, known S. invicta virus genes, and all other non-virus gene sequences, 61,944 unique singletons were identified with virus identity. These were assembled de novo yielding 171 contiguous sequences with significant identity to non-plant virus genes. Fifteen contiguous sequences exhibited very high expression rates and were detected in all four gene libraries. One contig (Contig_29) exhibited the highest expression level overall and across all four gene libraries. Random amplification of cDNA ends analyses expanded this contiguous sequence yielding a complete virus genome, which we have provisionally named Solenopsis invicta virus 5 (SINV-5). SINV-5 is a positive-sense, single-stranded RNA virus with genome characteristics consistent with insect-infecting viruses from the family Dicistroviridae. Moreover, the replicative genome strand of SINV-5 was detected in worker ants indicating that S. invicta serves as host for the virus. Many additional sequences were identified that are likely of viral origin. These sequences await further investigation to determine their origins and relationship with S. invicta. This study expands knowledge of the RNA virome diversity found within S. invicta populations.

Genes and genomic processes underpinning the social lives of ants

The >15000 ant species are all highly social and show great variation in colony organization, complexity and behavior. The mechanisms by which such sociality evolved, as well as those underpinning the elaboration of ant societies since their ∼140 million year old common ancestor, have long been pondered. Here, we review recent insights generated using various genomic approaches. This includes understanding the molecular mechanisms underlying caste differentiation and the diversity of social structures, studying the impact of eusociality on genomic evolutionary rates, and investigating gene expression changes associated with differences in lifespan between castes. Furthermore, functional studies involving RNAi and CRISPR have recently been successfully applied to ants, opening the door to exciting research that promises to revolutionize the understanding of the evolution and diversification of social living.

Single-stranded RNA viruses infecting the invasive Argentine ant, Linepithema humile

Social insects host a diversity of viruses. We examined New Zealand populations of the globally widely distributed invasive Argentine ant (Linepithema humile) for RNA viruses. We used metatranscriptomic analysis, which identified six potential novel viruses in the Dicistroviridae family. Of these, three contigs were confirmed by Sanger sequencing as Linepithema humile virus-1 (LHUV-1), a novel strain of Kashmir bee virus (KBV) and Black queen cell virus (BQCV), while the others were chimeric or misassembled sequences. We extended the known sequence of LHUV-1 to confirm its placement in the Dicistroviridae and categorised its relationship to closest relatives, which were all viruses infecting Hymenoptera. We examined further for known viruses by mapping our metatranscriptomic sequences to all viral genomes, and confirmed KBV, BQCV, LHUV-1 and Deformed wing virus (DWV) presence using qRT-PCR. Viral replication was confirmed for DWV, KBV and LHUV-1. Viral titers in ants were higher in the presence of honey bee hives. Argentine ants appear to host a range of' honey bee' pathogens in addition to a virus currently described only from this invasive ant. The role of these viruses in the population dynamics of the ant remain to be determined, but offer potential targets for biocontrol approaches.

Isolation and characterization of Nylanderia fulva virus 1, a positive-sense, single-stranded RNA virus infecting the tawny crazy ant, Nylanderia fulva

We report the discovery of Nylanderia fulva virus 1 (NfV-1), the first virus identified and characterized from the ant, Nylanderia fulva. The NfV-1 genome (GenBank accession KX024775) is 10,881 nucleotides in length, encoding one large open reading frame (ORF). Helicase, protease, RNA-dependent RNA polymerase, and jelly-roll capsid protein domains were recognized within the polyprotein. Phylogenetic analysis placed NfV-1 in an unclassified clade of viruses. Electron microscopic examination of negatively stained samples revealed particles with icosahedral symmetry with a diameter of 28.7±1.1nm. The virus was detected by RT-PCR in larval, pupal, worker and queen developmental stages. However, the replicative strand of NfV-1 was only detected in larvae. Vertical transmission did not appear to occur, but horizontal transmission was facile. The inter-colonial field prevalence of NfV-1 was 52±35% with some local infections reaching 100%. NfV-1 was not detected in limited samples of other Nylanderia species or closely related ant species.

Evolutionary history of the Azteca-like mariner transposons and their host ants

Three different complete mariner elements were found in the genome of the ant Tapinoma nigerrimum. One (Tnigmar-Mr) was interrupted by a 900-bp insertion that corresponded to an incomplete member of a fourth mariner element, called Azteca. In this work, we isolate and characterize full-length Tnigmar-Az elements in T. nigerrimum. The purpose of this study is to clarify the evolutionary history of Azteca elements and their hosts as well as the possible existence of horizontal transfer processes. For this, Azteca-like elements were also retrieved from the available sequences of various ant genomes, representing four different ant subfamilies: Dolichoderinae, Formicinae, Myrmicinae, and Ponerinae. The tree topology resulting for the Azteca-like elements bore very little resemblance to that of their respective hosts. The pervasive presence of Azteca-like elements in all ant genomes, together with the observation that extant copies are usually younger than the genomes that host them, could be explained either by lineage sorting or by recent horizontal transfer of active elements. However, the finding of closer genetic relationships between elements than between the ants that host them is consistent with the latter scenario. This is clearly observed in Sinvmar-Az, Tnigmar-Az, Acepmar-Az, and Cflomar-Az elements, suggesting the existence of horizontal transfer processes. On the contrary, some elements displayed more divergence than did the hosts harboring them. This may reflect either further horizontal transfer events or random lineage sorting.

Myrmecomorba nylanderiae gen. et sp. nov., a microsporidian parasite of the tawny crazy ant Nylanderia fulva

A new microsporidian genus and species, Myrmecomorba nylanderiae, is described from North American populations of the tawny crazy ant, Nylanderia fulva. This new species was found to be heterosporous producing several types of binucleate spores in both larval and adult stages and an abortive octosporoblastic sporogony in adult ants. While microsporidia are widespread arthropod parasites, this description represents only the fifth species described from an ant host. Molecular analysis indicated that this new taxon is phylogenetically closely allied to the microsporidian family Caudosporidae, a group known to parasitize aquatic black fly larvae. We report the presence of 3 spore types (Type 1 DK, Type 2 DK, and octospores) with infections found in all stages of host development and reproductive castes. This report documents the first pathogen infecting N. fulva, an invasive ant of considerable economic and ecological consequence.

A metatranscriptomic approach to the identification of microbiota associated with the ant Formica exsecta

Social insects live in cooperative colonies, often in high densities and with closely related individuals, and interact using social contact behaviours. Compared to solitary insects, social insects have evolved multi-level immunity that includes immune responses common to holometabolous insects, and social immunity, which is exclusive to social taxa. This suggests that social insects may be subject to high pathogen pressure, yet relatively little is known about the range of symbiotic and pathogenic microbial communities that associate with social insects. In this study we examined transcriptome data generated from the ant Formica exsecta for sequences identifying as microbes (or other organisms potentially of non-ant origin). Sequences showing homology to two viruses and several other potentially or obligate intracellular organisms, such as Wolbachia, Arsenophonus, Entomoplasmatales and Microsporidia, were present in the transcriptome data. These homologous sequence matches correspond to genera/species that have previously been associated with a variety of insects, including social insects. There were also sequences with identity to several other microbes such as common moulds and soil bacteria. We conclude that this sequence data provides a starting point for a deeper understanding of the biological interactions between a species of ant and the micro- and macrobiotic communities that it potentially encounters.

Mining RNA-seq data for infections and contaminations

RNA sequencing (RNA-seq) provides novel opportunities for transcriptomic studies at nucleotide resolution, including transcriptomics of viruses or microbes infecting a cell. However, standard approaches for mapping the resulting sequencing reads generally ignore alternative sources of expression other than the host cell and are little equipped to address the problems arising from redundancies and gaps among sequenced microbe and virus genomes. We show that screening of sequencing reads for contaminations and infections can be performed easily using ContextMap, our recently developed mapping software. Based on mapping-derived statistics, mapping confidence, similarities and misidentifications (e.g. due to missing genome sequences) of species/strains can be assessed. Performance of our approach is evaluated on three real-life sequencing data sets and compared to state-of-the-art metagenomics tools. In particular, ContextMap vastly outperformed GASiC and GRAMMy in terms of runtime. In contrast to MEGAN4, it was capable of providing individual read mappings to species and resolving non-unique mappings, thus allowing the identification of misalignments caused by sequence similarities between genomes and missing genome sequences. Our study illustrates the importance and potentials of routinely mining RNA-seq experiments for infections or contaminations by microbes and viruses. By using ContextMap, gene expression of infecting agents can be analyzed and novel insights in infection processes and tumorigenesis can be obtained.

RNA-seq-based metatranscriptomic and microscopic investigation reveals novel metalloproteases of Neobodo sp. as potential virulence factors for soft tunic syndrome in Halocynthia roretzi

Bodonids and trypanosomatids are derived from a common ancestor with the bodonids being a more primitive lineage. The Neobodonida, one of the three clades of bodonids, can be free-living, commensal or parasitic. Despite the ecological and evolutionary significance of these organisms, however, many of their biological and pathological features are currently unknown. Here, we employed metatranscriptomics using RNA-seq technology combined with field-emission microscopy to reveal the virulence factors of a recently described genus of Neobodonida that is considered to be responsible for ascidian soft tunic syndrome (AsSTS), but whose pathogenesis is unclear. Our microscopic observation of infected tunic tissues suggested putative virulence factors, enabling us to extract novel candidate transcripts; these included cysteine proteases of the families C1 and C2, serine proteases of S51 and S9 families, and metalloproteases grouped into families M1, M3, M8, M14, M16, M17, M24, M41, and M49. Protease activity/inhibition assays and the estimation of expression levels within gene clusters allowed us to identify metalloprotease-like enzymes as potential virulence attributes for AsSTS. Furthermore, a multimarker-based phylogenetic analysis using 1,184 concatenated amino acid sequences clarified the order Neobodo sp. In sum, we herein used metatranscriptomics to elucidate the in situ expression profiles of uncharacterized putative transcripts of Neobodo sp., combined these results with microscopic observation to select candidate genes relevant to pathogenesis, and used empirical screening to define important virulence factors.

Generation and analysis of a mouse intestinal metatranscriptome through Illumina based RNA-sequencing

With the advent of high through-put sequencing (HTS), the emerging science of metagenomics is transforming our understanding of the relationships of microbial communities with their environments. While metagenomics aims to catalogue the genes present in a sample through assessing which genes are actively expressed, metatranscriptomics can provide a mechanistic understanding of community inter-relationships. To achieve these goals, several challenges need to be addressed from sample preparation to sequence processing, statistical analysis and functional annotation. Here we use an inbred non-obese diabetic (NOD) mouse model in which germ-free animals were colonized with a defined mixture of eight commensal bacteria, to explore methods of RNA extraction and to develop a pipeline for the generation and analysis of metatranscriptomic data. Applying the Illumina HTS platform, we sequenced 12 NOD cecal samples prepared using multiple RNA-extraction protocols. The absence of a complete set of reference genomes necessitated a peptide-based search strategy. Up to 16% of sequence reads could be matched to a known bacterial gene. Phylogenetic analysis of the mapped ORFs revealed a distribution consistent with ribosomal RNA, the majority from Bacteroides or Clostridium species. To place these HTS data within a systems context, we mapped the relative abundance of corresponding Escherichia coli homologs onto metabolic and protein-protein interaction networks. These maps identified bacterial processes with components that were well-represented in the datasets. In summary this study highlights the potential of exploiting the economy of HTS platforms for metatranscriptomics.