Skip to main content

ORIGINAL RESEARCH article

Front. Psychiatry, 24 November 2020
Sec. Molecular Psychiatry
This article is part of the Research Topic In Search of Mechanisms: Genes, Brains, and Environment in Aggressive Behavior View all 9 articles

Variants of the Aggression-Related RBFOX1 Gene in a Population Representative Birth Cohort Study: Aggressiveness, Personality, and Alcohol Use Disorder

\nMariliis VahtMariliis Vaht1Kariina LaasKariina Laas1Nolia Fernndez-Castillo,,,Noèlia Fernàndez-Castillo2,3,4,5Triin KurrikoffTriin Kurrikoff6Margus KanarikMargus Kanarik1Stephen V. FaraoneStephen V. Faraone7Liina-Mai ToodingLiina-Mai Tooding6Toomas VeidebaumToomas Veidebaum8Barbara FrankeBarbara Franke9Andreas ReifAndreas Reif10Bru Cormand,,,Bru Cormand2,3,4,5Jaanus Harro
Jaanus Harro1*
  • 1Division of Neuropsychopharmacology, Department of Psychology, Estonian Center of Behavioral and Health Sciences, University of Tartu, Tartu, Estonia
  • 2Departament de Genètica, Microbiologia i Estadística, Facultat de Biologia, Universitat de Barcelona, Barcelona, Spain
  • 3Centro de Investigación Biomédica en Red de Enfermedades Raras (CIBERER), Instituto de Salud Carlos III, Madrid, Spain
  • 4Institut de Biomedicina de la Universitat de Barcelona (IBUB), Barcelona, Spain
  • 5Institut de Recerca Sant Joan de Déu (IR-SJD), Esplugues de Llobregat, Spain
  • 6Institute of Social Studies, University of Tartu, Tartu, Estonia
  • 7Departments of Psychiatry and of Neuroscience and Physiology, The State University of New York Upstate Medical University, Syracuse, NY, United States
  • 8National Institute for Health Development, Estonian Center of Behavioral and Health Sciences, Tallinn, Estonia
  • 9Department of Human Genetics, Radboud University Medical Center, Donders Institute for Brain, Cognition and Behavior, Nijmegen, Netherlands
  • 10Department of Psychiatry, Psychosomatic Medicine and Psychotherapy, University Hospital Frankfurt Goethe University, Frankfurt am Main, Germany

Background: Recently, RBFOX1, a gene encoding an RNA binding protein, has consistently been associated with aggressive and antisocial behavior. Several loci in the gene have been nominally associated with aggression in genome-wide association studies, the risk alleles being more frequent in the general population. We have hence examined the association of four RBFOX1 single nucleotide polymorphisms, previously found related to aggressive traits, with aggressiveness, personality, and alcohol use disorder in birth cohort representative samples.

Methods: We used both birth cohorts of the Estonian Children Personality Behavior and Health Study (ECPBHS; original n = 1,238). Aggressiveness was assessed using the Buss–Perry Aggression Questionnaire and the Lifetime History of Aggressiveness structured interview at age 25 (younger cohort) or 33 (older cohort). Big Five personality at age 25 was measured with self-reports and the lifetime occurrence of alcohol use disorder assessed with the MINI interview. RBFOX1 polymorphisms rs809682, rs8062784, rs12921846, and rs6500744 were genotyped in all participants. Given the restricted size of the sample, correction for multiple comparisons was not applied.

Results: Aggressiveness was not significantly associated with the RBFOX1 genotype. RBFOX1 rs8062784 was associated with neuroticism and rs809682 with extraversion. Two out of four analyzed RBFOX1 variants, rs8062784 and rs12921846, were associated with the occurrence of alcohol use disorder.

Conclusions: In the birth cohort representative sample of the ECPBHS, no association of RBFOX1 with aggressiveness was found, but RBFOX1 variants affected basic personality traits and the prevalence of alcohol use disorder. Future studies on RBFOX1 should consider the moderating role of personality and alcohol use patterns in aggressiveness.

Introduction

Globally, more than 1.3 million people worldwide die each year owing to aggressive behavior and violence (either self-directed, interpersonal, and collective), accounting for 2.5% of mortality (1). Nevertheless, from the evolutionary perspective, aggression can be described as adaptive. Aggression belongs to the behavioral repertoire of most species. Humans are a rather highly aggressive species compared with other animals. This may be related to the high benefit-to-cost ratio for intraspecific aggression (2). Aggression in mammals, including humans, has a high heritability (3, 4). Struggle for resources (e.g., territory, suitable mates, food) must have substantially contributed to the shaping of aggression by favoring gene variants that promote agonistic behavior. However, when humans express their inherent aggressiveness in an unfitting context, this may lead to social maladjustment and crime (5).

Variance in many genes has been associated with aggressiveness, with very small effects of each in the large population studies by GWAS (6, 7). Recently, a novel candidate gene, RBFOX1, has been linked to aggressive behavior by convergent evidence from GWAS, epigenetic analyses, and neuroimaging genetics, gene expression, and animal models (8). RBFOX1 (RNA binding protein, Fox-1 homolog 1; also known as ataxin 2-binding protein 1, A2BP1, or hexaribonucleotide-binding protein 1, HRNBP1) encodes for the Fox-1 protein and regulates alternative splicing that controls gene expression and, in turn, coordinates neuronal brain activity (9, 10).

