Jekyll or Hyde? The genome (and more) of Nesidiocoris tenuis, a zoophytophagous predatory bug that is both a biological control agent and a pest

The First FISH-Confirmed Non-Canonical Telomeric Motif in Heteroptera: Cimex lectularius Linnaeus, 1758 and C. hemipterus (Fabricius, 1803) (Hemiptera, Cimicidae) Have a 10 bp Motif (TTAGGGATGG)n

Fluorescence in situ hybridization (FISH) with two different probes, the canonical insect telomeric sequence (TTAGG)n and the sequence (TTAGGGATGG)n, was performed on meiotic chromosomes of two members of the true bug family Cimicidae (Cimicomorpha), the common bed bug Cimex lectularius Linnaeus, 1758 and the tropical bed bug C. hemipterus (Fabricius, 1803), whose telomeric motifs were not known. In both species, there were no hybridization signals with the first probe, but strong signals at chromosomal ends were observed with the second probe, indicating the presence of a telomeric motif (TTAGGGATGG)n. This study represents the first FISH confirmation of the presence of a non-canonical telomeric motif not only for the infraorder Cimicomorpha but also for the suborder Heteroptera (Hemiptera) as a whole. The present finding is of key significance for unraveling the evolutionary shifts in the telomeric sequences in this suborder.

Phenazine biosynthesis protein MoPhzF regulates appressorium formation and host infection through canonical metabolic and noncanonical signaling function in Magnaporthe oryzae

Microbe-produced secondary metabolite phenazine-1-carboxylic acid (PCA) facilitates pathogen virulence and defense mechanisms against competitors. Magnaporthe oryzae, a causal agent of the devastating rice blast disease, needs to compete with other phyllosphere microbes and overcome host immunity for successful colonization and infection. However, whether M. oryzae produces PCA or it has any other functions remains unknown. Here, we found that the MoPHZF gene encodes the phenazine biosynthesis protein MoPhzF, synthesizes PCA in M. oryzae, and regulates appressorium formation and host virulence. MoPhzF is likely acquired through an ancient horizontal gene transfer event and has a canonical function in PCA synthesis. In addition, we found that PCA has a role in suppressing the accumulation of host-derived reactive oxygen species (ROS) during infection. Further examination indicated that MoPhzF recruits both the endoplasmic reticulum membrane protein MoEmc2 and the regulator of G-protein signaling MoRgs1 to the plasma membrane (PM) for MoRgs1 phosphorylation, which is a critical regulatory mechanism in appressorium formation and pathogenicity. Collectively, our studies unveiled a canonical function of MoPhzF in PCA synthesis and a noncanonical signaling function in promoting appressorium formation and host infection.

High-quality genome of the zoophytophagous stink bug, Nesidiocoris tenuis, informs their food habit adaptation

The zoophytophagous stink bug, Nesidiocoris tenuis, is a promising natural enemy of micro-pests such as whiteflies and thrips. This bug possesses both phytophagous and entomophagous food habits, enabling it to obtain nutrition from both plants and insects. This trait allows us to maintain its population density in agricultural fields by introducing insectary plants, even when the pest prey density is extremely low. However, if the bugs' population becomes too dense, they can sometimes damage crop plants. This dual character seems to arise from the food preferences and chemosensation of this predator. To understand the genomic landscape of N. tenuis, we examined the whole genome sequence of a commercially available Japanese strain. We used long-read sequencing and Hi-C analysis to assemble the genome at the chromosomal level. We then conducted a comparative analysis of the genome with previously reported genomes of phytophagous and hematophagous stink bugs to focus on the genetic factors contributing to this species' herbivorous and carnivorous tendencies. Our findings suggest that the gustatory gene set plays a pivotal role in adapting to food habits, making it a promising target for selective breeding. Furthermore, we identified the whole genomes of microorganisms symbiotic with this species through genomic analysis. We believe that our results shed light on the food habit adaptations of N. tenuis and will accelerate breeding efforts based on new breeding techniques for natural enemy insects, including genomics and genome editing.

