Abstract
The present study examined the association of dopamine-related genes with short- and long-term neurobehavioral recovery, as well as neurobehavioral recovery trajectories over time, in children who had sustained early childhood traumatic brain injuries (TBI) relative to children who had sustained orthopedic injuries (OI). Participants were recruited from a prospective, longitudinal study evaluating outcomes of children who sustained a TBI (n = 68) or OI (n = 72) between the ages of 3 and 7 years. Parents completed ratings of child executive function and behavior at the immediate post-acute period (0–3 months after injury); 6, 12, and 18 months after injury; and an average of 3.5 and 7 years after injury. Thirty-two single nucleotide polymorphisms (SNPs) in dopamine-related genes (dopamine receptor D2 [DRD2], solute carrier family 6 member 3 [SLC6A3], solute carrier family 18 member A2 [SLC18A2], catechol-o-methyltransferase [COMT], and ankyrin repeat and kinase domain containing 1 [ANKK1]) were examined in association with short- and long-term executive function and behavioral adjustment, as well as their trajectories over time. After controlling for premorbid child functioning, genetic variation within the SLC6A3 (rs464049 and rs460000) gene was differentially associated with neurobehavioral recovery trajectories over time following TBI relative to OI, with rs464049 surviving multiple testing corrections. In addition, genetic variation within the ANKK1 (rs1800497 and rs2734849) and SLC6A3 (rs464049, rs460000, and rs1042098) genes was differentially associated with short- and long-term neurobehavioral recovery following TBI, with rs460000 and rs464049 surviving multiple testing corrections. The findings provide preliminary evidence that genetic variation in genes involved in DRD2 expression and density (ANKK1) and dopamine transport (SLC6A3) plays a role in neurobehavioral recovery following pediatric TBI.
Introduction
T
The unexplained heterogeneity in recovery and outcomes following TBI is one of the most significant barriers to the development of effective prognostic tools and therapeutic interventions. 9,10 The most well-established predictor of outcome following TBI is injury severity, most often measured using the Glasgow Coma Scale (GCS). 11 Injury factors, in combination with demographic factors (e.g., age), however, typically account for only ∼35% of the variance in outcomes following moderate to severe TBI. 12 Family environmental factors appear to account for additional variance in outcomes. 13,14 Despite the known contributions of injury and environmental factors to recovery, substantial unexplained heterogeneity exists in neurobehavioral outcomes following pediatric TBI. Injuries of similar severity in similar environmental contexts can result in significantly disparate outcomes. Individual factors, such as genetic variability, may influence the pathophysiological response to injury, recovery, and repair processes, and response to intervention. 15 Therefore, investigation of the influence of genetic factors on recovery is a critical step for advancing our understanding of variability in neurobehavioral outcomes following pediatric TBI.
Genes modulating the dopamine systems are of particular interest among those likely to influence neurobehavioral recovery following early TBI. Dopamine is a neurotransmitter that plays a key role in cognition and behavior, especially in the executive functions, including working memory, cognitive control, cognitive flexibility, and planning, 16,17 as well as in related behaviors, including hyperactivity, impulse control, and reward sensitivity. 18 –20 Executive functions are mediated by the mesocortical and mesolimbic dopamine pathways, which originate in the ventral tegmental area of the midbrain and project to the prefrontal cortex, anterior cingulate cortex, anterior temporal structures, and the basal forebrain. 17,21 In the developing brain, dopamine plays a role in myelination, synaptogenesis, and pruning processes. 22 –25 Throughout childhood, the development of executive functions and related behaviors is associated with maturation of the dopamine system. 26 Dopamine, dopamine transporter, and D1 and D2 receptor levels increase in the striatum until middle childhood, with more protracted development of dopaminergic modulation of the prefrontal cortex occurring into adolescence. 27 –29
Modulation of various components of the dopamine system by dopamine-related genes can result in disrupted neurobehavioral functioning. The “inverted-U” hypothesis suggests that optimal functioning of frontostriatal circuits underlying executive functions and related behaviors relies on optimal levels of dopamine, whereby too much or too little results in impaired function. 17,30 Differences in executive functions and related behaviors are associated with variation in dopamine-related genes in typically developing children, 31 –33 as well as in conditions such as attention-deficit/hyperactivity disorder (ADHD), 34,35 oppositional defiant disorder, and conduct disorder. 19 Evidence from in vivo TBI models suggests that functional reductions in dopamine can occur after TBI. 36 Moreover, pharmacotherapies targeting dopamine system pathways are beneficial for post-injury neurobehavioral function in both adults and children; however, the specific mechanisms of TBI pathology targeted by these therapies, the most effective time course, and which individuals will reap the most benefit remain unclear. 36 Recent studies of TBI in adults have hypothesized that the effects of low dopamine on neurobehavioral functioning may be related to genetic variants that modulate dopamine levels. 37 –39
Several sources of variation within genes that modulate dopamine signaling may influence neurobehavioral recovery following pediatric TBI. Dopamine receptors are distributed at various densities throughout the frontostriatal dopamine system and receive the transmission of dopamine across the synapse. Genetic variability associated with dopamine receptor D2 (DRD2) expression and density has been examined by assessing variation within the DRD2 gene as well as the ankyrin repeat and kinase domain containing 1 (ANKK1) gene immediately upstream of DRD2. 39 Studies have shown that genetic variation within polymorphisms rs1800497, rs11604671, and rs4938016 in ANKK1 and rs6279 in DRD2 is associated with neuropsychological performance within the 1st year post-injury in adults. 38 –41 Less is known regarding the influence of other dopamine receptor genes on neurobehavioral function following TBI.
In addition to dopamine receptors, dopamine's effects are regulated through complex processes of release, reuptake, and metabolism. The vesicular monoamine transporter (VMAT) facilitates dopamine transport into synaptic vesicles for release. A recent study found that the rs363226 genotype within the vesicular monoamine transporter 2 (VMAT2) gene, also known as solute carrier family 18 member A2 (SLC18A2), was associated with cognitive impairment at 6 months following severe TBI in adults. 42 In most brain regions, dopamine is removed from the synapse primarily by the dopamine transporter (DAT). Several studies have implicated a variable number tandem repeat (VNTR rs28363170) within the dopamine transporter gene, solute carrier family 6 member 3 (SLC6A3), in influencing both short-term 43 –45 and long-term 37 dopamine-relevant TBI pathology; however, a study in adults with severe TBI failed to find a significant association with composite cognitive scores at 6 or 12 months post-injury. 38
In the prefrontal cortex, dopamine transporter expression is lower and dopamine clearance is mediated by catechol-o-methyltransferase (COMT), an enzyme that metabolizes dopamine. 46,47 The COMT gene has received the most research attention in TBI, but has produced mixed results. Variation within the Val158Met (rs4680) single nucleotide polymorphism (SNP) has been associated with poorer neuropsychological performance in a cohort of adults with mild-severe TBI assessed within 1 year post-injury, 48 and in a cohort of adults with mild TBI examined 6 months post-injury; 49 however, the COMT gene showed no significant associations with post-injury performance in a cohort of adults with moderate-severe TBI relative to healthy controls assessed within 3 months post-injury, 50 or in the largest investigated cohort of adults with moderate-severe TBI assessed 1 month post-injury. 51 In the only prior study examining a dopamine-related gene in association with neurobehavioral function following TBI in childhood, variation within rs4680 was associated with poorer neuropsychological function across groups at 18 months following mild-severe TBI or orthopedic injury (OI). 52
Although initial investigations of the influence of dopamine-related genes on neurobehavioral function following TBI are promising, the present study addresses three important limitations of the existing literature. First, few studies have included a non-brain-injured comparison group for examining the differential association of dopamine-related genes with neurobehavioral recovery following injury to the brain. Without examination of the moderating effect of brain injury, previous studies have been limited in their ability to differentiate effects of dopamine genes on neurobehavioral recovery after TBI from differences that may reflect premorbid genetic influences on neurobehavioral function, as have been found in individuals without brain injuries. 49 Second, no studies have examined the association of dopamine-related genes with neurobehavioral outcomes beyond 18 months post-injury. The effects of dopamine on neurobehavioral function may vary over time post-injury, with different pathophysiological processes occurring during the acute and chronic phases of recovery and influencing recovery trajectories over time. Third, the relation of dopamine-related genes to neurobehavioral function following TBI has almost exclusively been studied in adults. Given the important role that dopamine systems play in brain development and development of executive functions and related behaviors throughout childhood and into adolescence, dopamine-related genes may play an even more prominent role in neurobehavioral recovery following TBI sustained during childhood.
Therefore, the primary aim of the present study was to examine the differential effects of dopamine-related genes on neurobehavioral recovery trajectories over time in children who sustained early childhood TBI relative to children who sustained OI. In addition, because we might expect different genetic effects during the short- and long-term phases of recovery, we also had an exploratory aim of examining the differential effects of dopamine-related genes on short-term (∼ 6 months post-injury) and long-term (∼7 years post-injury) neurobehavioral recovery. We hypothesized that variation within dopamine-related genes would be associated with changes in executive function and behavioral adjustment over time post-injury, as well as in the short and long term, and that these associations would be moderated by injury type (e.g., TBI vs. OI).
Materials and Methods
Participants
Participants were part of a prospective, longitudinal study evaluating outcomes of children who sustained a TBI or OI between the ages of 3 and 7 years. 5,6,53 They were recruited from three tertiary care children's hospitals and one tertiary care, general hospital in Ohio. Participants completed assessments at multiple visits, including the immediate post-acute period (0–3 months after injury); 6 (referred to as “short-term” throughout), 12, and 18 months after injury; and an average of 3.5 and 7 years after injury (referred to as “long-term” throughout).
Inclusion criteria included hospitalization overnight for traumatic injury (TBI or OI), no evidence of child abuse as the cause of the injury, no history of documented neurological problems or developmental delays pre-injury, and English as the primary language spoken in the home. Severity of TBI was characterized using the lowest post-resuscitation GCS score. 11 Severe TBI was defined as a GCS score ≤8. Moderate TBI was defined as a GCS score of 9–12 with or without abnormal neuroimaging (moderate TBI) or a higher GCS score with abnormal neuroimaging as defined by an intracranial or parenchymal injury or depressed skull fracture (complicated mild TBI). Mild TBI was defined as a GCS score ≥13 without abnormal neuroimaging. The OI group included children who had sustained a bone fracture (not including skull fractures), had an overnight stay in the hospital, and did not exhibit alterations in consciousness or other signs or symptoms of head trauma or brain injury. The study was approved by the institutional review boards at each of the participating medical centers and informed consent was obtained from participating caregivers.
DNA collection and genotyping
DNA was collected from saliva samples and purified using the United States Food and Drug Administration (FDA)-cleared Oragene OG-500 self-collection tubes (DNA Genotek, Ottawa, Canada). Oragene allows the collection, stabilization, and long-term storage of DNA from saliva at ambient temperature. Salivary DNA is a valid and reliable alternative to blood DNA for high-throughput genotyping.
54
The HumanExome v1.1 Bead Chip (Illumina, San Diego, CA) was used to perform genotyping using the Illumina iScan system to identify SNPs associated with the following dopamine-related genes: DRD2, SLC6A3, SLC18A2, COMT, and ANKK1. SNPs were selected within genes' boundaries plus 5 kilo bases upstream and downstream from the genes. These candidate genes were chosen based on their associations with outcomes following TBI in adults and their biological relevance in dopamine pathways in adults and children as described. Prior to analysis, quality of SNP calls from the chip was evaluated. SNPs that failed Hardy–Weinberg equilibrium (p < 0.0001) or had minor allele frequencies <10% were excluded. Thresholds for quality control for call rates at individual and SNP levels were 99% and 90%, respectively. Cryptic relatedness was checked using Graphical Representation of Relationships (GRR)
55
(
Outcome measures
Parents completed the age-appropriate form of the Behavior Rating Inventory of Executive Function (BRIEF) at all visits, the first of which was based on the child's behavior prior to injury, used in analysis to control for premorbid differences in child functioning. The BRIEF is a parent-report measure of child executive function as evident in everyday behavior. 56 –58 We analyzed the age-standardized global executive composite T score (BRIEF GEC) to assess global executive function behaviors. Higher scores reflect poorer executive functioning.
Parents completed the age-appropriate form of the Child Behavior Checklist (CBCL) 59 at all visits, the first of which was based on retrospective recall of the child's behavior prior to injury to control for premorbid differences in child functioning. The CBCL is a parent-report measure of child behavioral adjustment and possesses high test–retest reliability and criterion-related validity. We analyzed the age- and sex-standardized Total Problems T score to assess child behavioral adjustment. Higher scores reflect poorer behavioral adjustment.
Both outcome measures are National Institute of Neurological Disorders and Stroke (NINDS) recommended common data elements for pediatric TBI and are well validated for the pediatric TBI population. 60
Statistical analysis
Statistical analyses were conducted using SAS 9.4 (SAS, Cary, NC). Prior to analyses, child outcome data were reviewed for plausibility according to the following rules to reduce the potential influence of outliers: 1) participants with changes in outcome scores between the short-term and long-term time point that exceeded three standard deviations of change were excluded from respective analyses (n = 2 from BRIEF analysis; n = 0 from CBCL analysis); and 2) outcome scores of T > 90 (4 SDs above the normative mean) were winsorized to 90 (n = 5 from BRIEF analysis; n = 1 from CBCL analysis).
Our primary analysis examined genetic associations with outcome trajectories over time, using a random coefficients longitudinal mixed model including all time points. Using longitudinal analyses, we examined the association of each SNP with executive function and behavioral adjustment recovery trajectories over time and the moderating effect of injury group (TBI vs. OI) on these associations (e.g., SNP × injury group × time interaction). Mixed models included random intercepts, slopes, and subjects. Random slopes were modeled over time. Our exploratory analysis used multiple linear models to examine genetic associations with cross-sectional outcomes in the short-term (6 months post-injury; 0.62 ± 0.08 years) and long-term (6.82 ± 1.12 years post-injury) phases of recovery, corresponding to the second and sixth study visits, respectively. We employed multiple linear models to examine the association of each SNP with short-term and long-term executive functioning and behavioral adjustment, and the moderating effect of injury group (TBI vs. OI) on these associations (e.g., SNP × injury group interaction). In all models, higher level interactions were examined first, followed by each lower level interaction if nonsignificant at a nominal (p < 0.05) level. In all association tests, we used an additive genetic model in which major homozygotes were coded as 0, heterozygotes were coded as 1, and minor homozygotes were coded as 2.
Covariates initially included and then trimmed if nonsignificant were the child's pre-injury functioning on the measure of interest and socioeconomic status (SES) (defined by averaging sample z scores for maternal education and median income). 53 Continental ancestry was covaried in all models, dichotomized as European versus African, to account for potential confounding related to race. Results were evaluated after correcting for multiple testing using a Bonferroni adjustment taking into account average pairwise LD correlation between all SNPs considered in the analyses (r = 0.16) as a way to correct for the number of independent tests. 55 The results of the SNP times injury group interaction for each SNP are provided in Manhattan plots with both nominal (p < 0.05) and multiple testing-adjusted (p = 0.0027) significance thresholds indicated. For each interaction that was at least nominally significant, we plotted regression lines for each group by allele status to examine the direction of the interaction and identify the risk allele. When plotting the data, SNPs with <10 participants for minor allele homozygotes were combined with heterozygotes to ensure more robust estimates.
Results
Sample description
Of the 221 participants enrolled in the original study, 141 provided DNA samples. Participants with genetic data did not differ significantly from those without genetic data in demographic characteristics or on various study measures (see Table 1). Of those with genetic data, 140 of 141 participants had covariate and outcome data for at least one outcome visit to be included in the present analyses. Of those included in the present analyses, 10 had mild TBI, 42 had complicated mild to moderate TBI, 16 had severe TBI, and 72 had OI. The TBI and OI groups did not differ significantly in race, sex, age at injury, age at assessment, or SES (Table 2). As shown in Table 2, children with TBI had poorer executive functioning and poorer behavioral adjustment relative to children with OI during both the short-term and long-term phases of recovery. These results are consistent with our prior reports based on the larger sample from this cohort. 61,62
BRIEF, Behavior Rating Inventory of Executive Function; CBCL, Child Behavior Checklist; GCS, Glasgow Coma Scale.
Analysis of covariance (ANCOVA) controlling for premorbid child functioning on each measure and zSES.
BRIEF GEC, Behavior Rating Inventory of Executive Function; CBCL, Child Behavior Checklist; GCS, Glasgow Coma Scale; LSM, least squares mean; OI, orthopedic injury; TBI, traumatic brain injury; zSES, socioeconomic status z score.
SNPs investigated
Thirty-two SNPs across the five dopamine-related genes met quality control and inclusion criteria. The majority of SNPs (30 out of 32) were located within boundaries of selected genes, and two SNPs were located within the 5 kilo bases interval downstream from DRD2 and COMT genes. Table 3 contains descriptive information for each SNP. LD plots were generated for the ANKK1 and SLC6A3 genes based on their associations with study outcomes (see Fig. 1).