RNA binding proteins (RBPs) are key components in RNA metabolism (11) by influencing the structure and interactions of the RNAs and playing critical roles in their biogenesis, stability, function, transport, and cellular localization (12). Each RBP interacts with RNA in a unique sequence- or structure-specific manner. Alternative splicing is one of the central mechanisms regulating eukaryotic gene expression (13). RBPs coordinate elaborate networks of RNA–protein and protein–protein interactions that control RNA metabolism. Hence, alterations in their RNA-binding function could impact many genes and pathways, leading to complex, multifaceted phenotypes (11, 14). Mutations in RBPs cause and/or contribute to many human neurodevelopmental and neurologic disorders (11, 15, 16). Abnormalities in the Fox-1 encoding gene, RBFOX1, are associated with neurodevelopmental disorders (17). Variations in the RBFOX1 have been associated with anxiety (9), attention deficit/hyperactivity disorder (18), schizophrenia (19, 20), autism spectrum disorders (17), intellectual disability with epilepsy (21), and gray matter loss in Alzheimer's disease patients (22). Conclusively, potential effects of the RBFOX1 genotype are multiple and likely variable by sample specifics owing to environmental interactions.

Several variants of RBFOX1 have been associated with aspects of aggressiveness in a number of GWASs (23). The C allele of rs6500744, located within the first intron, was associated with conduct disorder symptoms in the interaction with mothers' warmth (24). Another SNP (rs8062784) in intron 1 of RBFOX1 was associated with anger in a GWAS assessing hostility (25), and a variant located in intron 3 of the gene (rs12921846) was associated with conduct disorder in a sample of ADHD trios (26). A meta-analysis of nine population-based GWASs including around 19,000 children provided evidence for the contribution of RBFOX1 to children's aggressive behavior (7). Four SNPs in the RBFOX1 gene (rs809682, rs12922093, rs12373031, and rs10521042, all located in the intron regions) showed suggestive associations. Of these, the rs809682 polymorphism was demonstrated to have the lowest association p-value, with the major T allele conferring the risk. Statistically significant association was found for rs809682 when comparing aggressive prisoners to controls and non-aggressive prisoners (23). It should be noticed though that the direction of the effect was not as in the original study.

We selected the SNPs showing the lowest association p-value from every single reported GWAS [rs6500744, rs8062784, rs809682, and rs12921846; see (23)] and assessed their association with aggressiveness and basic personality traits in a birth cohort representative sample of young adults. Given that antisocial behavior is the main predictor of alcohol (mis)use and the other way around [e.g., (27)], we also examined the occurrence of alcohol use disorder.

Materials and Methods

Study Population

The analysis was carried out on the Estonian Children Personality Behavior and Health Study (ECPBHS) sample, the original Estonian sample of the European Youth Heart Study (1998/99) that was subsequently incorporated into the longitudinal ECPBHS. All the subjects are of European descent. The principles of formation of the original sample and procedure of first data collection have been described in detail elsewhere (28). In brief, this is a representative birth cohort sample of the Tartu city and county with a school as the sampling unit. All schools of Tartu County, Estonia, that agreed to participate (54 of the total of 56) were included into the sampling using the probability proportional to the number of students of the respective age groups in the school, and 25 schools were selected. All children from grades 3 (younger birth cohort) and 9 (older birth cohort) were invited to participate. ECPBHS is population representative, while 79.1% of subjects of the randomized regional sample participated in the original sampling. Details on follow-ups have been described elsewhere [e.g., (29)], and the number of participants in the presented analyses is given below. The study was approved by the Ethics Review Committee on Human Research of the University of Tartu. Written informed consent was obtained from all participants and, in case of minors, also from their parents.

Measures

Aggressive Behavior

During the last data collection waves (at age 25 years in the younger cohort in 2014/2015 and at age 33 in the older cohort in 2016/2017), aggressive behavior was self-reported using the Buss-Perry Aggression Questionnaire (30). The 29-item self-report Buss-Perry Aggression Questionnaire (31) assesses four aspects of aggressive behavior: Physical aggression, Verbal aggression, Anger, and Hostility. Participants (n = 436 in the younger and n = 499 in the older cohort) rated each statement on a 5-point Likert Scale (uncharacteristic = 1, characteristic = 5). During the same data collection waves, the Life History of Aggression interview [LHA; (32)] was carried out by experienced clinical psychologists in order to assess dimensions of aggression (n = 427 in the younger and n = 495 in the older cohort). Items were scored only for the history of actual behavior (33). LHA has three subscales: Aggression (temper tantrums, physical fighting, verbal fighting, assaults on other people, and assaults on property); Consequences/Antisocial Behavior (school disciplinary problems, problems with supervisors at work, antisocial behavior not resulting in police involvement, and antisocial behavior involving the police); and Self-Directed Aggression (assaults on self and suicide attempts). Each item was rated on a 5-point scale, ranging from 0 = “no events” to 5 = “more events than can be counted.”

Personality

Personality traits of the five-factor model (34) were measured by self-reports at age 25 (n = 856) with EE.PIP-NEO (35), which is a semantically simplified 240-item version of the International Personality Item Pool (IPIP), which emulates the NEO-PI-R.

