Evidence for a new Graves disease susceptibility locus at chromosome 18q21

Graves disease (GD) is a common autoimmune thyroid disorder that is inherited as a complex multigenic trait. By using a single microsatellite marker at each locus, we screened the type 1 diabetes loci IDDM4, IDDM5, IDDM6, IDDM8, and IDDM10 and the fucosyltransferase-2 locus for linkage in sib pairs with GD. This showed a two-point nonparametric linkage (NPL) score of 1.57 (P=.06) at the IDDM6 marker D18S41, but NPL scores were <1.0 at the other five loci. Thus, the investigation of the IDDM6 locus was extended by genotyping 11 microsatellite markers spanning 48 cM across chromosome 18q12-q22 in 81 sib pairs affected with autoimmune thyroid disease (AITD). Multipoint analysis, designating all AITD sib pairs as affected, showed a peak NPL score of 3.46 (P=.0003), at the marker D18S487. Designation of only GD cases as affected (74 sib pairs) showed a peak NPL score of 3.09 (P=.001). Linkage to this region has been demonstrated in type 1 diabetes (IDDM6), rheumatoid arthritis, and systemic lupus erythematosus, which suggests that this locus may have a role in several forms of autoimmunity.     

A unified sampling approach for multipoint analysis of qualitative and quantitative traits in sib pairs

Recent advances in molecular biology have enhanced the opportunity to conduct multipoint mapping for complex diseases. Concurrently, one sees a growing interest in the use of quantitative traits in linkage studies. Here, we present a multipoint sib-pair approach to locate the map position (tau) of a trait locus that controls the observed phenotype (qualitative or quantitative), along with a measure of statistical uncertainty. This method builds on a parametric representation for the expected identical-by-descent statistic at an arbitrary locus, conditional on an event reflecting the sampling scheme, such as affected sib pairs, for qualitative traits, or extreme discordant (ED) sib pairs, for quantitative traits. Our results suggest that the variance about tau&d4;, the estimator of tau, can be reduced by as much as 60%-70% by reducing the length of intervals between markers by one half. For quantitative traits, we examine the precision gain (measured by the variance reduction in tau&d4;) by genotyping extremely concordant (EC) sib pairs and including them along with ED sib pairs in the statistical analysis. The precision gain depends heavily on the residual correlation of the quantitative trait for sib pairs but considerably less on the allele frequency and exact genetic mechanism. Since complex traits involve multiple loci and, hence, the residual correlation cannot be ignored, our finding strongly suggests that one should incorporate EC sib pairs along with ED sib pairs, in both design and analysis. Finally, we empirically establish a simple linear relationship between the magnitude of precision gain and the ratio of the number of ED pairs to the number of EC pairs. This relationship allows investigators to address issues of cost effectiveness that are due to the need for phenotyping and genotyping subjects.     

Evidence for a Language Quantitative Trait Locus on Chromosome 7q in Multiplex Autism Families

Autism is a syndrome characterized by deficits in language and social skills and by repetitive behaviors. We hypothesized that potential quantitative trait loci (QTLs) related to component autism endophenotypes might underlie putative or significant regions of autism linkage. We performed nonparametric multipoint linkage analyses, in 152 families from the Autism Genetic Resource Exchange, focusing on three traits derived from the Autism Diagnostic Interview: "age at first word," "age at first phrase," and a composite measure of "repetitive and stereotyped behavior." Families were genotyped for 335 markers, and multipoint sib pair linkage analyses were conducted. Using nonparametric multipoint linkage analysis, we found the strongest QTL evidence for age at first word on chromosome 7q (nonparametric test statistic [Z] 2.98; P=.001), and subsequent linkage analyses of additional markers and association analyses in the same region supported the initial result (Z=2.85, P=.002; chi(2)=18.84, df 8, P=.016). Moreover, the peak fine-mapping result for repetitive behavior (Z=2.48; P=.007) localized to a region overlapping this language QTL. The putative autism-susceptibility locus on chromosome 7 may be the result of separate QTLs for the language and repetitive or stereotyped behavior deficits that are associated with the disorder.     

