Differential subgenome expression underlies biomass accumulation in allotetraploid Pennisetum giganteum.

BACKGROUND: Pennisetum giganteum (AABB, 2n = 4x = 28) is a C4 plant in the genus Pennisetum with origin in Africa but currently also grown in Asia and America. It is a crucial forage and potential energy grass with significant advantages in yield, stress resistance, and environmental adaptation. However, the mechanisms underlying these advantageous traits remain largely unexplored. Here, we present a high-quality genome assembly of the allotetraploid P. giganteum aiming at providing insights into biomass accumulation. RESULTS: Our assembly has a genome size 2.03 Gb and contig N50 of 88.47 Mb that was further divided into A and B subgenomes. Genome evolution analysis revealed the evolutionary relationships across the Panicoideae subfamily lineages and identified numerous genome rearrangements that had occurred in P. giganteum. Comparative genomic analysis showed functional differentiation between the subgenomes. Transcriptome analysis found no subgenome dominance at the overall gene expression level; however, differentially expressed homoeologous genes and homoeolog-specific expressed genes between the two subgenomes were identified, suggesting that complementary effects between the A and B subgenomes contributed to biomass accumulation of P. giganteum. Besides, C4 photosynthesis-related genes were significantly expanded in P. giganteum and their sequences and expression patterns were highly conserved between the two subgenomes, implying that both subgenomes contributed greatly and almost equally to the highly efficient C4 photosynthesis in P. giganteum. We also identified key candidate genes in the C4 photosynthesis pathway that showed sustained high expression across all developmental stages of P. giganteum. CONCLUSIONS: Our study provides important genomic resources for elucidating the genetic basis of advantageous traits in polyploid species, and facilitates further functional genomics research and genetic improvement of P. giganteum.

A chromosome-level genome assembly of the beet armyworm Spodoptera exigua.

BACKGROUND: The beet armyworm Spodoptera exigua is a polyphagous caterpillar that causes serious damage to many species of crops and vegetables. To gain insight into how this polyphagous insect differs from less harmful oligophagous species, we generated a chromosome-level assembly and compared it to closely related species with the same or different feeding habits. RESULTS: Based on Illumina and Pacific Biosciences data and Hi-C technology, 425.6 Mb of genome sequences were anchored and oriented into 31 linkage groups, with an N50 length of 14.8 Mb. A total of 24,649 gene models were predicted, of which 97.4% were identified in the genome assembly. Chemosensory genes are vital for locating food: of the four main families, odorant-binding proteins, chemosensory proteins and olfactory receptors showed little difference, whereas gustatory receptors are greatly expanded in S. exigua. Examination of other polyphagous insects confirmed this difference from oligophagous congeners and further identified the bitter receptor subfamily as being particularly affected. CONCLUSION: Our high-quality genome sequence for beet armyworm identified a key expansion of the bitter gustatory receptor subfamily in this and other pests that differs crucially from more benign relatives and offers insight into the biology and possible future means of control for these economically important insects.

Chromosome-level genome assembly and population genomic analyses provide insights into adaptive evolution of the red turpentine beetle, Dendroctonus valens.

BACKGROUND: Biological invasions are responsible for substantial environmental and economic losses. The red turpentine beetle (RTB), Dendroctonus valens LeConte, is an important invasive bark beetle from North America that has caused substantial tree mortality in China. The lack of a high-quality reference genome seriously limits deciphering the extent to which genetic adaptions resulted in a secondary pest becoming so destructive in its invaded area. RESULTS: Here, we present a 322.41 Mb chromosome-scale reference genome of RTB, of which 98% of assembled sequences are anchored onto fourteen linkage groups including the X chromosome with a N50 size of 24.36 Mb, which is significantly greater than other Coleoptera species. Repetitive sequences make up 45.22% of the genome, which is higher than four other Coleoptera species, i.e., Mountain pine beetle Dendroctonus ponderosae, red flour beetle Tribolium castaneum, blister beetle Hycleus cichorii, and Colorado potato beetle Leptinotarsa decemlineata. We identify rapidly expanded gene families and positively selected genes in RTB, which may be responsible for its rapid environmental adaptation. Population genetic structure of RTB was revealed by genome resequencing of geographic populations in native and invaded regions, suggesting substantial divergence of the North American population and illustrates the possible invasion and spread route in China. Selective sweep analysis highlighted the enhanced ability of Chinese populations in environmental adaptation. CONCLUSIONS: Overall, our high-quality reference genome represents an important resource for genomics study of invasive bark beetles, which will facilitate the functional study and decipher mechanism underlying invasion success of RTB by integrating the Pinus tabuliformis genome.

