Genome-wide scans of genetic variants for psychophysiological endophenotypes: a methodological overview

Bistable Perception Discriminates Between Depressive Patients, Controls, Schizophrenia Patients, and Their Siblings

BACKGROUND AND HYPOTHESIS

Individuals with schizophrenia have less priors than controls, meaning they rely less upon their prior experiences to interpret the current stimuli. These differences in priors are expected to show as higher alternation rates in bistable perception tasks like the Structure-from-Motion (SfM) paradigm. In this paradigm, continuously moving dots in two dimensions are perceived subjectively as traveling along a three-dimensional sphere, which results in a direction of motion (left or right) that shifts approximately every few seconds.

STUDY DESIGN

Here, we tested healthy controls, patients with schizophrenia, siblings of patients with schizophrenia, and patients with depression with both the intermittent and continuous variants of the SfM paradigm.

STUDY RESULTS

In the intermittent variant of the SfM paradigm, depressive patients exhibited the lowest alternation rate, followed by unaffected controls. In contrast, patients with schizophrenia and their unaffected siblings displayed significantly higher alternation rates. In the continuous variant of the SfM paradigm, patients with schizophrenia showed the lowest mean percept durations, while there were no differences between the other three groups.

CONCLUSIONS

The intermittent SfM paradigm is a candidate endophenotype for schizophrenia. The aberrant processing in the patients may stem from alterations in adaptation and/or cross-inhibition mechanisms leading to changes in priors, as suggested by current models in the field. The intermittent SfM paradigm is, hence, a trait marker that offers the great opportunity to investigate perceptual abnormalities across the psychiatry spectrum, ranging from depression to psychosis.

A genome-wide association study of interhemispheric theta EEG coherence: implications for neural connectivity and alcohol use behavior

Aberrant connectivity of large-scale brain networks has been observed among individuals with alcohol use disorders (AUDs) as well as in those at risk, suggesting deficits in neural communication between brain regions in the liability to develop AUD. Electroencephalographical (EEG) coherence, which measures the degree of synchrony between brain regions, may be a useful measure of connectivity patterns in neural networks for studying the genetics of AUD. In 8810 individuals (6644 of European and 2166 of African ancestry) from the Collaborative Study on the Genetics of Alcoholism (COGA), we performed a Multi-Trait Analyses of genome-wide association studies (MTAG) on parietal resting-state theta (3-7 Hz) EEG coherence, which previously have been associated with AUD. We also examined developmental effects of GWAS findings on trajectories of neural connectivity in a longitudinal subsample of 2316 adolescent/young adult offspring from COGA families (ages 12-30) and examined the functional and clinical significance of GWAS variants. Six correlated single nucleotide polymorphisms located in a brain-expressed lincRNA (ENSG00000266213) on chromosome 18q23 were associated with posterior interhemispheric low theta EEG coherence (3-5 Hz). These same variants were also associated with alcohol use behavior and posterior corpus callosum volume, both in a subset of COGA and in the UK Biobank. Analyses in the subsample of COGA offspring indicated that the association of rs12954372 with low theta EEG coherence occurred only in females, most prominently between ages 25 and 30 (p < 2 × 10-9). Converging data provide support for the role of genetic variants on chromosome 18q23 in regulating neural connectivity and alcohol use behavior, potentially via dysregulated myelination. While findings were less robust, genome-wide associations were also observed with rs151174000 and parieto-frontal low theta coherence, rs14429078 and parieto-occipital interhemispheric high theta coherence, and rs116445911 with centro-parietal low theta coherence. These novel genetic findings highlight the utility of the endophenotype approach in enhancing our understanding of mechanisms underlying addiction susceptibility.

Using multivariate endophenotypes to identify psychophysiological mechanisms associated with polygenic scores for substance use, schizophrenia, and education attainment

BACKGROUND

To better characterize brain-based mechanisms of polygenic liability for psychopathology and psychological traits, we extended our previous report (Liu et al. Psychophysiological endophenotypes to characterize mechanisms of known schizophrenia genetic loci. Psychological Medicine, 2017), focused solely on schizophrenia, to test the association between multivariate psychophysiological candidate endophenotypes (including novel measures of θ/δ oscillatory activity) and a range of polygenic scores (PGSs), namely alcohol/cannabis/nicotine use, an updated schizophrenia PGS (containing 52 more genome-wide significant loci than the PGS used in our previous report) and educational attainment.

