Skip to main content

ORIGINAL RESEARCH article

Front. Behav. Neurosci., 20 August 2021
Sec. Motivation and Reward
This article is part of the Research Topic Motivation and Reward - Editors' Pick 2022 View all 10 articles

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

  • Laboratory of Centre for Preclinical Research, Department of Experimental and Clinical Physiology, Medical University of Warsaw, Warsaw, Poland

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.

Introduction

Ultrasonic vocalization is one of the most intensively studied components of social behavior in rodents. In rats, two main types of ultrasonic vocalization (USV) can be distinguished based on the dominant frequency of the signal: 22-kHz and 50-kHz (Sales and Pye, 1974; Barfield et al., 1979; Brudzynski, 2015, 2021). The first is characterized by long-lasting flat calls (up to 3 s in duration) with a relatively narrow frequency range around 22-kHz (Barfield and Geyer, 1972). In the frustration state (situation of the absence of expected appetitive reinforcement), these calls can be frequency modulated (20–35 kHz) with the preceding element at a frequency of about 45-kHz (Geyer et al., 1978; Bialy et al., 2019b). The so-called 50-kHz USV includes shorter ultrasounds (typically lasting up to several tens of milliseconds) with frequencies from 30–35 to 80-kHz. Signals of such calls, after undergoing Fourier transformations, display various shapes with 14 different subtypes (Wright et al., 2010). The rats emit 50-kHz USV during different elements of behavior related to their high arousal states (Bell, 1974; Berz et al., 2021) such as: socio-sexual interactions including copulation (Barfield et al., 1979; Bialy et al., 2000; Burgdorf et al., 2008), fighting (Sales, 1972; Burke et al., 2017), playing (Reinhold et al., 2019), and even tickling by the experimenter (Burgdorf et al., 2005; Panksepp, 2005). Additionally, pharmacologically induced high levels of general arousal [related to movement activity and sensory sensitization associated with increased level of wakefulness during activation of the gigantocellular reticular nucleus and associated structures (Pfaff, 2017)] result in the expression of 50-kHz vocalizations with a strong positive correlation between the number of 50-kHz ultrasounds emitted and the level of activation of the dopaminergic and noradrenergic pathways (Brudzynski, 2007, 2015; Wright et al., 2012; Simola, 2015; Hamed et al., 2016; Simola and Costa, 2018; Kuchniak et al., 2019).

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., 1996, 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 socio-sexual 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?

Materials and Methods

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, 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 females (N = 20, 10 in the first and 10 in the second experiment), were housed 2–3 animals per standard laboratory cage (55.6 cm × 33.4 cm × 19.5 cm).

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 males—one 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.

Results

Experiment 1

Number of PVs

The number of PVs (Figure 1) increased significantly in males (N = 16) during four consecutive sessions of the acquisition of sexual experience (Fr = 12.816, P < 0.01 with statistically significant differences between sessions 1 and 3 (P < 0.05), and sessions 1 and 4 (P < 0.01). Mount latency decreased significantly (Fr = 17.591, P < 0.001) with statistically significant differences between sessions 1 and 2 (P < 0.05), sessions 1 and 3 (P < 0.01) and 1 and 4 (P < 0.01. Intromission and ejaculation latencies display no significant differences (P = 0.0815 and P = 0.077 respectively).

FIGURE 1
www.frontiersin.org

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).

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).

FIGURE 2
www.frontiersin.org

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.

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.

TABLE 1
www.frontiersin.org

Table 1. PVs and ultrasonic subtypes during consecutive sessions.

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).

FIGURE 3
www.frontiersin.org

Figure 3. Means percentage of different ultrasonic subtypes during the acquisition of sexual experience. Subtypes are represent as: complex and composite (CC), trill, flat-trill and trill with jumps (Trills), split, step up (SU), step down (SD), multi-step (MS), flat, short, upward ramp (UR), downward ramp (DR), and inverted U (IU). Subtypes on the X-axis are shown from left to right of each column as sessions 1, 2, 3, and 4.

Experiment 2

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).

FIGURE 4
www.frontiersin.org

