{'fraction_non_alphanumeric': 0.034103410341034104, 'fraction_numerical': 0.024525982009965704, 'mean_word_length': 4.952812018489984, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 8, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Unusual Repertoire of Vocalizations in the BTBR T+tf/J Mouse Model of Autism

BTBR T+ tf/J (BTBR) is an inbred mouse strain that displays social abnormalities and repetitive behaviors analogous to the first and third diagnostic symptoms of autism. Here we investigate ultrasonic vocalizations in BTBR, to address the second diagnostic symptom of autism, communication deficits. As compared to the commonly used C57BL/6J (B6) strain, BTBR pups called more loudly and more frequently when separated from their mothers and siblings. Detailed analysis of ten categories of calls revealed an unusual pattern in BTBR as compared to B6. BTBR emitted high levels of harmonics, two-syllable, and composite calls, but minimal numbers of chevron-shaped syllables, upward, downward, and short calls. Because body weights were higher in BTBR than B6 pups, one possible explanation was that larger thoracic size was responsible for the louder calls and different distribution of syllable categories. To test this possibility, we recorded separation calls from FVB/NJ, a strain with body weights similar to BTBR, and 129X1/SvJ, a strain with body weights similar to B6. BTBR remained the outlier on number of calls, displaying low numbers of complex, upward, chevron, short, and frequency steps calls, along with high harmonics and composites. Further, developmental milestones and growth rates were accelerated in BTBR, indicating an unusual neurodevelopmental trajectory. Overall, our findings demonstrate strain-specific patterns of ultrasonic calls that may represent different lexicons, or innate variations in complex vocal repertoires, in genetically distinct strains of mice. Particularly intriguing is the unusual pattern of vocalizations and the more frequent, loud harmonics evident in the BTBR mouse model of autism that may resemble the atypical vocalizations seen in some autistic infants.

Abnormal reciprocal social interactions and communication deficits are two of the three diagnostic symptoms of autism [25]. We hypothesize that USVs may be a measure of social communication in mice [15,[26][27][28]. Reduced or unusual USVs in mice may offer a useful assay with reasonable face validity to the second diagnostic symptom of autism, impaired communication. A standard test for vocalizations in mice, the ultrasonic distress call of pups separated from the mother or removed from the nest [29][30][31], has been used to investigate the number of calls emitted by separated pups in mouse models of autism spectrum disorders [15,[32][33][34][35][36][37][38][39].
We previously reported on an inbred strain of mice, BTBR T+tf/J (BTBR), that displays several traits relevant to autism, including reduced social approach in adults, reduced reciprocal social interactions in juveniles and adults, decreased social transmission of food preference, and high repetitive self grooming when compared to C57BL/6J (B6), a standard inbred strain commonly used in behavioral genetics [40][41][42][43][44]. To address the hypothesis that BTBR also display communication deficits, and that autism-like phenotypes in mice can be detected at early developmental stages, we analyzed neonatal USV patterns (Cohort 1) and behavioral and somatic development (Cohort 2) in BTBR mice throughout the first two weeks of postnatal life. The primary goal of the present studies was to detect any unusual components of vocalizations in BTBR at infant stages, relevant to the absence of crying, and the unusual guttural grunts and squeals, reported for some babies that were later diagnosed with autism [45,46].
To define ''unusual'' BTBR vocalization patterns in the absence of normative sonographs of USVs across inbred strains of mice, we compared USVs in BTBR versus three control strains typically used in behavioral genetics and as backgrounds for breeding genetically modified mice. B6 represented a widely-used standard control with high levels of social behaviors [47] and low USV calling rate [10,[48][49][50][51][52]. FVB/NJ was selected as a strain with high body weights similar to BTBR, to control for the possibility that larger lung capacity could be responsible for unusual vocalizations. In addition, FVB/NJ displays high levels of sociability in the social approach apparatus, comparable to B6. 129X1/SvJ (129X1) was selected as genetically close to BTBR on the mouse phylogenetic tree [53] but with body weights similar to B6.
Our findings demonstrate a comparable use of ten different categories of calls that may represent a similar ''language'' in three inbred mouse strains. In contrast, the BTBR strain showed an unusual pattern of vocalizations, emitted only four out of ten categories of calls, and displayed more frequent, loud harmonics.  Table  S1 in the online Supporting Information).

Classification of ultrasonic vocalizations into distinct
categories. Examples of sonograms typical for each pattern category are shown in Figures 2 and S1. Audiofiles illustrating each call are available in the online Supporting Information (Sounds S1, S2, S3, S4, S5, S6, S7, S8, S9 and S10), as described at the end of the Methods section below. Figure 3 illustrates strain-dependent effects on both frequency [F (3,76) Figure 5, B6, 129X1 and FVB/NJ pups emitted a wide spectrum of call categories while BTBR pups emitted a narrower spectrum of call categories. BTBR displayed high prevalence in production of harmonic, frequency steps, composite and two-syllable calls, along with low prevalence in production of upward, chevron, short, and downward calls, indicating an apparently unusual pattern of vocalizations as compared to the three other strains tested in the present study. When analyzing each USV category separately, a straindependent effect was found on the probability of producing calls from nine of the ten call categories (number of calls in each category for each subject/total number of calls analyzed in each subject)      Figure 6A].
Righting reflex. The righting reflex, measured as latency to turn back onto all four paws when placed on the back, was tested at pnd 2 to 14 in all four strains in the first cohort used for ultrasonic vocalizations. Righting reflex latencies differed significantly between B6 and the other strains [strain6age interaction: F(18,216) = 6.88, p,.0001; Newman-Keuls p,.0001 for B6 vs all strains at pnd 2 and 4] ( Figure 6B). Righting reflex latencies were similarly measured in the second cohort of BTBR and B6, used for the full set of developmental milestones. A significant difference was detected between B6 and BTBR pups at pnd 2 and 4 [strain6age interaction: F(6,108) = 9.8, p,.0001; Newman-Keuls p,.01]( Figure 8A).     Homing test. A significant strain effect was found on latency to reach the area containing the nest litter [Mann-Whitney,

Discussion
BTBR is a minimally characterized inbred strain of mice that displays deficits in social behaviors and high levels of repetitive selfgrooming, traits relevant to the first and third diagnostic symptoms of autism [40][41][42][43][44]. We now report unusual properties of vocalizations in BTBR. BTBR pups separated from their mothers and siblings at postnatal days 2, 4, 6, 8, and 12 emitted significantly more calls, and those calls were of significantly longer duration than standard B6 mice. Total number of calls was higher in BTBR as compared to two other strains as well, 129X1, which is a close genetic relative, and FVB/NJ, which is close in body weight. In addition, BTBR calls showed smaller peak frequencies than the three other strains at pnd 6 and 8 and higher peak amplitudes than B6 at pnd 4, 6, 8 and 12, indicating overall qualitative differences as compared to the other strains. No differences were detected in body temperature after USV recording, excluding the possibility of unusual thermoregulation in BTBR.
We further investigated the specific types of calls emitted by each of the four inbred strains. Waveform patterns of the separation calls were classified into multiple categories of identified syllables emitted by BTBR, B6, 129X1 and FVB/NJ pups at postnatal age 8. Ten specific categories of calls were detected, designated as complex, harmonics, two-syllable, upward, downward, flat, chevron, short, composite and frequency steps, consistent with previous reports [21,22,24,54]. A detailed analysis of the probability of vocalizations revealed that 8-day-old B6, FVB/NJ and 129X1 pups emitted a wide repertoire of calls, including high numbers of frequency steps and complex calls. BTBR pups emitted a narrower repertoire of calls, which included high levels of harmonics and composites, but minimal numbers of chevron-shaped syllables, upward, downward, and short calls.
Mouse pup calls incorporate some properties that suggest they could serve some of the same functions as the crying of human babies, especially their ability to elicit parental retrieval behaviors [11,12,[55][56][57][58][59][60][61]. These studies indicated that most of the separation emissions carry the decisive acoustic and neurophysiological features for releasing maternal behavior, and thus are likely to be effective in communication [62][63][64]. Together with indications of left-hemisphere dominance of call perception in mice [65], these studies provide preliminary evidence for similar neuroanatomical and neurophysiological mechanisms mediating the perception of vocalizations in mice and humans.
Holy and Guo [23] first reported that vocalizations of adult male mice (C57BL/66DBA/2J F1 cross) are more complex than previously known, and share characteristics of birdsong. Mouse vocalizations consist of several categories, i.e. different syllable types, whose temporal sequencing includes the utterance of repeated phrases. Consistent with the study of song production by adult mice, Panksepp and coworkers [24] found that many USVs in adolescent B6 and BALB/cJ mice were frequencymodulated, occurring in repetitive bouts separated by periods of silence. These adolescent vocalizations were remarkably complex, with a significant effect of genotype on each distinct USV category that was classified (upward, downward, chevron, complex and punctuated calls). Further, in adult mice, call categories and frequencies varied across stages of sexual behaviors, and were influenced by amphetamine treatments and gene disruptions [8]. Inspired by these recent reports, we classified every call emitted by four strains of mouse pups into ten categories, to search for strain differences. Our 8 day old B6 mice showed similar percentages of upward-and downward-modulated calls as compared to 30 day old B6 mice [24] but higher percentages of chevrons and complex calls. Differences in age, social setting, analysis parameters, and the use of five versus ten distinct categories may explain the discrepancies in percentages reported by Panksepp and coworkers [24] and the present data. Moreover, the mean duration and peak frequencies of these types of calls are not identical in 8 day and 30 day old B6 mice, most likely due to the different ages of the subjects and differences in the test conditions. Branchi and coworkers [21] previously reported a spectrographic characterization of eight day old outbred CD-1 mice, isolated from their mothers and littermates. CD-1 pups emitted a wide spectrum of USVs classified into five categories (flat, complex, frequency steps, short and composite). As compared to the present results with B6, 129X1, and FVB/NJ, the eight day old CD-1 emitted a higher percentage of frequency steps and complex calls but low numbers of flat, short and complex calls.
Thus, vocalization categories or syllables in B6, F1 BXD, BALB/cJ, and CD-1 seem to be generally similar to each other, to the extent that comparisons can be made of results from different laboratories. BTBR appears to be an outlier strain on both higher number of calls and restricted number of categories of calls. BTBR has been proposed as a mouse model of autism, because of its low levels of sociability, juvenile play, and social transmission of food preference, and its high levels of repetitive self-grooming, resistance to change in a spatial habit in the Morris water maze, and poor shift performance in a holeboard task [40][41][42][43][44]66]. The present findings of unusual categories or syllables of vocalizations by BTBR may offer face validity to some forms of the second diagnostic symptom of autism, impaired communication. For example, some infants and young children later diagnosed with autism make atypical vocalizations. Instead of cooing and babbling, they may hum or grunt for extended periods, fail to add inflections into speech patterns, repeat ''pop up'' words out of context, squeal stereotypically, and laugh inappropriately [45,[67][68][69][70]. Others may be very irritable, cry for long periods of time, and be difficult to console [45]. More crying in these babies may be similar to the higher number of separation calls in BTBR pups, that were emitted at louder amplitudes than B6 beginning at pnd 4, and were primarily the harmonic syllable. In 81% of their 3633 calls, BTBR used only 4 of the 10 call types (harmonics, twosyllable, composite, and frequency steps). In contrast, each of the other three strains used a wider variety of the call types of their calls (2333 B6 calls, 1806 129X1 calls, 2575 FVB/NJ calls analyzed).
In rodent models of neurodevelopmental disorders, it is critical to conduct behavioral phenotyping during the early developmental period in order to document the precise onset of symptoms, identify transient signs, and provide a basis for the timing of early intervention [71,72]. BTBR reached many developmental milestones earlier than B6, an inbred strain that is commonly used in behavioral genetics and for breeding targeted gene mutations. Faster growth and development in BTBR was seen in a fully developed righting reflex response by pnd 2 as compared to pnd 6 in B6, the full negative geotaxis response acquired by pnd 4 in Figure 10. Open field activity was assayed at pnd 18. BTBR was initially more active higher than B6 on A) distance travelled and B) horizontal activity. There were no significant strain differences on C) center time and D) vertical activity in the open field arena. **p,.01 for B6 vs. BTBR. doi:10.1371/journal.pone.0003067.g010 BTBR versus pnd 8 in B6, and the more rapid growth of the body and tail length, earlier detachment of the pinnae, opening of the eyes, and incisor eruption in BTBR than B6. Faster acquisition of developmental milestones in BTBR may represent an unusual neurodevelopmental profile. Along with the larger head circumference at early ages in autism, accelerated growth of body length and higher weights were reported in three clinical studies [73][74][75]. However, it is important to note the possibility that B6 is the outlier on some measures of early neurodevelopment. The righting reflex developed most slowly in B6 as compared to the FVB/NJ and 129X1 strains which were more similar to BTBR. Body weights of B6 were similar to 129X1, while body weights of FVB/ NJ were similar to BTBR.
The homing test incorporates more complex measures of interest in social odors, cognitive abilities for differentiating own versus other home cage odors, and the sensory and motor abilities required to navigate toward litter containing the pup's own home cage odor. In the homing test, 9-day-old BTBR pups showed shorter latencies to reach the area containing nesting material from their own home cage, suggesting more rapid development of social olfactory and cognitive abilities as compared to B6 pups of same age. However, data relative to number of squares crossed during the homing test, and distance traveled during the open field test, showed that BTBR had higher general activity than agematched B6. It is possible that the shorter latencies on the homing test were related to the higher level of exploratory activity in BTBR.
From our results, it appears that qualitative as well as quantitative analyses will be useful in understanding the content of calls emitted during ultrasonic vocalization production in mice, to compare syllable categories from different strains or genotypes. It is interesting to speculate that variability in the types of syllable categories emitted by each strain represent ''dialects'' specific to each inbred strain. Similar experiments using many more inbred strains of mice will be needed, at different ages and in different environmental situations, to make definitive statements about strain variations in complex vocal repertoires. Our preliminary findings raise the notion that a large set of unusual syllable patterns in a mouse strain could be analogous to a distinct ''dialect'', ''vocabulary'' or ''language'' in that strain. However, an interpretation of ''language'' implies a communication function for the vocalizations. Detailed investigations of calls emitted by one mouse, call responses emitted by another mouse, and the resulting changes in behavior, are needed to conclude that communication has occurred. Playback experiments may be useful to discover whether mice actually communicate meaningful information to each other using ultrasonic vocalizations. If so, an accurate analysis of ultrasonic emissions could provide a reliable assay to model the second diagnostic symptom of autism, impaired communication, for use in identifying genetic and environmental causes of autism, and evaluating proposed treatments.

Behavioral observations
BTBR, B6, FVB/NJ, and 129X1 breeding pairs were purchased from The Jackson Laboratory (Bar Harbor, ME). Mice were housed in standard wire-topped Plexiglas cages (42 cm627 cm614 cm) in a temperature and humidity controlled vivarium, with food and water available ad libitum, with lights on from 19:00 to 7:00. Behavioral testing was conducted between 9.30 and 14.00 h, during the dark phase of the circadian cycle. All procedures followed NIH guidelines for the care and use of laboratory animals, and were approved by the NIMH Animal Care and Use Committee. Ten days after pairing for breeding, the females were individually housed and subsequently inspected daily at 9:30 for pregnancy and delivery. The day of birth was considered as postnatal day (pnd) 0.
Pups were tattooed on the paw with animal tattoo ink (Ketchum permanent Tattoo Inks green paste, Ketchum Manufacturing Inc., Brockville ON Canada) by loading the ink into a 30G hypodermic needle and inserting the ink subcutaneously through the needle tip into the center of the paw. The procedure was performed at two days of age, immediately after behavioral testing. The procedure causes only minor brief pain and distress and does not require the use of anesthesia.

Ultrasonic vocalizations in separated pups
Litters chosen for testing contained more than seven pups for BTBR (7.861.05), B6 (7.160.58) and FVB/NJ (7.560.34), and more than five pups for 129X1 (5.560.40; a strain known for small litters). One female and one male from each litter of BTBR, B6, FVB/NJ and 129X1 mice (n = 10 litters each strain) were used for baseline measurements of the ultrasonic vocalizations from pnd 2 to 12. Body weights and body temperatures of pups were measured after the ultrasonic vocalization test on pnd 2, 4, 6, 8 and 12. On each day of testing, each pup was placed into an empty plastic container (diameter, 5 cm; height 10 cm), located inside a sound-attenuating styrofoam box, and assessed for USVs during a five minute test. At the end of the five minute recording session, each pup was weighed and its axillary temperature measured by gentle insertion of the thermal probe in the skin pocket between upper foreleg and chest of the animal for about 30 seconds (Microprobe digital thermometer with mouse probe, Stoelting Co., Illinois, USA). No differences in patterns of calling were detected in a comparison of male and female pups, therefore data were collapsed across sex.
An Ultrasound Microphone (Avisoft UltraSoundGate condenser microphone capsule CM16, Avisoft Bioacoustics, Berlin, Germany) sensitive to frequencies of 10-180 kHz, recorded the pup vocalizations in the sound-attenuating chamber. The microphone was placed through a hole in the middle of the cover of the styrofoam sound-attenuating box, about 20 cm above the pup in its plastic container. The temperature of the room was maintained at 2261uC. Vocalizations were recorded using Avisoft Recorder software (Version 3.2). Settings included sampling rate at 250 kHz; format 16 bit. For acoustical analysis, recordings were transferred to Avisoft SASLab Pro (Version 4.40) and a fast Fourier transformation (FFT) was conducted. Spectrograms were generated with an FFT-length of 1024 points and a time window overlap of 75% (100% Frame, Hamming window). The spectrogram was produced at a frequency resolution of 488 Hz and a time resolution of 1 ms. A lower cut-off frequency of 15 kHz was used to reduce background noise outside the relevant frequency band to 0 dB. Call detection was provided by an automatic threshold-based algorithm and a hold-time mechanism (hold time: 0.01 s). An experienced user checked the accuracy of call detection, and obtained a 100% concordance between automated and observational detection. Parameters analyzed for each test day included number of calls, duration of calls, qualitative and quantitative analyses of sound frequencies measured in terms of frequency and amplitude at the maximum of the spectrum.
Waveform patterns of calls were examined in depth in twenty sonograms collected from every strain, one from each of the pups tested. The sonograms were one minute in length and selected from recordings at postnatal day 8. We classified 3633 BTBR calls, 2333 B6 calls, 1806 129X1 calls and 2575 FVB/NJ calls. Each call was identified as one of 10 distinct categories, based on internal pitch changes, lengths and shapes, using previously published categorizations [21,22,24]. Classification of USVs included ten waveform patterns described below, and illustrated visually in Figure 2 and S1 and as audiofiles (Sounds S1, S2, S3, S4, S5, S6, S7, S8, S9, S10) in Supporting Information.

1)
Complex calls displayed one syllable containing two or more directional changes in pitch, each $6.25 kHz.

