Escalated aggressive behavior with pathological signatures poses immense risk for violent and antisocial behavior and is a major challenge to human welfare (Miczek et al., 2007). Aggression is a behavioral adaptation to threat and competition, but when it becomes excessive, uncontrollable and out of contextd it is considered as maladaptive and pathological (Neves and Tudela, 2015; Waltes et al., 2016).

Animal aggression is considered escalated and pathological if it displays short attack latency, prolonged attack duration, attacks targeted on inappropriate partners, and body parts prone to serious injury; attacks not signaled by threats; or ignorance of signals of opponents (Bacq et al., 2020). In general, animal models of escalated aggression can largely explain abnormal aggressive behavior in humans (Miczek et al., 2015). Therefore, it is extremely important to understand the biological factors contributing to shift of normal adaptive aggression to escalated and pathological form.

Mounting epidemiological evidences link early life stressful experiences with deteriorating mental health (Nelson and Trainor, 2007; Duke et al., 2010; Haller et al., 2014; Hunt et al., 2019; Mitchell and Aamodt, 2005). In particular, stress exposures around puberty including childhood and adolescence including fear, maltreatment, physical and sexual abuse confers susceptibility to aggressive behavioral disorders (Veenema, 2009; Tzanoulinou and Sandi, 2017; Bounoua et al., 2020). Although, early life adversities are considered as one of the potential triggers for abnormal aggression, biological insights are obscure. More importantly, majority of research in the field of aggressive biology have focused on the adaptive form without really considering the escalated and pathological forms (de Boer, 2018). Márquez et al., 2013 developed a novel animal model which showed the effect of peripubertal fearful exposures on male pathological aggression at adulthood in Wistar Han rats. They primarily focused on neural circuits of aggression and on a single gene Maoa in isolation.

Considering aggressive behavior as a multidimensional trait, we rationalized that unbiased genome wide investigation would decipher key molecular pathways that can be exploited further as prediction and intervention targets. We modeled peripubertal stress (PPS) induced escalated aggression in laboratory bred Balb/c mice and screened the male cohort showing extremes of behavioral changes. Female mice showed resilience towards PPS induced escalated aggressive behavior as also reported previously (Konar et al., 2019).

Next, we performed a sex-specific transcriptome analysis in vulnerable brain regions of hypothalamus (Hypo) and prefrontal cortex (PFC). Hypothalamus is a key brain region for expression of aggressive behavior and neural circuit-specific manipulation experiments revealed that ventromedial hypothalamus is the crucial for inter-male aggression (Lin et al., 2011; Falkner et al., 2016). While Hypo is considered as the trigger center for aggression, PFC plays opposite regulatory role being involved in inhibition of threat provoked aggressive behavior. More importantly, direct neuronal projections from PFC to Hypo have been suggested to control both type and amplitude of aggressive behavior (Choy et al., 2018; Biro et al., 2018). Therefore, we primarily focused on hypothalamic molecular culprits of escalated aggression and also included PFC in our study to understand inter-brain regional molecular regulation if any.

We prioritized Ttr gene given its (i) top rank in Hypo transcriptome analysis and unique sex specific diametrically opposite expression pattern in Hypo and PFC and (iii) long-term gene expression changes from early peripubertal age till adulthood. Brain-region-specific changes in Ttr gene expression and promoter methylation and thyroid hormone (TH) availability, one of the key functions of TTR were evident in PPS-induced escalated aggressive male mice. Further, targeted gene manipulation revealed causal role for hypothalamic Ttr in development of escalated aggressive behavior. However, causal relationship of Ttr regulated hypothalamic TH availability and escalated aggression is still to be explored.

In the present study, we fill in the existent gap of knowledge about the molecular roots of escalated aggressive behavior. Based on unbiased transcriptome screening, we identify novel role of TTR in long-term programming of escalated aggression induced by PPS. Our findings also indicate the possible involvement of TTR dependent brain TH availability in such abnormal behavioral response. However, further studies are necessitated to identify the definite molecular mechanism.

TTR is a 55 kDa protein that is synthesized in choroid plexus epithelial cells of the brain (Alshehri et al., 2015) until recently it was identified in neurons (Li et al., 2011; Zawiślak et al., 2017) and astrocytes indicating wide expression of the protein in CNS. Interestingly, we observed TTR immunoreactivity in different cell types in Hypo and PFC. TTR expression in PFC was similar to brain endothelial cell marker proteins (Tang et al., 2017). Our cell type specificity analysis using publicly available databases also revealed highest expression Ttr mRNA in endothelial cells followed by neurons, microglia and astrocytes in mouse brain. Additional experiments mainly colocalization with specific markers are required to confirm the cell types.

We observed long term impact of PPS on Ttr expression in brain-region and sex-specific diametrically opposite manner. Such a unique pattern of the only top ranking gene from our transcriptome data intrigued us to explore the functional outcomes. TTR is involved in the uptake of T4 from blood to CSF and local distribution of TH in brain (Alshehri et al., 2015). TTR has also been assigned other functions including proteolysis of Neuropeptide Y (Nunes et al., 2006), neuroprotection and regeneration of damaged neurons (Sousa et al., 2004). We focused on circulating and brain TH levels based on multiple reasons. Our transcriptome data showed significant expression change in TH signaling genes primarily Trh and Nrgn in Hypo. Several other genes in our list (PFC-Rasd2, Fosl2, Nr4a3, Inf2, Arl4d, Dcn, Pcp4l1, Drd2, Syndig1l, Hspa1a, Spock3; Hypo- Oxt, Cdhr1, Col23a1, Dgkk, Dkk3, Cck, Ptpro) showed overlap with literature available on T3 responsive genes in primary cultured neurons (Gil-Ibáñez et al., 2014; Richard and Flamant, 2018). Also, TH action in developing brain have been considered as critical determinants of multiple neurological deficits (Préau et al., 2015).

In clinical settings, TH abnormalities are diagnosed by serum parameters. However, we show that alterations in Ttr gene expression paralleled perturbation in T4 and T3 content in brain tissues of Hypo and PFC without affecting the circulating levels of the hormone. Determining the effective concentration of T4 and T3 in brain tissues is difficult owing to multiple factors driving their synthesis, transport across blood-brain and blood-CSF barrier, intracellular distribution and activation/inactivation (Schroeder and Privalsky, 2014; Bárez-López and Guadaño-Ferraz, 2017). Therefore, future investigations on TH transporters, DIO2 and DIO3 enzymes are necessitated to determine the local TH availability in brain (Mayerl et al., 2014; Williams and Bassett, 2011).

Any change in brain TH state has a direct influence on TH responsive genes. TH deficiency during postnatal brain development causes irreversible neurological manifestations through target gene expression changes (Vallortigara et al., 2008). Nrgn is one such brain specific TH responsive gene that was also amongst the top ranking DEGs in our transcriptome data containing TRE elements in promoter and its transcription is dependent on TH in brain (Pak et al., 2000; Husson et al., 2004). Nrgn regulates synaptic plasticity by activating calmodulin kinase II (CaMKII) protein and spine density. We observed significant reduction in Nrgn transcript and NRGN protein levels in concordance with decrease in T4 and T3 levels in Hypo.

Until now we found a strong association between altered brain TTR expression and function with escalated aggressive behavior but causal relationship was yet to be established. Next, we showed that intra-hypothalamic Ttr gene knockdown evoked escalated aggression in unstressed males to a similar extent to that of PPS induced males. Ttr gene knockdown also led to decrease in Hypo T3 content and alteration in expression TH signaling genes, Nrgn and Trh. Ttr gene silencing also reduced hairless (Hr), a universal TH responsive gene that is studied to monitor the local TH status in brain (Herwig et al., 2014).

Escalated aggressive behavior, reduced Ttr mRNA expression and T3 content in Hypo was also evident in F1 male progenies of SEagg males. Previous studies suggest that TH changes in neonatal brain can elicit neuroendocrine abnormalities in their F1 progenies. Also, developmental exposure of thyroxine disrupting chemicals can affect gene expression and behavior in later generations (Morte et al., 2018). These data indicated that behavioral and molecular consequences are dependent on TTR. However, it is not clear whether the decrease of TTR expression causes the escalated aggression via the regulation of hypothalamic TH availability or the effect of altered TTR expression on the aggression and TH availability are independent.

Ttr promoter showed altered methylation pattern in Hypo and PFC of SEagg males. As DNA methylation is considered important for lasting influence of stressful experiences (Vukojevic et al., 2014; Gulmez Karaca et al., 2020), it might underlie long-term impact of PPS on Ttr expression and aggressive behavior. Further studies are required to determine the molecular pathways underlying involvement of TTR in aggressive behavior.

In conclusion, we delineate novel role of TTR in manifestation of early life stress induced escalated aggression. TTR signaling in brain can also serve as a valid molecular predictor as well as intervention target in excessive pathological aggression. Our findings have inherent limitations of investigations in animal models and therefore further studies are warranted in relevant human cohort to establish role of TTR in abnormal aggression and related psychopathologies. Our work also provides resource for investigating sexual dimorphism in behavioral disorders and deciphering susceptibility as well as protective pathways.

