We employed network science methods to systematically analyze the structural features of male and hermaphrodite neural networks (Cook et al., 2019). In order to simplify network size to emphasize the functional connectivity, we merged neurons with similar functions into a singular node. For instance, neurons ADAL (ADA Left) and ADAR (ADA Right) were consolidated as node ADA, and their total synaptic connections were integrated as the connections of ADA. Neurons with diverse functions were maintained as separate nodes, such as ASEL and ASER. We generated two simplified functional neural networks: the male neural network that comprises 161 neuron nodes and the hermaphrodite neural network that is composed of 122 nodes (Figure 1—figure supplement 1). In Figure 1—figure supplement 1, synaptic connections between neurons are denoted by clockwise lines. The line width corresponds to the synapse number, which is the sum of chemical synapses and electrical synapses (Cook et al., 2019). Given that chemical synapses are one-directional while electrical synapses are two-directional, we considered that one electrical synapse between two neurons is equivalent to two reciprocal chemical synapses to generate the connection matrix (Supplementary file 1). We found that 120 nodes are shared in these two neural networks (Figure 1—figure supplement 2), suggesting the dimorphic behaviors are largely mediated by male-specific neurons or synaptic connection differences among sex-shared neurons.

Given the directional nature and variable synapse numbers of the synaptic connections, we abstracted these neural networks as two directed-weighted graphs (Figure 1—figure supplement 1). We used graph theory methods designed for directed-weighted graphs to facilitate feature detection in both the male and hermaphrodite neural networks. One obvious observation is that the absence of connections between pharyngeal neurons and other neurons in the male neural network suggests a relatively independent functionality of pharyngeal neurons in males, though several studies imply the non-synaptic communications in C. elegans (Dag et al., 2023; Shen et al., 2016).

The node strength measures how strongly a node directly possesses with other nodes in the network. We calculated the node strength for each node and revealed that the distribution of node strengths in both hermaphrodite and male networks follows an exponential distribution (Figure 1A). Overall, neurons in the male network exhibit a larger node strength (Figure 1B). Comparing the node strengths of sex-specific neurons and sex-shared neurons, we found that the increased node strengths of sex-specific neurons in the male neural network, but not in the hermaphrodite network (Figure 1C and D), indicated the stronger connections of sex-specific neurons in the male neural network. A substantial node strength indicates a neuron’s potential to exert dominance in shaping behavioral outputs. By classifying behavioral outputs of the top 20 neurons based on node strength, we predict enhanced mating and suppressed locomotion and chemotaxis behaviors in males (Figure 1E). Several sex-shared neurons are in the top 20 list in both neural networks, including AVA, ALA, and RIA, highlighting their important roles in both sexes (Supplementary file 2).

We further analyzed the node strength by dividing it into node out-strength and in-strength, representing the synaptic outputs or inputs of individual neurons, respectively. Both out-strength and in-strength exhibit an exponential distribution (Figure 1F and J). It is observed that neurons in males have stronger out-strength compared to those in hermaphrodites, particularly for male-specific neurons (Figure 1G–I). Interestingly, although the overall in-strength in both male and hermaphrodite neural networks is similar (Figure 1K), the analysis shows a difference in how sex-specific and sex-shared neurons receive inputs. Male-specific neurons receive more substantial input signals compared to sex-shared neurons in males (Figure 1L), suggesting a heightened sensitivity or importance in processing inputs related to sex-specific behaviors in males. In contrast, female-specific neurons in hermaphrodites exhibit the pattern of decreased in-strength (Figure 1M), highlighting a sex-dependent variation in neural network architecture and function.

The node strength takes into account the total synapse number of individual neurons, regardless of the number of connections with other neurons. Next, we calculated edge weight that is defined by assigning a number to each edge, representing the synapse number of the individual connection between two neurons. Our analysis indicated that, on average, male neural networks tend to have larger edge weights compared to hermaphrodite networks (Figure 1N). This suggests that the male network may have stronger synaptic connections per connection on average. A closer look at the edge weights of sex-specific and sex-shared neuron-related edges reveals a nuanced picture. In males, the edge weights associated with sex-specific neurons are comparable to those associated with sex-shared neurons (Figure 1O). This indicates a uniformity in how connections are established within the male neural network, regardless of whether the neurons are involved in sex-specific or shared functions. However, in hermaphrodites, the edge weights of sex-specific neurons are significantly smaller than those of sex-shared neurons (Figure 1P), suggesting that synaptic connections in the hermaphrodite network may be organized differently, with decreased modification of sexual-related behaviors.