Alcohol Use Disorders

Assessment of lifetime occurrence of alcohol use disorders was based on DSM-IV and was carried out in both cohorts (n = 931) at age 25 by experienced clinical psychologists using the Mini-International Neuropsychiatric Interview [M.I.N.I.5.0.0; (36); Estonian version; (37)].

Genotyping

Genomic DNA was extracted from whole blood samples using Qiagen QIAamp® DNA Blood Midi Kit. The quantitative real-time polymerase chain reaction (qRT-PCR) for genotyping the four SNP polymorphisms was performed using TaqMan Pre-Designed SNP Genotyping Assays (Applied Biosystems; Foster City, CA, United States) containing primers and fluorescent probes. For rs809682, the Assay C___8926788_10 was utilized; for rs8062784, rs12921846, and rs6500744 polymorphisms, the Assays C__29081048_20, C__32104163_10, and C___3008571_10 were used, respectively. Genotyping reactions were performed in a total volume of 10 ml with ~25 ng of template DNA. QRT-PCR reaction components and final concentrations were as follows: 1:5 5 x HOT FIREPol® Probe qPCR Mix Plus (ROX) (Solis BioDyne) and 1:20 80 × TaqMan Primers Probe.

Context sequences [VIC/FAM] were as follows: rs809682–TAACAAACTACAGCCTAATTTAGTA[A/T]GTGAACTAAGTGAAAGCCTAACTTG, rs8062784–TTCTGTGAACCAACACTTTCTTTTC[A/T]TTGGGTTTGATACAGTGGCATCAAT,rs12921846–ATCTTGGAAAGCATTTGTTATTTCA[A/T]ACTCTTCAAATCTGCAAGTCTTACA, rs6500744–GCTTACCATTTATTTTATTTTCAGG[C/T]GGTTGTATTCATTATAATGCCATTA.

Reactions were performed on the Applied Biosystems ViiA™ 7 Real-Time PCR System. The amplification procedure consisted of an initial denaturation step at 95°C for 12 min and 40 cycles of 95°C for 15 s and 60°C for 1 min. Positive and negative controls were added to each reaction plate. No inconsistencies occurred. Genotyping was performed blind to all phenotypic data. Allele frequencies agreed with the National Center for Biotechnology Information database and published reports. Genotype frequencies were in Hardy-Weinberg equilibrium and are shown in Table 1.

TABLE 1
www.frontiersin.org

Table 1. RBFOX1 genotype frequencies in the ECPBHS sample.

Statistical Analysis

Birth cohorts were pooled for analysis. Categorical variable (genotype) relations to continuous variables were explored with analysis of variance (ANOVA) and presented as F-statistic, raw p-value and eta-squared (η2) as a measure of effect size. Fisher's least significance difference method (LSD) was used in all post hoc comparisons. Contrasts were calculated for significant model effects. Chi-square tests were conducted to assess the distribution of Alcohol Use Disorder by genotype, and by genotype and gender, and presented as χ2-statistic and raw p-value. For fitting of the path model, we used the AMOS package of structural equation modeling (SEM) by the MCMC (Markov chain Monte Carlo) method. Genotype was entered into the SEM model as a dichotomous variable: A/A homozygotes vs. T-allele carriers. All p-values are reported as two-tailed, and results are considered significant at the conventional p < 0.05 level; correction for multiple testing is not applied. Statistical analyses were performed using IBM SPSS Statistics, Version 25.

Results

The Selected RBFOX1 Polymorphisms and Aggressive Behavior

Aggressiveness assessed either using the self-report Buss-Perry Aggression Questionnaire (Table 2) or the Lifetime History of Aggression interview (Table 3) was not associated with any of the RBFOX1 polymorphisms. This was the case for total scores as well as subscales. No genotype by gender effect was found either.

TABLE 2
www.frontiersin.org

Table 2. RBFOX1 genotypes and the Buss-Perry Aggression Questionnaire at age 25.

TABLE 3
www.frontiersin.org

Table 3. RBFOX1 genotypes and the life history of aggression interview at age 25.

Personality and RBFOX1 Genotype

Neuroticism was associated with the RBFOX1 rs8062784, being lower in homozygotes for the aggressiveness risk allele A (Table 4). RBFOX1 rs809682 was associated with extraversion, while the homozygotes for the risk allele T having higher scores. No other polymorphism was statistically significantly associated with neuroticism or extraversion.

TABLE 4
www.frontiersin.org

Table 4. The effects of RBFOX1 genotypes on neuroticism and extraversion at age 25.

RBFOX1 Polymorphisms and the Occurrence of Alcohol Use Disorder

Alcohol use is a most salient mediator to aggressiveness, so the association of RBFOX1 polymorphisms with lifetime prevalence of alcohol use disorder by age 25 was examined. The overall genotype effect (χ2 = 4.14; p = 0.042) was revealed for RBFOX1 rs8062784 (Figure 1): This was largely based on male subjects carrying the less frequent low aggressiveness allele T who had almost twice higher risk of alcohol abuse. While analyzing males and females separately, we found that rs12921846 (Figure 2) was associated with alcohol use disorder in females (χ2 = 4.22; p = 0.045). Female RBFOX1 rs12921846 homozygotes for the less frequent T-allele, also has been related to lower aggressiveness, had higher alcohol abuse risk.

