Context
Identifying susceptibility genes for schizophrenia may be complicated by phenotypic heterogeneity, with some evidence suggesting that phenotypic heterogeneity reflects genetic heterogeneity.
Objective
To evaluate the heritability and conduct genetic linkage analyses of empirically derived, clinically homogeneous schizophrenia subtypes.
Design
Latent class and linkage analysis.
Setting
Taiwanese field research centers.
Participants
The latent class analysis included 1236 Han Chinese individuals with DSM-IV schizophrenia. These individuals were members of a large affected-sibling-pair sample of schizophrenia (606 ascertained families), original linkage analyses of which detected a maximum logarithm of odds (LOD) of 1.8 (z = 2.88) on chromosome 10q22.3.
Main Outcome Measures
Multipoint exponential LOD scores by latent class assignment and parametric heterogeneity LOD scores.
Results
Latent class analyses identified 4 classes, with 2 demonstrating familial aggregation. The first (LC2) described a group with severe negative symptoms, disorganization, and pronounced functional impairment, resembling “deficit schizophrenia.” The second (LC3) described a group with minimal functional impairment, mild or absent negative symptoms, and low disorganization. Using the negative/deficit subtype, we detected genome-wide significant linkage to 1q23-25 (LOD = 3.78, empiric genome-wide P = .01). This region was not detected using the DSM-IV schizophrenia diagnosis, but has been strongly implicated in schizophrenia pathogenesis by previous linkage and association studies.Variants in the 1q region may specifically increase risk for a negative/deficit schizophrenia subtype. Alternatively, these results may reflect increased familiality/heritability of the negative class, the presence of multiple 1q schizophrenia risk genes, or a pleiotropic 1q risk locus or loci, with stronger genotype-phenotype correlation with negative/deficit symptoms. Using the second familial latent class, we identified nominally significant linkage to the original 10q peak region.
Conclusion
Genetic analyses of heritable, homogeneous phenotypes may improve the power of linkage and association studies of schizophrenia and thus have relevance to the design and analysis of genome-wide association studies.
Schizophrenia (SZ) is a debilitating psychiatric illness for which the pathogenesis remains unclear. Although strongly implicating a number of genetic regions and promising candidate genes, genetic analyses have also highlighted the causal complexity of the disorder. Genomic variants and regions showing association or linkage often differ between studies, and statistical evidence is frequently unconvincing. Confirmed functional variants remain elusive. The presence of small-effect risk variants, genetic heterogeneity, incomplete penetrance, and phenocopies renders the phenotype immensely challenging for genetic dissection.
Genetic studies of SZ may also be confounded by clinical heterogeneity. Schizophrenia has no established biological markers, and diagnoses are based on the presence and duration of symptoms, the prominence of which can vary considerably between patients. Debate has long centered on whether this phenotypic heterogeneity reflects causal heterogeneity. Although the use of operationalized diagnoses has unarguably facilitated substantial progress in SZ research, it also assumes a single causal model for patients with distinct clinical profiles. If this clinical variation reflects the presence of different causal processes, the use of more homogeneous clinical phenotypes may facilitate the identification of genotype-phenotype correlations.1 Such approaches have yielded promising findings in other complex disorders, including attention-deficit/hyperactivity disorder2 and migraine.3
Some evidence suggests that the genetic cause of SZ may differ between patients with distinct clinical profiles. Kendler and colleagues4 demonstrated significant clinical differences between linked and unlinked families for a candidate region in 8p22-21. Also, in a large association study of the D-amino acid oxidase activator gene (DAOA) with SZ, significant association was detected in a subset of subjects (n = 112) who had experienced major mood episodes but not in the (much larger) subset (n = 597) who had not, and no association was detected by using the operationalized SZ diagnosis.5 However, few studies have used a fully data-driven approach to empirically derive clinical subtypes for genetic analyses of SZ; to our knowledge, there is just one. Fanous and colleagues6 recently reported linkage analyses of psychotic subtypes derived by means of latent class analysis (LCA) of 755 psychotic subjects from 270 Irish families. They detected 4 regions that showed suggestive linkage to 1 of the LCA subtypes but little linkage to the traditional clinical diagnosis.
Interestingly, 1 of the LCA subgroups identified in the study by Fanous et al resembled “deficit SZ” (DS), a well-characterized SZ subtype supported by taxometric analysis7 and showing good longitudinal stability.8 Compared with patients without DS, patients with DS exhibit primary, enduring negative symptoms9; poorer social functioning; less depression; less suicidal ideation; and fewer delusions of an exclusively social content while having similar overall severity of all delusions, hallucinations, and formal thought disorder.10 Family studies suggest utility of the DS subtype in genetic studies, showing significant correlation for deficit vs nondeficit subtypes in sibling pairs concordant for SZ11 and increased family risk of SZ in individuals with DS compared with patients without DS.12 However, few studies have reported genetic analyses of DS. Bakker and colleagues13 reported association of DS with the PIP5K2A gene on chromosome 10p12, whereas Fanous and colleagues6 identified suggestive linkage of their deficit class to a region in chromosome 20.
To further explore the utility of homogeneous subtypes in identifying SZ risk loci, we conducted genome-wide linkage analyses of empirically derived SZ subtypes in 1236 Taiwanese individuals with SZ.14 These individuals were members of the largest extant affected-sibling-pair sample for SZ, comprising 606 ascertained families. Initial genome-wide linkage analysis of the DSM-IV SZ diagnosis provided support for regions in 10q22.3, 2q14.1, 1p31.1, 15q14, and 4q21.23.14 However, no region achieved genome-wide significance, with a region in 10q22.3 showing the strongest evidence of linkage (z = 2.88, logarithm of odds [LOD] = 1.8). We aimed to investigate whether additional regions or increased significance could be detected by means of empirically derived clinical phenotypes.
This data set was accessed via the National Institute of Mental Health Center for Collaborative Genetic Studies on Mental Disorders (distribution 5.0) (http://zork.wustl.edu/nimh). The collection of this sample has been previously described.14 Briefly, affected-sibling-pair families containing at least 2 siblings with SZ were recruited from 6 centers in Taiwan. Each affected sibling underwent a diagnostic screen using supplemental medical records and a semistructured interview, followed by the Mandarin Chinese version of the Diagnostic Interview for Genetic Studies (DIGS).15,16 The DIGS interview data were supplemented with information from medical records and a semistructured interview with family members.17,18 Best-estimate final diagnoses were assigned by 2 research psychiatrists, with a third diagnostician resolving any disagreements. The sample was genotyped by the Center for Inherited Disease Research for 386 microsatellite markers spaced at 9-cM intervals (average) genome-wide.
The number and composition of distinct clinical groups were empirically determined by means of LCA, a statistical method for identifying subtypes of related cases from multivariate categorical data.19 On the basis of clinical experience, we selected 11 symptoms recorded during the DIGS interview that were deemed to capture the important aspects of clinical variation (Table 1). The LCA was conducted with Latent Gold software (version 4.5; Statistical Innovations Inc, Belmont, Massachusetts) and included 1236 Han Chinese participants with a DSM-IV diagnosis of SZ. Latent cluster models specifying from 1 to 10 latent classes were fitted, allowing up to 5000 iterations of the Expectation Maximization algorithm with a convergence criterion of 1 × 10−6. To avoid local maxima in parameter estimation, models were refitted multiple times with the use of different starting values. The best-fitting model was selected on the basis of the Bayesian information criterion, as previously recommended.20 The preferred model was that with the lowest Bayesian information criterion, providing optimal balance between fit and parsimony. In addition, we assessed the clinical interpretability of the identified classes and ensured that classes were of nontrivial size. Because many cases contained missing values for at least 1 indicator variable and to use all available data, cases containing missing values were included in the analysis. However, the impact of missing values was reduced by constructing variables from multiple, related DIGS items, providing a degree of redundancy. Latent Gold handles missing indicator values directly in the likelihood function by basing the likelihood contribution of each case on the observed indicators only. That is, parameters are estimated by using all available information for each of the cases.
An important assumption of LCA is that of local independence between observed indicator variables. That is, conditional on latent class membership, the indicator variables are statistically independent (uncorrelated). In cases in which the local independence assumption does not hold, LCA models may be modified to provide adequate fit.21 After initial model selection, we used the Latent Gold software to identify indicator variables with significant residual association, indicated by a bivariate residual (BVR) score greater than 3.84 (equivalent to P < .05 because the BVR statistic is distributed as χ21). Residual association was accounted for by incorporating direct effect parameters between significantly correlated variable pairs,22 starting with the highest BVR and proceeding iteratively until all BVRs were less than 3.84.
Familial aggregation analyses
Before performing genetic analyses of LCA subtypes, we sought evidence of a genetic contribution to latent class group membership, which should produce aggregation of group membership within families. On the basis of all 1236 affected individuals with SZ, the sample contained 758 possible affected relative pairs (ARPs). The proportion of ARPs concordant for membership in any latent class was compared with the expected value, which was derived by means of the marginal probabilities for each latent class. The expected and observed proportions were treated as binomial variables and compared by means of a normal approximation. Because the sample contained nonindependent relative pairs from families containing 3 or more affected individuals, the analysis was repeated with the use of an average of the number of concordant relative pairs within each family, in which each family was weighted to contribute a single relative pair.
On the basis of the posterior probabilities of latent class membership, affected individuals were assigned to their most likely latent class. Genome-wide linkage analyses were conducted for the 2 latent classes exhibiting familial aggregation. For each, individuals endorsing that class were considered affected and multipoint exponential LOD scores23 were calculated at 1-cM increments by means of the Spairs24 sharing statistic implemented in the Merlin program.25 Marker positions were determined by means of the sex-averaged Marshfield genetic map,26 as for the original analysis. This was considered the primary linkage analysis.
As a secondary analysis, we also calculated parametric heterogeneity LOD (HLOD) scores in the entire, original data set (606 families). Multipoint HLOD scores were calculated under simple dominant and recessive models, specifying disease allele frequencies of 0.01 and 0.1 and genotype penetrances of (0, 0.5, 0.5) and (0, 0, 0.5) for the dominant and recessive models, respectively, for 0, 1, or 2 copies of the disease allele.27,28
To estimate genome-wide significance levels for our primary analysis, we used an empiric method.29 For each of the 2 latent class groups used in genetic analyses, Merlin25 was used to generate 5000 data replicates under the null hypothesis of no linkage, which were analyzed at 1-cM increments by means of the exponential LOD score.23 We then derived the empiric distribution of the maximum LOD score at each position (LODmax) by recording, for each replicate, the highest LOD score obtained across the 2 latent class groups. This generated an empiric distribution of the highest score at each position that naturally accounts for both multiple testing and the correlation between statistics. The genome-wide significance level associated with the highest observed score was defined as the frequency of equivalent or higher peaks in the empiric distribution of LODmax. Linkage peaks separated by more than 40 cM were considered independent.30,31 Genome-wide significance and suggestive thresholds were defined as scores occurring with a probability of 0.05 and 1 in every LODmax replicate, respectively,32 and nominal P = .01 and P = .05 thresholds were defined as scores occurring with a probability of 0.01 and 0.05 at any map position.
The LCA included 1236 individuals with SZ (1175 affected siblings and 61 affected parents), with a mean (SD) age of 36.29 (9.58) years and a mean age at illness onset of 23.05 (6.92) years. Males composed 61.5% of the sample. The affected individuals were contained within 578 families, 561 of which contained 2 or more affected sibling pairs (28 of the original 606 families had no diagnostic and/or symptom data available and 17 had data for 1 sibling only). Of the 561 multiple-sibling families, 524 families contained 2 affected siblings, 36 contained 3 affected siblings, and 1 family contained 4 affected siblings. Genotype data were available for all individuals included in the LCA.
After fitting a single-class solution, a continuous improvement in fit was observed up to a 4-class solution; the relevant Bayesian information criterion values were 20463.27 (1 class), 19527.27 (2 classes), 19425.43 (3 classes), 19370.09 (4 classes), and 19401.08 (5 classes). To evaluate the robustness of our LCA model, we repeated the primary Latent Gold analysis by using the LCA procedure of SAS. Identical results were obtained, showing that our model was not sensitive to the particular software package used. We also repeated the Latent Gold analysis 20 times, each time randomly selecting 90% of the original sample for inclusion. Each run identified the same 4-class solution, with individuals being classified into the same latent class with 96.8% consistency across the 20 runs.
Assessment of Latent Gold's BVR output for the 4-class solution showed several BVRs greater than 3.84, suggesting a number of pairwise correlations between indicator variables remaining unexplained by the model. To ensure conditional independence of the LCA solution, we used the Latent Gold software to iteratively incorporate direct effects between all variable pairs with a BVR greater than 3.84. The results of the final, 4-class solution are shown in Table 2 and Figure 1.
Visual comparison of the 4 cluster profiles suggests qualitative differences between the identified phenotypic subgroups. A profile analysis of variance statistically validated these differences, rejecting the hypothesis of parallelism of the 4 plots (F57,1628 = 85.52, P < .001). This suggests an underlying phenotypic structure of 4 distinct categories rather than a single continuum of liability, which would produce classes differing primarily in symptom severity.33 The 4 classes were most differentiated by the levels of disorganization, functional impairment, negative symptoms (alogia and affective flattening), and overall symptom pattern. Negative symptoms best differentiated the 4 groups, with endorsement probabilities differing by approximately 90% between the least and most severely affected groups. All 4 classes exhibited high levels of delusions and hallucinations.
The largest subgroup, LC1 (γ = 0.43, n = 533) was not distinguished by marked presence or absence of any particular symptom, rather demonstrating moderate levels of delusions, hallucinations, disorganization, functional impairment, suicidality, and negative symptoms. The pattern of symptoms for this class was most likely to be a mixture of positive and negative symptoms.
The other 3 subgroups were more symptomatically distinctive. The second largest class, LC2 (γ = 0.29, n = 370) described a subtype characterized by moderate to severe negative symptoms (flat affect and alogia), prominent disorganization (speech and behavior), and marked to severe functional impairment. This class described a category of severely affected patients with SZ resembling DS. Of the 4 classes, individuals in LC2 were the least likely to experience delusions and hallucinations, although these symptoms were still present in a majority of individuals. Individuals in LC2 were also the least likely to have attempted suicide.
Of the 4 classes, subjects in LC3 (γ = 0.16, n = 213) were the least severely affected, demonstrating the highest level of global functioning, a virtual absence of negative symptoms, and the lowest levels of disorganization. This group had a pattern of symptoms characterized by predominantly positive symptoms when ill, but relatively few symptoms during remissions.
The fourth and smallest class (LC4; γ = 0.12, n = 120) displayed prominent delusions and hallucinations, moderate disorganization, and moderate to marked functional impairment. Individuals in this class were most likely to experience continuous, positive symptoms and were also the most likely to have attempted suicide. Negative symptoms were absent or mild for almost all individuals in this group.
Familial aggregation analyses
To justify linkage analysis of the identified latent classes, it was important to demonstrate familial aggregation of class membership, indicative of an underlying genetic influence. On the basis of the marginal probabilities of LC1, LC2, LC3, and LC4 (0.43, 0.29, 0.16, and 0.12, respectively), the expected number (proportion) of total concordant ARPs was 232 (0.31). The observed number (proportion) was 257 (0.339) (z = 1.97, P = .02), suggesting marginally significant overall familial aggregation. However, this aggregation was nonsignificant when each family was weighted to contribute a single relative pair (z = 1.38, P = .08). Among individual classes, a significant excess of concordant relative pairs was observed for only 2 classes, LC2 and LC3. For these classes, ARP class concordance was elevated on the basis of counts of both total relative pairs (LC2: z = 2.27, 1-sided P = .01; LC3: z = 2, P = .02) and independent pairs per family (counted as R − 1, where R is the number of affected individuals in the family) (LC2: z = 3.26, 1-sided P < .001; LC3: z = 2.58, P = .005). The excess was nonsignificant for LC1, and LC4 had fewer observed than expected concordant relative pairs. This supports a role of familial factors in the etiology of LC2 and LC3.
Incorporating correction for testing both LC2 and LC3, empiric genome-wide significant and suggestive linkage thresholds of 3.14 and 1.86 and nominal P = .01 and P = .05 thresholds of 1.45 and 0.84 were determined. Linkage analyses of LC2 included 257 individuals contained within 64 families, incorporating 67 total or 66 independent sibling pairs. Analyses of LC3 included 93 individuals contained within 23 families, with 25 total and 24 independent sibling pairs. Genome-wide multipoint LOD scores23 for the LC2 and LC3 subsets are shown in Figure 2.
For LC2, the strongest linkage evidence was observed at 1q23-25, 189 cM from the p-telomere (LOD = 3.78). This result easily surpassed the threshold for genome-wide significant linkage, with only 61 peaks of 3.78 or greater occurring in 5000 simulations of LODmax (genome-wide P = .01, incorporating correction for testing both LC2 and LC3). The 1-LOD drop (approximating a 95% confidence interval for the location of the peak) delimited a 15-cM region from 183 to 198 cM (approximately 166.6-180.6 megabases, National Center for Biotechnology Information Build 36.2). To explore the potential mode of inheritance of the putative chromosome 1 locus, we calculated multipoint HLOD scores under simple dominant and recessive models (see the “Methods” section for model parameters). These exploratory analyses demonstrated greater linkage evidence under a recessive (peak HLOD = 3.539) rather than a dominant (peak HLOD = 2.403) model.
Single-point LOD scores greater than 1 were observed for 4 markers in the 1q region: D1S1679 (LOD = 1.07), D1S1619 (LOD = 2.72), D1S1589 (LOD = 1.41), and D1S518 (LOD = 1.38), supporting the multipoint results. In the LC2 subset, an additional 2 regions achieved nominal P ≤ .01 and a total of 6 regions achieved nominal P ≤ .05 (Table 3).
For LC3, no region achieved genome-wide significant linkage, suggestive linkage, or nominal P ≤ .01. Ten regions achieved nominal P ≤ .05 (Table 3). One of these was located on 10q22.1 (peak LOD = 1.19, 94 cM from p-ter) and overlapped the maximum linkage peak detected in original analyses of these data using the DSM-IV SZ diagnosis (peak LOD = 1.8, z = 2.88, 100.9 cM from p-ter).
Parametric analyses of the original (entire) data set detected the strongest linkage evidence on 10q22.1 (HLOD-dominant = 2.402), in the same region as the original peak z score. Nine other peaks achieved HLOD scores greater than 1 under either the dominant or recessive model (Table 4). The 1q region detected with the use of LC2 showed little linkage evidence in the entire sample (peak HLOD = 0.738; recessive model, 216 cM from p-ter).
This is the second linkage study of SZ to use categorical phenotypes empirically derived from observed symptom profiles. We proposed that, if clinical heterogeneity reflects genetic heterogeneity, homogeneous clinical subtypes may allow the identification of genomic regions that the traditional diagnosis did not.
Our LCA identified a familial SZ subtype with prominent negative symptoms. The characteristics of this class showed striking similarity to a subtype also described by Fanous et al.6 In addition to prominent negative symptoms, individuals in both classes demonstrated high levels of social dysfunction and lower-than-average probabilities of delusions, hallucinations, and Schneiderian psychotic symptoms. Although we did not include depressive and manic items in our LCA, independent assessment of their probabilities in our negative subtype (data not shown) showed values lower than the sample median, consistent with Fanous et al.6 Individuals in our negative subtype also demonstrated prominent disorganization and low suicidality; however, these symptoms were not included in the Fanous et al study.
This negative subtype resembles DS. However, an important characteristic of DS is the presence of primary negative symptoms as opposed to those arising secondary to factors such as depression, paranoia, or preoccupation with psychotic symptoms.34 The available data did not allow us to distinguish between primary and secondary negative symptoms in our deficit (LC2) subset. Although LC2 demonstrated several features consistent with DS—including lower levels of depression,35 suicidality, and social/occupational function36 in the presence of similar (lower) overall levels of delusions and hallucinations—it also showed high levels of disorganization (incorporating formal thought disorder), which may influence the severity of negative symptoms. Thus, despite these similarities, our limitation in addressing the primary/secondary distinction precludes us from describing our LC2 subset as classic DS.
Linkage analyses of our negative subtype (LC2) identified a region of genome-wide significant linkage in 1q23-25. Both the original analysis and our secondary parametric heterogeneity analyses failed to detect linkage to this region, despite the latter's putative allowance for locus heterogeneity. This region has previously shown highly significant linkage to SZ in a Canadian sample37,38 and has also shown linkage to SZ in 2 Chinese studies39,40 and a European study.41 Furthermore, this region has been implicated in SZ susceptibility in 2 genome scan meta-analyses of SZ.42,43 Our results may suggest that variants in 1q23-25 specifically increase the risk for a negative/deficit subtype of SZ. In this case, although our LC2 subset was modest in size (64 families), its high phenotypic homogeneity may have increased the proportion of 1q-linked families and power to detect linkage. Alternatively, our detection of the 1q region may have resulted less from the specific phenotypic features of the LC2 subset than our use of a familial/heritable SZ subtype. Indeed, Fanous et al6 did not identify linkage of their deficit subtype to this 1q region, but rather detected linkage to a region in chromosome 20. However, these differences may reflect the presence of further genetic heterogeneity within negative/deficit subtypes or statistical variation relating to power, given the modest number of deficit sibling pairs included in our study (67 pairs) and the Fanous study (32 pairs).
Involvement of the 1q23-25 region in a negative/deficit SZ subtype may also be specific to the Han Chinese population or have resulted from this sample's ascertainment characteristics. A striking feature of this Taiwanese sample was a low rate of comorbid substance use/abuse, with just 5.2% and 3.7% of affected siblings reporting alcohol or drug abuse/dependence, respectively.44 In our negative class, the equivalent rates were 3.1% and 0.8%, which may have further increased phenotypic homogeneity and power to resolve SZ risk loci. Further research is required to clarify the structure, heritability, and relevant risk loci for negative/deficit subtypes of SZ and to determine the generalizability of our findings to other samples and populations.
Several genes in the 1q region have shown statistical association with SZ, although none of the evidence is compelling. Perhaps the best studied is the regulator of G-protein signaling 4 gene (RGS4), which is located in 1q23.3 and has shown modest association in multiple studies45-50 and one meta-analysis,51 although negative results have also been reported.52-54 Intriguingly, one recent study reported association of an RGS4 single-nucleotide polymorphism with severity of negative symptoms and neurocognitive performance,55 although another study reported negative association of DS with RGS4 single-nucleotide polymorphisms.13 Other genes in the 1q region showing some level of association with SZ include CAPON,56,57CHRNB2,58Cx50,59 and UHMK1.60,61 Interestingly, 2 recent, prominent studies detected association of SZ to a large, recurrent microdeletion in 1q21.1.62,63 However, the deletion is rare (frequency, <1% in cases) and thus likely accounts for little linkage to the 1q region. Other studies have shown association of SZ with a heterochromatin variant64 and a fragile site65 in 1q21. Comprehensive research in larger data sets—including genotype-phenotype analyses—may help to identify genetic variants in the 1q region involved in SZ susceptibility.
The other class demonstrating familiality in this study (LC3) described the least severely affected subgroup, demonstrating high global functioning, low disorganization, and a virtual absence of negative symptoms. Linkage analyses of LC3 detected evidence of linkage (nominal P ≤ .05) to a region in 10q overlapping the primary linkage peak detected in original analyses of these data.14 However, any phenotypic interpretation of these results should be undertaken cautiously because of the modest number of concordant sibling pairs included in the linkage analysis of LC3 (25 total pairs).
Our results should be interpreted in light of several considerations. First, the latent class typology we report is dependent on the response variables we selected for the LCA and the statistical method we used for selecting the best-fitting model. It remains to be shown whether these results differ with the use of different response variables, other rating instruments, or different diagnostic systems. The inclusion of additional symptoms may also result in the identification of further latent classes. Second, although the familial aggregation of LC2 and LC3 group membership suggests genetic influences, it could also reflect shared environmental effects. Twin studies could help to clarify the respective contributions of genetic and environmental effects to subtype liability. Third, the size and position of our linkage peaks may have been influenced by the number of sibling pairs concordant for each LCA subtype, with the number of total concordant sibling pairs for LC2 and LC3 being 67 and 25, respectively. It is possible that the detection of significant linkage with LC2, but not LC3, reflected power differences relating to sample size. However, the position of the 10q peak accords well with those reported in other studies, supporting the validity of its location.
In conclusion, our results suggest that the identification and analysis of homogeneous, heritable clinical subtypes may be an important avenue to progress in SZ genetics research. Reduced genetic heterogeneity in such subtypes may increase power to identify particular risk variants. Although genetic studies have achieved substantial progress by means of the traditional SZ diagnosis, exclusive reliance on DSM-IV diagnostic categories of SZ, albeit in ever-increasing samples, may not be sufficient for dissecting its genetic complexity. If the clinical SZ diagnosis comprises multiple, causally divergent subtypes or dimensions, additional value may be extracted from existing and future data sets by conducting empiric analyses of comprehensive phenotype data to identify homogeneous clinical subtypes and dimensions. The analysis of such alternative phenotypes may be an important, complementary approach for clarifying the genetic cause of SZ.
Correspondence: Bryan J. Mowry, MD, FRANZCP, Queensland Centre for Mental Health Research, The Park, Wacol, Queensland 4076, Australia (bryan_mowry@qcmhr.uq.edu.au).
Submitted for Publication: July 23, 2008; final revision received February 11, 2009; accepted February 27, 2009.
Author Contributions: Drs Holliday and Mowry had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.
Han Chinese Schizophrenia Linkage Study: Project leaders in Taiwan were Hai-Gwo Hwu, MD (Taiwan principal investigator, National Taiwan University Hospital), and Wei J. Chen, MD, SciD (Taiwan co–principal investigator). Other participants in Taiwan were Chih-Min Liu, MD, Shih-Kai Liu, MD, Ming-Hsien Shieh, MD, Tzung-Jeng Hwang, MD, MPH, Ming-Ming Tsuang, MD, Wen Chen OuYang, MD, PhD, Chun-Ying Chen, MD, Chwen-Cheng Chen, MD, PhD, Jin-Jia Lin, MD, Frank Huang-Chih Chou, MD, PhD, Ching-Mo Chueh, MD, Wei-Ming Liu, MD, Chiao-Chicy Chen, MD, Jia-Jiu Lo, MD, Jia-Fu Lee, MD, PhD, Seng Shen, MD, Yung Feng, MD, Shin-Pin Lin, MD, Shi-Chin Guo, MD, Ming-Cheng Kuo, MD, Liang-Jen Chuo, MD, Chih-Pin Lu, MD, Deng-Yi Chen, MD, Huan-Kwang Ferng, MD, Nan-Ying Chiu, MD, Wen-Kun Chen, MD, Tien-Cheng Lee, MD, Hsin-Pei Tang, MD, Yih-Dar Lee, MD, Wu-Shih Wang, MD, For-Wey Long, MD, PhD, Tiao-Lai Huang, MD, Jung-Kwang Wen, MD, Cheng-Sheng Chen, MD, Wen-Hsiang Huang, MD, Shu-Yu Yang, MD, Mei-Hua Hall, PhD, and Cheng-Hsing Chen, MD. Other participants in the United States were Stephen V. Faraone, PhD (co–principal investigator), Shao Zhu, MD (project coordinator), and Xingjia Cui, MD (project coordinator).
Financial Disclosure: None reported.
Funding/Support: This research was supported in part by Australian National Health and Medical Research Council (NHMRC) grants 143027 and 339454 (Dr Mowry). Dr Holliday was supported by an NHMRC Public Health Postgraduate Scholarship. Dr Nyholt is supported by NHMRC fellowship 339462 and was supported in part by NHMRC program grant 389892. Drs Holliday and Mowry and Mr McLean are supported by the Queensland Department of Health. The data from the Han Chinese Schizophrenia Linkage Study were collected with funding from grant R01 MH59624 from the US National Institute of Mental Health to Ming T. Tsuang, MD, PhD (principal investigator).
Additional Contributions: We appreciate the comments of 5 anonymous reviewers, which helped to considerably improve the manuscript.
1.Jablensky
A Subtyping schizophrenia: implications for genetic research.
Mol Psychiatry 2006;11
(9)
815- 836
PubMedGoogle Scholar 2.Todd
RDHuang
HSmalley
SLNelson
SFWillcutt
EGPennington
BFSmith
SDFaraone
SVNeuman
RJ Collaborative analysis of DRD4 and DAT genotypes in population-defined ADHD subtypes.
J Child Psychol Psychiatry 2005;46
(10)
1067- 1073
PubMedGoogle Scholar 3.Nyholt
DRMorley
KIFerreira
MAMedland
SEBoomsma
DIHeath
ACMerikangas
KRMontgomery
GWMartin
NG Genomewide significant linkage to migrainous headache on chromosome 5q21.
Am J Hum Genet 2005;77
(3)
500- 512
PubMedGoogle Scholar 4.Kendler
KSMyers
JMO’Neill
FAMartin
RMurphy
BMacLean
CJWalsh
DStraub
RE Clinical features of schizophrenia and linkage to chromosomes 5q, 6p, 8p, and 10p in the Irish Study of High-Density Schizophrenia Families.
Am J Psychiatry 2000;157
(3)
402- 408
PubMedGoogle Scholar 5.Williams
NMGreen
EKMacgregor
SDwyer
SNorton
NWilliams
HRaybould
RGrozeva
DHamshere
MZammit
SJones
LCardno
AKirov
GJones
IO’Donovan
MCOwen
MJCraddock
N Variation at the
DAOA/G30 locus influences susceptibility to major mood episodes but not psychosis in schizophrenia and bipolar disorder.
Arch Gen Psychiatry 2006;63
(4)
366- 373
PubMedGoogle Scholar 6.Fanous
AHNeale
MCWebb
BTStraub
REO’Neill
FAWalsh
DRiley
BPKendler
KS Novel linkage to chromosome 20p using latent classes of psychotic illness in 270 Irish high-density families.
Biol Psychiatry 2008;64
(2)
121- 127
PubMedGoogle Scholar 7.Blanchard
JJHoran
WPCollins
LM Examining the latent structure of negative symptoms: is there a distinct subtype of negative symptom schizophrenia?
Schizophr Res 2005;77
(2-3)
151- 165
PubMedGoogle Scholar 8.Amador
XFKirkpatrick
BBuchanan
RWCarpenter
WTMarcinko
LYale
SA Stability of the diagnosis of deficit syndrome in schizophrenia.
Am J Psychiatry 1999;156
(4)
637- 639
PubMedGoogle Scholar 9.Carpenter
WT
JrHeinrichs
DWWagman
AM Deficit and nondeficit forms of schizophrenia: the concept.
Am J Psychiatry 1988;145
(5)
578- 583
PubMedGoogle Scholar 10.Kirkpatrick
BRam
RBromet
E The deficit syndrome in the Suffolk County Mental Health Project.
Schizophr Res 1996;22
(2)
119- 126
PubMedGoogle Scholar 11.Ross
DEKirkpatrick
BKarkowski
LMStraub
REMacLean
CJO'Neill
FACompton
ADMurphy
BWalsh
DKendler
KS Sibling correlation of deficit syndrome in the Irish study of high-density schizophrenia families.
Am J Psychiatry 2000;157
(7)
1071- 1076
PubMedGoogle Scholar 12.Kirkpatrick
BRoss
DEWalsh
DKarkowski
LKendler
KS Family characteristics of deficit and nondeficit schizophrenia in the Roscommon Family Study.
Schizophr Res 2000;45
(1-2)
57- 64
PubMedGoogle Scholar 13.Bakker
SCHoogendoorn
MLHendriks
JVerzijlbergen
KCaron
SVerduijn
WSelten
JPPearson
PLKahn
RSSinke
RJ The
PIP5K2A and
RGS4 genes are differentially associated with deficit and non-deficit schizophrenia.
Genes Brain Behav 2007;6
(2)
113- 119
PubMedGoogle Scholar 14.Faraone
SVHwu
HGLiu
CMChen
WJTsuang
MMLiu
SKShieh
MHHwang
TJOu-Yang
WCChen
CYChen
CCLin
JJChou
FHChueh
CMLiu
WMHall
MHSu
JVan Eerdewegh
PTsuang
MT Genome scan of Han Chinese schizophrenia families from Taiwan: confirmation of linkage to 10q22.3.
Am J Psychiatry 2006;163
(10)
1760- 1766
PubMedGoogle Scholar 15.Nurnberger
JI
JrBlehar
MCKaufmann
CAYork-Cooler
CSimpson
SGHarkavy-Friedman
JSevere
JBMalaspina
DReich
T Diagnostic Interview for Genetic Studies: rationale, unique features, and training. NIMH Genetics Initiative.
Arch Gen Psychiatry 1994;51
(11)
849- 859
PubMedGoogle Scholar 16.Chen
WJLiu
SKChang
CJLien
YJChang
YHHwu
HG Sustained attention deficit and schizotypal personality features in nonpsychotic relatives of schizophrenic patients.
Am J Psychiatry 1998;155
(9)
1214- 1220
PubMedGoogle Scholar 17.Gershon
ESDeLisi
LEHamovit
JNurnberger
JI
JrMaxwell
MESchreiber
JDauphinais
DDingman
CW
IIGuroff
JJ A controlled family study of chronic psychoses: schizophrenia and schizoaffective disorder.
Arch Gen Psychiatry 1988;45
(4)
328- 336
PubMedGoogle Scholar 18.Maxwell
ME Family Interview for Genetic Studies (FIGS): A Manual for FIGS. Bethesda, MD Clinical Neurogenetics Branch, Intramural Research Program, National Institute of Mental Health1992;
19.Rindskopf
DRindskopf
W The value of latent class analysis in medical diagnosis.
Stat Med 1986;5
(1)
21- 27
PubMedGoogle Scholar 20.Sclove
S Application of model-selection criteria to some problems in multivariate analysis.
Psychometrika 1987;52
(3)
333- 343
Google Scholar 21.Magidson
JVermunt
JK Latent class models. Kaplan
D
The Sage Handbook of Quantitative Methodology for the Social Sciences. Thousand Oaks, CA Sage Publications2004;175- 198
Google Scholar 22.Hagenaars
JA Latent structure models with direct effects between indicators: local dependence models.
Sociol Methods Res 1988;16
(3)
379- 405
Google Scholar 23.Kong
ACox
NJ Allele-sharing models: LOD scores and accurate linkage tests.
Am J Hum Genet 1997;61
(5)
1179- 1188
PubMedGoogle Scholar 24.Whittemore
ASHalpern
J A class of tests for linkage using affected pedigree members.
Biometrics 1994;50
(1)
118- 127
PubMedGoogle Scholar 25.Abecasis
GRWigginton
JE Handling marker-marker linkage disequilibrium: pedigree analysis with clustered markers.
Am J Hum Genet 2005;77
(5)
754- 767
PubMedGoogle Scholar 26.Broman
KWMurray
JCSheffield
VCWhite
RLWeber
JL Comprehensive human genetic maps: individual and sex-specific variation in recombination.
Am J Hum Genet 1998;63
(3)
861- 869
PubMedGoogle Scholar 27.Greenberg
DAAbreu
PHodge
SE The power to detect linkage in complex disease by means of simple LOD-score analyses.
Am J Hum Genet 1998;63
(3)
870- 879
PubMedGoogle Scholar 28.Pal
DKDurner
MGreenberg
DA Effect of misspecification of gene frequency on the two-point LOD score.
Eur J Hum Genet 2001;9
(11)
855- 859
PubMedGoogle Scholar 29.Ferreira
MAO’Gorman
LLe Souëf
PBurton
PRToelle
BGRobertson
CFVisscher
PMMartin
NGDuffy
DL Robust estimation of experimentwise
P values applied to a genome scan of multiple asthma traits identifies a new region of significant linkage on chromosome 20q13.
Am J Hum Genet 2005;77
(6)
1075- 1085
PubMedGoogle Scholar 30.Ott
J Analysis of Human Genetic Linkage. 3rd ed. Baltimore, MD John Hopkins University Press1999;
31.Wiltshire
SCardon
LR McCarthy
MI Evaluating the results of genomewide linkage scans of complex traits by locus counting.
Am J Hum Genet 2002;71
(5)
1175- 1182
PubMedGoogle Scholar 32.Lander
EKruglyak
L Genetic dissection of complex traits: guidelines for interpreting and reporting linkage results.
Nat Genet 1995;11
(3)
241- 247
PubMedGoogle Scholar 33. McCutcheon
AL Latent Class Analysis. Beverly Hills, CA Sage Publications1987;
34.Kirkpatrick
BBuchanan
RWRoss
DECarpenter
WT
Jr A separate disease within the syndrome of schizophrenia.
Arch Gen Psychiatry 2001;58
(2)
165- 171
PubMedGoogle Scholar 35.Kirkpatrick
BBuchanan
RWBreier
ACarpenter
WT
Jr Depressive symptoms and the deficit syndrome of schizophrenia.
J Nerv Ment Dis 1994;182
(8)
452- 455
PubMedGoogle Scholar 36.Fenton
WS McGlashan
TH Antecedents, symptom progression, and long-term outcome of the deficit syndrome in schizophrenia.
Am J Psychiatry 1994;151
(3)
351- 356
PubMedGoogle Scholar 37.Brzustowicz
LMHodgkinson
KAChow
EWHoner
WGBassett
AS Location of a major susceptibility locus for familial schizophrenia on chromosome 1q21-q22.
Science 2000;288
(5466)
678- 682
PubMedGoogle Scholar 38.Logue
MWBrzustowicz
LMBassett
ASChow
EWVieland
VJ A posterior probability of linkage-based re-analysis of schizophrenia data yields evidence of linkage to chromosomes 1 and 17.
Hum Hered 2006;62
(1)
47- 54
PubMedGoogle Scholar 39.Hwu
HGLiu
CMFann
CSOu-Yang
WCLee
SF Linkage of schizophrenia with chromosome 1q loci in Taiwanese families.
Mol Psychiatry 2003;8
(4)
445- 452
PubMedGoogle Scholar 40.Zheng
YWang
XGu
NFeng
GZou
FQin
WZhang
JLin
WTao
RQian
XHe
L A two-stage linkage analysis of Chinese schizophrenia pedigrees in 10 target chromosomes.
Biochem Biophys Res Commun 2006;342
(4)
1049- 1057
PubMedGoogle Scholar 41.Shaw
SHKelly
MSmith
ABShields
GHopkins
PJLoftus
JLaval
SHVita
ADe Hert
MCardon
LRCrow
TJSherrington
RDeLisi
LE A genome-wide search for schizophrenia susceptibility genes.
Am J Med Genet 1998;81
(5)
364- 376
PubMedGoogle Scholar 42.Lewis
CMLevinson
DFWise
LHDeLisi
LEStraub
REHovatta
IWilliams
NMSchwab
SGPulver
AEFaraone
SVBrzustowicz
LMKaufmann
CAGarver
DLGurling
HMLindholm
ECoon
HMoises
HWByerley
WShaw
SHMesen
ASherrington
RO'Neill
FAWalsh
DKendler
KSEkelund
JPaunio
TLonnqvist
JPeltonen
LO’Donovan
MCOwen
MJWildenauer
DBMaier
WNestadt
GBlouin
JLAntonarakis
SEMowry
BJSilverman
JMCrowe
RRCloninger
CRTsuang
MTMalaspina
DHarkavy-Friedman
JMSvrakic
DMBassett
ASHolcomb
JKalsi
G McQuillin
ABrynjolfson
JSigmundsson
TPetursson
HJazin
EZoega
THelgason
T Genome scan meta-analysis of schizophrenia and bipolar disorder, part II: schizophrenia.
Am J Hum Genet 2003;73
(1)
34- 48
PubMedGoogle Scholar 43.Ng
MYMLevinson
DFFaraone
SVSuarez
BKDeLisi
LEArinami
TRiley
BPaunio
TPulver
AEIrmansyahHolmans
PAEscamilla
MWildenauer
DBWilliams
NMLaurent
CMowry
BJBrzustowicz
LMaziade
MSklar
PGarver
DLAbecasis
GLerer
BFallin
MDGurling
HMDGejman
PVLindholm
EMoises
HByerley
WWijsman
EMForabosco
PTsuang
MTHwu
H-GOkazaki
YKendler
KSWormley
BFanous
AWalsh
DO’Neill
APeltonen
LNestadt
GLasseter
VKLiang
KYPapadimitriou
GDikeos
DSchwab
SGOwen
MJO’Donovan
MCNorton
NHare
ERaventos
HNicolini
HAlbus
MMaier
WNimgaonkar
VLTerenius
LMallet
JJay
MGodard
SNertney
DAlexander
MCrowe
RRSilverman
JMBassett
ASRoy
M-AMérette
CPato
CNPato
MTLouw Roos
JKohn
YAmann-Zalcenstein
DKalsi
G McQuillin
ACurtis
DBrynjolfson
JSigmundsson
TPetursson
HSanders
ARDuan
JJazin
EMyles-Worsley
MKarayiorgou
MLewis
CM Meta-analysis of 32 genomewide linkage studies of schizophrenia [published online ahead of print December 30, 2008].
Mol Psychiatry 10.1038/mp.2008.135
Google Scholar 44.Hwu
HGFaraone
SVLiu
CMChen
WJLiu
SKShieh
MHHwang
TJTsuang
MMOuYang
WCChen
CYChen
CCLin
JJChou
FHChueh
CMLiu
WMHall
MH Tsuang
MT Taiwan schizophrenia linkage study: the field study.
Am J Med Genet B Neuropsychiatr Genet 2005;134B
(1)
30- 36
Google Scholar 45.Chowdari
KVMirnics
KSemwal
PWood
JLawrence
EBhatia
TDeshpande
SNB K
TFerrell
REMiddleton
FADevlin
BLevitt
PLewis
DANimgaonkar
VL Association and linkage analyses of
RGS4 polymorphisms in schizophrenia [published correction appears in
Hum Mol Genet. 2003;12(14):1781].
Hum Mol Genet 2002;11
(12)
1373- 1380
PubMedGoogle Scholar 46.Chen
XDunham
CKendler
SWang
XO’Neill
FAWalsh
DKendler
KS Regulator of G-protein signaling 4 (
RGS4) gene is associated with schizophrenia in Irish high density families.
Am J Med Genet B Neuropsychiatr Genet 2004;129B
(1)
23- 26
PubMedGoogle Scholar 47.Morris
DWRodgers
A McGhee
KASchwaiger
SScully
PQuinn
JMeagher
DWaddington
JLGill
MCorvin
AP Confirming
RGS4 as a susceptibility gene for schizophrenia.
Am J Med Genet B Neuropsychiatr Genet 2004;125B
(1)
50- 53
Google Scholar 48.Williams
NMPreece
ASpurlock
GNorton
NWilliams
HJ McCreadie
RGBuckland
PSharkey
VChowdari
KVZammit
SNimgaonkar
VKirov
GOwen
MJO’Donovan
MC Support for RGS4 as a susceptibility gene for schizophrenia.
Biol Psychiatry 2004;55
(2)
192- 195
PubMedGoogle Scholar 49.Zhang
FSt Clair
DLiu
XSun
XSham
PCCrombie
CMa
XWang
QMeng
HDeng
WYates
PHu
XWalker
NMurray
RMCollier
DALi
T Association analysis of the
RGS4 gene in Han Chinese and Scottish populations with schizophrenia.
Genes Brain Behav 2005;4
(7)
444- 448
PubMedGoogle Scholar 50.Fallin
MDLasseter
VKAvramopoulos
DNicodemus
KKWolyniec
PS McGrath
JASteel
GNestadt
GLiang
KYHuganir
RLValle
DPulver
AE Bipolar I disorder and schizophrenia: a 440-single-nucleotide polymorphism screen of 64 candidate genes among Ashkenazi Jewish case-parent trios.
Am J Hum Genet 2005;77
(6)
918- 936
PubMedGoogle Scholar 51.Talkowski
MESeltman
HBassett
ASBrzustowicz
LMChen
XChowdari
KVCollier
DACordeiro
QCorvin
APDeshpande
SNEgan
MFGill
MKendler
KSKirov
GHeston
LLLevitt
PLewis
DALi
TMirnics
KMorris
DWNorton
NO'Donovan
MCOwen
MJRichard
CSemwal
PSobell
JLSt Clair
DStraub
REThelma
BKVallada
HWeinberger
DRWilliams
NMWood
JZhang
FDevlin
BNimgaonkar
VL Evaluation of a susceptibility gene for schizophrenia: genotype based meta-analysis of RGS4 polymorphisms from thirteen independent samples.
Biol Psychiatry 2006;60
(2)
152- 162
PubMedGoogle Scholar 52.Guo
STang
WShi
YHuang
KXi
ZXu
YFeng
GHe
L
RGS4 polymorphisms and risk of schizophrenia: an association study in Han Chinese plus meta-analysis.
Neurosci Lett 2006;406
(1-2)
122- 127
PubMedGoogle Scholar 53.Ishiguro
HHoriuchi
YKoga
MInada
TIwata
NOzaki
NUjike
HMuratake
TSomeya
TArinami
T RGS4 is not a susceptibility gene for schizophrenia in Japanese: association study in a large case-control population.
Schizophr Res 2007;89
(1-3)
161- 164
PubMedGoogle Scholar 54.Rizig
MA McQuillin
APuri
VChoudhury
KDatta
SThirumalai
SLawrence
JQuested
DPimm
JBass
NLamb
GMoorey
HBadacsonyi
AKelly
KMorgan
JPunukollu
BKandasami
GCurtis
DGurling
H Failure to confirm genetic association between schizophrenia and markers on chromosome 1q23.3 in the region of the gene encoding the regulator of G-protein signaling 4 protein (
RGS4).
Am J Med Genet B Neuropsychiatr Genet 2006;141B
(3)
296- 300
PubMedGoogle Scholar 55.So
HCChen
RYChen
EYCheung
EFLi
TSham
PC An association study of RGS4 polymorphisms with clinical phenotypes of schizophrenia in a Chinese population.
Am J Med Genet B Neuropsychiatr Genet 2008;147B
(1)
77- 85
PubMedGoogle Scholar 56.Brzustowicz
LMSimone
JMohseni
PHayter
JEHodgkinson
KAChow
EWBassett
AS Linkage disequilibrium mapping of schizophrenia susceptibility to the
CAPON region of chromosome 1q22.
Am J Hum Genet 2004;74
(5)
1057- 1063
PubMedGoogle Scholar 57.Zheng
YLi
HQin
WChen
WDuan
YXiao
YLi
CZhang
JLi
XFeng
GHe
L Association of the carboxyl-terminal PDZ ligand of neuronal nitric oxide synthase gene with schizophrenia in the Chinese Han population.
Biochem Biophys Res Commun 2005;328
(4)
809- 815
PubMedGoogle Scholar 58.De Luca
VVoineskos
SWong
GKennedy
JL Genetic interaction between α4 and β2 subunits of high affinity nicotinic receptor: analysis in schizophrenia.
Exp Brain Res 2006;174
(2)
292- 296
PubMedGoogle Scholar 59.Ni
XValente
JAzevedo
MHPato
MTPato
CNKennedy
JL Connexin 50 gene on human chromosome 1q21 is associated with schizophrenia in matched case control and family-based studies.
J Med Genet 2007;44
(8)
532- 536
PubMedGoogle Scholar 60.Puri
V McQuillin
AChoudhury
KDatta
SPimm
JThirumalai
SKrasucki
RLawrence
JQuested
DBass
NMoorey
HMorgan
JPunukollu
BKandasami
GCurtis
DGurling
H Fine mapping by genetic association implicates the chromosome 1q23.3 gene UHMK1, encoding a serine/threonine protein kinase, as a novel schizophrenia susceptibility gene.
Biol Psychiatry 2007;61
(7)
873- 879
PubMedGoogle Scholar 61.Puri
V McQuillin
ADatta
SChoudhury
KPimm
JThirumalai
SKrasucki
RLawrence
JQuested
DBass
NCrombie
CFraser
GWalker
NMoorey
HRay
MKSule
ACurtis
DClair
DSGurling
H Confirmation of the genetic association between the U2AF homology motif (UHM) kinase 1 (UHMK1) gene and schizophrenia on chromosome 1q23.3.
Eur J Hum Genet 2008;16
(10)
1275- 1282
Google Scholar 62.Stefansson
HRujescu
DCichon
SPietilainen
OPIngason
ASteinberg
S Fossdal
RSigurdsson
ESigmundsson
TBuizer-Voskamp
JEHansen
TJakobsen
KDMuglia
PFrancks
CMatthews
PMGylfason
AHalldorsson
BVGudbjartsson
DThorgeirsson
TESigurdsson
AJonasdottir
AJonasdottir
ABjornsson
AMattiasdottir
SBlondal
THaraldsson
MMagnusdottir
BBGiegling
IMoller
HJHartmann
AShianna
KVGe
DNeed
ACCrombie
CFraser
GWalker
NLonnqvist
JSuvisaari
JTuulio-Henriksson
APaunio
TToulopoulou
TBramon
EDi Forti
MMurray
RRuggeri
MVassos
ETosato
SWalshe
MLi
TVasilescu
CMuhleisen
TWWang
AGUllum
HDjurovic
SMelle
IOlesen
JKiemeney
LAFranke
BSabatti
CFreimer
NBGulcher
JRThorsteinsdottir
UKong
AAndreassen
OAOphoff
RAGeorgi
ARietschel
MWerge
TPetursson
HGoldstein
DBNöthen
MMPeltonen
LCollier
DASt Clair
DStefansson
K Large recurrent microdeletions associated with schizophrenia.
Nature 2008;455
(7210)
232- 236
PubMedGoogle Scholar 63.International Schizophrenia Consortium, Rare chromosomal deletions and duplications increase risk of schizophrenia.
Nature 2008;455
(7210)
237- 241
PubMedGoogle Scholar 64.Kosower
NSGerad
LGoldstein
MParasol
NZipser
YRagolsky
MRozencwaig
SElkabetz
EAbramovitch
YLerer
BWeizman
A Constitutive heterochromatin of chromosome 1 and Duffy blood group alleles in schizophrenia.
Am J Med Genet 1995;60
(2)
133- 138
PubMedGoogle Scholar 65.Fañanás
LFuster
CGuillamat
RMiro
R Chromosomal fragile site 1q21 in schizophrenic patients.
Am J Psychiatry 1997;154
(5)
716
PubMedGoogle Scholar