METHOD

A large community-based twin/family sample (N = 4893) was genome-wide genotyped and imputed. PGSs were constructed for alcohol use, regular smoking initiation, lifetime cannabis use, schizophrenia, and educational attainment. Eleven endophenotypes were assessed: visual oddball task event-related electroencephalogram (EEG) measures (target-related parietal P3 amplitude, frontal θ, and parietal δ energy/inter-trial phase clustering), band-limited resting-state EEG power, antisaccade error rate. Principal component analysis exploited covariation among endophenotypes to extract a smaller number of meaningful dimensions/components for statistical analysis.

RESULTS

Endophenotypes were heritable. PGSs showed expected intercorrelations (e.g. schizophrenia PGS correlated positively with alcohol/nicotine/cannabis PGSs). Schizophrenia PGS was negatively associated with an event-related P3/δ component [β = -0.032, nonparametric bootstrap 95% confidence interval (CI) -0.059 to -0.003]. A prefrontal control component (event-related θ/antisaccade errors) was negatively associated with alcohol (β = -0.034, 95% CI -0.063 to -0.006) and regular smoking PGSs (β = -0.032, 95% CI -0.061 to -0.005) and positively associated with educational attainment PGS (β = 0.031, 95% CI 0.003-0.058).

CONCLUSIONS

Evidence suggests that multivariate endophenotypes of decision-making (P3/δ) and cognitive/attentional control (θ/antisaccade error) relate to alcohol/nicotine, schizophrenia, and educational attainment PGSs and represent promising targets for future research.

Male Ejaculatory Endophenotypes: Revealing Internal Inconsistencies of the Concept in Heterosexual Copulating Rats

Distinct manifestations of sexual behavior are conceived as separate phenotypes. Each sexual phenotype is assumed to be associated with a characteristic brain. These notions have justified the phenotyping of heterosexual copulator males based upon their ejaculation's latencies (EL) or frequencies (i.e., cumulative ejaculation number; EN). For instance, men and male rats showing premature, normal or retarded ejaculation are assumed to be distinctive endophenotypes. This concept, nonetheless, contradicts past and recent evidence that supports that sexual behavior is highly variable within each sex, and that the brain sexual functional morphology represents an intricate sexual phenotypic mosaic. Hence, for ejaculatory male endophenotypes to be considered as a valid biological concept, it must show internal consistency at various levels of organization (including genetic architectures), after being challenged by intrinsic and/or extrinsic factors. We then judged the internal consistency of the presumed ejaculatory endophenotypes by assessing whether copulatory behavior and the expression of copulation relevant genes and brain limbic structures are specific to each of the presumed EL- or EN-ejaculatory endophenotypes. To do this, copulating male rats were first phenotyped in groups consistently displaying short, average or long ejaculation latencies or very high, high, average, low or very low EN, based in their copulatory performance. Then, the internal consistency of the presumed EL- or EN-endophenotypes was tested by introducing as covariates of phenotyping other copulatory parameters (e.g., number of intromissions) in addition to EL or EN, or by analyzing the expression levels of genes encoding for estrogen receptor alpha, progesterone receptor, androgen receptor, aromatase, DNA methyl-transferase 3a and DNA methyl-transferase 1 in the amygdala, medial preoptic area, ventromedial hypothalamus and olfactory bulb. We found that even though there were group-level differences in all the variables that were studied, these differences did not add-up to create the presumed EL- or EN-ejaculatory endophenotypes. In fact, the extensive overlapping of copulatory parameters and expression levels of copulation relevant genes in limbic structures across EL- or EN-phenotyped copulating male rats, is not consistent with the hypothesis that distinct ejaculatory endophenotypes exist and that they are associated with specific brain characteristics.

Negative self-referential processing is associated with genetic variation in the serotonin transporter-linked polymorphic region (5-HTTLPR): Evidence from two independent studies

The current research examined whether carriers of the short 5-HTTLPR allele (in SLC6A4), who have been shown to selectively attend to negative information, exhibit a bias towards negative self-referent processing. The self-referent encoding task (SRET) was used to measure self-referential processing of positive and negative adjectives. Ratcliff's diffusion model isolated and extracted decision-making components from SRET responses and reaction times. Across the initial (N = 183) and replication (N = 137) studies, results indicated that short 5-HTTLPR allele carriers more easily categorized negative adjectives as self-referential (i.e., higher drift rate). Further, drift rate was associated with recall of negative self-referential stimuli. Findings across both studies provide further evidence that genetic variation may contribute to the etiology of negatively biased processing of self-referent information. Large scale studies examining the genetic contributions to negative self-referent processing may be warranted.