Microbiome of Zoophytophagous Biological Control Agent Nesidiocoris tenuis

Many insects are associated with endosymbionts that influence the feeding, reproduction, and distribution of their hosts. Although the small green mirid, Nesidiocoris tenuis (Reuter) (Hemiptera: Miridae), a zoophytophagous predator that feeds on plants as well as arthropods, is a globally important biological control agent, its microbiome has not been sufficiently studied. In the present study, we assessed the microbiome variation in 96 N. tenuis individuals from 14 locations throughout Japan, based on amplicon sequencing of the 16S ribosomal RNA gene. Nine major bacteria associated with N. tenuis were identified: Rickettsia, two strains of Wolbachia, Spiroplasma, Providencia, Serratia, Pseudochrobactrum, Lactococcus, and Stenotrophomonas. Additionally, a diagnostic PCR analysis for three typical insect reproductive manipulators, Rickettsia, Wolbachia, and Spiroplasma, was performed on a larger sample size (n = 360) of N. tenuis individuals; the most prevalent symbiont was Rickettsia (69.7%), followed by Wolbachia (39.2%) and Spiroplasma (6.1%). Although some symbionts were co-infected, their prevalence did not exhibit any specific tendency, such as a high frequency in specific infection combinations. The infection frequency of Rickettsia was significantly correlated with latitude and temperature, while that of Wolbachia and Spiroplasma was significantly correlated with host plants. The predominance of these bacteria and the absence of obligate symbionts suggested that the N. tenuis microbiome is typical for predatory arthropods rather than sap-feeding insects. Rickettsia and Wolbachia were vertically transmitted rather than horizontally transmitted from the prey. The functional validation of each symbiont would be warranted to develop N. tenuis as a biological control agent.

High-quality chromosome-level scaffolds of the plant bug Pachypeltis micranthus provide insights into the availability of Mikania micrantha control

BACKGROUND

The plant bug, Pachypeltis micranthus Mu et Liu (Hemiptera: Miridae), is an effective potential biological control agent for Mikania micrantha H.B.K. (Asteraceae; one of the most notorious invasive weeds worldwide). However, limited knowledge about this species hindered its practical application and research. Accordingly, sequencing the genome of this mirid bug holds great significance in controlling M. micrantha.

RESULTS

Here, 712.72 Mb high-quality chromosome-level scaffolds of P. micranthus were generated, of which 707.51 Mb (99.27%) of assembled sequences were anchored onto 15 chromosome-level scaffolds with contig N50 of 16.84 Mb. The P. micranthus genome had the highest GC content (42.43%) and the second highest proportion of repetitive sequences (375.82 Mb, 52.73%) than the three other mirid bugs (i.e., Apolygus lucorum, Cyrtorhinus lividipennis, and Nesidiocoris tenuis). Phylogenetic analysis showed that P. micranthus clustered with other mirid bugs and diverged from the common ancestor approximately 200 million years ago. Gene family expansion and/or contraction were analyzed, and significantly expanded gene families associated with P. micranthus feeding and adaptation to M. micrantha were manually identified. Compared with the whole body, transcriptome analysis of the salivary gland revealed that most of the upregulated genes were significantly associated with metabolism pathways and peptidase activity, particularly among cysteine peptidase, serine peptidase, and polygalacturonase; this could be one of the reasons for precisely and highly efficient feeding by the oligophagous bug P. micranthus on M. micrantha.

CONCLUSION

Collectively, this work provides a crucial chromosome-level scaffolds resource to study the evolutionary adaptation between mirid bug and their host. It is also helpful in searching for novel environment-friendly biological strategies to control M. micrantha.

Genome sequence and silkomics of the spindle ermine moth, Yponomeuta cagnagella, representing the early diverging lineage of the ditrysian Lepidoptera