Linkage disequilibrium (LD) plots for single nucleotide polymorphisms (SNPs) examined in the ankyrin repeat and kinase domain containing 1 (ANKK1) (top) and solute carrier family 6 member 3 (SLC6A3) (bottom) genes.
Single nucleotide polymorphism (SNP) is outside of gene region, within 5 kilo bases of extension downstream.
SLC6A3, solute carrier family 6 member 3; SLC18A2, solute carrier family 18 member A2; ANKK1, ankyrin repeat and kinase domain containing 1; DRD2, dopamine receptor D2; COMT, catechol-o-methyltransferase.
Associations of dopamine genes with longitudinal recovery trajectories
Figure 2 illustrates the SNP by injury group interactions that met nominal and multiple testing-adjusted significance thresholds for their effect on longitudinal recovery trajectories of executive functioning and behavioral adjustment. The three-way interaction of SNP by injury group by time was not significant for either child outcome, suggesting that differences in recovery trajectories among groups remained consistent across time since injury. Lower level analyses revealed significant interactions of SNP by injury group across time for rs464049 and rs460000 SNPs in SLC6A3 on both executive functioning (p = 0.005 and p = 0.005) and behavioral adjustment (p = 0.002 and p = 0.019) recovery trajectories. The interactions of SNP by injury group across time for rs464049 SNP in SLC6A3 on behavioral adjustment remained significant after multiple testing corrections.