Full-length codling moth transcriptome atlas revealed by single-molecule real-time sequencing.

Over the past decade, second-generation sequencing (SGS) has been widely used to elucidate the transcriptome across many organisms. However, the full-length (FL) transcripts and alternative splice (AS) isoforms could not be confidently and accurately defined with SGS. Pacific biosciences (PacBio) single-molecule real-time sequencing was conducted to obtain FL transcriptome data in the codling moth. In total, 25,940 high-quality FL isoforms were obtained and clustered to 14,099 nonredundant clusters. Interestingly, nearly 90% of nonredundant PacBio transcripts were novel compared to reference genes. Among them, 3389 transcripts potentially represented novel genes. Additionally, a large number of AS events were discovered, and most of the splice junctions in the PacBio isoforms could be supported by short reads in public datasets. Furthermore, 952 FL lncRNAs and 81 fusion transcripts were identified and validated using RT-PCR analysis. Overall, an atlas of FL transcripts was obtained in the codling moth, which will help provide further insights into the complexity of the transcriptome and facilitate improving genome annotations and functional studies in this insect.

A Genome-Wide Analysis of Serine Protease Inhibitors in Cydia pomonella Provides Insights into Their Evolution and Expression Pattern.

Serine protease inhibitors (serpins) appear to be ubiquitous in almost all living organisms, with a conserved structure and varying functions. Serpins can modulate immune responses by negatively regulating serine protease activities strictly and precisely. The codling moth, Cydia pomonella (L.), a major invasive pest in China, can cause serious economic losses. However, knowledge of serpin genes in this insect remain largely unknown. In this study, we performed a systematic analysis of the serpin genes in C. pomonella, obtaining 26 serpins from the C. pomonella genome. Subsequently, their sequence features, evolutionary relationship, and expression pattern were characterized. Comparative analysis revealed the evolution of a number of serpin genes in Lepidoptera. Importantly, the evolutionary relationship and putative roles of serpin genes in C. pomonella were revealed. Additionally, selective pressure analysis found amino acid sites with strong evidence of positive selection. Interestingly, the serpin1 gene possessed at least six splicing isoforms with distinct reactive-center loops, and these isoforms were experimentally validated. Furthermore, we observed a subclade expansion of serpins, and these genes showed high expression in multiple tissues, suggesting their important roles in C. pomonella. Overall, this study will enrich our knowledge of the immunity of C. pomonella and help to elucidate the role of serpins in the immune response.

High-quality Fagopyrum esculentum genome provides insights into the flavonoid accumulation among different tissues and self-incompatibility.

Common buckwheat (Fagopyrum esculentum) and Tartary buckwheat (Fagopyrum tataricum), the two most widely cultivated buckwheat species, differ greatly in flavonoid content and reproductive mode. Here, we report the first high-quality and chromosome-level genome assembly of common buckwheat with 1.2 Gb. Comparative genomic analysis revealed that common buckwheat underwent a burst of long terminal repeat retrotransposons insertion accompanied by numerous large chromosome rearrangements after divergence from Tartary buckwheat. Moreover, multiple gene families involved in stress tolerance and flavonoid biosynthesis such as multidrug and toxic compound extrusion (MATE) and chalcone synthase (CHS) underwent significant expansion in buckwheat, especially in common buckwheat. Integrated multi-omics analysis identified high expression of catechin biosynthesis-related genes in flower and seed in common buckwheat and high expression of rutin biosynthesis-related genes in seed in Tartary buckwheat as being important for the differences in flavonoid type and content between these buckwheat species. We also identified a candidate key rutin-degrading enzyme gene (Ft8.2377) that was highly expressed in Tartary buckwheat seed. In addition, we identified a haplotype-resolved candidate locus containing many genes reportedly associated with the development of flower and pollen, which was potentially related to self-incompatibility in common buckwheat. Our study provides important resources facilitating future functional genomics-related research of flavonoid biosynthesis and self-incompatibility in buckwheat.