The utility of twins in developmental cognitive neuroscience research: How twins strengthen the ABCD research design

The ABCD twin study will elucidate the genetic and environmental contributions to a wide range of mental and physical health outcomes in children, including substance use, brain and behavioral development, and their interrelationship. Comparisons within and between monozygotic and dizygotic twin pairs, further powered by multiple assessments, provide information about genetic and environmental contributions to developmental associations, and enable stronger tests of causal hypotheses, than do comparisons involving unrelated children. Thus a sub-study of 800 pairs of same-sex twins was embedded within the overall Adolescent Brain and Cognitive Development (ABCD) design. The ABCD Twin Hub comprises four leading centers for twin research in Minnesota, Colorado, Virginia, and Missouri. Each site is enrolling 200 twin pairs, as well as singletons. The twins are recruited from registries of all twin births in each State during 2006-2008. Singletons at each site are recruited following the same school-based procedures as the rest of the ABCD study. This paper describes the background and rationale for the ABCD twin study, the ascertainment of twin pairs and implementation strategy at each site, and the details of the proposed analytic strategies to quantify genetic and environmental influences and test hypotheses critical to the aims of the ABCD study.

Brain imaging genetics in ADHD and beyond - Mapping pathways from gene to disorder at different levels of complexity

Attention-deficit/hyperactivity disorder (ADHD) is a common and often persistent neurodevelopmental disorder. Beyond gene-finding, neurobiological parameters, such as brain structure, connectivity, and function, have been used to link genetic variation to ADHD symptomatology. We performed a systematic review of brain imaging genetics studies involving 62 ADHD candidate genes in childhood and adult ADHD cohorts. Fifty-one eligible research articles described studies of 13 ADHD candidate genes. Almost exclusively, single genetic variants were studied, mostly focussing on dopamine-related genes. While promising results have been reported, imaging genetics studies are thus far hampered by methodological differences in study design and analysis methodology, as well as limited sample sizes. Beyond reviewing imaging genetics studies, we also discuss the need for complementary approaches at multiple levels of biological complexity and emphasize the importance of combining and integrating findings across levels for a better understanding of biological pathways from gene to disease. These may include multi-modal imaging genetics studies, bioinformatic analyses, and functional analyses of cell and animal models.

Psychophysiological endophenotypes to characterize mechanisms of known schizophrenia genetic loci

BACKGROUND

Endophenotypes are laboratory-based measures hypothesized to lie in the causal chain between genes and clinical disorder, and to serve as a more powerful way to identify genes associated with the disorder. One promise of endophenotypes is that they may assist in elucidating the neurobehavioral mechanisms by which an associated genetic polymorphism affects disorder risk in complex traits. We evaluated this promise by testing the extent to which variants discovered to be associated with schizophrenia through large-scale meta-analysis show associations with psychophysiological endophenotypes.

METHOD

We genome-wide genotyped and imputed 4905 individuals. Of these, 1837 were whole-genome-sequenced at 11× depth. In a community-based sample, we conducted targeted tests of variants within schizophrenia-associated loci, as well as genome-wide polygenic tests of association, with 17 psychophysiological endophenotypes including acoustic startle response and affective startle modulation, antisaccade, multiple frequencies of resting electroencephalogram (EEG), electrodermal activity and P300 event-related potential.

RESULTS

Using single variant tests and gene-based tests we found suggestive evidence for an association between contactin 4 (CNTN4) and antisaccade and P300. We were unable to find any other variant or gene within the 108 schizophrenia loci significantly associated with any of our 17 endophenotypes. Polygenic risk scores indexing genetic vulnerability to schizophrenia were not related to any of the psychophysiological endophenotypes after correction for multiple testing.

CONCLUSIONS

The results indicate significant difficulty in using psychophysiological endophenotypes to characterize the genetically influenced neurobehavioral mechanisms by which risk loci identified in genome-wide association studies affect disorder risk.

Sensation seeking and impulsive traits as personality endophenotypes for antisocial behavior: Evidence from two independent samples

