A chromosome-level genome assembly provides new insights into paternal genome elimination in the cotton mealybug Phenacoccus solenopsis

Silencing PsASH2 affects embryo development in the cotton mealybug

Mealybugs are highly aggressive pests that infest various plants and cause substantial economic losses. Histone lysine methyltransferases (KMT) are evolutionarily conserved and proposed to be essential in early embryo development in animals. However, few KMTs have been reported in mealybugs. Here, we identified a novel KMT gene, PsASH2, in the cotton mealybug, Phenacoccus solenopsis Tinsley. This gene was highly expressed in the ovary of female adults. Through RNA interference (RNAi) of PsASH2 by dsRNA microinjection, we found a reduction in the number of male embryos and total embryos in the ovaries of pregnant females. Continuous downregulation of PsASH2 in mated females until their death resulted in few changes in sex ratio but significant decreases in the number of both male and female offspring. Therefore, we believe that PsASH2 plays essential roles in embryo survival for both sexes of the cotton mealybug which may provide a potential target gene for the management of cotton mealybug by disrupting embryo development.

The B Chromosome of Pseudococcus viburni: A Selfish Chromosome that Exploits Whole-Genome Meiotic Drive

Meiosis is generally a fair process: each chromosome has a 50% chance of being included into each gamete. However, meiosis can become aberrant with some chromosomes having a higher chance of making it into gametes than others. Yet, why and how such systems evolve remains unclear. Here, we study the unusual reproductive genetics of mealybugs, where only maternal-origin chromosomes are included in gametes during male meiosis, while paternal chromosomes are eliminated. One species-Pseudococcus viburni-has a segregating B chromosome that drives by escaping paternal genome elimination. We present whole genome and gene expression data from lines with and without B chromosomes. We identify B-linked sequences including 204 protein-coding genes and a satellite repeat that makes up a significant proportion of the chromosome. The few paralogs between the B and the core genome are distributed throughout the genome, arguing against a simple, or at least recent, chromosomal duplication of one of the autosomes to create the B. We do, however, find one 373 kb region containing 146 genes that appears to be a recent translocation. Finally, we show that while many B-linked genes are expressed during meiosis, most of these are encoded on the recently translocated region. Only a small number of B-exclusive genes are expressed during meiosis. Of these, only one was overexpressed during male meiosis, which is when the drive occurs: an acetyltransferase involved in H3K56Ac, which has a putative role in meiosis and is, therefore, a promising candidate for further studies.

Genome report: Genome sequence of tuliptree scale, Toumeyella liriodendri (Gmelin), an ornamental pest insect

Scale insects are of interest both to basic researchers for their unique reproductive biology and to applied researchers for their pest status. In spite of this interest, there remain few genomic resources for this group of insects. To begin addressing this lack of data, we present the genome sequence of tuliptree scale, Toumeyella liriodendri (Gmelin) (Hemiptera: Coccomorpha: Coccidae). The genome assembly spans 536Mb, with over 96% of sequence assembled into one of 17 chromosomal scaffolds. We characterize roughly 66% of this sequence as repetitive and annotate 16,508 protein coding genes. Then we use the reference genome to explore the phylogeny of soft scales (Coccidae) and evolution of karyotype within the family. We find that T. liriodendri is an early-diverging soft scale, less closely related to most sequenced soft scales than a species of the family Aclerdidae is. This molecular result corroborates a previous, morphology-based phylogenetic placement of Aclerdidae within Coccidae. In terms of genome structure, T. liriodendri has nearly twice as many chromosomes as the only other soft scale assembled to the chromosome level, Ericerus pela (Chavannes). In comparing the two, we find that chromosome number evolution can largely be explained by simple fissions rather than more complex rearrangements. These genomic natural history observations lay a foundation for further exploration of this unique group of insects.

Essential roles of histone lysine methyltransferases EZH2 and EHMT1 in male embryo development of Phenacoccus solenopsis