Focusing on the top edge weights in the male neural network, we found that the consistent high weights of ALA-related edges in both sexes underscore the importance of ALA-mediated harsh touch escape behavior for individuals of both males and hermaphrodites (Figure 1N, Figure 1—figure supplement 3). However, in males, the weights of edges associated with AVA and RIA are notably higher (Figure 1N, Figure 1—figure supplement 3), indicating the increased AVA- or RIA-mediated behaviors in males.

In graph theory, the shortest path problem involves finding a path between two nodes in a graph that minimizes the sum of the weights. For a neural network represented as a graph, the shortest path between two neurons can provide insight into the functional correlation between them, especially when they are not directly connected. A shorter shortest path in the male neural network suggests a more cohesive interconnection of functional relationships among neurons in males (Figure 1Q). When comparing the shortest paths between two sex-shared neurons or those involving at least one sex-specific neuron, it is observed that the latter shows a tighter connection in males compared to hermaphrodites (Figure 1R and S).

Betweenness centrality in a graph measures the number of shortest paths that pass through a particular node, thus highlighting the key functional nodes within the network. The top betweenness centrality of AVA in both male and hermaphrodite networks (Figure 1T–V) underscores the critical role of AVA in both sexes. By examining the top 20 neurons ranked by betweenness centrality, we infer that males exhibit enhanced mechanosensation, mate-searching, and mating behaviors (Figure 1W). These predictions are based on the understanding that nodes with high betweenness centrality are often involved in a significant number of shortest paths, indicating that they are integral to the network’s function and connectivity. AVA, in particular, may be central to various behaviors and processes, including mechanosensory responses and social interactions such as mate-searching and mating.

Acknowledging the crucial influence of sex-specific neurons and their interactions in driving behaviors unique to each sex (Fenk and de Bono, 2015; Hart and Hobert, 2015; Huang et al., 2023; Koo et al., 2011; Susoy et al., 2021), we generated the sex-specific sub-neural networks that include the sex-specific neurons along with their directed-connected sex-shared neurons. This approach allows us to delve into the intrinsic properties of neural networks that exhibit sexual dimorphism (Figure 2—figure supplement 1).

Our analysis of sex-specific sub-neural networks has revealed that they share similar characteristics with their respective entire neural networks. Node strengths, which include total node strength, node out-strength, and node in-strength, are found to follow exponential distributions within these sex-specific sub-networks for both males and hermaphrodites (Figure 2A, E and I). In males, the sex-specific sub-neural network exhibits enhanced node strengths when compared to the hermaphrodite’s, indicating a more robust connectivity within the male network (Figure 2B). And the male sex-specific neurons reveal larger total node strengths and node out-strengths (Figure 2C and G), while the node in-strengths remain similar to those of hermaphrodites (Figure 2K). In contrast, hermaphrodites show no significant difference in total strengths, out-strengths, and in-strengths between sex-specific neurons and sex-shared neurons (Figure 2D, H and L). The large edge weights and short length of the shortest path within the male sex-specific sub-neural network suggest a tighter functional connection, indicative of a more integrated network (Figure 2M–R). Specifically, the top edge weights of AIY→RIA and RIA→RMD_DV in the male sub-neural network, along with the role of RIA in sensorimotor integration for olfactory steering, suggest that this AIY→RIA→RMD_DV circuit plays a significant role in mediating chemotaxis-related behaviors in males (Figure 2M). Additionally, the high betweenness centrality of AVA and the top edge weights of AVA→DA and AVA→VA within the male sub-neural network point to AVA’s regulatory role in locomotion, which is essential for male’s sex-related behaviors (Figure 2M and S). These findings underscore the importance of understanding the connectivity and functional roles of neurons within sex-specific neural networks and how these factors contribute to the manifestation of sex-specific behaviors.

Previous studies have demonstrated that sex-shared neurons can form dimorphic circuits that are responsible for mediating behaviors that vary between males and hermaphrodites (Serrano-Saiz et al., 2017; Setty et al., 2022; White et al., 2007). Therefore, we generated the sex-shared sub-neural networks containing sex-shared neurons and the connections among them (Figure 3—figure supplement 1) and analyzed the characteristics of these sex-shared sub-neural networks for both sexes.

