Telomerase RNA in Hymenoptera (Insecta) switched to plant/ciliate-like biogenesis.

In contrast to the catalytic subunit of telomerase, its RNA subunit (TR) is highly divergent in size, sequence and biogenesis pathways across eukaryotes. Current views on TR evolution assume a common origin of TRs transcribed with RNA polymerase II in Opisthokonta (the supergroup including Animalia and Fungi) and Trypanosomida on one hand, and TRs transcribed with RNA polymerase III under the control of type 3 promoter, found in TSAR and Archaeplastida supergroups (including e.g. ciliates and Viridiplantae taxa, respectively). Here, we focus on unknown TRs in one of the largest Animalia order - Hymenoptera (Arthropoda) with more than 300 available representative genomes. Using a combination of bioinformatic and experimental approaches, we identify their TRs. In contrast to the presumed type of TRs (H/ACA box snoRNAs transcribed with RNA Polymerase II) corresponding to their phylogenetic position, we find here short TRs of the snRNA type, likely transcribed with RNA polymerase III under the control of the type 3 promoter. The newly described insect TRs thus question the hitherto assumed monophyletic origin of TRs across Animalia and point to an evolutionary switch in TR type and biogenesis that was associated with the divergence of Arthropods.

Immune system modulation & virus transmission during parasitism identified by multi-species transcriptomics of a declining insect biocontrol system.

BACKGROUND: The Argentine stem weevil (ASW, Listronotus bonariensis) is a significant pasture pest in Aotearoa New Zealand, primarily controlled by the parasitoid biocontrol agent Microctonus hyperodae. Despite providing effective control of ASW soon after release, M. hyperodae parasitism rates have since declined significantly, with ASW hypothesised to have evolved resistance to its biocontrol agent. While the parasitism arsenal of M. hyperodae has previously been investigated, revealing many venom components and an exogenous novel DNA virus Microctonus hyperodae filamentous virus (MhFV), the effects of said arsenal on gene expression in ASW during parasitism have not been examined. In this study, we performed a multi-species transcriptomic analysis to investigate the biology of ASW parasitism by M. hyperodae, as well as the decline in efficacy of this biocontrol system. RESULTS: The transcriptomic response of ASW to parasitism by M. hyperodae involves modulation of the weevil's innate immune system, flight muscle components, and lipid and glucose metabolism. The multispecies approach also revealed continued expression of venom components in parasitised ASW, as well as the transmission of MhFV to weevils during parasitism and some interrupted parasitism attempts. Transcriptomics did not detect a clear indication of parasitoid avoidance or other mechanisms to explain biocontrol decline. CONCLUSIONS: This study has expanded our understanding of interactions between M. hyperodae and ASW in a biocontrol system of critical importance to Aotearoa-New Zealand's agricultural economy. Transmission of MhFV to ASW during successful and interrupted parasitism attempts may link to a premature mortality phenomenon in ASW, hypothesised to be a result of a toxin-antitoxin system. Further research into MhFV and its potential role in ASW premature mortality is required to explore whether manipulation of this viral infection has the potential to increase biocontrol efficacy in future.

A supermatrix phylogeny of the world's bees (Hymenoptera: Anthophila).

The increasing availability of large molecular phylogenies has provided new opportunities to study the evolution of species traits, their origins and diversification, and biogeography; yet there are limited attempts to synthesise existing phylogenetic information for major insect groups. Bees (Hymenoptera: Anthophila) are a large group of insect pollinators that have a worldwide distribution, and a wide variation in ecology, morphology, and life-history traits, including sociality. For these reasons, as well as their major economic importance as pollinators, numerous molecular phylogenetic studies of family and genus-level relationships have been published, providing an opportunity to assemble a bee 'tree-of-life'. We used publicly available genetic sequence data, including phylogenomic data, reconciled to a taxonomic database, to produce a concatenated supermatrix phylogeny for the Anthophila comprising 4,586 bee species, representing 23% of species and 82% of genera. At family, subfamily, and tribe levels, support for expected relationships was robust, but between and within some genera relationships remain uncertain. Within families, sampling of genera ranged from 67 to 100% but species coverage was lower (17-41%). Our phylogeny mostly reproduces the relationships found in recent phylogenomic studies with a few exceptions. We provide a summary of these differences and the current state of molecular data available and its gaps. We discuss the advantages and limitations of this bee supermatrix phylogeny (available online at beetreeoflife.org), which may enable new insights into long standing questions about evolutionary drivers in bees, and potentially insects more generally.