Sensation seeking and impulsivity are personality traits that are correlated with risk for antisocial behavior (ASB). This paper uses two independent samples of twins to (a) test the extent to which sensation seeking and impulsivity statistically mediate genetic influence on ASB, and (b) compare this to genetic influences accounted for by other personality traits. In Sample 1, delinquent behavior, as well as impulsivity, sensation seeking and Big Five personality traits, were measured in adolescent twins from the Texas Twin Project. In Sample 2, adult twins from the Australian Twin Registry responded to questionnaires that assessed individual differences in Eysenck's and Cloninger's personality dimensions, and a structured telephone interview that asked participants to retrospectively report DSM-defined symptoms of conduct disorder. Bivariate quantitative genetic models were used to identify genetic overlap between personality traits and ASB. Across both samples, novelty/sensation seeking and impulsive traits accounted for larger portions of genetic variance in ASB than other personality traits. We discuss whether sensation seeking and impulsive personality are causal endophenotypes for ASB, or merely index genetic liability for ASB.

Displacement behaviors in chimpanzees (Pan troglodytes): A neurogenomics investigation of the RDoC Negative Valence Systems domain

The current study aimed to systematically investigate genetic and neuroanatomical correlates of individual variation in scratching behaviors, a well-validated animal-behavioral indicator of negative emotional states with clear links to the NIMH Research Domain Criteria (RDoC) response to potential harm ("anxiety") construct within the Negative Valence Systems domain. Utilizing data from a sample of 76 captive chimpanzees (Pan troglodytes), we (a) examined the association between scratching and presence or absence of the RS3-containing DupB element in the AVPR1A 5' flanking region, (b) utilized voxel-based morphometry (VBM) to identify gray matter (GM) voxel clusters that differentiated AVPR1A genotype, and (c) conducted a VBM-guided voxel-of-interest analysis to examine the association between GM intensity and scratching. AVPR1A evidenced sexually dimorphic associations with scratching. VBM analyses revealed significant differences in GM by genotype across twelve clusters largely in the frontal cortex. Regions differentiating AVPR1A genotype showed sex-specific associations with scratching. Results suggest that sexually dimorphic associations between AVPR1A and scratching may be explained by genotype-specific neuroanatomical variation. The current study provides an example of the way in which chimpanzee research is uniquely poised for multilevel, systematic investigations of psychopathology-relevant constructs within the context of the RDoC framework.

Endophenotype best practices

This review examines the current state of electrophysiological endophenotype research and recommends best practices that are based on knowledge gleaned from the last decade of molecular genetic research with complex traits. Endophenotype research is being oversold for its potential to help discover psychopathology relevant genes using the types of small samples feasible for electrophysiological research. This is largely because the genetic architecture of endophenotypes appears to be very much like that of behavioral traits and disorders: they are complex, influenced by many variants (e.g., tens of thousands) within many genes, each contributing a very small effect. Out of over 40 electrophysiological endophenotypes covered by our review, only resting heart, a measure that has received scant advocacy as an endophenotype, emerges as an electrophysiological variable with verified associations with molecular genetic variants. To move the field forward, investigations designed to discover novel variants associated with endophenotypes will need extremely large samples best obtained by forming consortia and sharing data obtained from genome wide arrays. In addition, endophenotype research can benefit from successful molecular genetic studies of psychopathology by examining the degree to which these verified psychopathology-relevant variants are also associated with an endophenotype, and by using knowledge about the functional significance of these variants to generate new endophenotypes. Even without molecular genetic associations, endophenotypes still have value in studying the development of disorders in unaffected individuals at high genetic risk, constructing animal models, and gaining insight into neural mechanisms that are relevant to clinical disorder.

What can time-frequency and phase coherence measures tell us about the genetic basis of P3 amplitude?

In a recent comprehensive investigation, we largely failed to identify significant genetic markers associated with P3 amplitude or to corroborate previous associations between P3 and specific single nucleotide polymorphisms (SNPs) or genes. In the present study we extended this line of investigation to examine time-frequency (TF) activity and intertrial phase coherence (ITPC) in the P3 time window, both of which are associated with P3 amplitude. Previous genome-wide research has reported associations between P3-related theta and delta activity and individual genetic variants. A large, population-based sample of 4211 subjects, comprising male and female adolescent twins and their parents, was genotyped for 527,828 single nucleotide polymorphisms (SNPs), from which over six million SNPs were accurately imputed. Heritability estimates were greater for TF energy than ITPC, whether based on biometric models or the combined influence of all measured SNPs (derived from genome-wide complex trait analysis). The magnitude of overlap in the specific SNPs associated with delta energy and ITPC and P3 amplitude was significant. A genome-wide analysis of all SNPs, accompanied by an analysis of approximately 17,600 genes, indicated a region of chromosome 2 around TEKT4 that was significantly associated with theta ITPC. Analysis of candidate SNPs and genes previously reported to be associated with P3 or related phenotypes yielded one association surviving correction for multiple tests: between theta energy and CRHR1. However, we did not obtain significant associations for SNPs implicated in previous genome-wide studies of TF measures. Identifying specific genetic variants associated with P3 amplitude remains a challenge.