Interestingly, our findings indicate that sex-shared sub-neural networks display distinct features when compared to the complete neural networks. Both of these sub-networks exhibit node strengths that follow a similar exponential distribution (Figure 3A, B, D, E, G and H). To quantify the differences in node strength between two sexes, we calculated the node strength difference for each neuron, defined as the node strength of the male neuron minus that of the hermaphrodite neuron (Figure 3C, F and I). Significantly, some neurons show an increase in node strength in males, including the commander neuron AVA, motor neurons VA and VB, and the interneuron RIA (Figure 3C, F and I). This increase in node strength suggests that the neural functions associated with AVA and RIA may be heightened in males. Conversely, the comparatively reduced node strength observed in several neurons in males could be indicative of behaviors that are specific to hermaphrodites, such as those related to neurons URY, PVC, and PVR (Figure 3C, F and I).

In examining the edge weight feature of the neural networks, we identified several edges that are specific to one sex (Supplementary file 3). These findings underscore the potential influence of these sex-specific edges on the generation of dimorphic behavioral outputs. For clarity, we focus on the top five weighted sex-specific edges and their edge counts, as indicated in Figure 3J. When comparing the weights of shared edges, it is observed that they tend to be more substantial in males (Figure 3K). This disparity in edge weight is particularly notable in circuits related to RIA, such as AIY→RIA→RMD_DV, and circuits associated with AVA, including AVA→DA and AVA→VA, which are markedly enriched in males. On the other hand, circuits like PVP-AQR and those involving FLP neurons, such as FLP→AVD and FLP→AVB, appear to be more dominant in hermaphrodites (Figure 3L). These circuit-specific differences suggest that they may contribute to distinct behavioral manifestations between two sexes.

Furthermore, we observed a decrease in the shortest paths among pharyngeal neurons, such as M3, M5, and NSM, in the male neural network (Figure 3M and N). Along with a lack of connections with other neurons (Figure 1—figure supplement 1), this observation indicates the formation of a functional unit consisting of pharyngeal neurons, which may be critical for male-specific behaviors. We have noted an increase in the shortest paths involving the RIP neuron, such as RIP→ASI, RIP→ADE, RIP→ALN, RIP→ASJ, and RIP→ASER (Figure 3N), along with a decrease in the betweenness centrality of the RIP neuron (Figure 3O and P). These changes strongly suggest an enhanced dimorphic function of the RIP neuron in hermaphrodites, which could be integral to hermaphrodite-specific behaviors.

To systemically study the dynamics of these two neural networks, we performed computational neuroscience methods to assess the functional role of individual neurons in the network. We utilized a non-spiking RC model to simulate neuronal activity because it is widely accepted that neurons in C. elegans are non-spiked, with the acknowledgment of action potentials in several neurons (Jiang et al., 2022; Liu et al., 2018b). To simplify the simulation, we standardized parameters for all neurons, maintaining a resting potential of –60 mV and the fixed linear correlation of the current input and synapse number (Source code 1). We stimulated individual neurons by shifting the membrane potential to +60 mV for a duration of 40 seconds and monitored the membrane potential changes of all the other neurons. Figure 4A illustrates the stimulation of AVA in hermaphrodite neural network as an example; the other neurons of hermaphrodite neural network exhibit diverse response patterns (Figure 4A). From these simulations, we identified the neuron that exhibited the strongest response when an individual neuron was stimulated, creating functional pairs. We found that only 30% of these functional pairs are conserved across both neural networks (Figure 4B and Supplementary file 4). To gauge the influence of a neuron as an upstream component within the network, we measured the cumulative changes in membrane potential across all other neurons when a particular neuron was stimulated. A substantial alteration in membrane potential across all other neurons would indicate that a particular neuron has a significant impact on the entire network as an upstream component. According to our results, neurons in the male network are more capable of activating the network than those in the hermaphrodite network (Figure 4C). Furthermore, sex-specific neurons, when acting as upstream neurons, have a greater effect on the network in males than do sex-shared neurons (Figure 4D), a trend not observed in hermaphrodites (Figure 4E). Similar patterns were observed when assessing the dynamic functions of neurons as part of the downstream component of the network (Figure 4G–I). By categorizing the top 20 upstream and top 20 downstream neurons, we attempted to predict sexually dimorphic behaviors. Both predictions highlighted mating behavior as a key aspect of male-specific behaviors (Figure 4F and J).

Given the structural and dynamical discrepancies identified between the male and hermaphrodite neural networks, several sexually dimorphic behaviors can be predicted. In particular, male behavior is primarily characterized by an emphasis on sexual-related behaviors. The AVA neuron emerges as pivotal to the function of the male network, participating in both sex-specific and sex-shared sub-neural networks. Considering that AVA regulates local search behavior in hermaphrodite (Gray et al., 2005; López-Cruz et al., 2019) and enhanced connections of AVA in the male neural network, therefore, we examined the hypothesis that AVA-related neural circuit could be a key regulatory component for dimorphic spontaneous local search behavior in males.