The landscape of lncRNAs in Cydia pomonella provides insights into their signatures and potential roles in transcriptional regulation.

BACKGROUND: Long noncoding RNAs (lncRNAs) have emerged as an important class of transcriptional regulators in cellular processes. The past decades have witnessed great progress in lncRNA studies in a variety of organisms. The codling moth (Cydia pomonella L.) is an important invasive insect in China. However, the functional impact of lncRNAs in this insect remains unclear. In this study, an atlas of codling moth lncRNAs was constructed based on publicly available RNA-seq datasets. RESULTS: In total, 9875 lncRNA transcripts encoded by 9161 loci were identified in the codling moth. As expected, the lncRNAs exhibited shorter transcript lengths, lower GC contents, and lower expression levels than protein-coding genes (PCGs). Additionally, the lncRNAs were more likely to show tissue-specific expression patterns than PCGs. Interestingly, a substantial fraction of the lncRNAs showed a testis-biased expression pattern. Additionally, conservation analysis indicated that lncRNA sequences were weakly conserved across insect species, though additional lncRNAs with homologous relationships could be identified based on synteny, suggesting that synteny could be a more reliable approach for the cross-species comparison of lncRNAs. Furthermore, the correlation analysis of lncRNAs with neighbouring PCGs indicated a stronger correlation between them, suggesting potential cis-acting roles of these lncRNAs in the regulation of gene expression. CONCLUSIONS: Taken together, our work provides a valuable resource for the comparative and functional study of lncRNAs, which will facilitate the understanding of their mechanistic roles in transcriptional regulation.

Inhibition of melanization by serpin-5 and serpin-9 promotes baculovirus infection in cotton bollworm Helicoverpa armigera.

Melanization, an important insect defense mechanism, is mediated by clip-domain serine protease (cSP) cascades and is regulated by serpins. Here we show that proteolytic activation of prophenoloxidase (PPO) and PO-catalyzed melanization kill the baculovirus in vitro. Our quantitative proteomics and biochemical experiments revealed that baculovirus infection of the cotton bollworm, Helicoverpa armigera, reduced levels of most cascade members in the host hemolymph and PO activity. By contrast, serpin-9 and serpin-5 were sequentially upregulated after the viral infection. The H. armigera serpin-5 and serpin-9 regulate melanization by directly inhibiting their target proteases cSP4 and cSP6, respectively and cSP6 activates PPO purified from hemolymph. Furthermore, serpin-5/9-depleted insects exhibited high PO activities and showed resistance to baculovirus infection. Together, our results characterize a part of the melanization cascade in H. armigera, and suggest that natural insect virus baculovirus has evolved a distinct strategy to suppress the host immune system.

Dynamics of the Interaction between Cotton Bollworm Helicoverpa armigera and Nucleopolyhedrovirus as Revealed by Integrated Transcriptomic and Proteomic Analyses.