Electrophysiological Endophenotypes for Schizophrenia

Endophenotypes are quantitative, heritable traits that may help to elucidate the pathophysiologic mechanisms underlying complex disease syndromes, such as schizophrenia. They can be assessed at numerous levels of analysis; here, we review electrophysiological endophenotypes that have shown promise in helping us understand schizophrenia from a more mechanistic point of view. For each endophenotype, we describe typical experimental procedures, reliability, heritability, and reported gene and neurobiological associations. We discuss recent findings regarding the genetic architecture of specific electrophysiological endophenotypes, as well as converging evidence from EEG studies implicating disrupted balance of glutamatergic signaling and GABAergic inhibition in the pathophysiology of schizophrenia. We conclude that refining the measurement of electrophysiological endophenotypes, expanding genetic association studies, and integrating data sets are important next steps for understanding the mechanisms that connect identified genetic risk loci for schizophrenia to the disease phenotype.

Heritability and molecular-genetic basis of the P3 event-related brain potential: a genome-wide association study

P3 amplitude is a candidate endophenotype for disinhibitory psychopathology, psychosis, and other disorders. The present study is a comprehensive analysis of the behavioral- and molecular-genetic basis of P3 amplitude and a P3 genetic factor score in a large community sample (N = 4,211) of adolescent twins and their parents, genotyped for 527,829 single nucleotide polymorphisms (SNPs). Biometric models indicated that as much as 65% of the variance in each measure was due to additive genes. All SNPs in aggregate accounted for approximately 40% to 50% of the heritable variance. However, analyses of individual SNPs did not yield any significant associations. Analyses of individual genes did not confirm previous associations between P3 amplitude and candidate genes but did yield a novel association with myelin expression factor 2 (MYEF2). Main effects of individual variants may be too small to be detected by GWAS without larger samples.

Heritability and molecular-genetic basis of resting EEG activity: a genome-wide association study

Several EEG parameters are potential endophenotypes for different psychiatric disorders. The present study consists of a comprehensive behavioral- and molecular-genetic analysis of such parameters in a large community sample (N = 4,026) of adolescent twins and their parents, genotyped for 527,829 single nucleotide polymorphisms (SNPs). Biometric heritability estimates ranged from .49 to .85, with a median of .78. The additive effect of all SNPs (SNP heritability) varied across electrodes. Although individual SNPs were not significantly associated with EEG parameters, several genes were associated with delta power. We also obtained an association between the GABRA2 gene and beta power (p < .014), consistent with findings reported by others, although this did not survive Bonferroni correction. If EEG parameters conform to a largely polygenic model of inheritance, larger sample sizes will be required to detect individual variants reliably.

Heritability and molecular genetic basis of electrodermal activity: a genome-wide association study

The molecular genetic basis of electrodermal activity (EDA) was analyzed using 527,829 single nucleotide polymorphisms (SNPs) in a large population-representative sample of twins and parents (N = 4,424) in relation to various EDA indices. Biometric analyses suggested that approximately 50% or more of variance in all EDA indices was heritable. The combined effect of all SNPs together accounted for a significant amount of variance in each index, affirming their polygenic basis and heritability. However, none of the SNPs were genome-wide significant for any EDA index. Previously reported SNP associations with disorders such as substance dependence or schizophrenia, which have been linked to EDA abnormalities, were not significant; nor were associations between EDA and genes in specific neurotransmitter systems. These results suggest that EDA is influenced by multiple genes rather than by polymorphisms with large effects.

The Power of Theory, Research Design, and Transdisciplinary Integration in Moving Psychopathology Forward