RNA sequencing data have been deposited in GEO under accession code GSE199844. All data generated or analyzed during this study are included in the manuscript and supplementary files. Source data files have been provided for Fig. 1C, Fig. 3D, Fig. 4B-4J, Fig. 4P, Fig.5B-5D, Fig.6B-6E, Fig.7 and Fig. 4-figure supplement 1.

1. 1. Konar A

   (2022) NCBI Gene Expression Omnibus

   ID GSE199844. Early life stressful experiences escalate aggressive behavior in adulthood via changes in transthyretin expression and function.

   http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE199844

1. 1. Alshehri B
   2. D’Souza DG
   3. Lee JY
   4. Petratos S
   5. Richardson SJ

   (2015) The diversity of mechanisms influenced by transthyretin in neurobiology: development, disease and endocrine disruption

   Journal of Neuroendocrinology 27:303–323.

   https://doi.org/10.1111/jne.12271

   • PubMed
   • Google Scholar
2. 1. Bacq A
   2. Astori S
   3. Gebara E
   4. Tang W
   5. Silva BA
   6. Sanchez-Mut J
   7. Grosse J
   8. Zanoletti O
   9. Maclachlan C
   10. Knott GW
   11. Gräff J
   12. Sandi C

   (2020) Amygdala glun2b-NMDAR dysfunction is critical in abnormal aggression of neurodevelopmental origin induced by st8sia2 deficiency

   Mol Psychiatry 25:2144–2161.

   https://doi.org/10.1038/s41380-018-0132-3

   • Google Scholar
3. 1. Bárez-López S
   2. Guadaño-Ferraz A

   (2017) Thyroid hormone availability and action during brain development in rodents

   Frontiers in Cellular Neuroscience 11:240.

   https://doi.org/10.3389/fncel.2017.00240

   • PubMed
   • Google Scholar
4. 1. Biro L
   2. Sipos E
   3. Bruzsik B
   4. Farkas I
   5. Zelena D
   6. Balazsfi D
   7. Toth M
   8. Haller J

   (2018) Task division within the prefrontal cortex: distinct neuron populations selectively control different aspects of aggressive behavior via the hypothalamus

   The Journal of Neuroscience 38:4065–4075.

   https://doi.org/10.1523/JNEUROSCI.3234-17.2018

   • PubMed
   • Google Scholar
5. 1. Bounoua N
   2. Miglin R
   3. Spielberg JM
   4. Sadeh N

   (2020) Childhood assaultive trauma and physical aggression: links with cortical thickness in prefrontal and occipital cortices

   NeuroImage. Clinical 27:102321.

   https://doi.org/10.1016/j.nicl.2020.102321

   • PubMed
   • Google Scholar
6. 1. Choy O
   2. Raine A
   3. Hamilton RH

   (2018) Stimulation of the prefrontal cortex reduces intentions to commit aggression: a randomized, double-blind, placebo-controlled, stratified, parallel-group trial

   The Journal of Neuroscience 38:6505–6512.

   https://doi.org/10.1523/JNEUROSCI.3317-17.2018

   • PubMed
   • Google Scholar
7. 1. de Boer SF

   (2018) Animal models of excessive aggression: implications for human aggression and violence

   Current Opinion in Psychology 19:81–87.

   https://doi.org/10.1016/j.copsyc.2017.04.006

   • PubMed
   • Google Scholar
8. 1. Duke NN
   2. Pettingell SL
   3. McMorris BJ
   4. Borowsky IW

   (2010) Adolescent violence perpetration: associations with multiple types of adverse childhood experiences

   Pediatrics 125:e778–e786.

   https://doi.org/10.1542/peds.2009-0597

   • PubMed
   • Google Scholar
9. 1. Falkner AL
   2. Grosenick L
   3. Davidson TJ
   4. Deisseroth K
   5. Lin D

   (2016) Hypothalamic control of male aggression-seeking behavior

   Nature Neuroscience 19:596–604.

   https://doi.org/10.1038/nn.4264

   • PubMed
   • Google Scholar
10. 1. Gil-Ibáñez P
    2. Bernal J
    3. Morte B

    (2014) Thyroid hormone regulation of gene expression in primary cerebrocortical cells: role of thyroid hormone receptor subtypes and interactions with retinoic acid and glucocorticoids

    PLOS ONE 9:e91692.

    https://doi.org/10.1371/journal.pone.0091692

    • PubMed
    • Google Scholar
11. 1. Gulmez Karaca K
    2. Kupke J
    3. Brito DVC
    4. Zeuch B
    5. Thome C
    6. Weichenhan D
    7. Lutsik P
    8. Plass C
    9. Oliveira AMM

    (2020) Neuronal ensemble-specific DNA methylation strengthens engram stability

    Nature Communications 11:639.

    https://doi.org/10.1038/s41467-020-14498-4

    • PubMed
    • Google Scholar
12. 1. Haller J
    2. Harold G
    3. Sandi C
    4. Neumann ID

    (2014) Effects of adverse early-life events on aggression and anti-social behaviours in animals and humans

    Journal of Neuroendocrinology 26:724–738.

    https://doi.org/10.1111/jne.12182

    • PubMed
    • Google Scholar
13. 1. Herwig A
    2. Campbell G
    3. Mayer C-D
    4. Boelen A
    5. Anderson RA
    6. Ross AW
    7. Mercer JG
    8. Barrett P

    (2014) A thyroid hormone challenge in hypothyroid rats identifies T3 regulated genes in the hypothalamus and in models with altered energy balance and glucose homeostasis

    Thyroid 24:1575–1593.

    https://doi.org/10.1089/thy.2014.0169

    • PubMed
    • Google Scholar
14. 1. Hu Z
    2. Dai Y
    3. Prentki M
    4. Chohnan S
    5. Lane MD

    (2005) A role for hypothalamic malonyl-CoA in the control of food intake

    The Journal of Biological Chemistry 280:39681–39683.

    https://doi.org/10.1074/jbc.C500398200

    • PubMed
    • Google Scholar
15. Book

    1. Hunt KS
    2. Iaconetti MJ
    3. Maass KT

    (2019)

    Recidivism among Federal Violent Offenders

    United States Sentencing Commission.

    • Google Scholar
16. 1. Husson M
    2. Enderlin V
    3. Alfos S
    4. Boucheron C
    5. Pallet V
    6. Higueret P

    (2004) Expression of neurogranin and neuromodulin is affected in the striatum of vitamin A-deprived rats

    Brain Research. Molecular Brain Research 123:7–17.

    https://doi.org/10.1016/j.molbrainres.2003.12.012

    • PubMed
    • Google Scholar
17. 1. Konar A
    2. Rastogi M
    3. Bhambri A

    (2019) Brain region specific methylation and SIRT1 binding changes in MAOA promoter is associated with sexual dimorphism in early life stress induced aggressive behavior

    Neurochemistry International 129:104510.

    https://doi.org/10.1016/j.neuint.2019.104510

    • PubMed
    • Google Scholar
18. 1. Koolhaas JM
    2. Coppens CM
    3. de Boer SF
    4. Buwalda B
    5. Meerlo P
    6. Timmermans PJA

    (2013) The resident-intruder paradigm: a standardized test for aggression, violence and social stress

    Journal of Visualized Experiments 4:e4367.

    https://doi.org/10.3791/4367

    • PubMed
    • Google Scholar
19. 1. Li X
    2. Masliah E
    3. Reixach N
    4. Buxbaum JN

    (2011) Neuronal production of transthyretin in human and murine Alzheimer’s disease: is it protective?

    The Journal of Neuroscience 31:12483–12490.

    https://doi.org/10.1523/JNEUROSCI.2417-11.2011

    • PubMed
    • Google Scholar
20. 1. Lin D
    2. Boyle MP
    3. Dollar P
    4. Lee H
    5. Lein ES
    6. Perona P
    7. Anderson DJ

    (2011) Functional identification of an aggression locus in the mouse hypothalamus

    Nature 470:221–226.

    https://doi.org/10.1038/nature09736

    • PubMed
    • Google Scholar
21. 1. Márquez C
    2. Poirier GL
    3. Cordero MI
    4. Larsen MH
    5. Groner A
    6. Marquis J
    7. Magistretti PJ
    8. Trono D
    9. Sandi C

    (2013) Peripuberty stress leads to abnormal aggression, altered amygdala and orbitofrontal reactivity and increased prefrontal MAOA gene expression

    Translational Psychiatry 3:e216.

    https://doi.org/10.1038/tp.2012.144

    • PubMed
    • Google Scholar
