Introduction

Reproductive behavior is critical for population sustenance and survival of species, and is influenced by a complex array of internal conditions, environmental factors, and social interactions. In insects, males typically initiate mating through courtship, while females determine acceptance by integrating external sensory stimuli with their internal state. Consequently, female reproductive activity is modulated by multiple signaling pathways ^1–8. In insects, successful copulation commonly induces profound physiological and behavioral changes in the females ^5,9–18.

Following their first mating, insect females typically reject further courtship attempts and initiate oviposition. This distinct shift in behavior and physiology is referred to as a post-mating response (PMR) ^{5–7,19–23}. A PMR has been documented across diverse insect taxa, including Drosophila melanogaster, Anopheles gambiae, Aedes aegypti and the brown planthopper (BPH), Nilaparvata lugens ^{6,18,21,24,25}.

Induction of a PMR commonly requires the transfer of specific seminal fluid components during copulation. In Drosophila, seminal fluid peptides, including sex peptide (SP) and DUP99B produced in male accessory glands (MAG), enter the female reproductive tract, bind to specific receptors on sensory neurons, and activate signaling pathways that drive post-mating behavior and physiology ^{6,13,14,18,26}. Additional factors such as Acp26Ab and ovulin (Acp26Aa) are also transferred with the seminal fluid and modulate female reproductive physiology in Drosophila ^27–30. Notably, it has been shown that the PMR-inducing substances, produced in the MAG of various insects, are highly species-specific and not even conserved across related taxonomic groups ^25,31.

While MAG-derived sex peptide (SP) and its receptor (SPR) are key regulators of the female post-mating response (PMR) in Drosophila ^5–7,16, additional signaling pathways intrinsic to females also modulate this process. Interneurons within abdominal neuromeres of the ventral nerve cord (VNC) promote receptivity to males by signaling via myoinhibitory peptide (MIP) to brain circuits ³². Furthermore, the neuropeptide diuretic hormone 44 (DH44) modulates female sexual receptivity through a sex-specific brain circuit; mating suppresses DH44 signaling, thereby reducing receptivity ¹. During sexual maturation, leucokinin (LK)-expressing neurons (ABLKs) in the abdominal ganglion, downstream of SP-activated circuitry, suppress female receptivity ³³. In addition to changes in receptivity, the Drosophila PMR includes a shift in diet preference toward yeast consumption. This altered feeding behavior is mediated by specific peptidergic neurons (ALKs) in the brain through LK signaling ¹².

Although the mechanisms and neuronal pathways underlying the post-mating response (PMR) in Drosophila are well characterized ^{3,4,6,7,34–37}, this process remains poorly understood in most non-model insects. A prominent example is the brown planthopper (BPH), N. lugens (Hemiptera). As a monophagous rice pest, the BPH exhibits robust reproductive capacity and remarkable environmental adaptability ³⁸. It has also developed high resistance to diverse chemical pesticides ^39–41. Consequently, elucidating reproductive physiology and behavior in BPHs is crucial for developing novel control strategies.

In BPHs, the recently identified MAG-derived peptide maccessin (macc) is transferred during copulation and triggers a PMR in females, although its receptor remains unknown ²⁵. Notably, the same study demonstrated that a second peptide, ion transport-like peptide 1 (ITPL-1), also modulates the PMR. Whereas macc is male-specific, ITPL-1 is produced both in the MAG and endogenously in females, where it reduces female receptivity to courting males ²⁵. Given its endogenous female expression, ITPL-1 likely functions as an intrinsic signaling component in the PMR of BPHs, analogous to the role of myoinhibitory peptide (MIP) and DH44 in Drosophila females ^1,32.

Here, we explored the role of another neuropeptide, corazonin (CRZ), in the reproductive behavior of BPHs. In Drosophila males, CRZ is released from interneurons of the ventral nerve cord (VNC) and acts on serotonergic projection neurons to regulate ejaculation and copulation duration ⁴². Additionally, CRZ release in males during mating mitigates mating-induced heart rate acceleration, to promote physiological recovery ⁴³. CRZ furthermore influences post-mating dietary preference in males ¹². Surprisingly, the role of CRZ signaling in female reproduction remains unexplored. To fill this gap, we examined the role of CRZ in female reproduction, including the PMR, in BPH and Drosophila melanogaster.

In this study, we demonstrate that CRZ induces a post-mating switch in receptivity in females of both BPH and Drosophila. Notably, CRZ injection in virgin females diminishes receptivity to courting males, while Crz knockdown or knockout in mated females renders them receptive to further mating. Importantly, CRZ is not produced in the MAG, and thus not transferred from males, but acts as an endogenously produced peptide in females. Furthermore, we detected strong expression of the CRZ receptor (CrzR) in the female reproductive tract of BPHs. We demonstrated that both CRISPR-Cas9-mediated CrzR knockout and RNAi-mediated CrzR knockdown disrupt the PMR, phenocopying Crz knockdown. Critically, we could show that the level of receptivity to mating in virgin females, whether by MAG extract (with seminal fluid proteins) injection or macc injection, is strictly dependent on presence of the CrzR receptor. Importantly, manipulations of CRZ signaling in female Drosophila yield similar effects on receptivity and oviposition. Interestingly, we find that knockdown or knockout of Crz in male N. lugens does not affect male copulation behavior, but Crz deficient males do not induce a full PMR in mated females. Notably, knockout of Crz in both male and female partners results in an almost complete abolishment of the PMR. Since CRZ is not produced in the MAG we performed a transcriptomic analysis of the MAG in Crz mutants to screen for other candidate CRZ-regulated genes affecting the PMR. We found altered expression of a large number of genes, containing seminal protein genes, but no change in the genes encoding macc and ITPL-1, suggesting presence of other MAG factors, yet to be identified. In summary, our findings strengthen the notion that the female PMR requires not only transfer of male-derived factors, but also endogenous female neuropeptide signaling components.

Results

CRZ signaling reduces receptivity to courting males and stimulates oviposition in N. lugens females

While numerous studies have established roles of the neuropeptide CRZ in regulating aspects of male reproduction across diverse insect species ^44–50, research addressing its potential functions in female reproduction remains conspicuously lacking.

