Using lod-score differences to determine mode of inheritance: a simple, robust method even in the presence of heterogeneity and reduced penetrance.

Determining the mode of inheritance is often difficult under the best of circumstances, but when segregation analysis is used, the problems of ambiguous ascertainment procedures, reduced penetrance, heterogeneity, and misdiagnosis make mode-of-inheritance determinations even more unreliable. The mode of inheritance can also be determined using a linkage-based method (maximized maximum lod score or mod score) and association-based methods, which can overcome many of these problems. In this work, we determined how much information is necessary to reliably determine the mode of inheritance from linkage data when heterogeneity and reduced penetrance are present in the data set. We generated data sets under both dominant and recessive inheritance with reduced penetrance and with varying fractions of linked and unlinked families. We then analyzed those data sets, assuming reduced penetrance, both dominant and recessive inheritance, and no heterogeneity. We investigated the reliability of two methods for determining the mode of inheritance from the linkage data. The first method examined the difference (delta) between the maximum lod scores calculated under the two mode-of-inheritance assumptions. We found that if delta was > 1.5, then the higher of the two maximum lod scores reflected the correct mode of inheritance with high reliability and that a delta of 2.5 appeared to practically guarantee a correct mode-of-inheritance inference. Furthermore, this reliability appeared to be virtually independent of alpha, the fraction of linked families in the data set, although the reliability decreased slightly as alpha fell below .50.(ABSTRACT TRUNCATED AT 250 WORDS)     

Direct power comparisons between simple LOD scores and NPL scores for linkage analysis in complex diseases.

Several methods have been proposed for linkage analysis of complex traits with unknown mode of inheritance. These methods include the LOD score maximized over disease models (MMLS) and the "nonparametric" linkage (NPL) statistic. In previous work, we evaluated the increase of type I error when maximizing over two or more genetic models, and we compared the power of MMLS to detect linkage, in a number of complex modes of inheritance, with analysis assuming the true model. In the present study, we compare MMLS and NPL directly. We simulated 100 data sets with 20 families each, using 26 generating models: (1) 4 intermediate models (penetrance of heterozygote between that of the two homozygotes); (2) 6 two-locus additive models; and (3) 16 two-locus heterogeneity models (admixture alpha = 1.0,.7,.5, and.3; alpha = 1.0 replicates simple Mendelian models). For LOD scores, we assumed dominant and recessive inheritance with 50% penetrance. We took the higher of the two maximum LOD scores and subtracted 0.3 to correct for multiple tests (MMLS-C). We compared expected maximum LOD scores and power, using MMLS-C and NPL as well as the true model. Since NPL uses only the affected family members, we also performed an affecteds-only analysis using MMLS-C. The MMLS-C was both uniformly more powerful than NPL for most cases we examined, except when linkage information was low, and close to the results for the true model under locus heterogeneity. We still found better power for the MMLS-C compared with NPL in affecteds-only analysis. The results show that use of two simple modes of inheritance at a fixed penetrance can have more power than NPL when the trait mode of inheritance is complex and when there is heterogeneity in the data set.     

The HLA system and leprosy in Thailand

Summary To investigate immunogenetics of leprosy, 205 leprosy patients (26 with tuberculoid, 57 with borderline-tuberculoid, 21 with borderline, 31 with borderline-lepromatous, and 70 with lepromatous leprosy) have been typed for HLA antigens, and compared with 183 healthy controls from the same region (Nothern Thailand). There was no significant difference between the overal group of leprosy patients or the three borderline classes and the controls. The two polar forms, tuberculoid and lepromatous leprosy, however, showed significant associations: HLA-A2 is decreased and HLA-Bw17 is increased in tuberculoid leprosy; HLA-B7 is increased in lepromatous leprosy. When both polar forms are compared with each other, HLA-A2 is significantly higher, HLA-Bw40 lower in patients with lepromatous than in those with tuberculoid leprosy. The results are discussed with respect to the different immune responsiveness in the two polar forms of leprosy.     

The modes of inheritance of insulin-dependent diabetes mellitus or the genetics of IDDM, no longer a nightmare but still a headache.

