Leptopilina heterotoma are obligate parasitoid wasps that develop in the body of their Drosophila hosts. During oviposition, female wasps introduce venom into the larval hosts’ body cavity. The venom contains discrete, 300 nm-wide, mixed-strategy extracellular vesicles (MSEVs), until recently referred to as virus-like particles. While the crucial immune suppressive functions of L. heterotoma MSEVs have remained undisputed, their biotic nature and origin still remain controversial. In recent proteomics analyses of L. heterotoma MSEVs, we identified 161 proteins in three classes: conserved eukaryotic proteins, infection and immunity related proteins, and proteins without clear annotation. Here we report 246 additional proteins from the L. heterotoma MSEV proteome. An enrichment analysis of the entire proteome supports vesicular nature of these structures. Sequences for more than 90% of these proteins are present in the whole-body transcriptome. Sequencing and de novo assembly of the 460 Mb-sized L. heterotoma genome revealed 90% of MSEV proteins have coding regions within the genomic scaffolds. Altogether, these results explain the stable association of MSEVs with their wasps, and like other wasp structures, their vertical inheritance. While our results do not rule out a viral origin of MSEVs, they suggest that a similar strategy for co-opting cellular machinery for immune suppression may be shared by other wasps to gain advantage over their hosts. These results are relevant to our understanding of the evolution of figitid and related wasp species.

The maize (Zea mays) mutant Unstable factor for orange1 (Ufo1) has been implicated in the epigenetic modifications of pericarp color1 (p1), which regulates the production of the flavonoid pigments phlobaphenes. Here, we show that the ufo1 gene maps to a genetically recalcitrant region near the centromere of chromosome 10. Transcriptome analysis of Ufo1-1 mutant and wild-type plants identified a candidate gene in the mapping region using a comparative sequence-based approach. The candidate gene, GRMZM2G053177, is overexpressed by >45-fold in multiple tissues of Ufo1-1, explaining the dominance of Ufo1-1 and its phenotypes. In the mutant stock, GRMZM2G053177 has a unique transcript originating within a CACTA transposon inserted in its first intron, and it is missing the first four codons of the wild-type transcript. GRMZM2G053177 expression is regulated by the DNA methylation status of the CACTA transposon, explaining the incomplete penetrance and poor expressivity of Ufo1-1. Transgenic overexpression lines of GRMZM2G053177 (Ufo1-1) phenocopy the p1-induced pigmentation in coleoptiles, tassels, leaf sheaths, husks, pericarps, and cob glumes. Transcriptome analysis of Ufo1 versus wild-type tissues revealed changes in several pathways related to abiotic and biotic stress. Thus, this study addresses the enigma of Ufo1 identity in maize, which had gone unsolved for more than 50 years.

Risks of postintroduction evolution in insects introduced to control invasive pests have been discussed for some time, but little is known about responses to selection or genetic architectures of host adaptation and thus about the likelihood or rapidity of evolutionary shifts. We report here results on the response to selection and genetic architecture of parasitism of a suboptimal, low‐preference host species by an aphid parasitoid, Aphelinus rhamni, a candidate for introduction against the soy bean aphid, Aphis glycines. We selected A. rhamni for increased parasitism of Rhopalsiphum padi by rearing the parasitoid on this aphid for three generations. We measured parasitism of R. padi at generations 2 and 3, and at generation 3, we crossed and backcrossed parasitoids from the populations reared on R. padi with those from populations reared on Aphis glycines and compared parasitism of both R. padi and Aphis glycines among F1 and backcross females. Aphelinus rhamni responded rapidly to selection for parasitism of R. padi. Selection for R. padi parasitism reduced parasitism of Aphis glycines, the original host of A. rhamni. However, parasitism of R. padi did not increase from generation 2 to generation 3 of selection, suggesting reduced variance available for selection, which was indeed found. We tested the associations between 184 single nucleotide polymorphisms (SNP) and increased parasitism of R. padi and found 28 SNP loci, some of which were associated with increased and others with decreased parasitism of R. padi. We assembled and annotated the A. rhamni genome, mapped all SNP loci to contigs and tested whether genes on contigs with SNP loci associated with parasitism were enriched for candidate genes or gene functions. We identified 80 genes on these contigs that mapped to 1.2 Mb of the 483 Mb genome of A. rhamni but found little enrichment of candidate genes or gene functions.

Allelic variation at the Zea mays (maize) pericarp color1 (p1) gene has been attributed to epigenetic gene regulation. A p1 distal enhancer, 5.2 kb upstream of the transcriptional start site, has demonstrated variation in DNA methylation in different p1 alleles/epialleles. In addition, DNA methylation of sequences within the 3’ end of intron 2 also plays a role in tissue-specific expression of p1 alleles. We show here a direct evidence for small RNAs’ involvement in regulating p1 that has not been demonstrated previously. The role of mediator of paramutation1 (mop1) was tested in the maintenance of somatic silencing at distinct p1 alleles: the non-paramutagenic P1-wr allele and paramutagenic P1-rr’ epiallele. The mop1-1 mutation gradually relieves the silenced phenotype after multiple generations of exposure; P1-wr;mop1-1 plants display a loss of 24-nt small RNAs and DNA methylation in the 3’ end of the intron 2, a region close to a Stowaway transposon. In addition, a MULE sequence within the proximal promoter of P1-wr shows depletion of 24nt siRNAs in mop1-1 plants. Release of silencing was not correlated with small RNAs at the distal enhancer region of the P1-wr allele. We found that the somatic silencing of the paramutagenic P1-rr’ is correlated with significantly reduced H3K9me2 in the distal enhancer of P1-rr’; mop1-1 plants, while symmetric DNA methylation is not significantly different. This study highlights that the epigenetic regulation of p1 alleles is controlled both via RdDM as well as non-RdDM mechanisms.

Safe, effective biological-control introductions against invasive pests depend on narrowly host-specific natural enemies with the ability to adapt to a changing environment. As part of a project on the genetic architectures of these traits, we assembled and annotated the genomes of two aphid parasitoids, Aphelinus atriplicis and Aphelinus certus. We report here several assemblies of A. atriplicis made with Illumina and PacBio data, which we combined into a meta-assembly. We scaffolded the meta-assembly with markers from a genetic map of hybrids between A. atriplicis and A. certus. We used this genetic-linkage scaffolded (GLS) assembly of A. atriplicis to scaffold a de novo assembly of A. certus. The de novo assemblies of A. atriplicis differed in contiguity, and the meta-assembly of these assemblies was more contiguous than the best de novo assembly. Scaffolding with genetic-linkage data allowed chromosomal-level assembly of the A. atriplicis genome and scaffolding a de novo assembly of A. certus with this GLS assembly, greatly increased the contiguity of the A. certus assembly to the point where it was also at the chromosomal-level. However, completeness of the A. atriplicis assembly, as measured by percent complete, single-copy BUSCO hymenopteran genes, varied little among de novo assemblies and was not increased by meta-assembly or genetic scaffolding. Furthermore, the greater contiguity of the meta-assembly and GLS assembly had little or no effect on the numbers of genes identified, the proportions with homologs or functional annotations. Increased contiguity of the A. certus assembly provided modest improvement in assembly completeness, as measured by percent complete, single-copy BUSCO hymenopteran genes. The total genic sequence increased, and while the number of genes declined, gene length increased, which together suggest greater accuracy of gene models. More contiguous assemblies provide uses other than gene annotation, for example, identifying the genes associated with quantitative trait loci and understanding of chromosomal rearrangements associated with speciation.