A chromosome‐level genome assembly of the parasitoid wasp Pteromalus puparum

Parasitoid wasps represent a large proportion of hymenopteran species. They have complex evolutionary histories and are important biocontrol agents. To advance parasitoid research, a combination of Illumina short‐read, PacBio long‐read and Hi‐C scaffolding technologies was used to develop a high‐quality chromosome‐level genome assembly for Pteromalus puparum, which is an important pupal endoparasitoid of caterpillar pests. The chromosome‐level assembly has aided in studies of venom and detoxification genes. The assembled genome size is 338 Mb with a contig N50 of 38.7 kb and a scaffold N50 of 1.16 Mb. Hi‐C analysis assembled scaffolds onto five chromosomes and raised the scaffold N50 to 65.8 Mb, with more than 96% of assembled bases located on chromosomes. Gene annotation was assisted by RNA sequencing for the two sexes and four different life stages. Analysis detected 98% of the BUSCO (Benchmarking Universal Single‐Copy Orthologs) gene set, supporting a high‐quality assembly and annotation. In total, 40.1% (135.6 Mb) of the assembly is composed of repetitive sequences, and 14,946 protein‐coding genes were identified. Although venom genes play important roles in parasitoid biology, their spatial distribution on chromosomes was poorly understood. Mapping has revealed venom gene tandem arrays for serine proteases, pancreatic lipase‐related proteins and kynurenine–oxoglutarate transaminases, which have amplified in the P. puparum lineage after divergence from its common ancestor with Nasonia vitripennis. In addition, there is a large expansion of P450 genes in P. puparum. These examples illustrate how chromosome‐level genome assembly can provide a valuable resource for molecular, evolutionary and biocontrol studies of parasitoid wasps.     

Genome-wide identification and analysis of genes encoding cuticular proteins in the endoparasitoid wasp Pteromalus puparum (Hymenoptera: Pteromalidae)

The multifunctional insect cuticle serves as the exoskeleton, determines body shape, restricts water loss, provides attachment sites for muscles and internal organs and is a formidable barrier to invaders. It is morphologically divided into three layers, including envelope, epicuticle, and procuticle and is composed of chitin and cuticular proteins (CPs). Annotation of CPs and their cognate genes may help understand the structure and functions of insect cuticles. In this paper, we interrogated the genome of Pteromalus puparum, an endoparasitoid wasp that parasitizes Pieris rapae and Papilio xuthus pupae, and identified 82 genes encoding CPs belonging to six CP families, including 62 in the CPR family, 8 in CPAP3, 5 in CPF/CPFL, 2 low complexity proteins, 2 in TWDL, and 3 in Apidermin. We used six RNA-seq libraries to determine CP gene expression profiles through development and compared the cuticle hydrophobicity between the P. puparum and the ectoparasitoid Nasonia vitripennis based on GRAVY values of CPR sequences. In the Nasonia-Pteromalus comparison, we found in both N. vitripennis and P. puparum, the peak of their CPR hydrophobicity displayed at their pupal stage, whereas their adult stage showed the lowest level. Except at the adult stage, the CPR hydrophobicity in N. vitripennis is always higher than P. puparum. Finally, we identified three novel Apidermin genes, a family found solely in Hymenoptera and revealed a new sequence feature of this family. This new information contributes to a broader understanding of insect CPs generally.     

EFFECT OF PARASITIZATION BY THE PUPAL ENDOPARISITOID,PTEROMALUS PUPARUM (HYMENOPTERA: PTEROMALIDAE) ON HUMORAL IMMUNE REACTIONS OF PIERIS RAPAE (LEPIDOPETERA: PIERIDAE)*

Abstract Studies on the effect of parasitization by the endoparasitoid on host humoral immune reactions are carried out with the pupal endoparisitic wasp, Pteromalus puparum, and its host, Pieris rapae. Phenoloxidase (PO) activity of parasitized hosts hemolymph increased significantly at 12 h, day four and day five after parasitization. Hem‐agglutination activity of parasitized hosts hemolymph was always higher than that of wounded and unparasitized ones. Moreover, antibacterial activity of parasitized hosts hemolymph became more and more stronger, whilst wounded and unparasitized pupae only owned a weak antibacterial activity. It suggested that activities of humoral immune factors of Pieris rapae could be influenced to some degrees by P. puparum.     

Comparative genomic study of arachnid immune systems indicates loss of beta‐1,3‐glucanase‐related proteins and the immune deficiency pathway

