Refining phenotype characterization in genetic linkage studies of schizophrenia

Abstract

A definitive replicable genetic linkage for a major locus underlying the susceptibility to schizophrenia has not been identified to date. Although there are several possible explanations for the failure to find linkage in schizophrenia, one major problem is that the range of phenotypic expressions of the genes for schizophrenia has not been clarified. A more refined understanding of the various phenotypic expressions of a gene related to schizophrenia would enhance the power of studies designed to detect a genetic linkage with a major chromosomal locus and would benefit other strategies for understanding the etiology of schizophrenia.

The genes for schizophrenia may be frequently expressed in relatives of schizophrenic patients, although with less severe symptoms than those of chronic schizophrenia. Two series of findings support this notion. Nonschizophrenic relatives of schizophrenic patients demonstrate an increased incidence of nonpsychotic schizophrenia‐like symptoms and traits, and they manifest deficit performances on several different laboratory tests of neurocognitive functioning. A more refined phenotypic expression of a schizophrenia‐related gene may thus be indicated by personality traits and subclinical neurocognitive deficits.

These personality traits and neurocognitive deficits are considered here as possible aids in the identification of affected cases in genetic linkage studies of schizophrenia. Terminology for different indicators of neurocognitive deficits is introduced, and the relative utility of personality traits and indicators of neurocognitive deficit for genetic linkage studies is discussed. As specific examples, schizophrenia‐related personality traits that are unrelated to affective symptoms and performance deficits on tasks of eye tracking and continuous attention are considered for strategies for broadening phenotype characterization without reducing the specificity of affected case identification.     

Refining phenotype characterization in genetic linkage studies of schizophrenia.

A definitive replicable genetic linkage for a major locus underlying the susceptibility to schizophrenia has not been identified to date. Although there are several possible explanations for the failure to find linkage in schizophrenia, one major problem is that the range of phenotypic expressions of the genes for schizophrenia has not been clarified. A more refined understanding of the various phenotypic expressions of a gene related to schizophrenia would enhance the power of studies designed to detect a genetic linkage with a major chromosomal locus and would benefit other strategies for understanding the etiology of schizophrenia. The genes for schizophrenia may be frequently expressed in relatives of schizophrenic patients, although with less severe symptoms than those of chronic schizophrenia. Two series of findings support this notion. Nonschizophrenic relatives of schizophrenic patients demonstrate an increased incidence of nonpsychotic schizophrenia-like symptoms and traits, and they manifest deficit performances on several different laboratory tests of neurocognitive functioning. A more refined phenotypic expression of a schizophrenia-related gene may thus be indicated by personality traits and subclinical neurocognitive deficits. These personality traits and neurocognitive deficits are considered here as possible aids in the identification of affected cases in genetic linkage studies of schizophrenia. Terminology for different indicators of neurocognitive deficits is introduced, and the relative utility of personality traits and indicators of neurocognitive deficit for genetic linkage studies is discussed. As specific examples, schizophrenia-related personality traits that are unrelated to affective symptoms and performance deficits on tasks of eye tracking and continuous attention are considered for strategies for broadening phenotype characterization without reducing the specificity of affected case identification.     

Programs,databases, and expert systems for human geneticists — a survey

We present an overview of the variety of databases and programs that offer substantial aid to medical and molecular geneticists. Databases and expert systems for genetic diseases and birth defects, programs for segregation and linkage analysis, certain DNA and protein sequence databases, and information resources in general for molecular biology are addressed. These systems cannot be used effectively without the newly developed techniques of information exchange based on international computer networks. A short introduction is given to the Internet and to European institutions and organizations that offer help with the aquisition and use of bioinformatic resources.     

A major lung cancer susceptibility locus maps to chromosome 6q23-25

Lung cancer is a major cause of death in the United States and other countries. The risk of lung cancer is greatly increased by cigarette smoking and by certain occupational exposures, but familial factors also clearly play a major role. To identify susceptibility genes for familial lung cancer, we conducted a genomewide linkage analysis of 52 extended pedigrees ascertained through probands with lung cancer who had several first-degree relatives with the same disease. Multipoint linkage analysis, under a simple autosomal dominant model, of all 52 families with three or more individuals affected by lung, throat, or laryngeal cancer, yielded a maximum heterogeneity LOD score (HLOD) of 2.79 at 155 cM on chromosome 6q (marker D6S2436). A subset of 38 pedigrees with four or more affected individuals yielded a multipoint HLOD of 3.47 at 155 cM. Analysis of a further subset of 23 multigenerational pedigrees with five or more affected individuals yielded a multipoint HLOD score of 4.26 at the same position. The 14 families with only three affected relatives yielded negative LOD scores in this region. A predivided samples test for heterogeneity comparing the LOD scores from the 23 multigenerational families with those from the remaining families was significant (P=.007). The 1-HLOD multipoint support interval from the multigenerational families extends from C6S1848 at 146 cM to 164 cM near D6S1035, overlapping a genomic region that is deleted in sporadic lung cancers as well as numerous other cancer types. Parametric linkage and variance-components analysis that incorporated effects of age and personal smoking also supported linkage in this region, but with somewhat diminished support. These results localize a major susceptibility locus influencing lung cancer risk to 6q23-25.     

