Anticipatory 50-kHz Precontact Ultrasonic Vocalizations and Sexual Motivation: Characteristic Pattern of Ultrasound Subtypes in an Individual Analyzed Profile

We verified the hypothesis of the existence of forms of individual-specific differences in the emission of anticipatory precontact vocalization (PVs) indicating individualization related to sexual experience and motivation in male rats. Long-Evans males were individually placed in a chamber and 50-kHz ultrasounds were recorded during 5-min periods. In experiment 1, PVs were recorded before the introduction of a female in four consecutive sessions during the acquisition of sexual experience. In experiment 2, PVs were analyzed in three groups of sexually experienced males: with the highest, moderate, and the lowest sexual motivation based on previous copulatory activity. In both experiments, the total number of ultrasounds, as well as 14 different specific subtypes, was measured. The ultrasound profiles for each male were created by analyzing the proportions of specific dominant subtypes of so-called 50-kHz calls. We decided that the dominant ultrasounds were those that represented more than 10% of the total recorded signals in a particular session. The number of PVs was positively correlated with the acquisition of sexual experience and previous copulatory efficiency (measured as the number of sessions with ejaculation). Furthermore, PVs showed domination of the frequency modulated signals (complex and composite) as well as flat and short with upward ramp ultrasounds with some individual differences, regardless of the level of sexual motivation. The results show a characteristic pattern of PVs and confirm the hypothesis that the number of PVs is a parameter reflecting the level of sexual motivation.