We placed individual worms onto distinct food patches and recorded the number of reversals over the initial ten minutes in the new environment; this number indicates local search behavior (Gray et al., 2005). Results showed that males display high levels of local search behavior on both thick and thin food patches, as well as patches without food (Figure 5A-D). The absence of food alters the initial local search behavior in hermaphrodites, with no impact on males (Figure 5C and D). The males’ enhanced local search behavior is dependent on sexual maturity, as evidenced by the similar reversal rates of late larval males and hermaphrodites (Figure 5E). To explore the potential impact of the prior experiences during development on male’s local search behavior, we separated L4 males on three distinct types of plates: a plate devoid of both food and hermaphrodites, a plate containing food without hermaphrodites or a plate with both food and hermaphrodites. There was no significant difference among males with different prior experiences, suggesting that developmental experiences do not influence the enhanced local search behavior in males (Figure 5F). Overall, we discovered a sexually dimorphic behavior in C. elegans, that sexually matured males exhibit enhanced local search behavior regardless of the environment.

To examine AVA’s role in local search behavior, we temporally suppressed AVA activity by expressing HisCl1 (histamine-gated chloride channel) in AVA and administrating histamine (Pokala et al., 2014). The significantly reduced reversal rates indicate the commanding role of AVA in driving reversals in both males and hermaphrodites (Figure 5G and H). Subsequently, we conducted in vivo calcium imaging studies on AVA neurons. Our findings revealed that AVA neurons exhibit spontaneous calcium transients, with adult males showing a significantly higher frequency of these events compared to hermaphrodites (Figure 5I and J). However, this difference is not observed in late larval males (Figure 5—figure supplement 1), which aligns with the behavioral outcomes.

Next, to identify male’s local search target, we assessed male’s behavioral preference towards hermaphrodites. We generated four food patches with equal distance to the center of the plate, then placed approximately 30 adult hermaphrodites on diagonal food patches for 1 hour, followed by introducing approximately 30 males to the center of the plate (Figure 5K). After 1 hour, we quantified the distribution of males across food patches with and without hermaphrodites. A greater proportion of wild-type males congregated on patches containing hermaphrodites. However, inhibiting AVA activity notably reduced the preference for hermaphrodites (Figure 5L). These findings suggest that AVA-mediated enhancement of male local search behavior aids in locating hermaphrodites for mating purposes.

In order to elucidate the AVA-related dimorphic circuit responsible for male-specific local search behavior, we quantified the differences in chemical synaptic connections with AVA neurons between males and hermaphrodites. Notably, PQR, RIC, and DVC neurons demonstrated the highest synaptic output to AVA (Figure 6A and B, Figure 6—figure supplement 1), suggesting their potential role as upstream neurons. Consequently, we investigated the effect of PQR, RIC, and DVC on male’s enhanced local search behavior by blocking their synaptic release using TeTx expression driven by cell-specific promoters (Figure 6—figure supplement 2). Our findings indicated that selectively blocking either RIC or DVC significantly reduced the local search behavior in males (Figure 6C and E), but not in hermaphrodites (Figure 6D and F). However, PQR is not involved in local search regulation in both male and hermaphrodite (Figure 6G, H and Figure 6—figure supplement 3). To further validate these results, we conducted AVA calcium imaging under conditions where PQR, DVC, or RIC synaptic activity was blocked. A decrease in the frequency of spontaneous calcium transients was observed specifically when DVC or RIC activity was suppressed, indicating that both DVC and RIC are indeed upstream neurons of AVA and play a crucial role in modulating the enhanced local search behavior observed in males (Figure 6I-N). These results are consistent with the behavioral observations.

DVC identified as a glutamatergic neuron (Pereira et al., 2015) and activates AVA calcium transient events for local search behavior in males (Figure 6E and K). Therefore, we tested whether NMDA or AMPA receptors on AVA regulate local search behavior. We generated transgenic strains that individually silence NMDA and AMPA receptor genes in AVA neurons (Figure 7A). By examining the reversal rates, we evaluated the impact on local search behavior. Our results indicated that individually silencing AMPA receptor gene glr-2, as well as NMDA receptor genes, nmr-1 or nmr-2, led to a significant decrease in reversal rates in males (Figure 7B, D and F). However, no such effect was observed in hermaphrodites (Figure 7C, E and G). This suggests that both AMPA and NMDA receptors are crucial for mediating local search behavior in males. Furthermore, since RIC is known to release octopamine (Pereira et al., 2015), we explored the influence of octopamine receptors SER-3 and SER-6 in AVA on local search behavior (Figure 7A). By silencing the ser-3 and ser-6 genes specifically in AVA neurons and assessing the reversal rates, we aimed to determine their contribution. Our findings demonstrated that silencing ser-3 in AVA significantly reduced the reversal rate in males, without affecting hermaphrodites (Figure 7H and I). In contrast, silencing ser-6 in AVA had no discernible effect on the reversal rate in either sex (Figure 7J and K, Figure 7—figure supplement 1). These results suggest that the octopamine receptor SER-3 plays a specific role in regulating the increased local search behavior observed in males.

