Parental experiences orchestrate locust egg hatching synchrony by regulating nuclear export of precursor miRNA.

Parental experiences can affect the phenotypic plasticity of offspring. In locusts, the population density that adults experience regulates the number and hatching synchrony of their eggs, contributing to locust outbreaks. However, the pathway of signal transmission from parents to offspring remains unclear. Here, we find that transcription factor Forkhead box protein N1 (FOXN1) responds to high population density and activates the polypyrimidine tract-binding protein 1 (Ptbp1) in locusts. FOXN1-PTBP1 serves as an upstream regulator of miR-276, a miRNA to control egg-hatching synchrony. PTBP1 boosts the nucleo-cytoplasmic transport of pre-miR-276 in a "CU motif"-dependent manner, by collaborating with the primary exportin protein exportin 5 (XPO5). Enhanced nuclear export of pre-miR-276 elevates miR-276 expression in terminal oocytes, where FOXN1 activates Ptbp1 and leads to egg-hatching synchrony in response to high population density. Additionally, PTBP1-prompted nuclear export of pre-miR-276 is conserved in insects, implying a ubiquitous mechanism to mediate transgenerational effects.

Chromosome-level genome assembly of the pygmy grasshopper Eucriotettix oculatus (Orthoptera: Tetrigoidea).

The pygmy grasshoppers, which belong to the superfamily Tetrigoidea, exhibit remarkable environmental adaptability. However, no study has yet reported a reference genome for this group. In this study, we assembled a high-quality chromosome-scale genome of Eucriotettix oculatus, which survive in the environment heavily polluted by heavy metals, achieved through Illumina and PacBio sequencing, alongside chromosome conformation capture techniques. The resulting genome spans 985.45 Mb across seven chromosomes (range: 71.55 to 266.65 Mb) and features an N50 length of 123.82 Mb. Chr5 is considered to be the single sex chromosome (X). This genome is composed of 46.42% repetitive elements and contains 14,906 predicted protein-coding genes, 91.63% of which are functionally annotated. Decoding the E. oculatus genome not only promotes future studies on environmental adaptation for the pygmy grasshopper, but also provides valuable resources for in-depth investigation on phylogeny, evolution, and behavior of Orthoptera.

A neural m6A pathway regulates behavioral aggregation in migratory locusts.

RNA N6-methyladenosine (m6A), as the most abundant modification of messenger RNA, can modulate insect behaviors, but its specific roles in aggregation behaviors remain unexplored. Here, we conducted a comprehensive molecular and physiological characterization of the individual components of the methyltransferase and demethylase in the migratory locust Locusta migratoria. Our results demonstrated that METTL3, METTL14 and ALKBH5 were dominantly expressed in the brain and exhibited remarkable responses to crowding or isolation. The individual knockdown of methyltransferases (i.e., METTL3 and METTL14) promoted locust movement and conspecific attraction, whereas ALKBH5 knockdown induced a behavioral shift toward the solitary phase. Furthermore, global transcriptome profiles revealed that m6A modification could regulate the orchestration of gene expression to fine tune the behavioral aggregation of locusts. In summary, our in vivo characterization of the m6A functions in migratory locusts clearly demonstrated the crucial roles of the m6A pathway in effectively modulating aggregation behaviors.

Limited genetic variability and spatial population structure in grasshoppers between natural and metal-contaminated areas in Egypt.

Pollutants in an environment can have long-term implications for the species living there, resulting in local adaptations with implications for their genetic structure. Heavy metal pollutants infiltrate soils and groundwater, bioaccumulate in food webs, and negatively impact biota. In this study, we investigated the degree to which the genetic structure and variability of the slender green-winged grasshopper (Aiolopus thalassinus (Fabricius) (Orthoptera: Acrididae)) were impacted by heavy metal pollution and distance. We used the random amplified polymorphic DNA-polymerase chain reaction (RAPD-PCR) method to examine the genetic variability of populations in 3 heavy metal-polluted and 3 unpolluted locations across varying geographical distances in Egypt. The heavy metal concentrations of cadmium, copper, lead, and zinc were measured from the grasshopper tissue and soils. Sixty-nine unique and polymorphic bands were produced by 4 primers. Cluster and principal component analyses separated the populations inside and outside Cairo into 2 main branches, which were further divided into smaller branches corresponding to their geographical regions. We found no differences in the Shannon genetic diversity index between populations or with increasing heavy metal concentrations in either the soil or the grasshopper tissue. Our results showed a greater genetic variation among populations than between populations within the same location, indicating populations within locations were less differentiated than those between locations. The moderate correlation between genetic similarity and spatial distance suggests geographical isolation influenced grasshopper population differentiation. Based on the RAPD analysis, environmental pollutants and geographical distances impact the A. thalassinus population structure, potentially restricting gene flow between sites even at small spatial scales.

