Structure of principal component analysis (PCA) of P300 amplitudes across electrodes showing factor loadings of P300 amplitudes at each electrode position on the temporoparietal and frontal P300 amplitude principal components and the associated eigenvalues and percentage of explained variance of each component. The PCAs did not include the 2 frontopolar electrodes Fp1 and Fp2 owing to their liability for artifact contamination. Using Scree-tests, an eigenvalue of 1 served as the minimum criteria for the determination of the number of factors. The PCAs were computed for the combined sample of schizophrenic patients, healthy siblings, and controls and for each group separately to obtain information about group differences of the factor structure. Because factor structure was similar between groups (data available on request), all subsequent analyses were performed based on the factor scores (weighted sum of standardized amplitude or latency values, weight = factor loading) from the combined sample (n = 270), which is depicted in this figure.
Oddball event-related potential group grand averages (targets). In schizophrenic patients, strongly reduced P300 (P3b) amplitudes are seen in the temporoparietal area, whereas in unaffected siblings, only slightly reduced P3b amplitudes are observed (strongest in the right posterior temporal region). In siblings, a double-peak configuration is seen, with an earlier peak between 250 and 300 milliseconds in the centroparietal region, which is less obvious in controls and absent in schizophrenic patients. In the frontal regions (F3 and F4) of schizophrenic patients and siblings, the slope of the P300 amplitude is shifted toward more positive values, and the P300 is generally less well expressed, as in the temporoparietal region. In schizophrenic patients, amplitudes at electrode positions Fp1 and Fp2 are possibly contaminated by eye movement artifacts, whereas there is no evidence of such contamination in siblings.
Group comparisons of mean frontal and temporoparietal P300 amplitude principal components in schizophrenic patients(n = 66), siblings (n = 115), and controls (n = 89). Error bars represent SEs. Multivariate analysis of covariance: Rao R4,514 = 3.94;P= .004.
Intraclass correlation coefficient (ICCU) values68 of temporoparietal (ICCU = 0.31; P= .006) and frontal (ICCU = 0.23; P= .03) P300 amplitude principal components in age- and sex-matched sib-pairs (n = 57). For comparison, ICCU values in age- and sex-matched controls (n = 57) (matched to the siblings) are also depicted. The temporoparietal P300 group difference between siblings and controls is less obvious in this smaller sample than in the large group comparison in Figure 3
Winterer G, Egan MF, Raedler T, Sanchez C, Jones DW, Coppola R, Weinberger DR. P300 and Genetic Risk for Schizophrenia. Arch Gen Psychiatry. 2003;60(11):1158-1167. doi:10.1001/archpsyc.60.11.1158
Copyright 2003 American Medical Association. All Rights Reserved. Applicable FARS/DFARS Restrictions Apply to Government Use.2003
We assessed the suitability of event-related potential frontal and temporoparietal P300 changes as intermediate phenotypes in genetic studies of schizophrenia. We applied a principal component analysis approach based on the notion that P300 abnormalities in siblings of schizophrenic patients may involve a widespread network of relatively weak cortical generators and because an earlier, smaller study that used a topographic analysis of covariance model did not show that localized P300 changes predict risk for schizophrenia.
P300 changes in 66 schizophrenic patients, 115 healthy siblings of schizophrenic patients, and 89 unrelated controls were studied during a standard auditory oddball paradigm. Principal components were calculated across electrodes, revealing frontal and temporoparietal components for latency and amplitude, respectively. For the frontal and temporoparietal P300 amplitude and latency components, the intraclass correlations (ICCs) between sib-pairs (pairs of unaffected siblings and schizophrenic index patients) and the relative risk ratios (λ) were determined.
Compared with controls, schizophrenic patients and their unaffected siblings showed significant reductions in the temporoparietal P300 amplitude component. Both groups were also characterized by a significantly higher frontal P300 amplitude component. Significant ICCs and increased relative risk ratios were found for the frontal (ICCU = 0.18; P =.04; λ = 3.4) and temporoparietal (ICCU = 0.24; P = .01; λ = 1.7) P300 amplitude components.
Temporoparietal P300 amplitude reduction and frontal P300 amplitude increase seem to be quantitative phenotypes associated with increased risk of schizophrenia. Both measures may be useful for increasing the statistical power of genetic studies of schizophrenia.
SCHIZOPHRENIA IS considered a heritable disorder with a complex genetic architecture interacting with environmental factors. Family-based linkage studies1 with noncoding markers spanning the genome have identified a few minor susceptibility loci that have been difficult to replicate. Because of this genetic complexity, positional cloning of these loci may be a daunting task, requiring large study samples.2 The goal of reducing the genetic complexity is the rationale for studying biological traits as intermediate phenotypes.3- 6 One example of biological traits is event-related potentials (ERPs).7 Event-related potentials tend to be relatively stable and heritable characteristics, and several studies have shown that the ERP P300, that is, its temporoparietal and frontal subcomponents, may be a useful electrophysiologic trait marker for schizophrenia.
The temporoparietal P300 (P3b) subcomponent has been thought to reflect the relative amount of resources allocated to processing an unpredicted stimulus.8- 10 P3b amplitude increases with greater expenditures of attentional resources11 and greater demand on memory updating.12 On the other hand, the frontal P300 subcomponent is considered to embody the cortical orienting response, and the amplitude increases with perceptual discrimination difficulties.13- 16
Temporoparietal and frontal P300 subcomponents show medium to high test-retest reliability for amplitude and low to medium test-retest reliability for latency.17- 20 Based on twin studies,7 heritability of P300 (P3b) measures was estimated to be 0.6 to 0.8, with latencies usually showing lower values than amplitudes. Heritability studies21 of the frontal P300 subcomponent suggest a heritability similar to that of the temporoparietal P3b subcomponent.
A few studies22- 26 have described P3b latency prolongation in schizophrenic patients. However, most studies report P3b amplitude reduction, which is relatively independent of medication status,23- 25,27- 36 suggesting that P3b amplitude reduction may reflect a trait feature of schizophrenia.19,37,38 Topographic and source analysis studies19,24,39- 48 have revealed that the P3b amplitude reduction in schizophrenia is most strongly expressed bilaterally over the temporoparietal areas, with some, but not all, studies showing a stronger left side amplitude reduction.
P3b amplitude reductions also have been described in nonschizophrenic family members of schizophrenic patients, including discordant twins49 and healthy siblings.50- 54 A few, but not all, studies46,52,55- 59 also reported delayed P3b latencies in family members of schizophrenic patients. However, several studies55- 59 did not replicate the finding of a decreased temporoparietal P300 amplitude reduction in family members of schizophrenic patients, including our own earlier study47 of a smaller subsample of the present study population. Most of these earlier studies were relatively underpowered because of small sample sizes and did not consider the possibility that the P300 data might reflect dysfunction in a disturbed network of weak source generators. The present P300 analysis of schizophrenic patients and their siblings differs from earlier studies in that the larger sample size now available made it possible to calculate intraclass correlations (ICCs) between sib-pairs (pairs of unaffected siblings and schizophrenic index patients) and relative risk ratios (λ) for P300. In addition, in the present study we applied a principal component analysis (PCA) model instead of a topographic analysis of variance model.
Data were collected between July 1, 1998, and October 31, 2000, as part of the CBDB Sibling Study. The details of participant recruitment, evaluation, and potential ascertainment biases are discussed elsewhere.60 Briefly, participants had to be aged 18 to 60 years, have a premorbid IQ score greater than 70, and be able to give informed consent. Applicants with alcohol or drug abuse in the past 6 months, dependence in the past year, or more than a 5-year history of abuse or dependence were excluded. Preliminary diagnosis was established by a research psychiatric social worker through review of psychiatric records and interviews with family members. Schizophrenic outpatients and their siblings were recruited from local and national sources as described elsewhere.60 Control subjects were recruited through the National Institutes of Health Normal Volunteer Office or from the community and were screened according to the same criteria, with the additional requirement that they not have a first-degree relative with schizophrenia. Participants were interviewed and videotaped by a research psychiatrist (masked to control or sibling status) using the Structured Clinical Interview for DSM-IV (SCID and SCID-II).61 All available psychiatric records were reviewed. A second psychiatrist reviewed all Structured Clinical Interview for DSM-IV and family data; disagreements with respect to diagnosis were reviewed by a third psychiatrist, who was the "tie breaker" (M.F.E.). Psychiatric diagnosis was established on the basis of DSM-IV criteria. All procedures were approved by the National Institute of Mental Health institutional review board.
All index patients fulfilled DSM-IV criteria for schizophrenia. Patients with schizoaffective disorder (n = 9) and individuals with fewer than 30 artifact-free electroencephalographic (EEG) trials (n = 13) were excluded, as were schizophrenic patients with an additional Axis I diagnosis or any other neurological or medical diagnosis. Excluded individuals are not listed in Table 1. Siblings of schizophrenic patients were only included in group comparisons if they were considered unaffected, that is, that they had no history of psychotic illness, including schizophrenia, schizoaffective disorder, mood disorder with psychotic symptoms, or any other disorder that at some point was accompanied by psychotic symptoms. Accordingly, affected siblings (n = 16) were excluded and are not listed in Table 1. Fourteen unaffected siblings who fulfilled the criteria for any present "nonpsychotic" psychiatric disorder such as major depression or personality disorder were required to be clinically stable. Thirty-nine unaffected siblings had a history of a nonpsychotic psychiatric disorder but did not fulfill diagnostic criteria at the time of testing. These individuals were included because we were interested in trait characteristics independent of the illness-related state.
Participants were presented 250 tones of 2 different pitches by loudspeaker at an approximately 80-dB sound pressure level. Stimulus duration was 50 milliseconds, and rise and decay were each 5 milliseconds. The pseudorandomized interstimulus interval was 1.0 to 1.5 seconds with all stimuli presented in pseudorandomized order (target: 1500 Hz; probability, 20%; nontarget: 1000 Hz; probability, 80%). Participants were asked to count the number of targets.
A more detailed description of ERP recording is given elsewhere.47 Using a Grass model (8-24D) EEG, P300 was recorded with eyes closed using gold electrodes.62,63 Electrode positions were defined according to the international 10/20 system, and impedance was kept below 5 kΩ. Eye movements were recorded on an additional channel across electrodes 1 cm lateral and either above or below the corner of the eye (data available on paper recordings only). The EEG was measured from 16 electrode positions (Fp1, Fp2, F3, F4, F7, F8, T3, T4, T5, T6, C3, C4, P3, P4, O1, and O2). Amplifiers were calibrated using a 50-µV square pulse. Sampling frequency was 256 Hz using 0.5- to 70.0-Hz filters.62,63 All channels were recorded against separate ears as a reference, that is, left hemispheric electrodes were referenced to the left ear and right hemispheric electrodes were referenced to the right ear. Before analysis of data, offline artifact detection was performed automatically (amplitude criterion, >70 µV), and it was manually edited afterward by a rater (G.W.) who was masked to diagnosis. Specifically, all EEG segments that contained obvious eye movement, head movement, or muscle artifacts were excluded. Individuals (n = 13) with fewer than 30 artifact-free target trials were excluded from further investigation. This procedure resulted in significant group differences with respect to the mean ± SD number of artifact-free target trials (controls vs schizophrenic patients: 46.5 ± 0.38 vs 48.1 ± 0.77; t95 = 1.99; P = .07; and controls vs siblings: 46.5 ± 0.38 vs 48.0 ± 0.56; t197 = 1.9; P = .02).
Electrophysiologic test-retest data were obtained from 17 individuals(11 controls, 1 unaffected sibling, and 5 schizophrenic inpatients), with a mean ± SD duration between sessions of 316.1 ± 66.9 days.
Generation of ERP grand averages across the montage and subsequent topographic amplitude and latency analysis were performed automatically using EEGSYS.64 P300 peak amplitudes against baseline and latencies were determined for the target condition of the ERPs by locating the peak within the specified time window for P300 (260-420 milliseconds after stimulus). Baseline was defined as the available 50 milliseconds before stimulus EEG.
Statistical analyses were performed using Statistica.65 Group comparisons of demographic and clinical data were performed using analysis of covariance (ANCOVA), the χ2 statistic, or the t test. All multivariate ANCOVAs (MANCOVAs) and ANCOVAs of neurophysiologic data were performed with the 2 factors diagnosis and sex; age and number of ERP trials were the covariates, and the P300 measures were the dependent variables if not otherwise indicated.
Descriptive MANCOVAs of P300 latencies and amplitudes were performed across electrode positions, with either P300 amplitudes or latencies at all 16 electrode positions as dependent variables (although these 16 variables are not completely independent of each other).66,67 These analyses were followed by post hoc honestly significance difference analyses for unequal participants per group (n) at selected electrode positions.
To maximize the ability to measure distributed and correlated small effects, for reasons of data reduction and to increase measurement stability, exploratory PCAs were performed. The PCAs (varimax rotation) were performed separately on P300 latencies and P300 amplitudes across electrode positions. Similar approaches have been applied in previous studies20,59 to analyze topographic P300 data from schizophrenic patients, siblings, and controls(Figure 1). Test-retest stability of principal components was calculated for principal components with ICCU.68 This was achieved by using factor loadings from the PCA as weight factors. These weights were multiplied by the raw amplitude or latency values at each electrode position and were subsequently summed for each individual. The resulting individual values for the test and retest set were then subjected to ICC calculations.
Group comparisons of the temporoparietal and frontal P300 amplitude and latency principal components were performed using MANCOVA and were followed by planned comparisons. These analyses were performed (1) under inclusion of all clinically unaffected siblings and schizophrenic index patients (multiple sibships) and (2) with only 1 randomly selected sibling from each family. Additional exploratory post hoc analyses were performed using the Tukey honest significance difference test. To test whether the observed effects were dependent on family membership, random-effects 2-way ANCOVAs (mixed model) were also performed, with family membership as the random factor under inclusion of those families with data available from at least 2 family members (unaffected siblings and schizophrenic patients).
Similarity between sib-pairs of selected phenotypes (ie, individual factor scores) was obtained by treating the electrophysiologic measures as quantitative variables and calculating ICCU values.68 These analyses were performed (1) for 1 randomly selected sib-pair for family and(2) for sib-pairs (1 per family) that were matched for sex and age as far as possible. For the matched sib-pairs, we also performed repeated-measures MANCOVAs with subsequent planned comparisons and the Pearson product moment correlation coefficient r to test the consistency of within-family slopes.69 Relative risk (λ) calculations were performed for the obtained ERP phenotypes based on "qualitative phenotype" definitions analogous to those previously reported for spectroscopic, neuropsychologic, and resting EEG data.60,70,71 Accordingly, a cutoff value for abnormality (ie, for the qualitative phenotype) was arbitrarily delimited based on the arithmetic mean of control subjects and subtracting 1 SD. Relative risk values (λ) were calculated for the phenotype definition of 1 SD under inclusion of 1 randomly selected unaffected sibling per family.72,73 A relative risk greater than 2.0 is considered a reasonable indication that a trait is familial and, except for the contribution of shared environmental factors, potentially heritable. A relative risk greater than 5.0 is considered high.5
No significant age differences were found between healthy controls and schizophrenic patients (t153 = 1.2; P = .23) and between controls and unaffected siblings (t201 = 1.6; P = .12)(Table 1). Sex differed significantly between controls and schizophrenic patients (χ21 =15.4; P<.001) and between controls and siblings(χ21 = 5.45; P = .02). Controls had more years of education than schizophrenic patients (t116 = 5.18; P<.001) but not than unaffected siblings (t168 = 0.33; P = .74). Controls were less frequently smokers than schizophrenic patients (χ21 = 16.3; P<.001) but not siblings (χ21 = 0.6; P = .44). With respect to handedness as measured by the Edinburgh Scale,74 no significant differences were found between controls and schizophrenic patients (χ21 = 0.1; P = .08) and between controls and siblings (χ21 = 2.1; P = .15). Comparing the percentage of correctly counted targets, controls did not differ significantly from schizophrenic patients (t80 = 0.53; P = .60) and siblings (t110 = −1.02; P = .31).
Figure 2 shows unadjusted grand averages of the raw ERPs across the montage. Adjusted 2-way MANCOVA analysis with diagnosis (controls, siblings, and schizophrenic patients) and sex as factors and P300 peak amplitudes across electrodes as dependent variables revealed a main effect of diagnosis (Rao R32,490 = 2.98; P<.001) and sex (Rao R16,245 = 2.90; P<.001) and a significant interaction between diagnosis and sex (Rao R32,490 = 1.71; P = .01). The equivalent multivariate analysis of variance for P300 peak latencies did not show significant main effects of diagnosis (Rao R32,490 =0.80; P = .77) or sex (Rao R16,245 = 0.68; P = .81). Results from additional post hoc group comparisons of P300 amplitudes at 2 frontal and 2 temporal electrode positions are given in Table 2. The numerical decrease of temporoparietal P300 amplitudes and the numerical increase of frontal P300 amplitudes is more obvious in the MANCOVA model than in the raw grand averages (Figure 2) after having adjusted for age, number of trials, and sex (2-way ANCOVA). However, in siblings, the P300 amplitude differences mostly do not achieve statistical significance. Overall, the findings at frontal and temporal electrode sites hardly differ from a previous study using a similar statistical approach of a subsample of the present cohort.47
The PCA of the P300 amplitudes across electrodes reveals a 2-factor structure in the combined sample comprising controls, siblings, and schizophrenic patients: a frontal P300 subcomponent and a temporoparietal P300 subcomponent(P3b) (Figure 1). The PCA of P300 latencies across electrode positions (data not shown) shows an almost identical 2-factor solution as found for the P300 amplitudes. The larger principal latency component is a temporoparietal P3b latency component, which includes temporal, parietal, and occipital electrode positions (eigenvalue, 6.1; variance, 43.5%). The smaller principal component is a frontal P300 latency component (eigenvalue, 2.2; variance, 16.0%).
Test-retest stability of the obtained principal components is low for the frontal latency component (ICCU = 0.25; P = .07) and temporoparietal latency component (ICCU = −0.23; P = .82), moderate for the temporoparietal amplitude component (ICCU = 0.61; P = .001), and high for the frontal amplitude component (ICCU = 0.88; P<.001).
Group comparisons based on the principal component values are depicted in Figure 3. The reader should recall that the factor scores (P300 amplitude principal component values) can be negative because these values were standardized. Overall MANCOVA of schizophrenic patients, siblings, and controls for the temporoparietal and frontal P300 amplitude principal components resulted in a significant main effect of diagnosis (Rao R4,514 = 3.94; P = .004) (Figure 3) and sex (Rao R2,257 = 11.16; P<.001). No significant diagnosis × sex interaction was seen (Rao R4,514 = 0.86; P = .49). The significance level increased when including the control group and randomly selecting only 2 individuals from each family (1 index patient and 1 clinically unaffected sibling): we observed a significant main effect of diagnosis (Rao R4,388 = 3.84; P = .005), a significant main effect of sex (Rao R2,194 = 9.04; P<.001), and, again, no significant interaction (Rao R4,388 = 0.86; P = .49). Finally, the MANCOVA of schizophrenic patients with diagnostic subtype (paranoid, undifferentiated) and sex as independent factors showed a significant effect of diagnosis (Rao R2,57 = 5.24; P = .008) and sex (Rao R2,57 = 6.23, P = .004) but no significant interaction (Rao R2,57 = 2.03; P = .14). An equivalent MANCOVA of siblings with the 2 subgroups (nonpsychotic, psychiatric lifetime diagnosis[DSM-IV Axis I and II] present vs not present) revealed no significant main effect of subgroup (Rao R2,100 = 1.80; P = .17) or sex (Rao R2,100 = 0.68; P = .51) and no interaction (Rao R2,100 = 1.05; P = .35). The MANCOVA of patients with neuroleptic drug treatment vs no medication showed no significant main effects of medication use (Rao R2,56 = 0.36; P = .70) or sex (Rao R2,56 = 2.03; P = .14) and no interaction (Rao R2,56 = 1.05; P = .86). Likewise, schizo phrenic patients with and without a family history of schizophrenia did not show significantly different P300 amplitude principal component values (Rao R2,54 = 1.72; P = .19), although a statistical trend interaction with sex was observed (Rao R2,54 = 2.96; P = .06).
Planned comparisons of the temporoparietal P300 amplitude principal components between groups under inclusion of all participants within the framework of the previous overall MANCOVA are shown in Figure 3. Similar results were obtained when randomly selecting only 1 sibling from each family, with siblings showing significantly lower temporoparietal P300 principal component amplitude mean values than controls(−0.02 vs 0.24; P = .01) and higher mean values than schizophrenic patients (−0.02 vs −0.29; P = .02). The temporoparietal P300 amplitude principal component values also differed significantly by sex, with males showing lower mean amplitudes in the whole group (−0.20 vs 0.18; P<.001), which was also seen when randomly selecting only 2 family members (−0.23 vs 0.18; P<.001). Random-effects 2-way ANCOVA results were consistent with the results of the previous analyses. A significant effect was found for diagnosis with respect to the dependent variable temporoparietal P300 amplitude component (F1,28.3 = 19.77; P<.001), without a significant effect of family membership (F59,2.1 = 3.54; P = .24).
Figure 3 also shows planned comparisons of frontal P300 amplitude principal component values within the framework of the previous overall MANCOVA under inclusion of all participants. When only 1 unaffected sibling and schizophrenic patient was randomly selected from each family, similar results were obtained. Siblings showed higher frontal P300 amplitude principal component mean values than controls (0.18 vs −0.23; P = .004) but did not significantly differ from schizophrenic patients (P = .60). No significant sex differences in frontal P300 amplitude principal component mean values were seen between males and females in the whole group (−0.03 vs −0.005; P = .85) or in the group with no more than 2 participants per family(0.07 vs −0.05; P = .32). The additional random-effects 2-way ANCOVA results were again consistent with the results of the previous analyses. No significant effect was seen for the frontal P300 principal component(F1,40.7 = 0.00; P = .95), again without a significant effect of family membership (F59,12.9 = 1.44; P = .24).
Relative risk values (λ) of unaffected siblings (1 randomly selected unaffected sibling per family) are presented for the principal component values given in Table 3. Separate relative risk calculations for males (7 of 29 siblings and 10 of 49 controls with the phenotype) and females (10 of 51 siblings and 2 of 40 controls with the phenotype) revealed an increased risk in females (λ = 3.92) but not males (λ= 1.18) with respect to the temporoparietal P300 amplitude principal component. With regard to the frontal P300 amplitude principal component, increased risk values were found in males (λ = 4.39, 95% confidence interval, 1.7-11.1) and females (λ = 2.94, 95% confidence interval, 1.1-8.2).
Sib-pair ICCU values (1 randomly selected sib-pair per family) are also presented in Table 3. Equivalent ICC analyses for P300 amplitude values at selected electrode positions(F3, F4, T5, and T6) were not significant (data not shown). Repeated-measures MANCOVAs, ICC calculations, and Pearson r values were calculated for sib-pairs, that is, 1 pair per family, with matching of sib-pairs for age and sex as much as possible (a separate analysis for male- or female-only sib-pairs was not possible because of small numbers). Repeated-measures MANCOVAs with frontal and temporoparietal P300 amplitude principal components as dependent variables resulted in significant effects of diagnosis (Rao R2,48 = 7.35; P = .002) and sex (Rao R6,96 = 3.26; P = .006), without any evidence of interaction (Rao R6,96 = 0.43; P = .86). Subsequent planned comparisons of matched siblings and index patients were consistent with the previously presented larger group comparisons (Figure 3), revealing a significant difference in temporoparietal P300 amplitude between schizophrenic patients and their related siblings (P<.001) and no significant difference in frontal P300 amplitude between the 2 groups (P = .74). When performing the same ICC calculations as depicted in Table 3 for the matched sib-pairs (Figure 4), ICCU values increased for the temporoparietal and frontal P300 amplitude components. For the temporoparietal P300 amplitude component, the Pearson r value was high (r = 0.531; P<.001), with a parallel slope for most sib-pairs (Figure 4). The Pearson r value for the frontal P300 amplitude component was considerably lower, with trend significance (r =0.282; P = .07). The lower correlation of the frontal component is reflected by the relatively steep slopes in those sib-pairs with the highest values, although, for the most part, schizophrenic patients with high frontal amplitude are related to unaffected siblings who also have relatively high amplitude values and vice versa.
The present study was performed to further explore the suitability of the ERP P300 as an intermediate phenotype for genetic studies of schizophrenia. To our knowledge, this is the first study to report similarity (ICC) and relative risk (λ) estimates of ERP P300 measures in clinically unaffected siblings of schizophrenic patients. This allows inferences about the association of P300 abnormalities and the genetic risk for schizophrenia, especially because electrophysiologic variables seem to be little affected by shared environment.7
The detection of statistically significant group differences between siblings and controls in the present and a previous study47 of a subsample of this investigation was unsuccessful when applying an analysis of variance model that only takes into account local information from single electrodes as dependent variables. Of note, any demonstration of group differences may have been difficult because of a possible ascertainment bias toward relatively healthy siblings.47,60 Therefore, we decided to use a PCA, which combines the weighted, correlated effects across the montage in a data-driven way and thus might offer greater statistical power. In this context, it is important to mention that intracortical investigations of P300 and source analyses of scalp-recorded data clearly indicate that frontal and temporoparietal P300 is composed of multiple, cortically widespread, and relatively weak cortical generators.16,41,48,75 Also, source analyses of scalp-recorded P300 have provided some evidence that frontal and temporoparietal P300 changes in schizophrenic patients involve an extended cortical "P300 network" rather than only affecting isolated P300 generators, for example, at the temporoparietal junction.48 Previous, equivalent PCAs of P300 have resulted in stable and neurophysiologically plausible frontal and temporoparietal components,20,59 which provides validity to this particular approach. The obtained simple 2-factor structure in the present study is, to a large extent, in accordance with these previous factor analyses.
Four potential phenotypes were analyzed: the temporoparietal and frontal P300 amplitude and latency principal components. Largely in accordance with previous literature,17- 20 the 2 latency components showed relatively low test-retest stability, whereas stability was moderate to high for the 2 amplitude components.
Based on earlier studies,47,49- 59 it was predicted that a reduction in the temporoparietal P300 amplitude is associated with increased risk of schizophrenia. The present study is in agreement with this prediction. Relative risk for reduced P3b amplitude is increased in clinically unaffected siblings, and group comparisons revealed that schizophrenic patients have lower P3b amplitudes than controls, with their siblings showing intermediate values. The random-effects analysis suggests that this finding can be generalized and is not explained by family membership. The similarity of P3b amplitudes between sib-pairs seems to be in the low to moderate range, which provides indirect evidence that P3b amplitude is heritable. The P3b amplitude distribution across groups suggests that reduced P3b amplitude is a quantitative rather than a qualitative risk factor, although a more conclusive answer to this question would have also required analysis of control sib-pairs with no schizophrenic relative in their family. Schizophrenic patients generally have lower amplitudes than their related sibling; however, many clinically affected-unaffected sib-pairs still have P300 amplitudes that extend well into the range of controls. Therefore, and because the relative risk (λ) value for temporoparietal P300 reduction in siblings is relatively low, this phenotype may have some value for increasing statistical power in large genetic studies but is obviously of limited value with respect to the prediction of the individual risk for schizophrenic illness. Also, sex-specific differences need to be taken into account.
Previous studies of frontal P300 amplitudes have found unchanged or slightly lower frontal P300 amplitudes in schizophrenic patients and siblings.23,24,36,39,40,43,45- 48,58,76- 78 However, in the present study, increased amplitudes were found in both groups. When exclusively regarding schizophrenic patients, the possibility of artifacts as an explanation cannot be discarded because an increased rate of eye movements is a frequent observation in schizophrenic patients and because the P300 amplitudes are highest over the frontopolar electrodes. In siblings, however, P300 amplitude differences are strongest over the frontal and not frontopolar electrodes, suggesting that the increased P300 amplitude is generated in the frontal cortex. An important observation in this context is that the ICCs between P300 amplitude principal components of clinically affected-unaffected sib-pairs are significant, which suggests that the increased P300 amplitude in schizophrenic patients is at least in part also cortically generated. Another possible explanation for the difference between our results and those of some earlier studies of frontal P300 may involve our method for measuring frontal P300. In the present study, an automatic peak detection procedure was used instead of manual peak detection, which usually has been used in previous and mostly smaller studies. This technical difference might be important because although the temporoparietal subcomponent of the P300 is a true ERP with a clear peak in the ERP tracings and cortical oscillations in selective frequency bands with phase locking and polarity reversal, the poorly expressed frontal P300 subcomponent can be considered to be a correlate of cortically generated, neuronal noise with some embedded signal-like activity.14,16,48,75 When peak detection of frontal P300 is performed manually, it is in many cases impossible to be objective about the real peak. When using automatic peak detection, the signal peak also can be missed because of the noisy structure of this subcomponent. Accordingly, it is conceivable that we have monitored with our approach an increase in cortically generated noise in patients and clinically unaffected siblings.
Another possibility that could explain the frontal P300 amplitude difference between studies has to do with statistical power to detect differences. Moreover, and in contrast to most earlier studies, we used a PCA, which is sensitive in the detection of relatively small but correlated activity changes. With respect to schizophrenic patients, it also might be relevant that most earlier studies of frontal P300 amplitude studied inpatients, whereas our study only included clinically stable outpatients. This could be important because it was shown that a frontal amplitude increase correlates with symptom improvement.19,26,36
Siblings of patients with high frontal P300 amplitudes also have relatively high amplitudes, and our findings show that familiality plays a role because relative risk is increased in the moderate range and ICCs are significant. As with the temporoparietal P300, the amplitude distribution suggests a quantitative rather than a qualitative trait. However, the steep slopes of frontal P300 amplitudes between sib-pairs within the upper portion of the amplitude factor score range suggest that as-yet-unidentified neurophysiologic mechanisms that might not directly be trait dependent also may play a role. One possible explanation could be that the high amplitude in some siblings and schizophrenic patients reflects a higher engagement of the frontal cortex to compensate for a cortical deficit either within the frontal cortex or elsewhere, for example, in the temporal lobe. The plausibility of this explanation would lie in the neurophysiologic observation that the orientation-related, frontal P300 amplitude component assists in probing widespread cortical areas for rapid evaluation of the biological significance of stimuli.14,16
Corresponding author and reprints: Georg Winterer, MD, Clinical Brain Disorders Branch, National Institute of Mental Health, National Institutes of Health, Building 10, Room 4s229A, Mail Stop Code 1379, Bethesda, MD 20892(e-mail: email@example.com).
Submitted for publication January 31, 2002; final revision received August 21, 2002; accepted November 18, 2002.