The overwhelming majority (99.9%) of the C. elegans population consists of hermaphrodites, which prefer self-fertilization over mating (Ebert and Bargmann, 2024; Garcia et al., 2007; Woodruff et al., 2014; Wu et al., 2023). Therefore, the evolution pressure has driven the rare males to develop multiple mechanisms to increase the mating efficiency.

We find that the male’s neural network contains about 25% (41/161) male-specific neurons, while the neural network in hermaphrodite only contains about 1.6% (2/122) hermaphrodite-specific neurons that regulates hermaphrodite-specific egg-laying behavior (Emmons, 2018). Structural analysis shows that male-specific neurons form strong connections within male neural networks, revealing the important roles of these neurons and setting the basis for males’ prioritized sexual-related behavior. Our functional simulation analysis confirms the prediction based on structural features, emphasizing the substantial impact of male-specific neurons on the male neural network.

Not only the dominant role of male-specific neurons, the distinct connections and functions of sex-shared neurons facilitate sexual-related behaviors through different mechanisms. Males express pheromone receptors in both sex-shared sensory neurons (Aprison and Ruvinsky, 2017; Fagan et al., 2018) and male-specific neurons (Narayan et al., 2016; Reilly et al., 2023) to increase pheromone sensitivity. Besides the increased pheromone receptors in the sensory neurons, males rewire the synaptic connections at the interneuron and motor neuron level to promote sex-related behavior. AVA neuron is the command neuron for driving reversals, which plays an essential role for hermaphrodite behaviors, including navigation, food searching, aversive response, learning, and memory, etc. (Choi et al., 2020; Gray et al., 2005; Liu et al., 2020). We find that AVA has strong connections with sex-specific neurons and shows the key position for the functional connection in the whole male network as well as the male-specific sub-network. It suggests that AVA plays an important role for mating behavior in males. Indeed, a previous study demonstrated that AVA integrates the signals for sex-specific neurons to complete a series of stereotyped steps during mating (Sherlekar et al., 2013). Beyond directly regulating mating behavior, we find that the AVA-related rewired circuit in males mediates increased local search behavior that indirectly enhances the mating efficiency. The highly increased connections between command neuron AVA and downstream motor neurons VA and DA indicate a strong motor control in males. Sex-shared neurons DVC and RIC have increased chemical synapses with AVA in male, and our behavioral assay and calcium imaging results show that the circuit involving glutamatergic neuron DVC, octopaminergic neuron RIC, and command neuron AVA regulate enhanced spontaneous local search behavior that facilitates mate-searching. For further prioritizing the value of hermaphrodites, males decrease food attraction to promote food lawn leaving (Ryan et al., 2014), and decrease aversive response to nociceptive stimuli (Pechuk et al., 2022).

In C. elegans, pharyngeal neurons that govern the pharyngeal pumping mechanism for feeding are organized into a specialized sub-network (Cook et al., 2020). Previous research has demonstrated that signals from non-pharyngeal neurons can modulate feeding behaviors through regulating pharyngeal neurons in hermaphrodites (Dallière et al., 2016; Trojanowski et al., 2016). Structural and functional studies characterized RIP neuron as a pivotal intermediary to connect pharyngeal neurons with others in hermaphrodite neural network (Albertson and Thomson, 1976; Bhatla et al., 2015). Interestingly, our current study reveals a distinct pattern in male worms. Unlike previous findings, we observed no direct connections between the pharyngeal neurons and other neurons, as depicted in Figure 1—figure supplement 1. Instead, we noted a robust interconnection among the pharynx neurons themselves and a loose functional linkage of the RIP neuron with sensory neurons responsible for internal and external environment changes (Davis et al., 2018; Guo et al., 2018; Otarigho et al., 2024; Sato et al., 2021) in neural network (Figure 1N). These findings suggest that males may possess an independent regulatory mechanism for feeding behaviors.