Phylogenetic insights into Central European Chorthippus and Pseudochorthippus (Orthoptera: Acrididae) species using ddRADseq data.

The evolution of several orthopteran groups, especially within the grasshopper family Acrididae, remains poorly understood. This is particularly true for the subfamily Gomphocerinae, which comprises cryptic sympatric and syntopic species. Previous mitochondrial studies have highlighted major discrepancies between taxonomic and phylogenetic hypotheses, thereby emphasizing the necessity of genome-wide approaches. In this study, we employ double-digest restriction site-associated DNA sequencing (ddRADseq) to reconstruct the evolution of Central European Chorthippus and Pseudochorthippus species, especially C.smardai, P.tatrae and the C.biguttulus group. Our phylogenomic analyses recovered deep discordance with mitochondrial DNA barcoding, emphasizing its unreliability in Gomphocerinae grasshoppers. Specifically, our data robustly distinguished the C.biguttulus group and confirmed the distinctiveness of C.eisentrauti, also shedding light on its presence in the Berchtesgaden Alps. Moreover, our results support the reclassification of C.smardai to the genus Pseudochorthippus and of P.tatrae to the genus Chorthippus. Our study demonstrates the efficiency of high-throughput genomic methods such as RADseq without prior optimization to elucidate the complex evolution of grasshopper radiations with direct taxonomic implications. While RADseq has predominantly been utilized for population genomics and within-genus phylogenomics, its application extends to resolve relationships between deeply-diverged clades representative of distinct genera.

Auditory robustness and resilience in the aging auditory system of the desert locust.

After overexposure to loud music, we experience a decrease in our ability to hear (robustness), which usually recovers (resilience). Here, we exploited the amenable auditory system of the desert locust, Schistocerca gregaria, to measure how robustness and resilience depend on age. We found that gene expression changes are dominated by age as opposed to noise exposure. We measured sound-evoked nerve activity for young and aged locusts directly, after 24 hours and 48 hours after noise exposure. We found that both young and aged locusts recovered their auditory nerve function over 48 hours. We also measured the sound-evoked transduction current in individual auditory neurons, and although the transduction current magnitude recovered in the young locusts after noise exposure, it failed to recover in the aged locusts. A plastic mechanism compensates for the decreased transduction current in aged locusts. We suggest key genes upregulated in young noise-exposed locusts that mediate robustness to noise exposure and find potential candidates responsible for compensatory mechanisms in the auditory neurons of aged noise-exposed locusts.

Comparative analysis of mitogenomes among three species of grasshoppers (Orthoptera: Acridoidea: Gomphocerinae) and their phylogenetic implications.

Whole mitochondrial genomes have been widely used in phylogenetic analysis, population genetics and biogeography studies. This study sequenced and characterized three complete mitochondrial genomes (Dasyhippus peipingensis, Myrmeleotettix palpalis, Aeropedellus prominemarginis) and determined their phylogenetic position in Acrididae. The length of the mitochondrial genomes ranged from 15,621-15,629 bp and composed of 13 PCGs, 2 rRNA, 22 tRNA genes and an AT control region. The arrangement and structure of the mitochondrial genomes were similar to those of other invertebrates. Comparative genomics revealed that the three mitochondrial genomes were highly conserved in terms of gene size, structure, and codon usage, all PCGs were purified selections with an ATN start codon and a TAN stop codon. All tRNAs could be folded into the typical clover-leaf structure, except tRNA Ser (AGN) that lacked a dihydrouridine (DHU) arm. Phylogenetic analysis based on 13 PCGs of 34 Acrididae species and seven outgroup species revealed that differences in the shape of antennae within the family Acrididae should be given less weight as a taxonomic character for higher-level classification. Moreover, the divergence time estimates indicates that in Gomphocerinae, the species with clubbed antennae were formed within the nearest 18 Mya, and Pacris xizangensis is more ancient.

An insect chemosensory protein facilitates locust avoidance to fungal pathogens via recognition of fungal volatiles.