Antennal transcriptome analysis of Psyttalia incisi (silvestri) (Hymenoptera: Braconidae): identification and tissue expression profiling of candidate odorant-binding protein genes.

BACKGROUND: Olfaction plays an important role in host-seeking by parasitoids, as they can sense chemical signals using sensitive chemosensory systems. Psyttalia incisi (Silvestri) (Hymenoptera: Braconidae) is the dominant parasitoid of Bactrocera dorsalis (Hendel) in fruit-producing regions of southern China. The olfactory behavior of P. incisi has been extensively studied; however, the chemosensory mechanisms of this species are not fully understood. RESULTS: Bioinformatics analysis of 64,515 unigenes from the antennal transcriptome of both male and female adults P. incisi identified 87 candidate chemosensory genes. These included 13 odorant-binding proteins (OBPs), seven gustatory receptors (GRs), 55 odorant receptors (ORs), 10 ionotropic receptors (IRs), and two sensory neuron membrane proteins (SNMPs). Phylogenetic trees were constructed to predict evolutionary relationships between these chemosensory genes in hymenopterans. Moreover, the tissue expression profiles of 13 OBPs were analyzed by quantitative real-time PCR, revealing high expression of seven OBPs (1, 3, 6, 7, 8, 12, and 13) in the antennae. CONCLUSION: This study represents the first identification of chemosensory genes and the determination of their expression patterns in different tissues of P. incisi. These results contribute to a better understanding of the function of the chemosensory system of this parasitoid species.

microRNAs in Syrista parreyssi (Hymenoptera) and Lepisma saccharina (Zygentoma) possibly involved in the mitochondrial function.

Mitochondria are essential organelles for maintaining vital cellular functions, and microRNAs (miRNAs) regulate gene expression posttranscriptionally. miRNAs exhibit tissue and time-specific patterns in mitochondria and specifically mitochondrial miRNAs (mitomiRs) can regulate the mRNA expression both originating from mitochondrial and nuclear transcription which affect mitochondrial metabolic activity and cell homeostasis. In this study, miRNAs of two insect species, Syrista parreyssi (Hymenoptera) and Lepisma saccharina (Zygentoma), were investigated for the first time. The known and possible novel miRNAs were predicted and characterized and their potential effects on mitochondrial transcription were investigated in these insect species using deep sequencing. The previously reported mitomiRs were also investigated and housekeeping miRNAs were characterized. miRNAs that are involved in mitochondrial processes such as apoptosis and signaling and that affect genes encoding the subunits of OXPHOS complexes have been identified in each species. Here, 81 and 161 novel mature miRNA candidates were bioinformatically predicted and 9 and 24 of those were aligned with reference mitogenomes of S. parreyssi and L. saccharina, respectively. As a result of RNAHybrid analysis, 51 and 69 potential targets of miRNAs were found in the mitogenome of S. parreyssi and L. saccharina, respectively. cox1 gene was the most targeted gene and cytB, rrnS, and rrnL genes were highly targeted in both of the species by novel miRNAs, hypothetically. We speculate that these novel miRNAs, originating from or targeting mitochondria, influence on rRNA genes or positively selected mitochondrial protein-coding genes. These findings may provide a new perspective in evaluating miRNAs for maintaining mitochondrial function and transcription.

Hymenoptera Genome Database: new genomes and annotation datasets for improved go enrichment and orthologue analyses.