Many lepidopteran species produce silk, cocoons, feeding tubes, or nests for protection from predators and parasites for caterpillars and pupae. Yet, the number of lepidopteran species whose silk composition has been studied in detail is very small, because the genes encoding the major structural silk proteins tend to be large and repetitive, making their assembly and sequence analysis difficult. Here we have analyzed the silk of Yponomeuta cagnagella, which represents one of the early diverging lineages of the ditrysian Lepidoptera thus improving the coverage of the order. To obtain a comprehensive list of the Y. cagnagella silk genes, we sequenced and assembled a draft genome using Oxford Nanopore and Illumina technologies. We used a silk-gland transcriptome and a silk proteome to identify major silk components and verified the tissue specificity of expression of individual genes. A detailed annotation of the major genes and their putative products, including their complete sequences and exon-intron structures is provided. The morphology of silk glands and fibers are also shown. This study fills an important gap in our growing understanding of the structure, evolution, and function of silk genes and provides genomic resources for future studies of the chemical ecology of Yponomeuta species.

The Adipokinetic Peptides of Hemiptera: Structure, Function, and Evolutionary Trends

The Hemiptera comprise the most species-rich order of the hemimetabolous insects. Members of a number of superfamilies, most notably especially the more basal ones such as white flies, psyllids and aphids, belong to the most destructive agricultural insects known worldwide. At the other end of the phylogenetic tree are hemipterans that are notorious medical pests (e.g. kissing bugs). Most of the hemipteran species are good flyers, and lipid oxidation plays a pivotal role to power the contraction of flight muscles and, in aquatic water bugs, also deliver the ATP for the extensive swimming action of the leg muscles. Mobilization of stored lipids (mostly triacylglycerols in the fat body) to circulating diacylglycerols in the hemolymph is regulated by a set of small neuropeptides, the adipokinetic hormones (AKHs). We searched the literature and publicly available databases of transcriptomes and genomes to present here AKH sequences from 191 hemipteran species. Only few of these peptides were sequenced via Edman degradation or mass spectrometry, and even fewer were characterized with molecular biology methods; thus, the majority of the AKHs we have identified by bioinformatics are merely predicted sequences at this stage. Nonetheless, a total of 42 AKH primary sequences are assigned to Hemiptera. About 50% of these structures occur also in other insect orders, while the remaining 50% are currently unique for Hemiptera. We find 9 novel AKHs not shown to be synthesized before in any insect. Most of the hemipteran AKHs are octapeptides (28) but there is an impressive number of decapeptides (12) compared to other speciose orders such as Diptera and Lepidoptera. We attempt to construct a hypothetical molecular peptide evolution of hemipteran AKHs and find quite a bit of overlapping with current phylogenetic ideas of the Hemiptera. Lastly, we discuss the possibility to use the sequence of the aphid AKH as lead peptide for the research into a peptide mimetic fulfilling criteria of a green insecticide.

Chromosome-level genome assembly of the mirid predator Cyrtorhinus lividipennis Reuter (Hemiptera: Miridae), an important natural enemy in the rice ecosystem

Though the genomes of many rice herbivorous pests have recently been well characterized, little is known about the genome of their natural enemies. Here, by using the Illumina and PacBio platforms, we sequenced and assembled the whole genome of the mirid species Cyrtorhinus lividipennis Reuter (Hemiptera: Miridae), which is an economically and ecologically important natural enemy in the rice ecosystem acting as a dominant predator for planthoppers and leafhoppers in the field. Through Hi-C scaffolding, 1615 scaffolds with a total size of 338.08 Mb were successfully anchored onto 13 chromosomes. The assembled genome size was 345.75 Mb with a final scaffold N50 of 27.58 Mb. Approximately 107.51 Mb of sequences accounting for 31.10% of the genome were identified as repeat elements, and 14,644 protein-coding genes were annotated. Phylogenetic analysis showed that C. lividipennis clustered with other Hemipteran species and diverged from Apolygus lucorum about 66.7 million years ago. Gene families related to detoxification, environmental adaptation and digestion were analysed comparatively with other Hemipteran species, but no significant expansion or contraction was found in C. lividipennis. We also observed male meiosis in C. lividipennis, which showed a typical post-reduction of sex chromosomes and a karyotype of 2n = 22 + XY. As the first natural-enemy genome in the rice ecosystem, the genomic resource of C. lividipennis not only expands our understanding of the multitrophic interactions (host plant-prey-predator), but also provides a genomic basis for better understanding this dominant predator and therefore promotes sustainable rice pest management and food grain production.