Locusts (Locusta migratoria) are one of the most destructive insect pests worldwide. Entomopathogenic fungi can infect and kill locusts, with Metarhizium acridum having evolved as a specialized acridid pathogen. However, locusts have evolved countermeasures to limit or avoid microbial pathogens, although the underlying molecular mechanisms behind these defenses remain obscure. Here, we demonstrate that L. migratoria exhibit avoidance behaviors towards M. acridum contaminated food via recognition of fungal volatiles, with locust perception of the volatile mediated by the LmigCSP60 chemosensory protein. RNAi-knockdown of LmigCSP60 lowered locust M. acridum avoidance behavior and increased infection and mortality. The fungal volatile, 2-phenylethanol (PEA), was identified to participate in locust behavioral avoidance. RNAi-knockdown of LmigCSP60 reduced antennal electrophysiological responses to PEA and impaired locust avoidance to the compound. Purified LmigCSP60 was able to bind a set of fungal volatiles including PEA. Furthermore, reduction of PEA emission by M. acridum via construction of a targeted gene knockout mutant of the alcohol dehydrogenase gene (ΔMaAdh strain) that contributes to PEA production reduced locust avoidance behavior towards the pathogen. These findings identify an olfactory circuit used by locusts to detect and avoid potential microbial pathogens before they are capable of initiating infection and highlight behavioral and olfactory adaptations affecting the co-evolution of host-pathogen interactions.

Phenotypic variation and genomic variation in insect virulence traits reveal patterns of intraspecific diversity in a locust-specific fungal pathogen.

Intraspecific pathogen diversity is crucial for understanding the evolution and maintenance of adaptation in host-pathogen interactions. Traits associated with virulence are often a significant source of variation directly impacted by local selection pressures. The specialist fungal entomopathogen, Metarhizium acridum, has been widely implemented as a biological control agent of locust pests in tropical regions of the world. However, few studies have accounted for natural intraspecific phenotypic and genetic variation. Here, we examine the diversity of nine isolates of M. acridum spanning the known geographic distribution, in terms of (1) virulence towards two locust species, (2) growth rates on three diverse nutrient sources, and (3) comparative genomics to uncover genomic variability. Significant variability in patterns of virulence and growth was shown among the isolates, suggesting intraspecific ecological specialization. Different patterns of virulence were shown between the two locust species, indicative of potential host preference. Additionally, a high level of diversity among M. acridum isolates was observed, revealing increased variation in subtilisin-like proteases from the Pr1 family. These results culminate in the first in-depth analysis regarding multiple facets of natural variation in M. acridum, offering opportunities to understand critical evolutionary drivers of intraspecific diversity in pathogens.

A chemical defense deters cannibalism in migratory locusts.

Many animals engage in cannibalism to supplement their diets. Among dense populations of migratory locusts, cannibalism is prevalent. We show that under crowded conditions, locusts produce an anticannibalistic pheromone called phenylacetonitrile. Both the degree of cannibalism and the production of phenylacetonitrile are density dependent and covary. We identified the olfactory receptor that detects phenylacetonitrile and used genome editing to make this receptor nonfunctional, thereby abolishing the negative behavioral response. We also inactivated the gene underlying phenylacetonitrile production and show that locusts that lack this compound lose its protection and are more frequently exposed to intraspecific predation. Thus, we reveal an anticannibalistic feature built on a specifically produced odor. The system is very likely to be of major importance in locust population ecology, and our results might therefore provide opportunities in locust management.

PLIN2-induced ectopic lipid accumulation promotes muscle ageing in gregarious locusts.

Ageing plasticity represents the flexibility of the ageing process in response to non-genetic factors, occurring commonly in animals. However, the regulatory mechanisms underlying ageing plasticity are largely unclear. The density-dependent polyphenism of locusts, Locusta migratoria, displays dramatic lifespan divergence between solitary and gregarious phases, providing a useful system for studying ageing plasticity. Here, we found that gregarious locusts displayed faster locomotor deficits and increased muscle degeneration on ageing than solitary locusts. Comparative transcriptome analysis in flight muscles revealed significant differences in transcriptional patterns on ageing between two phases. RNA interference screening showed that the knockdown of the upregulated PLIN2 gene significantly relieved the ageing-related flight impairments in gregarious locusts. Mechanistically, the gradual upregulation of PLIN2 could induce the accumulation of ectopic lipid droplets and triacylglycerols in flight muscles during the ageing process. Further experiments suggested that ectopic lipid accumulation led to an ageing-related ß-oxidation decline through limiting fatty acid transport and content. These findings reveal the key roles of lipid metabolism in the differences of muscle ageing between solitary and gregarious locusts and provide a potential mechanism underlying environment-induced muscle ageing plasticity.

