Svoboda | Graniru | BBC Russia | Golosameriki | Facebook
Skip to main content
Advertisement
  • Loading metrics

Paternally expressed gene 3 (Pw1/Peg3) promotes sexual dimorphism in metabolism and behavior

  • Karo Tanaka,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing

    Affiliation Stem Cells and Regenerative Medicine, Institute of Cardiometabolism and Nutrition (ICAN), INSERM U1166, University of Pierre and Marie Curie Paris VI, Paris, France

  • Vanessa Besson,

    Roles Investigation, Methodology, Resources

    Affiliation Stem Cells and Regenerative Medicine, Institute of Cardiometabolism and Nutrition (ICAN), INSERM U1166, University of Pierre and Marie Curie Paris VI, Paris, France

  • Manon Rivagorda,

    Roles Data curation, Formal analysis, Investigation, Methodology, Resources

    Current address: Integrative neurobiology, INSERM U1151, Institut Necker Enfants Malades, Paris, France

    Affiliation Hormonal Regulation of Brain Development and Functions, INSERM U1151, Institut Necker Enfants Malades, Paris, France

  • Franck Oury,

    Roles Conceptualization, Data curation, Formal analysis, Resources, Supervision, Writing – original draft, Writing – review & editing

    Current address: Integrative neurobiology, INSERM U1151, Institut Necker Enfants Malades, Paris, France

    Affiliation Hormonal Regulation of Brain Development and Functions, INSERM U1151, Institut Necker Enfants Malades, Paris, France

  • Giovanna Marazzi,

    Roles Conceptualization, Data curation, Formal analysis, Supervision, Validation, Writing – original draft, Writing – review & editing

    Current address: Paris Cardiovascular Research Center, INSERM U970, European Hospital Georges Pompidou, Paris, France

    Affiliation Stem Cells and Regenerative Medicine, Institute of Cardiometabolism and Nutrition (ICAN), INSERM U1166, University of Pierre and Marie Curie Paris VI, Paris, France

  • David A. Sassoon

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Writing – original draft, Writing – review & editing

    [email protected]

    Current address: Paris Cardiovascular Research Center, INSERM U970, European Hospital Georges Pompidou, Paris, France

    Affiliation Stem Cells and Regenerative Medicine, Institute of Cardiometabolism and Nutrition (ICAN), INSERM U1166, University of Pierre and Marie Curie Paris VI, Paris, France

Abstract

The paternally expressed gene 3 (Pw1/Peg3) is a mammalian-specific parentally imprinted gene expressed in stem/progenitor cells of the brain and endocrine tissues. Here, we compared phenotypic characteristics in Pw1/Peg3 deficient male and female mice. Our findings indicate that Pw1/Peg3 is a key player for the determination of sexual dimorphism in metabolism and behavior. Mice carrying a paternally inherited Pw1/Peg3 mutant allele manifested postnatal deficits in GH/IGF dependent growth before weaning, sex steroid dependent masculinization during puberty, and insulin dependent fat accumulation in adulthood. As a result, Pw1/Peg3 deficient mice develop a sex-dependent global shift of body metabolism towards accelerated adiposity, diabetic-like insulin resistance, and fatty liver. Furthermore, Pw1/Peg3 deficient males displayed reduced social dominance and competitiveness concomitant with alterations in the vasopressinergic architecture in the brain. This study demonstrates that Pw1/Peg3 provides an epigenetic context that promotes male-specific characteristics through sex steroid pathways during postnatal development.

Author summary

Pw1/Peg3 is under parental specific epigenetic regulation. We propose that Pw1/Peg3 confers a selective advantage in mammals by regulating sexual dimorphism. To address this question, we examined the consequences of Pw1/Peg3 loss of function in mice in an age- and sex-dependent context and found that Pw1/Peg3 mutants display reduced sexual dimorphism in growth, metabolism and behaviors. Our findings support the intralocus sexual conflict model of genomic imprinting where it contributes in sexual differentiation. Furthermore, our observations provide a unifying role of sex steroid signaling as a common property of Pw1/Peg3 expressing stem/progenitor cells and differentiated endocrine cells, both of which remain proliferative in response to gonadal hormones in adult life.

Introduction

Parental genomic imprinting is a form of epigenetic regulation by which one allele of a gene is expressed according to its parent-of-origin. In vertebrates, this form of imprinting is unique to placental mammals and its evolutionary advantage is still under active debate [13]. The parental conflict (or kinship) [4] and maternal-offspring coadaptation theories [5] are two widely recognized concepts to explain why parental genomic imprinting arose in mammals. Independently, Day and Bonduriansky proposed an ‘intralocus sexual conflict’ model [6] that predicts a physiological role for genomic imprinting in the genetic architecture of sexually dimorphic traits. This hypothesis is applicable to any species and traits under sex-specific selection pressure. However, empirical exploration of the role of imprinted genes in sexual differentiation is relatively limited [7,8].

Human diseases associated with deregulated genomic imprinting and gene knockout studies in mice have established pivotal roles of genomic imprinting in growth, metabolism, reproduction, and behavior [911]. In human and mouse, many imprinted genes are clustered in distinct chromosomal regions and are typically co-expressed in organs and tissues that regulate homeostasis of the whole-body energy metabolism, such as the brain hypothalamic region, liver, pancreas, muscle, fat, gonads, and placenta [9,12,13]. The generation of mutants corresponding to several imprinted genes in mice demonstrated global metabolic changes and their imprinting status (i.e. maternal or paternal) often correlates with inverse metabolic outcomes. Specifically, paternally expressed genes such as Magel2, Dlk1 and Zac1 promote growth and energy expenditure and restrict adiposity whereas the maternally expressed genes, H19 and Grb10, suppress growth and increase adiposity ([10] and references therein). Genome-wide transcriptome analyses have further demonstrated that inactivation or overexpression of a single imprinted gene alters the expression profile of multiple imprinted genes, suggesting that imprinted genes act in networks to coordinate cellular and organ development and functions [14].

The Pw1/paternally expressed gene 3 (hereafter referred to as Pw1) is a mammalian-specific, parentally imprinted gene that is widely expressed during early embryonic development and becomes restricted to subset of tissues in adulthood [15,16]. Using a Pw1 reporter transgenic mouse line (Pw1IRESnLacZ), we showed that Pw1 is expressed in a wide array of adult stem/progenitor cells [17]. Studies of different types of progenitor cells, all of which express high levels of Pw1, demonstrated that Pw1 dysfunction alters stem cell competence, self-renewal capacity, and cell cycle behaviors [1821]. At a molecular level, Pw1 modulates cell stress pathways including TNFα-NFκB signaling in cell growth and survival [22], p53 signaling in apoptosis [23,24], and decolin-induced autophagy [25]. The PW1 protein also acts as a transcription factor that is shown to regulate expression of mitochondrial genes in the brain [26] as well as oxytocin receptor [27]. To date, several lines of Pw1 mutant mice have been generated by different gene targeting strategies [2830]. Mice carrying a paternally inherited mutant allele for Pw1 consistently displayed reduced pre- and postnatal growth in all models. Pw1+/p- adult males were also shown to have altered energy homeostasis such as increased body fat and reduced thermogenesis, whereas metabolic phenotypes of female counterparts were not fully characterized in detail [31]. By contrast, a delayed onset of oestrus cycle and alterations in the reproductive physiology, such as smaller litter size and mature oocytes, were demonstrated [30,31]. It has been also reported that Pw1 deficient mice display deficits in adaptive traits, such as maternal care in females [28], and sexual experience-dependent olfactory learning in males [32]. All these findings indicate a significant involvement of Pw1 in sex-hormone dependent physiology, but the underlying mechanism by which this paternally expressed gene exerts such diverse biological functions remained unresolved.

In this study, we characterized paternally inherited Pw1 deficient phenotypes in male and female mice at different stages of postnatal development. We identified specific growth factor and hormonal axes that are deregulated at critical stages of postnatal development. At the cellular level, we demonstrate co-localization of Pw1 in sex-hormone dependent cell types in various organs. Our results point to a central role for Pw1 in establishing sexual dimorphism in mammals that regulates overall sex-specific physical traits, metabolism and behavior.

Results

Reduced masculinization of growth and metabolism in Pw1 deficient males

Mice carrying a paternally inherited mutant allele (Pw1+/pat) were distinguishable from their wildtype (WT) littermates (Pw1+/+) by a smaller size at birth and a reduced postnatal weight gain, as previously reported by our group and others [19,28,29]. Comparison of male and female littermates in the postnatal growth phase revealed that body weight was identical between males and females up to 3.5 weeks of age in both genotypes and Pw1+/p- mice were significantly smaller (Fig 1A). Sex differences emerge at 4 weeks of age in both genotypes, with slight delay in the Pw1+/pat littermates. Body weight at 7 weeks of age revealed a significant interaction between sex and genotype (p<0.05), and Pw1+/pat- are significantly smaller in both sexes (p<0.0001) while females are significantly smaller than males (p<0.0001). Multiple comparisons revealed all four groups are different (p<0.0001), however notably, there were no differences detected between Pw1+/+ males and Pw1+/pat- females. We observed a positive correlation between random-fed blood glucose and body weight during the postnatal growth phase regardless of sex and genotype (Fig 1B, left), and the Pw1+/pat- mice displayed reduced glucose levels up to 2 months of age (Fig 1B, right, p<0.001). Concomitantly, food consumption between 2 months and 2.5 months of age was reduced in a sex dependent manner (Fig 1C). We noted that for all sexually dimorphic parameters examined, Pw1+/pat- males were similar to Pw1+/+ female littermates, which indicated a role for Pw1 in the control of male-specific sexual differentiation during postnatal development.