A conserved role of the duplicated Masculinizer gene in sex determination of the Mediterranean flour moth, Ephestia kuehniella

Sex determination in the silkworm, Bombyx mori, is based on Feminizer (Fem), a W-linked Fem piRNA that triggers female development in WZ individuals, and the Z-linked Masculinizer (Masc), which initiates male development and dosage compensation in ZZ individuals. While Fem piRNA is missing in a close relative of B. mori, Masc determines sex in several representatives of distant lepidopteran lineages. We studied the molecular mechanisms of sex determination in the Mediterranean flour moth, Ephestia kuehniella (Pyralidae). We identified an E. kuehniella Masc ortholog, EkMasc, and its paralog resulting from a recent duplication, EkMascB. Both genes are located on the Z chromosome and encode a similar Masc protein that contains two conserved domains but has lost the conserved double zinc finger domain. We developed PCR-based genetic sexing and demonstrated a peak in the expression of EkMasc and EkMascB genes only in early male embryos. Simultaneous knock-down experiments of both EkMasc and EkMascB using RNAi during early embryogenesis led to a shift from male- to female-specific splicing of the E. kuehniella doublesex gene (Ekdsx), their downstream effector, in ZZ embryos and resulted in a strong female-biased sex-ratio. Our results thus confirmed the conserved role of EkMasc and/or EkMascB in masculinization. We suggest that the C-terminal proline-rich domain, we have identified in all functionally confirmed Masc proteins, in conjunction with the masculinizing domain, is important for transcriptional regulation of sex determination in Lepidoptera. The function of the Masc double zinc finger domain is still unknown, but appears to have been lost in E. kuehniella.

The sex-determining cascade in the silkworm, Bombyx mori, differs greatly from those of other insects. In B. mori, female development is initiated by Fem piRNA expressed from the W chromosome during early embryogenesis. Fem piRNA silences Masculinizer (Masc) thereby blocking the male pathway resulting in female development. It is currently unknown whether this cascade is conserved across Lepidoptera. In the Mediterranean flour moth, Ephestia kuehniella, we identified an ortholog of Masc and discovered its functional duplication on the Z chromosome, which has not yet been found in any other lepidopteran species. We provide two lines of evidence that the EkMasc and/or EkMascB genes play an essential role in masculinization: (i) they show a peak of expression during early embryogenesis in ZZ but not in WZ embryos and (ii) their simultaneous silencing by RNAi results in female-specific splicing of the E. kuehniella doublesex gene (Ekdsx) in ZZ embryos and in a female-biased sex ratio. Our results suggest a conserved role of the duplicated Masc gene in sex determination of E. kuehniella.

Large-scale comparative analysis of cytogenetic markers across Lepidoptera

Fluorescence in situ hybridization (FISH) allows identification of particular chromosomes and their rearrangements. Using FISH with signal enhancement via antibody amplification and enzymatically catalysed reporter deposition, we evaluated applicability of universal cytogenetic markers, namely 18S and 5S rDNA genes, U1 and U2 snRNA genes, and histone H3 genes, in the study of the karyotype evolution in moths and butterflies. Major rDNA underwent rather erratic evolution, which does not always reflect chromosomal changes. In contrast, the hybridization pattern of histone H3 genes was well conserved, reflecting the stable organisation of lepidopteran genomes. Unlike 5S rDNA and U1 and U2 snRNA genes which we failed to detect, except for 5S rDNA in a few representatives of early diverging lepidopteran lineages. To explain the negative FISH results, we used quantitative PCR and Southern hybridization to estimate the copy number and organization of the studied genes in selected species. The results suggested that their detection was hampered by long spacers between the genes and/or their scattered distribution. Our results question homology of 5S rDNA and U1 and U2 snRNA loci in comparative studies. We recommend the use of histone H3 in studies of karyotype evolution.