Transcriptome of the pygmy grasshopper Formosatettix qinlingensis (Orthoptera: Tetrigidae).

Formosatettix qinlingensis (Zheng, 1982) is a tiny grasshopper endemic to Qinling in China. For further study of its transcriptomic features, we obtained RNA-Seq data by Illumina HiSeq X Ten sequencing platform. Firstly, transcriptomic analysis showed that transcriptome read numbers of two female and one male samples were 25,043,314, 24,429,905, and 25,034,457, respectively. We assembled 65,977 unigenes, their average length was 1,072.09 bp, and the length of N50 was 2,031 bp. The average lengths of F. qinlingensis female and male unigenes were 911.30 bp, and 941.82 bp, and the N50 lengths were 1,745 bp and 1,735 bp, respectively. Eight databases were used to annotate the functions of unigenes, and 23,268 functional unigenes were obtained. Besides, we also studied the body color, immunity and insecticide resistance of F. qinlingensis. Thirty-nine pigment-related genes were annotated. Some immunity genes and signaling pathways were found, such as JAK-STAT and Toll-LIKE receptor signaling pathways. There are also some insecticide resistance genes and signal pathways, like nAChR, GST and DDT. Further, some of these genes were differentially expressed in female and male samples, including pigment, immunity and insecticide resistance. The transcriptomic study of F. qinlingensis will provide data reference for gene prediction and molecular expression study of other Tetrigidae species in the future. Differential genetic screening of males and females provides a basis for studying sex and immune balance in insects.

Long-term biogeographical processes dominate patterns of genetic diversity in a wingless grasshopper despite substantial recent habitat fragmentation.

Low-vagility species may hold strong genetic signatures of past biogeographical processes but are also vulnerable to habitat loss. Flightless grasshoppers of the morabine group were once widespread in southeastern Australia, including Tasmania, but are becoming restricted to remnant patches of vegetation, with local ranges impacted by agriculture and development as well as management. Habitat fragmentation can generate genetically differentiated "island" populations with low genetic variation. However, following revegetation, populations could be re-established, and gene flow increased. Here we characterize single nucleotide polymorphism-based genetic variation in a widespread chromosomal race of the morabine Vandiemenella viatica (race 19) to investigate the genetic health of remnant populations and to provide guidelines for restoration efforts. We update the distribution of this race to new sites in Victoria and Tasmania, and show that V. viatica populations from northern Tasmania and eastern Victoria have reduced genetic variation compared to other mainland populations. In contrast, there was no effect of habitat fragment size on genetic variation. Tasmanian V. viatica populations fell into two groups, one connected genetically to eastern Victoria and the other connected to southwestern Victoria. Mainland populations showed isolation by distance. These patterns are consistent with expectations from past biogeographical processes rather than local recent population fragmentation and emphasize the importance of small local reserves in preserving genetic variation. The study highlights how genomic analyses can combine information on genetic variability and population structure to identify biogeographical patterns within a species, which in turn can inform decisions on potential source populations for translocations.

New estimates of genome size in Orthoptera and their evolutionary implications.

Animal genomes vary widely in size, and much of their architecture and content remains poorly understood. Even among related groups, such as orders of insects, genomes may vary in size by orders of magnitude-for reasons unknown. The largest known insect genomes were repeatedly found in Orthoptera, e.g., Podisma pedestris (1C = 16.93 pg), Stethophyma grossum (1C = 18.48 pg) and Bryodemella holdereri (1C = 18.64 pg). While all these species belong to the suborder of Caelifera, the ensiferan Deracantha onos (1C = 19.60 pg) was recently found to have the largest genome. Here, we present new genome size estimates of 50 further species of Ensifera (superfamilies Gryllidea, Tettigoniidea) and Caelifera (Acrididae, Tetrigidae) based on flow cytometric measurements. We found that Bryodemella tuberculata (Caelifera: Acrididae) has the so far largest measured genome of all insects with 1C = 21.96 pg (21.48 gBp). Species of Orthoptera with 2n = 16 and 2n = 22 chromosomes have significantly larger genomes than species with other chromosome counts. Gryllidea genomes vary between 1C = 0.95 and 2.88 pg, and Tetrigidae between 1C = 2.18 and 2.41, while the genomes of all other studied Orthoptera range in size from 1C = 1.37 to 21.96 pg. Reconstructing ancestral genome sizes based on a phylogenetic tree of mitochondrial genomic data, we found genome size values of >15.84 pg only for the nodes of Bryodemella holdereri / B. tuberculata and Chrysochraon dispar / Euthystira brachyptera. The predicted values of ancestral genome sizes are 6.19 pg for Orthoptera, 5.37 pg for Ensifera, and 7.28 pg for Caelifera. The reasons for the large genomes in Orthoptera remain largely unknown, but a duplication or polyploidization seems unlikely as chromosome numbers do not differ much. Sequence-based genomic studies may shed light on the underlying evolutionary mechanisms.