The discovery of HLA antigen associations with juvenile-type insulin-dependent diabetes mellitus (IDDM) provided strong evidence separating this disorder, or group of disorders, from maturity-type noninsulin-dependent diabetes, as well as adding to the evidence for an immunologic pathogenesis. In addition, it was hoped that the use of these disease-marker associations in appropriate studies might clarify the genetics of IDDM. While these associations have provided a useful tool to further investigate the genetics and pathogenesis of IDDM, the mode or modes of inheritance of this group of disorders remain an area of great controversy. Susceptibility to IDDM is currently being proposed as being inherited as a single autosomal dominant, as a single autosomal recessive, as recessive and some dominant forms, in an intermediate gene dosage model, in a heterogeneous three-allele or two HLA loci model, and as a two-locus disorder. The arguments for each of these proposals is presented, as well as the problems of each. We surmise that the weight of evidence supports the heterogeneity hypothesis but that the modes of inheritance of IDDM will be fully resolved only when we can more reliably identify the diabetogenic genotype, rather than being limited in our investigations to the study of only full-blown clinical disease.     

The affected sib method. III. Selection and recombination.

The affected sib-pair method has been used to investigate the mode of inheritance, and to estimate the "disease" allele frequency, for a number of HLA-associated diseases. One of the assumptions of the original sib-pair method is that the disease confers no selective disadvantage on affected individuals. This is obviously not the situation for most diseases. We have determined the expected HLA haplotype-sharing distribution among affected sib-pairs when selection against individuals with the disease is taken into account. We have shown that if the mode of inheritance of the selectively disadvantageous disease is recessive or additive, the original affected sib-pair analysis, ignoring selection, still estimates the true mode of inheritance, but usually yields an underestimate of the "disease" allele frequency. For intermediate and dominant models of disease predisposition, both the estimates of the degree of penetrance of the "disease" genotypes, and the "disease" allele frequency, are altered if selection is ignored in the analysis. Similarly, allowing for recombination between the "disease" locus and the HLA region does not affect the determination of the mode of inheritance of the disease if it is recessive or additive; in other cases, however, the estimate of the mode of inheritance is affected. The "disease" allele frequency is overestimated when nonzero recombination is ignored for all the modes of inheritance that have been studied.     

Reproducibility and complications in gene searches: linkage on chromosome 6, heterogeneity, association, and maternal inheritance in juvenile myoclonic epilepsy

Evidence for genetic influences in epilepsy is strong, but reports identifying specific chromosomal origins of those influences conflict. One early study reported that human leukocyte antigen (HLA) markers were genetically linked to juvenile myoclonic epilepsy (JME); this was confirmed in a later study. Other reports did not find linkage to HLA markers. One found evidence of linkage to markers on chromosome 15, another to markers on chromosome 6, centromeric to HLA. We identified families through a patient with JME and genotyped markers throughout chromosome 6. Linkage analysis assuming equal male-female recombination probabilities showed evidence for linkage (LOD score 2.5), but at a high recombination fraction (theta), suggesting heterogeneity. When linkage analysis was redone to allow independent male-female thetas, the LOD score was significantly higher (4.2) at a male-female theta of.5,.01. Although the overall pattern of LOD scores with respect to male-female theta could not be explained solely by heterogeneity, the presence of heterogeneity and predominantly maternal inheritance of JME might explain it. By analyzing loci between HLA-DP and HLA-DR and stratifying the families on the basis of evidence for or against linkage, we were able to show evidence of heterogeneity within JME and to propose a marker associated with the linked form. These data also suggest that JME may be predominantly maternally inherited and that the HLA-linked form is more likely to occur in families of European origin.     

Inter- and intrafamilial heterogeneity: effective sampling strategies and comparison of analysis methods.

Heterogeneity, both inter- and intrafamilial, represents a serious problem in linkage studies of common complex diseases. In this study we simulated different scenarios with families who phenotypically have identical diseases but who genotypically have two different forms of the disease (both forms genetic). We examined the proportion of families displaying intrafamilial heterogeneity, as a function of mode of inheritance, gene frequency, penetrance, and sampling strategies. Furthermore, we compared two different ways of analyzing linkage in these data sets: a two-locus (2L) analysis versus a one-locus (SL) analysis combined with an admixture test. Data were simulated with tight linkage between one disease locus and a marker locus; the other disease locus was not linked to a marker. Our findings are as follows: (1) In contrast to what has been proposed elsewhere to minimize heterogeneity, sampling only "high-density" pedigrees will increase the proportion of families with intrafamilial heterogeneity, especially when the two forms are relatively close in frequency. (2) When one form is dominant and one is recessive, this sampling strategy will greatly decrease the proportions of families with a recessive form and may therefore make it more difficult to detect linkage to the recessive form. (3) An SL analysis combined with an admixture test achieves about the same lod scores and estimate of the recombination fraction as does a 2L analysis. Also, a 2L analysis of a sample of families with intrafamilial heterogeneity does not perform significantly better than an SL analysis. (4) Bilineal pedigrees have little effect on the mean maximum lod score and mean maximum recombination fraction, and therefore there is little danger that including these families will lead to a false exclusion of linkage.     