Selection in complex genetic systems IV. Multiple alleles and interactions between two loci

Summary A symmetric viability model for two loci with two alleles at one locus and m alleles at the other is suggested and analyzed. The analysis of the equilibria is complete if the two loci are absolutely linked, while if recombination is allowed the analysis is incomplete. The dynamics of the mode! resemble those of the two locus two allele model, namely that for loose linkage there will be no correlation between the loci and for tight linkage there may be strong correlation. The major caveats to this are: 1. The equilibria stable for tight linkage may belong to an array of different structures dependent on the selection and the number of alleles. 2. If both loci are overdominant in viability, the stable equilibria always contain all alleles segregating in the population; otherwise, the stable equilibria may only be two locus two allele high complementarity equilibria for tight linkage. 3. For intermediate linkage values and special selection values the boundary two locus two allele high complementarity equilibria may be stable simultaneously with the totally polymorphic central point at which there is no association between the loci.Dedicated to the memory of Ove Frydenberg.Research supported in part by a grant from the Danish Natural Science Research Council, a grant from National Science Foundation, U.S.A., and by USPHS grant NIH 10452-09-11.     

Linkage analysis of "necessary" disease loci versus "susceptibility" loci.

The association of some diseases with specific alleles of certain genetic markers has been difficult to explain. Several explanations have been proposed for the phenomenon of association, e.g. the existence of multiple, interacting genes (epistasis) or a disease locus in linkage disequilibrium with the marker locus. One might suppose that when marker data from families with associated diseases are analyzed for linkage, the existence of the association would assure that linkage will be found, and found at a tight recombination fraction. In fact, however, linkage analyses of some diseases associated with HLA, as well as diseases associated with alleles at other loci located throughout the genome, show significant evidence against linkage, and others show loose linkage, to the puzzlement of many researchers. In part, the puzzlement arises because linkage analysis is ideal for looking for loci that are necessary, even if not sufficient, for disease expression but may be much less useful for finding loci that are neither necessary nor sufficient for disease expression (so-called susceptibility loci). This work explores what happens when one looks for linkage to susceptibility loci. A susceptibility locus in this case means that the allele increases risk but is neither necessary nor sufficient for disease expression. It might be either an allele at the marker locus itself that is increasing susceptibility or an allele at a locus in linkage disequilibrium with the marker. This work uses computer simulation to examine how linkage analyses behave when confronted with data from such a model.(ABSTRACT TRUNCATED AT 250 WORDS)     

Comparison off genetic distance and order off DNA markers in five populations of rice.

A group of about 300 evenly distributed DNA markers from a high density RFLP linkage map of rice constructed using an F2 population derived from a japonica variety, Nipponbare, and an indica variety, Kasalath, were used to evaluate gene order and genetic distance in four other rice mapping populations. The purpose of this study was to determine the degree to which information gained from the high density linkage map could be applied to other mapping populations, particularly with regard to its utility in bridging quantitative traits and molecular and physical mapping information. The mapping populations consisted of two F2 populations derived from Dao Ren Qiao/Fl-1084 and Kinandangputi/Fl-1007, recombinant inbred lines from Asominori/IR24, and a backcross population from Sasanishiki/Habataki//Sasanishiki. All DNA markers commonly mapped in the four populations showed the same linkage groups as in the Nipponbare/Kasalath linkage map with conserved linkage order. The genetic distance between markers among the different populations did not vary to a significant level in any of the 12 chromosomes. The differences in some markers could be attributed to the size of the population used in the construction of the linkage maps. Furthermore, the conservation of linkage order found in the distal region of chromosomes 11 and 12 was also confirmed in the RFLP maps based on the four populations of rice. These suggest that any major genetic information from the Nipponbare/Kasalath map can be expected to be approximately the same in other crosses or populations. This high density RFLP linkage map, which is being utilized in constructing a physical map of rice, can be very useful in interpreting genome structure with great accuracy in other populations. Key words : linkage map, japonica, indica, gene order, genetic distance.     