Autosomal genomic scan for loci linked to obesity and energy metabolism in Pima Indians.

An autosomal genomic scan to search for linkage to obesity and energy metabolism was completed in Pima Indians, a population prone to obesity. Obesity was assessed by percent body fat (by hydrodensitometry) and fat distribution (the ratio of waist circumference to thigh circumference). Energy metabolism was measured in a respiratory chamber as 24-h metabolic rate, sleeping metabolic rate, and 24-h respiratory quotient (24RQ), an indicator of the ratio of carbohydrate oxidation to fat oxidation. Five hundred sixteen microsatellite markers with a median spacing of 6.4 cM were analyzed, in 362 siblings who had measurements of body composition and in 220 siblings who had measurements of energy metabolism. These comprised 451 sib pairs in 127 nuclear families, for linkage analysis to obesity, and 236 sib pairs in 82 nuclear families, for linkage analysis to energy metabolism. Pointwise and multipoint methods for regression of sib-pair differences in identity by descent, as well as a sibling-based variance-components method, were used to detect linkage. LOD scores >=2 were found at 11q21-q22, for percent body fat (LOD=2.1; P=.001), at 11q23-q24, for 24-h energy expenditure (LOD=2.0; P=.001), and at 1p31-p21 (LOD=2.0) and 20q11.2 (LOD=3.0; P=.0001), for 24RQ, by pointwise and multipoint analyses. With the variance-components method, the highest LOD score (LOD=2.3 P=.0006) was found at 18q21, for percent body fat, and at 1p31-p21 (LOD=2.8; P=.0003), for 24RQ. Possible candidate genes include LEPR (leptin receptor), at 1p31, and ASIP (agouti-signaling protein), at 20q11.2.     

Evidence for Linkage of Bipolar Disorder to Chromosome 18 with a Parent-of-Origin Effect

A susceptibility gene on chromosome 18 and a parent-of-origin effect have been suggested for bipolar affective disorder (BPAD). We have studied 28 nuclear families selected for apparent unilineal transmission of the BPAD phenotype, by using 31 polymorphic markers spanning chromosome 18. Evidence for linkage was tested with affected-sib-pair and LOD score methods under two definitions of the affected phenotype. The affected-sib-pair analyses indicated excess allele sharing for markers on 18p within the region reported previously. The greatest sharing was at D18S37: 64% in bipolar and recurrent unipolar (RUP) sib pairs (P = .0006). In addition, excess sharing of the paternally, but not maternally, transmitted alleles was observed at three markers on 18q: at D18S41, 51 bipolar and RUP sib pairs were concordant for paternally transmitted alleles, and 21 pairs were discordant (P = .0004). The evidence for linkage to loci on both 18p and 18q was strongest in the 11 paternal pedigrees, i.e., those in which the father or one of the father's sibs is affected. In these pedigrees, the greatest allele sharing (81%; P = .00002) and the highest LOD score (3.51; θ = 0.0) were observed at D18S41. Our results provide further support for linkage of BPAD to chromosome 18 and the first molecular evidence for a parent-of-origin effect operating in this disorder. The number of loci involved, and their precise location, require further study.     

Equivalence between Haseman-Elston and variance-components linkage analyses for sib pairs

The Haseman-Elston regression method offers a simpler alternative to variance-components (VC) models, for the linkage analysis of quantitative traits. However, even the "revisited" method, which uses the cross-product--rather than the squared difference--in sib trait values, is, in general, less powerful than VC models. In this report, we clarify the relative efficiencies of existing Haseman-Elston methods and show how a new Haseman-Elston method can be constructed to have power equivalent to that of VC models. This method uses as the dependent variable a linear combination of squared sums and squared differences, in which the weights are determined by the overall trait correlation between sibs in a population. We show how this method can be used for both the selection of maximally informative sib pairs for genotyping and the subsequent analysis of such selected samples.     

Effects of misspecification of allele frequencies on the type I error rate of model-free linkage analysis

