Assemblies of the genomes of parasitic wasps using meta-assembly and scaffolding with genetic linkage.

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.

Response to selection for parasitism of a suboptimal, low-preference host in an aphid parasitoid.

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 F 1 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.

Maize unstable factor for orange1 is essential for endosperm development and carbohydrate accumulation.

Maize (Zea mays L.) Ufo1-1 is a spontaneous dominant mutation of the unstable factor for orange1 (ufo1). We recently cloned ufo1, which is a Poaceae-specific gene highly expressed during seed development in maize. Here, we have characterized Ufo1-1 and a loss-of-function Ds insertion allele (ufo1-Dsg) to decipher the role of ufo1 in maize. We found that both ufo1 mutant alleles impact sugars and hormones, and have defects in the basal endosperm transfer layer (BETL) and adjacent cell types. The Ufo1-1 BETL had reduced cell elongation and cell wall ingrowth, resulting in cuboidal shaped transfer cells. In contrast, the ufo1-Dsg BETL cells showed a reduced overall size with abnormal wall ingrowth. Expression analysis identified the impact of ufo1 on several genes essential for BETL development. The overexpression of Ufo1-1 in various tissues leads to ectopic phenotypes, including abnormal cell organization and stomata subsidiary cell defects. Interestingly, pericarp and leaf transcriptomes also showed that as compared with wild type, Ufo1-1 had ectopic expression of endosperm development-specific genes. This study shows that Ufo1-1 impacts the expression patterns of a wide range of genes involved in various developmental processes.

Immune Suppressive Extracellular Vesicle Proteins of Leptopilina heterotoma Are Encoded in the Wasp Genome.

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 Dominant and Poorly Penetrant Phenotypes of Maize Unstable factor for orange1 Are Caused by DNA Methylation Changes at a Linked Transposon.

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.

Overlapping RdDM and non-RdDM mechanisms work together to maintain somatic repression of a paramutagenic epiallele of maize pericarp color1.

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.