thumbnail
Fig 1. Reduced masculinization of metabolisms in Pw1+/pat- males.

(A) Postnatal growth of Pw1+/pat- compared with Pw1+/+ littermates showing Pw1+/pat- animals are smaller in both sexes throughout the growing phase (p<0.001), and the onset of sexual dimorphism in body weight is delayed in Pw1+/pat- at four weeks of age. Two-way ANOVA test revealed significant interaction between genotype and sex after 6 weeks of age. (B) Random-fed glucose levels at 2 month-old and its positive correlation with body weight (r = 0.690, p<0.0001). Each symbol represents independent measurement. (C) Sex dimorphic food consumption at 2 month-old. Two way ANOVA test with multiple comparisons demonstrated significant interaction between sex and genotype. (D) Sexual dimorphisms in body composition in young adults (n = 8–14 each group). Male-specific increase in lean mass and decrease in fat mass in Pw1+/+ males at 2–3 month of age was less prominent in Pw1+/pat- males, while the females Pw1+/pat- are proportionally smaller in lean and fat mass. * in black represents comparison in males and * in red represents comparison in females. (E) An inverse correlation between lean mass and fat mass at 10 month of age was only found in the Pw1+/+ males (r = 0.757, p<0.001) and not in the Pw1+/p- males (r = 0.076, p = 0.771) nor in the females (r = 0.006, p = 0.981). Comparison between genotypes and sexes was performed using two-way ANOVA with Tukey’s multiple comparisons. Correlation was determined with simple linear regression analysis. *P<0.05, **P<0.01, and ***P<0.001. NS: non-significant. Values are mean ± SEM. Each symbol represents individual animals in (B) and (E).

https://doi.org/10.1371/journal.pgen.1010003.g001

In adult mammals, including humans and mice, males are typically larger with an increased skeletal mass, whereas females are smaller with higher adiposity. To further monitor the sex-dependent postnatal development and maturation, we performed a longitudinal analysis of body composition (lean/fat mass) of Pw1+/pat- males and females in comparison to Pw1+/+ littermates using non-invasive NMR imaging. During secondary sexual maturation, both male and female Pw1+/pat- animals showed reduced lean mass at all time points analyzed, but the difference became more marked in males (by 20%, p<0.001) than in females (by 15%, p<0.001) (Fig 1D, top left). In contrast, the fat mass development was highly sex-dependent. Pw1+/+ males manifested transient reduction of fat mass at 3 months of age whereas Pw1+/pat- males did not undergo this transition resulting in an accelerated fat accumulation in later adulthood (Fig 1D, bottom left). When expressed in percentage, the composition of lean mass shows a steady increase up to three months of age in Pw1+/+ males (Fig 1D, top right), while fat mass decrease comparatively (Fig 1D, bottom right). These male-specific changes in body composition were less prominent in Pw1+/pat- males. Furthermore, the phenotype is highly specific to male, as lean and fat mass were proportionally reduced in the Pw1+/pat- females as compared to their Pw1+/+ littermates. Therefore, there was no difference in % of body composition in females (Fig 1D, right). Two-way ANOVA test revealed an interaction between genotype and sex with age, indicating that loss of Pw1 has significantly different impacts on body composition in males and females. There was an inverse correlation between lean mass and fat mass in mature age specifically in Pw1+/+ adult males (r = 0.757, p<0.001) (Fig 1E). Taken together, Pw1+/pat- mice displayed a significant reduction in male-specific body growth.

Pw1 deficient mice have altered GH/IGF signaling that reduces body size and sexual dimorphism during postnatal development

The growth hormone (GH)/insulin-like growth factor-1 (IGF-1) axis plays a pivotal role in directing postnatal growth and regulates fat metabolism [33], whereas gonadal androgens stimulate the male-specific pulsatile secretion of GH in early puberty [34,35] thereby promoting sexually dimorphic patterns of somatic growth and body composition. The anabolic effect of GH is exerted by the stimulation of endocrine IGF-1 production primarily in the liver, and the circulating IGF-1 levels are considered as an indicator of GH activity in the postnatal growth phase [33]. Therefore, we examined IGF-1 activity in the Pw1 mutant mice during postnatal development. Plasma IGF-1 levels correlated with body weight at 3 weeks old, as commonly expected (Fig 2A), and IGF-1 levels were reduced in the Pw1+/pat- mice as compared to Pw1+/+ littermates (Fig 2B). We further monitored the circulating levels of IGF-1 in the same animals weekly up to 6 weeks of age corresponding to the period when the circulating IGF-1 levels dynamically change in a sex-dependent manner [36]. At five weeks of age, IGF-1 levels were significantly different between sexes (p<0.001) when IGF-1 levels decline in females corresponding to an earlier cessation of growth and increase in males to further promote their growth. Therefore, the levels of IGF-1 in Pw1+/+ males were significantly higher as compared to Pw1+/+ and Pw1+/pat- females (p<0.05 and p<0.01, respectively). Remarkably, no statistical differences were found between Pw1+/pat- males and females of both genotypes. We conclude that the Pw1+/pat- animals display reduced sex-specific regulation in IGF-1 secretion compared to the wild-type littermates.

thumbnail
Fig 2. Deregulated GH/IGF axis in Pw1+/pat- youngs and insulin homeostasis in Pw1+/pat- adult males.

(A) Circulating IGF-1 levels at 3 weeks of age and its positive correlation with body weight. Correlation was determined with simple linear regression analysis. (B) Blunted sexual dimorphism of circulating IGF-1 dynamics in Pw1+/pat- mice as compared to that of Pw1+/+. N = 4–8 each group from 5 litters. Two-way ANOVA showed that the IGF-1 levels at 3 weeks old are significantly lower in Pw1+/pat- mice, with no difference between sexes. **P<0.01 in black represent comparison between genotypes, whereas ***P<0.001 in red represent comparison between sex. (C) mRNA expression of growth hormone (Gh) in the pituitary gland, and GH receptor (Ghr) and insulin-like growth factor (Igf1) in the liver at 3 weeks of age. Values were normalized with Tbp and presented relative to the Pw1+/+ male littermates (n = 6–12 each group). (D) Random-fed blood insulin levels at 3 months and 6 months of age in Pw1+/pat- and Pw1+/+ littermates (n = 4–6). (E) Representative insulin tolerance test (ITT) on Pw1+/pat- and Pw1+/+ males from a single litter (n = 3 each genotype) at 3 months and 6 months of age. Similar results were obtained from two other litters. (F) Representative oral glucose tolerance test (OGTT) with insulin secretion measurement on the same set of mice as in ITT. (G) Liver size and mRNA expression of lipogenic genes in the liver at 3 months and 6 months of age showing an age-dependent development of hepatic steatosis in Pw1+/pat- mice as compared to Pw1+/+ littermates. Srebp1, sterol regulatory element binding protein 1; Acaca, acetyl-CoA carboxylase alpha; Fasn, fatty acid synthese; Scd1, Stearoyl-CoA desaturase; Pparg1 & 2, peroxisome proliferator activated receptor gamma 1 & 2. (H) Fat deposition revealed by Oil Red-O staining in 8-month-old livers (n = 4–6 for each group). Lipid droplets were quantified in number and in size using particle analysis tool in Image-J software. Values are mean ± SEM; *, P < 0.05; **, P < 0.01; ***, P < 0.001 by two-way ANOVA with Tukey’s multiple comparisons.

https://doi.org/10.1371/journal.pgen.1010003.g002

Based on the observation that circulating IGF-1 levels are reduced in Pw1+/pat- mice at 3 weeks of age, we performed gene expression analysis on pituitary gland and liver of different sets of littermates. Consistently, the expression of Gh in the pituitary gland was significantly reduced in Pw1+/pat- males as compared to Pw1+/+ littermates, whereas growth hormone receptor (Ghr) and Igf1 expression levels in the liver also showed a trend of down-regulation (Fig 2C). Taken together, these results show a global suppression of GH/IGF-1 activity during postnatal growth in Pw1+/pat- mice in a sex-dependent manner.

Pw1 deficiency deregulates insulin sensitivity and increases adiposity in adult males

Insulin is a key regulator of energy and fat metabolism throughout life. Its anabolic action promotes postnatal growth after weaning [37], however, chronically elevated insulin levels are associated with obesity and abnormal fat metabolism [38]. To evaluate the steady state insulin levels, we measured blood insulin levels in fed animals. At 3 months of age, circulating insulin levels were lower in females than males (p<0.001) corresponding to their lower levels of glycemia (Fig 2D). The insulin levels of Pw1+/pat- were also reduced at this age although this trend is only confirmed with Fisher’s LSD test. By contrast, the Pw1+/pat- genotype exhibited a male specific increase in insulin levels at 6 months of age.

Linear regression between insulin levels and body composition revealed a positive correlation between plasma insulin levels and lean mass at 3 months of age in Pw1+/+ and Pw1+/pat- males (r = 0.714, p<0.01, and r = 0.699, p<0.01), respectively (S1A Fig). At 6 months of age, on the other hand, the insulin levels correlated better with fat mass in Pw1+/+ and Pw1+/pat- males (r = 0.875, **p<0,01 and r = 0.644, p = 0.118, respectively).