Linkage analyses of simulated quantitative trait data were performed using the Haseman-Elston (H-E) sib pair regression test to investigate the effects of inaccurate allele frequency estimates on the type I error rates of this test. Computer simulations generating a quantitative trait in nuclear families were performed using GASP [1]. Assuming no linkage, several data sets were simulated; they differed in marker allele numbers and frequencies, number of sib pairs and number of sibships. Each set of simulated data was analyzed using (1) all parental marker data, (2) half of the parental marker data, and (3) no parental marker data, using both correct and incorrect allele frequencies in the latter 2 cases. The H-E sib pair linkage method was found to be robust to misspecification of marker allele frequencies regardless of the number of alleles.     

A quantitative-trait locus on chromosome 6p influences different aspects of developmental dyslexia.

Recent application of nonparametric-linkage analysis to reading disability has implicated a putative quantitative-trait locus (QTL) on the short arm of chromosome 6. In the present study, we use QTL methods to evaluate linkage to the 6p25-21.3 region in a sample of 181 sib pairs from 82 nuclear families that were selected on the basis of a dyslexic proband. We have assessed linkage directly for several quantitative measures that should correlate with different components of the phenotype, rather than using a single composite measure or employing categorical definitions of subtypes. Our measures include the traditional IQ/reading discrepancy score, as well as tests of word recognition, irregular-word reading, and nonword reading. Pointwise analysis by means of sib-pair trait differences suggests the presence, in 6p21.3, of a QTL influencing multiple components of dyslexia, in particular the reading of irregular words (P=.0016) and nonwords (P=.0024). A complementary statistical approach involving estimation of variance components supports these findings (irregular words, P=.007; nonwords, P=.0004). Multipoint analyses place the QTL within the D6S422-D6S291 interval, with a peak around markers D6S276 and D6S105 consistently identified by approaches based on trait differences (irregular words, P=.00035; nonwords, P=.0035) and variance components (irregular words, P=.007; nonwords, P=.0038). Our findings indicate that the QTL affects both phonological and orthographic skills and is not specific to phoneme awareness, as has been previously suggested. Further studies will be necessary to obtain a more precise localization of this QTL, which may lead to the isolation of one of the genes involved in developmental dyslexia.     

Pleiotropic effects of a chromosome 3 locus on speech-sound disorder and reading

Speech-sound disorder (SSD) is a complex behavioral disorder characterized by speech-sound production errors associated with deficits in articulation, phonological processes, and cognitive linguistic processes. SSD is prevalent in childhood and is comorbid with disorders of language, spelling, and reading disability, or dyslexia. Previous research suggests that developmental problems in domains associated with speech and language acquisition place a child at risk for dyslexia. Recent genetic studies have identified several candidate regions for dyslexia, including one on chromosome 3 segregating in a large Finnish pedigree. To explore common genetic influences on SSD and reading, we examined linkage for several quantitative traits to markers in the pericentrometric region of chromosome 3 in 77 families ascertained through a child with SSD. The quantitative scores measured several processes underlying speech-sound production, including phonological memory, phonological representation, articulation, receptive and expressive vocabulary, and reading decoding and comprehension skills. Model-free linkage analysis was followed by identification of sib pairs with linkage and construction of core shared haplotypes. In our multipoint analyses, measures of phonological memory demonstrated the strongest linkage (marker D3S2465, P=5.6 x 10(-5), and marker D3S3716, P=6.8 x 10(-4)). Tests for single-word decoding also demonstrated linkage (real word reading: marker D3S2465, P=.004; nonsense word reading: marker D3S1595, P=.005). The minimum shared haplotype in sib pairs with similar trait values spans 4.9 cM and is bounded by markers D3S3049 and D3S3045. Our results suggest that domains common to SSD and dyslexia are pleiotropically influenced by a putative quantitative trait locus on chromosome 3.     

A Score-Statistic Approach for the Mapping of Quantitative-Trait Loci with Sibships of Arbitrary Size