Our previous study demonstrated that the AIY→RIA circuit is fundamental for sensorimotor integration in olfactory steering in hermaphrodites (Liu et al., 2022b; Liu et al., 2018a). AIY receives inputs from various sensory neurons, including AWC and AWA, in both sexes (Figure 1—figure supplement 4). Notably, AWC and AWA in males are involved in sexual attraction (Wan et al., 2019; White et al., 2007). Here, we discover an enhanced circuit at the inter-motor neuron level, AIY→RIA→RMD_DV circuit (Figure 1K and L); therefore, we predict that enhanced AIY→RIA→RMD_DV circuit facilitates pheromone-guided olfactory steering behavior in males. A recent study showed that AFD expresses a dimorphic level of the insulin-like peptide INS-39 in males, suggesting a dimorphic function of AFD (Haque et al., 2024). Our analysis shows that neuron AFD exhibits the highest number of synaptic connections with AIY (Figure 1—figure supplement 4); therefore, further exploring the function of the AFD→AIY circuit in males is of significant interest.

Previous functional studies have demonstrated that AVA, as the command neuron, drives locomotion to regulate diverse behaviors in both sexes (Bayer and Hobert, 2018; Choi et al., 2020; Gat et al., 2023; Gray et al., 2005; Lin et al., 2024; Pechuk et al., 2022). Here, we present structural evidence for the essential role of AVA in the neural networks of males and hermaphrodites. We find that AVA has a large number of synapses, including both synaptic inputs and outputs, in two neural networks (Figure 1B, G and K). Interestingly, with respect to the difference in synapse numbers among sex-shared neurons, AVA exhibits a dimorphic connection pattern in males and hermaphrodites, as well as playing a key role in the sex-specific sub-neural network of males. Therefore, it is reasonable to predict that a wide range of behaviors mediated by AVA are sexually dimorphic. Here, we have demonstrated that AVA and the upstream neurons DVC and RIC regulate enhanced local search behavior in males (Figures 5—7). However, PQR, which reveals the largest synapse number difference with AVA between the two sexes, is not involved. Thus, it is intriguing to investigate the function of the PQR-AVA circuit in males to discover how this dimorphic circuit mediates sexually specific behavior.

We have identified several synaptic connections that exist only in one sex (Supplementary file 3 and Figure 3J). One male-specific circuit involving PHA and two hermaphrodite-specific circuits comprising PHB and URY are highlighted according to the sex-specific connections (Figure 2—figure supplement 2). Indeed, studies have demonstrated the mechanisms of sex-specific pruning of PHA and PHB neuronal synapses and the corresponding sexually dimorphic behaviors (Barrios et al., 2012; Oren-Suissa et al., 2016; Setty et al., 2022). Sensory neuron URY is involved in pathogen avoidance in hermaphrodites (Gallrein et al., 2021; Pujol et al., 2001). Therefore, it is intriguing to examine whether URY activates downstream motor neurons RMD_DV and SMD or interneuron AVE, to generate avoidance behavior in a sex-specific manner.

Overall, this study elucidates the neural circuitry underlying prioritized sexual behaviors in males, also predicting potential dimorphic behaviors and the associated circuits.

C. elegans strains maintenance followed the protocol from wormbook (Stiernagle, 2006). The Bristol N2 strain was used as wild-type. Males were generated through heat shock and maintained in a mixed population with hermaphrodites (Anderson et al., 2010). The following strains were employed in this study, with their detailed genotypes or microinjection mixtures described within square brackets:

The plasmids used to generate the aforementioned transgenic strains were constructed via Gibson assembly (Gibson et al., 2009), utilizing the Gibson Assembly Master Mix (NEB, E2611) according to the manufacturer’s protocol. Gibson assembly is an isothermal, single-reaction method that assembles multiple overlapping DNA fragments through the concerted action of a 5’ exonuclease, DNA polymerase, and DNA ligase. The DNA fragments (listed in Supplementary file 5) used for Gibson assembly were generated via PCR with the 2×High-Fidelity PCR Master Mix (MedChemExpress, Cat. No.: HY-K0533). Primer and template details are provided in Supplementary file 5. The identity of each plasmid was confirmed by sequencing.

For cell-specific expression, the following promoters were employed: AVA neuron expression: Driven by an overlapping strategy using the flp-18 promoter (Pflp-18, 1527 bp upstream of the flp-18 gene) and the nmr-1s promoter (Pnmr-1s, 990 bp upstream of the nmr-1s gene) Zhang and Zhang, 2012; RIC neuron-specific expression: Driven by the tbh-1 promoter (Ptbh-1, 2819 bp upstream of the tbh-1 gene) Suo et al., 2006; DVC neuron-specific expression: Driven by the ceh-63 promoter (Pceh-63, 646 bp upstream of the ceh-63 gene) Oh et al., 2019; PQR neuron expression: Driven by the gcy-36 promoter (Pgcy-26, 1089 bp upstream of the gcy-36 gene) (Yu et al., 2017).