Manhattan plots for all single nucleotide polymorphisms (SNPs) examined in association with longitudinal recovery trajectories of executive function and behavioral adjustment. Triangles represent SNPs located in 5 Kb extensions from dopamine receptor D2 (DRD2) and catechol-o-methyltransferase (COMT). The dashed line indicates nominal significance at p < 0.05. The dotted line indicates multiple testing-adjusted significance (p = 0.0027), taking into account linkage disequilibrium (LD) correlation between SNPs.
Associations of dopamine genes with short- and long-term recovery
Figure 3 illustrates the SNP by injury group interactions that met nominal and multiple testing significance thresholds for their effect on short- and long-term executive functioning and behavioral adjustment. The rs464049 and rs460000 SNPs in SLC6A3 exhibited interactions with injury group on both short- (p = 0.010 and p < 0.001) and long-term (p = 0.025 and p = 0.014) executive functioning. The rs1800497 SNP (p = 0.029) in ANKK1 also exhibited interactions with injury group on short-term executive functioning. The rs464049 SNP in SLC6A3 also exhibited interactions with injury group for both short- (p = 0.015) and long-term (p = 0.002) behavioral adjustment. In addition, the rs2734849 SNP (p = 0.040) in ANKK1 exhibited an interaction with injury group on short-term behavioral adjustment. Finally, the rs1042098 SNP (p = 0.027) in SLC6A3 exhibited an interaction with injury group on long-term behavioral adjustment. The interactions of SNP by injury group for rs460000 in SLC6A3 on short-term executive functioning and for rs464049 SNP in SLC6A3 on long-term behavioral adjustment remained significant after multiple testing corrections.