We verified the hypothesis of the existence of forms of individual-specific differences in the emission of anticipatory precontact vocalization (PVs) indicating individualization related to sexual experience and motivation in male rats. Long-Evans males were individually placed in a chamber and 50-kHz ultrasounds were recorded during 5-min periods. In experiment 1, PVs were recorded before the introduction of a female in four consecutive sessions during the acquisition of sexual experience. In experiment 2, PVs were analyzed in three groups of sexually experienced males: with the highest, moderate, and the lowest sexual motivation based on previous copulatory activity. In both experiments, the total number of ultrasounds, as well as 14 different specific subtypes, was measured. The ultrasound profiles for each male were created by analyzing the proportions of specific dominant subtypes of so-called 50-kHz calls. We decided that the dominant ultrasounds were those that represented more than 10% of the total recorded signals in a particular session. The number of PVs was positively correlated with the acquisition of sexual experience and previous copulatory efficiency (measured as the number of sessions with ejaculation). Furthermore, PVs showed domination of the frequency modulated signals (complex and composite) as well as flat and short with upward ramp ultrasounds with some individual differences, regardless of the level of sexual motivation. The results show a characteristic pattern of PVs and confirm the hypothesis that the number of PVs is a parameter reflecting the level of sexual motivation.
Contrary to 50-kHz USV and according to the arousal hypothesis, 22-kHz USV occurs during abrupt decreases of arousal (Bell, 1974) and in this case can reflect a relaxation state after ejaculation (Bialy et al., 2016) or a safety signal during aversive conditioning (Jelen et al., 2003).
Although USV is accompanied by certain behavioral situations, the question about the function of 50-kHz vocalizations is still open. One hypothesis assuming a non-semiotic communicative character of vocalization was proposed by Sales and Pye in 1974. An alternative hypothesis postulated that USV can simply be artifacts associated with breathing patterns in rodents (Blumberg and Alberts, 1991). Some researchers have also suggested that it is possible to accurately assess a rat's emotional state based on the profile of the signals emitted (Brudzynski, 2015;Simola and Granon, 2019).
In this context, the phenomenon of USV can be considered as a complex central response containing an easily measurable autonomic reflex component (objective component) along with a co-occurring mental constituent (subjective component) that is difficult to measure. While the precise division between subjective and objective emotional components is not easy to establish, an accurate behavioral analysis may assess the proportion between them. To categorize some aspects of behavior as a reflection of subjective mental processes, it is necessary to demonstrate individual differences in USV expression. This form of individualization in the case of USV emission was shown in mice (Holy and Guo, 2005;Arriaga and Jarvis, 2013;Asaba et al., 2014;Chabout et al., 2015;Zala et al., 2020).
Moreover, Matsumoto et al. (2016) revealed that rats have the potential ability to discriminate calls emitted by their own vs. conspecific. Based on electrophysiological measurements, authors showed an independent reaction of different groups of neurons in the dorsal amygdala (area of the lateral amygdaloid nucleus) responded specifically to own vs. conspecifics ultrasounds calls.
Furthermore, rats seem to be able to diminish the emission of USV while playing hide-and-seek with humans when hiding themselves which requires being quiet even when having an elevated level of arousal (Reinhold et al., 2019). Such results provide a possible approach for the conscious use of USV by rats.
In sexual behavior, 50-kHz vocalization can be observed during anticipation, during the initiation of copulation, during the copulatory performance, and during the late phase of the postejaculatory interval (Barfield et al., 1979;Bialy et al., 1996Bialy et al., , 2019b. Furthermore, the medial preoptic area and the nucleus accumbens-the neuronal circuits important in sexual activity (Hull et al., 2002) are involved in 50-kHz emissions (Fu and Brudzynski, 1994), also during sexual activity (Gao et al., 2019;Karigo et al., 2021), and anticipatory sexual behavior (Bialy et al., 2010). The anticipatory precontact 50-kHz vocalization (PVs) is convenient for the analysis of individual male vocalizations related to general arousal along with sexual motivation (Bialy et al., 2019a). The usefulness of this model derives, from the possibility, to separately analyze the ultrasounds emitted by the male in a conveniently short period (absence of the female provides a single source of USV). It allows describing the factors that trigger and modify individual vocalization in an appetitive state. Furthermore, it can also be a convenient parameter in the experimental models based on the acquisition of sociosexual experience and extinction reactions, facilitate recognition of crucial cues regulating acquisition/extinction reactions. It is due to the direct relationship of PVs with the memory of emotional state, copulatory efficiency, and reward value of preceding socio-sexual interactions (Bialy et al., 2000).
However, this relationship concerns the total number of emitted ultrasounds, while the question of the relationship between the level of sexual motivation and emitted subtypes remains unknown. In the present experiments, we recorded PVs emitted during the acquisition of sexual experience as well as those emitted by sexually experienced males with different levels of copulatory activity. We have classified 50-kHz calls into one of 14 subtypes (Wright et al., 2010) and created individual USV profiles for each male. We compared the profiles thus obtained with the copulatory history of each rat. We would answer whether PVs subtypes profiles have the traits suggesting the uniqueness for each animal. If they have, is there any relation between this emerging inter-individual differentiation and the level of sexual motivation for each animal?
The purpose of this study was to answer the question: are there individual-specific distinctions of 50-kHz vocalization associated with the level of sexual experience and sexual motivation in the ultrasound recordings? If not, what is the general physiological pattern of anticipatory precontact vocalizations in the sexual behavior of the male rat?

Animals
Long-Evans rats, 6-7 months old, were the subjects in this study. The choice of animals' age was dictated by the stabilized profile of sexual parameters occurring between 150 and 500 days of their life (allows to omit the first life period with significant fluctuation in levels of sex hormones and corresponding tissue sensitivity) (Larsson, 1956(Larsson, , 1967. Males and females were housed in separate rooms with a reversed 12-h light-dark cycle (lights switched off at 09:00 h) and at a temperature maintained at 22 ± 1 • C. All of the animals had food and water freely available. The food consisted of standard laboratory chow (experiment 1 and 2) with some enrichment in experiment 2 (cereals, fresh vegetables, and fruits) due to the effort to achieve more natural-like conditions conducive to the emergence of social hierarchies. The males arrived from the Department of Experimental Medicine, Medical University of Silesia, Katowice, from different cohorts for experiments 1 and 2. The males were sexually naïve at the beginning of the experiments and have been described separately in the corresponding experiments.
The housing conditions of males are described in corresponding "experiment 1 and 2" subsections.
All of the cages were provided with wood shavings and dedicated plastic tubes as enrichment. The ovariectomized females were brought into estrus with a single injection of estradiol benzoate (50 µg/rat s.c., Sigma-Aldrich) and progesterone (500 µg/rat s.c., Sigma-Aldrich). Hormones were dissolved in sesame oil and administered at a dose of 0.05 ml per individual. Hormonal injections were given 48-72 h before the test for estradiol and 4-8 h before the test for progesterone. Estrus was induced not more often than once a week and not less often than once every 2 weeks. The females were sexually experienced at the beginning of experiments. During an experimental day, a female copulated with up to two males.

Behavioral Procedures
All of the behavioral tests were conducted between 13:00 and 17:00 h during the dark phase of the light-dark cycle. We maintained at least a 1-week interval between tests to counteract the influence of the sexual exhaustion phenomenon on the sexual parameters, which is particularly important in less active groups of males (Larsson, 1956).
The test chamber was a transparent Plexiglas container (50 cm × 25 cm × 30 cm) in experiment 1 and a container made of transparent reinforced polyethylene (39 cm × 59 cm × 37 cm) in experiment 2.
Before an experiment, all of the males were acclimated 3-4 times to the experimental chamber for 10 min the first time and then for 5 min in consecutive acclimating sessions.

Acquisition of Sexual Experience
A male was introduced into the experimental chamber and a female was introduced 5 min later. Ultrasounds were recorded during the 5-min period between the introduction of the male and the female to the experimental chamber (precontact anticipatory ultrasonic vocalizations-PVs). The session was conducted until the first ejaculation and ended just after the male resumed copulatory activity. The maximum duration of a single session was 30 min and the session was terminated after this time.

Ultrasound Analysis and Behavior Recording
Visual recording of behavior was made using the Noldus EthoVision system.
Simultaneous to visual recording, ultrasounds were recorded on the same computer using the Metris Sonotrack system.
The microphone was placed 50 cm above the floor during ultrasonic recording. Spectral analysis of ultrasounds was performed using the Sonotrack software. Each ultrasound was analyzed manually using the Sonotrack cursor and was assigned to one of the 14 subtypes according to the classification proposed by Wright et al. (2010) based on its characteristic shape, complexity, and average frequency.
Each session was independently analyzed by two experienced observers to minimize bias in classification. The most inconclusive results were found between the complex and composite categories. For this reason, we decided to modify the classification and distinguish only 13 subtypes. Therefore, the complex and the composite subtypes were combined into one complex/composite category.
In this classification, subtypes of ultrasounds are defined as: complex/composite, trill, flat-trill, trill with jumps, split, step up, step down, multi-step, flat, short, upward ramp, downward ramp, and inverted U. In ambiguous cases, the ultrasounds were listened to at slow speed to classify them correctly.
Subsequently, we counted the total number of ultrasounds and the proportion of their subtypes emitted by each male. The identification of dominant subtypes in the recording was used to create a simplified code characteristic for each individual. Subtypes accounting for at least 10% of the total number of signals were understood to be dominant. Additionally, a low proportion (less than 10%) of the flat vocalization, which is usually frequently emitted by male rats, was considered as a distinguishing feature. In experiment 1, a characteristic vocalization profile for each male was created based on the dominant subtypes visible in all sessions with vocalizations above 9 USV.

Experiment 1
Males (N = 17) were housed 2-3 animals per standard laboratory cage (55.6 cm × 33.4 cm × 19.5 cm). All were sexually naïve at the beginning of the experiment. They acquired sexual experience during four consecutive copulatory sessions. Our previous results showed that during four sessions changes in sexual parameters related to the acquisition of sexual experience are most relevant and stabilized around the fifth session (Bialy et al., 2000).
All of the precontact 50-kHz vocalizations (PVs), quantification of every subtype and their percentage share in the total spectrum were analyzed. The sum of PVs was additionally correlated with copulatory efficiency measured as the sum of sessions when the male achieved ejaculation (during sessions 2-4 after the male had his first sexual experience).
One male was excluded from the experiment due to extremely aggressive behavior.

Experiment 2
Males (N = 20) were housed in the enriched environment as the group contained up to five animals which supported hierarchization within the group. These special larger cages had a vertical structure with wire sidewalls and a solid base (base dimensions 48 cm × 80 cm and 142 cm high) and were equipped with three levels of wooden platforms. The males before the experimental test achieved socio-sexual experience during five copulatory training sessions.
After such training, rats achieved a relatively stable level of sexual performance related to sexual experience (Bialy et al., 2000). During each session, the male was placed individually with a receptive female in one cage. Training sessions were terminated after 30 min. As the indicator of copulatory efficiency, we used the number of sessions when males achieved ejaculation(s) (or lack of them) from all five training sessions. Then we divided each male into one of three groups based on their copulatory history.
Alpha males (N = 8) ejaculated at least during two sessions, beta males (N = 7) ejaculated only once and gamma males (N = 5) never ejaculated during the five copulatory sessions.
One week after the last training session, precontact anticipatory ultrasonic vocalizations were recorded. A male was introduced for 5 min into the familiar experimental/training chamber where odor cues from receptive females were present. The scent stimuli were provided by the earlier placement of three randomly selected receptive females for 5 min into the experimental cage. Prior to the introduction of the male, the females were taken out.
Males that emitted less than 10 ultrasounds were excluded from statistical analysis. This criterion was adopted because, in the case of few ultrasounds, a single vocalization showed a disproportionately high percentage of a given subtype of vocalization, which led to a significant error in the interpretation of the results. This criterion was met by a total of three malesone male from each examined group.

Statistics
Data from the acquisition phase (the number of ultrasounds, mount, intromission, and ejaculation latencies) were analyzed by non-parametric Friedman repeated measure ANOVA and Dunn's post hoc tests (Experiment 1). The number of ultrasounds in Experiment 2 was analyzed by non-parametric Kruskal-Wallis ANOVA and Mann-Whitney for independent groups. Additionally, in experiment 1, the Spearman test was used to correlate the number of the session with ejaculation(s) with the total number of PVs.
A significant correlation (Spearman r = 0.5608, P < 0.05, N = 16) was found between the number of sessions with ejaculation (sessions 2-4) and the total number of PVs emitted by a male during these sessions (Figure 2).

PVs Subtypes
We counted the percentage of different subtypes when a male vocalized 10 times or more in a session. An analysis of USV subtypes was performed and more than 10 vocalizations occurred the most often in three sessions for the same male (six males) and less frequently: in four sessions (five males), in two sessions (four males), and one male vocalized less than 10 times in all four sessions. Those subtypes that appeared at least 10% in the total pool of male vocalizations were marked as the dominant calls. Table 1 shows those subtypes that predominated in all sessions (100%) and those that occurred in at least two sessions out of three or four (50-75%) sessions and in only one session (25%). The complex/composite and flat (or short or upward ramp) subtypes occurred most frequently. The complex/composite and flat (or short or upward ramp) profile was typical to anticipatory precontact vocalizations regardless of the male's copulatory efficiency as measured by the number of ejaculations achieved during sessions 2-4.
Subtypes of PVs did not differ as a function of the acquisition of sexual experience. The percentage of different subtypes showed a similar level during all four copulatory sessions (Figure 3).

Number and Subtypes of PVs
The number of PVs differed significantly between the alpha, beta, and gamma groups [H(2), N = 17] = 6.0396, P < 0.05. The alpha group vocalized more frequently compared with the gamma group (P < 0.05). The beta group did not differ significantly compared with the alpha and gamma groups (Figure 4).
There were no statistically significant differences in the percentage of selected subtypes in the total pool of ultrasounds emitted by males when these three groups were compared (Figure 5). Table 2 shows the main ultrasonic subtypes emitted by each male.

DISCUSSION
The results indicate a significant positive relationship between the total number of anticipatory precontact vocalizations and the number of sessions with achieved ejaculation(s). It reveals FIGURE 1 | The means number (SE) of PVs during the acquisition of sexual experience. Numbers on the Y axis represent the means value of PVs emitted during subsequent sessions. Statistically significant differences compared with session 1 ( * P < 0.05 and * * P < 0.01).
FIGURE 2 | Total number of PVs with relation to the sum of ejaculations during sessions 2-4. Numbers on the X-axis represent the sum of ejaculation from 0 to 3 achieved by the male during all three sessions 2-4. 0-means that male did not achieve any ejaculation during three sessions, respectively, 1-only one ejaculation during three sessions, 2-two ejaculations, and 3-ejaculated in each session. Numbers on the Y-axis represent the sum of PVs emitted by the male during sessions 2-4. The diamond represents an individual male. r = 0.5608, P < 0.05.
an association between the level of socio-sexual motivation and the expression of PVs. Shortening of mount latency-parameter related to sexual motivation (Hull et al., 2002) during the acquisition of sexual experience additionally supports the relation between PVs and sexual motivation.
In two independent groups of rats (experiments 1 and 2), similar patterns of emitted ultrasounds were observed. In all of the males, regardless of their level of sexual motivation or sexual experience, the dominant signals were the frequency modulated complex/composite (CC) signals and the unmodulated flat (F)  and, less frequently, the slightly modulated short (Sh) and upward ramp (UR) ultrasounds appeared as co-occurring subtypes. Trills were emitted on an elevated level in experiment 1 only.
Other signals were seen sporadically and rarely exceeded 10% of the total PVs.

PVs and Sexual Activity
Anticipatory precontact vocalizations depend on the level of sexual activity and reward value of the contacts. In experiment 1, the number of PVs increased during the acquisition of sexual experience. Moreover, the total number of PVs relates to the number of ejaculations. Ejaculation has the highest reward value compared with other elements of copulatory behavior (Tenk et al., 2009) and, in this context, PVs are related to the reward value of socio-sexual contacts. In experiment 2, the active males vocalized on a significantly higher level compared with much less active or sexually inactive males.
A more detailed explanation of the processes related to anticipatory ultrasonic precontact vocalizations follows. The number of PVs depends on the acquisition of sexual experience, conditioning to odor and background cues, the reward value of contact, NMDA (Bialy et al., 2000), and D1 receptor activity (Bialy et al., 2010).
Blocking of the NMDA receptor (receptor important in neuronal plasticity processes) inhibited changes in PVs during the acquisition of sexual experience but have no significant effect on PVs in sexually experienced males (Bialy et al., 2000). Acquisition of sexual experience and PVs are also related to D1 receptor activity. Surprisingly, peripheral repeated administration of both: antagonist and agonist of D1 receptor inhibited changes in PVs during the acquisition   phase (Bialy et al., 2010) but, D1 receptor agonist had a minor effect on PVs in sexually experienced males (Beck et al., 2002). Repeated D1 receptor agonist injections into the nucleus accumbens similarly diminished PVs during the acquisition phase and this effect prolonged at least 4 weeks.  Previously, were identified the specific, dopaminergic neuronal group in the nucleus accumbens responding either to appetitive unconditioned and conditioned stimuli that regulating, reward related behavior (Fiorillo et al., 2003). Thus, probably any significant disturbance of their activity: inhibition or overstimulation can be responsible for such attenuation of the acquisition phase and diminishing PVs.
In this context, PVs reflect learning and emotional memory processes. The number of PVs seems to be related to the level of general arousal and sexual motivation rather than to sexual arousal (Bialy et al., 2019a). In addition, social motivation should also be taken into account, as contacts with an anestrus female also has a triggering effect on PVs although at a significantly lower level (Bialy et al., 2000). We can assume that in the present experiments, PVs in low sexually active males were related to social motivation and general arousal rather than sexual motivation. All of these data seem to be consistent and repeatable, enabling the changing number of PVs to be used as a parameter for the functioning of the memory circuits related to the reward system.

PVs and Individualization
We observed individual differences in the percentage of nondominant (less than 10% of the total spectrum) subtypes as well as in the occasionally very low level of the flat subtype. The composite signals contained within the CC group consisted of a modulated sound of the complex subtype with a directly accompanying sound of the short subtype. The distinction between the two ultrasounds was often only possible after listening to the signal. In some studies, this type of signal is listed as a "dual-type." However, due to the high subjectivity in the evaluation between different observers and the inability to listen on different recording systems other than our system, we decided on a common category. The composite subtype calls, consisting of calls other than complex and short components, occurred very rarely. Furthermore, the individual specificity of vocalizations was manifested in different proportions in the short, upward ramp, and flat subtypes. In general, the presented form of individualization seems to be subtle and not easy to use as a behavioral parameter.
This limitation is also related to the fact that a vast majority of inter-individual differentiation can occur within the so-called frequency-modulated non-trill group (split, step down, step up, multistep, downward, and inverted U) Despite their lack of visibility in the profiles after applying our quantitative criterion, variation within these types exists. Each of them is characterized by its specific shape, although their acoustic parameters, such as duration, mean frequencies, or modulation ranges, are heterogeneous. In summary, the observed individualization additionally occurs within the types of ultrasound that occur sporadically, in their exact parameters, and their potential combinations (like syllables).
Nevertheless, given the existence of the presented pattern, we can assume that it reflects a relatively constant emotional state. In this context, large, comparatively easy to visualize deviations from it may indicate significant central changes. The question of a behaviorally relevant reason for the existence of certain forms of individualization remains open.

USV Profiles
The main question, and aim of our experiments, was to answer the question: do changes in the level of PVs related to elevated general arousal and sexual motivation correspond to changes in the profile of the ultrasound subtypes. Qualitative analysis of ultrasound subtypes revealed an existing characteristic vocalization pattern with complex/composite subtypes combined with flat, short, and upward ramp subtypes dominating in the recordings. This model of vocalization appears similar in both sexually inactive and active males. Moreover, it seems constant during consecutive sessions of sexual experience acquisition and sexual motivation influenced the number of PVs but not the subtype profile. In experiment 2, more escalated social hierarchization was evident in the form of an increase in the non-copulatory number of rats (the gamma group). Regardless of this, the pattern of the profiles remained similar.
In our experiments, trills were visible in some sexually moderately active males. They did not change as a function of sexual experience. Probably they are related to the general arousal level but not sexual motivation per se. This hypothesis can be supported by the results of experiments with the repeated administration of psychostimulant drug agents (amphetamine and derivatives, apomorphine, cocaine) which evokes the expression of the trills subtype (Wright et al., 2012(Wright et al., , 2013Mulvihill and Brudzynski, 2018). However, the biological context of the expression of this type of vocalization is not entirely clear. For example, Mahler et al. (2013) observed no difference in the number of trill emissions during the methamphetamine self-administration period and the subsequent extinction period, obtaining an increase only after re-exposure of the animals to the drug.
Furthermore, many authors, due to the simplicity of ultrasonic classification, often cluster all frequency modulated subtypes (Taracha et al., 2012;Brenes and Schwarting, 2015;Mulvihill and Brudzynski, 2018). The lack of correlation between sexual motivation and the profile of PVs could indicate that the specific type of frequency modulated calls presented in our experiments is an innate pattern of general anticipatory activity. In this context, the degree of modulation could represent a continuous behavioral spectrum, with complex/composite subtype domination associated with more physiological situations and the trills subtype domination reflecting highly aroused states. Further research is needed to confirm this hypothesis. Also, the neurobehavioral context of flat, unmodulated subtype is not clear. In these experiments, we observed flat call as one of the dominant types in most rats. However, in few animals, we observed a noticeable decrease of flats calls in profiles (Tables 1, 2). This reduction wasn't related to copulatory status (activity), indicating the involvement of a factor other than socio-sexual motivation. Recently, some authors proposed that this type can be associated with social-coordinating function and alimentary appetitive reactions, including the food approach (Brudzynski, 2021). According to our results, additional data are needed to clarify this issue.

Further Implications
We have shown that the pattern of physiological anticipatory vocalization was relatively constant and persisted from one experimental session to the next. It is worth noting, however, that vocalization can change significantly in some situations. We observed rapid changes in vocalization occurring during a single session with characteristic transitions between the 50-kHz and 22-kHz frequencies during the frustrated states (Bialy et al., 2019b). Using an analogous non-contact model, changes in the vocalization profile across consecutive sessions have also been demonstrated in mice (Zala et al., 2020). Similarly, significant changes during consecutive experimental sessions have been shown repeatedly with addiction models, e.g., morphine (Covington and Miczek, 2003;Hamed and Kursa, 2018), cocaine (Ma et al., 2010), and methamphetamine (Mahler et al., 2013). Another example of such an application is the change in vocalization profiles that have been observed in induced Parkinson's disease in rats where ultrasound flattening occurs as the disease progresses (Ciucci et al., 2009). Individual vocalization and own vs. foreign recognition can be disturbed in schizophrenic-like symptoms in the rat model (Matsumoto et al., 2016). Furthermore, changes in the vocalization profile occur in the rat model of affective disorders and psychotic states (Nikiforuk et al., 2013;Wendler et al., 2019;Wöhr, 2021) as well as in suggested models of autism (Caruso et al., 2020).

CONCLUSION
We proposed a relatively simple method to discriminate individual characters of ultrasonic vocalization in rats based on dominant subtype ultrasonic vocalizations. Our results indicate a persistent similar vocalization pattern during anticipatory behavior, regardless of the level of socio-sexual motivation or experience. The number of ultrasonic vocalizations but not the different subtypes seems to be the most related to sexual motivation.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

ETHICS STATEMENT
The animal study was reviewed and approved by First Local Ethical Committee in Warsaw.

AUTHOR CONTRIBUTIONS
WB-R and MB designed and performed experiments, analyzed data, and prepared the manuscript. MR analyzed USV data in experiments 2. All authors contributed to the article and approved the submitted version.