Figure 4. The means number (SE) of PVs emitted by alpha, beta, and gamma groups. 1—Alpha group, 2—Beta group, and 3—Gamma group. Statistically significant differences (P < 0.05) compared with the alpha group.

The three most frequent subtypes of ultrasounds detected were: complex/composite (alpha 39.9%, beta 42.4%, and gamma 38% of total ultrasounds), flat (alpha 23.1%, beta 23.8%, gamma 24% of total ultrasounds) and upward ramp (alpha 10.2%, beta 11.5%, gamma 9.2% of total ultrasounds).

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.

FIGURE 5
www.frontiersin.org

Figure 5. The dominant means percentage of different subtypes in the alpha, beta and gamma groups. Subtypes are shown from left to right as: complex and composite (CC), trill, flat-trill and trill with jumps (Trills), split, step up (SU), step down (SD), multi-step (MS), flat, short, upward ramp (UR), downward ramp (DR), and inverted U (IU).

TABLE 2
www.frontiersin.org

Table 2. The individual profiles of emitted ultrasounds.

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 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 non-dominant (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, 2013; Mulvihill 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.

Funding

This work was supported by the Medical University of Warsaw and the MUW Young Researcher Grant: 1MA/5/M/MB/N/20/20. The study was carried out with the use of the CePT infrastructure financed by the European Union—the European Regional Development Fund within the Operational Program “Innovative Economy” for 2007–2013.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Arriaga, G., and Jarvis, E. D. (2013). Mouse vocal communication system: are ultrasounds learned or innate? Brain Lang. 124, 96–116. doi: 10.1016/j.bandl.2012.10.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Asaba, A., Hattori, T., Mogi, K., and Kikusui, T. (2014). Sexual attractiveness of male chemicals and vocalizations in mice. Front. Neurosci. 8:231. doi: 10.3389/fnins.2014.00231

PubMed Abstract | CrossRef Full Text | Google Scholar

Barfield, R. J., and Geyer, L. A. (1972). Sexual behavior: ultrasonic postejaculatory song of the male rat. Science 176, 1349–1350. doi: 10.1126/science.176.4041.1349

PubMed Abstract | CrossRef Full Text | Google Scholar

Barfield, R. J., Auerbach, P., Geyer, L. A., and McIntosh, T. K. (1979). Ultrasonic Vocalizations in Rat Sexual Behavior. Am. Zool. 19, 469–480. doi: 10.1093/icb/19.2.469

PubMed Abstract | CrossRef Full Text | Google Scholar

Beck, J., Biały, M., and Kostowski, W. (2002). Effects of D(1) receptor agonist SKF 38393 on male rat sexual behavior and postcopulatory departure in the goal compartment-runway paradigm. Physiol. Behav. 76, 91–97. doi: 10.1016/s0031-9384(02)00678-9

CrossRef Full Text | Google Scholar

Bell, R. W. (1974). Ultrasounds in small rodents: arousal-produced and arousal-producing. Dev. Psychobiol. 7, 39–42. doi: 10.1002/dev.420070107

PubMed Abstract | CrossRef Full Text | Google Scholar

Berz, A., Pasquini, de Souza, C., Wöhr, M., and Schwarting, R. K. W. (2021). Limited generalizability, pharmacological modulation, and state-dependency of habituation towards pro-social 50-kHz calls in rats. Science 24:102426. doi: 10.1016/j.isci.2021.102426

PubMed Abstract | CrossRef Full Text | Google Scholar

Bialy, M., Beck, J., Abramczyk, P., Trzebskj, A., and Przybylski, J. (1996). Sexual behavior in male rats after nitric oxide synthesis inhibition. Physiol. Behav. 60, 139–143. doi: 10.1016/0031-9384(95)02272-4

CrossRef Full Text | Google Scholar

Bialy, M., Bogacki-Rychlik, W., Kasarello, K., Nikolaev, E., and Sajdel-Sulkowska, E. M. (2016). Modulation of 22-khz postejaculatory vocalizations by conditioning to new place: Evidence for expression of a positive emotional state. Behav. Neurosci. 130, 415–421. doi: 10.1037/bne0000153

PubMed Abstract | CrossRef Full Text | Google Scholar

Bialy, M., Bogacki-Rychlik, W., Przybylski, J., and Zera, T. (2019a). The Sexual Motivation of Male Rats as a Tool in Animal Models of Human Health Disorders. Front. Behav. Neurosci. 13:257. doi: 10.3389/fnbeh.2019.00257

PubMed Abstract | CrossRef Full Text | Google Scholar

Bialy, M., Kalata, U., Nikolaev-Diak, A., and Nikolaev, E. (2010). D1 receptors involved in the acquisition of sexual experience in male rats. Behav. Brain Res. 206, 166–176. doi: 10.1016/j.bbr.2009.09.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Bialy, M., Podobinska, M., Barski, J., Bogacki-Rychlik, W., and Sajdel-Sulkowska, E. M. (2019b). Distinct classes of low frequency ultrasonic vocalizations in rats during sexual interactions relate to different emotional states. Acta Neurobiol. Exp. 79, 1–12. doi: 10.21307/ane-2019-001

CrossRef Full Text | Google Scholar

Bialy, M., Rydz, M., and Kaczmarek, L. (2000). Precontact 50-kHz vocalizations in male rats during acquisition of sexual experience. Behav. Neurosci. 114, 983–990. doi: 10.1037/0735-7044.114.5.983

PubMed Abstract | CrossRef Full Text | Google Scholar

Blumberg, M. S., and Alberts, J. R. (1991). On the significance of similarities between ultrasonic vocalizations on infant and adult rats. Neuroscie. Biobehav. Rev. 15, 383–390. doi: 10.1016/S0149-7634(05)80031-4

CrossRef Full Text | Google Scholar

Brenes, J. C., and Schwarting, R. K. W. (2015). Individual differences in anticipatory activity to food rewards predict cue-induced appetitive 50-kHz calls in rats. Physiol. Behav. 149, 107–118. doi: 10.1016/j.physbeh.2015.05.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Brudzynski, S. M. (2007). Ultrasonic calls of rats as indicator variables of negative or positive states: acetylcholine-dopamine interaction and acoustic coding. Behav. Brain Res. 182, 261–273. doi: 10.1016/j.bbr.2007.03.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Brudzynski, S. M. (2015). Pharmacology of Ultrasonic Vocalizations in adult Rats: Significance, Call Classification and Neural Substrate. Curr. Neuropharmacol. 13, 180–192. doi: 10.2174/1570159x13999150210141444

PubMed Abstract | CrossRef Full Text | Google Scholar

Brudzynski, S. M. (2021). Biological Functions of Rat Ultrasonic Vocalizations, Arousal Mechanisms, and Call Initiation. Brain Sci. 11:605. doi: 10.3390/brainsci11050605

PubMed Abstract | CrossRef Full Text | Google Scholar

Burgdorf, J., Kroes, R. A., Moskal, J. R., Pfaus, J. G., Brudzynski, S. M., and Panksepp, J. (2008). Ultrasonic vocalizations of rats (Rattus norvegicus) during mating, play, and aggression: Behavioral concomitants, relationship to reward, and self-administration of playback. J. Comp. Psychol. 122, 357–367. doi: 10.1037/a0012889

PubMed Abstract | CrossRef Full Text | Google Scholar

Burgdorf, J., Panksepp, J., Brudzynski, S. M., Kroes, R., and Moskal, J. R. (2005). Breeding for 50-kHz positive affective vocalization in rats. Behav. Genet. 35, 67–72. doi: 10.1007/s10519-004-0856-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Burke, C. J., Kisko, T. M., Pellis, S. M., and Euston, D. R. (2017). Avoiding escalation from play to aggression in adult male rats: The role of ultrasonic calls. Behav. Process. 144, 72–81. doi: 10.1016/j.beproc.2017.09.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Caruso, A., Ricceri, L., and Scattoni, M. L. (2020). Ultrasonic vocalizations as a fundamental tool for early and adult behavioral phenotyping of Autism Spectrum Disorder rodent models. Neurosci. Biobehav. Rev. 116, 31–43. doi: 10.1016/j.neubiorev.2020.06.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Chabout, J., Sarkar, A., Dunson, D. B., and Jarvis, E. D. (2015). Male mice song syntax depends on social contexts and influences female preferences. Front. Behav. Neurosci. 9:76. doi: 10.3389/fnbeh.2015.00076

PubMed Abstract | CrossRef Full Text | Google Scholar

Ciucci, M. R., Ahrens, A. M., Ma, S. T., Kane, J. R., Windham, E. B., Woodlee, M. T., et al. (2009). Reduction of dopamine synaptic activity: degradation of 50-kHz ultrasonic vocalization in rats. Behav. Neurosci. 123, 328–336. doi: 10.1037/a0014593

PubMed Abstract | CrossRef Full Text | Google Scholar

Covington, H. E., and Miczek, K. A. (2003). Vocalizations during withdrawal from opiates and cocaine: possible expressions of affective distress. Eur. J. Pharmacol. 467, 1–13. doi: 10.1016/S0014-2999(03)01558-9

CrossRef Full Text | Google Scholar

Fiorillo, C. D., Tobler, P. N., and Schultz, W. (2003). Discrete coding of reward probability and uncertainty by dopamine neurons. Science 299, 1898–1902. doi: 10.1126/science.1077349

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu, X. W., and Brudzynski, S. M. (1994). High-frequency ultrasonic vocalization induced by intracerebral glutamate in rats. Pharmacol. Biochem. Behav. 49, 835–841. doi: 10.1016/0091-3057(94)90231-3

CrossRef Full Text | Google Scholar

Gao, S. C., Wei, Y. C., Wang, S. R., and Xu, X. H. (2019). Medial Preoptic Area Modulates Courtship Ultrasonic Vocalization in Adult Male Mice. Neurosci. Bull. 35, 697–708. doi: 10.1007/s12264-019-00365-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Geyer, L. A., Barfield, R. J., and McIntosh, T. K. (1978). Influence of gonadal hormones and sexual behavior on ultrasonic vocalization in rats: II. Treatment of males. J. Comparat. Physiol. Psychol. 92, 447–456. doi: 10.1037/h0077487

CrossRef Full Text | Google Scholar

Hamed, A., and Kursa, M. B. (2018). Inter-individual differences in serotonin and glutamate co-transmission reflect differentiation in context-induced conditioned 50-kHz USVs response after morphine withdrawal. Brain Struct. Funct. 223, 3149–3167. doi: 10.1007/s00429-018-1683-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Hamed, A., Daszczuk, P., Kursa, M. B., Turzyńska, D., Sobolewska, A., Lehner, M., et al. (2016). Non-parametric analysis of neurochemical effects and Arc expression in amphetamine-induced 50-kHz ultrasonic vocalization. Behav. Brain Res. 312, 174–185. doi: 10.1016/j.bbr.2016.05.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Holy, T. E., and Guo, Z. (2005). Ultrasonic Songs of Male Mice. PLoS Biol. 3:e386. doi: 10.1371/journal.pbio.0030386

PubMed Abstract | CrossRef Full Text | Google Scholar

Hull, E. M., Meisel, R. L., and Sachs, B. D. (2002). “1 - Male Sexual Behavior,” in Hormones, Brain and Behavior, eds D. W. Pfaff, A. P. Arnold, S. E. Fahrbach, A. M. Etgen, and R. T. Rubin (San Diego: Academic Press), 3–137.

Google Scholar

Jelen, P., Soltysik, S., and Zagrodzka, J. (2003). 22-kHz ultrasonic vocalization in rats as an index of anxiety but not fear: behavioral and pharmacological modulation of affective state. Behav. Brain Res. 141, 63–72. doi: 10.1016/s0166-4328(02)00321-2

CrossRef Full Text | Google Scholar

Karigo, T., Kennedy, A., Yang, B., Liu, M., Tai, D., Wahle, I. A., et al. (2021). Distinct hypothalamic control of same- and opposite-sex mounting behaviour in mice. Nature 589, 258–263. doi: 10.1038/s41586-020-2995-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuchniak, K., Wyszogrodzka, E., Chrapusta, S. J., Czarna, M., Michalak, M., Płaźnik, A., et al. (2019). Using anticipatory and drug-evoked appetitive ultrasonic vocalization for monitoring the rewarding effect of amphetamine in a rat model of drug self-administration. Behav. Brain Res. 376:112187. doi: 10.1016/j.bbr.2019.112187

PubMed Abstract | CrossRef Full Text | Google Scholar

Larsson, K. (1956). Conditioning and sexual behavior in the male albino rat. Stockholm: Almqvist & Wiksell.

Google Scholar

Larsson, K. (1967). Testicular hormone and developmental changes in mating behavior of the male rat. J. Comparat. Physiol. Psychol. 63, 223–230. doi: 10.1037/h0024358

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, S. T., Maier, E. Y., Ahrens, A. M., Schallert, T., and Duvauchelle, C. L. (2010). Repeated intravenous cocaine experience: development and escalation of pre-drug anticipatory 50-kHz ultrasonic vocalizations in rats. Behav. Brain Res. 212, 109–114. doi: 10.1016/j.bbr.2010.04.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Mahler, S. V., Moorman, D. E., Feltenstein, M. W., Cox, B. M., Ogburn, K. B., Bachar, M., et al. (2013). A rodent “self-report” measure of methamphetamine craving? Rat ultrasonic vocalizations during methamphetamine self-administration, extinction, and reinstatement. Behav. Brain Res. 236, 78–89. doi: 10.1016/j.bbr.2012.08.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Matsumoto, J., Nishimaru, H., Takamura, Y., Urakawa, S., Ono, T., and Nishijo, H. (2016). Amygdalar Auditory Neurons Contribute to Self-Other Distinction during Ultrasonic Social Vocalization in Rats. Front. Neurosci. 10:399. doi: 10.3389/fnins.2016.00399

PubMed Abstract | CrossRef Full Text | Google Scholar

Mulvihill, K. G., and Brudzynski, S. M. (2018). Non-pharmacological induction of rat 50 kHz ultrasonic vocalization: Social and non-social contexts differentially induce 50 kHz call subtypes. Physiol. Behav. 196, 200–207. doi: 10.1016/j.physbeh.2018.09.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Nikiforuk, A., Kos, T., Fijał, K., Hołuj, M., Rafa, D., and Popik, P. (2013). Effects of the Selective 5-HT7 Receptor Antagonist SB-269970 and Amisulpride on Ketamine-Induced Schizophrenia-like Deficits in Rats. PLoS One 8:e66695. doi: 10.1371/journal.pone.0066695

PubMed Abstract | CrossRef Full Text | Google Scholar

Panksepp, J. (2005). Beyond a Joke: From Animal Laughter to Human Joy? Science 308, 62–63. doi: 10.1126/science.1112066

PubMed Abstract | CrossRef Full Text | Google Scholar

Pfaff, D. (ed.) (2017). “Central nervous system arousal fueling instinctive behaviors,” in How the Vertebrate Brain Regulates Behavior, (Cambridge: Harvard University Press), 215–226. doi: 10.4159/9780674978751-009

CrossRef Full Text | Google Scholar

Reinhold, A. S., Sanguinetti-Scheck, J. I., Hartmann, K., and Brecht, M. (2019). Behavioral and neural correlates of hide-and-seek in rats. Science 365, 1180–1183. doi: 10.1126/science.aax4705

PubMed Abstract | CrossRef Full Text | Google Scholar

Sales, G. D. (1972). Ultrasound and aggressive behaviour in rats and other small mammals. Anim. Behav. 20, 88–100. doi: 10.1016/S0003-3472(72)80177-5

CrossRef Full Text | Google Scholar

Sales, G. D., and Pye, J. D. (1974). Ultrasonic communication by animals. London: Chapman and Hall.

Google Scholar

Simola, N. (2015). Rat Ultrasonic Vocalizations and Behavioral Neuropharmacology: From the Screening of Drugs to the Study of Disease. Curr. Neuropharmacol. 13, 164–179. doi: 10.2174/1570159x13999150318113800

PubMed Abstract | CrossRef Full Text | Google Scholar

Simola, N., and Costa, G. (2018). Emission of categorized 50-kHz ultrasonic vocalizations in rats repeatedly treated with amphetamine or apomorphine: Possible relevance to drug-induced modifications in the emotional state. Behav. Brain Res. 347, 88–98. doi: 10.1016/j.bbr.2018.02.041

PubMed Abstract | CrossRef Full Text | Google Scholar

Simola, N., and Granon, S. (2019). Ultrasonic vocalizations as a tool in studying emotional states in rodent models of social behavior and brain disease. Neuropharmacology 159:107420. doi: 10.1016/j.neuropharm.2018.11.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Taracha, E., Hamed, A., Krza̧ścik, P., Lehner, M., Skórzewska, A., Płaźnik, A., et al. (2012). Inter-individual diversity and intra-individual stability of amphetamine-induced sensitization of frequency-modulated 50-kHz vocalization in Sprague-Dawley rats. Psychopharmacology 222, 619–632. doi: 10.1007/s00213-012-2658-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Tenk, C. M., Wilson, H., Zhang, Q., Pitchers, K. K., and Coolen, L. M. (2009). Sexual reward in male rats: Effects of sexual experience on conditioned place preferences associated with ejaculation and intromissions. Hormones Behav. 55, 93–97. doi: 10.1016/j.yhbeh.2008.08.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Wendler, E., de Souza, C. P., Dornellas, A. P. S., Santos, L. E., Ferreira, S. T., Galduróz, J. C. F., et al. (2019). Mania-like elevated mood in rats: Enhanced 50-kHz ultrasonic vocalizations after sleep deprivation. Prog. Neuropsychophar. Biol. Psychiatry 88, 142–150. doi: 10.1016/j.pnpbp.2018.07.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Wöhr, M. (2021). Measuring mania-like elevated mood through amphetamine-induced 50-kHz ultrasonic vocalizations in rats. Br. J. Pharmacol. [Preprint]. doi: 10.1111/bph.15487

PubMed Abstract | CrossRef Full Text | Google Scholar

Wright, J. M., Dobosiewicz, M. R. S., and Clarke, P. B. S. (2012). α- and β-Adrenergic Receptors Differentially Modulate the Emission of Spontaneous and Amphetamine-Induced 50-kHz Ultrasonic Vocalizations in Adult Rats. Neuropsychopharmacology 37, 808–821. doi: 10.1038/npp.2011.258

PubMed Abstract | CrossRef Full Text | Google Scholar

Wright, J. M., Dobosiewicz, M. R. S., and Clarke, P. B. S. (2013). The role of dopaminergic transmission through D1-like and D2-like receptors in amphetamine-induced rat ultrasonic vocalizations. Psychopharmacology 225, 853–868. doi: 10.1007/s00213-012-2871-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Wright, J. M., Gourdon, J. C., and Clarke, P. B. (2010). Identification of multiple call categories within the rich repertoire of adult rat 50-kHz ultrasonic vocalizations: effects of amphetamine and social context. Psychopharmacology 211, 1–13. doi: 10.1007/s00213-010-1859-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Zala, S. M., Nicolakis, D., Marconi, M. A., Noll, A., Ruf, T., Balazs, P., et al. (2020). Primed to vocalize: Wild-derived male house mice increase vocalization rate and diversity after a previous encounter with a female. PLoS One 15:e0242959. doi: 10.1371/journal.pone.0242959

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: 50-kHz ultrasonic vocalizations, anticipatory behavior, male rats, sexual experience, sexual motivation, ultrasound subtypes

Citation: Bogacki-Rychlik W, Rolf M and Bialy M (2021) Anticipatory 50-kHz Precontact Ultrasonic Vocalizations and Sexual Motivation: Characteristic Pattern of Ultrasound Subtypes in an Individual Analyzed Profile. Front. Behav. Neurosci. 15:722456. doi: 10.3389/fnbeh.2021.722456

Received: 08 June 2021; Accepted: 03 August 2021;
Published: 20 August 2021.

Edited by:

Fabrizio Sanna, University of Cagliari, Italy

Reviewed by:

Nicola Simola, University of Cagliari, Italy
Wendy Portillo, National Autonomous University of Mexico, Mexico

Copyright © 2021 Bogacki-Rychlik, Rolf and Bialy. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Michal Bialy, michal.bialy@wum.edu.pl; michalbialy@yahoo.com

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.