Manhattan plots for all single nucleotide polymorphisms (SNPs) examined in association with short- and long-term executive function and behavioral adjustment. Triangles represent SNPs located in 5 Kb extensions from dopamine receptor D2 (DRD2) and catechol-o-methyltransferase (COMT). The dashed line indicates nominal significance at p < 0.05. The dotted line indicates multiple testing-adjusted significance (p = 0.0027), taking into account linkage disequilibrium (LD) correlation between SNPs.
Identification of risk alleles
Given the presence of significant SNP by injury group interactions, we performed stratified, descriptive analyses to understand the direction of the association of the variants with outcomes within each injury group and to identify the risk allele. Figures 4 and 5 illustrate the direction of associations of genetic variation in SLC6A3 and ANKK1 with executive function (Fig. 4) and behavioral adjustment (Fig. 5) outcomes within each injury group at the short- and long-term phases of recovery, as well as on longitudinal recovery trajectories. Homozygotes for the minor allele were combined with heterozygotes because of small cell sizes for rs1042098, rs460000, and rs1800497. Results revealed the following risk alleles: G for rs1042098, G for rs464049, A for rs460000, A for rs2734849, and A for rs1800497.

Regression lines for each injury group by allele status plotted to examine the direction of the interaction and identify the risk allele. An additive model was used to plot regression lines for all single nucleotide polymorphisms (SNPs). SNPs with <10 participants for minor allele homozygotes were combined with heterozygotes.