2)
Harmonics displayed one main call, resembling the complex calls, but with additional calls of different frequencies surrounding the main call.

3)
Two-syllable calls consisted of two components: a main call (flat or downward) with an additional punctuated component towards the end.

4)
Upward-modulated calls exhibited a continuous increase in pitch that was $12.5 kHz, with a terminal dominant frequency at least 6.25 kHz more than the pitch at the beginning of the vocalization.

5)
Downward-modulated calls exhibited a continuous decrease in pitch that was $12.5 kHz, with a terminal dominant frequency at least 6.25 kHz less than the pitch at the beginning of the vocalization.

6)
Flat calls displayed a constant beginning and the ending of the pitch frequency remained constant (#3 kHz of each other).

7)
Chevron calls resembled an 'inverted-U', which was identified by a continuous increase in pitch $12.5 kHz followed by a decrease that was $6.25 kHz.

8)
Short calls were punctuated and shorter than 5 ms.

9)
Composite calls were formed by two harmonically independent components, emitted simultaneously. 10) Frequency steps were instantaneous frequency changes appearing as a vertically discontinuous ''step'' on a spectrogram, but with no interruption in time.
Inter-rater reliability in scoring the call categories was 98%. Call category data were subjected to two different analyses: a) straindependent effects on the frequency and duration of the vocalizations emitted by each subject at pnd 8 b) strain-dependent effects on the probability of producing calls from each of the ten categories of USV, as described below under Statistical analysis.