22. 1. Mayerl S
    2. Müller J
    3. Bauer R
    4. Richert S
    5. Kassmann CM
    6. Darras VM
    7. Buder K
    8. Boelen A
    9. Visser TJ
    10. Heuer H

    (2014) Transporters MCT8 and OATP1C1 maintain murine brain thyroid hormone homeostasis

    The Journal of Clinical Investigation 124:1987–1999.

    https://doi.org/10.1172/JCI70324

    • PubMed
    • Google Scholar
23. 1. Miczek KA
    2. de Almeida RMM
    3. Kravitz EA
    4. Rissman EF
    5. de Boer SF
    6. Raine A

    (2007) Neurobiology of escalated aggression and violence

    The Journal of Neuroscience 27:11803–11806.

    https://doi.org/10.1523/JNEUROSCI.3500-07.2007

    • PubMed
    • Google Scholar
24. 1. Miczek KA
    2. Takahashi A
    3. Gobrogge KL
    4. Hwa LS

    (2015) Escalated aggression in animal models: shedding new light on mesocorticolimbic circuits

    Curr Opin Behav Sci 3:90–95.

    https://doi.org/10.1016/j.cobeha.2015.02.007

    • Google Scholar
25. 1. Mitchell H
    2. Aamodt MG

    (2005) The incidence of child abuse in serial killers

    Journal of Police and Criminal Psychology 20:40–47.

    https://doi.org/10.1007/BF02806705

    • Google Scholar
26. 1. Mohn F
    2. Weber M
    3. Schübeler D
    4. Roloff TC

    (2009) Methylated DNA immunoprecipitation (medip)

    Methods in Molecular Biology 507:55–64.

    https://doi.org/10.1007/978-1-59745-522-0_5

    • PubMed
    • Google Scholar
27. 1. Morte B
    2. Gil-Ibáñez P
    3. Bernal J

    (2018) Regulation of gene expression by thyroid hormone in primary astrocytes: factors influencing the genomic response

    Endocrinology 159:2083–2092.

    https://doi.org/10.1210/en.2017-03084

    • PubMed
    • Google Scholar
28. 1. Mukhopadhyay A
    2. Deplancke B
    3. Walhout AJM
    4. Tissenbaum HA

    (2008) Chromatin immunoprecipitation (CHIP) coupled to detection by quantitative real-time PCR to study transcription factor binding to DNA in Caenorhabditis elegans

    Nature Protocols 3:698–709.

    https://doi.org/10.1038/nprot.2008.38

    • PubMed
    • Google Scholar
29. 1. Nelson RJ
    2. Trainor BC

    (2007) Neural mechanisms of aggression

    Nature Reviews. Neuroscience 8:536–546.

    https://doi.org/10.1038/nrn2174

    • PubMed
    • Google Scholar
30. 1. Neves S
    2. Tudela J

    (2015) Aggression and violence – “ a journey through the human brain.”

    European Psychiatry 30:1243.

    https://doi.org/10.1016/S0924-9338(15)30976-7

    • Google Scholar
31. 1. Nunes AF
    2. Saraiva MJ
    3. Sousa MM

    (2006) Transthyretin knockouts are a new mouse model for increased neuropeptide Y

    FASEB Journal 20:166–168.

    https://doi.org/10.1096/fj.05-4106fje

    • PubMed
    • Google Scholar
32. 1. Pak JH
    2. Huang FL
    3. Li J
    4. Balschun D
    5. Reymann KG
    6. Chiang C
    7. Westphal H
    8. Huang KP

    (2000) Involvement of neurogranin in the modulation of calcium/calmodulin-dependent protein kinase II, synaptic plasticity, and spatial learning: a study with knockout mice

    PNAS 97:11232–11237.

    https://doi.org/10.1073/pnas.210184697

    • PubMed
    • Google Scholar
33. 1. Préau L
    2. Fini JB
    3. Morvan-Dubois G
    4. Demeneix B

    (2015) Thyroid hormone signaling during early neurogenesis and its significance as a vulnerable window for endocrine disruption

    Biochimica et Biophysica Acta 1849:112–121.

    https://doi.org/10.1016/j.bbagrm.2014.06.015

    • PubMed
    • Google Scholar
34. 1. Richard S
    2. Flamant F

    (2018) Regulation of T3 availability in the developing brain: the mouse genetics contribution

    Frontiers in Endocrinology 9:265.

    https://doi.org/10.3389/fendo.2018.00265

    • PubMed
    • Google Scholar
35. 1. Schroeder AC
    2. Privalsky ML

    (2014) Thyroid hormones, T3 and T4, in the brain

    Frontiers in Endocrinology 5:40.

    https://doi.org/10.3389/fendo.2014.00040

    • PubMed
    • Google Scholar
36. 1. Sousa JC
    2. Grandela C
    3. Fernández-Ruiz J
    4. de Miguel R
    5. de Sousa L
    6. Magalhães AI
    7. Saraiva MJ
    8. Sousa N
    9. Palha JA

    (2004) Transthyretin is involved in depression-like behaviour and exploratory activity

    Journal of Neurochemistry 88:1052–1058.

    https://doi.org/10.1046/j.1471-4159.2003.02309.x

    • PubMed
    • Google Scholar
37. 1. Takahashi A
    2. Miczek KA

    (2014) Neurogenetics of aggressive behavior: studies in rodents

    Current Topics in Behavioral Neurosciences 17:3–44.

    https://doi.org/10.1007/7854_2013_263

    • PubMed
    • Google Scholar
38. 1. Tang M
    2. Gao G
    3. Rueda CB
    4. Yu H
    5. Thibodeaux DN
    6. Awano T
    7. Engelstad KM
    8. Sanchez-Quintero M-J
    9. Yang H
    10. Li F
    11. Li H
    12. Su Q
    13. Shetler KE
    14. Jones L
    15. Seo R
    16. McConathy J
    17. Hillman EM
    18. Noebels JL
    19. De Vivo DC
    20. Monani UR

    (2017) Brain microvasculature defects and GLUT1 deficiency syndrome averted by early repletion of the glucose transporter-1 protein

    Nature Communications 8:14152.

    https://doi.org/10.1038/ncomms14152

    • PubMed
    • Google Scholar
39. 1. Tzanoulinou S
    2. Sandi C

    (2017) The programming of the social brain by stress during childhood and adolescence: from rodents to humans

    Current Topics in Behavioral Neurosciences 30:411–429.

    https://doi.org/10.1007/7854_2015_430

    • PubMed
    • Google Scholar
40. 1. Vallortigara J
    2. Alfos S
    3. Micheau J
    4. Higueret P
    5. Enderlin V

    (2008) T3 administration in adult hypothyroid mice modulates expression of proteins involved in striatal synaptic plasticity and improves motor behavior

    Neurobiology of Disease 31:378–385.

    https://doi.org/10.1016/j.nbd.2008.05.015

    • PubMed
    • Google Scholar
41. 1. Veenema AH

    (2009) Early life stress, the development of aggression and neuroendocrine and neurobiological correlates: what can we learn from animal models?

    Frontiers in Neuroendocrinology 30:497–518.

    https://doi.org/10.1016/j.yfrne.2009.03.003

    • PubMed
    • Google Scholar
42. 1. Vukojevic V
    2. Kolassa I-T
    3. Fastenrath M
    4. Gschwind L
    5. Spalek K
    6. Milnik A
    7. Heck A
    8. Vogler C
    9. Wilker S
    10. Demougin P
    11. Peter F
    12. Atucha E
    13. Stetak A
    14. Roozendaal B
    15. Elbert T
    16. Papassotiropoulos A
    17. de Quervain DJ-F

    (2014) Epigenetic modification of the glucocorticoid receptor gene is linked to traumatic memory and post-traumatic stress disorder risk in genocide survivors

    The Journal of Neuroscience 34:10274–10284.

    https://doi.org/10.1523/JNEUROSCI.1526-14.2014

    • PubMed
    • Google Scholar
43. 1. Waltes R
    2. Chiocchetti AG
    3. Freitag CM

    (2016) The neurobiological basis of human aggression: A review on genetic and epigenetic mechanisms

    American Journal of Medical Genetics. Part B, Neuropsychiatric Genetics 171:650–675.

    https://doi.org/10.1002/ajmg.b.32388

    • PubMed
    • Google Scholar
44. 1. Williams GR
    2. Bassett JHD

    (2011) Deiodinases: the balance of thyroid hormone: local control of thyroid hormone action: role of type 2 deiodinase

    The Journal of Endocrinology 209:261–272.

    https://doi.org/10.1530/JOE-10-0448

    • PubMed
    • Google Scholar
45. 1. Zawiślak A
    2. Jakimowicz P
    3. McCubrey JA
    4. Rakus D

    (2017) Neuron-Derived transthyretin modulates astrocytic glycolysis in hormone-independent manner

    Oncotarget 8:106625–106638.

    https://doi.org/10.18632/oncotarget.22542

    • PubMed
    • Google Scholar

Decision letter after peer review:

Thank you for submitting your article "Early life trauma leads to escalated aggressive behavior and its inheritance by impairing thyroid hormone availability in brain" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, one of whom is a member of our Board of Reviewing Editors, and the evaluation has been overseen by Catherine Dulac as the Senior Editor. The following individual involved in review of your submission has agreed to reveal their identity: Petr N Menshanov (Reviewer #2).

The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.

Essential revisions:

1) The role of the decreased hypothalamic T3 availability in the development of aggression is not causally established. In the current form of the paper, it is not clear whether the decrease of TTR expression causes the aggression via the regulation of hypothalamic T3 availability or the effect of altered TTR expression on the aggression and T3 availability are independent. Therefore, the authors should either focus on the role of TTR or they should provide additional evidence that lack of hypothalamic T3 is the cause of the development of aggression.

2) The studies on epigenetic mediated inheritance of the aggression phenotype are not well supported and should be removed. The data on methylation of the TTR promoter are an important contribution and should remain but should not be considered as a diagnostic for aggressive behavior or the basis for trans-generational inheritance.

3) Statistical analyses should be redone with factorial analyses to make the proper comparisons between sexes, stress condition and brain region.

4) RNA-Seq data should include greater in depth analyses.

5) The images of immunohistochemistry staining are not of sufficient quality and should be replaced.

Reviewer #1 (Recommendations for the authors):

Deiodinase 3 (Dio3) plays critical role in the regulation of hypothalamic T3 availability, therefore, the D3 expression should also be determined along the other TH regulated genes.

According to my knowledge, there is no reliable antibody against dio2. This fact highly questions the validity of the dio2 ELISA data. The only accepted method for the determination of DIO2 protein level is the DIO2 enzyme assay (PMID: 24001133).

The quality of the images illustrating immunocytochemistry is very weak. Better images would be necessary. In addition, it should be described in more details where the TTR-immunoreactive cells were observed in the hypothalamus. In addition, colocalization study should be performed at least with neuronal, glial and endothelial markers to determine what kind of cells express TTR in the hypothalamus.

Reviewer #2 (Recommendations for the authors):

Suggestions for improvement

1.1. In the present form submitted to GEO database, the transcriptomic data on the levels of individual transcripts are reported in aggregate per each experimental group, with no possibility to estimate variance between individual biological replicates. Reporting raw transcriptomic data on the transcript levels in each biological replicate will be beneficial for the readers interested in reanalyses of the data.

It is indispensable to report key elements of the transcriptomic analysis in full. Please update the GEO data with a single file on all levels of individual transcripts in each individual biological replicate.

1.2. The personal experience of the reviewer evidences that old Cufflink-Cuffdiff Tuxedo bioinformatic pipeline for RNA-Seqs (done after Tophat or STAR aligners), if applied properly, is highly sensitive to individual mRNA levels with minor shares. Nevertheless, in the present form of the description of bioinformatics pipeline in the "Methods" section, with missed basic options applied by the authors for Tophat, Cufflink, Cuffdiff utilits, it is hard to estimate the validity of the transcriptomic analyses done by the authors. Please consider improving the description of the bioinformatics methods applied.

1.3. Technical figures like volcano plots are important to control the quality of RNA-Seq, but are too uninformative to demonstrate the results of the experiment. Please consider illustrating the differences between PFC and HPT changes in DEGs on a single figure panel by depicting the parallel changes in the transcript levels identified in PFC (for example, x-axis) and HPT (for example, y-axis). At the same time, it is possible to move the volcano plots to the supplementary figures.

1.4. Preliminary analysis done by the reviewer by applying tissue cell deconvolution methods evidenced in favour of possible specific trends in cell numbers in the limbic brain regions studied. In particular, it cannot be excluded that the juvenile stress episodes suppressed microglia numbers in both the PFC and the HPT.

Please consider providing the data on tissue cell composition in RNA-Seq individual samples with respect to the experimental groups.

(For details, see Sutton et al. 2022 https://doi.org/10.1038/s41467-022.28655-4 or other).

Such an analysis might be illustrative that stress itself is necessary but not sufficient to induce changes in the aggressive behaviours in affected male mice demonstrated excessive violence.

1.5. Please consider additional analyse of the data on the individual transcripts levels reported in the RNA-Seq analysis in a way similar to 2-way or 3-way (when appropriate) factorial ANOVA, to identify possible additive and non-additive patterns of changes in the levels of transcripts.

1.6. It will be also interesting to know on cell specificity of DEGs identified in PFC. Please consider checking the prevalence of individual DEGs with the help of Allen brain transcriptomic atlas (cited in Yao et al.,2021) or another one.

The above mentioned atlas can be assessed by the following link (https://portal.brain-map.org/atlases-and-data/rnaseq)

1.7. Please consider redrawing all figure panels depicted as traditional bars with the (Median, IQR, SD) box plots with individual data dots depicted. In particular, it is possible to achieve this at ease with the JASP freeware (https://jasp-stats.org/)

1.8. It is arguable to use GAPDH as a reference for qPCR assay in the experiments with stress exposures, since GAPDH levels might be affected and even programmed by stress experienced by animals.

Please consider providing a rationale on the applying of GAPDH as an internal standard for mRNA levels instead of β-actin or other mRNAs conventional for stress studies.

1.9. For future studies.

In the absence of cross-fostering experimental schedule for F1 experiments, it is hard to delineate the origins of epigenetic changes identified in the F1 descendants, whether these changes were transmitted directly or indirectly, via mother's specific behaviours on the descendants.

Recommendations for improving the writing and presentation

2.1. The present description of animal procedures in the "Methods" section does not provide enough details on the environments, in which the experimental mice were grown up. Please provide all details on animal housing procedures that might be stressful (social isolation events, social crowding events, numbers of animals per cell etc)

2.2. The description of several methods must be improved to clarify details critical for data comprehension.

For example,

Please indicate the details of screening tests done with subjects (females and anesthetized intruder) that were attacked by mice with pathological aggressive behaviours.

The experimental schedules must be reported in a clear way also. In particular, a brief statement must be done on how the distinct populations of "adult control" and "PPS adult male" mice screened by the Authors were originated from. A similar clarification must be done for female mice groups also.

2.3. The initial two paragraphs in the "Introduction" section do not provide the linear story tale on violence, "escalated" violence and a difference of these two concepts of aggression. The logic of this section must be improved.

2.3a. In the introduction, it looks like that the authors consider an aggression trait as a behavioural continuum between "zero-level" aggression to appropriate violence, and then to "escalated" aggressive behaviour. This point of view is arguable since it cannot be excluded that aggression is a multidimensional trait. Please consider revising and clarifying.

2.3b. It is possible to criticize "escalated" aggressive behaviour as unproductive. However, please do not make generalized negative statements on the nature of general violence in the Introduction and in the Discussion. Such statements might blackmail protective types of aggression critical for survival in mammals and humans.

2.4. Please avoid to made generalized statements that "A can lead to B". Such statements with a strong modal verb "can" are highly misleading since the development of behavioural traits is not linear and depends on both genotype and environmental context often. Better to speak that "A might lead to B" under certain circumstances.

Several examples:

2.4a."Escalated aggressive behavior … can lead to antisocial and criminal activities" [P3-S1-LL1-2]. This is misleading for numerous animals and for specific types of aggressive behaviours in H. sapiens. Please consider revising.

2.4b. "…pathological aggression has emerged as a consequence of early life adversities…" [P3-S2.LL9-10]. This statement is misleading, as it blames unrightfully all children affected by harsh life. Please consider revising.

2.4c. "…we inferred that Ttr promoter methylation could serve as a predictor of … behavioral deficits." [P24-S1-LL14-15] – Better to talk about "possible behavioural deficits", not about "behavioural deficits".

The above-mentioned list of examples is not exhaustive.

Reviewer #3 (Recommendations for the authors):

1) Examination of some of the DEGs from the other pathways identified in the KEGG analysis would either strengthen the argument that changes are specific to the TH pathway or highlight that changes are more wide spread.

2) Some discussion of how a change in the amount of T3 and T4 could lead to aberrant aggression would enhance the manuscript.

3) Overall the sex difference, which is profound, is not given much attention. There are more DEGs identified in females subject to peripubertal stress than males, yet there is no change in aggressive behavior, so what does this tell us? Also, does testosterone play a role in the sex difference in both the transcriptome and the behavioral changes? If males were gonadectomized, would the same transcriptional profile be apparent and would the behavior also be there? Or would the two endpoints diverge, belying the noting that there is a casual connect between them.

4) The peri-pubertal stress was conducted on 7 random days from PN28 to PN42. The timing of puberty is different in males and females, being earlier in females. Where measures taken to determine the stage of puberty in each animal (i.e. vaginal opening, preputial separation)? Did the stress impact the timing of puberty?

5) The transgenerational assertions should either be dropped or the study carried out to the F2 generation.

6) How was the use of MeDIP specific to the promoter for Ttr?

7) What are the circulating androgen levels in the males from the various groups? Could the PP have altered the HPGA that then in turns alters behavior?