The Haseman-Elston method is widely used for the mapping of quantitative-trait loci. However, this method does not use all the information in the data, because it only considers the sib-pair trait-value difference. In addition, the Haseman-Elston method was developed for independent sib pairs; its generalization to nonindependent sib pairs is not straightforward. Here we introduce a score test statistic derived from a normal likelihood based on multiplex sibship data, conditional on identical-by-descent sharing statuses. This score test is asymptotically equivalent to the corresponding likelihood-ratio test, but it is much easier to implement. Because the proposed test uses all of the trait values, it makes more efficient use of the data than does the Haseman-Elston method. The proposed test is naturally applicable to sibships of arbitrary size. The finite-sample properties of the proposed score statistic are evaluated via simulations.     

A whole-genome scan in 164 Dutch sib pairs with attention-deficit/hyperactivity disorder: suggestive evidence for linkage on chromosomes 7p and 15q

A genome scan was performed on 164 Dutch affected sib pairs (ASPs) with attention-deficit/hyperactivity disorder (ADHD). All subjects were white and of Dutch descent and were phenotyped according to criteria set out in the Diagnostic and Statistical Manual Of Mental Disorders, 4th edition. Initially, a narrow phenotype was defined, in which all the sib pairs met the full ADHD criteria (117 ASPs). In a broad phenotype, additional sib pairs were included, in which one child had an autistic-spectrum disorder but also met the full ADHD criteria (164 ASPs). A set of 402 polymorphic microsatellite markers with an average intermarker distance of 10 cM was genotyped and analyzed using the Mapmaker/sibs program. Regions with multipoint maximum likelihood scores (MLSs) >1.5 in both phenotypes were fine mapped with additional markers. This genome scan indicated several regions of interest, two of which showed suggestive evidence for linkage. The most promising chromosome region was located at 15q, with an MLS of 3.54 under the broad phenotype definition. This region was previously implicated in reading disability and autism. In addition, MLSs of 3.04 and 2.05 were found for chromosome regions 7p and 9q in the narrow phenotype. Except for a region on chromosome 5, no overlap was found with regions mentioned in the only other independent genome scan in ADHD reported to date.     

A genome-wide search for genes predisposing to familial psoriasis by using a stratification approach

We have performed a genome scan, using markers spaced by 10 cM, in the search for psoriasis-susceptibility loci. The family material of 134 affected sibling pairs was ascertained on the basis of a population genetic study in which 65% of the probands had two healthy parents. Genotyping results were analyzed for non-random excessive allele-sharing between sib pairs by using GENEHUNTER ver 1.1. A stratification approach was applied to increase the homogeneity of the material by means of an operational definition of joint complaints among affected individuals. Significant linkage to the human leukocyte antigen region on chromosome 6p in a cohort including 42 families without joint complaints (nonparametric linkage score of 2.83, P=0.002) strongly supported the validity of this operational definition as it replicated results from an earlier linkage report with similar stratification criteria. New candidate regions on chromosomes 3 and 15 were identified. The highest non-parametric linkage values in this study, 2.96 (P=0.0017) and 2.89 (P=0.0020), were reached on chromosome 15 in a subgroup with joint complaints and on chromosome 3 in a subgroup without joint complaints. In addition, confirmation of previously reported loci was established on chromosomes 4q, 6p, and 17q. This study indicates that distinct disease loci might be involved in psoriasis etiology for various phenotypes.     

Efficiency robust tests for mapping quantitative trait loci using extremely discordant sib pairs

In 1972, Haseman and Elston proposed a pioneering regression method for mapping quantitative trait loci using randomly selected sib pairs. Recently, the statistical power of their method was shown to be increased when extremely discordant sib pairs are ascertained. While the precise genetic model may not be known, prior information that constrains IBD probabilities is often available. We investigate properties of tests that are robust against model uncertainty and show that the power gain from further constraining IBD probabilities is marginal. The additional linkage information contained in the trait values can be incorporated by combining the Haseman-Elston regression method and a robust allele sharing test.     

Association of the chromosomal region 2q35 with type 1 diabetes mellitus in the Russian patients from Moscow