An insulin tolerance test (ITT) was performed in a set of Pw1+/+ and Pw1+/pat- male littermates, which revealed no differences between genotypes at 3 months of age, whereas the Pw1+/pat- males developed a modest insulin resistance at 6 months of age as compared to Pw1+/+ males (Fig 2E). Oral glucose tolerance test (OGTT) on the animals of the same litter showed that insulin secretion and glucose clearance were slightly lower in Pw1+/pat- at 3 months of age (Fig 2F). Notably, these patterns were inverted at 6 months of age and Pw1+/pat- males displayed a higher insulin secretion and clearance.

Pw1 reporter gene expression was high in pancreatic β-cells and in hepatocytes (S2A and S2B Fig), both of which were characterized by the presence of sex steroid hormone receptors [39,40]. Co-localization of Pw1 with ERα in various endocrine cells indicates a pivotal role of Pw1 in these cell types via sex steroid signaling.

Paternal loss of Pw1 has been shown to lead to increased β-cells cycling in Pw1+/pat- males at 3 months of age [41]. The increase of proliferation at a younger age may result in increased insulin production in later adulthood. We analyzed the insulin content of pancreas in mature adult males and observed that the pancreatic insulin is slightly elevated in the Pw1+/pat- males (S1B Fig, top). In addition, random-fed glycemia was significantly elevated in Pw1+/pat- males, in agreement with their insulin resistance in adulthood (S1B Fig, bottom).

Pw1+/pat- animals also exhibited age- and sex-dependent hepatic phenotypes: liver size was significantly higher in Pw1+/pat- males at 6 months of age as compared to Pw1+/+ males and females (Fig 2G, left). Gene expression of major adipogenic genes in these animals demonstrated significant changes in the 6-month-old Pw1+/pat- livers in a sex-dependent manner (Fig 2G, right). Notably, the two major isoforms of Pparg1 and Pparg2, differentially expressed between males and females [42], were differently affected by Pw1 loss of function. While the Pparg1 is similarly expressed in Pw1+/+ and Pw1+/pat- livers, Pparg2 expression levels were significantly increased in 6-month-old male livers. In contrast, Pparg2 levels were significantly lower in female livers as compared to male livers at 3 months of age and no increase was observed in Pw1+/pat- female livers at 6 months of age.

PPARG2 is selectively increased in human obesity [43] and is specifically elevated in the steatotic livers of ob/ob mice [44]. We therefore performed hepatic histology using Oil Red-O staining on the 8-month-old livers of both sexes (Figs 2H and S1C). Pw1+/+ livers revealed multiple small lipid droplets in both sexes. In contrast, Pw1+/pat- mice showed abundant, large lipid droplets that were more marked in males. Digital quantification revealed that the total number of lipid droplets and Oil Red-O positive area size were significantly increased in Pw1+/pat- livers (p<0.001) in mature adulthood (Fig 2H). We note that smaller droplets are more abundant in Pw1+/+, whereas larger droplets increased by age in Pw1+/pat- livers, and that this trend was more pronounced in males (S1C Fig).

Taken together, our findings demonstrated that paternal Pw1 deficiency affected multiple stages of early life that proceeded to an age- and sex-dependent global shift of body metabolism towards accelerated adiposity, diabetic-like insulin resistance, and fatty liver in later adulthood, and the impact is more profound in males.

Paternal Pw1 deficiency reduces aggressive behavior and social dominance in males

During routine handling of the Pw1 mutant colony, we observed that adult Pw1+/pat- males seldom display typical aggressive behavior as compared to their Pw1+/+ littermates. We scored incidents of spontaneous fights among Pw1+/+ (n = 75) and Pw1+/pat- (n = 57) male offspring that were group-caged with littermates, and found that Pw1+/pat- males were significantly less aggressive (S3A Fig). When male offspring were separated according to genotype at the time of weaning (Pw1+/+ or Pw1+/pat-), we observed little incidents of fight in the Pw1+/pat- cages, suggesting that the reduced aggressive behavior is, at least in part, intrinsic to the paternal Pw1 loss of function.

To quantitatively assess the competitive ability of Pw1+/pat- males, we used a social-confrontation tube test [45] on adult offspring derived from Pw1+/pat- breeder males. The first test was to examine whether Pw1 is involved in establishing social hierarchy among littermates by using litters consisting of two genotypes. We observed a typical social dominance pattern in which Pw1+/+ males dominate the siblings in a given cage at 10 months of age (S3B Fig, squared in red). Each animal was ranked within each litter by the number of wins and the score was compared between genotypes. This ranking revealed that the Pw1+/+ males rank higher and there is a significant difference between genotype (***p<0.001) (Fig 3A). Notably, the same analysis on younger litters at 3 months of age revealed no significant difference in the inter-litter rank between genotypes, suggesting that younger males have not yet established social rank at this age. A second test was performed in the context of stranger encounter as described by Garfield et. al. [46] in which animals were tested against unfamiliar opponents from different cages with mixed genotypes. The winning rate was determined by the percentage of win in all matches against unknown opponents (S3B Fig). This test demonstrated that the Pw1+/+ males have a greater likelihood of winning in a forced encounter (*P<0.05) (Fig 3B).

thumbnail
Fig 3. Altered social behavior and brain architecture in Pw1+/pat- males.

(A) Interlitter social rank by tube test in Pw1+/+ and Pw1+/pat- males from mixed genotypes at 10 months of age (n = 11 vs n = 10, from 5 litters) and at 3 months of age (n = 7 vs n = 10, from 4 litters). (B) Assessment of social dominance in the stranger encounter tube test. Animals used were listed in S3B Fig. The winning rate was calculated from 17–18 matches against unfamiliar opponents. (C) Pw1 reporter expression (β-gal+) is observed in the vasopressinergic neurons of PVN (top), whereas β-gal signals are strongly co-localizing with ERα receptor in SON (bottom) in the hypothalamus (x400). (D) Representative images of AVP expressing neurons in the PVN of hypothamamus in the Pw1+/+ and Pw1+/pat- male brains (x40). Coronal sections at 120μm intervals through PVN from anterior to posterior axis were immunostained with an anti-AVP antibody. AVP positive area size (dotted line) and cell count were quantified from five sequential sections. E. Digital quantification of AVP+ area size and the total cell count in PVN and their positive correlation in the Pw1+/+ male brain. Columns, mean; bars, SEM; *, P < 0.05, Mann-Whitney U test.

https://doi.org/10.1371/journal.pgen.1010003.g003

In the female Pw1+/p- brains, the oxytocinergic architecture appeared under-developed concomitant with alteration in maternal care [28], a female specific behavior that is acquired at pregnancy. On the other hand, the aggressive behavior commonly observed in laboratory mice is male specific and develop during postnatal growth period. Oxytocin and arginine-vasopressin (AVP) are the two major neuropeptide that regulates sex-specific mammalian behaviors (reviewed in [47,48]). In particular, the AVP system is androgen-dependent [49] and central AVP plays a pivotal role in inter-male aggressive behavior [50,51]. Pw1 is shown to be expressed in both oxytocinergic and vasopressinergic neurons [52]. Therefore, we hypothesized that Pw1 plays a pivotal role in regulating the function of these cell types through sex hormone signaling. We first examined Pw1 expression in the brain using the Pw1 reporter transgenic mouse line Pw1IRESnLacZ [17]. As predicted, we found high levels of reporter gene expression in brain nuclei known to be sexually dimorphic and express sex steroid hormone receptors [53,54], including paraventricular nucleus (PVN) of hypothalamus, the bed nucleus of stria terminalis (BnST), the medial preoptic area (mPOA), and the medial amygdala (MeA) (S3C Fig). These regions are primary sites of AVP production and vasopressinergic neuronal projections [50,55]. We therefore examined the brains from Pw1 reporter mice by immunofluorescence using anti-β-gal and anti-AVP antibodies and found that the vasopressinergic cells in the PVN and SON are the sites of high Pw1 reporter gene expression (Fig 3C, top) which co-express ERα (Fig 3C, bottom), suggesting a role of Pw1 in this cell type. We next examined the architecture of AVP+ cells in the Pw1+/+ and Pw1+/pat- males whose competitive ability had been already established by the tube test (Litter 1–6 in S3B Fig). Using anti-AVP antibody, we immunostained the coronal sections of entire brain and the total AVP+ cell number in the PVN and its area size were determined (Fig 3D). Concordant with the reduced social competitiveness, the AVP+ PVN area size was significantly reduced in the Pw1+/+ brains at 10 months of age (Fig 3E). In the Pw1+/+ brain, we found a strong positive correlation between the area size and the cell number (Fig 3E). Remarkably, this correlation is abolished in the Pw1+/pat- brain, implying that the proliferation and/or expansion of the AVP+ cells are deregulated. Finally, we examined the correlation between the AVP+ cell structure and social behavior in the litter 1 which consists of four Pw1+/+ males. Both area size and cell count in the PVN showed positive correlation with the winning rate in this set of animals (S3D Fig).

These data suggests that Pw1 promotes acquired social dominance and aggressive behavior by modulating the AVP+ neuroendocrine architecture in male mice.

Pw1 promotes testosterone production in young male mice