8) It does not seem appropriate to refer to "donut shaped cells".

9) Figure 4J – appears mislabeled, has Hypo twice and no PFC.

10) Figure 4L-P – why aren't the individual points plotted for the mRNA and protein.

11) Whenever both brain areas are considered the statistics should be 2-way ANOVA with brain region and treatment as factors?

12) The word "trauma" in the context used here connotes an emotional interpretation of stressful or fearful events. We do not know if the mice are experiencing trauma, instead we know they are being subject to fearful and stress-inducing experiences. It is suggested that the word trauma be removed throughout and replaced with more precise terminology.

[Editors' note: further revisions were suggested prior to acceptance, as described below.]

Thank you for resubmitting your work entitled "Pivotal role of TTR in early life stressful experiences induced escalated aggressive behavior" for further consideration by eLife. Your revised article has been evaluated by Catherine Dulac (Senior Editor) and a Reviewing Editor.

The manuscript has been improved but there are some remaining issues that need to be addressed, as outlined in detail below:

The following additional revisions are essential:

1) The conclusions on large sex differences between male and female hypothalamic transcriptomes are not supported by the data due to the problematic sample GSM5988437 (Female Hypo Experimental BiologicalReplicate1). This sample must be withdrawn from the analysis, and the male-female comparisons for hypothalamic transcriptomes must be re-estimated without this problematic sample. This will help to dismiss any spurious claims about the "augmented" male-female hypothalamic differences.

2) Data about the specificity of the TTR immunocytochemistry and D2 elisa would be absolutely necessary. The authors describe the changes of 2 TH transmitters and Dio2. Dio3 is at least an important regulator of TH availability in the brain as Dio2. So either data about Dio3 expression should be added or data about the expression of TH transporters and Dio2 should be removed.

3) Determination of the cell types expressing TTR should be very fast and easy with double-labeling immunocytochemistry and would increase the value of the paper, however, this is not absolutely necessary to support the conclusions of the paper.

Reviewer #1 (Recommendations for the authors):

The manuscript was highly improved, but there are some points that need further change.

The immunocytochemical images of Figure 3 were highly improved. On E, the TTR immunoreactivity seems to be neuronal, while the TTR seems to be present in perivascular localization suggesting endothelial or astrocytic expression. The cell type expressing TTR in the studied brain regions should be determined using double-labeling immunofluorescence.

The specificity of the TTR immunostaining should be proven.

It should be described which hypothalamic region is shown in the images.

The red immunofluorescent signal is visible inside the ventricle on F. This questions the specificity of the staining.

The length of the scale bars is hardly visible. In addition, while the magnification of the left and right panels are obviously different, the scale bars are the same on all images. This should be corrected.

As the balance of DIO2 and DIO3 determines the TH availability, in addition to the expression level of transporters and DIO2, the expression of DIO3 should be also presented.

As there is no specific DIO2 antibody, the validity of the DIO2 protein data is highly questionable. Therefore, data about the specificity of the assay should be provided.

Dio2 is negatively regulated by TH in most parts of the brain, but changes in TH level have no effect on Dio2 activity in the hypothalamus suggesting that the observed changes of Dio2 are not compensatory.

Reviewer #2 (Recommendations for the authors):

1. It is still unclear whether Tophat, Cufflinks, Cuffmerge, and Cuffdiff tools were applied with the default settings or specific settings.

Please indicate clearly the settings applied for Tophat, Cuffdiff, and other bioinformatic tools used, within the "RNA-sequencing" subsection of the "Materials and methods" section.

2. Individual data provided for RNA-seq evidenced that Cuffdiff2 "misbehaved" during bioinformatics analysis while estimating FPKM levels for some abundant transcripts and filtering biological noise. As a result, Cuffdiff2 nullifies FPKM levels for some abundant transcripts like mitochondrial mt-Nd6 non-uniformly, only in specific samples. This issue should not be considered as a large problem, as Cuffdiff2's misbehaviour in such situations is already known.

To deal with this issue, the authors must include in the "Discussion" section a short sentence on this possible limitation of the study arising when comparing levels of highly abundant transcripts caused by Cuffdiff2 due to normalizations enacted.

3. The sample GSM5988437 (Female Hypo Experimental BiologicalReplicate1) failed to pass the Replicate concordance test (R=0.59-0.63), while all other replicates demonstrated R>0.8. To understand the Replicate concordance test, please consult with the ENCODE guidelines on RNA-Seq (https://www.encodeproject.org/about/experiment-guidelines/).

Since all reads were found to be above phred score 22 with no adapter contamination or over-represented sequences P46-LL882.883. , the FPKM levels in the sample GSM5988437 (Female Hypo Experimental BiologicalReplicate1) were nullified somehow (by Cuffdiff2 or in another way) for a large number of transcripts with 1.to-10 FPKM levels, as compared to other samples.

This might result in a spurious identification of DEGs in the female hypothalami of exposed animals and in a spurious identification of sex differences between male and female hypothalamic DEGs (please see Supplementary File 1 with FPKMs, and the Figures 1D, 2A, and 2D).

To deal with these complications,

3a. The authors must repeat the concordance test for their FPKM estimates in the individual samples, providing the results of the concordance test in the supplement Table in accordance with the ENCODE guidelines.

3b. It is necessary to repeat a Cuffdiff analysis for hypothalamic RNA-seq samples without the affected replicate GSM5988437 (Female Hypo Experimental BiologicalReplicate1), to ensure unbiased data for the Figure 2A, Figure 2D, DEG lists and the "Transcriptome analyses identify brain region-specific gene signatures…" subsection.

3c. The subsection "Transcriptome analyses identify brain region-specific gene signatures…" must be improved after a Cuffdiff-mediated reanalysis of the RNA-Seq data.

4. ANOVA results gained by the authors are an important addon improving the text comprehensibility. Nevertheless, there is no need to incorporate ANOVA Tables into the manuscript text. Please move these Tables to Supplement.

5. Please identify units of measurements for X- and Y-axes on the subpanels of Figures1D and 2A.

[Editors' note: further revisions were suggested prior to acceptance, as described below.]

Thank you for resubmitting your work entitled "Pivotal role of TTR in early life stressful experiences induced escalated aggressive behavior" for further consideration by eLife. Your revised article has been evaluated by Catherine Dulac (Senior Editor) and a Reviewing Editor.

The manuscript has been improved but there are some remaining issues that need to be addressed, as outlined below:

Thank you for the additional revisions to your manuscript. Upon further examination, it appears that some of the wording continues to be imprecise and is at times inappropriate in over-interpreting the significance of the findings to human behavior. Some specific suggestions are provided below but it would be useful to revisit the entire manuscript with an eye towards assuring there are no hidden biases and that are statements are both factually and grammatically correct.

Specific Comments:

1) The change in the title now makes it both awkward and inaccessible to the general reader. One suggestion is: "Early life stressful experiences escalate adult mouse aggressive behavior via changes in transthyretin expression and function" or some variation thereof.

2) There are many examples of statements regarding violence and criminality that are inappropriate or incorrectly attribute to causality. For instance, the second sentence of the Abstract states "Early life stressful experiences triggers adulthood violence and criminality". This suggests that anyone that experiences early life stress will grow up to be a violent criminal. It is critically important not to make blanket statements that can be misinterpreted as sound science.

3) It is unclear how there can be an "escalated aggressive phenotype", a phenotype cannot be escalated, perhaps instead state "enhanced aggressive phenotype" or "resulted in escalated aggressive behavior".

4) Second sentence of the Introduction – "Such aberrant behavioral patterns are also manifested in patients of multiple psychiatric disorders including schizophrenia and bipolar disorder (2, 3) necessitating the identification of predisposing factors and early intervention strategies." – this seems strongly stigmatizing of individuals with mental illness, the vast majority of whom are not violent. It also is not important to the current findings and I would recommend removing such references to mental illness entirely.

5) Introduction – "Brain region-specific long-term changes in Ttr gene expression and thyroid hormone (TH) availability was evident in PPS induced escalated aggressive male mice, circulating TH being unaltered". – should read "were evident".

6) Introduction – "….it is extremely important to understand the biological culprits underlying brutal shift of normal adaptive aggression to escalated and pathological form." – the words "culprits" and "brutal" are emotionally laden terms that are inappropriate when discussing research findings.

7) Introduction – "We selected the extreme phenotypes for better understanding of the behavior observed in human violent offenders and psychopathy" – it is important to limit the conclusions of your study to what you observed which was changes in mouse behavior. Given the enormous complexity and multifactorial nature of violence in humans, it behooves you to not try and make direct connections between your studies and humans in the absence of any evidence that similar mechanisms are at play in human violent offenders.

Essential revisions:

1) The role of the decreased hypothalamic T3 availability in the development of aggression is not causally established. In the current form of the paper, it is not clear whether the decrease of TTR expression causes the aggression via the regulation of hypothalamic T3 availability or the effect of altered TTR expression on the aggression and T3 availability are independent. Therefore, the authors should either focus on the role of TTR or they should provide additional evidence that lack of hypothalamic T3 is the cause of the development of aggression.