To map human chromosome 2 region associated with type 1 diabetes mellitus, 89 families with concordant and discordant sib pairs were analyzed. Linkage and association with type 1 diabetes were examined using polymorphic microsatellite markers spanning the region of about 4 Mb. The linkage plot was constructed, and association of the five microsatellite markers within the chromosomal region 2q35 was examined. Polymorphic marker D2S137 (Z' = 3.225, p(c) = 0.0048) demonstrated maximum linkage and association with type 1 diabetes.     

A susceptibility locus for bipolar affective disorder on chromosome 4q35.

Bipolar affective disorder (BAD) affects approximately 1% of the population and shows strong heritability. To identify potential BAD susceptibility loci, we undertook a 15-cM genome screen, using 214 microsatellite markers on the 35 most informative individuals of a large, statistically powerful pedigree. Data were analyzed by parametric two-point linkage methods under several diagnostic models. LOD scores >1.00 were obtained for 21 markers, with four of these >2.00 for at least one model. The remaining 52 individuals in the family were genotyped with these four markers, and LOD scores remained positive for three markers. A more intensive screen was undertaken in these regions, with the most positive results being obtained for chromosome 4q35. Using a dominant model of inheritance with 90% maximum age-specific penetrance and including bipolar I, II, schizoaffective/mania, and unipolar individuals as affected, we obtained a maximum two-point LOD score of 2.20 (theta = .15) at D4S1652 and a maximum three-point LOD score of 3.19 between D4S408 and D4S2924. Nonparametric analyses further supported the presence of a locus on chromosome 4q35. A maximum score of 2.62 (P=.01) was obtained between D4S1652 and D4S171 by use of the GENEHUNTER program, and a maximum score of 3.57 (P=.0002) was obtained at D4S2924 using the affected pedigree member method. Analysis of a further 10 pedigrees suggests the presence of this locus in at least one additional family, indicating a possible predisposing locus and not a pedigree-specific mutation. Our results suggest the presence of a novel BAD susceptibility locus on chromosome 4q35.     

A new gene for Tourette's syndrome: a window into causal mechanisms?

Gilles de la Tourette syndrome (GTS) is a neurodevelopmental disorder characterized by impairing motor-vocal tics. Locating genetic loci by associating the phenotype with DNA translocations, inversions, gain or losses, State et al. identified SLITRK1 as a candidate gene in an individual with GTS and inv(13) (q31.1; q33.1). This gene was also associated with abnormal axonal-dendritic development in embryonic mouse cells. Although SLITRK1 is not a major causal gene for GTS, it can shed light on our understanding of the gene-based neural correlates of this disease.     

Model-free linkage analysis with covariates confirms linkage of prostate cancer to chromosomes 1 and 4

As with many complex genetic diseases, genome scans for prostate cancer have given conflicting results, often failing to provide replication of previous findings. One factor contributing to the lack of consistency across studies is locus heterogeneity, which can weaken or even eliminate evidence for linkage that is present only in a subset of families. Currently, most analyses either fail to account for locus heterogeneity or attempt to account for it only by partitioning data sets into smaller and smaller portions. In the present study, we model locus heterogeneity among affected sib pairs with prostate cancer by including covariates in the linkage analysis that serve as surrogate measures of between-family linkage differences. The model is a modification of the Olson conditional logistic model for affected relative pairs. By including Gleason score, age at onset, male-to-male transmission, and/or number of affected first-degree family members as covariates, we detected linkage near three locations that were previously identified by linkage (1q24-25 [HPC1; LOD score 3.25, P=.00012], 1q42.2-43 [PCAP; LOD score 2.84, P=.0030], and 4q [LOD score 2.80, P=.00038]), near the androgen-receptor locus on Xq12-13 (AR; LOD score 3.06, P=.00053), and at five new locations (LOD score > 2.5). Without covariates, only a few weak-to-moderate linkage signals were found, none of which replicate findings of previous genome scans. We conclude that covariate-based linkage analysis greatly improves the likelihood that linked regions will be found by incorporation of information about heterogeneity within the sample.     

Linkage analysis of extremely discordant and concordant sibling pairs identifies quantitative trait loci influencing variation in human menopausal age

Age at natural menopause may be used as parameter for evaluating the rate of ovarian aging. Environmental factors determine only a small part of the large variation in menopausal age. Studies have shown that genetic factors are likely to be involved in variation in menopausal age. To identify quantitative-trait loci for this trait, we performed a genomewide linkage study with age at natural menopause as a continuous quantitative phenotype in Dutch sister pairs, through use of a selective sampling scheme. A total of 165 families were ascertained using extreme selected sampling and were genotyped for 417 markers. Data were analyzed by Haseman-Elston regression and by an adjusted variance-components analysis. Subgroup analyses for early and late menopausal age were conducted by Haseman-Elston regression. In the adjusted variance-components analysis, 12 chromosomes had a LOD score of > or =1.0. Two chromosomal regions showed suggestive linkage: 9q21.3 (LOD score 2.6) and Xp21.3 (LOD score 3.1). Haseman-Elston regression showed rather similar locations of the peaks but yielded lower LOD scores. A permutation test to obtain empirical P values resulted in a significant peak on the X chromosome. To our knowledge, this is the first study to attempt to identify loci responsible for variability in menopausal age and in which several chromosomal regions were identified with suggestive and significant linkage. Although the finding of the region on the X chromosome comes as no surprise, because of its widespread involvement in premature ovarian failure, the definition of which particular gene is involved is of great interest. The region on chromosome 9 deserves further consideration. Both findings require independent confirmation.     

Recursive partitioning models for linkage in COGA data

We have developed a recursive-partitioning (RP) algorithm for identifying phenotype and covariate groupings that interact with the evidence for linkage. This data-mining approach for detecting gene x environment interactions uses genotype and covariate data on affected relative pairs to find evidence for linkage heterogeneity across covariate-defined subgroups. We adapted a likelihood-ratio based test of linkage parameterized with relative risks to a recursive partitioning framework, including a cross-validation based deviance measurement for choosing optimal tree size and a bootstrap sampling procedure for choosing robust tree structure. ALDX2 category 5 individuals were considered affected, categories 1 and 3 unaffected, and all others unknown. We sampled non-overlapping affected relative pairs from each family; therefore, we used 144 affected pairs in the RP model. Twenty pair-level covariates were defined from smoking status, maximum drinks, ethnicity, sex, and age at onset. Using the all-pairs score in GENEHUNTER, the nonparametric linkage tests showed no regions with suggestive linkage evidence. However, using the RP model, several suggestive regions were found on chromosomes 2, 4, 6, 14, and 20, with detection of associated covariates such as sex and age at onset.     

Indications of linkage and association of Gilles de la Tourette syndrome in two independent family samples: 17q25 is a putative susceptibility region

Gilles de la Tourette syndrome (GTS) is characterized by multiple motor and phonic tics and high comorbidity rates with other neurobehavioral disorders. It is hypothesized that frontal-subcortical pathways and a complex genetic background are involved in the etiopathogenesis of the disorder. The genetic basis of GTS remains elusive. However, several genomic regions have been implicated. Among them, 17q25 appears to be of special interest, as suggested by various independent investigators. In the present study, we explored the possibility that 17q25 contributes to the genetic component of GTS. The initial scan of chromosome 17 performed on two large pedigrees provided a nonparametric LOD score of 2.41 near D17S928. Fine mapping with 17 additional microsatellite markers increased the peak to 2.61 (P=.002). The original families, as well as two additional pedigrees, were genotyped for 25 single-nucleotide polymorphisms (SNPs), with a focus on three genes in the indicated region that could play a role in the development of GTS, on the basis of their function and expression profile. Multiple three-marker haplotypes spanning all three genes studied provided highly significant association results (P<.001). An independent sample of 96 small families with one or two children affected with GTS was also studied. Of the 25 SNPs, 3 were associated with GTS at a statistically significant level. The transmission/disequilibrium test for a three-site haplotype moving window again provided multiple positive results. The background linkage disequilibrium (LD) of the region was studied in eight populations of European origin. A complicated pattern was revealed, with the pairwise tests producing unexpectedly high LD values at the telomeric TBCD gene. In conclusion, our findings warrant the further investigation of 17q25 as a candidate susceptibility region for GTS.