Thitarodes armoricanus is a medicinal and economically important Lepidopteran insect species. The larvae infected by Paecilomyces hepiali survive no more than four days, while those infected by Ophiocordyceps sinensis can survive for several months before mummification. This provides a valuable comparative system to study interactions between an insect host and different pathogenic fungi. By using the T. armoricanus genome, a time-course transcriptome analysis of the whole larvae without guts was performed to explore the larvae response to P. hepiali and O. sinensis infection. A total of 3106 differentially expressed genes in five clusters were identified. The genes involved in coagulation and multiple metabolisms were both suppressed after P. hepiali or O. sinensis infection, whereas those related to environmental information responses, cell processes, biotic stimulus, and immunity (such as cecropin (CEC)) were elevated. The rapid death of T. armoricanus after P. hepiali infection might be caused by osmotic imbalance, immunocompromise (such as DEFs and GLVs), and nervous system dysfunction (glutamatergic synapse). Up-regulation of the genes related to cuticle structure, nervous system (such as neurotrophin signal pathway and dopaminergic synapse) and immune effectors (such as attacin (ATT) and proline-rich antimicrobial peptide 1 (PRAMP1)) in T. armoricanus, may contribute to the co-existence of T. armoricanus and O. sinensis. This study provides a global view and potential key genes of the interaction between T. armoricanus and two fungal entomopathogens.

Aggressive behavior in animals is important for survival and reproduction. It is well studied in adult insects, such as flies, ants, honey bees, and crickets. However, the larvae of Lepidopteran insects are also aggressive, studies of which are still lacking. Here, RNA-seq was used to generate a high-quality database for the aggressive behavior of Thitarodes xiaojinensis toward conspecifics and heterospecifics. Although there was similar aggressive behavior between the conspecific group and heterospecific group, significant differences were identified at the transcriptional level. When there was aggressive behavior toward conspecifics, T. xiaojinensis trended toward higher expression at the respiratory chain, while cuticle development and metabolism may have interfered. On the other hand, when there was aggressive behavior toward H. armigera, genes related to neuron and cuticle development, cellular processes, and its regulated signaling pathways were significantly upregulated, while the genes associated with oxidation-reduction and metabolism were downregulated. Weighted gene co-expression networks analysis (WGCNA) was performed, and two modules with properties correlating to the aggressive behavior of T. xiaojinensis were identified. Several hub genes were predicted and confirmed by qRT-PCR, such as CLTC, MYH, IGF2BP1, and EMC. This study provides a global view and potential key genes for the aggressive behavior of T. xiaojinensis toward conspecifics and heterospecifics. Further investigation of the hub genes would help us to better understand the aggressive behavior of insects.

The Chinese cordyceps, a parasitic Ophiocordyceps sinensis fungus–Thitarodes/Hepialus larva complex, is a valuable biological resource endemic to the Tibetan Plateau. Protection of the Plateau environment and huge market demand make it necessary to culture this complex in an artificial system. A method for the large-scale artificial rearing of the Thitarodes/Hepialus insect host has been established. However, the deterioration of the insect rearing population and low mummification of the infected larvae by the fungus constrain effective commercial cultivation. Hybridization of Thitarodes/Hepialus populations may be needed to overcome this problem. The species T. shambalaensis (GG♂ × GG♀) and an undescribed Thitarodes species (SD♂ × SD♀) were inbred or hybridized to evaluate the biological parameters, larval sensitivity to the fungal infection and mitochondrial genomes of the resulting populations. The two parental Thitarodes species exhibited significant differences in adult fresh weights and body lengths but not in pupal emergence rates. Hybridization of T. shambalaensis and Thitarodes sp. allowed producing a new generation. The SD♂ × GG♀ population showed a higher population trend index than the SD♂ × SD♀ population, implying increased population growth compared with the male parent. The sensitivity of the inbred larval populations to four fungal isolates of O. sinensis also differed. This provides possibilities to create Thitarodes/Hepialus populations with increased growth potential for the improved artificial production of the insect hosts. The mitochondrial genomes of GG♂ × GG♀, SD♂ × SD♀ and SD♂ × GG♀ were 15,612 bp, 15,389 bp and 15,496 bp in length, with an A + T content of 80.92%, 82.35% and 80.87%, respectively. The A + T-rich region contains 787 bp with two 114 bp repetitive sequences, 554 bp without repetitive sequences and 673 bp without repetitive sequences in GG♂ × GG♀, SD♂ × SD♀ and SD♂ × GG♀, respectively. The hybrid population (SD♂ × GG♀) was located in the same clade with GG♂ × GG♀, based on the phylogenetic tree constructed by 13 PCGs, implying the maternal inheritance of mitochondrial DNA.