Perinatal androgen secretion leads to changes in the CNS and underlie sexual dimorphism in the brain [56]. The testis produce androgens in adolescence that contribute to the development of adult male characteristics in metabolism and behavior that typically underlies reproductive success including male mating and aggression [5759]. We therefore measured levels of testosterone in young males during peripubertal development when testosterone secretion peaks in postnatal growth (between 6–9 weeks)[60]. First, we measured testosterone levels of male Pw1+/+ and Pw1+/pat- littermates derived from 3 litters and found that the Pw1+/pat- mice display a delayed pubertal surge (Fig 4A). Based on this timing, we further analyzed different young males from multiple cages containing mixed genotypes in the entire colony at the age between 2.5 and 3 month old. These analyses revealed that Pw1+/pat- males have significantly lower levels of testosterone as compared to Pw1+/+ males (Fig 4B). In agreement with these observations, mRNA expression of steroidogenic genes encoding the rate limiting enzymes for testosterone biosynthesis (Cyp17a, Cyp11a, and 3β-HSD), as well as that of luteinizing hormone receptor (LHR) in the testes are significantly reduced in Pw1+/pat- males (Fig 4C). Notably, there was no significant difference in the expression of aromatase Cyp19 or 17β-HSD3.

thumbnail
Fig 4. Reduced testosterone production in the Pw1+/pat- males.

(A) Plasma testosterone measurements at peripubertal age in Pw1+/+ and Pw1+/pat- males (n = 7 and 5, respectively, from 3 litters). Two-way repeated measures ANOVA with Sidak’s multiple comparisons test revealed a significant difference between genotypes at 9 weeks of age (*P<0.05). (B) Plasma testosterone levels in Pw1+/+ and Pw1+/pat- non-breeder littermates (n = 47 and n = 42) at 2.5–3 months of age. Each symbol represents independent animals. Bars represent mean ± SEM; **, P<0.01, by Mann Whitney U test. (C) mRNA expression of steroidogenic genes in the 2.5–3 month-old testes of Pw1+/+ and Pw1+/pat- non-breeder littermates (n = 15 versus n = 11). StAR, Steroidogenic acute regulatory protein; Cyp17a1, cytochrome P-450 17a; Cyp11a1, cholesterol side-chain cleavage enzyme; Hsd3b1, 3-b-hydroxysteroid dehydrogenase; Lhr, luteinizing hormone receptor; Cyp19a1, aromatase enzyme; Hsd17b3, 17-b-hydroxysteroid dehydrogenase. (D) Pw1 reporter expression in the 2.5-month-old testis of Pw1IRESnLacZ mice, showing X-gal staining (x100) and β-gal immunofluorescence (x400) in the interstitial compartment and in the epithelium of seminiferous tubules. The β-gal immunofluorescence overlapped with that of androgen receptor (AR) in the Leydig and Sertoli cells (arrows). L, interstitial Leydig cells; S, seminiferous tubule.

https://doi.org/10.1371/journal.pgen.1010003.g004

A previous study from our laboratory using the Pw1 reporter transgenic mouse line (Pw1IRESnLacZ) demonstrated that Pw1 is highly expressed in the peritubular cells near the basement membrane of seminiferous tubules, a part of which were identified as Bmi1+ spermatogonia [17]. Since the Pw1+/pat- mice displayed reduced steroidogenesis, we were interested in whether Pw1 is expressed in the cell types responsible for the testosterone production in the testis. We performed histological analyses on the testis of Pw1IRESnLacZ mice using an anti-β-gal antibody and antibodies against sex steroid hormone receptors and found that reporter gene expression colocalized with androgen receptor expression (Fig 4D). In the testis, Leidig cells and Sertoli cells are the two cell types predominantly express androgen receptor (AR) [61]. Therefore, we conclude that Pw1 reporter gene is abundantly expressed in the testosterone producing Leydig cells in the interstitial compartment and in the supporting Sertoli cells near the basement membrane of seminiferous tubules. We note that other imprinted genes are highly expressed in the Leidig cells and that a number of imprinted genes are simultaneously deregulated in Leidig cells of human patients with idiopathic germ cell aplasia [62].

Taken together, our data demonstrate that male sex-hormone signaling during secondary sexual maturation is suppressed in Pw1+/pat- male mice, concomitant with reduced gene expression for testosterone biosynthesis in the testis, which may account for the reduced masculinization in metabolism and social behavior during the postnatal development.

Discussion

We demonstrate that Pw1 plays a key role at specific stages of postnatal development for sex-specific growth, metabolism and behavior. Pw1 deficient male mice exhibit a smaller body size and a reduced masculinization of body composition followed by a global shift of metabolism leading to early onset obesity and related metabolic changes. A marked reduction of growth promoting GH/IGF-1 and insulin characterized early postnatal life and a decrease in testosterone activity at puberty coupled with insulin-resistance lead to a male-specific deficiency by adulthood (Fig 5). These observations are concomitant with a significantly reduced sexual dimorphism in body composition. In addition, our data suggests a model in which a paternally expressed gene promotes male-specific brain development and behavior through a sex steroid pathway that contribute to reproductive success in mammals. Our finding support the intralocus sexual conflict model of genomic imprinting [6] for the control of mammalian sexual dimorphism.

thumbnail
Fig 5. Hormonal cascades are deregulated in Pw1+/pat- mice.

Pw1 is abundantly expressed in the gonads and target organs and colocalizes with steroid hormone receptors in hormone secreting cells. GHRH, growth hormone-releasing hormone; GH, growth hormone; IGF-1, insulin-like growth factor 1; GnRH, gonadotropin-releasing hormone; LH, luteinizing hormone; FSH, follicle-stimulating hormone; AVP, arginine vasopressin; +, stimulation; -, inhibition by Pw1.

https://doi.org/10.1371/journal.pgen.1010003.g005

The physiopathological alterations in Pw1 deficient mice identified in this study are similar to the phenotypic spectrum of Prader-Willi syndrome (PWS), a disorder with hypothalamic dysfunction and hypopituitarism due to imprinting errors [63]. PWS is characterized by biphasic clinical manifestations, i.e., reduced growth velocity, hypoglycemia, and hypotonia at infancy, followed by hyperphagia, extreme obesity, and hypogonadism in childhood through adulthood. In addition, patients present low levels of GH and IGF-1, gonadotropins, and gonadal sex steroids, hence, the established treatment includes GH administration at infancy and sex hormone replacement in young adult that improve their growth and metabolism by enhancing muscle mass and reducing fat mass [64].

We found that the expression and function of Pw1 are particularly concentrated in the hypothalamus, pituitary gland, and gonads. Accordingly, our findings suggest that a primary Pw1 deficiency involves hypothalamic-pituitary-gonadal dysfunction at critical periods of postnatal development, leading to changes in sex steroids levels and GH/IGF signaling. Our observations are in line with several studies reporting that mouse mutants for paternally expressed genes typically display hypothalamo-pituitary phenotypes [65,66]. For instance, a loss of Grf1 or Dlk1 resulted in reduced GH content and secretion [67,68], while disruption of the imprinting domain encompassing Dlk1, Rtl1, and Dio3 resulted in transient perturbation in the GH/IGF-1 growth pathway with hypothyroidism [69]. Together, these observations pointed that imprinted genes play pivotal roles for the regulation of hypothalamo-pituitary-gonadal (HPG) axis governing growth, metabolism and reproduction.

Sexual dimorphism in mammalian growth is regulated by the sex steroid-dependent GH/IGF-1 pathway. During postnatal development, testosterone and 17-estradiol evoke pituitary production of GH in discrete manners [70,71]. GH in turn stimulates systemic IGF-1 through the JAK2/Stat5b pathways [72,73]. Estrogen-bound estrogen receptor (ER) also induces transcription of IGF-1 and IGF receptor through estrogen-responsive elements at their gene loci [74]. The activation of IGF-1/IGF1R signaling further activates transcription factors including ERs, forming a complex crosstalk between IGF-1 and sex steroid signaling that amplifies somatic growth in puberty [75]. Notably, sex differences during postnatal growth are diminished in the Stat5b deficient mice [76], the GH receptor (Ghr)- and Igf1-deficient mice [77]. Hence, reduced sexual dimorphism in Pw1 deficient young is likely due to an attenuated hypothalamic regulation of sex hormones and downstream GH/IGF-1 signaling.

Androgens and their cognate receptors regulate male secondary sexual differentiation [78], and testosterone deficiency is related to metabolic syndromes in men [79,80]. In androgen receptor deficient (ARKO) mice, late-onset obesity and fatty liver, as well as insulin- and leptin-resistance have been reported which is restricted to males [8183]. Elevated levels of circulating leptin and leptin-resistance were also reported previously in the Pw1+/pat- adult males [31]. Leptin was shown to reciprocally regulate gonadal hormone production at central and peripheral levels [84] and the steroidogenesis in the Leydig cells was impaired in the leptin-deficient ob/ob males [85] concomitant with loss of skeletal masculinization [86]. Therefore, the defect in leptin signaling in Pw1+/pat- males may causally relate to further reduction in gonadal steroidogenesis. We note, however, that the Pw1 deficient phenotypes were not restricted to males, as both Pw1+/pat- males and females displayed delayed postnatal growth until weaning. Smaller body size of Pw1+/pat- mice before weaning may originate, at least in part, from changes arising during prenatal development, as observed by the lower body weight and myofiber numbers at birth [19]. Sex specificity of resulting phenotypes upon Pw1 loss is presumably due to sex dimorphic postnatal dependency to sex steroid signaling in males and females. Sex-dependent phenotypes in Pw1 deficiency have also been demonstrated in placental function, where female placenta compensates the loss of Pw1 better than male placenta [30,87]. This indicates that the sex-dependent Pw1 function starts much earlier than the offspring’s postnatal manifestation of sexually dimorphism in growth.