This study investigates the role of CRZ in female reproduction, starting with the brown planthopper (BPH), Nilaparvata lugens, a major agricultural pest. Initially, the Crz gene (GenBank accession: AB817247.1) was cloned based on transcriptome data ⁵¹. The deduced mature CRZ sequence in BPH, pQTFQYSRGWTNamide, is identical to that reported for most studied insect species (Figure 1-figure supplement 1A). Consistent with other known CRZ precursors, the BPH precursor encodes a single CRZ peptide starting at the first amino acid of the propeptide (Figure 1-figure supplement 1A and B). A comparison of the organization of selected CRZ precursor genes is shown in Figure 1-figure supplement 1B. Genomic analysis revealed that the CRZ precursor is encoded by two exons separated by a single intron (Figure 1-figure supplement 1B). The mature CRZ peptide and a scrambled control peptide (sCRZ) were synthesized for pharmacological experiments.

First, we investigated the role of CRZ signaling in the reproductive behavior and PMR of female BPHs. The PMR in BPHs, is characterized by reduced receptivity to courting males, including display of rejection behaviors, and increased oviposition ²⁵. This response is partially induced by the male-derived peptide maccessin (macc), transferred during copulation ²⁵. We found that injection of CRZ significantly reduced receptivity of virgin female BPHs to courting males (Figure 1A and B). Furthermore, virgin females injected with CRZ actively rejected courting males, exhibiting a distinct rejection behavior characterized by ovipositor extrusion and kicking or fleeing. This behavior was not observed in control females injected with sCRZ (Figure 1 - Supplementary Video 1). Notably, the CRZ-induced mating refusal behavior in virgin females was no longer observed 24 hours after injection (Figure 1-figure supplement 2A). While injection of CRZ peptide did not stimulate egg production in virgin BPHs (Figure 1-figure supplement 2B), it significantly promoted oviposition in mated females (Figure 1C).

Figure 1.

Figure 1 figure supplement 1.

Figure 1 figure supplement 2.

Figure 1 figure supplement 3.

Figure 1 figure supplement 4.

To further investigate CRZ action, we performed RNA interference (RNAi) by injecting double-stranded Crz RNA (dsNlcrz) into female BPHs, which efficiently reduced Crz expression (Figure 1-figure supplement 2). Knockdown of Crz resulted in increased receptivity with 10% of the mated females mating again, and significantly reduced oviposition to approximately 60% of control levels (Figure 1E and F). In contrast, Crz knockdown had no effect on the receptivity of virgin females (Figure 1E).

Since RNAi-mediated knockdown of Crz in female BPHs left a residual level of gene expression (Figure 1-figure supplement 2), we generated two Crz deletion alleles, ΔCrz¹ and ΔCrz², using CRISPR/Cas9 (Figure 1G-G’) to achieve a complete ablation. Both alleles harbor deletions eliminating the mature CRZ peptide coding sequence (Figure 1G’’). Sanger sequencing confirmed that the deletions compromised essential triplet codons within the Crz open reading frame (Figure 1-figure supplement 4), thereby efficiently disrupting CRZ peptide expression (Figure 1G’’ and H). Immunohistochemical analysis further validated the mutants, revealing a complete absence of anti-CRZ immunolabeling in homozygous ΔCrz¹ and ΔCrz²females (Figure 1H). Homozygous ΔCrz¹ and ΔCrz² virgin females are fully viable and fertile. They exhibit a courtship vigor similar to controls and their initial mating success rate with wild-type males exceeds 60% (Figure 1I and J). However, 25% of the mated mutant females were receptive to further mating (Figure 1J), and oviposition was significantly reduced to approximately 50% of the wild-type levels (Figure 1K). Collectively, these results demonstrate that diminished CRZ signaling severely perturbs the PMR in female BPHs, seen as increased post-mating receptivity and reduced egg-laying. Importantly, the Crz deletion abolished the onset of typical post-mating rejection behavior. We, thus, conclude that endogenous CRZ signaling is important for the induction of a PMR in female BPHs.

CRZ is expressed in similar neurons in both sexes of BPHs and is absent from the male accessory gland

To further understand how CRZ regulates the PMR in BPHs, we examined the expression of CRZ and Crz across different tissues in males and females. The Crz mRNA expression was examined in various BPH tissues using both RT-PCR and qPCR (Figure 2A). We detected Crz mRNA expression primarily in the central nervous system (CNS), regardless of the primer set or detection method used, while expression was nearly undetectable in other tissues, including reproductive organs of both sexes (Figure 2B and C).

Figure 2.

Figure 2 figure supplement 1.

A previous study reported CRZ expression in four male-specific interneurons within the abdominal ganglia of Drosophila ⁴⁴. To determine whether CRZ exhibits sexually dimorphic expression in BPHs, we performed immunolabeling using an anti-CRZ antiserum. This labeled four pairs of neurons in the protocerebrum of the brain (Figure 2D and E) and one pair of neurons in the subesophageal ganglion in both males and females (Figure 2D and E). Since CRZ immunolabeled neurons in brains of most, if not all, insects studied include three or more pairs of lateral neurosecretory cells ^52,53, we suggest that CRZ signaling regulating the PMR may be hormonal (at least in part). Importantly, no CRZ-immunopositive neurons were detected in the abdominal ganglia of either sex (Figure 2D and E), and thus there is no evidence for sex-specific expression of CRZ neurons in BPHs. This is in contrast to Drosophila, where we confirmed, as previously reported ⁴², the presence of CRZ-immunolabeled neurons specifically in the abdominal ganglia of males (Figure 2 - figure supplement 1A).

To exclude that CRZ acts as a seminal fluid peptide transferred to females during copulation, similar to sex peptide (SP) in Drosophila and macc in BPHs ²⁵, we specifically examined CRZ expression in reproductive tissues of both sexes. However, we observed no CRZ immunolabeling in the male reproductive organs (Figure 2 - figure supplement 1B1 - B3) or female reproductive organs (Figure 2 - figure supplement 1B4 – B6). This aligns with our transcriptional Crz data shown above (Figure 2B and C). Furthermore, Crz mRNA and CRZ protein were absent from transcriptomic and proteomic datasets of the male accessory glands (MAGs) in BPHs ²⁵. Collectively the available data suggest, that CRZ is not produced by the MAG and is consequently unlikely to be transferred to females as a seminal fluid component.

The CRZ receptor is essential for mediating the post mating response in female BPHs

Given our findings implicating CRZ in modulating the PMR in female BPHs, we identified the CRZ receptor (CrzR) gene in this species; it encodes a structurally conserved GPCR orthologous to other insect CrzRs (Figure 3-supplemental figure 1). We next explored the effect of knocking down the CrzR on the PMR in BPHs. RNAi-mediated knockdown of the receptor (dsCrzR injection) resulted in 15% of the mated females re-mating, a phenotype not seen in dsgfp-injected controls (Figure 3A and B) and in a significantly reduced oviposition (Figure 3C). Importantly, CRZ injection failed to diminish receptivity in CrzR-knockdown virgins (Figure 3D and E). Furthermore, MAG extract (with seminal fluid proteins) injection into dsCrzR-treated virgins did not affect receptivity (Figure 3F).