Developmental Milestones
To avoid potential confounds from using previously handled animals, a second cohort of B6 and BTBR pups was bred for the assays of developmental milestones, homing, and open field activity. One female and one male from each litter of BTBR and B6 mice (n = 10 litters each strain) were tested from pnd 2 to 18. Pups were transferred to a cage filled with clean bedding placed over a heating pad to maintain the temperature at 35uC. Every other day from pnd 2 to 14, pups were weighed to the nearest 0.01 g and their body and tail lengths were measured. Fur development, day of eyelid opening, pinnae detachment and incisor eruption were also recorded.
Pups were tested according to a slightly modified Fox battery [76][77][78]. The tests were conducted during the light phase of the circadian cycle, between 10:00 and 15:00 h. Each subject was tested at approximately the same time of day. Reflexes and responses were scored in the following order: 1. Righting reflex: pup turns over with all four feet on the ground when placed on its back.
2. Negative geotaxis: pup turns approximately 180u to either side when placed head down on a wire mesh screen (464 mm) held at a 45u angle. 3. Cliff aversion: pup withdraws from the edge of a flat surface when its snout and forepaws are placed over the edge of a table. 4. Forelimb grasping reflex: pup grasps the shaft of a toothpick when the forepaw is stroked. 5. Forelimb placing reflex: pup raises and places its forepaw on the surface of the edge of an object when stroked on the dorsum of the paw. 6. Vibrissa placing reflex: pup places its forepaw onto a cotton swab stroked across its vibrissae. 7. Auditory startle: pup reacts to acoustic stimuli (snapping of the fingers) by a startle response. 8. Level screen tests: pup holds onto a wire mesh screen when it is dragged across it horizontally by the tail. 9. Screen climbing test: pup climbs up the vertical wire mesh screen (90u angle) using both fore-and hind-paws. Maximal response is scored when the subject reaches the top of the vertical screen (10610 cm), which usually takes about 5 s. 10. Bar holding: pup grasps a small wire bar by its forelimbs. Bar holding was scored as present if the pup is able to hang suspended for 10 s.
Latencies were measured in seconds, using a stopwatch for righting reflex, negative geotaxis and bar holding. Other somatic and behavioral variables were rating semi-quantitatively: 0 = no response/reaction, 1 = slight response/reaction, 2 = incomplete response/reaction, and 3 = a complete adult-like response. Investigators were trained until the inter-observer reliability was greater than 95%. Unless otherwise noted, absence of a milestone was scored as zero if the mouse did not exhibit the behavior within 60 s.
Homing test (pnd 9). On pnd 9, the litter was separated from the dam and kept for 30 min in one holding cage. The cage was placed on a heating pad set at a temperature of 35uC to maintain normal body temperature of the pups in the nest. Individual pups were then transferred to a Plexiglas cage (36 cm622.5 cm, walls 10 cm high). Wood shavings from the home cage were evenly spread on one side (14 cm622.5 cm, nest area) while the rest of the cage was covered with clean bedding. The pup was placed in the middle of the Plexiglas cage and videorecorded for four minutes. The floor of the arena was virtually subdivided into squares of 7 cm67 cm each, to enhance scoring of locomotor activity from the video digital DVDs, using Noldus Observer 5.0 software (Noldus Information Technology, Leesburg, VA). Homing performance was scored for latency to reach the area containing nest litter, time spent in the area containing nesting litter and locomotor activity by square crossings.
Open field test (pnd 18). One male and one female from each litter, which were not used in the developmental milestones and homing assessments, were tested for locomotor activity on pnd 18. General exploratory locomotion in a novel environment was tested by placing pups in a VersaMax Animal Activity Monitoring System (AccuScan Instruments, Columbus, OH, USA) for a 15-min test session. To compensate for the relatively small size of 18-day-oldmale mice, the VersaMax vertical sensor was adjusted to the lowest setting of 7 cm, and the floor of the open field arena was elevated by 1.0 cm, so that the final height of the vertical sensor was 6.0 cm above the floor of the arena. The testing room was illuminated with a single 25-W red lamp and kept at a similar temperature as the colony room. The test was carried out between 10.00 and 13.00 h.
Statistical analysis. In the neonatal studies, statistical analysis based on litters as statistical units and pups as repeated trials within each litter was performed. Since no sex differences were detected, data were collapsed across sex. A mixed-model Analysis of Variance (ANOVA) with Repeated Measures was used to analyze a) neonatal USVs and developmental milestone responses, with the strain as factor and postnatal day as the repeated measures, b) homing responses, frequency and duration, c) open field activity, with session minutes as the repeated measure. Nonparametric analysis (Mann Whitney test) was used to analyze latencies in the developmental milestones and homing tests. Post-hoc comparisons were performed using Newman Keuls test only when a significant F-value was determined. For all comparisons, significance was set at P,0.05.
An ANOVA with repeated measures was performed to analyze strain-dependent effects on the vocalizations data, with the strain as factor and call categories as the repeated measures. Probability of vocalizations within strain was calculated as number of calls in each category for each subject/total number of calls analyzed in each subject and standardized by angular transformation. To better analyze strain differences within call categories, an ANOVA was performed for each category and Newman-Keuls comparison was followed by Bonferroni correction. Figure S1 Complex vocal repertoire of mouse pup separation calls. Audioclips representing examples of the ten distinct categories of calls are provided (Sounds S1, S2, S3, S4, S5, S6, S7, S8, S9 and S10). Recordings were collected and converted from mouse's ultrasonic range to the human hearing range by Avisoft software. In order to appreciate the waveforms in greater detail, a selected call of each subtype was converted from the sample rate of 250 kHz in the original wav file to 11.025 kHz, resulting in a slower speed for human listening. Found at: doi:10.1371/journal.pone.0003067.s001 (0.20 MB PPT) Sound S1 Chevron vocalization from a male B6 pup at pnd 8, corresponding to the sonogram in Figure S1  Sound S7 Two-Syllable vocalization from a male B6 pup at pnd 8, corresponding to the sonogram in Figure S1  Sound S10 Frequency step vocalization from a male B6 pup at pnd 8, corresponding to the sonogram in Figure S1