Testosterone exposure at specific phases of development, such as perinatal and pubertal periods, is crucial for the brain masculinization [56,58]. This study shows that Pw1 expression and sex steroid receptors highly overlap in key endocrine and neural tissues that regulate sex-specific behaviors. Specifically, we found that paternal loss of Pw1 function attenuated social dominance and correlates with the decreased expansion of the testosterone-sensitive vasopressinergic nuclei in the hypothalamus. In general, the sexually dimorphic nuclear volume/neuron number is ascribed to the sex steroid control of proliferation and cell death [88]. Given that Pw1 regulates cell survival and/or apoptosis at a cellular level, it is conceivable that Pw1 participates in sex-specific brain formation and neuronal proliferation/differentiation for various acquired behaviors. In contrast, aggression in males is ERα dependent [89] and the inter-male competitive ability reflects the level of ERα expression in BnST [45]. Pw1 may contribute to masculinization by regulating the expression of cognate receptors in sex-hormone dependent cells, thereby determining the dependency of organisms on sex steroid signaling.

Our finding suggests that Pw1 plays a role in regulating the proliferative and/or cell number in various tissues through sex steroid signaling during tissue remodeling in adult life. We demonstrate further that Pw1 expressing adult stem/progenitor cells are characterized by the presence of sex steroid receptors. Consistent with this observation, the number of Pw1 expressing cells in the heart were shown to increase during pregnancy [90], while the pregnancy-induced expansion of pancreatic β-cells inversely correlates with the levels of Pw1 expression [41]. Accumulating evidence also demonstrates the interaction of sex hormones in stem cell behaviors. Kim et. al. demonstrated that sex steroid hormones promote the establishment of the adult quiescent satellite cell pool through the HPG axis [91]. Likewise, ERα signaling has been shown to drive pancreatic β-cell replication in development as well as in post-injury regeneration [92]. How Pw1 loss of function alters adult stem/progenitor cell behaviors or steroid hormones per se affects pluripotency is beyond the scope of this study, however, we have shown previously that loss of Pw1 function leads to a loss in self-renewal capacity in adult muscle stem cells [19] as well as skin stem cell [17]. The link between self-renewal, stem cell competence and sex steroid hormones presents an emerging concept in stem cell biology.

This study uncovers a key role of the parentally imprinted gene Pw1 in establishing male-specific characteristics in metabolism, brain structure and behavior by modulating sex steroid pathways. We note that similar metabolic phenotypes observed in Pw1 deficient mice were reported in the progeny of human populations that have been subjected to overnutrition or prolonged nutritional deprivation [93], coinciding with a deregulation of several parentally imprinted genes [94]. Our findings show that Pw1 provides an epigenetic context that links energy metabolism to the generation of male-specific traits in mice. Analogous mechanisms may also underlie the sex differences in the metabolic syndrome in humans.

Materials and methods

Ethics statement

All experiments were in adherence to the institutional guideline for experiment and husbandry of laboratory animals. Approval for the animal (mouse) work performed in this study was obtained through review by the French Ministry of Education, Agreement #A751320.

Animals

The Pw1IRESnLacZ reporter mice [17] and the Pw1 deficient mice carrying targeted mutation Pw1Δ9 [29] were generated in our laboratory and maintained in C57BL/6J genetic background. Eight Pw1+/pat- males (two generation after the first deletion of Pw1 allele, designated F3) were used as breeders. In order to compare Pw1+/pat- offspring to their Pw1+/+ littermates, Pw1+/pat- breeder males under the age of 8 months old were crossed with young C57BL/6J females (purchased from Janvier Laboratory, Le Genest St Isle, France), and large litters with both genotypes were pooled for each analysis whenever available. All metabolic phenotypes presented in this study were obtained from offspring of F3-F4, and behavior phenotypes from F5 generations. Offspring were weaned at the age between 3–4 weeks, and housed in groups with littermates unless otherwise stated (maximum 6 per cage). Mice were kept on 12-hour light/dark cycle at 24°C with ad libitum access to water and standard chow diet.

Postnatal growth and body composition

Postnatal growth of Pw1+/pat- mice was measured along with Pw1+/+ littermates, which were nursed together by C57BL/6J females. Longitudinal body composition analysis was realized using Brucker minispec nuclear magnetic resonance (NMR) imaging (Brucker, USA), starting from 1 month of age. The lean and fat mass expressed in percentage (%) was converted to weight (g) in respect to the total body mass (g).

Blood glucose and food intake

Blood was sampled from the tail tip and glycemia was measured using an Accu-Check glucometer with disposable test strips (Roche Diagnostics, Basel, Switzerland). Food intake was assessed on N = 5–7 mice per genotype per sex, individually housed at 8 weeks of age and food consumption was measured twice a week for 3 weeks.

Plasma hormone measurement

Blood was collected either from the tail tip or from the facial vein of randomly fed animals at dormant period of circadian rhythms (11am-16pm), using Microvette CB300 capillary action blood collection tubes (Sarstedt, Nümbrecht, Germany). Plasma IGF, insulin, and testosterone levels were determined using ELISA kits for mouse (Chrystal Chem. Ltd, NJ, USA).

Glucose and insulin tolerance test

Oral glucose tolerance test (OGTT) and insulin tolerance test (ITT) were performed on 5–6 hours fasted mice. D-glucose (2 g/kg) was orally administered, and human insulin (0.25–0.55 units/kg) was intraperitoneally injected. For both tests, tail blood was sampled at 0, 15, 30, 60, and 90 minutes after the administration, and glycemia was immediately measured as described above. For monitoring insulin secretion upon OGTT, tail blood was collected at 0, 15, 30 minutes and subjected for ELISA protocol.

Quantification of fatty liver

To visualize lipid droplets, Oil red O-staining was performed in the 10um frozen sections of 8-month-old liver. Number of lipid droplets and the total oil Red positive area were quantified using ‘particle measurement’ in the Image-J software [95].

Gene expression analysis

Total RNA was isolated using TRIzol reagent (Ambion) and subjected to cDNA synthesis (Lifetechnology), followed by quantitative PCR analysis (LightCycler, Roche). The expression levels of target genes were normalized with the levels of a housekeeping gene coding TATA binding protein (Tbp). Primer sequences are available upon request.

Assessment of social dominance by confrontation tube test

We applied the social confrontation tube test previously described by Garfield et. al. on Grb10 deficient mice [46] to determine social dominance among Pw1+/+ and Pw1+/p- males. All mice used in the test were group-housed with littermates consist of both Pw1+/+ and Pw1+/p- genotypes (n = 4–6) at mature age between 10–12 months. Independent sets of mice were analyzed at 3 months old for comparison.

To determine the social hierarchy within a given litter, animals from a given litter were challenged with their littermate. On subsequent days, the same animals were also challenged with unfamiliar opponents from different litters. On the day of experiment, the animals were removed from home cages, and isolated for the duration of the test. All the tests were performed during daylight hours. The test apparatus was a 32-cm semi-transparent tube with adjustable internal diameter to prevent crossing over or turning of animals. Test animals were placed opposite ends of the tube and released simultaneously. Winners and losers were scored as one animal retreated from the tube completely. Tests in which the animal remained in the tube (five minutes maximum) were scored as draw. Inter-litter rank was determined by the number of winning within littermates (most win to least win = from 1, 2, 3, etc.).

Histology and Immunohistochemistry

X-Gal staining and immunohistochemistry were performed following standard protocols described previously [17]. All tissues were fixed in 4% paraformaldehyde and embedded in Optimal Temperature Compound (OTC) after cryoprotection in 20–30% sucrose. 12 μm cryosections were used for X-gal staining and co-immunofluorescences. The Pw1nLacZ reporter gene expression faithfully identified the Pw1 expressing cells [17], therefore, mouse anti-β-galactocidase antibody (1:500, Z3781: Promega) was used to represent Pw1 expression in combination with either rabbit anti-AR (1: 1000, Santa-cruz, sc-816), anti-ERα (1:1000, Santa-cruz, sc-542), or AVP (1:1000, AB1565: Millipore) in co-localization study. For brain structure analysis, mice were anesthetized with Ketamine/Xylazine and transcardially perfused with 4% paraformaldehyde. Brains were post-fixed overnight and 30μm serial sections were obtained from the entire brains using a vibratome. Comparable regions consist of 13–16 coronal cross-sections. The brain section images were collected from the same brain coordinates (Bregma: between -0.70 to -0.94 mm). Typically the AVP+ PVN structure spanning in four sections were subjected for analysis. Floating sections were incubated with rabbit anti-AVP (1, 1000, AB1565: Millipore) and detected using Vectastain ABC HRP kit (Vector Labs). Number of AVP+ cells and the total area were quantified using ImageJ software [95].

Statistical analysis

All figures and statistical analyses were generated using Prism 5 (GraphPad). Comparison between genotypes and sexes was performed using two-way ANOVA with Tukey’s multiple comparison test. For male specific variables, Mann-Whitney U-test was applied for independent measurement. Correlation between two variables was tested by a simple linear regression analysis. P value < 0.05 was considered statistically significant.

Supporting information

S1 Fig. Insulin action in Pw1+/pat- adult males.

(A) Positive correlation of insulin levels with lean mass (top) and with body fat (bottom) was observed at at 3 months and 6 months of age, respectively. Correlation was determined by simple linear regression analysis. (B) Pancreatic insulin content and blood glucose levels at 10 months of age. Pancreatic insulin content was measured using acid-ethanol extraction protocol, followed by the insulin ELISA and normalized with total protein content. Columns, mean; bars, SEM; **, P< 0.01 by Mann-Whitney U test. (C) Typical images of Oil Red-O stained 8-month-old livers (top) and the quantification of lipid droplet number and size distribution (bottom), illustrating the differences in lipid content between the groups. P<0.001 by two-way ANOVA with Tukey’s multiple comparisons. Original magnifications: x200.