While the past few decades have seen much work in psychopathology research that has yielded provocative insights, relatively little progress has been made in understanding the etiology of mental disorders. We contend that this is due to an overreliance on statistics and technology with insufficient attention to adequacy of experimental design, a lack of integration of data across various domains of research, and testing of theoretical models using relatively weak study designs. We provide a conceptual discussion of these issues and follow with a concrete demonstration of our proposed solution. Using two different disorders - depression and substance use - as examples, we illustrate how we can evaluate competing theories regarding their etiology by integrating information from various domains including latent variable models, neurobiology, and quasi-experimental data such as twin and adoption studies, rather than relying on any single methodology alone. More broadly, we discuss the extent to which such integrative thinking allows for inferences about the etiology of mental disorders, rather than focusing on descriptive correlates alone. Greater scientific insight will require stringent tests of competing theories and a deeper conceptual understanding of the advantages and pitfalls of methodologies and criteria we use in our studies.

Heritability and molecular genetic basis of antisaccade eye tracking error rate: a genome-wide association study

Antisaccade deficits reflect abnormalities in executive function linked to various disorders including schizophrenia, externalizing psychopathology, and neurological conditions. We examined the genetic bases of antisaccade error in a sample of community-based twins and parents (N = 4,469). Biometric models showed that about half of the variance in the antisaccade response was due to genetic factors and half due to nonshared environmental factors. Molecular genetic analyses supported these results, showing that the heritability accounted for by common molecular genetic variants approximated biometric estimates. Genome-wide analyses revealed several SNPs as well as two genes-B3GNT7 and NCL-on Chromosome 2 associated with antisaccade error. SNPs and genes hypothesized to be associated with antisaccade error based on prior work, although generating some suggestive findings for MIR137, GRM8, and CACNG2, could not be confirmed.

Heritability and molecular genetic basis of acoustic startle eye blink and affectively modulated startle response: a genome-wide association study

Acoustic startle responses have been studied extensively in relation to individual differences and psychopathology. We examined three indices of the blink response in a picture-viewing paradigm-overall startle magnitude across all picture types, and aversive and pleasant modulation scores-in 3,323 twins and parents. Biometric models and molecular genetic analyses showed that half the variance in overall startle was due to additive genetic effects. No single nucleotide polymorphism was genome-wide significant, but GRIK3 produced a significant effect when examined as part of a candidate gene set. In contrast, emotion modulation scores showed little evidence of heritability in either biometric or molecular genetic analyses. However, in a genome-wide scan, PARP14 produced a significant effect for aversive modulation. We conclude that, although overall startle retains potential as an endophenotype, emotion-modulated startle does not.

Genetic associations of nonsynonymous exonic variants with psychophysiological endophenotypes

We mapped ∼85,000 rare nonsynonymous exonic single nucleotide polymorphisms (SNPs) to 17 psychophysiological endophenotypes in 4,905 individuals, including antisaccade eye movements, resting EEG, P300 amplitude, electrodermal activity, affect-modulated startle eye blink. Nonsynonymous SNPs are predicted to directly change or disrupt proteins encoded by genes and are expected to have significant biological consequences. Most such variants are rare, and new technologies can efficiently assay them on a large scale. We assayed 247,870 mostly rare SNPs on an Illumina exome array. Approximately 85,000 of the SNPs were polymorphic, rare (MAF < .05), and nonsynonymous. Single variant association tests identified a SNP in the PARD3 gene associated with theta resting EEG power. The sequence kernel association test, a gene-based test, identified a gene PNPLA7 associated with pleasant difference startle, the difference in startle magnitude between pleasant and neutral images. No other single nonsynonymous variant, or gene-based group of variants, was strongly associated with any endophenotype.

In search of rare variants: preliminary results from whole genome sequencing of 1,325 individuals with psychophysiological endophenotypes

Whole genome sequencing was completed on 1,325 individuals from 602 families, identifying 27 million autosomal variants. Genetic association tests were conducted for those individuals who had been assessed for one or more of 17 endophenotypes (N range = 802-1,185). No significant associations were found. These 27 million variants were then imputed into the full sample of individuals with psychophysiological data (N range = 3,088-4,469) and again tested for associations with the 17 endophenotypes. No association was significant. Using a gene-based variable threshold burden test of nonsynonymous variants, we obtained five significant associations. These findings are preliminary and call for additional analysis of this rich sample. We argue that larger samples, alternative study designs, and additional bioinformatics approaches will be necessary to discover associations between these endophenotypes and genomic variation.