Paternal genome elimination (PGE) is an intriguing but poorly understood reproductive strategy in which females are typically diploid, but males lose paternal genomes. Paternal genome heterochromatin (PGH) occurs in arthropods with germline PGE, such as the mealybug, coffee borer beetles, and booklice. Here, we present evidence that PGH initially occurs during early embryo development at around 15 h post-mating (hpm) in the cotton mealybug, Phenacoccus solenopsis Tinsley. Transcriptome analysis followed by qPCR validation indicated that six histone lysine methyltransferase (KMT) genes are predominantly expressed in adult females. We knocked down these five genes through dsRNA microinjection. We found that downregulation of two KMT genes, PsEZH2-X1 and PsEHMT1, resulted in a decrease of heterochromatin-related methylations, including H3K27me1, H3K27me3, and H3K9me3 in the ovaries, fewer PGH male embryos, and reduced male offspring. For further confirmation, we obtained two strains of transgenic tobacco highly expressing dsRNA targeting PsEZH2-X1 and PsEHMT1, respectively. Similarly, fewer PGH embryos and fewer male offspring were observed when feeding on these transgenic tobacco plants. Overall, we present evidence that PsEZH2-X1 and PsEHMT1 have essential roles in male embryo survival by regulating PGH formation in cotton mealybugs.

Chromosome-level genome assembly of the cottony cushion scale Icerya purchasi

The cottony cushion scale, Icerya purchasi, a polyphagous pest, poses a significant threat to the global citrus industry. The hermaphroditic self-fertilization observed in I. purchasi is an exceptionally rare reproductive mode among insects. In this study, we successfully assembled a chromosome-level genome sequence for I. purchasi using PacBio long-reads and the Hi-C technique, resulting in a total size of 1,103.38 Mb and a contig N50 of 12.81 Mb. The genome comprises 14,046 predicted protein-coding genes, with 462,722,633 bp occurrence of repetitive sequences. BUSCO analysis revealed a completeness score of 93.20%. The genome sequence of I. purchasi serves as a crucial resource for comprehending the reproductive modes in insects, with particular emphasis on hermaphroditic self-fertilization.

Chromosome-level genome assembly of the Asian citrus psyllid, Diaphorina citri

Diaphorina citri is a global citrus pest. As a vector insect, it can transmit the causative agents of citrus huanglongbing, causing irreversible losses to the citrus industry. The acquisition of genomic information can provide a molecular genetic basis for effective control of D. citri. Here, the DNBSEQ™ , Oxford Nanopore Technologies, and Hi-C technologies are applied to generate a high-quality chromosome-level genome of D. citri. The genome size of D. citri was 523.78 Mb with a scaffold N50 of 47.05 Mb distributed on 13 chromosomes. A total of 250.64 Mb (47.85%) repeat sequences and 24 048 protein-coding genes were predicted. Genome resequencing of female and male individuals indicated that the sex chromosome system of D. citri is XO. Phylogenetic analysis demonstrated that D. citri and Pachypsylla venusta, which separated from their most recent common ancestor about 336.62 million years ago, were the most closely related. Additionally, we identified genes potentially involved in detoxification metabolism, pathogen transmission, and honeydew secretion for further investigation. The high-quality genome provides an important reference for developing effective management strategies of D. citri.

Endosymbiont Tremblaya phenacola influences the reproduction of cotton mealybugs by regulating the mechanistic target of rapamycin pathway

The intricate evolutionary dynamics of endosymbiotic relationships result in unique characteristics among the genomes of symbionts, which profoundly influence host insect phenotypes. Here, we investigated an endosymbiotic system in Phenacoccus solenopsis, a notorious pest of the subfamily Phenacoccinae. The endosymbiont, "Candidatus Tremblaya phenacola" (T. phenacola PSOL), persisted throughout the complete life cycle of female hosts and was more active during oviposition, whereas there was a significant decline in abundance after pupation in males. Genome sequencing yielded an endosymbiont genome of 221.1 kb in size, comprising seven contigs and originating from a chimeric arrangement between betaproteobacteria and gammaproteobacteria. A comprehensive analysis of amino acid metabolic pathways demonstrated complementarity between the host and endosymbiont metabolism. Elimination of T. phenacola PSOL through antibiotic treatment significantly decreased P. solenopsis fecundity. Weighted gene coexpression network analysis demonstrated a correlation between genes associated with essential amino acid synthesis and those associated with host meiosis and oocyte maturation. Moreover, altering endosymbiont abundance activated the host mechanistic target of rapamycin pathway, suggesting that changes in the amino acid abundance affected the host reproductive capabilities via this signal pathway. Taken together, these findings demonstrate a mechanism by which the endosymbiont T. phenacola PSOL contributed to high fecundity in P. solenopsis and provide new insights into nutritional compensation and coevolution of the endosymbiotic system.