FIGURE 1
www.frontiersin.org

Figure 1. Association of RBFOX1 rs8062784 with alcohol use disorder in males. Males, χ2 = 7.01; p = 0.008 (n = 411); females, χ2 = 0.45; p = 0.509 (n = 520). Lifetime diagnosis is based on MINI interview at age 25. AUD in 158 males out of 413 and 54 females out of 523 (in total, 212 out of 936).

FIGURE 2
www.frontiersin.org

Figure 2. Association of RBFOX1 rs12921846 with alcohol use disorder in females. Males, χ2 = 0.41; p = 0.815 (n = 411); females, χ2 = 6.22; p = 0.045 (n = 520). Lifetime diagnosis is based on MINI interview at age 25. AUD in 158 males out of 413 and 54 females out of 523 (in total, 212 out of 936).

Alcohol Use Disorder, Sociodemographic Factors, Personality, and the RBFOX1 rs8062784 Genotype

Next we examined the association of several factors potentially related to alcohol use disorder at age 25, which is the last observation for both cohorts. AUD was more prevalent in males: 10% of women and 38% of men had experienced AUD by age 25 (χ2 = 102.8; df = 1; p < 0.001). In these birth cohorts, lifetime alcohol use disorder was not associated with income in either male or female subjects (data not shown), probably owing to early age. At age 25, 58% of the participants were living with a spouse, but there was no difference in lifetime occurrence of AUD between subjects who lived with their spouse (21%) compared with subjects who did not (25%). However, the lifetime prevalence of AUD was higher (15%; n = 29) in females who did not live with a spouse compared with those who did (8%; χ2 = 6.3; df = 1; p = 0.01). There was a higher prevalence of AUD among subjects with lower education level at age 25 (32 vs. 17%; χ2 = 28.99; df = 1; p < 0.001; males, 48 vs. 30%; χ2 = 15.35; df = 1; p < 0.001; females, 14 vs. 9%; χ2 = 3.85; df = 1; p = 0.050); but RBFOX1 genotype was not associated with education (data not shown).

No interaction between the gender factor and lifetime alcohol use disorder was found for the Big Five personality traits (Table 5), but participants with AUD had significantly higher neuroticism [F(1, 781) = 30.01; p < 0.001; η2 = 0.037], lower agreeableness [F(1, 774) = 9.39; p = 0.002; η2 = 0.012], and lower conscientiousness [F(1, 779) = 18.15; p < 0.001; η2 = 0.023]. Because the RBFOX1 rs8062784 genotype and occurrence of alcohol use disorder were both associated with neuroticism, we examined whether or not the association of the RBFOX1 genotype could be mediated by this personality trait. According to a Bayesian model by the Markov chain Monte Carlo method, the RBFOX1 rs8062784 direct effect on AUD was not significant (regression weight 0.19; Bayes' credible interval −0.069…0.455; probability 95%), the RBFOX1 rs8062784 effect on Neuroticism was 7.69 (Bayes' credible interval 1.84…13.64; probability 95%), and the neuroticism effect on AUD was 0.006 (0.0002…0.009), suggestive of the neuroticism mediated association of RBFOX1 with alcohol use disorder (Figure 3).

TABLE 5
www.frontiersin.org

Table 5. Personality traits at age 25 by gender and lifetime alcohol use disorder (AUD).

FIGURE 3
www.frontiersin.org

Figure 3. Mediation by neuroticism of the association of RBFOX1 rs8062784 and alcohol use disorder. Regression weights with credible intervals (95%) are shown.

Discussion

Four polymorphisms of the RBFOX1 gene, previously linked to aspects of aggressive behavior, were in the ECPBHS sample not associated with aggressiveness in young adulthood, neither by self-report nor interview measurement. Genome-wide association studies indicate a very small effect of each gene on aggressiveness in human population, so given the size of the sample, this may not be surprising. Nevertheless, when performing subgroup analyses by ethnicity, age of participants, sample characteristics, and outcome measures, significant associations emerge (6). Several variants in a number of candidate genes have been found strongly associated with aggression measures in the ECPBHS sample [see (29, 33) and references therein]. It should hence be concluded that in this specific population of which the ECPBHS sample is highly representative, other genes than RBFOX1 play a role in aggressiveness. Thus, other genetic variants than the ones inspected play a role in aggressiveness in this sample, although we cannot exclude other variants in RBFOX1 not investigated here.

Neuroimaging genetic studies support the notion that the RBFOX1 gene contributes to brain function and structure. The expression levels of RBFOX1 in the brain are highest in the anterior cingulate cortex that regulates emotions and social behavior (3841). The thickness of the ACC has been associated with aggressiveness (41, 42), and in individuals prone to aggression, its activity in response to provocation is attenuated (43). Furthermore, higher activity of ACC has been found in aggressive adolescents with disruptive behavior disorders (44). Yet anterior cingulate is involved in the whole complexity of emotion regulation, and the meaning of any alteration at the molecular level likely depends on the overall state of the whole complex [e.g., (45)]. Because previous studies that have implicated these specific variants of RBFOX1 were either on selected groups of participants or relied on specific measures of phenotype, we took the approach to study the association of RBFOX1 with basic personality traits that could provide insight into which factors play a role in the potential effect of RBFOX1. Two of the polymorphisms were nominally associated with such basic personality traits as neuroticism and extraversion.