Figure 3.

Figure 3 figure supplement 1.

Figure 3 figure supplement 2.

To further confirm the role of the CrzR in the PMR of BPHs, we applied CRISPR/Cas9 genome editing to mutate the CrzR gene and constructed a homozygous mutant strain of CrzR (CrzR^M) by genetic crosses (Figure 3G-3G’’, Fig. 3 - figure supplement 2A). We found that the developmental time and life span of the female CrzR mutants, like the Crz mutants, were slightly extended compared to wild type animals (Figure 3 - supplement figure 2D and 2E). Homozygous CrzR^M females displayed increased remating rates (30%) (Figure 3I and J) and the number of eggs laid by CrzR^M females was also significantly decreased (Figure 3K). Consistent with the RNAi results, neither CRZ nor MAG extract injections reduced receptivity in CrzR^Mvirgins (Figure 3K - M). Finally, macc peptide injection also failed to decrease receptivity in CrzR-mutant virgins (70% acceptance rate), whereas it significantly suppressed mating in wild-type controls (Figure 3N). These results collectively demonstrate that intact CRZ/CrzR signaling in female BPHs is required for transducing the PMR-inducing effects of male-derived seminal factors, including macc ²⁵.

The CrzR is highly expressed in the reproductive organs of female BPHs

To further elucidate the role of the CRZ signaling in regulation of the PMR, we analyzed the CrzR expression in N. lugens tissues. We first monitored CrzR expression by means of RT-PCR and qPCR. Both methods revealed predominant CrzR transcript expression in the female reproductive tract (Figure 4A-C). qPCR further demonstrated significantly higher CrzR expression in adult females than in males (Figure 4 - figure supplement 1A), particularly localized to female abdomens (Figure 4 - figure supplement 1B).

Figure 4.

Figure 4 figure supplement 1.

As we failed to generate a functional NlCrzR antibody, to determine the tissue expression of the receptor in more detail, we used in situ hybridization and confirmed strong CrzR expression within the female reproductive system (Figure 4D). No hybridization signal was detected in the lateral or common oviducts (Figure 4D1). However, specific labeling occurred in the sperm storage organs: the spermathecae and pouched glands of the lower reproductive tract (Figure 4D2 and 4D2’). These structures are the primary sites for long-term sperm storage. Specificity was validated by absence of signal in sense-probe controls (Figure 4D3 - D4’). This enrichment in sperm storage organs suggests that CRZ signaling may modulate post-mating physiology in part by regulating the mobilization of stored sperm and seminal fluid.

Notably, while CrzR is highly expressed in the female reproductive tract, we acknowledge that RT-PCR and qPCR analyses also detected low but detectable CrzR expression in other tissues, including the CNS and fat body (Figure 4B-C). Thus, we cannot exclude the possibility that CrzR in the brain (e.g., by modulating neural circuits governing reproductive behavior) or fat body (e.g., by affecting post-mating metabolism and nutrient/energy allocation) may also contribute to the PMR. These non-reproductive tract sites of CrzR expression warrant further investigation to fully delineate the systemic roles of CRZ/CrzR signaling in female PMR regulation.

CRZ signaling in male BPHs impacts the male accessory glands and secondarily modulates the female PMR

Although our primary focus is on CRZ signaling in the female PMR, we also investigated its role in male BPHs, given its established importance in male reproductive biology in Drosophila ^44,46,48. In Drosophila, silencing CRZ-expressing neurons in the male abdominal ganglia disrupts sperm/seminal fluid transfer and prolongs copulation, while their activation, or CRZ peptide injection, induces rapid ejaculation ⁴⁴. To test whether this function is conserved in BPHs, we first injected synthetic CRZ peptide into male BPHs. In contrast to Drosophila, CRZ injection did not trigger precocious ejaculation in male BPHs (Figure 5-figure supplement 1A) nor alter the morphology of the male accessory gland (MAG) (Figure 5-figure supplement 1B). However, CRZ-injected males exhibited a significantly reduced mating success rate (defined as the proportion of all tested individuals that successfully court and complete copulation), compared to controls injected with scrambled CRZ (sCRZ) (Figure 5A and B), but maintained unaltered courtship rates (courtship drive) and copulation durations (Figure 5-figure supplement 2A and B).

Figure 5.

Figure 5 figure supplement 1:

Figure 5 figure supplement 2.

Since the male mating success rate in some cases was affected by CRZ, we made sure that the female PMR defects do not stem from incomplete male mating: all females included in our PMR analyses were mated with males that had completed full, successful copulation (i.e., completed all mating steps, including seminal fluid transfer, as verified by behavioral observation and subsequent validation of sperm presence in female reproductive tracts). Strikingly, even with fully successful mating, females mated with CRZ-injected males displayed a severely impaired PMR: nearly 30% of these females re-mated within the test window (vs. 0% of females mated with sCRZ-injected controls) (Figure 5B), and their oviposition was significantly reduced (Figure 5C). Although CRZ injection did not quantitatively change the total protein content of the MAG (Figure 5-figure supplement 1C), we hypothesize that it alters the composition of seminal fluid proteins, consistent with findings in the oriental fruit moth Grapholita molesta ⁴⁹ and thereby reducing the ability to induce a robust female PMR. In male BPHs, RNAi-mediated knockdown of the Crz transcript (via dsNlCrz injection) had no significant effect on key male reproductive behavior: mating success rate (proportion of males that complete copulation), courtship rate (proportion of males exhibiting courtship behaviors such as following or wing display), and copulation duration (Figure 5D and E, Figure 5—figure supplement 2C and D). However, females mated with dsNlCrz-injected males exhibited a defective PMR: about 10% remated (Figure 5E) and laid approximately 60% fewer eggs (Figure 5F). Consistent with the RNAi results, males carrying Crz mutant alleles (ΔCrz¹ or ΔCrz²) also failed to elicit a full PMR in mated females, despite displaying normal courtship behavior and mating success (Figure 5—figure supplement 2E and F). Specifically, ∼25% of the females mated with Crz mutant males accepted re-mating (Figure 5G and H), and their egg production was halved compared to females mated with wild-type males (Figure 5I).