Some mass-spectral and n.m.r. analytical studies of a glutathione conjugate of aflatoxin B1.

A system for the formation of an aflatoxin B1-reduced glutathione conjugate in vitro was developed, capable of yielding 80% conversion of aflatoxin B1 into the conjugate. A reverse-phase high-pressure-liquid-chromatography system was also devised that not only facilitates improved resolution of the compound but that, by manipulation of the pH, is also capable of an extensive purification of the compound from other aflatoxin B1 metabolites in a single step. Material produced by these techniques, after further purification, has been used in 1H-n.m.r. and mass-spectroscopic studies. Results were obtained that support the proposed linkage of the aflatoxin B1 to reduced glutathione in a 1:1 molar ratio via a thioether linkage. Amino acid analyses were also consistent with this structure. The absence of a Schiff-base linkage of aflatoxin B1 8,9-dihydrodiol to glutamate was further demonstrated by the presence of a gamma-glutamyltransferase-catalysed-transferable glutamate moiety. These data are consistent with the structure 8,9-dihydro-8-(S-glutathionyl)-9-hydroxy-aflatoxin B1.     

A symmetric two locus model with viability and fertility selection

A two locus deterministic population genetic model is analysed. One locus is under viability selection, the other under fertility selection with both forms of selection completely symmetric. It is shown that linkage equilibrium may occur at two different equilibrium points. For a two-locus polymorphism to be stable, it is necessary that the viability locus be overdominant but not necessary that the fertility locus, considered separately, be able to support a stable polymorphism. The overlaps in stability are not as complex as under two locus symmetric fertilities, but considerably more complex than with symmetric viabilities. Extensions of the analysis for the central linkage equilibrium point with multiple viability and fertility loci are indicated.Research supported in part by NIH grants GM 28106 and GM 10452     

Gait reaction reconstruction and a heel strike algorithm

A mathematical model of gait ground loading is presented. The model allows the ground reactions produced by any particular single- or multiple-footfall pattern to be constructed, given a sufficient variety of other measured ground reactions. An algorithm which uses center of vertical pressure data only to determine the instants of successive heel strikes on a large force plate is then presented. Experiments show the high accuracy of the heel strike algorithm and show that reconstructions of the vertical component of ground reactions are typically within 3% of corresponding measured reactions. The techniques presented allow certain problems associated with small force plates and other problems associated with large force plates to be largely overcome.     

Genome scan meta-analysis of schizophrenia and bipolar disorder,part III: Bipolar disorder

Genome scans of bipolar disorder (BPD) have not produced consistent evidence for linkage. The rank-based genome scan meta-analysis (GSMA) method was applied to 18 BPD genome scan data sets in an effort to identify regions with significant support for linkage in the combined data. The two primary analyses considered available linkage data for “very narrow” (i.e., BP-I and schizoaffective disorder–BP) and “narrow” (i.e., adding BP-II disorder) disease models, with the ranks weighted for sample size. A “broad” model (i.e., adding recurrent major depression) and unweighted analyses were also performed. No region achieved genomewide statistical significance by several simulation-based criteria. The most significant P values (<.01) were observed on chromosomes 9p22.3-21.1 (very narrow), 10q11.21-22.1 (very narrow), and 14q24.1-32.12 (narrow). Nominally significant P values were observed in adjacent bins on chromosomes 9p and 18p-q, across all three disease models on chromosomes 14q and 18p-q, and across two models on chromosome 8q. Relatively few BPD pedigrees have been studied under narrow disease models relative to the schizophrenia GSMA data set, which produced more significant results. There was no overlap of the highest-ranked regions for the two disorders. The present results for the very narrow model are promising but suggest that more and larger data sets are needed. Alternatively, linkage might be detected in certain populations or subsets of pedigrees. The narrow and broad data sets had considerable power, according to simulation studies, but did not produce more highly significant evidence for linkage. We note that meta-analysis can sometimes provide support for linkage but cannot disprove linkage in any candidate region.     

Osteopenia in 37 members of seven families: analysis based on a model of dominant inheritance.