Over the past decades, Helicoverpa armigera nucleopolyhedrovirus (HearNPV) has been widely used for biocontrol of cotton bollworm, which is one of the most destructive pest insects in agriculture worldwide. However, the molecular mechanism underlying the interaction between HearNPV and host insects remains poorly understood. In this study, high-throughput RNA-sequencing was integrated with label-free quantitative proteomics analysis to examine the dynamics of gene expression in the fat body of H. armigera larvae in response to challenge with HearNPV. RNA sequencing-based transcriptomic analysis indicated that host gene expression was substantially altered, yielding 3,850 differentially expressed genes (DEGs), whereas no global transcriptional shut-off effects were observed in the fat body. Among the DEGs, 60 immunity-related genes were down-regulated after baculovirus infection, a finding that was consistent with the results of quantitative real-time RT-PCR. Gene ontology and functional classification demonstrated that the majority of down-regulated genes were enriched in gene cohorts involved in energy, carbohydrate, and amino acid metabolic pathways. Proteomics analysis identified differentially expressed proteins in the fat body, among which 76 were up-regulated, whereas 373 were significantly down-regulated upon infection. The down-regulated proteins are involved in metabolic pathways such as energy metabolism, carbohydrate metabolism (CM), and amino acid metabolism, in agreement with the RNA-sequence data. Furthermore, correlation analysis suggested a strong association between the mRNA level and protein abundance in the H. armigera fat body. More importantly, the predicted gene interaction network indicated that a large subset of metabolic networks was significantly negatively regulated by viral infection, including CM-related enzymes such as aldolase, enolase, malate dehydrogenase, and triose-phosphate isomerase. Taken together, transcriptomic data combined with proteomic data elucidated that baculovirus established systemic infection of host larvae and manipulated the host mainly by suppressing the host immune response and down-regulating metabolism to allow viral self-replication and proliferation. Therefore, this study provided important insights into the mechanism of host-baculovirus interaction.

A Critical Role for CLSP2 in the Modulation of Antifungal Immune Response in Mosquitoes.

Entomopathogenic fungi represent a promising class of bio-insecticides for mosquito control. Thus, detailed knowledge of the molecular mechanisms governing anti-fungal immune response in mosquitoes is essential. In this study, we show that CLSP2 is a modulator of immune responses during anti-fungal infection in the mosquito Aedes aegypti. With a fungal infection, the expression of the CLSP2 gene is elevated. CLSP2 is cleaved upon challenge with Beauveria bassiana conidia, and the liberated CLSP2 CTL-type domain binds to fungal cell components and B. bassiana conidia. Furthermore, CLPS2 RNA interference silencing significantly increases the resistance to the fungal challenge. RNA-sequencing transcriptome analysis showed that the majority of immune genes were highly upregulated in the CLSP2-depleted mosquitoes infected with the fungus. The up-regulated immune gene cohorts belong to melanization and Toll pathways, but not to the IMD or JAK-STAT. A thioester-containing protein (TEP22), a member of α2-macroglobulin family, has been implicated in the CLSP2-modulated mosquito antifungal defense. Our study has contributed to a greater understanding of immune-modulating mechanisms in mosquitoes.

High throughput profiling of the cotton bollworm Helicoverpa armigera immunotranscriptome during the fungal and bacterial infections.

BACKGROUND: Innate immunity is essential in defending against invading pathogens in invertebrates. The cotton bollworm, Helicoverpa armigera (Hübner) is one of the most destructive lepidopteran pests, which causes enormous economic losses in agricultural production worldwide. The components of the immune system are largely unknown in this insect. The application of entomopathogens is considered as an alternative to the chemical insecticides for its control. However, few studies have focused on the molecular mechanisms of host-pathogen interactions between pest insects and their pathogens. Here, we investigated the immunotranscriptome of H. armigera larvae and examined gene expression changes after pathogen infections. This study provided insights into the potential immunity-related genes and pathways in H. armigera larvae. RESULTS: Here, we adopted a high throughput RNA-seq approach to determine the immunotranscriptome of H. armigera larvae injected with buffer, fungal pathogen Beauveria bassiana, or Gram-negative bacterium Enterobacter cloacae. Based on sequence similarity to those homologs known to participate in immune responses in other insects, we identified immunity-related genes encoding pattern recognition receptors, signal modulators, immune effectors, and nearly all members of the Toll, IMD and JAK/STAT pathways. The RNA-seq data indicated that some immunity-related genes were activated in fungus- and bacterium-challenged fat body while others were suppressed in B. bassiana challenged hemocytes, including the putative IMD and JAK-STAT pathway members. Bacterial infection elevated the expression of recognition and modulator genes in the fat body and signal pathway genes in hemocytes. Although fat body and hemocytes both are important organs involved in the immune response, our transcriptome analysis revealed that more immunity-related genes were induced in the fat body than that hemocytes. Furthermore, quantitative real-time PCR analysis confirmed that, consistent with the RNA-seq data, the transcript abundances of putative PGRP-SA1, Serpin1, Toll-14, and Spz2 genes were elevated in fat body upon B. bassiana infection, while the mRNA levels of defensin, moricin1, and gloverin1 were up-regulated in hemocytes. CONCLUSIONS: In this study, a global survey of the host defense against fungal and bacterial infection was performed on the non-model lepidopteran pest species. The comprehensive sequence resource and expression profiles of the immunity-related genes in H. armigera are acquired. This study provided valuable information for future functional investigations as well as development of specific and effective agents to control this pest.