Linkage analysis of idiopathic generalized epilepsy (IGE) and marker loci on chromosome 6p in families of patients with juvenile myoclonic epilepsy: no evidence for an epilepsy locus in the HLA region.

Evidence for a locus (EJM1) in the HLA region of chromosome 6p predisposing to idiopathic generalized epilepsy (IGE) in the families of patients with juvenile myoclonic epilepsy (JME) has been obtained in two previous studies of separately ascertained groups of kindreds. Linkage analysis has been undertaken in a third set of 25 families including a patient with JME and at least one first-degree relative with IGE. Family members were typed for eight polymorphic loci on chromosome 6p: F13A, D6S89, D6S109, D6S105, D6S10, C4B, DQA1/A2, and TCTE1. Pairwise and multipoint linkage analysis was carried out assuming autosomal dominant and autosomal recessive inheritance and age-dependent high or low penetrance. No significant evidence in favor of linkage was obtained at any locus. Multipoint linkage analysis generated significant exclusion data (lod score < -2.0) at HLA and for a region 10-30 cM telomeric to HLA, the extent of which varied with the level of penetrance assumed. These observations indicate that genetic heterogeneity exists within this epilepsy phenotype.     

Evidence for recessive and against dominant inheritance at the HLA-"linked" locus in coeliac disease.

It has been proposed that gluten sensitive enteropathy (GSE) results from the interaction of two loci: one locus linked to HLA and associated with dominant inheritance, and the other, a non-HLA-linked GSE-associated B-cell alloantigen, exhibiting recessive inheritance. We have shown in previous analyses that a two-locus, dominant-recessive model is less compatible with the existing population prevalence and observed familial segregation data than is a recessive-recessive two-locus model. Here we present additional analyses of reported population and familial HLA data that support the recessive mode of inheritance for the HLA-linked disease locus. Reported data from HLA typing of affected sib pairs, the association of GSE with DR3 and DR7 in different populations, and the proportions of different HLA phenotypes and genotypes were compared with expected data derived by three different methods. The HLA data analyses consistently reject a dominant mode of inheritance for the presumed HLA-linked disease allele but do not reject a recessive model. The affected sib-pair data also support a recessive model. These analyses are consistent with our previous prediction that the HLA-"linked" disease allele in GSE is recessive inherited.     

The genetic epidemiology of leprosy in a Brazilian population.

Data on leprosy patients have been obtained from the Dispensary of Leprosy of Campinas, São Paulo, where records on practically all cases of leprosy in the Campinas area during the period 1960-70 are filed. The whole sample comprises 10,886 individuals, distributed among 1,568 families. Complex segregation analysis was utilized to determine the nature of the genetic factors that may operate on leprosy and its subtypes. The results suggest the presence of a recessive major gene controlling susceptibility to leprosy per se, with frequency of approximately .05, although there are deviations from the expected Mendelian segregation proportions. Possible etiologic heterogeneity was examined by considering two subtypes separately: for lepromatous leprosy and tuberculoid leprosy there are suggestions for a segregating major effect; however, Mendelian transmission could not be demonstrated in either case. Therefore, there is no evidence to suggest unique genetic determinants for leprosy subtypes.     

Genetic analysis of multiply-affected families of insulin-dependent diabetes mellitus (IDDM) probands

The present study combines segregation and linkage information on 30 families ascertained through a proband and a first degree relative affected with insulin-dependent diabetes mellitus (IDDM). An autosomal dominant model with incomplete penetrance was much more likely to fit the family data than a recessive model, whether or not linkage to HLA was assumed. The lod scores for linkage to HLA were 2.46 at theta M = theta F = 0.00 for dominant and 1.45 at theta M = theta F = 0.22 for a recessive model. The results are discussed in light of heterogeneity in likelihood and lod scores when the families are grouped by familial types, which indicate that the increase in likelihood of a dominant hypothesis can be attributed to the parent-child families and not the sib-sib families.     

HLA and the inheritance of multiple sclerosis: linkage analysis of 72 pedigrees.

Linkage analysis of 72 pedigrees by the maximum-likelihood method provides evidence of linkage between HLA and the hypothesized multiple sclerosis susceptibility gene (MSSG) for both the dominant and recessive models of inheritance and for penetrance values ranging from 5%--65% (or higher). This MSSG, if it exists, is most likely located at 15%--20% recombination units from the HLA complex, probably on the B-D side. The analysis also shows that there is no heterogeneity in the estimates of linkage, and lod scores are not artifically inflated because of the association of multiple sclerosis (MS) with HLA-B7.     

Complex segregation analysis of Parkinson's disease in the Finnish population