BACKGROUND: The genetic factors involved in determining bone mineral density (BMD) have not been fully elucidated. We have begun genetic linkage analysis of seven families in which many members are osteopenic, in order to identify chromosomal loci that are potentially involved in determining BMD. MATERIALS AND METHODS: Spine BMD was measured in 143 members of seven kindred with familial osteopenia. The absolute BMD values for the spine (L2-L4) were converted to the age-, gender-, and weight-adjusted Z scores, and this corrected value was used as the quantitative trait on which to base subsequent genetic analyses. Simulations of linkage were performed in order to determine the information content of the pedigree set, and actual linkage analysis was conducted using polymorphic markers either within or near three candidate loci: COL1A1, COL1A2, and vitamin D receptor (VDR). RESULTS: The distribution of the corrected Z scores was bimodal (p = 0.001) suggesting a monogenic mode of inheritance of the low BMD trait. Simulation of linkage analysis suggested that the family data set was sufficient to detect linkage under a single major gene model. Actual linkage analysis did not support linkage to the three candidate loci. In addition, the VDR genotype was not statistically associated with low bone density at the spine. CONCLUSIONS: Loci other than COL1A1, COL1A2 and VDR are very likely responsible for the low BMD trait observed in these families. These families are suitable for a genome-wide screen using microsatellite repeats in order to identify the loci that are involved in osteopenia.     

A framework for structural equation models in general pedigrees

Background/Aims: Structural Equation Modeling (SEM) is an analysis approach that accounts for both the causal relationships between variables and the errors associated with the measurement of these variables. In this paper, a framework for implementing structural equation models (SEMs) in family data is proposed. Methods: This framework includes both a latent measurement model and a structural model with covariates. It allows for a wide variety of models, including latent growth curve models. Environmental, polygenic and other genetic variance components can be included in the SEM. Kronecker notation makes it easy to separate the SEM process from a familial correlation model. A limited information method of model fitting is discussed. We show how missing data and ascertainment may be handled. We give several examples of how the framework may be used. Results: A simulation study shows that our method is computationally feasible, and has good statistical properties. Conclusion: Our framework may be used to build and compare causal models using family data without any genetic marker data. It also allows for a nearly endless array of genetic association and/or linkage tests. A preliminary Matlab program is available, and we are currently implementing a more complete and user-friendly R package.     

Fine-mapping QTLs in advanced intercross lines and other outbred populations

Quantitative genetic studies in model organisms, particularly in mice, have been extremely successful in identifying chromosomal regions that are associated with a wide variety of behavioral and other traits. However, it is now widely understood that identification of the underlying genes will be far more challenging. In the last few years, a variety of populations have been utilized in an effort to more finely map these chromosomal regions with the goal of identifying specific genes. The common property of these newer populations is that linkage disequilibrium spans relatively short distances, which permits fine-scale mapping resolution. This review focuses on advanced intercross lines (AILs) which are the simplest such population. As originally proposed in 1995 by Darvasi and Soller, an AIL is the product of intercrossing two inbred strains beyond the F₂ generation. Unlike recombinant inbred strains, AILs are maintained as outbred populations; brother–sister matings are specifically avoided. Each generation of intercrossing beyond the F₂ further degrades linkage disequilibrium between adjacent makers, which allows for fine-scale mapping of quantitative trait loci (QTLs). Advances in genotyping technology and techniques for the statistical analysis of AILs have permitted rapid advances in the application of AILs. We review some of the analytical issues and available software, including QTLRel, EMMA, EMMAX, GEMMA, TASSEL, GRAMMAR, WOMBAT, Mendel, and others.     

Ethnic variation in genetic disease: possible roles of hitchhiking and epistasis.

The high incidence of some genetic diseases in certain ethnic groups is important in planning of medical genetic programs. Simple interaction models predict that at least some lethal recessive alleles will have "hitchhiked" to increased frequencies because of linkage to genes whose alleles have been favored by selection for other reasons in certain populations. In the absence of linkage or epistasis with a gene favored by selection, heterozygote advantage for a recessive lethal may produce the same phenomenon. In the hitchhiking model (linkage), the increase in the gene frequency is temporary, but the length of time that the increased gene frequency is at least double the base frequency may be quite long. Changes in gene frequency for the unlinked epistatic model result in a new equilibrium with a possibly higher gene frequency. The most likely chromosomal regions in which hitchhiked lethal recessives would be found are in the vicinity of genes whose allelic frequencies vary substantially among human racial groups (e.g., Gm, Rh, Duffy, lactose tolerance, or HL-A). There will be a hitchhiking effect if recombination distance is less than the selective advantage. The closer the linkage of two loci, the easier hitchhiking effects will be to detect. Hitchhiking is suggested by nonrandom association of the recessive disease and one of the selected markers, as in the case of Gm and cystic fibrosis. However, there is so far insufficient evidence of linkage between them. More pedigree information is necessary than is now available.     

Molecular genetics and heterogeneity in manic depression