A near-complete genome assembly of Monochamus alternatus a major vector beetle of pinewood nematode.

The Japanese sawyer beetle, Monochamus alternatus, is not only one of the most important wood boring pest itself, but also a major vector of the invasive pinewood nematode (PWN), which is the causal agent of the devastative pine wilt disease (PWD) and threats the global pine forest. Here, we present a near-complete genome of M. alternatus at the chromosome level. The assembled genome was 792.05 Mb with contig N50 length of 55.99 Mb, which is the largest N50 size among the sequenced Coleoptera insects currently. 99.57% of sequence was anchored onto ten pseudochromosomes (one X-chromosome and nine autosomes), and the final genome harbored only 13 gaps. BUSCO evaluation revealed the presence of 99.0% of complete core genes. Thus, our genome assembly represented the highest-contiguity genome assembly as well as high completeness in insects so far. We identified 20,471 protein-coding genes, of which 20,070 (98.04%) were functionally annotated. The genome assembly of M. alternatus provides a valuable resource for exploring the evolution of the symbiosis between PWN and the vector insects.

The haplotype-resolved genome assembly of autotetraploid rhubarb Rheum officinale provides insights into its genome evolution and massive accumulation of anthraquinones.

Rheum officinale, a member of the Polygonaceae family, is an important medicinal plant that is widely used in traditional Chinese medicine. Here, we report a 7.68-Gb chromosome-scale assembly of R. officinale with a contig N50 of 3.47 Mb, which was clustered into 44 chromosomes across four homologous groups. Comparative genomics analysis revealed that transposable elements have made a significant contribution to its genome evolution, gene copy number variation, and gene regulation and expression, particularly of genes involved in metabolite biosynthesis, stress resistance, and root development. We placed the recent autotetraploidization of R. officinale at ∼0.58 mya and analyzed the genomic features of its homologous chromosomes. Although no dominant monoploid genomes were observed at the overall expression level, numerous allele-differentially-expressed genes were identified, mainly with different transposable element insertions in their regulatory regions, suggesting that they functionally diverged after polyploidization. Combining genomics, transcriptomics, and metabolomics, we explored the contributions of gene family amplification and tetraploidization to the abundant anthraquinone production of R. officinale, as well as gene expression patterns and differences in anthraquinone content among tissues. Our report offers unprecedented genomic resources for fundamental research on the autopolyploid herb R. officinale and guidance for polyploid breeding of herbs.

The genome of okra (Abelmoschus esculentus) provides insights into its genome evolution and high nutrient content.

Okra (Abelmoschus esculentus) is an important vegetable crop with high nutritional value. However, the mechanism underlying its high nutrient content remains poorly understood. Here, we present a chromosome-scale genome of okra with a size of 1.19 Gb. Comparative genomics analysis revealed the phylogenetic status of A. esculentus, as well as whole-genome duplication (WGD) events that have occurred widely across the Malvaceae species. We found that okra has experienced three additional WGDs compared with the diploid cotton Gossypium raimondii, resulting in a large chromosome number (2n = 130). After three WGDs, okra has undergone extensive genomic deletions and retained substantial numbers of genes related to secondary metabolite biosynthesis and environmental adaptation, resulting in significant differences between okra and G. raimondii in the gene families related to cellulose synthesis. Combining transcriptomic and metabolomic analysis, we revealed the relationship between gene expression and metabolite content change across different okra developmental stages. Furthermore, the sinapic acid/S-lignin biosynthesis-related gene families have experienced remarkable expansion in okra, and the expression of key enzymes involved in the sinapic acid/S-lignin biosynthesis pathway vary greatly across developmental periods, which partially explains the differences in metabolite content across the different stages. Our study gains insights into the comprehensive evolutionary history of Malvaceae species and the genetic basis that underlies the nutrient content changes in okra, which will facilitate the functional study and genetic improvement of okra varieties.