Analyses of arthropod genomes have shown that the genes in the different innate humoral immune responses are conserved. These genes encode proteins that are involved in immune signalling pathways that recognize pathogens and activate immune responses. These immune responses include phagocytosis, encapsulation of the pathogen and production of effector molecules for pathogen elimination. So far, most studies have focused on insects leaving other major arthropod groups largely unexplored. Here, we annotate the immune‐related genes of six arachnid genomes and present evidence for a conserved pattern of some immune genes, but also evolutionary changes in the arachnid immune system. Specifically, our results suggest that the family of recognition molecules of beta‐1,3‐glucanase‐related proteins (βGRPs) and the genes from the immune deficiency (IMD) signalling pathway have been lost in a common ancestor of arachnids. These findings are consistent with previous work suggesting that the humoral immune effector proteins are constitutively produced in arachnids in contrast to insects, where these have to be induced. Further functional studies are needed to verify this. We further show that the full haemolymph clotting cascade found in the horseshoe crab is retrieved in most arachnid genomes. Tetranychus lacks at least one major component, although it is possible that this cascade could still function through recruitment of a different protein. The gel‐forming protein in horseshoe crabs, coagulogen, was not recovered in any of the arachnid genomes; however, it is possible that the arachnid clot consists of a related protein, spätzle, that is present in all of the genomes.     

昆虫天然免疫相关基因研究进展

在长期进化过程中,昆虫形成了强大的天然免疫防御系统,即体液免疫和细胞免疫。体液免疫主要包括Toll、IMD和JAK/STAT 3条信号通路,通过信号转导及免疫途径调控免疫相关基因的表达,诱导产生抗菌肽和其他效应分子。细胞免疫由血细胞介导,主要完成对病原物的包裹、吞噬和集结等。近年来,昆虫基因组学快速发展,通过生物信息学等方法从昆虫基因组数据中已鉴定到大量免疫相关基因,对这些基因的研究加深了人们对昆虫天然免疫分子机制的认识和理解。根据基因功能,免疫相关基因分为识别、信号转导、调制器、效应分子、黑化反应、RNA干扰和其他基因等7类,这些基因通过互作来调控体液免疫和细胞免疫。本文对昆虫免疫相关基因的分类、功能及家族进化等方面的研究成果进行总结,并对今后昆虫免疫的研究重点进行了展望,以期为昆虫免疫分子机制的研究及开发新的害虫防治策略提供依据。     

玫烟色棒束孢诱导的小菜蛾免疫响应表达谱分析

【目的】本研究旨在筛选小菜蛾Plutella xylostella应对玫烟色棒束孢Isaria fumosorosea侵染的免疫应答及其网络调控基因,以进一步探讨小菜蛾对玫烟色棒束孢的防御机制。【方法】采用第二代高通量测序技术RNA-seq,对处理后12 h的感染玫烟色棒束孢和健康(平行对照)的小菜蛾4龄幼虫转录组进行了测序,通过生物信息学分析对得到的差异基因进行了功能注释和分类,对差异基因参与的信号通路进行了分析。【结果】在感菌和健康小菜蛾幼虫转录组两个表达谱里,分别获得了12 346 987和12 315 210个clean reads,有60.93%和61.26%的reads分别能比对到参考基因库里,其中完美匹配(perfect match)的比例分别为32.15%和32.73%。共得到351个显著差异表达基因(differentially expressed unigenes,DEUs),上调表达基因有275个,下调表达基因有76个,与免疫防御反应潜在相关的基因有156个。GO(Gene Ontology)富集分析有102个DEUs分布到46个GO term里,KEGG(Kyoto Encyclopedia of Genes and Genomes)pathway富集分析结果显示,有132个DEUs显著富集在13个代谢通路(pathway)里。【结论】这些差异表达基因中,大部分编码潜在的与免疫识别及调控相关的基因,集中在能量代谢、疾病反应和防御反应等相关通路。研究结果为挖掘与玫烟色棒束孢侵染相关的小菜蛾免疫应答候选基因提供了重要数据库,也为阐明小菜蛾对玫烟色棒束孢的免疫机制奠定理论基础。     

Isolation and characterization of an immunosuppressive protein from venom of the pupa-specific endoparasitoid Pteromalus puparum