Recent research has shown that there are X-linked and possibly chromosome 11-linked forms of manic depression as well as at least one other autosomal form. Segregation analyses of large affected families and the finding of genetic linkage between chromosome specific markers and manic depression mutations provide strong evidence that bipolar as well as unipolar forms of manic depression (MD) within the same family are inherited as a dominant gene disorder. This clarification of the etiology of certain types of depression should bring changed attitudes within psychiatry and may serve to stimulate discussion of the role of evolutionary mechanisms. From a clinical point of view, it has now become possible to determine whether clinical (phenotypic) variation reflects the underlying genotypic heterogeneity of linkage. A preliminary analysis of data from four recent studies shows that there is no clear correlation between such clinical features as the ratio of unipolar to bipolar cases and the genotypic form of manic depression. Further recombinant DNA research, proven to be successful in other genetic diseases, can soon be applied to manic depression. The specific problems posed by manic depression for these techniques are discussed.     

Interventional neuroradiology.

A wide variety of diseases affecting the central nervous system and head and neck can be treated using interventional neuroradiologic techniques. These new treatments have depended on advances in radiologic imaging, catheter technology, and the development of new embolic agents. These procedures may be an adjunct to other therapy, palliative or curative. Diseases for which interventional neuroradiologic techniques have been major advances in treatment include cerebral aneurysms, vasospasm after subarachnoid hemorrhage, cerebral arteriovenous malformations, dural arteriovenous fistulas, dural sinus thrombosis, atherosclerosis, scalp arteriovenous fistulas, carotid-cavernous fistulas, and stroke. This field is rapidly evolving as advances are made in catheter technology and new embolic agents are developed.     

Syntactical Organization of the Gargle Vocalization of the Black-capped Chickadee,Parus atricapillus

The gargle is a complex call of the black-capped chickadee Parus atricapillus associated with agonistic interactions. A large sample of calls was recorded from five sites in southeastern Wisconsin. Three sites were near one another and the other two with a different dialect pattern were more distant (all sites were within 15 km of each other). A typical call consists of about 6 notes selected from an array of about 15 different note-types that are shared among individuals in the same population. Note-types are given in many combinations, and any individual may deliver many different variations. Analysis of two- and three-note sequences shows a syntactical organization with certain note-types much more likely to precede or follow than others. In some cases linkage between note-types was nearly complete. Some differences occurred among individuals and sites, but general organizational rules were similar. One effect on ordering was acoustical structure of the notes. Higher pitched, shorter duration notes were more likely at the beginning than the end of the call which often consisted of trills. Organization of the gargle was compared to that of the chick-a-dee call of the same species. Similarities include recombinant properties, general acoustic rules and syntax. A major difference may be that in chick-a-dee calls various note combinations constituting a call have different meanings, while the various note combinations comprising an individual's repertoire of gargles may be functionally equivalent.     

Mapping loci influencing blood pressure in the Framingham pedigrees using model-free LOD score analysis of a quantitative trait

This paper presents a method of performing model-free LOD-score based linkage analysis on quantitative traits. It is implemented in the QMFLINK program. The method is used to perform a genome screen on the Framingham Heart Study data. A number of markers that show some support for linkage in our study coincide substantially with those implicated in other linkage studies of hypertension. Although the new method needs further testing on additional real and simulated data sets we can already say that it is straightforward to apply and may offer a useful complementary approach to previously available methods for the linkage analysis of quantitative traits.     

Simultaneous maximum likelihood estimation of linkage and linkage phases in outcrossing species

With the advent of new molecular marker technologies, it is now feasible to initiate genome projects for outcrossing plant species, which have not received much attention in genetic research, despite their great agricultural and environmental value. Because outcrossing species typically have heterogeneous genomes, data structure for molecular markers representing an entire genome is complex: some markers may have more alleles than others, some markers are codominant whereas others are dominant, and some markers are heterozygous in one parent but fixed in the other parent whereas the opposite can be true for other markers. A major difficulty in analyzing these different types of marker at the same time arises from uncertainty about parental linkage phases over markers. In this paper, we present a general maximum-likelihood-based algorithm for simultaneously estimating linkage and linkage phases for a mixed set of different marker types containing fully informative markers (segregating 1:1:1:1) and partially informative markers (or missing markers, segregating 1:2:1, 3:1, and 1:1) in a full-sib family derived from two outbred parent plants. The characterization of linkage phases is based on the posterior probability distribution of the assignment of alternative alleles at given markers to two homologous chromosomes of each parent, conditional on the observed phenotypes of the markers. Two- and multi-point analyses are performed to estimate the recombination fraction and determine the most likely linkage phase between different types of markers. A numerical example is presented to demonstrate the statistical properties of the model for characterizing the linkage phase between markers.