Whole-genome sequence of the Cooley spruce gall adelgid, Adelges cooleyi (Hemiptera: Sternorrhyncha: Adelgidae)

The adelgids (Adelgidae) are a small family of sap-feeding insects, which, together with true aphids (Aphididae) and phylloxerans (Phylloxeridae), make up the infraorder Aphidomorpha. Some adelgid species are highly destructive to forest ecosystems such as Adelges tsugae, Adelges piceae, Adelges laricis, Pineus pini, and Pineus boerneri. Despite this, there are no high-quality genomic resources for adelgids, hindering advanced genomic analyses within Adelgidae and among Aphidomorpha. Here, we used PacBio continuous long-read and Illumina RNA-sequencing to construct a high-quality draft genome assembly for the Cooley spruce gall adelgid, Adelges cooleyi (Gillette), a gall-forming species endemic to North America. The assembled genome is 270.2 Mb in total size and has scaffold and contig N50 statistics of 14.87 and 7.18 Mb, respectively. There are 24,967 predicted coding sequences, and the assembly completeness is estimated at 98.1 and 99.6% with core BUSCO gene sets of Arthropoda and Hemiptera, respectively. Phylogenomic analysis using the A. cooleyi genome, 3 publicly available adelgid transcriptomes, 4 phylloxera transcriptomes, the Daktulosphaira vitifoliae (grape phylloxera) genome, 4 aphid genomes, and 2 outgroup coccoid genomes fully resolves adelgids and phylloxerans as sister taxa. The mitochondrial genome is 24 kb, among the largest in insects sampled to date, with 39.4% composed of noncoding regions. This genome assembly is currently the only genome-scale, annotated assembly for adelgids and will be a valuable resource for understanding the ecology and evolution of Aphidomorpha.

Unraveling the complex evolutionary history of lepidopteran chromosomes through ancestral chromosome reconstruction and novel chromosome nomenclature

BACKGROUND

Lepidoptera is one of the most species-rich animal groups, with substantial karyotype variations among species due to chromosomal rearrangements. Knowledge of the evolutionary patterns of lepidopteran chromosomes still needs to be improved.

RESULTS

Here, we used chromosome-level genome assemblies of 185 lepidopteran insects to reconstruct an ancestral reference genome and proposed a new chromosome nomenclature. Thus, we renamed over 5000 extant chromosomes with this system, revealing the historical events of chromosomal rearrangements and their features. Additionally, our findings indicate that, compared with autosomes, the Z chromosome in Lepidoptera underwent a fast loss of conserved genes, rapid acquisition of lineage-specific genes, and a low rate of gene duplication. Moreover, we presented evidence that all available 67 W chromosomes originated from a common ancestor chromosome, with four neo-W chromosomes identified, including one generated by fusion with an autosome and three derived through horizontal gene transfer. We also detected nearly 4000 inter-chromosomal gene movement events. Notably, Geminin is transferred from the autosome to the Z chromosome. When located on the autosome, Geminin shows female-biased expression, but on the Z chromosome, it exhibits male-biased expression. This contributes to the sexual dimorphism of body size in silkworms.

CONCLUSIONS

Our study sheds light on the complex evolutionary history of lepidopteran chromosomes based on ancestral chromosome reconstruction and novel chromosome nomenclature.

Paternal genome elimination: patterns and mechanisms of drive and silencing