Notably, when Crz mutant males were crossed with Crz mutant females, the PMR defects were dramatically exacerbated: the re-mating rate of females surged to nearly 50%, and egg production plummeted by ∼75% relative to wild-type male-female pairs (Figure 5J - L). This near-complete ablation of the PMR underscores that intact CRZ signaling in both sexes is indispensable for triggering a robust PMR in female BPHs.

CRZ signaling in male BPHs alters seminal fluid composition

The observation that manipulating Crz (via knockdown or mutation) in male BPHs affects the female PMR is particularly intriguing, prompting us to investigate the underlying mechanism. Since we had ruled out production of CRZ in the male reproductive system, including the male accessory gland (MAG), we set out to test whether males with impaired CRZ signaling exhibit defects in basic reproductive performance (e.g., mating success, sperm production/transfer).

To eliminate confounding effects of abnormal sperm production or transfer, we first assessed some key parameters: spermathecal size in females mated to wild-type (WT) or ΔCrz¹ males showed no significant differences (Figure 6A), and sperm counts in the testes of unmated ΔCrz¹ males (vs. WT) and in the spermathecae of females mated to ΔCrz¹ males (vs. WT) were comparable (Figure 6B). These data confirm intact sperm production and transfer, pointing to potential alterations in seminal fluid composition as the driver of female PMR defects. This led us to focus on the protein production in the MAG, encouraged by prior work in the oriental fruit moth (G. molesta) where it was shown that CRZ signaling modulates biosynthetic activity in the MAG ⁴⁹.

Figure 6.

To approach this, we performed RNA-sequencing (RNA-seq) on MAGs from WT and ΔCrz¹ males. Comparative transcriptomics identified 2,547 differentially expressed genes (DEGs) in ΔCrz¹ MAGs (|logLJFC| > 1, FDR < 0.05), with 2,114 upregulated and 431 downregulated (Figure 6C). Notably, 16 genes encoding seminal fluid proteins (SFPs) were dysregulated (7 downregulated, 9 upregulated). Gene Ontology (GO) enrichment analysis of downregulated DEGs revealed significant overrepresentation of terms related to serine-type endopeptidase activity (GO: 0004252; Padj < 0.05)—a function critical for SFP processing (Figure 6D).

Consistent with transcriptional dysregulation, the total protein content of ΔCrz¹ male MAGs was significantly lower than that of WT MAGs, as quantified by the Bicinchoninic Acid (BCA) assay (Figure 6E). Targeted qPCR of specific genes in the MAG confirmed significant downregulation of two genes (Figure 6F and G), and also confirmed unchanged transcript levels of the MAG peptides, macc and ITPL-1, in mutants (Figure 6H and I). In summary, the transcriptome analysis indicates that CRZ influences the expression of MAG genes and this results in a defective PMR in mated females. Since macc and ITPL-1 are not among the CRZ-regulated genes, some other gene(s) encoding PMR inducing factors need to be identified in future work.

Endogenous CRZ signaling modulates aspects of the PMR in female

Drosophila melanogaster

Having established the critical role of CRZ signaling in regulating the PMR of the female BPH, we next investigated its potential function in female reproduction, including the PMR, of the genetically tractable fly Drosophila melanogaster.

We first assessed re-mating frequency in Crz mutant female flies (Crz^attp). Approximately 40% of the Crz^attpfemales that had experienced a first successful copulation mated again, compared to only 3% of mated wild-type control females (Figure 7A and B). Importantly, the Crz mutation did not affect the first mating of females (Figure 7B). Furthermore, mated Crz^attp females laid approximately 45% fewer eggs than mated wild-type females (Figure 7C). These results demonstrate that CRZ signaling is necessary for aspects of the PMR in female Drosophila, specifically inhibiting re-mating and stimulating post-mating egg production. Immunohistochemical labeling confirmed the absence of CRZ protein in Crz^attpmutants (Figure 7D).

Figure 7.

To determine whether release of CRZ from Crz neurons is responsible for inhibition of remating and stimulation of oviposition, we knocked down Crz expression specifically in CRZ-producing neurons using the GAL4-UAS system (Crz-GAL4 > UAS-Crz-RNAi). Consistent with the mutant phenotype, 22% of the mated Crz-GAL4 > UAS-Crz-RNAi females mated again after their first successful mating with wild-type males, whereas only 2% of the control females (Crz-GAL4/+) did so (Figure 7E). Virgin receptivity was again unaffected by Crz knockdown (Figure 7E). Additionally, mated Crz-GAL4 > UAS-Crz-RNAi females laid approximately 20% fewer eggs than mated controls (Crz-GAL4/+) (Figure 7F). To further confirm the role of CRZ neuropeptide in modulation of the PMR in females, we injected this peptide in virgin females. Consistent with the above results, CRZ injected virgins displayed a reduced response to mating attempts (Figure 7G). In conclusion, like in BPHs, the neuropeptide CRZ endogenously modulates post-mating behavior also in female Drosophila.

Discussion

Our study unveils a novel function of CRZ in an endogenous female signaling pathway that modulates the post-mating reproductive behavior and physiology (collectively PMR) in BPHs. It is, thus, the first report linking CRZ to a PMR in insects. In contrast to SP and DUP99B in Drosophila ^5,18, HP-1 in Aedes aegypti ¹⁵ and maccessin (macc) recently identified in BPHs ²⁵, CRZ is not produced in the MAG, and thus not transferred from males during copulation. Instead, it is acting as an endogenous peptide derived from neuroendocrine cells in the female brain. Moreover, this CRZ-mediated endogenous modulation in females is seen also in Drosophila. When Crz is knocked down or knocked out in flies, mated females likewise exhibit an altered PMR.

Interestingly, we found that in BPH males, endogenous CRZ acts on the MAG to affect biosynthesis of seminal fluid proteins (SFP), similar to another insect, the lepidopteran Grapholita molesta ⁴⁹. That study revealed that silencing the Crz gene results in a downregulation of multiple genes in the MAG, including several serine-type endopeptidases, however, the authors did not analyze the effect on the female PMR. Our experiments demonstrate that with impaired CRZ signaling in males the seminal fluid appears to lack some component(s) necessary for inducing a PMR since females mated by Crz knockdown males display a reduced PMR. Our gene expression data reveals that neither macc nor ITP-L1 are among the MAG genes affected by Crz knockout. These are encoding peptides that are known to induce a partial PMR in female BPHs ²⁵. Our results therefore suggest that there are additional PMR-inducing factors produced by the MAG of BPHs and that these are regulated by CRZ. Further work is required to identify these factors, but from the present data we can conclude that in BPHs, CRZ signaling influences the female PMR both by action in males and females, but in distinct ways. Similar to our CRZ findings here, a possible dual role of ITPL-1 was also seen in BPHs ²⁵. This peptide is produced in the MAG, but also in tissues of female BPHs, and knockdown in males results in a reduced PMR in females ²⁵.