As suggested, the focus of the paper has been shifted to the role of TTR in aggression and accordingly title has been modified to “Pivotal role of TTR in early life stressful experiences induced escalated aggressive behavior” in the revised version of the manuscript.

2) The studies on epigenetic mediated inheritance of the aggression phenotype are not well supported and should be removed. The data on methylation of the TTR promoter are an important contribution and should remain but should not be considered as a diagnostic for aggressive behavior or the basis for trans-generational inheritance.

As suggested, studies on epigenetic mediated inheritance of the aggression phenotype that is Ttr promoter methylation changes in F1 generation (Figure 7B) has been removed. Behavioural and other molecular studies (Ttr and thyroid hormone gene expression) in F1 generation has been kept as an important observation (Figure 6) without claiming anything on epigenetic inheritance. However, if the reviewers wish the behavioural and molecular studies in F1 generation also to be removed, we are open to do the modifications.

3) Statistical analyses should be redone with factorial analyses to make the proper comparisons between sexes, stress condition and brain region.

Two way ANOVA with stress condition and brain region as factors (Revised Figure 3D, Figure 4F-4P and Figure 7) and three way ANOVA with sex, stress condition and brain region as factors (Revised Figure 3B and Figure 3C) has been performed and described in the result and figure legend section of the revised manuscript. ANOVA summary tables (Table 1 to Table 15) have also been incorporated in the main manuscript and details of ANOVA analyses for all the above mentioned figures have been given in Supplementary File 2.

4) RNA-Seq data should include greater in depth analyses.

As suggested, RNA-Seq data has been analysed in depth, addressing the following points which are also incorporated in revised manuscript.

• GEO data has been updated with FPKM values of individual transcripts in individual biological replicate of control and experimental groups. Private token for access of reviewers has been provided.

• Single figure panel depicting parallel changes in transcripts in PFC and Hypo has been included (Figure 1D and 2A)

• Cell type specificity analyses of Male DEGs by Barres Lab database (brainrnaseq.org) and Allen brain transcriptomic atlas (https://portal.brain-map.org/atlases-and-data/rnaseq) has been included as Figure 1—figure supplement 2.10.

5) The images of immunohistochemistry staining are not of sufficient quality and should be replaced.

As suggested, IHC images of better quality and clarity has been included in the revised manuscript (Figure 3E-3H)

Reviewer #1 (Recommendations for the authors):

Deiodinase 3 (Dio3) plays critical role in the regulation of hypothalamic T3 availability, therefore, the D3 expression should also be determined along the other TH regulated genes.

According to my knowledge, there is no reliable antibody against dio2. This fact highly questions the validity of the dio2 ELISA data. The only accepted method for the determination of DIO2 protein level is the DIO2 enzyme assay (PMID: 24001133).

It is very true that Deiodinase 3 (Dio3) plays crucial role in the regulation of hypothalamic T3 availability and therefore, determining the D3 expression would provide further insights into TH signalling in brain. In view of the comment 1 under “Essential revision” section, we have now shifted the focus of the manuscript on TTR and limited our claims on TH availability. However, authors appreciate suggestion of the reviewer and future studies will include both Dio2 and Dio3 with more advanced methods in the context of aggressive behaviour.

The quality of the images illustrating immunocytochemistry is very weak. Better images would be necessary. In addition, it should be described in more details where the TTR-immunoreactive cells were observed in the hypothalamus. In addition, colocalization study should be performed at least with neuronal, glial and endothelial markers to determine what kind of cells express TTR in the hypothalamus.

As suggested, IHC images of better quality and clarity has been included in the revised manuscript (Figure 3E-H). Hypothalamic regions in which TTR immunoreactive cells were observed, have been described in result section of Figure 3E.

Cells expressing TTR in hypothalamus have been highlighted in Figure 3E. Barres Lab database (brainrnaseq.org) analyses revealed presence of Ttr mRNA in different mouse brain cell types (endothelial cells> neurons> microglia> astrocytes> oligodendrocytyes). We have given this information in Figure 1—figure supplement 4. However, co-localization studies with appropriate neuronal, glial and endothelial markers for TTR protein are warranted to confirm the cell type specificity.

Reviewer #2 (Recommendations for the authors):

Suggestions for improvement

1.1. In the present form submitted to GEO database, the transcriptomic data on the levels of individual transcripts are reported in aggregate per each experimental group, with no possibility to estimate variance between individual biological replicates. Reporting raw transcriptomic data on the transcript levels in each biological replicate will be beneficial for the readers interested in reanalyses of the data.

It is indispensable to report key elements of the transcriptomic analysis in full. Please update the GEO data with a single file on all levels of individual transcripts in each individual biological replicate.

As suggested GEO data has been updated with FPKM values of all transcripts in individual biological replicate of all the control and experimental groups.

In addition, Supplementary File 1 containing FPKM values of individual transcripts in individual biological replicate of control and experimental groups has been incorporated in the revised manuscript.

1.2. The personal experience of the reviewer evidences that old Cufflink-Cuffdiff Tuxedo bioinformatic pipeline for RNA-Seqs (done after Tophat or STAR aligners), if applied properly, is highly sensitive to individual mRNA levels with minor shares. Nevertheless, in the present form of the description of bioinformatics pipeline in the "Methods" section, with missed basic options applied by the authors for Tophat, Cufflink, Cuffdiff utilits, it is hard to estimate the validity of the transcriptomic analyses done by the authors. Please consider improving the description of the bioinformatics methods applied.

As suggested, description of the bioinformatics methods applied in RNA seq data analyses has been improved and incorporated in Materials and methods section of the revised manuscript.

1.3. Technical figures like volcano plots are important to control the quality of RNA-Seq, but are too uninformative to demonstrate the results of the experiment. Please consider illustrating the differences between PFC and HPT changes in DEGs on a single figure panel by depicting the parallel changes in the transcript levels identified in PFC (for example, x-axis) and HPT (for example, y-axis). At the same time, it is possible to move the volcano plots to the supplementary figures.

As suggested a single figure panel (Figure 1D) depicting the parallel changes in transcript levels (log2 Fold change) identified in PFC (x-axis) and Hypo (y-axis) has been incorporated and volcano plots has been moved to Figure 1—figure supplement 1 (Male DEGs) and Figure 2—figure supplement 1 (Female DEGs) of the revised manuscript.

1.4. Preliminary analysis done by the reviewer by applying tissue cell deconvolution methods evidenced in favour of possible specific trends in cell numbers in the limbic brain regions studied. In particular, it cannot be excluded that the juvenile stress episodes suppressed microglia numbers in both the PFC and the HPT.

Please consider providing the data on tissue cell composition in RNA-Seq individual samples with respect to the experimental groups.