Regression lines for each injury group by allele status plotted to examine the direction of the interaction and identify the risk allele. An additive model was used to plot regression lines for all single nucleotide polymorphisms (SNPs). SNPs with <10 participants for minor allele homozygotes were combined with heterozygotes.
Discussion
The present study sought to examine the association of dopamine-related genes with neurobehavioral recovery trajectories over time, as well as short- and long-term neurobehavioral recovery, in children who had sustained early childhood TBI relative to children with OI. Genetic variation within the SLC6A3 (rs464049 and rs460000) gene was differentially associated with neurobehavioral recovery trajectories over time following TBI relative to OI. In addition, genetic variation within ANKK1 (rs1800497 and rs2734849) and SLC6A3 (rs464049, rs460000, and rs1042098) was differentially associated with short- and long-term neurobehavioral recovery following TBI. The findings provide preliminary evidence that genetic variation in genes involved in DRD2 expression and density (ANKK1) and dopamine transport (SLC6A3) plays a role in neurobehavioral recovery following pediatric TBI.
Our most robust and consistent findings were for SNPs within the SLC6A3 gene that encodes expression of the dopamine transporter. The rs1042098 SNP exhibited an interaction with injury group on long-term behavioral adjustment, and the rs464049 and rs460000 SNPs exhibited interactions with injury group across time on both executive functioning and behavioral adjustment recovery trajectories. In the cross-sectional analyses, the rs464049 and rs460000 SNPs exhibited interactions with injury group on both short- and long-term executive functioning, and the rs464049 SNP exhibited interactions with injury group for both short- and long-term behavioral adjustment. Risk alleles were A for rs460000, G for rs464049, and G for rs1042098. These results provide the first evidence for an association of variation within the SLC6A3 gene with clinical outcomes following TBI.
Although little is known about the functional impact of the significant SNPs associated with recovery in the present findings, other polymorphisms within this gene have shown associations with dopamine transporter expression and with pathology in adult TBI. Experimental TBI studies have demonstrated reduced maximal evoked dopamine overflow, altered dopamine clearance kinetics, and dopamine transporter protein expression changes post-injury. 63 Prior studies in adult TBI have shown reduced striatal dopamine transporter binding in association with the VNTR rs28363170 in SLC6A3, 37 although examination of clinical outcomes failed to find a significant association of this polymorphism with composite cognitive scores at 6 or 12 months following severe TBI. 38,39,64 As we did not directly assess this VNTR variant, it is not possible to determine if our variants are contributing independently to TBI recovery. Nevertheless, the significance of SLC6A3 for both trajectories of neurobehavioral recovery over time and long-term recovery in the present study might suggest that genetic variation in dopamine transporter expression may underlie the chronic pathology associated with persistent neurobehavioral impairment years after early childhood TBI.
Our exploratory cross-sectional analyses revealed additional differential associations with outcomes following TBI involving the ANKK1 gene. The rs1800497 SNP in ANKK1 exhibited interactions with injury group on short-term executive functioning, with post-hoc analyses showing poorest executive functioning in children with TBI who were carriers of the A allele, referred to as A1 in prior literature. Although this association failed to survive multiple comparison corrections, confidence in the findings is bolstered by its consistency with prior research in adult TBI. Several studies have shown associations of rs1800497 with cognitive outcomes following TBI in adulthood, 38 –41 although findings have been mixed with regard to the risk allele or genotype most closely associated with these outcomes. 39,64 Because our study included a non-brain-injured control group, our findings not only generalize the involvement of the ANKK1 gene to outcomes following TBI in childhood, but also lend evidence for a specific role of this gene in recovery from injury to the brain. The results for ANKK1 rs1800497 suggest a role for DRD2 expression and density in neurobehavioral recovery following TBI in childhood. Studies of healthy individuals suggest that A1 carriers show lower receptor densities, reduced receptor binding, and increased dopamine synthesis rate, especially in the striatum. 65 –67 In view of the mixed findings on the role of the rs1800497 in neurobehavioral outcomes of adult TBI, further study is needed to clarify the effects of genetic variation within ANKK1 on dopamine neurotransmission following brain injury.
Also within ANKK1, our results showed an interaction between SNP rs2734849 and injury group on short-term behavioral adjustment, with post-hoc analyses indicating poorest behavioral adjustment in children with TBI who were homozygous for the A allele. In contrast to the previous evidence for a role of rs1800497 in outcome following TBI, no previous studies have implicated rs2734849 in recovery from TBI. Moreover, little is known regarding the role of this SNP aside from possible involvement in dyskinesia in Parkinson's disease 68 and in nicotine dependence. 69 This SNP may serve as a new target of investigation in future studies of pediatric TBI.
Given that our exploratory cross-sectional analyses identified an additional gene (ANKK1) that was not significantly associated with outcomes in our longitudinal models, these findings might suggest that genetic variation related to DRD2 expression and density is more salient in determining neurobehavioral outcomes in the first several months post-injury, whereas genetic variation related to dopamine transport (SLC6A3) has more lasting influences on neurobehavioral recovery over time. Our longitudinal mixed models may have been underpowered to detect a three-way interaction among our ANKK1 SNPs, injury group, and time. Viewed from a more conservative perspective, however, it is possible that our associations of ANKK1 with short-term outcomes reflect type I error, as these findings were not identified in the more conservative longitudinal mixed models. Studies with larger samples are needed to replicate the present results.
The findings of significant associations of dopamine-related genes with neurobehavioral recovery following pediatric TBI may add to a conceptual model that incorporates the various risk and protective factors contributing to neurobehavioral phenotypes following pediatric TBI. In addition to being implicated in recovery from TBI, the SLC6A3 gene and other dopamine genes are also implicated in risk for ADHD 70,71 and related endophenotypes, including poor inhibitory control, 72 mood instability, 73 and reduced working memory. 74 The substantial overlap of such traits in children who have sustained TBI and in those with ADHD, combined with potential commonalities in genetic risk factors, suggests that shared pathological mechanisms may contribute to shared clinical presentations. 71 This hypothesis is supported by the frequent diagnosis of secondary ADHD in children with TBI, 75,76 as well as evidence for poorer recovery and exacerbation of symptoms following TBI in children with premorbid neurobehavioral difficulties. 77 –80 Conceptualizing traumatic injury to the brain as an environmental event that interacts with an individual's genetic predisposition, outcome, and recovery following TBI in childhood would be consistent with a diathesis-stress model in which children with genetic vulnerabilities for dopamine-related neurobehavioral difficulties experience poorer outcomes because of an emergence or exacerbation of these problems following traumatic disruption to the dopamine system.
The other side of the risk perspective discussed is the potential protective effects that certain genotypes may confer following TBI. Interestingly, cognitive neuroscience research into learning and cognitive training has also implicated roles for both striatal dopamine D2 receptor density/binding and dopamine transporter expression, as well as their genetic determinants. Dopamine may be important for plasticity by enhancing neural sprouting and synaptogenesis. 81 Positron emission tomography (PET) studies have shown changes in striatal D2 binding following working memory training, 82 and lower D2 densities have been associated with better motor sequence learning 83 in healthy adults. Recently, children and adolescents who were T/A1 carriers of ANKK1 rs1800497 demonstrated greater improvements during working memory training. 84 Similarly, variation in SLC6A3 influences improvements in working memory following cognitive training in preschool and school-age children. 85,86 Taken together, these studies suggest that the ANKK1 and SLC6A3 genes may confer higher cognitive or neural plasticity 84,87 that may be especially relevant to recovery from injury to the brain. Clinically, individuals with certain genotypes may be more responsive to cognitive or medical interventions following TBI. Identifying those who are more or less likely to respond to specific interventions is critical in treatment selection and optimizing recovery.
The present findings should be considered in light of several study limitations. Although the present study is one of the largest investigations of the association of genetic factors after TBI in childhood, larger sample sizes are needed to replicate the findings. Because the study population consisted primarily of children with moderate TBI and all injuries were sustained during early childhood, sample size limitations precluded investigation of potential differences in genetic effects in relation to injury severity or age at injury. Another significant limitation is that measurement of child outcomes relied entirely on parent report, which may be subject to bias. Nevertheless, parent reports, especially of executive function behaviors, may be more ecologically valid than paper-and-pencil tests administered in a controlled environment. 57 Future studies that also include objective, performance-based measures of neuropsychological functioning, however, would be informative. Related, our reliance on retrospective recall of premorbid child functioning is another limitation. Given our focus on dopamine-related genes, it would have also been informative to know how many children were taking dopaminergic agents during their recovery and whether these medications influenced the results. Unfortunately, however, sufficiently reliable and valid data on medication use was unavailable for analysis. Finally, the exome genotyping chip that we used did not capture all common variation in our genes of interest, hampering the identification of the most promising variants for future follow- up. Because our candidate gene approach evaluated only preselected dopamine genes, we did not explore the potential role of other polymorphisms or gene–gene interactions. Given the known moderating role of the family environment on recovery following TBI, 53,88 examination of gene–environment interactions in relation to dopamine genes is underway.
This work aimed to advance understanding of unexplained heterogeneity in neurobehavioral outcomes following pediatric TBI through the identification of genetically influenced biological processes affecting recovery. The present study is the first to provide preliminary evidence for the influence of variation in dopamine-related genes on neurobehavioral recovery following TBI in children. Our inclusion of the OI group was a significant strength, allowing us to identify genetic effects related specifically to recovery from brain injury. The fact that genetic associations were specific to the TBI group and were evident even when controlling for premorbid child functioning further increases our confidence that these associations are related to recovery processes rather than reflecting premorbid expression of genetic differences. The longitudinal study design also permitted examination of genetic effects on recovery trajectories over time as well as on both short- and long-term recovery. Elucidation of biological processes underlying neurobehavioral recovery following pediatric TBI has the potential to substantially advance the field by facilitating more individualized prognosis, the development of novel targets for intervention, and a reduction in public health burden.
Footnotes
Acknowledgments
Funding for this study was supported in part by the Rehabilitation Medicine Scientist Training Program (RMSTP) K-12 HD001097-16, National Institute for Child Health and Human Development K23HD074683-01A1, R01 HD42729, and Grant 8 UL1 TR000077 from the National Center for Advancing Translational Sciences (NCATS) of the National Institutes of Health (NIH). The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH or other supporting agencies.
Author Disclosure Statement
No competing financial interests exist.
