昆虫学报 ›› 2024, Vol. 67 ›› Issue (3): 431-442.doi: 10.16380/j.kcxb.2024.03.013

• 简 报 • 上一篇    下一篇

高通量测序分析褐飞虱circRNA表达谱及其潜在功能

于鑫滢1,3, 赵哲奇1,3, 闸雯俊1,2,4,*   

  1. (1. 湖北省农业科学院粮食作物研究所, 粮食作物种质创新与遗传改良湖北省重点实验室, 武汉 430064; 2. 湖北洪山实验室, 武汉 430070; 3. 武汉大学生命科学学院, 武汉 430072; 4. 农业农村部作物分子育种重点实验室, 武汉 430070)
  • 出版日期:2024-03-20 发布日期:2024-04-17

Expression profiles and potential functions of circRNAs in Nilaparvata lugens (Hemiptera: Delphacidae) analyzed by high-throughput sequencing

YU Xin-Ying1,3, ZHAO Zhe-Qi1,3, ZHA Wen-Jun1,2,4,*   

  1.  (1. Hubei Key Laboratory of Food Crop Germplasms and Genetic Improvement, Institute of Food Crops, Hubei Academy of Agricultural Sciences, Wuhan 430064, China; 2. Hubei Hongshan Laboratory, Wuhan 430070, China; 3. College of Life Sciences, Wuhan University, Wuhan 430072, China; 4. Key Laboratory of Crop Molecular Breeding, Ministry of Agriculture and Rural Affairs, Wuhan 430072, China)
  • Online:2024-03-20 Published:2024-04-17

摘要: 【目的】探究褐飞虱Nilaparvata lugens如何适应抗性水稻品种YHY15,并为研究环状RNA(circular RNA, circRNA)在褐飞虱对抗性水稻适应机制中的作用奠定基础。【方法】利用高通量测序技术鉴定生物型1和生物型Y(能致害抗性水稻YHY15的生物型)褐飞虱的circRNA,并统计其类型和分布。分析生物型1和生物型Y褐飞虱间circRNA表达差异,并对circRNA来源基因进行GO功能注释和KEGG富集分析。【结果】生物型1和生物型Y褐飞虱共鉴定到19个circRNA,分布在9条染色体上,其中17个circRNA由剪接位点之间的外显子构成,2个circRNA由剪接位点之间的所有碱基构成,17个circRNA长度在200~800 bp。生物型Y褐飞虱的circRNA表达丰度与表达量比生物型1褐飞虱的更高。分析结果表明,在生物型1褐飞虱克服抗性水稻YHY15形成新的生物型Y的过程中,在自然选择压力下,circRNA的表达出现变化,推测这种变化影响了褐飞虱新生物型的形成。KEGG分析结果表明,鉴定到的生物型1和生物型Y褐飞虱circRNA来源基因主要富集到自噬相关通路上。【结论】褐飞虱长期取食抗性水稻,导致褐飞虱消化、解毒和代谢能力进化,能降解抗性水稻为抵御褐飞虱而产生的次生代谢物。鉴定到的生物型1和生物型Y褐飞虱的circRNA来源基因与细胞自噬相关,表明褐飞虱通过细胞自噬过程应答抗虫水稻,这种应答反应促进褐飞虱适应抗虫水稻,致害性增强和生物型形成。

关键词:  褐飞虱, circRNA, 细胞自噬, 抗虫水稻, 高通量测序

Abstract: 【Aim】 To explore how the brown planthopper (Nilaparvata lugens) adapts to the resistant rice variety YHY15 and lay a foundation for studying the role of circular RNA (circRNA) in the adaptation mechanism of N. lugens to resistant rice. 【Methods】 The circRNAs of biotype 1 and biotype Y (biotype that can cause damage to resistant rice YHY15) of N. lugens were identified by high-throughput sequencing technology, and the type and distribution of circRNAs were counted. The expression differences of circRNAs between biotype 1 and biotype Y of N. lugens were analyzed, and GO functional annotation and KEGG enrichment analysis of the circRNA source genes were performed. 【Results】 A total of 19 circRNAs were identified in biotype 1 and biotype Y of N. lugens, and distributed on nine chromosomes. Among them, 17 circRNAs were composed of exons between splicing sites, two circRNAs were composed of all bases between splicing sites, and 17 circRNAs were 200-800 bp in length. The expression abundance and expression level of circRNAs in biotype Y of N. lugens were higher than those in biotype 1 of N. lugens. The analysis results showed that in the process of biotype 1 of N. lugens overcoming the resistance of rice YHY15 to form a new biotype Y of N. lugens, under the pressure of natural selection, the expression of circRNA changed resulting in the effects on the formation of a new biotype of N. lugens. The KEGG analysis results indicated that the identified circRNA source genes of biotype 1 and biotype Y of N. lugens were mainly enriched in autophagy-related pathways. 【Conclusion】 Long-term feeding of resistant rice by N. lugens leads to the evolution of the digestion, detoxification, and metabolic abilities of N. lugens, which can degrade the secondary metabolites produced by resistant rice to resist N. lugens. The circRNA source genes of biotype 1 and biotype Y of N. lugens are enriched in autophagy-related pathways, indicating that N. lugens respond to insect-resistant rice through the process of autophagy. This response promotes the adaptation of N. lugens to insect-resistant rice, enhances pathogenicity, and forms biotypes.

Key words: Nilaparvata lugens; circRNA, autophagy, insect-resistant rice, high-throughput sequencing