https://doi.org/10.1371/journal.pgen.1010003.s001

(TIF)

S2 Fig. Pw1 colocalisation with sex hormone receptors.

Pw1 reporter expression predominantly co-localizing with the nuclear expression of sex hormone receptors in diverse cell types. (A) pancreatic islets, (B) mono- and dinucleated hepatocytes, (C) anterior pituitary cells, (D) adipocytes, and (E) skeletal muscle.

https://doi.org/10.1371/journal.pgen.1010003.s002

(TIF)

S3 Fig. Reduced aggressive behavior and social rank of Pw1+/pat- adult males.

(A) Spontaneous fights among group-housed littermates were monitored per cage and individuals involved in fights were identified. ***, P< 0.001 by Fisher’s exact test. (B) Result of social confrontation tube test in Pw1+/+ (+/+) and Pw1+/pat- (pat-) littermates from six litters. Pw1+/pat- males (n = 10) and Pw1+/+ littermates (n = 13) were subjected for confrontation against each other and their i) intra-litter rank and ii) winning rate against unfamiliar opponents was determined for each animal. Squared in red on the diagonal line show matches within littermates. Mice were derived from C57B6 x Pw1+/p- breeding except for Litter 1, which was derived from C57B6 x Pw1+/+ breeding). (C) Coronal sections of Pw1IRESnLacZ transgenic brain at 2.5 months of age revealing Pw1 reporter gene expression (X-gal staining) in sexually dimorphic brain regions. BnST, the bed nucleus of the stria terminalis; mPOA, medial preoptic area; PVN, paraventricular nucleus of hypothalamus; MeA, medial amygdala; PIR, piriform cortex (x40). (D) Size comparisons of AVP+ PVN area in four Pw1+/+ male siblings from C57B6 x Pw1+/+ breeding (Litter 1). Positive correlation was found between winning rate and AVP+ cells.

https://doi.org/10.1371/journal.pgen.1010003.s003

(TIF)

Acknowledgments

We thank Drs. AL. Denizot, RM. Correra, D. Ollitrault, K. Kyrylkova, S. Kyryachenko and JR. Courbard for scientific discussions, and Ms I. Lopez for technical assistance. We thank Ms. A. Lacombe and Dr. T. Huby in the PreclinICAN, Institute of Cardiometabolism and Nutrition (IHU-ICAN, ANR-10-IAHU-05) for providing facility and technical advice for the metabolic analysis. We also thank the animal facility of Pierre and Marie Curie University for their excellent management and maintenance of laboratory animals.