Cell-specific gene knockdown was carried out based on the methodology outlined in a previous publication (Esposito et al., 2007) with some modifications. In this approach, double-stranded RNAs (dsRNAs), consisting of both sense and antisense RNA strands, are transcribed using cell-specific promoters. The DNA fragments responsible for generating the promoter::sense or promoter::antisense RNA were subcloned into a plasmid, which was sequenced to confirm the accuracy of the dsRNA constructs.

For AVA neuron-specific RNA interference (RNAi), the expression of sense and antisense RNA strands was achieved using the aforementioned overlapping promoter strategy, where the flp-18 promoter drove the expression of the sense RNA, and the nmr-1s promoter drove the expression of the antisense RNA (Zhang and Zhang, 2012). This promoter-overlapping approach facilitated the generation of dsRNA specifically in AVA neurons. Details of the DNA fragments used to produce the dsRNAs are provided in Supplementary file 5.

Plasmids designed to express sense and antisense RNA were constructed using Gibson assembly. These plasmids, along with a microinjection marker plasmid (CC::GFP or CC::DsRed2), were injected into the germline of C. elegans to create RNAi transgenic lines, specifically PSC132, PSC136, PSC149, PSC233, PSC320, and PSC321.

In the graph theory analysis, several parameters were calculated. The definitions for each parameter are as follows:

We generated an artificial network following the same principles as the neural networks in C. elegans (Figure 1—figure supplement 5) and used it as an example to demonstrate the calculation of each parameter.

Taking neuron D as an example, there are six shortest paths that pass through neuron D. Therefore, the betweenness centrality of neuron D is 6/12=0.5.

We model the dynamics of neural networks in both sexes of Caenorhabditis elegans using the RC neuron model. The RC neuron model is represented as follows:cdVdt=−gNa(V−ENa)−gK(V−EK)−gL(V−EL)+Iext\begin{document}$$\displaystyle c\frac{\mathrm{d}V}{\mathrm{d}t}=-g_{Na}(V-E_{Na})-g_{K}(V-E_{K})-g_{L}(V-E_{L })+I_{ext}$$\end{document}

In this equation, V\begin{document}$V$\end{document} is the membrane potential, c=1.0\begin{document}$c=1.0$\end{document} is the membrane capacitance, gNa=120.0\begin{document}$g_{Na}=120.0$\end{document}, gK=36.0\begin{document}$g_{K}=36.0$\end{document}, and gL=0.03\begin{document}$g_{L}=0.03$\end{document} are the conductance of sodium ion channel, potassium ion channel, and leakage channel, respectively. ENa=50.0mV\begin{document}$E_{Na}=50.0\,\rm mV$\end{document}, EK=−77.0mV\begin{document}$E_{K}=-77.0\, \rm mV$\end{document}, and EL=−54.4mV\begin{document}$E_{L}=-54.4\,\rm mV$\end{document} are the reversal potentials of the sodium ion, potassium ion, and leakage current, respectively. Iext\begin{document}$I_{ext}$\end{document} is the external input current. A neuron receives synaptic currents from all upstream neurons, and the sum of all synaptic currents is its external input current Iext\begin{document}$I_{ext}$\end{document}. If there are N\begin{document}$N$\end{document} synapses connecting neuron A to neuron B, then neuron B will receive a synaptic current IA→B\begin{document}$I_{A\rightarrow B}$\end{document},IA→B=kN(VA−VB)\begin{document}$$\displaystyle I_{A\rightarrow B}=kN(V_{A}-V_{B})$$\end{document}

where k=100.0\begin{document}$k=100.0$\end{document} is the coefficient, VA\begin{document}$V_{A}$\end{document} and VB\begin{document}$V_{B}$\end{document} are the membrane potentials of neuron A and neuron B.

The details of the computation are provided in the Python code included in the supplementary materials (Source code 1).

A 10 uL of OP50 bacterial liquid culture was transferred onto a 3 cm NGM plate and allowed to dry at room temperature for either 0.5 hours or 16 hours to differentiate food thickness. Individual one-day-old adult male or hermaphrodite C. elegans were then placed on the assay plate. The plate was positioned under a high-speed camera (MV-CA060-11GM, Hikvision) to record for 10 minutes at a frame rate of 10 Hz. Reversal behaviors were manually counted for further analysis. In the subsequent local search behavioral assay with transgenic animals, only the last 3 minutes of the 10-minute search period were recorded and analyzed.