Antiviral function of peptidoglycan recognition protein in Spodoptera exigua (Lepidoptera: Noctuidae).

Peptidoglycan recognition proteins (PGRPs) are a class of molecules that play a critical role in insect immunity. Understanding the function of PGRPs is important to improve the efficiency of microbial insecticides. In this study, we investigated the role of PGRP-LB (a long type PGRP) in insect immunity against viruses using Spodoptera exigua and Spodoptera exigua multiple nucleopolyhedrovirus (SeMNPV) as an insect-virus model. We cloned and identified a PGRP-LB gene from S. exigua; the gene consisted of 7 exons that encoded a polypeptide of 234 amino acids with a signal peptide and a typical amidase domain. Expression analysis revealed that the abundance of SePGRP-LB transcripts in the fat body was greater than in other tissues. Overexpression of SePGRP-LB resulted in a significant decrease of 49% in the rate of SeMNPV-infected cells. In addition, the multiplication of SeMNPV was significantly decreased: a decrease of 79% in the production of occlusion-derived virion (ODV), and a maximum decrease of 50% in the production of budded virion (BV). In contrast, silencing of SePGRP-LB expression by RNA interference resulted in a significant 1.65-fold increase in the rate of SeMNPV-infected cells, a significant 0.54-fold increase in ODV production, a maximum 1.57-fold increase in BV production, and the larval survival dropped to 21%. Our findings show that SePGRP-LB has an antiviral function against SeMNPV, and therefore this gene may provide a target for lepidopteran pest control using virus insecticides.

Plant pan-genomics: recent advances, new challenges, and roads ahead.

Pan-genomics can encompass most of the genetic diversity of a species or population and has proved to be a powerful tool for studying genomic evolution and the origin and domestication of species, and for providing information for plant improvement. Plant genomics has greatly progressed because of improvements in sequencing technologies and the rapid reduction of sequencing costs. Nevertheless, pan-genomics still presents many challenges, including computationally intensive assembly methods, high costs with large numbers of samples, ineffective integration of big data, and difficulty in applying it to downstream multi-omics analysis and breeding research. In this review, we summarize the definition and recent achievements of plant pan-genomics, computational technologies used for pan-genome construction, and the applications of pan-genomes in plant genomics and molecular breeding. We also discuss challenges and perspectives for future pan-genomics studies and provide a detailed pipeline for sample selection, genome assembly and annotation, structural variation identification, and construction and application of graph-based pan-genomes. The aim is to provide important guidance for plant pan-genome research and a better understanding of the genetic basis of genome evolution, crop domestication, and phenotypic diversity for future studies.

Identification of Pinelliae Rhizoma and its counterfeit species based on enzymatic signature peptides from toxic proteins.

BACKGROUND: Pinelliae Rhizoma (PR), a toxic medication, with long history, is commonly used for eliminating phlegm. Due to the shortage of wild resources and the relative lacking of cultivation technology, it is often confused with its counterfeit species in the market, such as Typhonii Rhizoma (TR), Arisaematis Rhizoma (AR) and tubers of Typhonium flagelliforme (TF) and Pinellia pedatisecta (PP). PURPOSE: It was aimed to screen signature enzymatic peptides from toxic proteins to identify PR and its four counterfeit species. STUDY DESIGN: A comparative proteogenomics strategy based on open-source transcriptome data was applied for screening signature peptides from toxic proteins, which were applied for species authentication of PR and its counterfeit species. METHODS: Firstly, the open-source transcriptome data was used for constructing the annotated protein database, which was used for peptides identification. Secondly, the toxicity of different fractions of PR were evaluated by the rat peritoneal inflammation model. Furthermore, sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) were used to profile the main proteins bands of five species, whose sequences were identified based on the in-gel digestion experiment by using ultra-high-performance liquid chromatography/quadrupole-Orbitrap mass spectrometry. Finally, the label-free proteomic analysis was performed to character the proteins and screen the signature peptides of five species, which were validated in commercially available products by dynamic multi reaction monitor (DMRM). RESULTS: The results in this study confirmed that protein was the main toxic components of PR. Both Pinellia ternata agglutinin (PTA) and trypsin inhibitor (TI) like proteins are the main proteins, which were characterized by proteomic analysis based on four annotated protein database. Meanwhile, seven signature peptides from toxic proteins were screened and validated with good repeatability and specificity in commercial products. CONCLUSION: Seven signature enzymatic peptides from toxic protein screened by the comparative proteogenomics strategy based on open-source transcriptome data achieved good identification ability of PR and its four counterfeit species.