References

  1. 1. Renfree MB, Suzuki S, Kaneko-Ishino T. The origin and evolution of genomic imprinting and viviparity in mammals. Philos Trans R Soc Lond B Biol Sci. 2013;368(1609):20120151. pmid:23166401; PubMed Central PMCID: PMC3539366.
  2. 2. Spencer HG, Clark AG. Non-conflict theories for the evolution of genomic imprinting. Heredity (Edinb). 2014;113(2):112–8. pmid:24398886; PubMed Central PMCID: PMC4105448.
  3. 3. Patten MM, Ross L, Curley JP, Queller DC, Bonduriansky R, Wolf JB. The evolution of genomic imprinting: theories, predictions and empirical tests. Heredity (Edinb). 2014;113(2):119–28. pmid:24755983; PubMed Central PMCID: PMC4105453.
  4. 4. Haig D, Graham C. Genomic imprinting and the strange case of the insulin-like growth factor II receptor. Cell. 1991;64(6):1045–6. pmid:1848481.
  5. 5. Curley JP, Barton S, Surani A, Keverne EB. Coadaptation in mother and infant regulated by a paternally expressed imprinted gene. Proc Biol Sci. 2004;271(1545):1303–9. pmid:15306355; PubMed Central PMCID: PMC1691726.
  6. 6. Day T, Bonduriansky R. Intralocus sexual conflict can drive the evolution of genomic imprinting. Genetics. 2004;167(4):1537–46. pmid:15342496; PubMed Central PMCID: PMC1470977.
  7. 7. Gregg C, Zhang J, Butler JE, Haig D, Dulac C. Sex-specific parent-of-origin allelic expression in the mouse brain. Science. 2010;329(5992):682–5. pmid:20616234; PubMed Central PMCID: PMC2997643.
  8. 8. Faisal M, Kim H, Kim J. Sexual differences of imprinted genes’ expression levels. Gene. 2014;533(1):434–8. pmid:24125951; PubMed Central PMCID: PMC3872986.
  9. 9. Peters J. The role of genomic imprinting in biology and disease: an expanding view. Nat Rev Genet. 2014;15(8):517–30. pmid:24958438.
  10. 10. Cleaton MA, Edwards CA, Ferguson-Smith AC. Phenotypic outcomes of imprinted gene models in mice: elucidation of pre- and postnatal functions of imprinted genes. Annu Rev Genomics Hum Genet. 2014;15:93–126. pmid:24898037.
  11. 11. Tucci V, Isles AR, Kelsey G, Ferguson-Smith AC, Erice Imprinting G. Genomic Imprinting and Physiological Processes in Mammals. Cell. 2019;176(5):952–65. pmid:30794780.
  12. 12. Babak T, DeVeale B, Tsang EK, Zhou Y, Li X, Smith KS, et al. Genetic conflict reflected in tissue-specific maps of genomic imprinting in human and mouse. Nat Genet. 2015;47(5):544–9. pmid:25848752; PubMed Central PMCID: PMC4414907.
  13. 13. Baran Y, Subramaniam M, Biton A, Tukiainen T, Tsang EK, Rivas MA, et al. The landscape of genomic imprinting across diverse adult human tissues. Genome Res. 2015;25(7):927–36. pmid:25953952; PubMed Central PMCID: PMC4484390.
  14. 14. Patten MM, Cowley M, Oakey RJ, Feil R. Regulatory links between imprinted genes: evolutionary predictions and consequences. Proc Biol Sci. 2016;283(1824). pmid:26842569; PubMed Central PMCID: PMC4760173.
  15. 15. Kuroiwa Y, Kaneko-Ishino T, Kagitani F, Kohda T, Li LL, Tada M, et al. Peg3 imprinted gene on proximal chromosome 7 encodes for a zinc finger protein. Nat Genet. 1996;12(2):186–90. pmid:8563758.
  16. 16. Relaix F, Weng X, Marazzi G, Yang E, Copeland N, Jenkins N, et al. Pw1, a novel zinc finger gene implicated in the myogenic and neuronal lineages. Dev Biol. 1996;177(2):383–96. pmid:8806818.
  17. 17. Besson V, Smeriglio P, Wegener A, Relaix F, Nait Oumesmar B, Sassoon DA, et al. PW1 gene/paternally expressed gene 3 (PW1/Peg3) identifies multiple adult stem and progenitor cell populations. Proc Natl Acad Sci U S A. 2011;108(28):11470–5. pmid:21709251; PubMed Central PMCID: PMC3136256.
  18. 18. Nicolas N, Marazzi G, Kelley K, Sassoon D. Embryonic deregulation of muscle stress signaling pathways leads to altered postnatal stem cell behavior and a failure in postnatal muscle growth. Dev Biol. 2005;281(2):171–83. pmid:15893971.
  19. 19. Correra RM, Ollitrault D, Valente M, Mazzola A, Adalsteinsson BT, Ferguson-Smith AC, et al. The imprinted gene Pw1/Peg3 regulates skeletal muscle growth, satellite cell metabolic state, and self-renewal. Sci Rep. 2018;8(1):14649. pmid:30279563; PubMed Central PMCID: PMC6168517.
  20. 20. Bonfanti C, Rossi G, Tedesco FS, Giannotta M, Benedetti S, Tonlorenzi R, et al. PW1/Peg3 expression regulates key properties that determine mesoangioblast stem cell competence. Nat Commun. 2015;6:6364. pmid:25751651; PubMed Central PMCID: PMC4366533.
  21. 21. Besson V, Kyryachenko S, Janich P, Benitah SA, Marazzi G, Sassoon D. Expression Analysis of the Stem Cell Marker Pw1/Peg3 Reveals a CD34 Negative Progenitor Population in the Hair Follicle. Stem Cells. 2017;35(4):1015–27. pmid:27862634.
  22. 22. Relaix F, Wei XJ, Wu X, Sassoon DA. Peg3/Pw1 is an imprinted gene involved in the TNF-NFkappaB signal transduction pathway. Nat Genet. 1998;18(3):287–91. pmid:9500555.
  23. 23. Relaix F, Wei X, Li W, Pan J, Lin Y, Bowtell DD, et al. Pw1/Peg3 is a potential cell death mediator and cooperates with Siah1a in p53-mediated apoptosis. Proc Natl Acad Sci U S A. 2000;97(5):2105–10. pmid:10681424; PubMed Central PMCID: PMC15761.
  24. 24. Deng Y, Wu X. Peg3/Pw1 promotes p53-mediated apoptosis by inducing Bax translocation from cytosol to mitochondria. Proc Natl Acad Sci U S A. 2000;97(22):12050–5. pmid:11050235; PubMed Central PMCID: PMC17292.
  25. 25. Buraschi S, Neill T, Goyal A, Poluzzi C, Smythies J, Owens RT, et al. Decorin causes autophagy in endothelial cells via Peg3. Proc Natl Acad Sci U S A. 2013;110(28):E2582–91. pmid:23798385; PubMed Central PMCID: PMC3710796.
  26. 26. Thiaville MM, Huang JM, Kim H, Ekram MB, Roh TY, Kim J. DNA-binding motif and target genes of the imprinted transcription factor PEG3. Gene. 2013;512(2):314–20. pmid:23078764; PubMed Central PMCID: PMC3513644.
  27. 27. Frey WD, Sharma K, Cain TL, Nishimori K, Teruyama R, Kim J. Oxytocin receptor is regulated by Peg3. PLoS One. 2018;13(8):e0202476. Epub 2018/08/15. pmid:30106994; PubMed Central PMCID: PMC6091971.
  28. 28. Li L, Keverne EB, Aparicio SA, Ishino F, Barton SC, Surani MA. Regulation of maternal behavior and offspring growth by paternally expressed Peg3. Science. 1999;284(5412):330–3. pmid:10195900.
  29. 29. Denizot AL, Besson V, Correra RM, Mazzola A, Lopes I, Courbard JR, et al. A Novel Mutant Allele of Pw1/Peg3 Does Not Affect Maternal Behavior or Nursing Behavior. PLoS Genet. 2016;12(5):e1006053. pmid:27187722; PubMed Central PMCID: PMC4871489.
  30. 30. Kim J, Frey WD, He H, Kim H, Ekram MB, Bakshi A, et al. Peg3 mutational effects on reproduction and placenta-specific gene families. PLoS One. 2013;8(12):e83359. Epub 2014/01/07. pmid:24391757; PubMed Central PMCID: PMC3877027.
  31. 31. Curley JP, Pinnock SB, Dickson SL, Thresher R, Miyoshi N, Surani MA, et al. Increased body fat in mice with a targeted mutation of the paternally expressed imprinted gene Peg3. FASEB J. 2005;19(10):1302–4. pmid:15928196.
  32. 32. Swaney WT, Curley JP, Champagne FA, Keverne EB. Genomic imprinting mediates sexual experience-dependent olfactory learning in male mice. Proc Natl Acad Sci U S A. 2007;104(14):6084–9. pmid:17389373; PubMed Central PMCID: PMC1851620.
  33. 33. Le Roith D, Scavo L, Butler A. What is the role of circulating IGF-I? Trends Endocrinol Metab. 2001;12(2):48–52. pmid:11167121.
  34. 34. Gatford KL, Egan AR, Clarke IJ, Owens PC. Sexual dimorphism of the somatotrophic axis. J Endocrinol. 1998;157(3):373–89. pmid:9691971.
  35. 35. Lichanska AM, Waters MJ. How growth hormone controls growth, obesity and sexual dimorphism. Trends Genet. 2008;24(1):41–7. pmid:18063438.
  36. 36. Callewaert F, Venken K, Kopchick JJ, Torcasio A, van Lenthe GH, Boonen S, et al. Sexual dimorphism in cortical bone size and strength but not density is determined by independent and time-specific actions of sex steroids and IGF-1: evidence from pubertal mouse models. J Bone Miner Res. 2010;25(3):617–26. pmid:19888832.
  37. 37. Menon RK, Sperling MA. Insulin as a growth factor. Endocrinol Metab Clin North Am. 1996;25(3):633–47. pmid:8879990.
  38. 38. Czech MP. Insulin action and resistance in obesity and type 2 diabetes. Nat Med. 2017;23(7):804–14. pmid:28697184.
  39. 39. Mauvais-Jarvis F. Role of Sex Steroids in beta Cell Function, Growth, and Survival. Trends Endocrinol Metab. 2016;27(12):844–55. pmid:27640750; PubMed Central PMCID: PMC5116277.
  40. 40. Shen M, Shi H. Sex Hormones and Their Receptors Regulate Liver Energy Homeostasis. Int J Endocrinol. 2015;2015:294278. pmid:26491440; PubMed Central PMCID: PMC4600502.
  41. 41. Sojoodi M, Stradiot L, Tanaka K, Heremans Y, Leuckx G, Besson V, et al. The zinc finger transcription factor PW1/PEG3 restrains murine beta cell cycling. Diabetologia. 2016;59(7):1474–9. pmid:27130279; PubMed Central PMCID: PMC4901110.
  42. 42. Rando G, Wahli W. Sex differences in nuclear receptor-regulated liver metabolic pathways. Biochim Biophys Acta. 2011;1812(8):964–73. pmid:21211563.
  43. 43. Vidal-Puig AJ, Considine RV, Jimenez-Liñan M, Werman A, Pories WJ, Caro JF, et al. Peroxisome proliferator-activated receptor gene expression in human tissues. Effects of obesity, weight loss, and regulation by insulin and glucocorticoids. J Clin Invest. 1997;99(10):2416–22. Epub 1997/05/15. pmid:9153284; PubMed Central PMCID: PMC508081.
  44. 44. Rahimian R, Masih-Khan E, Lo M, van Breemen C, McManus BM, Dubé GP. Hepatic over-expression of peroxisome proliferator activated receptor gamma2 in the ob/ob mouse model of non-insulin dependent diabetes mellitus. Molecular and cellular biochemistry. 2001;224(1–2):29–37. Epub 2001/11/06. pmid:11693197.
  45. 45. Greenberg GD, Howerton CL, Trainor BC. Fighting in the home cage: Agonistic encounters and effects on neurobiological markers within the social decision-making network of house mice (Mus musculus). Neurosci Lett. 2014;566:151–5. pmid:24602985; PubMed Central PMCID: PMC4019314.
  46. 46. Garfield AS, Cowley M, Smith FM, Moorwood K, Stewart-Cox JE, Gilroy K, et al. Distinct physiological and behavioural functions for parental alleles of imprinted Grb10. Nature. 2011;469(7331):534–8. pmid:21270893; PubMed Central PMCID: PMC3031026.
  47. 47. McCall C, Singer T. The animal and human neuroendocrinology of social cognition, motivation and behavior. Nat Neurosci. 2012;15(5):681–8. pmid:22504348.
  48. 48. Caldwell HK. Oxytocin and Vasopressin: Powerful Regulators of Social Behavior. Neuroscientist. 2017;23(5):517–28. pmid:28492104.
  49. 49. De Vries GJ, Buijs RM, Van Leeuwen FW, Caffe AR, Swaab DF. The vasopressinergic innervation of the brain in normal and castrated rats. J Comp Neurol. 1985;233(2):236–54. pmid:3882778.
  50. 50. Caldwell HK, Lee HJ, Macbeth AH, Young WS 3rd. Vasopressin: behavioral roles of an "original" neuropeptide. Prog Neurobiol. 2008;84(1):1–24. pmid:18053631; PubMed Central PMCID: PMC2292122.
  51. 51. Veenema AH, Neumann ID. Central vasopressin and oxytocin release: regulation of complex social behaviours. Prog Brain Res. 2008;170:261–76. pmid:18655888.
  52. 52. Yamashita M, Glasgow E, Zhang BJ, Kusano K, Gainer H. Identification of cell-specific messenger ribonucleic acids in oxytocinergic and vasopressinergic magnocellular neurons in rat supraoptic nucleus by single-cell differential hybridization. Endocrinology. 2002;143(11):4464–76. pmid:12399444.
  53. 53. Shah NM, Pisapia DJ, Maniatis S, Mendelsohn MM, Nemes A, Axel R. Visualizing sexual dimorphism in the brain. Neuron. 2004;43(3):313–9. pmid:15294140.
  54. 54. Wu MV, Manoli DS, Fraser EJ, Coats JK, Tollkuhn J, Honda S, et al. Estrogen masculinizes neural pathways and sex-specific behaviors. Cell. 2009;139(1):61–72. pmid:19804754; PubMed Central PMCID: PMC2851224.
  55. 55. Dumais KM, Veenema AH. Vasopressin and oxytocin receptor systems in the brain: Sex differences and sex-specific regulation of social behavior. Front Neuroendocrinol. 2016;40:1–23. pmid:25951955; PubMed Central PMCID: PMC4633405.
  56. 56. Sato T, Matsumoto T, Kawano H, Watanabe T, Uematsu Y, Sekine K, et al. Brain masculinization requires androgen receptor function. Proc Natl Acad Sci U S A. 2004;101(6):1673–8. pmid:14747651; PubMed Central PMCID: PMC341816.
  57. 57. Mauvais-Jarvis F. Sex differences in metabolic homeostasis, diabetes, and obesity. Biology of sex differences. 2015;6:14. Epub 2015/09/05. pmid:26339468; PubMed Central PMCID: PMC4559072.
  58. 58. Sato SM, Schulz KM, Sisk CL, Wood RI. Adolescents and androgens, receptors and rewards. Horm Behav. 2008;53(5):647–58. pmid:18343381; PubMed Central PMCID: PMC2435368.
  59. 59. Juntti SA, Coats JK, Shah NM. A genetic approach to dissect sexually dimorphic behaviors. Horm Behav. 2008;53(5):627–37. pmid:18313055; PubMed Central PMCID: PMC2464277.
  60. 60. Hill RA, Wu YW, Kwek P, van den Buuse M. Modulatory effects of sex steroid hormones on brain-derived neurotrophic factor-tyrosine kinase B expression during adolescent development in C57Bl/6 mice. J Neuroendocrinol. 2012;24(5):774–88. pmid:22221196.
  61. 61. Domanskyi A, Zhang FP, Nurmio M, Palvimo JJ, Toppari J, Janne OA. Expression and localization of androgen receptor-interacting protein-4 in the testis. Am J Physiol Endocrinol Metab. 2007;292(2):E513–22. pmid:17003240.
  62. 62. Alfano M, Tascini AS, Pederzoli F, Locatelli I, Nebuloni M, Giannese F, et al. Aging, inflammation and DNA damage in the somatic testicular niche with idiopathic germ cell aplasia. Nat Commun. 2021;12(1):5205. Epub 2021/09/03. pmid:34471128; PubMed Central PMCID: PMC8410861.
  63. 63. Angulo MA, Butler MG, Cataletto ME. Prader-Willi syndrome: a review of clinical, genetic, and endocrine findings. J Endocrinol Invest. 2015;38(12):1249–63. pmid:26062517; PubMed Central PMCID: PMC4630255.
  64. 64. Goldstone AP, Holland AJ, Hauffa BP, Hokken-Koelega AC, Tauber M, speakers contributors at the Second Expert Meeting of the Comprehensive Care of Patients with PWS. Recommendations for the diagnosis and management of Prader-Willi syndrome. J Clin Endocrinol Metab. 2008;93(11):4183–97. pmid:18697869.
  65. 65. Davies W, Lynn PM, Relkovic D, Wilkinson LS. Imprinted genes and neuroendocrine function. Front Neuroendocrinol. 2008;29(3):413–27. pmid:18206218.
  66. 66. Ivanova E, Kelsey G. Imprinted genes and hypothalamic function. J Mol Endocrinol. 2011;47(2):R67–74. pmid:21798993.
  67. 67. Itier JM, Tremp GL, Leonard JF, Multon MC, Ret G, Schweighoffer F, et al. Imprinted gene in postnatal growth role. Nature. 1998;393(6681):125–6. pmid:9603515.
  68. 68. Cheung LY, Rizzoti K, Lovell-Badge R, Le Tissier PR. Pituitary phenotypes of mice lacking the notch signalling ligand delta-like 1 homologue. J Neuroendocrinol. 2013;25(4):391–401. pmid:23279263; PubMed Central PMCID: PMC3664429.
  69. 69. Charalambous M, da Rocha ST, Hernandez A, Ferguson-Smith AC. Perturbations to the IGF1 growth pathway and adult energy homeostasis following disruption of mouse chromosome 12 imprinting. Acta Physiol (Oxf). 2014;210(1):174–87. pmid:24034272; PubMed Central PMCID: PMC3992899.
  70. 70. Hassan HA, Enright WJ, Tucker HA, Merkel RA. Estrogen and androgen elicit growth hormone release via dissimilar patterns of hypothalamic neuropeptide secretion. Steroids. 2001;66(2):71–80. pmid:11146085.
  71. 71. Chowen JA, Frago LM, Argente J. The regulation of GH secretion by sex steroids. Eur J Endocrinol. 2004;151 Suppl 3:U95–100. pmid:15554893.
  72. 72. D’Ercole AJ, Stiles AD, Underwood LE. Tissue concentrations of somatomedin C: further evidence for multiple sites of synthesis and paracrine or autocrine mechanisms of action. Proc Natl Acad Sci U S A. 1984;81(3):935–9. pmid:6583688; PubMed Central PMCID: PMC344954.
  73. 73. Mathews LS, Norstedt G, Palmiter RD. Regulation of insulin-like growth factor I gene expression by growth hormone. Proc Natl Acad Sci U S A. 1986;83(24):9343–7. pmid:3467309; PubMed Central PMCID: PMC387134.
  74. 74. Hewitt SC, Li Y, Li L, Korach KS. Estrogen-mediated regulation of Igf1 transcription and uterine growth involves direct binding of estrogen receptor alpha to estrogen-responsive elements. J Biol Chem. 2010;285(4):2676–85. pmid:19920132; PubMed Central PMCID: PMC2807324.
  75. 75. Mauras N, Rogol AD, Haymond MW, Veldhuis JD. Sex steroids, growth hormone, insulin-like growth factor-1: neuroendocrine and metabolic regulation in puberty. Horm Res. 1996;45(1–2):74–80. pmid:8742123.
  76. 76. Udy GB, Towers RP, Snell RG, Wilkins RJ, Park SH, Ram PA, et al. Requirement of STAT5b for sexual dimorphism of body growth rates and liver gene expression. Proc Natl Acad Sci U S A. 1997;94(14):7239–44. pmid:9207075; PubMed Central PMCID: PMC23803.
  77. 77. Lupu F, Terwilliger JD, Lee K, Segre GV, Efstratiadis A. Roles of growth hormone and insulin-like growth factor 1 in mouse postnatal growth. Dev Biol. 2001;229(1):141–62. pmid:11133160.
  78. 78. Giagulli VA, Castellana M, Pelusi C, Triggiani V. Androgens, Body Composition, and Their Metabolism Based on Sex. Frontiers of hormone research. 2019;53:18–32. Epub 2019/09/10. pmid:31499499.
  79. 79. Rao PM, Kelly DM, Jones TH. Testosterone and insulin resistance in the metabolic syndrome and T2DM in men. Nat Rev Endocrinol. 2013;9(8):479–93. pmid:23797822.
  80. 80. Dimopoulou C, Goulis DG, Corona G, Maggi M. The complex association between metabolic syndrome and male hypogonadism. Metabolism. 2018;86:61–8. pmid:29656047.
  81. 81. Sato T, Matsumoto T, Yamada T, Watanabe T, Kawano H, Kato S. Late onset of obesity in male androgen receptor-deficient (AR KO) mice. Biochem Biophys Res Commun. 2003;300(1):167–71. pmid:12480537.
  82. 82. Fan W, Yanase T, Nomura M, Okabe T, Goto K, Sato T, et al. Androgen receptor null male mice develop late-onset obesity caused by decreased energy expenditure and lipolytic activity but show normal insulin sensitivity with high adiponectin secretion. Diabetes. 2005;54(4):1000–8. pmid:15793238.
  83. 83. Lin HY, Xu Q, Yeh S, Wang RS, Sparks JD, Chang C. Insulin and leptin resistance with hyperleptinemia in mice lacking androgen receptor. Diabetes. 2005;54(6):1717–25. pmid:15919793.
  84. 84. Tena-Sempere M, Pinilla L, Gonzalez LC, Dieguez C, Casanueva FF, Aguilar E. Leptin inhibits testosterone secretion from adult rat testis in vitro. J Endocrinol. 1999;161(2):211–8. pmid:10320818.
  85. 85. Martins FF, Aguila MB, Mandarim-de-Lacerda CA. Impaired steroidogenesis in the testis of leptin-deficient mice (ob/ob -/-). Acta Histochem. 2017;119(5):508–15. pmid:28506466.
  86. 86. Wang X, Rundle CH, Wergedal JE, Srivastava AK, Mohan S, Lau KH. Loss of sex-specific difference in femoral bone parameters in male leptin knockout mice. Calcif Tissue Int. 2007;80(6):374–82. pmid:17516018.
  87. 87. Tunster SJ, Boqué-Sastre R, McNamara GI, Hunter SM, Creeth HDJ, John RM. Peg3 Deficiency Results in Sexually Dimorphic Losses and Gains in the Normal Repertoire of Placental Hormones. Frontiers in cell and developmental biology. 2018;6:123. Epub 2018/10/16. pmid:30320110; PubMed Central PMCID: PMC6170603.
  88. 88. Forger NG, Strahan JA, Castillo-Ruiz A. Cellular and molecular mechanisms of sexual differentiation in the mammalian nervous system. Front Neuroendocrinol. 2016;40:67–86. pmid:26790970; PubMed Central PMCID: PMC4897775.
  89. 89. Ogawa S, Lubahn DB, Korach KS, Pfaff DW. Behavioral effects of estrogen receptor gene disruption in male mice. Proc Natl Acad Sci U S A. 1997;94(4):1476–81. pmid:9037078; PubMed Central PMCID: PMC19816.
  90. 90. Gospocic J. Role of epicardium-derived Pw1 cells in heart homeostasis and growth during pregnancy. PhD thesis. 2013.
  91. 91. Kim JH, Han GC, Seo JY, Park I, Park W, Jeong HW, et al. Sex hormones establish a reserve pool of adult muscle stem cells. Nat Cell Biol. 2016;18(9):930–40. pmid:27548913.
  92. 92. Yuchi Y, Cai Y, Legein B, De Groef S, Leuckx G, Coppens V, et al. Estrogen Receptor alpha Regulates beta-Cell Formation During Pancreas Development and Following Injury. Diabetes. 2015;64(9):3218–28. pmid:26015547.
  93. 93. Pullar J, Wickramasinghe K, Demaio AR, Roberts N, Perez-Blanco KM, Noonan K, et al. The impact of maternal nutrition on offspring’s risk of non-communicable diseases in adulthood: a systematic review. Journal of global health. 2019;9(2):020405. Epub 2019/10/28. pmid:31656604; PubMed Central PMCID: PMC6790233 members of the World Health Organisation.
  94. 94. Sales VM, Ferguson-Smith AC, Patti ME. Epigenetic Mechanisms of Transmission of Metabolic Disease across Generations. Cell Metab. 2017;25(3):559–71. Epub 2017/03/09. pmid:28273478; PubMed Central PMCID: PMC5404272.
  95. 95. Mandarim-de-Lacerda CA, Fernandes-Santos C, Aguila MB. Image analysis and quantitative morphology. Methods in molecular biology (Clifton, NJ). 2010;611:211–25. Epub 2009/12/05. pmid:19960334.