The risk of Parkinson's disease (PD) is higher among relatives of affected individuals than among other members of the population, and most family studies have suggested autosomal dominant inheritance, although both autosomal dominant and recessive susceptibility genes have recently been identified. We carried out a complex segregation analysis with POINTER to assess the mode of inheritance of PD in the population of northern Finland. Nuclear families (n=265) were identified through a proband with idiopathic PD. The analysis was first carried out for the total data set, and then the heterogeneity between early-onset (proband under 55 years at onset) and late-onset families was examined. Finally, families with more than one affected individual were analyzed separately. The sporadic model was rejected (P<0.0001). Significant heterogeneity was found between the early-onset and late-onset families, suggesting that major genes have a greater role in early-onset PD than in late-onset PD and that the etiology of idiopathic PD is heterogeneous, even in the Finnish population, which has evolved from a small group of founders. The analysis of familial PD supported the hypothesis that a major locus was present in this subset, but it was not possible to distinguish between a recessive model with a high penetrance and a dominant model with lower penetrance.     

Nonparametric evaluation of quantitative traits in population-based association studies when the genetic model is unknown

Statistical association between a single nucleotide polymorphism (SNP) genotype and a quantitative trait in genome-wide association studies is usually assessed using a linear regression model, or, in the case of non-normally distributed trait values, using the Kruskal-Wallis test. While linear regression models assume an additive mode of inheritance via equi-distant genotype scores, Kruskal-Wallis test merely tests global differences in trait values associated with the three genotype groups. Both approaches thus exhibit suboptimal power when the underlying inheritance mode is dominant or recessive. Furthermore, these tests do not perform well in the common situations when only a few trait values are available in a rare genotype category (disbalance), or when the values associated with the three genotype categories exhibit unequal variance (variance heterogeneity). We propose a maximum test based on Marcus-type multiple contrast test for relative effect sizes. This test allows model-specific testing of either dominant, additive or recessive mode of inheritance, and it is robust against variance heterogeneity. We show how to obtain mode-specific simultaneous confidence intervals for the relative effect sizes to aid in interpreting the biological relevance of the results. Further, we discuss the use of a related all-pairwise comparisons contrast test with range preserving confidence intervals as an alternative to Kruskal-Wallis heterogeneity test. We applied the proposed maximum test to the Bogalusa Heart Study dataset, and gained a remarkable increase in the power to detect association, particularly for rare genotypes. Our simulation study also demonstrated that the proposed non-parametric tests control family-wise error rate in the presence of non-normality and variance heterogeneity contrary to the standard parametric approaches. We provide a publicly available R library nparcomp that can be used to estimate simultaneous confidence intervals or compatible multiplicity-adjusted p-values associated with the proposed maximum test.     

A susceptibility locus for migraine with aura, on chromosome 4q24

Migraine is a complex neurovascular disorder with substantial evidence supporting a genetic contribution. Prior attempts to localize susceptibility loci for common forms of migraine have not produced conclusive evidence of linkage or association. To date, no genomewide screen for migraine has been published. We report results from a genomewide screen of 50 multigenerational, clinically well-defined Finnish families showing intergenerational transmission of migraine with aura (MA). The families were screened using 350 polymorphic microsatellite markers, with an average intermarker distance of 11 cM. Significant evidence of linkage was found between the MA phenotype and marker D4S1647 on 4q24. Using parametric two-point linkage analysis and assuming a dominant mode of inheritance, we found for this marker a maximum LOD score of 4.20 under locus homogeneity (P=.000006) or locus heterogeneity (P=.000011). Multipoint parametric (HLOD = 4.45; P=.0000058) and nonparametric (NPL(all) = 3.43; P=.0007) analyses support linkage in this region. Statistically significant linkage was not observed in any other chromosomal region.     

Genetics of Type I diabetes mellitus: a single, recessive predisposition gene mapping between HLA-B and GLO. With an appendix on the estimation of selection bias.

Three different published sets of HLA-typed families of juvenile diabetes mellitus (JDM) patients have been analyzed. There was no significant genetic heterogeneity between them according to the criterion of Morton, and the total material was analyzed on the assumption of a single recessive (JDM-P) gene with incomplete penetrance. The analysis, carried out with the NYLIP program modified to account for penetrance less than 1 and for selection bias, yields highly significant lod scores for linkage between HLA and JDM-P, with a maximum value of 7.40 at theta = .05 +/- .03. The segregation of HLA and GLO in five affected sib pairs, in which one of the sibs carries an HLA/GLO recombinant, places JDM-P closer to HLA than the GLO locus: four of these five pairs are HLA-identical and GLO-different, in agreement with the conclusions of the formal linkage analysis. The data from these three independent sets of families are therefore consistent with our earlier claim that JDM is inherited as a recessive trait closely linked to HLA with reduced penetrance, and its analysis does not require more complicated genetic models.     