(For details, see Sutton et al. 2022 https://doi.org/10.1038/s41467-022.28655-4 or other)

Such an analysis might be illustrative that stress itself is necessary but not sufficient to induce changes in the aggressive behaviours in affected male mice demonstrated excessive violence.

Applying tissue cell deconvolution methods in our RNA-seq data is indeed a very good suggestion and will definitely provide deeper insights. Due to technical constraints, we could not perform the exact analyses asked for, though we have provided the cell type data for PFC and Hypo DEGs in Figure 1—figure supplement 2-10 using other databases Barres Lab database (brainrnaseq.org) and Allen Brain Map website (https://celltypes.brain-map.org/).

1.5. Please consider additional analyse of the data on the individual transcripts levels reported in the RNA-Seq analysis in a way similar to 2-way or 3-way (when appropriate) factorial ANOVA, to identify possible additive and non-additive patterns of changes in the levels of transcripts.

The RNA sequencing experiment was performed using three replicates for each condition. All these replicates were then used for analysis using the Tuxedo pipeline which uses Cuffdiff2 to determine differentially expressed genes.

Cuffdiff2 combines a β distribution model to account for the count uncertainty arising from the high sequencing depth while it uses a negative binomial distribution to take care of the over dispersion problem that may exist during sequencing. The algorithm counts the fragments aligned per kb of exon per million reads which is then used to calculate p value. Moreover, it performs a Benjamini-Hochberg correction on the p value to account for false-discovery rate. After the whole procedure, the corrected p value (also labelled as q value in the output) is a statistically trusted significance measure. Therefore, ANOVA may not be required on this value as it would not account for FDR. Further, it will not be able to account for biological noise which do not show up in the FPKM values of individual replicates. The noise is accounted for by analysis of raw alignment data which Cuffdiff2 performs. The algorithm is used widely in the field and currently has more than 3000 citations.

The authors completely agree with reviewers that ANOVA is essential when analyzing qRT-PCR data and therefore, it has been included in relevant graphs of the revised manuscript.

1.6. It will be also interesting to know on cell specificity of DEGs identified in PFC. Please consider checking the prevalence of individual DEGs with the help of Allen brain transcriptomic atlas (cited in Yao et al.,2021) or another one.

The above mentioned atlas can be assessed by the following link (https://portal.brain-map.org/atlases-and-data/rnaseq)

Investigating cell type specificity of DEGs is indeed a very good idea. However, not all genes from our RNA seq analysis can be searched in publicly available data since we performed Total RNA sequencing while majority of researchers only perform polyA+ RNA sequencing. We have extracted the data on cell type of top ranking PFC and Hypo DEGs from Barres Lab database (brainrnaseq.org) which was made by bulk sequencing of separated cells. Also, we checked cell type specificity of some top ranking PFC DEGs in Allen Brain Map website (https://celltypes.brain-map.org/) and generated heat map. The figures are provided as Figure 1 —figure supplement 2.10 and mentioned in result section.

1.7. Please consider redrawing all figure panels depicted as traditional bars with the (Median, IQR, SD) box plots with individual data dots depicted. In particular, it is possible to achieve this at ease with the JASP freeware (https://jasp-stats.org/)

As suggested, individual data points have been depicted in figure panels.

1.8. It is arguable to use GAPDH as a reference for qPCR assay in the experiments with stress exposures, since GAPDH levels might be affected and even programmed by stress experienced by animals.

Please consider providing a rationale on the applying of GAPDH as an internal standard for mRNA levels instead of β-actin or other mRNAs conventional for stress studies.

We did not get any change in GAPDH mRNA levels in our transcriptome data and as per previous literature related to our experimental regime and other stress exposure studies (some references given below) in rodents, GAPDH is widely being used as an internal control. Therefore, we used GAPDH as an internal standard for mRNA levels. Ttr mRNA expression was also normalized with 18s rRNA but the results did not vary (data not shown).

References

Kathleen E. Morrison et al. (2016) Peripubertal Stress With Social Support Promotes Resilience in the Face of Aging, Endocrinology. 2016 May; 157(5): 2002–2014.

Jordi Tomas‐Roig et al. (2018). Effects of repeated long‐term psychosocial stress and acute cannabinoid exposure on mouse corticostriatal circuitries: Implications for neuropsychiatric disorders. CNS Neurosci Ther 24(6): 528–538.

Du, P et al. (2020). Chronic stress promotes EMT-mediated metastasis through activation of STAT3 signaling pathway by miR-337-3p in breast cancer. Cell Death Dis 11, 761.

1.9. For future studies.

In the absence of cross-fostering experimental schedule for F1 experiments, it is hard to delineate the origins of epigenetic changes identified in the F1 descendants, whether these changes were transmitted directly or indirectly, via mother's specific behaviours on the descendants.

We have limited our claims in the revised version of the manuscript and studies on epigenetic mediated inheritance of the aggression phenotype has been removed. However, further studies with cross-fostering experimental schedule till F2 generation are in progress to elucidate the molecular basis of inheritance in behaviour.

Recommendations for improving the writing and presentation

2.1. The present description of animal procedures in the "Methods" section does not provide enough details on the environments, in which the experimental mice were grown up. Please provide all details on animal housing procedures that might be stressful (social isolation events, social crowding events, numbers of animals per cell etc)

ARRIVE guidelines for animal research has already been uploaded as supporting document. Also, we have included details in the “Methods” section describing animal housing environments.

2.2. The description of several methods must be improved to clarify details critical for data comprehension.

For example,

Please indicate the details of screening tests done with subjects (females and anesthetized intruder) that were attacked by mice with pathological aggressive behaviours.

The experimental schedules must be reported in a clear way also. In particular, a brief statement must be done on how the distinct populations of "adult control" and "PPS adult male" mice screened by the Authors were originated from. A similar clarification must be done for female mice groups also.

As suggested, methods have been described in detail for better understanding of the readers.

2.3. The initial two paragraphs in the "Introduction" section do not provide the linear story tale on violence, "escalated" violence and a difference of these two concepts of aggression. The logic of this section must be improved.

2.3a. In the introduction, it looks like that the authors consider an aggression trait as a behavioural continuum between "zero-level" aggression to appropriate violence, and then to "escalated" aggressive behaviour. This point of view is arguable since it cannot be excluded that aggression is a multidimensional trait. Please consider revising and clarifying.

2.3b. It is possible to criticize "escalated" aggressive behaviour as unproductive. However, please do not make generalized negative statements on the nature of general violence in the Introduction and in the Discussion. Such statements might blackmail protective types of aggression critical for survival in mammals and humans.

Introduction and Discussion section has been revised in light of the above comments.

2.4. Please avoid to made generalized statements that "A can lead to B". Such statements with a strong modal verb "can" are highly misleading since the development of behavioural traits is not linear and depends on both genotype and environmental context often. Better to speak that "A might lead to B" under certain circumstances.

Several examples:

2.4a. "Escalated aggressive behavior … can lead to antisocial and criminal activities" [P3-S1-LL1-2]. This is misleading for numerous animals and for specific types of aggressive behaviours in H. sapiens. Please consider revising.

The statement has been revised in light of the above comment.

2.4b. "…pathological aggression has emerged as a consequence of early life adversities…" [P3-S2-LL9-10]. This statement is misleading, as it blames unrightfully all children affected by harsh life. Please consider revising.

The statement has been revised in light of the above comment.

2.4c. "…we inferred that Ttr promoter methylation could serve as a predictor of … behavioral deficits." [P24-S1-LL14-15]. – Better to talk about "possible behavioural deficits", not about "behavioural deficits".

As suggested in the essential revision section, the entire statement has now been removed in text of revised manuscript.

The above-mentioned list of examples is not exhaustive.

Reviewer #3 (Recommendations for the authors):

1) Examination of some of the DEGs from the other pathways identified in the KEGG analysis would either strengthen the argument that changes are specific to the TH pathway or highlight that changes are more wide spread.

Considering the key role of TTR in TH transport and KEGG analyses showing TH as the topmost ranking pathway, we explored TH in detail. However, DEGs from other pathways are worth exploring in the future. It is also possible that TH directly or indirectly influences the other pathways as well.

2) Some discussion of how a change in the amount of T3 and T4 could lead to aberrant aggression would enhance the manuscript.

Considering the comments mentioned in the “Essential revision” section, the focus of the manuscript has now been shifted to TTR and it is important to mention that further studies are necessary to identify the precise molecular pathway responsible for aberrant aggression. Therefore detailed discussion on role of T3 and T4 has been not been incorporated in the revised manuscript.

3) Overall the sex difference, which is profound, is not given much attention. There are more DEGs identified in females subject to peripubertal stress than males, yet there is no change in aggressive behavior, so what does this tell us? Also, does testosterone play a role in the sex difference in both the transcriptome and the behavioral changes? If males were gonadectomized, would the same transcriptional profile be apparent and would the behavior also be there? Or would the two endpoints diverge, belying the noting that there is a casual connect between them.

The authors really appreciate the suggestion to investigate sex differences in detail, though this was not the prime focus of this manuscript. However, work with similar objectives suggested by reviewer is in progress.

4) The peri-pubertal stress was conducted on 7 random days from PN28 to PN42. The timing of puberty is different in males and females, being earlier in females. Where measures taken to determine the stage of puberty in each animal (i.e. vaginal opening, preputial separation)? Did the stress impact the timing of puberty?

The timing of peripubertal stress exposure from PN28-PN42 is well established both for male and female mice {Marquez et al. (2013) Transl Psychiatry,15;3(1):e216; Kathleen E. Morrison et al. (2016) Endocrinology.157(5): 2002–2014;Morató L et al. (2022) Sci Adv, 8(9):eabj9109.}

However, we also assessed pubertal onset in males through preputial separation and in females by observing vaginal opening and determining start of estrous (sexual) by crystal violet staining based vaginal cytology. We did not observe any conspicuous changes in the above mentioned visible signs of pubertal onset in both sexes, though it is indeed a very interesting point and worth exploring in detail in future studies.

5) The transgenerational assertions should either be dropped or the study carried out to the F2 generation.

We have not made any transgenerational assertions, rather referred to the observations as “intergenerational inheritance” as we conducted experiments only in F1 generation. However, studies related to epigenetic mediated inheritance of aggressive phenotype has been removed in the revised version of manuscript

6) How was the use of MeDIP specific to the promoter for Ttr?

Ttr proximal promoter (-184 to -33 bp from TSS) was amplified with specific primers (Supplementary File 3) generating a 151 bp product in MedIP-qPCR and this has been mentioned in methods section.

7) What are the circulating androgen levels in the males from the various groups? Could the PP have altered the HPGA that then in turns alters behavior?

The primary objective of the study was to identify brain region specific transcriptional responses and as we did not find androgen signalling genes in our top DEGs, we did not dig down further. However, high testosterone and testosterone/cortisol ratio being considered important for aggression, circulating androgens and PPS triggered changes in HPGA axis is worth exploring in future studies.

8) It does not seem appropriate to refer to "donut shaped cells".

