Spindoc

Models

All models for this gene with cis-heritability P-value <0.01 were tested and are shown here. In each tissue, a gene can have zero, one, or multiple RNA phenotypes for each RNA modality. After the heritability statistics, subsequent columns show cross-validation performance (R2 and P-value) of each modeling method. The best-performing method for each models is used for TWAS.

# Tissue Modality RNA phenotype h2 h2 SE h2 P xQTL R2 BLUP R2 ENET R2 LASSO R2 xQTL P BLUP P ENET P LASSO P
1 pVTA isoform ratio NM_001399680.1 0.05 0.03 5.0e-03 0 0.01 -0 0 1.6e-01 9.2e-02 6.5e-01 1.7e-01
2 Adipose mRNA stability Spindoc 0.03 0.02 3.6e-03 0.02 0.01 0 0.01 1.7e-03 1.7e-02 8.8e-02 1.3e-02
3 BLA gene expression Spindoc 0.07 0.05 8.4e-03 0.01 0 0.01 0.01 8.7e-02 2.9e-01 1.1e-01 1.1e-01
4 NAcc mRNA stability Spindoc 0.03 0.02 3.7e-04 0.02 0.03 0.01 0.03 9.4e-04 5.8e-05 6.3e-03 6.1e-05

Trait associations

The numbered columns correspond to the numbered models in the table above. They contain the TWAS Z-score for each trait-model pair. Colors indicate strength of Z-scores: cells with lower negative Z-scores are darker blue, cells with higher positive Z-scores are darker red. Significant associations have white and bolded text.