Multiple introgression events during the diversification history of the edible Mexican grasshopper genus Sphenarium (Orthoptera: Pyrgomorphidae).

Speciation with gene flow often leads to ambiguous phylogenetic reconstructions, reticulate patterns of relatedness and conflicting nuclear versus mitochondrial (mt) lineages. Here we employed a fragment of the COI mtDNA gene and nuclear genome-wide data (3RAD) to assess the diversification history of Sphenarium, an orthopteran genus of great economic importance in Mexico that is presumed to have experienced hybridisation events in some of its species. We carried out separate phylogenetic analyses to evaluate the existence of mito-nuclear discordance in the species relationships, and also assessed the genomic diversity and population genomic structure and investigated the existence of interspecific introgression and species limits of the taxa involved based on the nuclear dataset. The species delineation analyses discriminated all the currently recognised species, but also supported the existence of four undescribed species. The mt and nuclear topologies had four discordant species relationships that can be explained by mt introgression, where the mt haplotypes of S. purpurascens appear to have replaced those of S. purpurascens A and B, S. variabile and S. zapotecum. Moreover, our analyses supported the existence of nuclear introgression events between four species pairs that are distributed in the Sierra Madre del Sur province in southeast Mexico, with three of them occurring in the Tehuantepec Isthmus region. Our study highlights the relevance of genomic data to address the relative importance of allopatric isolation versus gene flow in speciation.

Tandem Repeat DNA Provides Many Cytological Markers for Hybrid Zone Analysis in Two Subspecies of the Grasshopper Chorthippus parallelus.

Recent advances in next generation sequencing (NGS) have greatly increased our understanding of non-coding tandem repeat (TR) DNA. Here we show how TR DNA can be useful for the study of hybrid zones (HZ), as it serves as a marker to identify introgression in areas where two biological entities come in contact. We used Illumina libraries to analyse two subspecies of the grasshopper Chorthippus parallelus, which currently form a HZ in the Pyrenees. We retrieved a total of 152 TR sequences, and used fluorescent in situ hybridization (FISH) to map 77 families in purebred individuals from both subspecies. Our analysis revealed 50 TR families that could serve as markers for analysis of this HZ, using FISH. Differential TR bands were unevenly distributed between chromosomes and subspecies. Some of these TR families yielded FISH bands in only one of the subspecies, suggesting the amplification of these TR families after the geographic separation of the subspecies in the Pleistocene. Our cytological analysis of two TR markers along a transect of the Pyrenean hybrid zone showed asymmetrical introgression of one subspecies into the other, consistent with previous findings using other markers. These results demonstrate the reliability of TR-band markers for hybrid zone studies.

Mwhitiwhiti Aotearoa: Phylogeny and synonymy of the silent alpine grasshopper radiation of New Zealand (Orthoptera: Acrididae).

Aotearoa New Zealand has a fauna of endemic alpine grasshoppers, consisting of thirteen species distributed among four genera. The many re-classifications of species within this group and the presence of species complexes highlight the uncertainty that surrounds relationships within and between these genera. High-throughput Next Generation Sequencing was used to assemble the complete mitochondrial genomes, 45S ribosomal cassettes and histone sequences of New Zealands four endemic alpine genera: Alpinacris, Brachaspis, Paprides and Sigaus. Phylogenetic analysis of these molecular datasets, as individual genes, partitions and combinations returned a consistent topology that is incompatible with the current classification. The genera Sigaus, Alpinacris, and Paprides all exhibit paraphyly. A consideration of the pronotum, epiphallus and terminalia of adult specimens reveals species-specific differences, but fails to provide compelling evidence for species groups justifying distinct genera. In combination with phylogenetic, morphological and spatial evidence we propose a simplified taxonomy consisting of a single genus for the mwhitiwhiti Aotearoa species radiation.

Identification and analysis of odorant receptors expressed in the two main olfactory organs, antennae and palps, of Schistocerca americana.