Image analysis to quantify S100-positive dendritic cells in leprosy-affected skin

A morphometric analysis of skin dendritic cells was done on biopsies of patients with different forms of leprosy. An anti S100 antibody was used to determine dendritic cell quantity and extension. Patients with a better immune response to the bacillus showed a greater number of dendritic cells in the cases of dimorphic tuberculoid leprosy and tuberculoid leprosy. This result contrasted with that from patients with dimorphic lepromatous leprosy and lepromatous leprosy.     

Linkage analysis identifies the thyroglobulin gene region as a major locus for familial congenital hypothyroidism

Congenital hypothyroidism affects 1/3000-4000 newborns and it has been estimated that 10-20% are familial cases with an autosomal recessive mode of inheritance. Previous studies of mostly individual cases have led to the identification of mutations in a number of genes, indicating that it is a genetically heterogeneous disease, but no major gene has been identified. In the present investigation, a population-based sample of 23 families with autosomal recessive congenital hypothyroidism, but no signs of goitre, were subject to linkage analysis. When markers located close to the thyroglobulin gene on chromosome 8q24 were used in a two-point analysis allowing for heterogeneity, a Z(max) of 4.10 was obtained with the microsatellite marker D8S557, indicating heterogeneity with 43% of the families being linked. A multipoint analysis using the markers D8S557 and D8S1835 gave a Z(max) of 3.51, assuming homogeneity. There was significant evidence of heterogeneity with 44.5% of the families being linked. The results indicate that a gene in 8q24 is a common cause of familial congenital hypothyroidism. Since thyroglobulin is essential for thyroid physiology, the gene encoding this protein is the obvious candidate for mutation analysis in the linked families.     

Genetic interrelationship between insulin-dependent diabetes mellitus, the autoimmune thyroid diseases, and rheumatoid arthritis.

To investigate the possible coinheritance of autoimmune diseases that are associated with the same HLA antigen, we studied 70 families in which at least two siblings had either type I diabetes mellitus (IDDM), autoimmune thyroid disease (ATD), rheumatoid arthritis (RA), or a combination of these diseases. HLA-A, B, and C typing was performed on all affected sibs in one generation or more. First, we estimated by sib-pair analysis the disease allele frequency (pD) and the mode of inheritance for each disease. According to the method of ascertainment entered into the analysis, the pD for ATD ranged from .120 to .180, for an additive (dominant) mode of inheritance. For RA, the pD ranged from .254 to .341, also for additive inheritance, although recessive inheritance could not be excluded. For IDDM, the pD ranged from .336 to .337 for recessive inheritance; additive inheritance was rejected. Second, we examined the distribution of shared parental haplotypes in pairs of siblings that were discordant for their autoimmune diseases. The results suggested that the same haplotype may predispose to both IDDM and ATD, or IDDM and RA, but not to both RA and ATD. Analysis of pedigrees supported this hypothesis. In 16 families typed for HLA-DR also, the haplotype predisposing to both IDDM and ATD was assigned from pedigree information to DR3 (44%), DR4 (39%), or DR5, DR6, or DR7 (5.5% each). In some families, these haplotypes segregated over several generations with ATD only (either clinical or subclinical), suggesting that in such families, ATD was a marker for a susceptibility to IDDM. In several families, an IDDM haplotype segregated with RA but not with ATD. This suggests that ATD- and RA-associated susceptibilities to IDDM may be biologically different and thus independently increase the risk of IDDM.     

The detection of linkage and heterogeneity in nuclear families for complex disorders: one versus two marker loci.

Using exact expected likelihoods, we have computed the average number of phase-unknown nuclear families needed to detect linkage and heterogeneity. We have examined the case of both dominant and recessive inheritance with reduced penetrance and phenocopies. Most of our calculations have been carried out under the assumption that 50% of families are linked to a marker locus. We have varied both the number of offspring per family and the sampling scheme. We have also investigated the increased power when the disease locus is midway between two marker loci 10 cM apart. For recessive inheritance, both linkage and heterogeneity can be detected in clinically feasible sample sizes. For dominant inheritance, linkage can be detected but heterogeneity cannot be detected unless larger sibships (four offspring) are sampled or two linked markers are available. As expected, if penetrance is reduced, sampling families with all sibs affected is most efficient. Our results provide a basis for estimating the amount of resources needed to find genes for complex disorders under conditions of heterogeneity.