Microbiota modulates gut immunity and promotes baculovirus infection in Helicoverpa armigera.

Baculoviruses are natural enemies of agricultural and forest insect pests and play an important role in biological pest control. Oral infection by baculovirus in the insect midgut is necessary for establishing systemic infection and eventually killing the insect. Since the insect midgut continuously encounters microbiota, the gut microbiota could affect baculovirus infection. Here, we demonstrated that gut microbiota modulates immune responses and promotes baculovirus infection in the cotton bollworm, Helicoverpa armigera. After oral infection, numerous host immunity-related genes including genes encoding Toll and immune deficiency (IMD) pathway components were upregulated in the midgut. Elimination of the gut microbiota significantly increased the resistance to viral infection in H. armigera. Quantitative real-time reverse transcription polymerase chain reaction and proteomic analysis showed that downregulation of the antiviral factor prophenoloxidase (PPO) could be mediated by microbiota during infection. It implied that midgut microbiota diminishes the expression of PPO to facilitate viral infection in H. armigera. Our findings revealed that the microbiota plays an important role in modulating the resistance of H. armigera to baculovirus infection, providing new insights in applying biopesticide.

Gene expression patterns of Cydia pomonella granulovirus in codling moth larvae revealed by RNAseq analysis.

The Cydia pomonella granulovirus (CpGV) has been used as a biological control agent of codling moth (Cydia pomonella), a severe global pest on pome fruit. Despite the economic importance, our knowledge of its molecular biology is still limited and a detailed picture of its gene expression is still missing. Here, we sequenced the transcriptome of codling moth larvae infected with the Mexican isolate CpGV-M and analyzed the expression of viral genes at 12, 48, and 96 h post infection (hpi). The results showed that two genes (p6.9 and pp31/39K) related to DNA binding of virus production, were highly expressed at 48 and 96 hpi. From 48 to 96 hpi, the expression of genes associated with virus replication and dissemination decreased, whereas the expression of genes related to infectious virion production and per os infectivity increased. This study provides a comprehensive view of CpGV gene expression patterns in host larvae.

Comparative Genomics Provide Insights Into Function and Evolution of Odorant Binding Proteins in Cydia pomonella.

Insect olfaction is vital for foraging, mating, host-seeking, and avoidance of predators/pathogens. In insects, odorant binding proteins (OBPs) are involved in transporting hydrophobic odor molecules from the external environment to receptor neurons. The codling moth, Cydia pomonella, one of the most destructive insect fruit pests, causes enormous economic losses. However, little is known about the number, variety, gains and losses, and evolution of OBP genes in C. pomonella. Here we report the identification of 40 OBPs in C. pomonella, most (75%) of which are classic OBPs, using genomic and transcriptomic analyses. Two OBP genes were lost in C. pomonella relative to possible distant ancestor in Lepidoptera lineage based on an analysis of gene gains and losses. The phylogenetic tree and chromosome location showed that the expansion of OBP genes mainly resulted from tandem duplications, as the CpomGOBP2 gene was duplicated twice along with loss of CpomPBPB. Two positive selection sites of the CpomGOBP1 gene were identified while other OBP genes evolved under purifying selection. Our results provide fundamental knowledge of OBP genes allowing further study of their function in C. pomonella.