Neuroticism was lower with two A-alleles of RBFOX1 rs8062784. Previously, this allele was associated with anger in a GWAS that was targeted at hostility (25). Anger or hostility is indeed a component of neuroticism in the five-factor model. For this reason, we separately analyzed the six facets of neuroticism and found four of them, including anger/hostility significantly associated with the genotype (data not shown). Interestingly, the weakest of the subscale associations was with impulsivity, often a facilitative factor in aggressive behavior. While this association of genotype and neuroticism could well-be a mere chance finding, its direction also supports the hypothesis that low neuroticism in the risk genotype in the present sample is a mechanism preventing overt aggressiveness: In population-derived samples, aggressive subjects appear to have higher neuroticism (4648).

RBFOX1 rs809682 was associated with extraversion: Homozygotes for the T-allele had higher scores. The meta-analysis of Pappa et al. (7) of nine population-based GWASs of 19,000 children suggested the major T allele to carry the aggression risk. Aggressiveness rating in these studies was based on maternal scores, and this may differ from self-assessment and interviews; moreover, levels of continuity of aggression from childhood to early adulthood have been found to be only moderate (49). Furthermore, in a study of adult aggressive prisoners, a significant association was found for rs809682; however, the direction of the effect was opposite as in another (7) study, this time with the minor A allele linked to aggression (8), although the sample size was very limited. These conflicting results may also find an explanation in the dual potential of extraversion to moderate aggressiveness. Extraversion is the tendency toward interaction with others, assertiveness, liveliness, and action-orientation (50, 51). While aspects of extraversion are defined as deriving from positive emotionality, another central part of extraversion is assertiveness. Frost (52) made the early notion that extraverts tend to be more assertive and lacking in submissiveness and self-criticism. Assertiveness and verbal aggression have been demonstrated to correlate positively (53), and extraversion has recently been shown to be positively associated with physical aggression (54). Indeed, in our exploratory analysis of the facets of extraversion, the RBFOX1 rs809682 T/T homozygocity was related rather to assertiveness, activity, gregariousness, and positive emotions facets but less to friendliness and excitement seeking (data not shown). Thus, somewhat speculatively, lower prevalence of rs809682 T-allele among less aggressive subjects in some studies may indirectly suggest that this allele is promoting assertive behavior, and the absence of it mitigates aggressive behavioral choices.

Two of the RBFOX1 variants were associated with alcohol use disorder, rs8062784 in males and rs12921846 in females. In both occasions, higher prevalence of alcohol use disorder was present with the minor variant previously associated with lower level aggression. Problematic alcohol use is often predictive of violent behavior, but in this sample, the risk alleles were not associated with aggressiveness. This result rather converges with the personality findings in that in the risk allele carriers, who constitute large majority of subjects, aggressiveness is not common owing to the generally positive side of these variants. According to this scenario, the RBFOX1 aggressiveness-related variants owe this association to other coinciding genetic or environmental factors. In the present sample, alcohol use disorder was associated with both rs8062784 genotype and neuroticism, and the genotype also with neuroticism. A path analysis supported the possibility of mediation by neuroticism between the genotype and alcohol use disorder. Speculatively, the major rs8062784 allele, promoting lower neuroticism and protective of alcohol use disorder in this sample, may in general be associated with proactive behaviors and thus become associated with anger in different types of environment where the positive side of agonistic behavior cannot be properly channeled.

Analysis simultaneously involving four gene polymorphisms and many comparisons inflates the possibility of false positive findings, and conventional correction for multiple testing would render all associations non-significant. This is a major limitation of the present longitudinal study, in that it is restricted in its sample size. RBFOX1 has previously been associated with aggression-related phenotypes in a number of GWAS. The present sample is much smaller, but nevertheless, any large direct effect of the genotype would have been detected. We do, however, suggest a few potential mediating mechanisms by which the variants of the RBFOX1 gene may exert an indirect and therefore small effect on aggression that becomes revealed in large samples. The strengths of the study are its standard of data collection performed in the uniform conditions of a laboratory, the strong representation of the regional population, and the solid rationale behind the selection of RBFOX1 as a target. Thus, the findings that RBFOX1 variants appear to be associated with personality traits and alcohol use disorder merit attention in further studies.

Data Availability Statement

The datasets generated for this study are available on request to the corresponding author.

Ethics Statement

The studies involving human participants were reviewed and approved by Ethics Review Committee on Human Research of the University of Tartu. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

NF-C, SF, BF, AR, BC, and JH contributed to the conception and design of the study. KL, TK, and JH selected and prepared the questionnaires, and KL conducted interviews. MV, KL, TK, MK, TV, and JH collected the data. MV, KL, and L-MT analyzed the data. MV performed genotyping. MV and JH wrote the manuscript. All authors were involved in the final editing and have approved the manuscript.

Funding

This work was supported by grants from the Estonian Research Council (IUT20-40) and European Community's Seventh Framework Programme (FP7/2007–2013) under grant agreement no 602805 (Aggressotype), European Union's Horizon 2020 research and innovation program under grant agreement no 667302 (CoCA), and European Union's Horizon 2020 research and innovation program under grant agreement no 728018 (Eat2beNICE).

Conflict of Interest

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

