Castor1

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 Adipose gene expression Castor1 0.09 0.04 1.4e-09 0.08 0.07 0.06 0.07 3.6e-09 6.0e-08 5.9e-07 3.4e-08
2 NAcc gene expression Castor1 0.05 0.03 2.2e-06 0.06 0.05 0.06 0.06 2.4e-07 1.1e-06 1.0e-07 9.4e-08
3 BLA gene expression Castor1 0.15 0.1 9.2e-03 0.03 0.01 0.03 0.03 1.3e-02 1.2e-01 1.2e-02 1.4e-02

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
Bone: apparent density0.61.11.81.0-0.6-1.3
Area under the delay curve7.17.59.7-3.13.1-1.8
BMI with tail0.41.41.7-0.91.31.3
BMI without tail0.10.20.4-0.20.40.6
Body weight1.25.17.7-2.42.8-1.3
Weight adjusted by age0.60.71.2-0.91.1-0.3
Bone surface0.30.61.2-0.40.6-1.1
Bone volume0.30.51.3-0.20.3-1.1
Breakpoint value in progressive ratio session0.50.50.9-0.70.9-0.5
Number of licking bursts0.91.12.1-1.40.9-0.5
Context. condit. distance diff. score2.42.84.12.0-1.81.0
Locomotion velocity, session 10.10.10.30.10.00.6
Locomotion distance, session 11.31.73.71.9-1.1-0.3
Locomotion velocity, session 20.10.20.4-0.0-0.3-0.6
Locomotion distance, session 21.31.82.71.6-1.50.8
Locomotion velocity, session 30.91.12.31.5-1.0-0.2
Locomotion distance, session 30.40.40.8-0.50.9-0.3
Stereotopy head waving bouts, day 31.21.42.51.6-1.20.5
Stereotopy head waving duration, day 32.53.04.42.1-1.71.3
Locomotion velocity, session 70.10.10.30.00.50.2
Locomotion distance, session 70.60.71.31.2-0.50.7
Stereotopy head waving bouts, day 70.00.00.0-0.00.00.1
Stereotopy head waving duration, day 70.10.10.2-0.30.40.0
Locomotion velocity, session 81.51.82.51.2-1.61.3
Locomotion distance, session 86.910.114.13.8-3.52.1
Degree of sensitization distance4.74.75.72.4-2.41.6
Degree of sensitization stereotypy0.30.40.7-0.80.6-0.3
Conditioned locomotion0.20.30.50.7-0.5-0.1
Cocaine response after cond. corrected1.41.52.6-0.81.1-1.6
Cocaine response after cond. not corrected3.23.56.7-1.11.6-2.6
Cocaine response before conditioning0.30.41.10.20.2-1.0
Saline control response0.70.91.61.3-0.80.6
Bone: connectivity density0.60.91.31.1-0.5-1.1
Bone: cortical apparent density0.20.30.9-0.30.10.9
Bone: cortical area2.13.14.7-1.92.2-1.1
Bone: cortical porosity1.01.02.6-0.40.6-1.6
Bone: cortical porosity0.20.30.7-0.10.1-0.9
Bone: cortical thickness1.21.93.0-1.41.7-0.9
Bone: cortical thickness0.20.30.70.9-0.20.2
Bone: cortical tissue density0.30.30.9-0.30.10.9
Condit. Reinf. active minus inactive responses0.20.20.40.10.50.6
Condit. Reinf. active-inactive response ratio1.21.33.70.2-0.21.9
Condit. Reinf. active responses0.20.30.50.50.20.7
Condit. Reinf. inactive responses0.10.20.40.20.2-0.6
Incentive salience index mean0.40.51.3-0.51.10.1
Condit. Reinf. lever presses0.40.51.2-0.61.1-0.0
Condit. Reinf. active minus inactive responses0.60.71.31.2-0.60.7
Condit. Reinf. active-inactive response ratio0.10.10.20.10.40.0
Condit. Reinf. active responses1.52.13.71.9-1.21.1
Condit. Reinf. inactive responses2.53.13.82.0-1.81.5
Condit. Reinf. lever presses0.50.71.21.1-0.60.7
Condit. Reinf. lever reinforcers received1.31.83.61.9-1.00.9
Conditioned reinforcement - actives0.90.91.5-1.21.1-0.4
Conditioned reinforcement - inactives2.32.63.4-1.71.8-1.2
Conditioned locomotion0.50.60.7-0.70.7-0.9
Patch foraging indifference point 0 sec0.10.20.3-0.60.40.2
Patch foraging indifference point AUC0.10.10.30.50.10.2
Delay disc. indifference point, 0s delay4.75.26.5-2.42.6-1.9
Delay disc. indifference point, 16s delay4.64.76.7-2.62.4-1.4
Delay disc. indifference point, 2s delay4.65.16.3-2.22.5-2.1
Delay disc. indifference point, 24s delay3.33.14.3-2.12.0-1.0
Delay disc. indifference point, 4s delay5.15.46.3-2.02.5-2.4
Delay disc. indifference point, 8s delay1.61.62.6-1.61.5-0.0
Indifference point function ln k0.50.71.80.5-0.2-1.3
Indifference point function log k0.50.71.80.5-0.2-1.3
Patch foraging total patch changes 0 sec0.20.30.5-0.70.50.3
Patch foraging total patch changes 12 sec0.10.10.3-0.30.50.1
Patch foraging total patch changes 18 sec0.70.91.5-1.20.90.6
Patch foraging total patch changes 24 sec2.83.55.71.6-1.52.4
Patch foraging total patch changes 6 sec0.10.20.3-0.50.50.1
Patch foraging time to switch 0 sec1.42.25.5-0.3-1.02.3
Patch foraging water rate 0 sec0.60.91.8-0.51.3-0.9
Patch foraging water rate 12 sec0.20.30.50.60.2-0.7
Patch foraging water rate 18 sec1.11.53.21.8-1.1-0.5
Patch foraging water rate 24 sec0.10.20.40.4-0.30.6
Patch foraging water rate 6 sec0.10.20.40.40.1-0.6
Seeking ratio, delayed vs. immediate footshock0.40.41.30.0-0.0-1.1
Distance traveled before self-admin0.30.30.7-0.20.4-0.8
Distance traveled after self-admin0.30.40.6-0.60.8-0.4
Extensor digitorum longus weight2.28.813.2-3.03.6-2.1
Bone: elastic displacement0.60.61.00.9-1.00.3
Bone: elastic work0.00.10.10.2-0.3-0.2
Bone: endosteal estimation0.10.10.3-0.30.1-0.5
Bone: endosteal perimeter0.81.61.8-1.21.2-1.3
Epididymis fat weight0.10.20.3-0.5-0.20.4
Delta time in closed arm before/after self-admin3.33.25.72.4-1.8-0.8
Delta time in open arm before/after self-admin4.44.28.1-2.82.2-0.2
Diff in mean of infusions in LGA sessions0.10.10.2-0.50.3-0.3
Delay discount exponential model param6.56.99.03.0-3.01.7
Extinction: sum of active levers before priming1.71.72.4-1.21.6-1.1
Active lever presses in extinction session 62.42.34.7-0.91.2-2.2
Bone: final force0.20.30.7-0.50.30.8
Bone: final moment0.10.20.4-0.60.40.2
Food consumed during 24 hour testing period0.30.30.4-0.60.40.7
Glucose11.216.928.2-5.34.7-0.6
Heart weight0.40.60.70.6-0.90.8
Delay discount hyperbolic model param6.87.29.53.0-3.11.8
Intermittent access intake day 1-15 change0.70.71.1-1.01.00.0
Intermittent access intake escalation0.40.40.6-0.80.70.1
Intermittent access intake escalation 20.30.30.7-0.3-0.00.9
Intermitt. access day 1 inactive lever presses0.50.51.11.0-0.4-0.6
Intermitt. access day 15 inactive lever presses0.70.70.9-0.71.00.9
Intermitt. access escalation Index0.50.51.20.4-0.21.1
Intermittent access day 1 total infusions0.40.51.1-1.10.40.3
Intermittent access day 15 total infusions1.31.42.9-1.71.10.3
Intermittent access terminal intake (last 3 days)0.50.61.2-1.10.60.4
Intermittent access total infusions0.60.71.3-1.20.7-0.3
Intermittent access day 1 locomotion0.91.01.40.7-1.21.0
Intermittent access day 15 total locomotion1.61.73.1-1.81.40.3
Intermittent access total locomotion0.60.71.3-1.10.9-0.3
Intermittent access standard deviation1.01.01.70.8-0.91.3
Left kidney weight0.51.32.11.3-1.40.4
Right kidney weight0.71.93.01.6-1.7-0.1
Cocaine induced anxiety0.91.11.7-1.11.3-0.6
Post-drug Anxiety0.70.80.9-1.01.0-0.7
Baseline Anxiety0.30.30.8-0.90.30.3
Tail length1.23.68.5-0.81.2-2.9
Length with tail1.24.47.0-1.81.8-2.7
Length without tail0.92.13.3-1.81.5-0.7
Times rat made contact with spout0.40.51.2-0.6-0.01.1
Lifetime Intake0.50.61.0-1.00.7-0.5
Liver weight, left3.14.98.22.9-2.6-0.1
Liver weight, right4.47.210.03.2-2.91.9
Locomotor activity0.50.81.2-1.10.8-0.7
Locomotor testing distance0.60.92.0-0.40.8-1.4
Locomotor testing rearing1.52.43.2-1.51.8-1.3
Locomotion in novel chamber0.50.60.90.8-1.00.3
Locomotion in novel chamber post-restriction1.72.02.51.3-1.41.6
Fecal boli incidents, locomotor time 10.20.20.50.7-0.4-0.2
Change in fecal boli incidents, locomotor task0.80.81.70.2-0.81.3
Time >=10cm from walls, locomotor time 11.31.33.3-0.0-0.71.8
Bouts of movement, locomotor time 10.20.30.6-0.30.30.8
Total resting periods, locomotor time 10.40.51.1-0.50.41.0
Rest time, locomotor task time 10.10.20.40.2-0.10.7
Distance moved, locomotor task time 10.00.00.1-0.30.0-0.3
Vertical activity count, locomotor time 10.10.10.1-0.3-0.3-0.3
Time >=10cm from walls, locomotor time 20.30.30.40.6-0.70.5
Bouts of movement, locomotor time 20.40.51.10.50.2-1.0
Total resting periods, locomotor time 20.30.30.70.50.4-0.8
Rest time, locomotor task time 20.70.81.50.7-0.61.2
Distance moved, locomotor task time 20.10.10.10.0-0.20.4
Vertical activity count, locomotor time 20.20.20.5-0.10.2-0.7
Light reinforcement 10.20.30.50.20.7-0.5
Light reinforcement 20.10.10.10.3-0.10.3
Bone: marrow area0.00.10.2-0.30.1-0.5
Bone: maximum diameter3.05.47.9-2.82.8-0.8
Bone: maximum force0.10.10.2-0.40.20.4
Bone: maximum moment0.10.10.3-0.50.2-0.2
Mean time between licks in bursts0.40.51.40.2-0.31.2
Mean num. licks in bursts1.01.21.71.0-1.01.3
Cd content in liver1.41.62.31.1-1.01.5
Co content in liver0.50.61.2-0.20.8-1.1
Cu content in liver0.10.10.2-0.00.4-0.5
Fe content in liver0.30.40.50.5-0.60.7
K content in liver1.41.42.7-0.41.1-1.6
Mg content in liver2.12.23.0-1.11.7-1.5
Mn content in liver1.11.23.60.10.4-1.9
Na content in liver5.55.87.9-2.72.8-1.6
Rb content in liver1.01.22.0-1.11.4-0.7
Se content in liver1.82.22.7-1.41.4-1.7
Sr content in liver1.01.01.6-1.30.60.9
Zn content in liver0.10.10.2-0.40.20.1
Bone: minimum diameter0.10.30.70.2-0.1-0.8
Time in familiar zone, hab. session 12.22.04.00.6-1.32.0
Time in novel zone, hab. session 11.81.73.9-0.61.0-2.0
Total zone transitions, hab. session 10.40.50.60.7-0.60.8
Total locomotion distance, hab. session 10.10.20.50.20.1-0.7
Locomotion velocity, hab. session 10.10.10.20.4-0.0-0.2
Time in familiar zone, hab. session 23.83.95.21.9-2.31.8
Time in novel zone, hab. session 22.52.53.8-1.81.9-0.8
Total zone transitions, hab. session 20.10.20.5-0.20.70.1
Total locomotion distance, hab. session 20.30.50.8-0.70.90.2
Locomotion velocity, hab. session 20.60.81.4-0.91.2-0.2
Novel to familiar place preference ratio1.71.72.9-1.41.7-0.6
Novelty place preference1.61.62.7-1.31.7-0.7
Time in novel zone, NPP test1.31.32.3-1.21.5-0.6
Total zone transitions, NPP test0.40.50.90.9-0.2-0.8
Total locomotion distance, NPP test0.20.30.70.30.2-0.8
Locomotion velocity, NPP test0.40.51.30.00.6-1.1
Intraocular pressure1.01.43.8-0.60.02.0
Delta distance traveled before/after self-admin0.10.10.3-0.40.20.5
Delta ambulatory time before/after self-admin0.10.10.10.1-0.30.2
Parametrial fat weight1.01.93.2-1.81.50.3
Pavlov. Cond. lever latency0.40.50.6-0.80.6-0.7
Pavlov. Cond. magazine entry latency0.81.11.41.2-1.10.7
Pavlov. Cond. change in total contacts0.60.71.1-0.91.00.1
Pavlov. Cond. index score1.01.31.91.4-1.20.8
Pavlov. Cond. latency score1.01.31.71.3-1.20.9
Pavlov. Cond. lever contacts0.60.70.80.8-0.70.9
Pavlov. Cond. magazine entry number0.81.01.5-1.20.9-0.7
Pavlov. Cond. intertrial magazine entries1.01.23.10.30.6-1.8
Pavlov. Cond. lever-magazine prob. diff.1.11.42.11.5-1.20.8
Pavlov. Cond. response bias0.10.10.30.6-0.20.0
Pavlov. Cond. lever latency1.41.92.2-1.51.4-1.2
Pavlov. Cond. magazine entry latency0.00.00.0-0.20.1-0.0
Pavlov. Cond. change in total contacts0.81.11.60.9-0.91.3
Pavlov. Cond. index score0.81.01.11.0-1.00.9
Pavlov. Cond. latency score0.70.91.21.1-1.00.8
Pavlov. Cond. lever contacts1.72.42.61.6-1.61.4
Pavlov. Cond. magazine entry number0.10.10.30.30.1-0.6
Pavlov. Cond. intertrial magazine entries0.20.20.40.60.1-0.5
Pavlov. Cond. lever-magazine prob. diff.0.81.01.11.0-1.10.9
Pavlov. Cond. response bias0.81.11.50.8-1.01.2
Time in closed arm before self-admin0.60.61.30.3-0.71.1
Time in closed arm after self-admin2.62.74.82.2-1.80.2
Time in open arm before self-admin0.20.20.30.60.0-0.5
Time in open arm after self-admin0.91.01.5-1.21.10.2
Bone: periosteal estimation0.61.31.6-1.31.3-0.8
Bone: periosteal perimeter2.13.54.4-2.02.1-1.5
Bone: post-yield work0.30.50.8-0.90.2-0.8
Progressive ratio test 1 active lever presses2.02.12.2-1.31.5-1.5
Progressive ratio test 1 breakpoint1.51.52.1-1.01.2-1.4
Progressive ratio test 1 inactive lever presses3.63.95.1-1.72.3-1.9
Progressive ratio test 2 active lever presses1.01.02.2-1.50.9-0.0
Progressive ratio test 2 breakpoint0.91.02.1-1.51.00.1
Progressive ratio test 2 inactive lever presses2.73.15.4-2.31.71.0
Incentive Sensitization - Responses1.11.12.0-1.40.51.0
Incentive Sensitization - Breakpoint1.11.12.1-1.40.51.1
Sum of active levers in priming session hrs 5-60.90.91.4-0.91.2-0.7
Progressive ratio0.40.41.3-1.1-0.10.0
Food seeking constrained by brief footshock2.83.44.2-2.11.7-1.8
Seeking ratio, punishment vs. effort1.01.21.9-0.01.3-1.4
Retroperitoneal fat weight0.31.13.1-0.3-0.11.8
Reaction time number correct0.20.20.40.3-0.60.5
Reaction time mean minus median0.10.10.2-0.4-0.1-0.3
Reaction time mean minus median AUC0.10.20.2-0.50.4-0.4
Reaction time num false alarms0.50.71.0-0.81.0-0.6
Reaction time num false alarms AUC0.30.40.5-0.60.7-0.5
Reaction time trials correct on left0.20.20.40.3-0.60.5
Reaction time trials on left0.20.20.40.3-0.60.5
Reaction time mean0.00.00.10.1-0.3-0.0
Reaction time mean AUC0.10.10.20.3-0.4-0.1
Median of all reaction times0.10.10.20.3-0.40.0
Reaction time omissions0.20.30.40.6-0.5-0.4
Reaction time false alarm rate0.50.80.9-0.81.0-0.9
Reaction time premature initiation rate3.44.68.22.9-2.11.1
Reaction time premature initiations2.63.75.92.4-1.91.2
Std. dev. reaction times0.10.10.2-0.5-0.10.4
Reaction time trials completed0.20.20.40.3-0.60.5
Reaction time trials AUC0.10.10.20.2-0.50.3
Run reversals in cocaine runway, females0.70.71.5-0.2-0.71.2
Run reversals in cocaine runway, males0.00.00.1-0.20.2-0.0
Latency to leave start box in cocaine runway0.30.30.6-0.10.5-0.8
Latency to leave start box in cocaine runway, F0.00.00.1-0.20.30.1
Latency to leave start box in cocaine runway, M1.31.23.40.10.6-1.8
Std. dev. time between licks in bursts0.70.81.11.1-0.70.8
Total sessions with >9 infusions0.10.10.20.5-0.1-0.3
Short access day 1 total inactive lever presses1.61.72.5-1.11.6-1.2
Short access day 10 total inactive lever presses0.91.01.81.4-0.4-1.0
Short access day 1 total infusions0.30.30.90.30.1-1.0
Short access day 10 total infusions0.10.10.10.2-0.4-0.2
Short access total infusions0.10.10.30.2-0.1-0.6
Short access day 1 locomotion1.31.42.41.6-1.3-0.1
Short access day 10 total locomotion1.41.62.41.5-1.50.1
Short access total locomotion3.94.06.62.6-2.30.6
One hour access (0.3 mA shock)0.50.51.1-1.00.7-0.3
Compulsive drug intake3.93.96.42.5-2.30.2
One hour access (shock baseline)3.84.46.7-2.62.3-1.1
Social response proportion0.60.61.2-0.70.21.1
Social responses0.50.61.0-1.00.50.8
Social time0.40.51.10.1-1.00.7
Soleus weight0.72.03.1-1.31.1-1.8
Bone: stiffness0.10.10.20.40.00.4
Tibialis anterior weight1.34.67.4-2.72.3-1.0
Time to tail flick, vehicle, before self-admin0.20.20.3-0.5-0.20.5
Time to tail flick, vehicle, after self-admin0.10.10.30.1-0.50.2
Time to tail flick, test, before self-admin0.80.81.31.2-1.00.2
Time to tail flick, test, after self-admin1.01.01.6-1.01.3-0.7
Delta time to tail flick, vehicle, before/after SA0.10.10.10.3-0.40.1
Delta time to tail flick, test, before/after SA1.31.11.8-1.11.4-0.5
Tibia length0.10.30.60.6-0.8-0.2
Ambulatory time before self-admin0.40.41.2-0.20.3-1.1
Ambulatory time after self-admin0.30.30.5-0.50.5-0.7
Bone: tissue strength0.50.50.90.1-0.90.8
Total heroin consumption0.10.10.40.2-0.0-0.6
Bone: trabecular number0.60.92.20.7-0.2-1.5
Bone: trabecular spacing1.21.52.2-1.40.61.5
Bone: trabecular thickness1.72.12.61.5-1.1-1.6
Bone: trabecular tissue density0.30.40.8-0.1-0.50.9
Water consumed over 24 hours0.40.51.3-0.3-0.51.1
Number of responses in last shaping day0.30.30.60.5-0.50.7