Long noncoding RNAs are potentially involved in the degeneration of virulence in an aphid-obligate pathogen, (Entomophthoromycotina).

Virulence attenuation frequently occurs in culturing of pathogenic microbes. In this study, we investigated the total putative long noncoding RNAs (lncRNAs) in an aphid-obligate pathogen, , and screened the differentially expressed (DE) lncRNAs and protein-coding genes involved in the virulence decline. The virulence was significantly attenuated after eight subculturing events, in which the median lethal concentration of the conidia ejected from mycelial mats relative to the bamboo aphid, , increased from 36.1 to 126.1 conidia mm, four days after inoculation. In total, 1,252 lncRNAs were identified based on the genome-wide transcriptional analysis. By characterizing their molecular structures and expression patterns, we found that the lncRNAs possessed shorter transcripts, lower expression, and fewer exons than did protein-coding genes in . A total of 410 DE genes of 329 protein-coding genes and 81 lncRNAs were identified. The functional enrichment analysis showed the DE genes were enriched in peptidase activity, protein folding, autophagy, and metabolism. Moreover, target prediction analysis of the 81 lncRNAs revealed 3,111 -regulated and 23 -regulated mRNAs, while 121 DE lncRNA-mRNA pairs were possibly involved in virulence decline. Moreover, the DE lncRNA-regulated target genes mainly encoded small heat shock proteins, secretory proteins, transporters, autophagy proteins, and other stress response-related proteins. This implies that the decline in virulence regulated by lncRNAs was likely associated with the environmental stress response of . Hence, these findings can provide insights into the lncRNA molecules of Entomophthoromycotina, with regards to virulence regulators of entomopathogens.