Acknowledgments

We are grateful to all participants of the ECPBHS and to the whole ECPBHS Study Team.

References

1. UN World Health Organization. Global Status Report on Violence Prevention. (2014). Available online at: http://www.refworld.org/docid/54aa8de14.html (accessed March 11, 2018).

Google Scholar

2. Georgiev AV, Klimczuk ACE, Traficonte DM, Maestripieri D. When violence pays: a cost-benefit analysis of aggressive behavior in animals and humans. Evol Psychol. (2013) 11:678–99. doi: 10.1177/147470491301100313

CrossRef Full Text | Google Scholar

3. Gómez JM, Verdú M, González-Megías A, Méndez M. The phylogenetic roots of human lethal violence. Nature. (2016) 538:233–7. doi: 10.1038/nature19758

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Veroude K, Zhang-James Y, Fernàndez -Castillo N, Bakker MJ, Cormand B, Faraone SV. Genetics of aggressive behavior: an overview. Am J Med Genet B Neuropsychiatr Genet. (2016) 171B:3–43. doi: 10.1002/ajmg.b.32364

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Asherson P, Cormand B. The genetics of aggression: where are we now? Am J Med Genet B. (2016) 171B:559–61. doi: 10.1002/ajmg.b.32450

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Vassos E, Collier DA, Fazel S. Systematic meta-analyses and field synopsis of genetic association studies of violence and aggression. Mol Psychiatry. (2014) 19:471–7. doi: 10.1038/mp.2013.31

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Pappa I, St Pourcain B, Benke K, Cavadino A, Hakulinen C, Nivard MG, et al. A genome-wide approach to children's aggressive behavior: the EAGLE consortium. Am J Med Genet Part B Neuropsychiatr Genet. (2016) 171:562–72. doi: 10.1002/ajmg.b.32333

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Fernàndez-Castillo N, Cormand B. Aggressive behavior in humans: Genes and pathways identified through association studies. Am J Med Genet B Neuropsychiatr Genet. (2016) 171:676–96. doi: 10.1002/ajmg.b.32419

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Davies MN, Verdi S, Burri A, Trzaskowski M, Lee M, Hettema JM, et al. Generalised anxiety disorder – a twin study of genetic architecture, genome-wide association and differential gene expression. PLoS ONE. (2015) 10:e0134865. doi: 10.1371/journal.pone.0134865

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Fogel BL, Wexler E, Wahnich A, Friedrich T, Vijayendran C, Gao F, et al. RBFOX1 regulates both splicing and transcriptional networks in human neuronal development. Hum Mol Genet. (2012) 21:4171–86. doi: 10.1093/hmg/dds240

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Lukong KE, Chang KW, Khandjian EW, Richard S. RNA-binding proteins in human genetic disease. Trends Genet. (2008) 24:416–25. doi: 10.1016/j.tig.2008.05.004

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Glisovic T, Bachorik JL, Yong J, Dreyfuss G. RNA-binding proteins and post-transcriptional gene regulation. FEBS Lett. (2008) 582:1977–86. doi: 10.1016/j.febslet.2008.03.004

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Jin Y, Suzuki H, Maegawa S, Endo H, Sugano S, Hashimoto K, et al. A vertebrate RNA-binding protein Fox-1 regulates tissue-specific splicing via the pentanucleotide GCAUG. EMBO J. (2003) 22:905–12. doi: 10.1093/emboj/cdg089

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Anantharaman V, Koonin EV, Aravind L. Comparative genomics and evolution of proteins involved in RNA metabolism. Nucleic Acids Res. (2002) 30:1427–64. doi: 10.1093/nar/30.7.1427

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Darnell JC, Richter JD. Cytoplasmic RNA-binding proteins and the control of complex brain function. Cold Spring Harb Perspect Biol. (2012) 4:a012344. doi: 10.1101/cshperspect.a012344

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Lee JA, Damianov A, Lin CH, Fontes M, Parikshak NN, Anderson ES, et al. Cytoplasmic Rbfox1 regulates the expression of synaptic and autism-related genes. Neuron. (2015) 89:113–28. doi: 10.1016/j.neuron.2015.11.025

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Hamada N, Ito H, Nishijo T, Iwamoto I, Morishita R, Tabata H, et al. Essential role of the nuclear isoform of RBFOX1, a candidate gene for autism spectrum disorders, in the brain development. Sci Rep. (2016) 6:30805. doi: 10.1038/srep30805

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Elia J, Gai X, Xie HM, Perin JC, Geiger E, Glessner JT, et al. Rare structural variants found in attention-deficit hyperactivity disorder are preferentially associated with neurodevelopmental genes. Mol Psychiatry. (2010) 15:637–46. doi: 10.1038/mp.2009.57

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Hamshere ML, Green EK, Jones IR, Jones L, Moskvina V, Kirov G, et al. Genetic utility of broadly defined bipolar schizoaffective disorder as a diagnostic concept. Br J Psychiatry. (2009) 195:23–9. doi: 10.1192/bjp.bp.108.061424

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Xu B, Roos JL, Levy S, van Rensburg EJ, Gogos JA, Karayiorgou M. Strong association of de novo copy number mutations with sporadic schizophrenia. Nat Genet. (2008) 40:880–5. doi: 10.1038/ng.162

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Bhalla K, Phillips HA, Crawford J, McKenzie OL, Mulley JC, Eyre H, et al. The de novo chromosome 16 translocations of two patients with abnormal phenotypes (mental retardation and epilepsy) disrupt the A2BP1 gene. J Hum Genet. (2004) 49:308–11. doi: 10.1007/s10038-004-0145-4

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Vounou M, Janousova E, Wolz R, Stein JL, Thompson PM, Rueckert D, et al. Sparse reduced-rank regression detects genetic associations with voxel-wise longitudinal phenotypes in Alzheimer's disease. Neuroimage. (2012) 60:700–16. doi: 10.1016/j.neuroimage.2011.12.029

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Fernàndez-Castillo N, Gan G, van Donkelaar MMJ, Vaht M, Weber H, Retz W, et al. RBFOX1, encoding a splicing regulator, is a candidate gene for aggressive behavior. Eur Neuropsychopharmacol. (2020) 30:44–55. doi: 10.1016/j.euroneuro.2017.11.012

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Sonuga-Barke EJS, Lasky-Su J, Neale BM, Oades R, Chen W, Franke B, et al. Does parental expressed emotion moderate genetic effects in ADHD? an exploration using a genome wide association scan. Am J Med Genet B Neuropsychiatr Genet. (2008) 147B:1359–68. doi: 10.1002/ajmg.b.30860

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Merjonen P, Keltikangas-Järvinen L, Jokela M, Seppälä I, Lyytikäinen LP, Pulkki-Råback L, et al. Hostility in adolescents and adults: a genome-wide association study of the Young Finns. Transl Psychiatry. (2011) 1:e11. doi: 10.1038/tp.2011.13

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Anney RJL, Lasky-Su J, Ó'Dúshláine C, Kenny E, Neale BM, Mulligan A, et al. Conduct disorder and ADHD: evaluation of conduct problems as a categorical and quantitative trait in the international multicenter ADHD genetics study. Am J Med Genet B Neuropsychiatr Genet. (2008) 147:1369–78. doi: 10.1002/ajmg.b.30871

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Young R, Sweeting H, West P. A longitudinal study of alcohol use and antisocial behaviour in young people. Alcohol Alcohol. (2008) 43:204–14. doi: 10.1093/alcalc/agm147

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Harro M, Eensoo D, Kiive E, Merenäkk L, Alep J, Oreland L, et al. Platelet monoamine oxidase in healthy 9- and 15-years old children: the effectof gender, smoking and puberty. Prog Neuropsychopharmacol Biol Psychiatry. (2001) 25:1497–511. doi: 10.1016/S0278-5846(01)00212-3

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Harro J, Laas D, Eensoo D, Kurrikoff T, Sakala K, Vaht M, et al. Orexin/hypocretin receptor gene (HCRTR1) variation is associated with aggressive behaviour. Neuropharmacology. (2019) 156:107527. doi: 10.1016/j.neuropharm.2019.02.009

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Laas K, Kiive E, Mäestu J, Vaht M, Veidebaum T, Harro J. Nice guys: homozygocity for the TPH2−703G/T (rs4570625) minor allele promotes low aggressiveness and low anxiety. J Affect Disord. (2017) 215:230–6. doi: 10.1016/j.jad.2017.03.045

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Buss AH, Perry M. The aggression questionnaire. J Pers Soc Psychol. (1992) 63:259–452. doi: 10.1037/0022-3514.63.3.452

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Coccaro EF, Berman ME, Kavoussi RJ. Assessment of life history of aggression: development and psychometric characteristics. Psychiatry Res. (1997) 73:147–57. doi: 10.1016/S0165-1781(97)00119-4