Focusing on the endogenous CRZ signaling in female BPHs, our main findings using mutants and RNAi of both Crz and CrzR, demonstrate that CRZ signaling decreases the female responsiveness to courting males, especially in already mated females. Also, oviposition in mated BPHs is increased by CRZ signaling. CRZ immunolabeling revealed a small set of neurons in the brain and subesophageal ganglion and no CRZ was detected in reproductive organs. Thus, the likely source of CRZ affecting the PMR is the CRZ-producing neuroendocrine cells in the brain. In insects studied so far (including other hemipterans), there are at least three pairs of CRZ expressing lateral neurosecretory cells with axon terminations in the corpora cardiaca ^52,53, suggesting that in BPHs and Drosophila hormonal CRZ may contribute to the modulation of the PMR. However local CRZ signaling in brain circuits cannot be excluded since the CRZ neurons arborize extensively within the brain.

We found that the CRZ receptor, CrzR, is predominantly expressed in the female reproductive tract in BPHs, more specifically in the lower tract including the spermatheca and pouched gland, which are sites of sperm and seminal fluid storage. Thus, in BPHs, CRZ may act to regulate the release of stored sperm and seminal fluid in the female. This regulation would also gate the flow of male-derived factors that induce and maintain the PMR, such as macc and ITPL-1. As mentioned, we cannot exclude additional action of CRZ in neuronal circuits in the CNS that regulate female sexual arousal and reproductive behavior, since our PCR data indicate expression of the CrzR also in the CNS. Additionally, the fat body of BPHs, express CrzR, similar to Drosophila ⁵⁴, possibly suggesting a role of CRZ in post-mating metabolism.

The brain CRZ neurons of BPHs are likely to be activated by sensory signals from the female reproductive tract upon copulation and hence release CRZ into the circulation (and maybe within the brain). In Drosophila, sensory neurons in the reproductive tract are activated by SP (via the SP receptor), and signal via a chain of ascending axons of the VNC to sex-specific neurons in the brain to induce the PMR switch in behavior and physiology ^7,35,37,55. Thus, in Drosophila, these sex-specific circuits may also activate the CRZ neurons. Possibly similar sensory and sex-specific neurons exist in BPHs and respond either to mechanical stimuli or factors in the seminal fluid leading to activation of CRZ neurons and other brain neurons that regulate the PMR.

The female PMR is complex in Drosophila and not only receptivity and oviposition are affected, but also feeding, metabolism, sleep patterns and aggression ^{5,16,18,19,56,57}. In BPHs only the receptivity (and associated rejection behavior) and oviposition have been analyzed in any detail. We found that the MAG-derived peptides transferred with the semen, macc and ITPL-1, only affected the receptivity ²⁵, whereas we show here that CRZ affects also oviposition. Thus, we have identified three different peptides CRZ, macc and ITPL-1 that all influence the PMR in different ways, suggesting complex regulatory mechanisms also in BPHs.

It is interesting to note that CRZ is an evolutionarily conserved neuropeptide that has been found in most, but not all, insect species and it displays a wide array of functions ^53,54,58,59. The CRZ function varies across different insect species, but a common role seems to be in regulation of stress ^59–61. Other functions are for example: in the American cockroach, CRZ accelerates the heart rate ⁶², in locusts, it triggers the formation of pigments that darken the body color in gregarious insects ⁶³, in the moth Bombyx mori, CRZ affects developmental speed and reduces the rate of silk spinning behavior ⁶⁴; in Drosophila larvae, CRZ regulates molting motor behavior ⁶⁵ and growth ⁶⁶; in adult Drosophila CRZ regulates stress responses, food intake and metabolism ^54,67, and protein food preference ¹², as well as sperm ejection and mating duration in males ⁴²; and in social ants, it regulates the social hierarchy and inhibits egg-laying in the queen ⁴⁵. Finally, a recent study found that CRZ mediates a diapause-induced suppression of reproduction in female bean bugs ⁶⁸. To all these functions can now be added the role of CRZ in reproductive behavior and physiology of BPHs, especially in the female PMR of BPHs and Drosophila. However, we also found a role of CRZ signaling in regulating the production of seminal fluid proteins in the MAG. Finally, building on the role of CRZ in metabolism in Drosophila ⁵⁴, it is possible that CRZ also regulates post-mating metabolism in female BPHs. Thus, there may be a link between CRZ and the altered feeding and metabolism commonly seen in mated female insects.

Key resources table

Method details

Experiment insects

All the brown planthoppers used in this research were collected from Hangzhou, Zhejiang Province in 1995. The insects were kept indoors on fresh rice seedlings (Taichung Native 1, TN1) in a growth chamber at 27 ± 1 °C and 70 ± 10% relative humidity, under a 16 h:8 h (Light: Dark) photoperiod ⁶⁹.

All Drosophila stocks were raised in standard cornmeal–molasses–agar medium maintained at 25° C, 60% humidity, and under 12 hr:12 hr light:dark conditions. We used Canton s flies as the wild-type strain. The relevant information of fruit flies used in the experiments are shown in key resources table.

Behavior assays of Nilaparvata lugens

Mating and receptivity

All BPHs were raised in incubators with a 16 h:8 h dark:light cycle. Unmated female and male BPHs were collected within 12 hours of emergence. For RNAi or neuropeptide injection treatment, the protocol is shown in Figure 1A and Figure 1D. The mating behavior (or female receptivity behavior) test in this study usually conducted the first mating between 3 and 5 days after emergence. The mating experiment was completed in a plastic tube (2.5cm in diameter and 10cm in height). Fresh rice seedlings of appropriate length were placed in each tube, and then a female and male BPHs were added. The mating conditions and other reproductive behavior parameters during this period were observed in a constant temperature room at 25℃ for 30min, and the males were taken out after 30 min. For the experiment requiring a second mating, females successfully mated for the first time were retained. After 48 hours, wild-type males were put in for the re-mating test. The males in all secondary mating experiments were wild type. The monitoring time for each mating behavior was 30 minutes; the receptivity rate refers to the proportion of successful matings of female individuals within the 30-minute measurement period.

Egg laying

