Profiling Diuraphis noxia (Hemiptera: Aphididae) Transcript Expression of the Biotypes SA1 and SAM Feeding on Various Triticum aestivum Varieties

Whole-body transcriptome mining for candidate effectors from Diuraphis noxia

Background

Proteins within aphid saliva play a crucial role as the molecular interface between aphids and their host plants. These salivary effectors modulate plant responses to favour aphid feeding and facilitate infestation. The identification of effectors from economically important pest species is central in understanding the molecular events during the aphid-plant interaction. The Russian wheat aphid (Diuraphis noxia, Kurdjumov) is one such pest that causes devastating losses to wheat and barley yields worldwide. Despite the severe threat to food security posed by D. noxia, the non-model nature of this pest and its host has hindered progress towards understanding this interaction. In this study, in the absence of a salivary gland transcriptome, whole-body transcriptomics data was mined to generate a candidate effector catalogue for D. noxia.

Results

Mining the transcriptome identified 725 transcripts encoding putatively secreted proteins amongst which were transcripts specific to D. noxia. Six of the seven examined D. noxia putative effectors, termed DnE’s (Diuraphis noxia effectors) exhibited salivary gland-specific expression. A comparative analysis between whole-body D. noxia transcriptome data versus the head and body transcriptomes from three other aphid species allowed us to define a catalogue of transcripts putatively upregulated in D. noxia head tissue. Five of these were selected for RT-qPCR confirmation, and were found to corroborate the differential expression predictions, with a further three confirmed to be highly expressed in D. noxia salivary gland tissue.

Conclusions

Determining a putative effector catalogue for D. noxia from whole-transcriptome data, particularly the identification of salivary-specific sequences potentially unique to D. noxia, provide the basis for future functional characterisation studies to gain further insight into this aphid-plant interaction. Furthermore, due to a lack of publicly available aphid salivary gland transcriptome data, the capacity to use comparative transcriptomics to compile a list of putative effector candidates from whole-body transcriptomics data will further the study of effectors in various aphid species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08712-4.

DNA Methylation and Demethylation Are Regulated by Functional DNA Methyltransferases and DnTET Enzymes in Diuraphis noxia

Aphids are economically important insect pests of crops worldwide. Despite resistant varieties being available, resistance is continuously challenged and eventually broken down, posing a threat to food security. In the current study, the epigenome of two related Russian wheat aphid (Diuraphis noxia, Kurdjumov) biotypes (i.e., SA1 and SAM) that differ in virulence was investigated to elucidate its role in virulence in this species. Whole genome bisulfite sequencing covered a total of 6,846,597,083 cytosine bases for SA1 and 7,397,965,699 cytosine bases for SAM, respectively, of which a total of 70,861,462 bases (SA1) and 74, 073,939 bases (SAM) were methylated, representing 1.126 ± 0.321% (SA1) and 1.105 ± 0.295% (SAM) methylation in their genomes. The sequence reads were analyzed for contexts of DNA methylation and the results revealed that RWA has methylation in all contexts (CpG, CHG and CHH), with the majority of methylation within the CpG context (± 5.19%), while the other contexts show much lower levels of methylation (CHG - ± 0.27%; CHH - ± 0.34%). The top strand was slightly (0.02%) more methylated than the bottom strand. Of the 35,493 genes that mapped, we also analyzed the contexts of methylation of each of these and found that the CpG methylation was much higher in genic regions than in intergenic regions. The CHG and CHH levels did not differ between genic and intergenic regions. The exonic regions of genes were more methylated (±0.56%) than the intronic regions. We also measured the 5mC and 5hmC levels between the aphid biotypes, and found little difference in 5mC levels between the biotypes, but much higher levels of 5hmC in the virulent SAM. RWA had two homologs of each of the DNA methyltransferases 1 (DNMT1a and DNMT1b) and DNMT3s (DNMT3a and DNMT3b), but only a single DNMT2, with only the expression of DNMT3 that differed significantly between the two RWA biotypes. RWA has a single ortholog of Ten eleven translocase (DnTET) in the genome. Feeding studies show that the more virulent RWA biotype SAM upregulate DnDNMT3 and DnTET in response to wheat expressing antibiosis and antixenosis.

Gene Expression and Diet Breadth in Plant-Feeding Insects: Summarizing Trends

Transcriptomic studies lend insights into the role of transcriptional plasticity in adaptation and specialization. Recently, there has been growing interest in understanding the relationship between variation in herbivorous insect gene expression and the evolution of diet breadth. We review the studies that have emerged on insect gene expression and host plant use, and outline the questions and approaches in the field. Many candidate genes underlying herbivory and specialization have been identified, and a few key studies demonstrate increased transcriptional plasticity associated with generalist compared with specialist species. Addressing the roles that transcriptional variation plays in insect diet breadth will have important implications for our understanding of the evolution of specialization and the genetic and environmental factors that govern insect-plant interactions.