In thousands of arthropod species, males inherit, but subsequently eliminate the entire haploid genome of their father. However, why this peculiar reproductive strategy evolved repeatedly across diverse species and what mechanisms are involved in paternal genome elimination (PGE) remains largely unknown. In this review, we summarize what we know about the patterns of paternal chromosome elimination during various stages of development in the diverse taxa that have been studied. We also discuss some other unusual features often associated with PGE, such as the transcriptional silencing of paternally derived chromosomes in males and sex determination through the early embryonic elimination of X chromosomes. Little is known about the molecular mechanisms underlying the parent-of-origin-dependent chromosome elimination and silencing under PGE, but we discuss the insight of several studies that are pioneering this work and highlight directions for future research.

The transcriptome of Icerya aegyptiaca (Hemiptera: Monophlebidae) and comparison with neococcoids reveal genetic clues of evolution in the scale insects

BACKGROUND

Scale insects are worldwide sap-sucking parasites, which can be distinguished into neococcoids and non-neococcoids. Neococcoids are monophyletic with a peculiar reproductive system, paternal genome elimination (PGE). Different with neococcoids, Iceryini, a tribe in non-neococcoids including several damaging pests, has abdominal spiracles, compound eyes in males, relatively abundant wax, unique hermaphrodite system, and specific symbionts. However, the current studies on the gene resources and genomic mechanism of scale insects are mainly limited in the neococcoids, and lacked of comparison in an evolution frame.

RESULT

We sequenced and de novo assembled a transcriptome of Icerya aegyptiaca (Douglas), a worldwide pest of Iceryini, and used it as representative of non-neococcoids to compare with the genomes or transcriptomes of other six species from different families of neococcoids. We found that the genes under positive selection or negative selection intensification (simplified as "selected genes" below) in I. aegyptiaca included those related to neurogenesis and development, especially eye development. Some genes related to fatty acid biosynthesis were unique in its transcriptome with relatively high expression and not detected in neococcoids. These results may indicate a potential link to the unique structures and abundant wax of I. aegyptiaca compared with neococcoids. Meanwhile, genes related to DNA repair, mitosis, spindle, cytokinesis and oogenesis, were included in the selected genes in I. aegyptiaca, which is possibly associated with cell division and germ cell formation of the hermaphrodite system. Chromatin-related process were enriched from selected genes in neococcoids, along with some mitosis-related genes also detected, which may be related to their unique PGE system. Moreover, in neococcoid species, male-biased genes tend to undergo negative selection relaxation under the PGE system. We also found that the candidate horizontally transferred genes (HTGs) in the scale insects mainly derived from bacteria and fungi. bioD and bioB, the two biotin-synthesizing HTGs were exclusively found in the scale insects and neococcoids, respectively, which possibly show potential demand changes in the symbiotic relationships.

CONCLUSION

Our study reports the first I. aegyptiaca transcriptome and provides preliminary insights for the genetic change of structures, reproductive systems and symbiont relationships at an evolutionary aspect. This will provide a basis for further research and control of scale insects.

Fatty acyl-CoA reductase influences wax biosynthesis in the cotton mealybug, Phenacoccus solenopsis Tinsley

Mealybugs are highly aggressive to a diversity of plants. The waxy layer covering the outermost part of the integument is an important protective defense of these pests. However, the molecular mechanisms underlying wax biosynthesis in mealybugs remain largely unknown. Here, we analyzed multi-omics data on wax biosynthesis by the cotton mealybug, Phenacoccus solenopsis Tinsley, and found that a fatty acyl-CoA reductase (PsFAR) gene, which was highly expressed in the fat bodies of female mealybugs, contributed to wax biosynthesis by regulating the production of the dominant chemical components of wax, cuticular hydrocarbons (CHCs). RNA interference (RNAi) against PsFAR by dsRNA microinjection and allowing mealybugs to feed on transgenic tobacco expressing target dsRNA resulted in a reduction of CHC contents in the waxy layer, and an increase in mealybug mortality under desiccation and deltamethrin treatments. In conclusion, PsFAR plays crucial roles in the wax biosynthesis of mealybugs, thereby contributing to their adaptation to water loss and insecticide stress.

The role of a fatty acyl-CoA reductase (PsFAR) in wax biosynthesis of cotton mealybug is investigated, RNAi against PsFAR resulted in insects with lower generation of waxy filaments and higher mortality under desiccation and deltamethrin treatments.