After completing the first mating, the male is sucked out and a sufficient number of rice seedlings are added to the tube, and the female is retained to lay eggs in the tube. After 72 hours, the spawning female is transferred out of the tube. After the eggs have matured, about 4-5 days later, we used tweezers to gently manipulate the rice seedlings and counted the number of fertilized eggs produced by the females. Note that 4 to 5 days after the egg is produced, the successfully fertilized egg will display red eye-spots. (Oviposition assay protocol: For quantitative assessment of post-mating responses, egg-laying activity was monitored over a 72-hour observation window following verified initial copulation. To eliminate confounding effects from sequential mating trials, experimental subjects undergoing oviposition measurement were permanently excluded from subsequent re-mating evaluations.)

To determine the oviposition profile of females injected with CRZ peptide. virgin female or mated female BPHs that had emerged 3-5 days earlier were transferred into plastic tubes with sufficient amount of rice seedlings 6h after recovery from CRZ injection, and the BPHs were taken out 24 hours later, using tweezers to gently dissect the rice seedlings the number of fertilized eggs produced by the females were counted.

Courtship rate

During the mating assessment, an unmated (virgin) female and an unmated male BPH was placed in a plastic tube containing rice seedlings. The occurence of courtship behavior in the male individuals was recorded, which included a series of courtship actions such as following, wing display, licking, and attempting to mate. The proportion of individuals displaying courtship behavior among all experimental individuals is defined as the courtship rate.

Copulation duration

During the mating assessment, an unmated (virgin) female and an unmated male BPH were placed in a plastic tube containing rice seedlings. Timing begins when mating behavior is observed between the female and male, and stops immediately after completed copulation. The duration from the start of the mating behavior to its conclusion is defined as the mating duration.

Behavior assays of Drosophila melanogaster

Mating and re-mating : During the receptivity assessment, all experimental fruit flies were virgin flies that emerged within 12 hours in a 20°C incubator in a 12 h:12 h dark:light cycle. After collection, they were kept at 25°C incubator in a 12 h:12 h dark:light cycle for 4 days before behavioral testing. Receptivity behavior measurement : a perforated plate was placed on top of a non-perforated plate. Male fruit flies were briefly anesthetized on ice and gently transferred into the holes of the perforated plate, with one male per hole. Once the male transfer was completed, a plastic film was gently placed over the perforated plate, followed by covering it with another perforated plate. Then, the virgin female fruit fly, also anesthetized on ice, was placed in each hole. After the transfer was completed, a non-perforated plate was placed on top. This setup was then transferred to a 25°C incubator to acclimate for 30 minutes. The plastic film separating the male and female flies was quickly removed, and video recording began for 30 minutes. The recorded videos were analyzed to calculate the acceptance rate of the females. Successful mating females were collected and individually transferred into plastic tubes containing sufficient artificial food for fruit flies to lay eggs. After 24 hours of laying eggs, the female flies were transferred to a new finger-shaped tube with ample artificial feed for another 24 hours of laying. After 48 hours, these egg-laying females were subjected to a secondary mating assessment with wild-type, unmated males, using the same method as described above, to evaluate the acceptance rate during the second mating process. Additionally, the number of eggs laid by each female during the 48 hours was recorded.

Egg laying for mated females : First, collect virgin female fruit flies that eclosion within 12 hours at 20°C incubator in a 12 h:12 h dark:light cycle and then keep them in a 25°C incubator in a 12 h:12 h dark:light cycle for 4 days. After that, unmated males of the wild type and virgin females were mixed 1:1.5 times. 24 hours after full mating, female flies were lightly anesthetized with ice and transferred to plastic vials containing enough artificial fruit fly food for them to lay eggs for 24 hours. After 24 hours, the flies were transferred again to a new bottle full of artificial food for another 24 hours of egg laying. The number of eggs in each vial is calculated under a microscope, and the total number of eggs in the two vials of each female fruit fly represents the egg-laying capacity of that individual.

Developmental duration statistics

Nymphs of BPHs that hatched synchronously from distinct genotypes were collected and reared in plastic cups containing fresh rice seedlings; the seedlings were replaced regularly to maintain optimal conditions. From the day the nymphs reached the fifth instar, cups were inspected daily for adult emergence, and the number of newly emerged adults was recorded. After all surviving individuals had completed eclosion, the total nymph-to-adult developmental duration was calculated.

BPHs survival rate statistics

For neuropeptide injection: After 4-day emergence, virgin BPHs were injected with CRZ neuropeptide or scrambled CRZ, and after 6 hours of recovery, these insects were placed into a plastic tube with sufficient quantity of rice seedlings. Survival was monitored every 24 hours until all specimens were dead.

For mutant: Wild-type and mutant BPHs adults that eclosed on the same day were individually confined in plastic tubes, each supplied with ample fresh rice seedlings and water renewed regularly to ensure optimal living and nutritional conditions. The number of surviving individuals of each genotype was recorded every 24 h until all had died, and survival curves were generated from these data.

Neuropeptide synthesis and injection