Locusts depend upon their sense of smell and provide useful models for understanding olfaction. Extending this understanding requires knowledge of the molecular and structural organization of the olfactory system. Odor sensing begins with olfactory receptor neurons (ORNs), which express odorant receptors (ORs). In insects, ORNs are housed, in varying numbers, in olfactory sensilla. Because the organization of ORs within sensilla affects their function, it is essential to identify the ORs they contain. Here, using RNA sequencing, we identified 179 putative ORs in the transcriptomes of the two main olfactory organs, antenna and palp, of the locust Schistocerca americana. Quantitative expression analysis showed most putative ORs (140) are expressed in antennae while only 31 are in the palps. Further, our analysis identified one OR detected only in the palps and seven ORs that are expressed differentially by sex. An in situ analysis of OR expression suggested ORs are organized in non-random combinations within antennal sensilla. A phylogenetic comparison of OR predicted protein sequences revealed homologous relationships among two other Acrididae species. Our results provide a foundation for understanding the organization of the first stage of the olfactory system in S. americana, a well-studied model for olfactory processing.

Transposable element expansion and low-level piRNA silencing in grasshoppers may cause genome gigantism.

BACKGROUND: Transposable elements (TEs) have been likened to parasites in the genome that reproduce and move ceaselessly in the host, continuously enlarging the host genome. However, the Piwi-interacting RNA (piRNA) pathway defends animal genomes against the harmful consequences of TE invasion by imposing small-RNA-mediated silencing. Here we compare the TE activity of two grasshopper species with different genome sizes in Acrididae (Locusta migratoria manilensis♀1C = 6.60 pg, Angaracris rhodopa♀1C = 16.36 pg) to ascertain the influence of piRNAs. RESULTS: We discovered that repetitive sequences accounted for 74.56% of the genome in A. rhodopa, more than 56.83% in L. migratoria, and the large-genome grasshopper contained a higher TEs proportions. The comparative analysis revealed that 41 TEs (copy number > 500) were shared in both species. The two species exhibited distinct "landscapes" of TE divergence. The TEs outbreaks in the small-genome grasshopper occurred at more ancient times, while the large-genome grasshopper maintains active transposition events in the recent past. Evolutionary history studies on TEs suggest that TEs may be subject to different dynamics and resistances in these two species. We found that TE transcript abundance was higher in the large-genome grasshopper and the TE-derived piRNAs abundance was lower than in the small-genome grasshopper. In addition, we found that the piRNA methylase HENMT, which is underexpressed in the large-genome grasshopper, impedes the piRNA silencing to a lower level. CONCLUSIONS: Our study revealed that the abundance of piRNAs is lower in the gigantic genome grasshopper than in the small genome grasshopper. In addition, the key gene HENMT in the piRNA biogenesis pathway (Ping-Pong cycle) in the gigantic genome grasshopper is underexpressed. We hypothesize that low-level piRNA silencing unbalances the original positive correlation between TEs and piRNAs, and triggers TEs to proliferate out of control, which may be one of the reasons for the gigantism of grasshopper genomes.

Sexual repurposing of juvenile aposematism in locusts.

Adaptive plasticity requires an integrated suite of functional responses to environmental variation, which can include social communication across life stages. Desert locusts (Schistocerca gregaria) exhibit an extreme example of phenotypic plasticity called phase polyphenism, in which a suite of behavioral and morphological traits differ according to local population density. Male and female juveniles developing at low population densities exhibit green- or sand-colored background-matching camouflage, while at high densities they show contrasting yellow and black aposematic patterning that deters predators. The predominant background colors of these phenotypes (green/sand/yellow) all depend on expression of the carotenoid-binding "Yellow Protein" (YP). Gregarious (high-density) adults of both sexes are initially pinkish, before a YP-mediated yellowing reoccurs upon sexual maturation. Yellow color is especially prominent in gregarious males, but the reason for this difference has been unknown since phase polyphenism was first described in 1921. Here, we use RNA interference to show that gregarious male yellowing acts as an intrasexual warning signal, which forms a multimodal signal with the antiaphrodisiac pheromone phenylacetonitrile (PAN) to prevent mistaken sexual harassment from other males during scramble mating in a swarm. Socially mediated reexpression of YP thus adaptively repurposes a juvenile signal that deters predators into an adult signal that deters undesirable mates. These findings reveal a previously underappreciated sexual dimension to locust phase polyphenism, and promote locusts as a model for investigating the relative contributions of natural versus sexual selection in the evolution of phenotypic plasticity.