In hymenopteran parasitoids devoid of symbiotic viruses, venom proteins appear to play a major role in host immune suppression and host regulation. Not much is known about the active components of venom proteins in these parasitoids, especially those that have the functions involved in the suppression of host cellular immunity. Here, we report the isolation and characterization of a venom protein Vn.11 with 24.1 kDa in size from Pteromalus puparum, a pupa-specific endoparasitoid of Pieris rapae. The Vn.11 venom protein is isolated with the combination of ammonium sulfate precipitation and anion exchange chromatography, and its purity is verified using SDS-PAGE analysis. Like crude venom, the Vn.11 venom protein significantly inhibits the spreading behavior and encapsulation ability of host hemocytes in vitro. It is suggested that this protein is an actual component of P. puparum crude venom as host cellular-immune suppressive factor.     

Biogenic amine biosynthetic and transduction genes in the endoparasitoid wasp Pteromalus puparum (Hymenoptera: Pteromalidae)

Biogenic amines (BAs), such as octopamine, tyramine, dopamine, serotonin, and acetylcholine regulate various behaviors and physiological functions in insects. Here, we identified seven genes encoding BA biosynthetic enzymes and 16 genes encoding BA G protein-coupled receptors in the genome of the endoparasitoid wasp, Pteromalus puparum. We compared the genes with their orthologs in its host Pieris rapae and the related ectoparasitic wasp Nasonia vitripennis. All the genes show high (>90%) identity to orthologs in N. vitripennis. P. puparum and N. vitripennis have the smallest number of BA receptor genes among the insect species we investigated. We then analyzed the expression profiles of the genes, finding those acting in BA biosynthesis were highly expressed in adults and larvae and those encoding BA receptors are highly expressed in adults than immatures. Octα1R and 5-HT₇ genes were highly expressed in salivary glands, and a high messenger RNA level of 5-HT_1A was found in venom apparatuses. We infer that BA signaling is a fundamental component of the organismal organization, homeostasis and operation in parasitoids, some of the smallest insects.     

Fungal incompatibility: Evolutionary origin in pathogen defense?

In fungi, cell fusion between genetically unlike individuals triggers a cell death reaction known as the incompatibility reaction. In Podospora anserina, the genes controlling this process belong to a gene family encoding STAND proteins with an N‐terminal cell death effector domain, a central NACHT domain and a C‐terminal WD‐repeat domain. These incompatibility genes are extremely polymorphic, subject to positive Darwinian selection and display a remarkable genetic plasticity allowing for constant diversification of the WD‐repeat domain responsible for recognition of non‐self. Remarkably, the architecture of these proteins is related to pathogen‐recognition receptors ensuring innate immunity in plants and animals. Here, we hypothesize that these P. anserina incompatibility genes could be components of a yet‐unidentified innate immune system of fungi. As already proposed in the case of plant hybrid necrosis or graft rejection in mammals, incompatibility could be a by‐product of pathogen‐driven divergence in host defense genes.     

Role of the penetration‐resistance genes PEN1, PEN2 and PEN3 in the hypersensitive response and race‐specific resistance in Arabidopsis thaliana

Plants are highly capable of recognizing and defending themselves against invading microbes. Adapted plant pathogens secrete effector molecules to suppress the host's immune system. These molecules may be recognized by host‐encoded resistance proteins, which then trigger defense in the form of the hypersensitive response (HR) leading to programmed cell death of the host tissue at the infection site. The three proteins PEN1, PEN2 and PEN3 have been found to act as central components in cell wall‐based defense against the non‐adapted powdery mildew Blumeria graminis fsp. hordei (Bgh). We found that loss of function mutations in any of the three PEN genes cause decreased hypersensitive cell death triggered by recognition of effectors from oomycete and bacterial pathogens in Arabidopsis. There were considerable additive effects of the mutations. The HR induced by recognition of AvrRpm1 was almost completely abolished in the pen2 pen3 and pen1 pen3 double mutants and the loss of cell death could be linked to indole glucosinolate breakdown products. However, the loss of the HR in pen double mutants did not affect the plants' ability to restrict bacterial growth, whereas resistance to avirulent isolates of the oomycete Hyaloperonospora arabidopsidis was strongly compromised. In contrast, the double and triple mutants demonstrated varying degrees of run‐away cell death in response to Bgh. Taken together, our results indicate that the three genes PEN1, PEN2 and PEN3 extend in functionality beyond their previously recognized functions in cell wall‐based defense against non‐host pathogens.