PubMed Abstract | CrossRef Full Text | Google Scholar

33. O'Leary A, Laas K, Vaht M, Kiive M, Veidebaum T, Reif A, et al. Nitric oxide synthase genotype interacts with stressful life events to increase aggression in male subjects in a population-representative sample. Eur Neuropsychopharmacol. (2020) 30:56–65. doi: 10.1016/j.euroneuro.2019.07.241

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Costa PT Jr, McCrae RR. The NEO-PI/NEO-FFI Manual Supplement. Odessa, FL: Psychological Assessment Resources (1989).

Google Scholar

35. Mõttus R, Pullmann H, Allik J. Toward more readable big five personality inventories. Eur J Psychol Assess. (2006) 22:149–57. doi: 10.1027/1015-5759.22.3.149

CrossRef Full Text | Google Scholar

36. Sheehan DV, Lecrubier Y, Sheehan KH, Amorim P, Janavs J, Weiller E, et al. The Mini-International Neuropsychiatric Interview (MINI): the development and validation of a structured diagnostic psychiatric interview for DSM-IV and ICD-10. J Clin Psychiatry. (1998) 59(Suppl. 20):22–33.

PubMed Abstract | Google Scholar

37. Shlik J, Aluoja A, Kihl E. MINI 5.0.0. Mini Rahvusvaheline Neuropsühhiaatriline Intervjuu DSM –IV. Estonian version of MINI international neuropsychiatric interview. Tartu (1999).

Google Scholar