The term "donut shaped cells" has been removed from the text.

9) Figure 4J – appears mislabeled, has Hypo twice and no PFC.

Mislabelling has been corrected in revised Figure 4J.

10) Figure 4L-P – why aren't the individual points plotted for the mRNA and protein.

As suggested individual points have been plotted for the mRNA and protein in revised Figure 4L-P.

11) Whenever both brain areas are considered the statistics should be 2-way ANOVA with brain region and treatment as factors?

As suggested 2-way ANOVA has been performed with brain region and treatment as factors (Revised Figure 3D, Figure 4F-4P and Figure 7) and 3-way ANOVA has been performed with sex, brain region and treatment as factors (Revised Figure 3B and Figure 3C).

12) The word "trauma" in the context used here connotes an emotional interpretation of stressful or fearful events. We do not know if the mice are experiencing trauma, instead we know they are being subject to fearful and stress-inducing experiences. It is suggested that the word trauma be removed throughout and replaced with more precise terminology.

As suggested the word “trauma” has been replaced with “stressful experiences” throughout the text of the revised manuscript.

[Editors' note: further revisions were suggested prior to acceptance, as described below.]

The manuscript has been improved but there are some remaining issues that need to be addressed, as outlined in detail below:

The following additional revisions are essential:

1) The conclusions on large sex differences between male and female hypothalamic transcriptomes are not supported by the data due to the problematic sample GSM5988437 (Female Hypo Experimental BiologicalReplicate1). This sample must be withdrawn from the analysis, and the male-female comparisons for hypothalamic transcriptomes must be re-estimated without this problematic sample. This will help to dismiss any spurious claims about the "augmented" male-female hypothalamic differences.

We really appreciate the suggestion of reviewers and thus we have withdrawn Female Hypo Experimental Biological Replicate1 sample, redone the analysis and revised the relevant figure (Figure 2) in the manuscript. We have included the updated list of DEGs in females along with earlier male DEGs in supplementary file 1. Also, we mentioned in the RNA seq method section that replicates that passed the concordance test (R> 0.8) were included for analysis of DEGs.

2) Data about the specificity of the TTR immunocytochemistry and D2 elisa would be absolutely necessary. The authors describe the changes of 2 TH transmitters and Dio2. Dio3 is at least an important regulator of TH availability in the brain as Dio2. So either data about Dio3 expression should be added or data about the expression of TH transporters and Dio2 should be removed.

No primary TTR antibody control was used to determine the specificity of TTR immunofluorescence and included as supplementary data (Figure 3—figure supplement 1) Moreover, the same TTR primary antibody has been used for TTR western blotting that gave the precise molecular weight band further confirming the specificity.

Considering the need for timely publication of our main finding that TTR is involved in stress driven aggression, we are removing the data about the expression of TH transporters and Dio2 and accordingly revised figures have been included (Figure 4. Figure 5 and Figure 6).

However, future studies focussing on role of brain TH availability in aggression and more importantly determining Dio2 and Dio3 levels with more valid and advanced methods is definitely on priority.

3) Determination of the cell types expressing TTR should be very fast and easy with double-labeling immunocytochemistry and would increase the value of the paper, however, this is not absolutely necessary to support the conclusions of the paper.

We absolutely agree that determination of the cell types expressing TTR would increase the value of the paper. However, we have shifted our laboratory to a new place and the project involving this study is also completed. Therefore, doing the co-localization experiments will take more time than anticipated. In our future planned experiments co-localization of TTR protein with appropriate neuronal, glial and endothelial markers in mouse brain is definitely on the priority list.

[Editors' note: further revisions were suggested prior to acceptance, as described below.]

The manuscript has been improved but there are some remaining issues that need to be addressed, as outlined below:

Thank you for the additional revisions to your manuscript. Upon further examination, it appears that some of the wording continues to be imprecise and is at times inappropriate in over-interpreting the significance of the findings to human behavior. Some specific suggestions are provided below but it would be useful to revisit the entire manuscript with an eye towards assuring there are no hidden biases and that are statements are both factually and grammatically correct.

Specific Comments:

1) The change in the title now makes it both awkward and inaccessible to the general reader. One suggestion is: "Early life stressful experiences escalate adult mouse aggressive behavior via changes in transthyretin expression and function" or some variation thereof.

As suggested the title has now been modified to “Early life stressful experiences escalate aggressive behavior in adulthood via changes in transthyretin expression and function”

2) There are many examples of statements regarding violence and criminality that are inappropriate or incorrectly attribute to causality. For instance, the second sentence of the Abstract states "Early life stressful experiences triggers adulthood violence and criminality". This suggests that anyone that experiences early life stress will grow up to be a violent criminal. It is critically important not to make blanket statements that can be misinterpreted as sound science.

As suggested the statement has been modified to “. Early life stressful experiences might increase the risk of developing pathological aggressive behavior in adulthood

3) It is unclear how there can be an "escalated aggressive phenotype", a phenotype cannot be escalated, perhaps instead state "enhanced aggressive phenotype" or "resulted in escalated aggressive behavior".

As suggested the above mentioned phrase of “escalated aggressive phenotype in the abstract has been changed to “resulted in escalated aggressive behavior”. In addition escalated aggressive phenotype has been changed to escalated aggressive behavior throughout the text.

4) Second sentence of the Introduction – "Such aberrant behavioral patterns are also manifested in patients of multiple psychiatric disorders including schizophrenia and bipolar disorder (2, 3) necessitating the identification of predisposing factors and early intervention strategies." – this seems strongly stigmatizing of individuals with mental illness, the vast majority of whom are not violent. It also is not important to the current findings and I would recommend removing such references to mental illness entirely.

As suggested the above references have been removed.

5) Introduction – "Brain region-specific long-term changes in Ttr gene expression and thyroid hormone (TH) availability was evident in PPS induced escalated aggressive male mice, circulating TH being unaltered". – should read "were evident".

As suggested the statement has been corrected to Brain region-specific long-term changes in Ttr gene expression and thyroid hormone (TH) availability were evident in PPS induced escalated aggressive male mice, circulating TH being unaltered.

6) Introduction – "….it is extremely important to understand the biological culprits underlying brutal shift of normal adaptive aggression to escalated and pathological form." – the words "culprits" and "brutal" are emotionally laden terms that are inappropriate when discussing research findings.

As suggested the statement has been modified as “it is extremely important to understand the biological factors contributing to shift of normal adaptive aggression to escalated and pathological form"

7) Introduction – "We selected the extreme phenotypes for better understanding of the behavior observed in human violent offenders and psychopathy" – it is important to limit the conclusions of your study to what you observed which was changes in mouse behavior. Given the enormous complexity and multifactorial nature of violence in humans, it behooves you to not try and make direct connections between your studies and humans in the absence of any evidence that similar mechanisms are at play in human violent offenders.

The comment is highly appreciated and accordingly all the sentences including “We selected the extreme phenotypes for better understanding of the behavior observed in human violent offenders and psychopathy” making direct connections with human violence behavior have been removed.

Author details

1. Rohit Singh Rawat

   CSIR-Institute of Genomics & Integrative Biology, New Delhi, India

   Contribution

   Data curation, Formal analysis, Validation, Investigation, Methodology, Writing - original draft

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0001-6329-6636

2. Aksheev Bhambri

   1. CSIR-Institute of Genomics & Integrative Biology, New Delhi, India
   2. Academy of Scientific and Innovative Research (AcSIR), New Delhi, India

   Contribution

   Data curation, Software, Formal analysis, Validation, Investigation, Visualization, Methodology

   Competing interests

   No competing interests declared

3. Muneesh Pal

   CSIR-Institute of Genomics & Integrative Biology, New Delhi, India

   Contribution

   Data curation, Formal analysis, Validation, Investigation, Visualization, Methodology

   Competing interests

   No competing interests declared

4. Avishek Roy

   All India Institute of Medical Sciences, New Delhi, India

   Contribution

   Data curation, Investigation, Methodology

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0003-2633-487X

5. Suman Jain

   All India Institute of Medical Sciences, New Delhi, India

   Contribution

   Resources, Supervision

   Competing interests

   No competing interests declared

6. Beena Pillai

   1. CSIR-Institute of Genomics & Integrative Biology, New Delhi, India
   2. Academy of Scientific and Innovative Research (AcSIR), New Delhi, India

   Contribution

   Resources, Funding acquisition, Writing - review and editing

   Competing interests

   No competing interests declared

7. Arpita Konar

   1. CSIR-Institute of Genomics & Integrative Biology, New Delhi, India
   2. Institute of Health Sciences, Presidency University, Kolkata, India

   Contribution

   Conceptualization, Resources, Data curation, Formal analysis, Supervision, Funding acquisition, Validation, Investigation, Visualization, Methodology, Writing - original draft, Project administration, Writing - review and editing

   For correspondence

   konar.arpita24@gmail.com

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-1761-1065

• 1,919

  Page views

• 327

  Downloads

• 0

  Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.