- Open Access
Female rats are not more variable than male rats: a meta-analysis of neuroscience studies
Biology of Sex Differencesvolume 7, Article number: 34 (2016)
Not including female rats or mice in neuroscience research has been justified due to the variable nature of female data caused by hormonal fluctuations associated with the female reproductive cycle. In this study, we investigated whether female rats are more variable than male rats in scientific reports of neuroscience-related traits.
PubMed and Web of Science were searched for the period from August 1, 2010, to July 31, 2014, for articles that included both male and female rats and that measured diverse aspects of brain function. Only empirical articles using both male and female gonad-intact adult rats, written in English, and including the number of subjects (or a range) were included. This resulted in 311 articles for analysis. Data were extracted from digital images from article PDFs and from manuscript tables and text. The mean and standard deviation (SD) were determined for each data point and their quotient provided a coefficient of variation (CV) as a measure of trait-specific variability for each sex. Additionally, the results were coded for the type of research being measured (behavior, electrophysiology, histology, neurochemistry, and non-brain measures) and for the strain of rat. Over 6000 data points were extracted for both males and females. Subsets of the data were coded for whether male and female mean values differed significantly and whether animals were grouped or individually housed.
Across all traits, there were no sex differences in trait variability, as indicated by the CV, and there were no sex differences in any of the four neuroscience categories, even in instances in which mean values for males and females were significantly different. Female rats were not more variable at any stage of the estrous cycle than male rats. There were no sex differences in the effect of housing conditions on CV. On one of four measures of non-brain function, females were more variable than males.
We conclude that even when female rats are used in neuroscience experiments without regard to the estrous cycle stage, their data are not more variable than those of male rats. This is true for behavioral, electrophysiological, neurochemical, and histological measures. Thus, when designing neuroscience experiments to include both male and female rats, power analyses based on variance in male measures are sufficient to yield accurate numbers for females as well, even when the estrous cycle is not taken into consideration.
Twenty years ago, the NIH began requiring all clinical research to include both women and men in research and to report research findings for both sexes. While women are now included in research, fewer women than men are used as subjects, and findings are still not being reported by sex . As a result, any chances for achieving personalized medicine for women in the near future seem remote, as the scientific basis upon which medical decisions are being based are still on data primarily derived from men.
Accounting for sex as a biological variable in all biomedical research is considered fundamental for enhancing rigor and reproducibility in preclinical research [2, 3]. Yet there is considerable concern among preclinical scientists that including female animals will increase costs and variability in data collected . This bias remains entrenched in spite of evidence demonstrating that there are important fundamental biological differences between the sexes, and failure to elucidate these differences is impeding progress in both basic and clinical research [2, 3]. As summarized by Clayton [3, p. 522], “A continual growth in knowledge of the influence of sex at molecular, cellular, and biochemical levels and the various ways that sex exerts influence will inform the design and conduct of additional biomedical research, which is imperative to the NIH mission of turning discovery into health. Understanding scientific findings in the context of sex—be they similarities, differences, and/or complex nuances—is crucial for correctly applying research-derived knowledge toward achieving our ultimate objectives”.
Nevertheless, there is substantial bias in biomedical research to not study female rats or mice and/or to not report the sex of the subjects at all [5–7]. Not including female rats or mice in neuroscience research has been justified due to the variable nature of female data caused by hormonal fluctuations associated with the female’s reproductive cycle, in spite of lack of data in support of this position. A recent meta-analysis reported that female mice are not inherently more variable than male mice across diverse physiological traits . Similar results have been obtained for measures of gene expression in mice and humans .
In this study, we investigated whether female and male rats differ in their variability in studies that focused on neuroscience outcomes. We chose to focus on one field for this study in order to examine a dataset that is relatively homogenous, so that failure to find a sex difference in variability would not be due to heterogeneity of the measures being examined. We examined studies that included intact adult male and female rats. The majority of the studies used female rats without regard to the stage of the estrous cycle, but we also examined 26 studies that included male and female rats at specific stages of the estrous cycle. We now report that female rats are not more variable than male rats on studies of neuroscience-related traits. This is true when females are used without regard to the estrous cycle or when studied at specific days of the estrous cycle.
PubMed and Web of Science were searched for the period from August 1, 2010 to July 31, 2014. The PubMed search terms used were as follows: (1) (rat AND gender differences) AND (brain OR neuroscience OR neuron) = 411 articles and (2) (rat AND sex differences) AND (brain OR neuroscience OR neuron) = 525 articles. When these lists were manually combined, this yielded 543 unique articles. On Web of Science, the search terms were TS = (male and female) AND TS = (neuro* AND rat) NOT TS = (adolescent) NOT TS = (mice). These articles were then filtered by neuroscience, behavior, article (not review) and 2010–2014. The Web of Science search generated 743 references; these were manually curated to identify 151 unique additional relevant references using the titles and abstracts (manually eliminated January 1, 2010–July 31, 2010, and any in August 2014). When combined with the PubMed search there were a total of 562 articles. These articles were manually reviewed to determine appropriateness for inclusion. Only empirical articles using both male and female gonad-intact adult rats, written in English, and describing the number of subjects (or a range) were included—resulting in 311 articles for analysis. A list of the articles used is included in the supplemental information for this article (see Additional file 1).
Data were extracted from digital image files generated from high-resolution screenshots of article PDFs and from manuscript tables and text. Vector graphics software (Adobe Illustrator) was used to quantify the mean and standard deviation (STDEV) or standard error of the mean (SEM) values directly from figure images (in mm), which provided a relative measure of the mean and STDEV/SEM for each data point as described in . Briefly, figures were imported into Adobe Illustrator, and for each data point used, rectangles were positioned on the graphs over the SEM/SD bar from the middle of the data point or bar to the end of the error bar. A rectangle was also positioned from the X-axis to the middle of the data point or bar (with corrections if the scale was discontinuous), and the length of each of these rectangles in millimeters (determined by the graphics software) was used as a relative measure of the mean and error reported. Data were only used if the mean and STDEV or SEM could be extracted from the article. Data presented in tables were transcribed directly from the table. For line graphs with more than three time points, values were obtained from the beginning, middle, and end of the time course, so that no one study contributed a disproportionate number of data points to the overall analysis. When a range for the number of subjects was given, the lowest number in the range was used. Data were collected by 10 undergraduate students with an inter-rater reliability coefficient >0.96.
Results were coded for the type of research (behavior, electrophysiology, histology, neurochemistry, and non-brain measures). Behavior was any behavioral measure (N = 2245 data points). Electrophysiology included measures of electrical neural activity (LTP, unit activity, cell clamp electrophysiology, etc.; N = 364 data points). Neurochemistry was any measure of neurotransmitter or neurotransmitter receptor amount, protein amount, synthesis, second messengers, or neurotransmitter release (N = 1809 data points). Most of the molecular studies were included this category. Histology was measures of cellular location, dendritic/axonal branching, brain regions, and brain region activity, including c-fos; measures that quantify physical structure in the brain (N = 1233 data points). Non-brain measures (N = 601 data points) were any measures of non-central nervous system biology including body weight (N = 127 data points), blood/serum hormone measures (N = 214 data points), cardio measures (heart rate, blood pressure, etc.; N = 54 data points), and blood or organ measurement of exogenous compounds or organ weights (“organ” N = 207 data points).
For histology and neurochemistry measures, each pair of data points was also coded for whether male and female values were significantly different from each other. For the histology data, the number of data points each for males and females was as follows: no sex difference = 648 data points; sex difference = 585 data points. For the neurochemistry data, the number of data points each for males and females was as follows: no sex difference = 1177 data points; sex difference = 451 data points; not measured = 181 data points.
In a subset of manuscripts, one or more estrous cycle stages were recorded (n = 26 manuscripts). Analysis was without respect to subject category. Not all studies examined all phases of the estrous cycle. We obtained the following number of values: males = 343 data points; diestrus = 330 data points; proestrus = 151 data points; estrus = 241 data points.
For neurochemistry and behavior measures (n = 4137 data points, in total), we also evaluated whether the animals were housed individually (N = 872 data points; 29 studies), in pairs or two to three/cage (N = 1311 data points; 57 studies), three or more per cage (N = 1062 data points; 47 studies), or not reported (N = 892 data points; 39 studies or 22.6 % of the studies). Housing conditions were the same for males and females in all studies. Thus, the number of data points is the same for both males and females.
The strain of rat was coded when it was indicated in the article (Sprague-Dawley: N = 2871 data points; Long-Evans: N = 1053 data points; Wistar: N = 2221 data points; Norway Brown: N = 50 data points).
The coefficient of variation (CV) was calculated as the standard deviation divided by the mean (STDEV/mean) for each data point. Male-female differences were analyzed by paired t tests (pairing by data points for male and female collected in an individual study) or analysis of variance (ANOVA; depending on whether individual traits or multiple traits were being compared, respectively). The ANOVAs were followed by pairwise comparisons with Tukey’s multiple comparisons test.
Female to male ratios of CV were calculated to determine if the distribution of variation differed by sex. To calculate, the female to male ratio = [CV female/(CV female + CV male)]. The theoretical mean for the ratios would be 0.5 if males and females did not differ in the coefficient of variability. The CV ratios for each trait were tested for each sex against the theoretical mean by t test to examine whether each differed from 0.5.
Inter-rater reliability was determined by Pearson r correlation to be 0.960–0.997.
Female and male trait variability
There were no sex differences in the coefficient of trait variability (CV = STDEV/mean) for any of the neuroscience measures when the CVs for data points obtained from males and females for a given measure from each study were compared with paired t tests (Table 1). For behavior, electrophysiology, histology, and neurochemistry data, we found that females were not more variable than males (Fig. 1).
There were, however, differences among the traits in the extent of variability. On a two-way ANOVA (sex X trait), there was no main effect of sex (F(1, 12,500) = 1.927; p = 0.1651) and no significant sex by trait interaction (F(4, 12,500) = 1.574; p = 0.1787). There was a main effect of trait (F(4, 12,500) = 18.98; p < 0.0001) indicating that the CVs for some traits were more variable than other traits. Using Tukey’s multiple comparisons test, the CV for behavior for males was greater than that of histology or neurochemistry (p < 0.001), and the CV for histology was lower than that for neurochemistry or non-brain measures (p < 0.01). For females, the CV for behavior was also was greater than that of histology or neurochemistry (p < 0.001), the CV for histology was lower than the CV for neurochemistry (p < 0.01), and the CV for non-brain measures was greater than that for electrophysiology (p < 0.05), histology (p < 0.001), and neurochemistry (p < 0.01). This indicates that even though males and females do not differ from each other, behavioral measures were more variable for both males and females than were neurochemistry and histology measures. On the other hand, histology CV data were less variable for both males and females than neurochemistry or the non-brain measures.
For “non-brain measures,” there was a significant difference when females and males were compared with a paired t test (Table 1). The non-brain measures included measures where the mean would be expected to vary with the estrous cycle (body weight, heart rate, blood pressure, organ weights, serum gonadal and adrenal hormones, etc.). To further investigate the source of the variance, we further assigned these measures to sub-categories. These categories were as follows: (1) body weight—body weight/fat weight (N = 127); (2) endo—hormone measures (N = 214); (3) cardio—blood pressure, heart rate, and other cardiac measures (N = 54); and (4) blood/organ—measures of organ weight, organ or blood proteins or exogenous substances, and other organ-specific measures (N = 207). As illustrated in Fig. 2, the measures of blood/organ were the primary source of the sex difference in the non-brain measures (t = 1.952; DF = 412; p = 0.0516; Mann-Whitney U test p = 0.036).
Distribution of CV ratios
There was a trend for the distribution of CV ratios (female CV/(female CV + male CV)) to vary by trait on ANOVA (F = 2.594, DF = 3, 5650; p = 0.0509). As the variance was not normally distributed, the Brown-Forsythe test was considered appropriate to apply and there the analysis indicated that the distribution of CV ratios varied by trait (F = 11.91, DF = 3, 5650, p < 0.0001; Fig. 3).
We then went on to examine whether there were sex differences in the CV ratios for the different traits. The theoretical mean for the ratios would be 0.5 if males and females did not differ in the CV ratio. When the CV ratios for each trait are tested for each sex against this theoretical mean by t test, there was no sex difference in the CV ratio on the behavior (mean = 0.4943 ± 0.0057; t = 1.893, DF = 2243) or histology (mean = 0.5050 ± 0.005; t = 1.130, DF = 1232) trait categories, and males were more variable than females on the electrophysiology (mean = 0.4863 ± 0.014; t = 2.092, DF = 363, p = 0.037) and the neurochemistry (mean = 0.4916 ± 0.0084; t = 2.336, DF = 1824, p = 0.0196) trait categories. Females were more variable than males on the non-brain measures (mean = 0.5308 ± 0.0308, t = 4.316, DF = 600, p < 0.0001).
CV values when there is a sex difference in the value
We went on to examine whether there were sex differences in CV values if the data points being compared differed significantly between males and females. This analysis examined two trait categories: the histology measures (where the CV ratio distribution did not differ between females and males) and the neurochemistry measures (where CV ratios indicated greater variability in males). As illustrated in Fig. 4, there was no effect of whether a given data point was significantly different between the sexes on the CV values for either histology or neurochemistry. However, in the neurochemistry category, CVs were greater in females when the mean did not differ significantly from those of males as compared to females whose means differed from those of males (Sidak’s multiple comparisons test, p < 0.05).
Impact of estrous cycle on trait variability
There was no significant effect of sex/estrous cycle stage on CV with a one-way ANOVA (F(3, 1061) = 2.199, p = 0.0865; Fig. 5). Females did not differ from males on any day of the estrous cycle nor did any of the female groups differ from each other.
Impact of housing on trait variability
For the neurochemistry and behavior values, we also examined whether the housing conditions contributed to the variability in trait data. As can be seen in Fig. 6, there was no sex difference in the effect of housing conditions on trait (F(1, 8266) = 0.4282, p = 0.5139). Overall, there was a main effect of housing conditions on trait (F(3, 8266) = 6.175; p = 0.0003), but no sex by trait interaction (F(3, 8266) = 0.4282; p = 0.5129). These effects do not change if the data from studies where housing conditions were not reported re excluded.
Impact of rat strain on trait variability
Lastly, we examined whether the strain of rat contributed to variability in data and whether there were effects of sex on the CV; however, there was no effect of strain on sex differences in CV (F(1, 12,382) = 0.0889, p = 0.765). Overall, male Sprague-Dawley rats were more variable than male Wistar rats (two-way ANOVA; main effect of strain: F(3, 12,382) = 3.941; p = 0.008; subsequent Tukey’s multiple comparisons test, p < 0.05) (see Fig. 7). There were no other effects of strain.
These results indicate that among diverse traits relevant to neuroscience, female rats are no more variable than male rats. When the data are categorized by type of information reported, some types of data have greater intrinsic variability than others: behavioral data are more variable than histology or neurochemistry data, for example; but females and males did not differ in this regard. Thus some types of neuroscience tests may yield more precise, or less variable, data values, but this does not differ by sex. An important and novel aspect of this analysis is that, there were no sex differences evident when males were compared with (1) either randomly cycling females or (2) females at specific, defined stages of the estrous cycle. Moreover, females did not exhibit greater variability at any stage of the estrous cycle, compared with males or with females at other estrous cycle stages.
It is important to note that trait variability was not greater for females or males even when there was a significant sex difference in the mean value reported in the studies analyzed. A significant difference between the sexes on a given measure does not mean that females are more variable than males. What our findings mean is that it is possible to see sex differences in neuroscience studies when equal numbers of male and female rats are used.
There was greater variability among females in the “non-brain” category. Upon further analysis, three of the four defined sub-categories of “non-brain” exhibited no difference whatsoever between males and females. For one indistinct sub-category with a relatively small sample size, there was greater variability in females. Thus, there will be instances where females are more variable than males.
Recently, Itoh and Arnold  conducted a meta-analysis of 103 human microarray datasets and 190 mouse microarray datasets to examine gene expression variability in males and females. The results indicated that variability was similar for females and males in humans and in mice and no evidence that female gene expression was more variable than male gene expression in either species. The present report extends the study of sex differences in variability to a species that is widely used in neuroscience and documents the overall absence of sex differences in variability across diverse traits of interest to neuroscientists.
In conclusion, female rats are not more variable than male rats in neuroscience research. Across a substantially large sampling of research, the data indicate that on average, females exhibit the same (or less) variability on a given trait that male rats do. One implication of these data is that for those investigators initiating research on female rats, power calculations based on data from males would likely be sufficient to determine the number of female subjects needed in order to see a sex difference. There will be particular topics where well-documented effects of the estrous cycle should be considered by investigators in the experimental design in order to get meaningful results. In all datasets, there exists a distribution of CV ratios; thus one single trait may be more variable in males than females (or vice versa). On the other hand, for topics where females have not been studied, these data suggest that inclusion of intact females, without regard to estrous cycle, and intact males is a valid approach to learn about females in neuroscience research.
Analysis of variance
Coefficient of variation
Standard error of the mean
Mazure CM, Jones DP. Twenty years and still counting: including women as participants and studying sex and gender in biomedical research. BMC Womens Health. 2015;15:94.
Clayton JA, Collins FS. Policy: NIH to balance sex in cell and animal studies. Nature. 2014;509:282–3.
Clayton JA. Studying both sexes: a guiding principle for biomedicine. FASEB J. 2016;30:519–24.
Fields RD. NIH policy: mandate goes too far. Nature. 2014;510:340.
Beery AK, Zucker I. Sex bias in neuroscience and biomedical research. Neurosci Biobehav Rev. 2011;35:565–72.
Flórez-Vargas O, Brass A, Karystianis G, Bramhall M, Stevens R, Cruickshank S, Nenadic G: Bias in the reporting of sex and age in biomedical research on mouse models. Elife. 2016 Mar 3;5. pii: e13615. doi: 10.7554/eLife.13615.
Yoon DY, Mansukhani NA, Stubbs VC, Helenowski IB, Woodruff TK, Kibbe MR. Sex bias exists in basic science and translational surgical research. Surgery. 2014;156:508–16.
Prendergast BJ, Onishi KG, Zucker I. Female mice liberated for inclusion in neuroscience and biomedical research. Neurosci Biobehav Rev. 2014;40:1–5.
Itoh Y, Arnold AP. Are females more variable than males in gene expression? Meta-analysis of microarray datasets. Biol Sex Differ. 2015;6:18.
We gratefully acknowledge the support from the National Institute on Drug Abuse R01 DA039952 to JBB. The funding agency was not involved in the design, analysis, or interpretation of the data reported. The authors would like to thank the undergraduates who helped to collect the data: Rachel Moore, Sam Gieseker, Qisi Yao, Cosette Kathawa, Jennifer Veith, Will Zech, John Kruszewski, Nikki Koll, and Krisitn Soreide.
Availability of data and materials
The articles used for collection of the data are listed in the supplementary material for this article (see Additional file 1). Spreadsheets are available at the following URL (https://umich.box.com/s/en23t90uem280fuf8kwk9nnsnd0kcj06 ).
JBB supervised the undergraduates collecting the data from the articles, collected the data, collated and analyzed the data, and wrote the first draft of the manuscript. BJP helped with the methods for collection of the data, collected the data, and provided the feedback and editing of the manuscript. JWL collected the data, collated the data, and provided the feedback and editing of the manuscript. All authors read and approved the final manuscript.
The authors declare they have no competing interests.
Pubmed references used. (DOCX 240 kb)