The protocol for temporally inhibiting AVA neurons with histamine treatment in AVA::hisCl1 transgenic animals was adapted from a previous study (Pokala et al., 2014). A 1 M histamine-dihydrochloride (H7125, Sigma-Aldrich) stock solution was prepared and stored in a –4°C freezer. For experimental use, the stock solution was diluted with double-distilled water (ddH₂O) to a final concentration of 200 mM. To prepare the histamine-containing food lawn, 10 uL of the 200 mM histamine-dihydrochloride solution was mixed with 190 uL of E. coli OP50 culture to achieve a final concentration of 10 mM histamine. Assay plates were prepared by seeding 10 uL of this bacterial mixture onto the center of 3 cm NGM plates and allowing the plates to dry for 0.5–3 hours before the experiment.

A 9 cm NGM assay plate pretreated with 10 uM histamine solution (H7125, Sigma-Aldrich) was divided into four equal sectors. A food lawn was created at the center of each sector by drying 50 uL of OP50 bacterial liquid culture for 2 hours at room temperature. Approximately 30 one-day-old adult GR1333 hermaphrodites were transferred onto two diagonal food lawns, initiating a 2-hour pheromone-releasing period. After this period, about 30 one-day-old adult males were quickly placed at the center of the assay plate and allowed to freely select a sector for 1 hour. Following the selection period, the number of males in each sector was manually counted. The Choice Index was calculated based on the distribution of males across the sectors:Choiceindex(CI)=n(lawnswithpheromone)−n(lawnwithoutpheromone)n(lawnswithpheromone)+n(lawnwithoutpheromone)\begin{document}$$\displaystyle Choice\, index(CI)=\frac{n(lawns\, with\, pheromone)-n(lawn\, without\, pheromone)} {n(lawns\, with\, pheromone)+n(lawn\, without\, pheromone)}$$\end{document}

One-day-old males or hermaphrodites were paralyzed with 10 mM levamisole (A606644-0100, Sango Biotech) and placed on a glass slide with a thin layer of NGM agar. Images were captured for ten minutes at a frame rate of 4.84 frames per second using a Zeiss LSM980 confocal microscope equipped with an EMCCD camera and a ×60 oil immersion objective. The images were processed using ImageJ software. The dynamic tracking plugin “TurboReg” was used to pre-process the raw image gallery and correct any frame misalignment. Regions of interest (ROIs) were manually drawn around the AVA soma, and the gray values of these ROIs were measured to obtain the time series of calcium signals. Approximately 144 frames (about 5% of the total 2876–2879 frames) were averaged to establish the F0\begin{document}$F_{0}$\end{document} baseline value. △F\begin{document}$\bigtriangleup F$\end{document} was calculated as △F=F−F0\begin{document}$\bigtriangleup F=F-F_{0}$\end{document}, where F represents the average fluorescence intensity of the ROI in each frame. The frequency of calcium events was determined by counting the number of signal intervals with △F/F0\begin{document}$\bigtriangleup F/F_{0}$\end{document} exceeding a threshold value of 1 during the 10-minute recording period.

No explicit statistical power analysis was used to pre-determine sample sizes. Instead, sample sizes were chosen based on standard practices widely accepted in the C. elegans field to ensure adequate statistical power. All experiments were performed independently at least two to three times (for technical replicates on distinct days). Unless otherwise specified, the sample size (n\begin{document}$n$\end{document}) represents the number of independent biological replicates (e.g., independent synchronized worm populations or individual animals) rather than technical replicates. Age-synchronized populations of C. elegans were randomly allocated to experimental and control groups. Due to the obvious phenotypic differences or the nature of the experimental setup, investigators were not blinded to group allocation during data collection or analysis. Pre-established criteria for data exclusion were not required, and no outliers, data points, or animals were omitted from the final statistical analyses. All statistical analyses were performed using GraphPad Prism version 9.5. Data distribution was assumed to be normal, and variances were similar between groups being compared. Differences between two distinct groups were evaluated using an unpaired two-tailed Student’s t-test (or Mann–Whitney U test for non-parametric data). For comparisons among three or more groups, a one-way analysis of variance (ANOVA) followed by Tukey’s post hoc test was applied. Summary statistics are presented as mean ± standard error of the mean (SEM) In accordance with eLife guidelines, exact p values, exact sample sizes (n\begin{document}$n$\end{document}), and specific statistical tests used for each experiment are explicitly detailed in the respective figures.