Gene Editing and Genetic Control of Hemipteran Pests: Progress, Challenges and Perspectives

The origin of the order Hemiptera can be traced to the late Permian Period more than 230 MYA, well before the origin of flowering plants 100 MY later in during the Cretaceous period. Hemipteran species consume their liquid diets using a sucking proboscis; for phytophagous hemipterans their mouthparts (stylets) are elegant structures that enable voracious feeding from plant xylem or phloem. This adaptation has resulted in some hemipteran species becoming globally significant pests of agriculture resulting in significant annual crop losses. Due to the reliance on chemical insecticides for the control of insect pests in agricultural settings, many hemipteran pests have evolved resistance to insecticides resulting in an urgent need to develop new, species-specific and environmentally friendly methods of pest control. The rapid advances in CRISPR/Cas9 technologies in model insects such as Drosophila melanogaster, Tribolium castaneum, Bombyx mori, and Aedes aegypti has spurred a new round of innovative genetic control strategies in the Diptera and Lepidoptera and an increased interest in assessing genetic control technologies for the Hemiptera. Genetic control approaches in the Hemiptera have, to date, been largely overlooked due to the problems of introducing genetic material into the germline of these insects. The high frequency of CRISPR-mediated mutagenesis in model insect species suggest that, if the delivery problem for Hemiptera could be solved, then gene editing in the Hemiptera might be quickly achieved. Significant advances in CRISPR/Cas9 editing have been realized in nine species of Hemiptera over the past 4 years. Here we review progress in the Hemiptera and discuss the challenges and opportunities for extending contemporary genetic control strategies into species in this agriculturally important insect orderr.

Genome and transcriptome analysis of the mealybug Maconellicoccus hirsutus: Correlation with its unique phenotypes

Mealybugs are aggressive pests with world-wide distribution and are suitable for the study of different phenomena like genomic imprinting and epigenetics. Genomic approaches facilitate these studies in absence of robust genetics in this system. We sequenced, de novo assembled, annotated Maconellicoccus hirsutus genome. We carried out comparative genomics it with four mealybug and eight other insect species, to identify expanded, specific and contracted gene classes that relate to pesticide and desiccation resistance. We identified horizontally transferred genes adding to the mutualism between the mealybug and its endosymbionts. Male and female transcriptome analysis indicates differential expression of metabolic pathway genes correlating with their physiology and the genes for sexual dimorphism. The significantly lower expression of endosymbiont genes in males relates to the depletion of endosymbionts in males during development.

Males That Silence Their Father's Genes: Genomic Imprinting of a Complete Haploid Genome

Genetic conflict is considered a key driver in the evolution of reproductive systems with non-Mendelian inheritance, where parents do not contribute equally to the genetic makeup of their offspring. One of the most extraordinary examples of non-Mendelian inheritance is paternal genome elimination (PGE), a form of haplodiploidy which has evolved repeatedly across arthropods. Under PGE, males are diploid but only transmit maternally inherited chromosomes, while the paternally inherited homologues are excluded from sperm. This asymmetric inheritance is thought to have evolved through an evolutionary arms race between the paternal and maternal genomes over transmission to future generations. In several PGE clades, such as the mealybugs (Hemiptera: Pseudococcidae), paternal chromosomes are not only eliminated from sperm, but also heterochromatinized early in development and thought to remain inactive, which could result from genetic conflict between parental genomes. Here, we present a parent-of-origin allele-specific transcriptome analysis in male mealybugs showing that expression is globally biased toward the maternal genome. However, up to 70% of somatically expressed genes are to some degree paternally expressed, while paternal genome expression is much more restricted in the male reproductive tract, with only 20% of genes showing paternal contribution. We also show that parent-of-origin-specific gene expression patterns are remarkably similar across genotypes, and that genes with completely biparental expression show elevated rates of molecular evolution. Our results provide the clearest example yet of genome-wide genomic imprinting in insects and enhance our understanding of PGE, which will aid future empirical tests of evolutionary theory regarding the origin of this unusual reproductive strategy.