Trait Avg chi2 ratio Avg chi2 Max chi2 1 2 3 4
Bone: apparent density1.32.32.71.61.51.5-1.3
Area under the delay curve0.30.30.60.40.60.5-0.8
BMI with tail6.121.222.94.64.54.5-4.8
BMI without tail1.23.85.52.11.61.7-2.3
Body weight3.917.223.64.73.33.5-4.9
Weight adjusted by age0.71.12.21.50.70.9-1.0
Bone surface4.68.514.63.82.62.5-2.4
Bone volume3.05.37.62.82.42.3-1.7
Breakpoint value in progressive ratio session0.91.13.31.80.10.1-0.9
Number of licking bursts0.10.10.30.50.10.3-0.2
Context. condit. distance diff. score0.20.20.6-0.80.30.40.1
Locomotion velocity, session 10.00.00.1-0.00.1-0.1-0.3
Locomotion distance, session 11.01.31.80.91.21.0-1.3
Locomotion velocity, session 21.71.93.01.71.00.8-1.7
Locomotion distance, session 20.20.30.70.50.30.0-0.8
Locomotion velocity, session 32.93.35.10.72.32.1-1.8
Locomotion distance, session 30.20.20.50.10.60.3-0.7
Stereotopy head waving bouts, day 31.31.52.70.21.51.7-1.0
Stereotopy head waving duration, day 30.30.40.9-0.20.80.9-0.2
Locomotion velocity, session 70.80.92.30.11.51.0-0.7
Locomotion distance, session 70.30.40.90.00.90.5-0.8
Stereotopy head waving bouts, day 70.50.60.80.80.90.9-0.2
Stereotopy head waving duration, day 71.01.22.21.50.90.9-0.8
Locomotion velocity, session 80.30.41.00.3-0.5-1.0-0.0
Locomotion distance, session 80.40.61.20.21.10.9-0.4
Degree of sensitization distance1.01.01.50.51.21.2-0.9
Degree of sensitization stereotypy0.80.81.61.30.90.8-0.5
Conditioned locomotion0.50.51.11.00.20.4-0.9
Cocaine response after cond. corrected0.30.31.1-1.1-0.1-0.30.2
Cocaine response after cond. not corrected0.10.10.5-0.70.2-0.10.0
Cocaine response before conditioning0.30.30.91.00.50.2-0.2
Saline control response1.62.14.12.00.80.3-1.9
Bone: connectivity density3.14.97.12.72.02.0-2.1
Bone: cortical apparent density0.60.81.7-0.5-1.3-1.00.6
Bone: cortical area1.21.82.21.51.31.2-1.3
Bone: cortical porosity0.50.61.1-0.2-0.8-1.10.7
Bone: cortical porosity0.30.30.90.2-0.6-0.90.3
Bone: cortical thickness0.50.91.51.0-0.6-0.7-1.2
Bone: cortical thickness0.30.30.6-0.1-0.5-0.8-0.6
Bone: cortical tissue density0.60.81.9-0.4-1.4-1.00.5
Condit. Reinf. active minus inactive responses0.30.30.50.60.40.6-0.7
Condit. Reinf. active-inactive response ratio0.10.10.3-0.1-0.5-0.40.2
Condit. Reinf. active responses0.81.01.70.71.01.3-0.9
Condit. Reinf. inactive responses0.91.12.20.51.21.5-0.6
Incentive salience index mean1.11.41.81.01.01.2-1.3
Condit. Reinf. lever presses1.72.32.61.41.61.6-1.5
Condit. Reinf. active minus inactive responses0.10.10.30.10.60.10.5
Condit. Reinf. active-inactive response ratio0.60.71.20.21.11.10.5
Condit. Reinf. active responses0.10.10.30.30.50.0-0.1
Condit. Reinf. inactive responses0.20.30.90.1-0.3-1.0-0.5
Condit. Reinf. lever presses1.62.23.81.91.51.1-1.2
Condit. Reinf. lever reinforcers received0.10.20.40.60.50.1-0.2
Conditioned reinforcement - actives0.10.20.50.50.30.20.7
Conditioned reinforcement - inactives0.50.71.5-0.31.21.10.5
Conditioned locomotion1.62.05.3-2.30.30.41.6
Patch foraging indifference point 0 sec0.20.30.50.30.70.5-0.6
Patch foraging indifference point AUC0.60.92.81.70.40.3-0.6
Delay disc. indifference point, 0s delay0.20.20.30.10.40.5-0.5
Delay disc. indifference point, 16s delay0.50.51.2-0.11.10.9-0.4
Delay disc. indifference point, 2s delay0.30.30.50.40.70.7-0.3
Delay disc. indifference point, 24s delay0.10.10.2-0.10.20.1-0.5
Delay disc. indifference point, 4s delay0.40.41.20.20.40.4-1.1
Delay disc. indifference point, 8s delay0.10.10.4-0.2-0.4-0.60.0
Indifference point function ln k1.83.29.6-3.1-0.9-0.81.3
Indifference point function log k1.83.29.6-3.1-0.9-0.81.3
Patch foraging total patch changes 0 sec0.00.10.3-0.00.30.1-0.5
Patch foraging total patch changes 12 sec2.22.95.72.41.10.6-2.1
Patch foraging total patch changes 18 sec0.20.31.01.00.30.1-0.4
Patch foraging total patch changes 24 sec2.44.16.42.11.91.4-2.5
Patch foraging total patch changes 6 sec1.01.52.41.31.10.8-1.5
Patch foraging time to switch 0 sec1.73.64.82.22.01.8-1.4
Patch foraging water rate 0 sec1.42.96.7-2.6-1.3-1.11.4
Patch foraging water rate 12 sec2.75.713.9-3.7-1.7-1.61.8
Patch foraging water rate 18 sec0.61.03.5-1.9-0.6-0.20.2
Patch foraging water rate 24 sec1.62.24.8-2.2-1.0-1.11.3
Patch foraging water rate 6 sec2.13.18.0-2.8-1.3-1.31.1
Seeking ratio, delayed vs. immediate footshock0.10.10.2-0.30.40.4-0.4
Distance traveled before self-admin1.92.52.91.71.51.6-1.4
Distance traveled after self-admin0.00.10.10.3-0.3-0.1-0.1
Extensor digitorum longus weight4.316.829.75.42.52.6-4.9
Bone: elastic displacement1.82.02.9-1.7-1.5-1.40.9
Bone: elastic work0.30.41.40.4-0.5-0.2-1.2
Bone: endosteal estimation1.63.67.30.82.42.7-0.8
Bone: endosteal perimeter1.73.57.10.82.52.7-0.3
Epididymis fat weight0.62.03.8-1.01.71.90.4
Delta time in closed arm before/after self-admin0.50.51.9-0.10.30.3-1.4
Delta time in open arm before/after self-admin0.90.91.51.21.00.9-0.6
Diff in mean of infusions in LGA sessions0.30.41.10.60.1-0.1-1.0
Delay discount exponential model param0.30.30.6-0.4-0.6-0.50.8
Extinction: sum of active levers before priming2.22.43.71.91.11.2-1.8
Active lever presses in extinction session 61.81.84.62.10.71.0-1.1
Bone: final force2.64.26.32.51.71.8-2.1
Bone: final moment2.64.26.12.51.71.8-2.1
Food consumed during 24 hour testing period0.00.00.10.20.20.00.3
Glucose0.50.71.3-0.81.10.90.2
Heart weight0.61.01.40.11.11.2-1.1
Delay discount hyperbolic model param0.30.30.5-0.5-0.6-0.50.7
Intermittent access intake day 1-15 change0.50.92.7-0.8-0.4-0.51.6
Intermittent access intake escalation0.61.13.3-0.7-0.6-0.61.8
Intermittent access intake escalation 20.30.41.3-0.40.10.01.1
Intermitt. access day 1 inactive lever presses0.31.02.3-0.6-0.9-0.91.5
Intermitt. access day 15 inactive lever presses1.02.04.00.12.01.9-0.5
Intermitt. access escalation Index0.91.33.2-1.8-0.7-0.70.9
Intermittent access day 1 total infusions2.13.36.50.22.52.5-0.7
Intermittent access day 15 total infusions0.91.31.8-0.91.11.11.3
Intermittent access terminal intake (last 3 days)1.92.55.10.02.32.20.2
Intermittent access total infusions1.92.34.8-0.12.22.10.0
Intermittent access day 1 locomotion1.72.24.20.02.02.0-0.7
Intermittent access day 15 total locomotion1.01.93.8-0.21.91.90.2
Intermittent access total locomotion3.94.79.50.03.13.0-0.3
Intermittent access standard deviation0.20.40.9-0.70.30.11.0
Left kidney weight0.21.12.60.11.41.6-0.3
Right kidney weight0.10.71.4-0.01.11.2-0.2
Cocaine induced anxiety0.00.00.10.3-0.2-0.10.3
Post-drug Anxiety0.10.10.50.70.10.2-0.3
Baseline Anxiety0.10.20.60.10.1-0.0-0.8
Tail length0.20.71.40.2-1.2-1.0-0.4
Length with tail0.52.14.32.1-0.10.1-2.1
Length without tail1.23.78.42.30.70.9-2.9
Times rat made contact with spout0.30.30.60.0-0.8-0.8-0.1
Lifetime Intake1.21.53.1-0.21.81.7-0.1
Liver weight, left0.30.51.2-1.1-0.6-0.50.6
Liver weight, right0.30.51.9-1.4-0.00.10.1
Locomotor activity1.72.64.5-2.1-1.5-1.51.2
Locomotor testing distance1.22.24.7-2.2-1.3-1.30.9
Locomotor testing rearing0.81.63.9-2.0-0.8-1.10.8
Locomotion in novel chamber0.70.91.8-1.30.30.61.1
Locomotion in novel chamber post-restriction0.91.53.1-1.6-0.4-0.41.8
Fecal boli incidents, locomotor time 10.20.30.40.60.60.30.5
Change in fecal boli incidents, locomotor task0.40.40.8-0.7-0.9-0.60.3
Time >=10cm from walls, locomotor time 10.60.71.00.41.01.00.7
Bouts of movement, locomotor time 11.61.83.8-0.2-1.6-2.01.1
Total resting periods, locomotor time 11.01.12.4-0.1-1.2-1.50.9
Rest time, locomotor task time 11.41.72.8-1.5-0.9-1.01.7
Distance moved, locomotor task time 10.50.81.01.00.70.8-1.0
Vertical activity count, locomotor time 10.10.10.4-0.6-0.1-0.10.4
Time >=10cm from walls, locomotor time 20.10.20.60.3-0.2-0.10.8
Bouts of movement, locomotor time 20.10.20.40.50.30.1-0.6
Total resting periods, locomotor time 20.00.10.10.30.30.0-0.4
Rest time, locomotor task time 21.01.32.3-1.5-0.4-0.71.4
Distance moved, locomotor task time 20.30.30.50.7-0.7-0.4-0.2
Vertical activity count, locomotor time 20.71.02.1-0.8-1.4-1.10.5
Light reinforcement 10.20.40.7-0.6-0.4-0.80.7
Light reinforcement 20.60.91.8-0.11.21.40.5
Bone: marrow area1.53.36.80.72.32.6-0.7
Bone: maximum diameter1.32.43.21.61.81.8-0.9
Bone: maximum force5.610.012.53.53.03.0-3.0
Bone: maximum moment5.49.511.73.43.03.0-3.0
Mean time between licks in bursts0.20.40.7-0.2-0.7-0.6-0.8
Mean num. licks in bursts0.40.60.8-0.8-0.7-0.90.5
Cd content in liver0.40.51.6-0.50.20.11.3
Co content in liver3.64.25.2-2.2-2.1-2.31.5
Cu content in liver0.20.30.4-0.4-0.5-0.60.5
Fe content in liver0.10.10.20.50.20.1-0.1
K content in liver1.21.22.01.01.00.8-1.4
Mg content in liver0.40.50.7-0.10.70.7-0.9
Mn content in liver2.22.44.22.11.31.3-1.5
Na content in liver0.50.60.70.70.80.6-0.9
Rb content in liver1.21.42.9-1.7-0.8-0.91.0
Se content in liver0.40.51.10.40.81.00.0
Sr content in liver1.10.91.91.40.60.4-1.1
Zn content in liver0.91.01.90.80.90.9-1.4
Bone: minimum diameter2.03.76.70.82.62.6-0.9
Time in familiar zone, hab. session 10.30.30.5-0.10.70.5-0.6
Time in novel zone, hab. session 10.80.80.8-0.9-0.9-0.90.8
Total zone transitions, hab. session 10.10.10.30.60.30.0-0.3
Total locomotion distance, hab. session 10.50.71.40.11.21.0-0.4
Locomotion velocity, hab. session 10.40.50.90.30.90.9-0.6
Time in familiar zone, hab. session 20.90.91.9-0.7-1.1-1.40.2
Time in novel zone, hab. session 21.11.11.80.81.31.3-0.7
Total zone transitions, hab. session 21.72.13.41.41.41.1-1.8
Total locomotion distance, hab. session 21.01.32.00.51.41.1-1.3
Locomotion velocity, hab. session 21.01.32.20.41.51.3-1.1
Novel to familiar place preference ratio1.51.72.81.71.31.1-1.0
Novelty place preference1.01.11.71.31.21.0-0.7
Time in novel zone, NPP test1.11.21.61.31.21.0-0.8
Total zone transitions, NPP test0.91.22.61.60.90.4-1.2
Total locomotion distance, NPP test0.30.40.60.30.70.5-0.7
Locomotion velocity, NPP test0.10.10.20.10.40.2-0.4
Intraocular pressure0.91.32.9-1.1-0.8-0.81.7
Delta distance traveled before/after self-admin1.81.92.8-1.4-1.7-1.50.9
Delta ambulatory time before/after self-admin1.51.92.8-1.2-1.7-1.50.9
Parametrial fat weight0.61.52.7-0.91.61.50.6
Pavlov. Cond. lever latency1.11.41.8-1.1-1.0-1.21.4
Pavlov. Cond. magazine entry latency1.82.42.91.51.51.7-1.4
Pavlov. Cond. change in total contacts1.01.01.51.01.21.2-0.5
Pavlov. Cond. index score0.60.70.90.90.70.9-0.8
Pavlov. Cond. latency score0.60.81.00.90.71.0-1.0
Pavlov. Cond. lever contacts0.91.11.50.90.91.2-1.2
Pavlov. Cond. magazine entry number0.30.40.7-0.8-0.4-0.60.4
Pavlov. Cond. intertrial magazine entries0.50.61.1-0.9-0.7-0.41.0
Pavlov. Cond. lever-magazine prob. diff.0.81.01.41.20.91.0-0.9
Pavlov. Cond. response bias0.60.80.90.90.81.0-0.8
Pavlov. Cond. lever latency1.11.52.0-1.4-1.4-1.00.9
Pavlov. Cond. magazine entry latency2.73.36.22.51.41.1-2.0
Pavlov. Cond. change in total contacts1.62.64.12.01.81.3-1.1
Pavlov. Cond. index score2.63.55.82.41.81.5-1.6
Pavlov. Cond. latency score2.73.65.62.41.91.6-1.6
Pavlov. Cond. lever contacts0.81.11.71.31.20.7-0.7
Pavlov. Cond. magazine entry number4.15.19.8-3.1-1.7-1.42.3
Pavlov. Cond. intertrial magazine entries0.91.42.4-1.2-0.9-0.91.6
Pavlov. Cond. lever-magazine prob. diff.2.83.76.42.51.81.4-1.7
Pavlov. Cond. response bias2.02.74.42.11.61.1-1.6
Time in closed arm before self-admin0.20.30.50.30.60.70.2
Time in closed arm after self-admin0.20.20.6-0.8-0.2-0.10.4
Time in open arm before self-admin1.41.52.3-1.5-1.4-1.30.2
Time in open arm after self-admin0.10.10.20.1-0.4-0.4-0.4
Bone: periosteal estimation2.04.36.91.62.42.6-1.4
Bone: periosteal perimeter1.32.33.21.41.81.8-1.0
Bone: post-yield work5.78.413.42.53.73.3-1.7
Progressive ratio test 1 active lever presses0.10.20.40.00.60.50.2
Progressive ratio test 1 breakpoint0.20.20.40.10.60.50.3
Progressive ratio test 1 inactive lever presses0.50.71.00.51.01.0-0.7
Progressive ratio test 2 active lever presses0.80.91.7-0.31.31.20.7
Progressive ratio test 2 breakpoint0.80.91.8-0.31.41.20.6
Progressive ratio test 2 inactive lever presses0.91.41.7-1.31.31.20.9
Incentive Sensitization - Responses1.62.23.91.12.01.7-0.9
Incentive Sensitization - Breakpoint1.62.03.80.92.01.7-0.8
Sum of active levers in priming session hrs 5-60.30.30.60.8-0.3-0.20.6
Progressive ratio1.42.23.51.91.41.5-1.1
Food seeking constrained by brief footshock0.00.00.00.00.0-0.2-0.2
Seeking ratio, punishment vs. effort1.11.53.6-1.9-1.0-1.10.5
Retroperitoneal fat weight1.45.918.1-4.30.30.42.3
Reaction time number correct0.10.10.20.2-0.5-0.3-0.4
Reaction time mean minus median0.41.14.12.00.40.1-0.2
Reaction time mean minus median AUC0.20.51.31.1-0.3-0.80.2
Reaction time num false alarms1.31.82.5-1.6-1.6-1.50.5
Reaction time num false alarms AUC0.81.32.2-1.0-1.4-1.50.2
Reaction time trials correct on left0.10.10.20.2-0.5-0.3-0.4
Reaction time trials on left0.10.20.30.3-0.5-0.4-0.4
Reaction time mean0.10.10.30.40.0-0.20.6
Reaction time mean AUC0.10.10.30.40.1-0.20.6
Median of all reaction times0.00.10.10.00.20.10.4
Reaction time omissions0.20.30.60.8-0.4-0.70.0
Reaction time false alarm rate0.60.81.2-1.1-0.9-1.00.4
Reaction time premature initiation rate2.02.87.41.40.71.1-2.7
Reaction time premature initiations1.32.36.61.50.10.5-2.6
Std. dev. reaction times0.30.50.70.6-0.5-0.90.9
Reaction time trials completed0.10.20.30.3-0.5-0.4-0.4
Reaction time trials AUC0.10.10.20.2-0.5-0.3-0.5
Run reversals in cocaine runway, females1.81.83.20.81.31.2-1.8
Run reversals in cocaine runway, males0.50.50.6-0.6-0.8-0.80.7
Latency to leave start box in cocaine runway0.40.40.90.10.91.0-0.0
Latency to leave start box in cocaine runway, F0.20.20.40.10.60.6-0.3
Latency to leave start box in cocaine runway, M0.10.10.20.10.20.40.4
Std. dev. time between licks in bursts0.40.50.8-0.8-0.7-0.9-0.2
Total sessions with >9 infusions0.10.10.20.20.40.40.1
Short access day 1 total inactive lever presses0.40.51.1-1.00.60.80.0
Short access day 10 total inactive lever presses0.41.01.4-1.2-0.9-0.81.0
Short access day 1 total infusions0.20.30.7-0.80.40.50.4
Short access day 10 total infusions1.31.73.10.61.81.7-0.7
Short access total infusions0.60.71.40.11.21.1-0.4
Short access day 1 locomotion0.40.51.0-0.30.81.00.1
Short access day 10 total locomotion1.01.83.40.41.81.8-0.7
Short access total locomotion1.11.42.80.11.71.6-0.4
One hour access (0.3 mA shock)1.11.32.1-1.01.51.30.6
Compulsive drug intake0.60.71.30.0-1.0-1.1-0.5
One hour access (shock baseline)1.01.51.9-1.11.21.21.4
Social response proportion0.10.10.40.70.1-0.0-0.3
Social responses0.30.40.90.1-0.9-0.90.0
Social time0.40.61.40.0-0.9-1.20.2
Soleus weight1.64.810.8-3.3-1.3-1.02.4
Bone: stiffness2.03.26.12.51.61.6-1.4
Tibialis anterior weight4.114.030.15.51.91.9-4.3
Time to tail flick, vehicle, before self-admin1.61.72.21.31.51.4-1.0
Time to tail flick, vehicle, after self-admin1.62.12.81.71.21.2-1.6
Time to tail flick, test, before self-admin1.92.03.20.01.71.8-1.4
Time to tail flick, test, after self-admin0.10.10.40.60.10.20.1
Delta time to tail flick, vehicle, before/after SA0.20.20.50.60.10.1-0.7
Delta time to tail flick, test, before/after SA1.51.62.40.2-1.5-1.51.3
Tibia length1.13.15.42.31.71.7-1.1
Ambulatory time before self-admin1.32.02.71.61.41.5-1.1
Ambulatory time after self-admin0.10.10.20.00.30.5-0.0
Bone: tissue strength0.10.10.2-0.3-0.3-0.40.4
Total heroin consumption1.31.54.02.0-0.10.2-1.4
Bone: trabecular number2.43.94.52.12.12.0-1.6
Bone: trabecular spacing2.02.63.3-1.7-1.8-1.81.1
Bone: trabecular thickness0.10.10.4-0.10.40.60.1
Bone: trabecular tissue density5.27.710.9-2.9-2.3-2.43.3
Water consumed over 24 hours0.20.30.40.5-0.7-0.4-0.4
Number of responses in last shaping day0.10.10.1-0.3-0.3-0.1-0.4