38. Botvinick MM. Conflict monitoring and decision making: reconciling two perspectives on anterior cingulate function. Cogn Affect Behav Neurosci. (2007) 7:356–66. doi: 10.3758/CABN.7.4.356

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Devinsky O, Morrell MJ, Vogt BA. Contributions of anterior cingulate cortex to behaviour. Brain. (1995) 118 (Pt. 1):279–306. doi: 10.1093/brain/118.1.279

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Carretié L, Albert J, López-Martín S, Tapia M. Negative brain: an integrative review on the neural processes activated by unpleasant stimuli. Int J Psychophysiol. (2009) 71:57–63. doi: 10.1016/j.ijpsycho.2008.07.006

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Yang Y, Joshi SH, Jahanshad N, Thompson PM, Baker LA. Neural correlates of proactive and reactive aggression in adolescent twins. Aggress Behav. (2017) 43:230–40. doi: 10.1002/ab.21683

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Ducharme S, Hudziak JJ, Botteron KN, Ganjavi H, Lepage C, Collins DL, et al. Right anterior cingulate cortical thickness and bilateral striatal volume correlate with child behavior checklist aggressive behavior scores in healthy children. Biol Psychiatry. (2011) 70:283–90. doi: 10.1016/j.biopsych.2011.03.015

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Bufkin JL, Luttrell VR. Neuroimaging studies of aggressive and violent behavior: current fndings and implications for criminology and criminal justice. Trauma Violence Abuse. (2005) 6:176–91. doi: 10.1177/1524838005275089

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Cohn MD, Popma A, van den Brink W, Pape LE, Kindt M, van Domburgh L, et al. Fear conditioning, persistence of disruptive behavior and psychopathic traits: an fMRI study. Transl Psychiatry. (2013) 3:e319. doi: 10.1038/tp.2013.89

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Harro J, Marcusson J, Oreland L. Alterations in brain cholecystokinin receptors in suicide victims. Eur Neuropsychopharmacol. (1992) 2:57–63. doi: 10.1016/0924-977X(92)90037-9

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Caprara GV, Barbaranelli C, Zimbardo P. Understanding the complexity of human aggression: Affective, cognitive and social dimensions of individual differences in propensity toward aggression. Eur J Pers. (1996) 10:133–55. doi: 10.1002/(SICI)1099-0984(199606)10:2<133::AID-PER252>3.0.CO;2-E

CrossRef Full Text | Google Scholar

47. Tremblay PF, Ewart LA. the buss and perry aggression questionnaire and its relations to values, the Big Five, provoking hypothetical situations, alcohol consumption patterns, and alcohol expectancies. Pers Indiv Diff. (2005) 38:337–48. doi: 10.1016/j.paid.2004.04.012

CrossRef Full Text | Google Scholar

48. KodŽopeljić J, Smederevac S, Mitrović D, Dinić B, Colović P. School bullying in adolescence and personality traits: a person-centered approach. J Interpers Violence. (2014) 29:736–57. doi: 10.1177/0886260513505216

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Huesmann LR, Dubow EF, Boxer P. Continuity of aggression from childhood to early adulthood as a predictor of life outcomes: implications for the adolescent-limited and life-course-persistent models. Aggress Behav. (2009) 35:136–49. doi: 10.1002/ab.20300

PubMed Abstract | CrossRef Full Text | Google Scholar

50. McCrae RR, Costa PT. Handbook of personality: Theory and research. 2nd Edn. New York, NY: Guilford Press (1999).

Google Scholar

51. McCrae RR, John OP. An introduction to the five-factor model and its applications. J Pers. (1992) 60:175–215. doi: 10.1111/j.1467-6494.1992.tb00970.x

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Frost BP. On the relationship between extraversion and aggression. Psychol Rep. (1981) 49:1009–10. doi: 10.2466/pr0.1981.49.3.1009

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Galassi JP, Galassi MD. Relationship between assertiveness and aggressiveness. Psychol Rep. (1975) 36:352–4. doi: 10.2466/pr0.1975.36.2.352

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Cavalcanti JG, Pimentel CE. Personality and aggression: a contribution of the general aggression model. Estud Psicol. (2016) 33:443–51. doi: 10.1590/1982-02752016000300008

CrossRef Full Text | Google Scholar

Keywords: RBFOX1, A2BP1, aggressiveness, neuroticism, extraversion, alcohol use disorder, gender

Citation: Vaht M, Laas K, Fernàndez-Castillo N, Kurrikoff T, Kanarik M, Faraone SV, Tooding LM, Veidebaum T, Franke B, Reif A, Cormand B and Harro J (2020) Variants of the Aggression-Related RBFOX1 Gene in a Population Representative Birth Cohort Study: Aggressiveness, Personality, and Alcohol Use Disorder. Front. Psychiatry 11:501847. doi: 10.3389/fpsyt.2020.501847

Received: 30 September 2019; Accepted: 09 October 2020;
Published: 24 November 2020.

Edited by:

Ju Wang, Tianjin Medical University, China

Reviewed by:

Yoshihiko Matsumoto, Yamagata University, Japan
Yanling Zhou, Guangzhou Medical University, China

Copyright © 2020 Vaht, Laas, Fernàndez-Castillo, Kurrikoff, Kanarik, Faraone, Tooding, Veidebaum, Franke, Reif, Cormand and Harro. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Jaanus Harro, amFhbnVzLmhhcnJvJiN4MDAwNDA7dXQuZWU=

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