We report an update of the Hymenoptera Genome Database (HGD; http://HymenopteraGenome.org), a genomic database of hymenopteran insect species. The number of species represented in HGD has nearly tripled, with fifty-eight hymenopteran species, including twenty bees, twenty-three ants, eleven wasps and four sawflies. With a reorganized website, HGD continues to provide the HymenopteraMine genomic data mining warehouse and JBrowse/Apollo genome browsers integrated with BLAST. We have computed Gene Ontology (GO) annotations for all species, greatly enhancing the GO annotation data gathered from UniProt with more than a ten-fold increase in the number of GO-annotated genes. We have also generated orthology datasets that encompass all HGD species and provide orthologue clusters for fourteen taxonomic groups. The new GO annotation and orthology data are available for searching in HymenopteraMine, and as bulk file downloads.

Blood Transcriptomics Identifies Multiple Gene Expression Pathways Associated with the Clinical Efficacy of Hymenoptera Venom Immunotherapy.

Allergen-specific venom immunotherapy (VIT) is a well-established therapy for Hymenoptera venom allergy (HVA). However, the precise mechanism underlying its clinical effect remains uncertain. Our study aimed to identify the molecular mechanisms associated with VIT efficiency. We prospectively included 19 patients with HVA undergoing VIT (sampled before the beginning of VIT, after reaching the maintenance dose, one year after finishing VIT, and after a sting challenge) and 9 healthy controls. RNA sequencing of whole blood was performed on an Illumina sequencing platform. Longitudinal transcriptomic profiling revealed the importance of the inhibition of the NFκB pathway and the downregulation of DUX4 transcripts for the early protection and induction of tolerance after finishing VIT. Furthermore, successful treatment was associated with inhibiting Th2, Th17, and macrophage alternative signalling pathways in synergy with the inhibition of the PPAR pathway and further silencing of the Th2 response. The immune system became activated when reaching the maintenance dose and was suppressed after finishing VIT. Finally, successful VIT restores the immune system's balance to a state similar to that of healthy individuals. Our results underline the important role of the inhibition of four pathways in the clinical effect of VIT: Th2, Th17, NFκB, and macrophage signalling. Two biomarkers specific for successful VIT, regardless of the time of sampling, were C4BPA and RPS10-NUDT3 and should be further tested as potential biomarkers.

The complete mitochondrial genome and novel gene arrangement in Nesodiprion zhejiangensis Zhou & Xiao (Hymenoptera: Diprionidae).

The complete mitochondrial genome (mitogenome) of the sawfly, Nesodiprion zhejiangensis Zhou & Xiao, was sequenced, assembled, and deposited in GenBank (Accession Number: OM501121). The 15,660 bp N. zhejiangensis mitogenome encodes for 2 ribosomal RNAs (rrnL and rrnS), 22 transfer RNAs (tRNAs), 13 protein-coding genes (PCGs), and an AT-rich region of 450 bp in length. The nucleotide composition is biased toward adenine and thymine (A + T = 81.8%). Each PCG is initiated by an ATN codon, except for cox2, which starts with a TTG. Of 13 PCGs, 9 have a TAA termination codon, while the remainder terminate with a TAG or a single T. All tRNAs have the classic cloverleaf structure, except for the dihydrouridine (DHU) arm of tRNAval, which forms a simple loop. There are 49 helices belonging to 6 domains in rrnL and 30 helices belonging to 4 domains in rrnS. In comparison to the ancestral architecture, N. zhejiangensis has the most rearranged mitogenome in Symphyta, in which rearrangement events of local inversion and transposition are identified in three gene clusters. Specifically, the main hotspot of gene rearrangement occurred between rrnS and trnY, and rearranged from rrnS-(AT-rich region)-I-Q-M-nd2-W-C-Y to rrnS-Q-W-C-nd2-I-M-(AT-rich region)-Y, involving a local inversion event of a large gene cluster and transposition events of some tRNAs. Transposition of trnA and trnR (rearranged from A-R to R-A) was observed at the nd3-nd5 gene junction while shuffling of trnP and trnT (rearranged from T-P to P-T) occurred at the nd4l-nd6 gene junction. While illegitimate inter-mtDNA recombination might explain the opposite orientations of transcription between rrnS and trnY, transposition events of tRNA in some gene blocks can be accounted for by the tandem duplication/random loss (TDRL) model. Our phylogenetic analysis suggests that N. zhejiangensis is closely related to congeneric species N. biremis and N. japonicus, which together form a sister lineage with the European pine sawfly, Neodiprion sertifer.

Genomic evidence for contrasting patterns of host-associated genetic differentiation across shared host-plant species in leaf- and bud-galling sawflies.

Resource specialization and ecological speciation arising through host-associated genetic differentiation (HAD) are frequently invoked as an explanation for the high diversity of plant-feeding insects and other organisms with a parasitic lifestyle. While genetic studies have demonstrated numerous examples of HAD in insect herbivores, the rarity of comparative studies means that we still lack an understanding of how deterministic HAD is, and whether patterns of host shifts can be predicted over evolutionary timescales. We applied genome-wide single nucleotide polymorphism and mitochondrial DNA sequence data obtained through genome resequencing to define species limits and to compare host-plant use in population samples of leaf- and bud-galling sawflies (Hymenoptera: Tenthredinidae: Nematinae) collected from seven shared willow (Salicaceae: Salix) host species. To infer the repeatability of long-term cophylogenetic patterns, we also contrasted the phylogenies of the two galler groups with each other as well as with the phylogeny of their Salix hosts estimated based on RADseq data. We found clear evidence for host specialization and HAD in both of the focal galler groups, but also that leaf gallers are more specialized to single host species compared with most bud gallers. In contrast to bud gallers, leaf gallers also exhibited statistically significant cophylogenetic signal with their Salix hosts. The observed discordant patterns of resource specialization and host shifts in two related galler groups that have radiated in parallel across a shared resource base indicate a lack of evolutionary repeatability in the focal system, and suggest that short- and long-term host use and ecological diversification in plant-feeding insects are dominated by stochasticity and/or lineage-specific effects.

Why do Hymenopteran workers drift to non-natal groups? Generalized reciprocity and the maximization of group and parental success.

Eusocial Hymenoptera are often characterized by having facultatively or obligately sterile worker castes. However, findings across an increasing number of species are that some workers are non-natal-they have 'drifted' away from where they were born and raised. Moreover, drifters are often indistinguishable from natal workers in the work and benefits provided to joined groups. This seems an evolutionary paradox of providing benefits to potentially unrelated individuals over close kin. Rather than being mistakes, drifting is proposed to be adaptive if joiners either gain inclusive fitness by preferentially moving to other kin groups or through generalized reciprocity in which exchanging workers across groups raises group-level genetic diversity and creates social heterosis. It is unclear, however, if reciprocity is unlikely because of a susceptibility to cheating. In resolving this question, a series of evolutionary simulations show: (1) Reciprocity can persist under a range of genetic assumptions and scenarios of cheating, (2) cheating almost always evolves, but can be expressed in a variety of ways that are not always predictable, (3) the inclusive fitness hypothesis is equally or more susceptible to cheating. Moreover, existing data in Hymenoptera (although not extensive) are more consistent with generalized reciprocity. This supports a hypothesis that drifting, as a phenomenon, may more often reflect maximization of group and parental fitness rather than fitness gains for the individual drifters.

Half-chromatid mutation as a possible cause of mosaic males and females in Hymenoptera and rare fertile male tortoiseshell cats.

Half-chromatid mutations occur when a single base change in a gamete is transmitted to the zygote, which, after DNA replication and cleavage, will result in a mosaic individual. These mutations will be passed on through the germ plasm and also may be expressed somatically. Half-chromatid mutation has been suggested to account for the observed lower frequency of males than expected for lethal X-linked recessive disorders in humans, such as Lesch-Nyhan syndrome, incontinentia pigmenti, and Duchene muscular dystrophy. Although attention has been paid to half-chromatid mutation in humans, it otherwise has been ignored. Here I show that half-chromatid mutation in haplodiploid organisms, such as Hymenoptera, has some interesting and important consequences: (i) since all genes follow the X-linked pattern of inheritance, half-chromatid mutations should be relatively easier to detect; (ii) recessive mutations of all viabilities may be expected; (iii) mosaics of both sexes are expected in haplodiploids with half-chromatid mutation; (iv) gynandromorphs could result from half-chromatid mutation at the sex-determination locus, in species with single-locus complementary sex-determination. Finally, half-chromatid mutation can account for the rare fertile male tortoiseshell phenotype observed in the domestic cat, Felis catus, and which still has not been fully accounted for by other mechanisms.

Differential sense and antisense expression profiles of Syrista parreyssi (Hymenoptera: Cephidae) mitochondrial transcripts.

The transcription of the mitogenome shows a unique pattern that is both similar to and different from the nuclear and bacterial patterns. Mitochondrial transcription generates five polycistronic units from three promoters in Drosophila melanogaster, and different expression levels of genes were observed in both different and, interestingly, the same polycistronic units in D. melanogaster. This study was conducted to test this phenomenon in the mitogenome of Syrista parreyssi (Hymenoptera: Cephidae). RNA isolation and DNase digestion were performed using only one whole individual, and real-time polymerase chain reaction analyses were performed with complementary DNAs of 11 gene regions using gene-specific primers. It was found that the expression level of each gene exhibited differences from each other, and some genes (e.g., cox genes, and rrnS) were interestingly expressed at significant levels in the corresponding antisense chain. Additionally, the mitogenome of S. parreyssi was found to have the capacity to encode 169 additional peptides from 13 known protein-coding genes, most of which were encoded in antisense transcript units. One of the unique findings was a potential open reading frame sequence that was potentially encoded in the antisense rrnL gene and included a conserved cox3 domain.

Sex biased expression and co-expression networks in development, using the hymenopteran Nasonia vitripennis.

Sexual dimorphism requires regulation of gene expression in developing organisms. These developmental differences are caused by differential expression of genes and isoforms. The effect of expressing a gene is also influenced by which other genes are simultaneously expressed (functional interactions). However, few studies have described how these processes change across development. We compare the dynamics of differential expression, isoform switching and functional interactions in the sexual development of the model parasitoid wasp Nasonia vitripennis, a system that permits genome wide analysis of sex bias from early embryos to adults. We find relatively little sex-bias in embryos and larvae at the gene level, but several sub-networks show sex-biased functional interactions in early developmental stages. These networks provide new candidates for hymenopteran sex determination, including histone modification. In contrast, sex-bias in pupae and adults is driven by the differential expression of genes. We observe sex-biased isoform switching consistently across development, but mostly in genes that are already differentially expressed. Finally, we discover that sex-biased networks are enriched by genes specific to the Nasonia clade, and that those genes possess the topological properties of key regulators. These findings suggest that regulators in sex-biased networks evolve more rapidly than regulators of other developmental networks.

Meta-Analysis of Public RNA Sequencing Data Revealed Potential Key Genes Associated with Reproductive Division of Labor in Social Hymenoptera and Termites.

Eusociality in insects has evolved independently many times. One of the most notable characteristics of eusociality is the reproductive division of labor. In social insects, the reproductive division of labor is accomplished by queens and workers. Transcriptome analyses of queens and workers have been conducted for various eusocial species. However, the genes that regulate the reproductive division of labor across all or multiple eusocial species have not yet been fully elucidated. Therefore, we conducted a meta-analysis using publicly available RNA-sequencing data from four major groups of social insects. In this meta-analysis, we collected 258 pairs (queen vs. worker) of RNA-sequencing data from 34 eusocial species. The meta-analysis identified a total of 20 genes that were differentially expressed in queens or workers. Out of these, 12 genes have not previously been reported to be involved in the reproductive division of labor. Functional annotation of these 20 genes in other organisms revealed that they could be regulators of behaviors and physiological states related to the reproductive division of labor. These 20 genes, revealed using massive datasets of numerous eusocial insects, may be key regulators of the reproductive division of labor.

Phylogenomics of braconid wasps (Hymenoptera, Braconidae) sheds light on classification and the evolution of parasitoid life history traits.

The parasitoid lifestyle is largely regarded as a key innovation that contributed to the evolutionary success and extreme species richness of the order Hymenoptera. Understanding the phylogenetic history of hyperdiverse parasitoid groups is a fundamental step in elucidating the evolution of biological traits linked to parasitoidism. We used a genomic-scale dataset based on ultra-conserved elements and the most comprehensive taxon sampling to date to estimate the evolutionary relationships of Braconidae, the second largest family of Hymenoptera. Based on our results, we propose Braconidae to comprise 41 extant subfamilies, confirmed a number of subfamilial placements and proposed subfamily-level taxonomic changes, notably the restoration of Trachypetinae stat. rev. and Masoninae stat. rev. as subfamilies of Braconidae, confirmation that Apozyx penyai Mason belongs in Braconidae placed in the subfamily Apozyginae and the recognition of Ichneutinae sensu stricto and Proteropinae as non-cyclostome subfamilies robustly supported in a phylogenetic context. The correlation between koinobiosis with endoparasitoidism and idiobiosis with ectoparasitoidism, long thought to be an important aspect in parasitoid life history, was formally tested and confirmed in a phylogenetic framework. Using ancestral reconstruction methods based on both parsimony and maximum likelihood, we suggest that the ancestor of the braconoid complex was a koinobiont endoparasitoid, as was that of the cyclostome sensu lato clade. Our results also provide strong evidence for one transition from endo- to ectoparasitoidism and three reversals back to endoparasitoidism within the cyclostome sensu stricto lineage. Transitions of koino- and idiobiosis were identical to those inferred for endo- versus ectoparasitoidism, except with one additional reversal back to koinobiosis in the small subfamily Rhysipolinae.

Antennal transcriptome analysis of olfactory genes and characterizations of odorant binding proteins in two woodwasps, Sirex noctilio and Sirex nitobei (Hymenoptera: Siricidae).

BACKGROUND: The woodwasp Sirex noctilio Fabricius is a major quarantine pest worldwide that was first discovered in China in 2013 and mainly harms Pinus sylvestris var. mongolica Litv.. S. nitobei Matsumura is a native species in China and is closely related to S. noctilio. Recently, the two woodwasps species were found attacking the P. sylvestris var. mongolica Litv in succession. The olfactory system is the foundation of insect behavior. Olfactory genes were identified through antennal transcriptome analysis. The expression profiles odorant binding proteins (OBPs) were analyzed with RT-qPCR. RESULTS: From our transcriptome analysis, 16 OBPs, 7 chemosensory proteins (CSPs), 41 odorant receptors (ORs), 8 gustatory receptors (GRs), 13 ionotropic receptors (IRs), and one sensory neuron membrane protein (SNMP) were identified in S. noctilio, while 15 OBPs, 6 CSPs, 43 ORs, 10 GRs, 16 IRs, and 1 SNMP were identified in S. nitobei. Most of the olfactory genes identified in two species were homologous. However, some species-specific olfactory genes were identified from the antennal transcriptomes, including SnocOBP13, SnocCSP6, SnocOR26, SnocGR2, SnocIR7 in S. noctilio and SnitGR9, SnitGR11, SnitIR17 in S. nitobei. In total, 14 OBPs were expressed primarily in the antennae. SnocOBP9 and SnitOBP9, identified as PBP homologues, were sex-biased expression in two siricid, but with different pattern. SnocOBP11 and SnitOBP11 were highly expressed in antennae and clearly expressed in external genitalia. SnocOBP7 and SnitOBP7 were highly expressed in male genitalia. SnocOBP3 and SnocOBP10 were highly expressed in female genitalia and male heads, while SnitOBP3 and SnitOBP10 did not show obvious tissue bias. CONCLUSION: We analyzed 86 and 91 olfactory genes from S. noctilio and S. nitobei, respectively. Most of the olfactory genes identified were homologous, but also some species-specific olfactory genes were identified, which indicated the similarities and differences of the molecular mechanisms between the two closely-related species. Different expression in the antennae, external genitals or heads, exhibiting an obvious sex bias, suggested their different role in recognizing sex pheromones or plant volatiles. Species-specific expression for several OBPs genes may suggest that they strengthened or lost their original function during species differentiation, resulting in olfactory differences between the two species.

Phylogenomics of Ichneumonoidea (Hymenoptera) and implications for evolution of mode of parasitism and viral endogenization.

Ichneumonoidea is one of the most diverse lineages of animals on the planet with >48,000 described species and many more undescribed. Parasitoid wasps of this superfamily are mostly beneficial insects that attack and kill other arthropods and are important for understanding diversification and the evolution of life history strategies related to parasitoidism. Further, some lineages of parasitoids within Ichneumonoidea have acquired endogenous virus elements (EVEs) that are permanently a part of the wasp's genome and benefit the wasp through host immune disruption and behavioral control. Unfortunately, understanding the evolution of viral acquisition, parasitism strategies, diversification, and host immune disruption mechanisms, is deeply limited by the lack of a robust phylogenetic framework for Ichneumonoidea. Here we design probes targeting 541 genes across 91 taxa to test phylogenetic relationships, the evolution of parasitoid strategies, and the utility of probes to capture polydnavirus genes across a diverse array of taxa. Phylogenetic relationships among Ichneumonoidea were largely well resolved with most higher-level relationships maximally supported. We noted codon use biases between the outgroups, Braconidae, and Ichneumonidae and within Pimplinae, which were largely solved through analyses of amino acids rather than nucleotide data. These biases may impact phylogenetic reconstruction and caution for outgroup selection is recommended. Ancestral state reconstructions were variable for Braconidae across analyses, but consistent for reconstruction of idiobiosis/koinobiosis in Ichneumonidae. The data suggest many transitions between parasitoid life history traits across the whole superfamily. The two subfamilies within Ichneumonidae that have polydnaviruses are supported as distantly related, providing strong evidence for two independent acquisitions of ichnoviruses. Polydnavirus capture using our designed probes was only partially successful and suggests that more targeted approaches would be needed for this strategy to be effective for surveying taxa for these viral genes. In total, these data provide a robust framework for the evolution of Ichneumonoidea.

Robertsonian rearrangements in Neotropical Meliponini karyotype evolution (Hymenoptera: Apidae: Meliponini).

Genome changes, evidenced through karyotype or nuclear genome size data, can result in reproductive isolation, diversification and speciation. The aim of this study was to understand how changes in the karyotype such as chromosome number and nuclear genome size accompanied the evolution of neotropical stingless bees, and to discuss these data in a phylogenetic context focusing on the karyotype evolution of this clade. We sampled 38 species representing the three Neotropical Meliponini groups; 35 for karyotype analyses and 16 for 1C value measurement. The chromosome number varied from 2n = 16 to 2n = 34, with distinct karyotypic formulae and the presence of a few polymorphisms, such as B chromosomes in one species and arm size differences between homologous chromosomes in two species. The mean 1C value varied from 0.31 pg to 0.92 pg. We associated empirical data on chromosome number and mean 1C value to highlight the importance of Robertsonian fusion rearrangements, leading to a decrease in chromosome number during the Neotropical Meliponini evolution. These data also allowed us to infer the independent heterochromatin amplification in several genera. Although less frequent, Melipona species with 2n = 22 represent evidence of Robertsonian fissions. We also pointed out the importance of chromosomal rearrangements that did not alter chromosome number, such as inversions and heterochromatin amplification.

Molecular cytogenetics in the study of repetitive sequences helping to understand the evolution of heterochromatin in Melipona (Hymenoptera, Meliponini).

The eukaryote genome is enriched by different types of repetitive DNA sequences and is most abundant in heterochromatin regions. Historically, no function has been assigned to these sequences, which makes them the target of studies that have demonstrated their structural and functional importance in the genome. Despite having a constant chromosome number, the genus Melipona has species with wide variation in heterochromatin content, from 8 to 73%, which is an important feature to be investigated regarding its origin and evolution. In the present study, a repetitive DNA sequence of Melipona mondury was isolated by restriction enzyme digestion. This sequence was used to hybridize chromosomes of eight Melipona species that include representatives of the four subgenera and present divergent characteristics in relation to the heterochromatin content. Considering that rDNA localization has shown differences in Melipona, 16 species of this genus were analyzed with 18S rDNA probe. Our data suggest that heterochromatin growth occurred independently in the Michmelia and Melikerria subgenera, considering that the isolated repetitive DNA sequence was shared only by the Michmelia species. Amplification possibly occurred from the centromeric region, causing the displacement of the rDNA sites to the ends of the chromosomes. The repetitive DNA sequence used is a constituent of Michmelia heterochromatin, which that arose from the common ancestor of the species of this subgenus.

Vitellogenin-like A-associated shifts in social cue responsiveness regulate behavioral task specialization in an ant.

Division of labor and task specialization explain the success of human and insect societies. Social insect colonies are characterized by division of labor, with workers specializing in brood care early and foraging later in life. Theory posits that this task switching requires shifts in responsiveness to task-related cues, yet experimental evidence is weak. Here, we show that a Vitellogenin (Vg) ortholog identified in an RNAseq study on the ant T. longispinosus is involved in this process: using phylogenetic analyses of Vg and Vg-like genes, we firstly show that this candidate gene does not cluster with the intensively studied honey bee Vg but falls into a separate Vg-like A cluster. Secondly, an experimental knockdown of Vg-like A in the fat body caused a reduction in brood care and an increase in nestmate care in young ant workers. Nestmate care is normally exhibited by older workers. We demonstrate experimentally that this task switch is at least partly based on Vg-like A-associated shifts in responsiveness from brood to worker cues. We thus reveal a novel mechanism leading to early behavioral maturation via changes in social cue responsiveness mediated by Vg-like A and associated pathways, which proximately play a role in regulating division of labor.