Corazonin (CRZ) and scrambled CRZ (sCRZ) were synthesized by Genscript Co. Ltd (Nanjing, China). Peptide masses were confirmed by Mass spectrometry and the amount of peptide was quantified by amino acid analysis. The amino acid sequences of the peptides used in this study are: N. lugens corazonin: pQTFQYSRGWNamide; scrambled corazonin: pQSRYTTFGQWTNamide. On the fourth day after adult emergence, virgin female BPHs were injected with different concentrations of synthetic CRZ peptide and control scrambled peptide. The scrambled peptide control was used to monitor for potential non-specific effects on behavior due to the process of injection. Each mg peptide was dissolved in 100 μl DMSO solution (Acmec, #D54264), and then diluted to 1.5 mM, 0.15 mM, 0.015 mM and 0.0015 mM with 1×PBS as solvent. The BPHs were anesthetized with CO₂ for about 30 seconds and placed abdomen up on 2% agarose (Solarbia, #A8201) plate and 50 nl neuropeptide solution injected into each female, 30 nl neuropeptide solution injected into each male. The male virgins were injected on the third day of emergence. Six hours after neuropeptide injection the insects were used for behavior assay.

Protein concentration determination of male accessory glands

MAGs of 3-day-old unmated males (n = 48) were dissected in 1×PBS( Solarbio, #P1020 ) and then gently crushed in 60 μL of 1×PBS. The tissues were homogenized for 60 s and centrifuged at 10,000g for 15 min at 4° C and the supernatant removed. The protein concentration in the supernatant was determined using the BCA Protein Assay Kit (CWBIO, #CW0014S).

Gene cloning and sequence analysis

We used the NCBI database with BLAST programs to carry out sequence alignment and analysis. Then we predicted Open Reading Frames (ORFs) with EditSeq. The primers were designed by Primer designing tool in NCBI. Total RNA Extraction was using the TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instructions. The cDNA template used for cloning was synthesized using the Biotech M-MLV reverse transcription kit and the synthesized cDNA template was stored at -20°C. The transmembrane segments and topology of proteins were predicted by TMHMM v2.0 (http://www.cbs.dtu.dk/services/TMHMM-2.0/) ⁷⁰. Multiple alignments of the complete amino acid sequences were performed with Clustal Omega (http://www.ebi.ac.uk/Tools/msa/clustalo).

qRT-PCR for spatial expression pattern and RNAi efficiency

The nervous system (NS), fat body (Fb), digestive system (DS), epidermis (EP), male reproductive organs (MRO) and female reproductive organs (FRO) from 3-5 days old adults were dissected in 1×PBS, to determine the tissue-specific expression pattern. With the exception of the reproductive system, the rest of the tissue comes from both sexes equally. The adult BPHs were placed on ice under slight anesthesia and dissected under a stereoscopic microscope (ZEISS). At least 50 individuals were dissected as one sample for tissue-specific analysis with three replicates. All tissues were quickly transferred to Trizol (Invitrogen, Carlsbad, CA, USA) solution at 0° C after dissection.

For RNAi efficiency testing, we collected the whole bodies of BPHs and transferred to a 1.5 mL EP tube containing Trizol reagent. Total RNA extraction was using Trizol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instructions. The cDNA template used for cloning was synthesized using the Biotech M-MLV reverse transcription kit and the synthesized cDNA template was stored at -20°C.

Real-time qPCRs of the various samples used the UltraSYBR Mixture (with ROX) Kit (CWBIO, Beijing, China). The PCR was performed in 20μl reaction including 2 μl of 10-fold diluted cDNA, 1 μl of each primer (10 μM), 10μl 2 × UltraSYBR Mixture, and 6 μl RNase-free water. The PCR conditions used were as follows: initial incubation at 95°C for 10 min, followed by 40 cycles of 95°C for 10 s and 60°C for 45 s. N. lugens 18S rRNA and actin were used as an internal control (primers listed in Table1). Relative quantification was performed via the comparative 2^−ΔΔCT method ⁷¹.

RNAi and RNAi efficiency determination

For lab-synthesized dsRNA, gfp, NlCrzR and NlCrz were amplified by PCR using specific primers conjugated with the T7 RNA polymerase promoter (primers listed in Table1). dsRNA was synthesized by the MEGAscript T7 transcription kit (Ambion, Austin, TX, USA) according to the manufacturer’s instructions. Finally, the quality and size of the dsRNA products were verified by 1% agarose gel electrophoresis and the Nanodrop 1000 spectrophotometer and kept at -70°C until use.

3-5 days old virgin female and virgin male BPHs were used for injection of 50nl / 30nl of 5 μg/μl dsRNA per insect. Injection of an equal volume of dsgfp was used as negative control. After 48 hours, the insects were used behavior assays and gene relative expression analysis.

Neuropeptide injection after RNAi

48 hours after the dsCrzR injection, each insect was injected again with 10 ng of peptide dissolved in 50 nL of 1 × PBS balanced salt solution. Recovery was for six hours after injection. These insects were used for behavior assays. RNAi efficiency was examined by qPCR using a pool of ten individuals 48 hours after dsRNA injections.

In vitro synthesis of sgRNA and Cas9 mRNA

The single guide RNA (sgRNA) was designed as previously reported ⁷². Briefly, sgRNA was designed by manually searching genomic sequence around the region of the Crz and CrzR exon for the sequences corresponding to 5’-N_17-20NGG-3’, where NGG is the protospacer-adjacent motif (PAM) of SpCas9 and N is any nucleotide. For in vitro transcription of sgRNA, were synthesized in vitro using the GeneArt™Precision gRNA Synthesis Kit (Thermo Fisher Scientific, Vilnius, Lithuania) according to the instructions of manufacturer. The Cas9 protein (TrueCut™Cas9 Protein v2, Cat. NO. A36497) was purchased from Thermo Fisher Scientific (Shanghai, China).

Embryo microinjection, crossing and genotyping

The embryonic injection was performed as previously reported ⁷³. The female BPHs, which had been mated for 3-6 days after emergence, were selected to lay eggs on fresh rice leaves. After 30 minutes, the embryos were transferred onto double-sided tape attached to a microscope slide along the long edge by gently pressing the slide onto the dorsal surface of embryos. A few drops of halocarbon oil 700 was added to these eggs. About 0.5–1 nL mixture of Cas9 protein (200 ng/μL) and sgRNA (100 ng/μL) was injected into each egg using a FemtoJet and Inject Man NI 2 microinjection system (Eppendorf, Germany). After the injection, the eggs were placed in Petri dishes covered with moist filter paper, and the Petri dishes were placed in plant growth chamber at 27±1 °C, 70±10% RH for hatching. These eggs hatched into G0 generations.

Crossing and genotyping (Re-edit)

To establish homozygous mutants of Nilaparvata lugens through CRISPR-mediated gene editing, newly emerged G0 adults were outcrossed with wild-type counterparts to generate G1 progeny. Following successful mating and oviposition, genomic DNA was extracted from G0 adults for mutation screening. The target region surrounding NlCrzR was amplified by PCR using specific primers NlCrzR-check-F and NlCrzR-check-R (primer sequences listed in Table 1), followed by Sanger sequencing to identify G0 individuals carrying mutations.

Table 1.

Progeny (G1) from mutation-positive G0 parents were maintained through sibling crosses to establish G2 populations. Subsequent genotyping via Sanger sequencing was performed on G1 individuals to select breeding pairs where both male and female parents carried heterozygous mutations. The resulting G2 offspring from these validated pairs were then subjected to intercrossing to ultimately generate homozygous mutant lines.

Quantitative RT-PCR

The first-strand cDNA was synthesized with HiScript II Q RT SuperMix for qPCR (+gDNA wiper) kit (Vazyme, Nanjing, China) using an oligo (dT)18 primer and 500 ng total RNA template in a 10 μl reaction, following the instructions. Real-time qPCRs in the various samples used the UltraSYBR Mixture (with ROX) Kit (CWBIO, Beijing, China). The PCR was performed in 20 μl reaction including 4 μl of 10-fold diluted cDNA, 1 μl of each primer (10 μM), 10 μl 2 × UltraSYBR Mixture, and 6 μl RNase-free water. The PCR conditions used were as follows: initial incubation at 95°C for 10 min, followed by 40 cycles of 95°C for 10 s and 60°C for 45 s. N. lugens 18S rRNA or Drosophila rp49 were used as an internal control (Table 1). Relative quantification was performed via the comparative 2^−ΔΔCT method (Livak, K. J,2001).

Immunohistochemistry

Adult BPHs, 3-5 days old, were dissected under 1x phosphate-buffered saline (PBS; pH 7.4) in Schneider’s insect medium (S2). We dissected the brain, ventral nerve cord and reproductive system of the BPHs. These tissues were fixed in 4% paraformaldehyde in PBS for 30 min at room temperature. After extensive washing with PTX (0.5% Triton X100, 0.5% bovine serum albumin in PBS), blocked in 3% normal goat serum. Then, the tissues were incubated in Anti - CRZ for 12 hours at 4℃ and in secondary antibody for 12 h at 4℃. Primary antibodies used were: rabbit anti - CRZ (1:1000, Anti - CRZ was gift from Jan A. Veenstra), mouse anti-Bruchpilot (1:1000, Developmental Studies Hybridoma Bank nc82). Secondary antibodies used: donkey anti-rabbit IgG conjugated to Alexa 488 (1:500, R37118, Thermo Fisher Scientific ) and donkey anti-mouse IgG conjugated to Alexa Alexa 555 (1:500, R37115, Thermo Fisher Scientific). The samples were mounted in Vectorshield (Vector Laboratory). Images were acquired with Zeiss LSM 700 confocal microscopes, and were processed with Image J software. All antibodies were diluted in PTX solution. The primary antibody was diluted at a ratio of 1:1000, while the secondary antibody was diluted at a ratio of 1:500.

Adult female Drosophila melanogaster, 3-5 days old, were dissected under 1 x phosphate-buffered saline (PBS; pH 7.4) in Schneider’s insect medium (S2). We dissected the brain, ventral nerve cord and reproductive system. These tissues were fixed in 4% paraformaldehyde in PBS for 30 min at room temperature. After extensive washing with PTX (0.5% Triton X100, 0.5% bovine serum albumin in PBS), blocked in 3% normal goat serum. Then, the tissues were incubated in Anti - CRZ and Anti - nc82 for 12 hours at 4℃ and in secondary antibody for 12 h at 4℃. Primary antibodies used were: rabbit anti-CRZ (1:1000, A11122, Provided by Jan A. Veenstra). Secondary antibodies used: donkey anti-rabbit IgG conjugated to Alexa 488 and anti-mouse IgG conjugated to Alexa 555 (R37118, Thermo Fisher Scientific). The samples were mounted in Vectorshield (Vector Laboratory). Images were acquired with Zeiss LSM 700 confocal microscopes, and were processed with Image J software. All antibodies were diluted in PTX solution.

Fluorscent in situ hybridization (FISH)

FISH was performed as previously reported by Yan et al. ⁷⁴. Briefly, the reproductive system of the female BPH dissection was performed under 1x Phosphate Buffered Saline (PBS) supplemented with a protease and phosphatase inhibitor cocktail (Shyuanye, #R40012) on a clean, clear rubber pad using two fine forceps (Dumont, #0108-5-po). Dissected tissues were fixed using 4% paraformaldehyde (PFA) (Shyuanye, #22039) diluted in 1 x PBS with 0.1% Tween20 (0.1% PBST) at room temperature (RT) for 20 min on a shaker. Fixed tissues were washed three times with 0.1% PBST for 15 min each on a shaker. The tissues were dehydrated in 25%, 50%, 90% and 100% methanol, and then rehydrated in reverse order, performed at room temperature for 10 minutes at a time. The tissue was then washed twice with 0.1% PBST and shaken for 5 minutes. A post fixation was performed by incubating the tissue in 4% PFA in 1 x PBS on a shaker for 20 minutes, followed by two washes with PBST for 5 minutes each while shaking. The tissue was permeabilized with proteinase K (10 µg/ml) at room temperature for 2 minutes. After incubation, the tissue was washed with 0.1% PBST for 5 minutes while shaking. The PBST was removed, and pre-hybridization solution (Boster, #AR0152) was applied at 55°C for approximately 2 hours. The DIG probes were denatured in a metal bath at 80°C for 5 minutes and then immediately placed on ice. The denatured DIG-labeled probes were added to the hybridization solution, and hybridization was carried out overnight at 55°C. Then, the hybridized tissue was washed twice at 55°C for 40 minutes with warm pre-hybridization solution, followed by four washes at room temperature with 0.1% PBST for 10 minutes each. The detection of the DIG-labeled probes was performed using an anti-DIG conjugated fluorescent antibody (1:100 dilution, Jackson, #200-542-156), incubated at room temperature in 0.1% PBST for 2 hours, followed by three washes with 0.1% PBST. Finally, the tissue was mounted in mounting medium and imaged using a Zeiss confocal microscope.

Imaging of reproductive organs

Male BPHs of different genotypes and different treatments, and female BPHs mated with males of different genotypes, were used for dissecting reproductive organs on the fourth day after emergence. Motic Images Devices and Motic Images Plus 3.0 were used to capture images of these reproductive organs.

Quantification and statistical analysis

All graphs were generated using Prism 9 software (GraphPad Software, La Jolla, CA). Data presented in this study were first verified for normal distribution by D’Agostino– Pearson normality test. If normally distributed, Student’s t test was used for pairwise comparisons, and one-way ANOVA or two-way ANOVA was used for comparisons among multiple groups. If not normally distributed, Mann–Whitney test was used for pairwise comparisons, and Kruskal–Wallis test was used for comparisons among multiple groups, followed by Dunn’s multiple comparisons. All data are presented as mean ± s.e.m.

Data availability

All of the RNA sequence data in this article have been deposited in the China National Center for Bioinformation database and are accessible in PRJCA045433. All data generated or analyzed during this study are included in the manuscript and supporting files; source data files have been provided for all figures.

Acknowledgements

We thank Dr. Jan A. Veenstra for providing CRZ antibody. This project was supported by National Key R&D Program of China (2022YFD1700200) to S.-F.W., the National Natural Science Foundation of China (No. 32022011 & 32472542) to S.-F.W..

Additional files

Figure 1 - Supplementary Video 1

Additional information

Funding

National Key R&D Program of China (2022YFD1700200)

National Natural Science Foundation of China (32022011)

National Natural Science Foundation of China (32472542)