Segregation analysis reveals evidence of a major gene for Alzheimer disease.

In an attempt to resolve the relative influences of major genes, multifactorial heritability, and cohort effects on the susceptibility to Alzheimer disease (AD), complex segregation analysis was performed on 232 nuclear families. All families were consecutively ascertained through a single proband who was referred for diagnostic evaluation of a memory disorder. The results suggest that susceptibility to AD is determined, in part, by a major autosomal dominant allele with an additional multifactorial component. Single-locus, polygenic, sporadic, and no-transmission models, as well as recessive inheritance of the major effect, were significantly rejected. Excess transmission from the heterozygote was marginally significant and probably reflects the presence of phenocopies or perhaps the existence of two or more major loci for AD. The frequency of the AD susceptibility allele was estimated to be .038, but the major locus accounts for only 24% of the transmission variance, indicating a substantial role for other genetic and nongenetic mechanisms in the causation of AD.     

Complex segregation analysis of non-myoclonic idiopathic generalized epilepsy in families ascertained from probands affected with idiopathic epilepsy with tonic-clonic seizures in Antioquia,Colombia

In an attempt to identify the possible role of major genes, multifactorial inheritance, and cohort effects in the susceptibility to idiopathic epilepsy with generalized tonic-clonic seizures of the awakening type (GTCS), complex segregation analysis was performed in 196 nuclear families ascertained through affected probands with idiopathic epilepsy with GTCS belonging to the Paisa community of Antioquia (Colombia). Models postulating no transmission, single major locus (dominant and recessive) only, and multifactorial component only, were rejected. Since the codominant single major locus model could not be rejected and models that assign no major locus to transmission, no polygenic component to transmission, and no transmission of the major effect were rejected, complex segregation analysis suggested that a major autosomal codominant allele together with a multifactorial component (mixed model) best explained clustering of idiopathic epilepsy with GTCS in families of the Paisa community. The deficit of transmission of heterozygotes (0.17) is compatible with the existence of epistasis acting on a major gene whose frequency was estimated to be 0.0211. Its transmission variance accounts for 81% of the susceptibility to idiopathic epilepsy with GTCS. The complementary variance (19%) is due to the polygenic component. Received: 19 January 1996 / Revised: 11 March 1996     

Detection of major gene for Gilles de la Tourette syndrome

The families of 250 consecutive, unselected patients with Tourette syndrome (TS) were analyzed. If the parents had either motor or vocal tics, but not both, there was an increased risk of both TS and tics in the offspring. The mode of inheritance of the combined tic-Tourette trait was evaluated in both nuclear families and extended pedigrees. Complex segregation analysis was carried out allowing for possible contributions from both a major autosomal locus and multifactorial inheritance of variation in the background of each genotype. The most likely mode of inheritance was a major semidominant gene, Ts, with low heritability of the multifactorial background variation. This was true regardless of assumptions about the prevalence of the disorder. The hypothesis of strict multifactorial inheritance could not be rejected with nuclear family data alone. However, the hypothesis of no major gene effect was rejected using data on 3 generations for any estimate of lifetime risk less than 12 per 1,000 in the general population. A pure recessive major gene effect was also rejected. With a gene frequency of approximately .5%, the penetrance was estimated to be about 94% in abnormal Ts/Ts homozygotes, 50% in Ts/ts heterozygotes, and less than 0.3% in normal ts/ts homozygotes. More than two of every three cases are heterozygotes, and nearly all other cases are phenocopies or new mutations. This is the first demonstration by segregation analysis of a major gene in a human neuropsychiatric disorder with a frequency approaching 1% of the population.     

A combined segregation and linkage analysis of insulin-dependent diabetes mellitus.

In an effort to clarify the mode of inheritance of insulin-dependent diabetes mellitus (IDDM), a total of 230 nuclear families with pointers were analyzed using the computer program COMBIN. Each family was ascertained without deliberate selection for multiplex families, and most families were completely typed for HLA-B, HLA-DR, and properdin factor B (Bf). There were 186 families with normal parents, 44 families with one affected parent, and no families with two affected parents. The computer program COMBIN evaluates evidence for a major locus of disease susceptibility, linkage of the major locus to a known genetic marker locus, linkage disequilibrium between the marker haplotypes and disease susceptibility, pleiotropic effects, and presence of an unlinked modifier. The parameters of COMBIN are T, Q, and D, representing the displacement, gene frequency of the IDDM allele, and dominance, respectively, of the major locus--and TM, QM, and DM being the analogous parameters of the modifier. In addition, the recombination fraction, theta, between the IDDM locus and HLA as well as the coupling frequencies are estimated. Finally, COMBIN simultaneously performs segregation and linkage analysis, with the optimal model being adjusted by the fit to the haplotype sharing distribution of IDDM. The results of these analyses indicated that the best-fitting genetic model of diabetic susceptibility appears to be a single major locus with near recessivity on a scale of standardized genetic liability, with gene frequency of the IDDM susceptibility allele of approximately 14%. In addition, the recombination fraction between the major locus and HLA is zero in all models; that is, for the B-BF-DR haplotype, the IDDM locus is tightly linked, probably (according to data from previous studies) to HLA-DR. Information determined by magnitude of coupling frequencies indicated that there is significant positive linkage disequilibrium with the haplotypes B8-BfS-DR4 and B15-BfS-DR4, significant negative linkage disequilibrium with B7-BfS-DR2, and intermediate disequilibrium for B8-BfS-DR3, B18-BfF1-DR3, and B40-BfS-DR4. Significant evidence in favor of an unlinked (to HLA) modifier (either single major locus or polygenes) could not be demonstrated. In conclusion, genetic susceptibility to IDDM appears to be most consistent with a single major locus with near recessivity that is tightly linked to HLA.     

Evidence, from family studies, for linkage disequilibrium between TGFA and a gene for nonsyndromic cleft lip with or without cleft palate.

The inheritance of alleles of the transforming growth factor alpha (TGFA) locus has been studied in families affected with cleft lip with or without cleft palate (CL/P), by using the transmission/disequilibrium test described by Spielman and colleagues. Only heterozygous parents with an affected child can be included in this test, but within such families a significantly greater frequency of C2 alleles were transmitted to affected children than would be expected by chance. There was no evidence that the total number of C2 alleles transmitted to affected and unaffected children differed significantly from random segregation. These data provide evidence from within families that a gene for susceptibility to CL/P is in significant linkage disequilibrium with the C2 allele of the TGFA locus.     

How do homozygous parents affect TDT as a test for association?

The traditional transmission disequilibrium test (TDT) (Spielman et al., 1993) is a powerful test for association only in the presence of linkage. Since allele transmissions from homozygous parents do not carry any information on linkage, the TDT statistic uses data only on heterozygous parents. However, homozygous parents carry information on association between alleles at a marker locus and a disease locus. In this article, we explore whether inclusion of homozygous parents increases the power to detect association. The resultant test statistic follows a chi(2) distribution with 2 degrees of freedom. Monte-Carlo simulations are included to compare the performance of this test with the traditional TDT under different disease models.     

Segregation Analysis of Body Mass Index in an Unselected French-Canadian Sample: The Québec Family Study

Interest in a single gene etiology for obesity, as assessed by the body mass index (BMI), has been spurred recently by reports of a putative recessive major gene for extreme values, which accounts for as much as 40% of the variance. The major gene hypothesis was evaluated here in the Québec Family Study, a random sample of 375 French-Canadian volunteer families. This report represents one component in a more complete investigation of obesity in these families. In contrast to the recent studies, a major gene hypothesis for BMI was not verified here. Although there was a major effect, it did not conform to a Mendelian pattern of transmission. A multifactorial component (i.e., polygenic and/or common environmental factors) accounted for 42% of the phenotypic variance. In addition, evidence of heterogeneity between the generations was found. The heterogeneity was traced to the major non-Mendelian component (which accounted for 0.01% of the variance in parents and over 40% in offspring) rather than to the multifactorial one. These results would suggest that a simple recessive gene mixed model may not be sufficient to explain the familial distribution of the BMI. Several factors which may have contributed to these results include temporal trends and surrogate effects such as those related to variation in body composition and energy balance components. (OBESITY RESEARCH 1993; 1:288–294)     

Platelet monoamine oxidase (MAO) activity: evidence for a single major locus.

Monoamine oxidase (MAO), a mitochondrial enzyme involved in the degradation of biogenic amines, has been associated with psychiatric morbidity. Although twin and family studies have indicated that MAO activity is familial, the exact mode of transmission is unclear. We performed segregation analysis on 154 nuclear families containing 419 individuals using the mixed model, which allows for a single major locus with a polygenic background. We were able to reject a dominant and additive locus with or without a heritable background and a recessive locus without background. The acceptable models were: (1) a codominant model without background where the mean of the heterozygote distribution was 30% of the distance from the low to the high homozygote distributions, and (2) a recessive locus with heritable background. In both cases, the gene frequency for the high-MAO allele is approximately .25--at odds with suggestions that low-MAO represents a genetic marker for a disorder such as schizophrenia with a lifetime risk of only 0.85%. To ensure that results were not artifacts from a familial, skewed distribution, the data were also analyzed after power transformation. In addition, hypotheses were tested using both the joint and conditional likelihoods to examine for possible misspecification of the model with respect to intergenerational differences. Finally, we allowed for non-Mendelian transmission probabilities to provide another class of alternatives against which to test the hypothesis of a major locus. All these approaches provided additional confirmation for the presence of a major locus segregating within these families.     

Erythrocyte phosphoglucomutase: A family study of a PGM1 deficient allele

Summary We have observed a large Mexican American family segregating for a low activity allele at the phosphoglucomutase-1 locus. The deficient allele is detectable by starch gel electrophoresis and by direct activity determination. The presence of the deficient allele in either the homozygous or heterozygous condition is not associated with any other phenotypic finding.     

Fine mapping of IGAD1 in IgA deficiency and common variable immunodeficiency: identification and characterization of haplotypes shared by affected members of 101 multiple-case families

To limit the region containing a mutation predisposing to selective IgA deficiency (IgAD) and common variable immunodeficiency (CVID), 554 informative members of 101 multiple-case families were haplotyped at the IGAD1 candidate locus in the MHC. Microsatellite markers were placed onto the physical map of IGAD1 to establish their order and permit rapid haplotype analyses. Linkage analysis of this extended family set provided additional support for a strong susceptibility locus at IGAD1 with a maximum multipoint nonparametric linkage score in excess of 3. Although the transmission of maternal IGAD1 haplotypes from unaffected heterozygous parents to the affected offspring was in excess, this was not apparent in multiple-case families with a predominance of affected mothers, suggesting that this parental bias is influenced by the affection status of transmitting parents and supporting a maternal effect in disease susceptibility. Of 110 haplotypes shared by 258 affected family members, a single haplotype (H1) was found in 44 pairs of affected relatives, accounting for the majority of the IGAD1 contribution to the development of IgAD/CVID in our families. The H1 allelic variability was higher in the telomeric part of the class III region than in the distal part of the class II region in both single- and multiple-case families. Incomplete H1 haplotypes had most variant alleles in the telomeric part of the analyzed region in homozygous IgAD/CVID patients, whereas this was not observed in unaffected homozygotes. These data suggest that a telomeric part of the class II region or centromeric part of the class III region is the most likely location of IGAD1.     

Germline integration of moloney murine leukemia virus at the Mov13 locus leads to recessive lethal mutation and early embryonic death

Thirteen mouse substrains genetically transmitting the exogenous Moloney murine leukemia virus (M-MuLV) at a single locus (Mov locus) have been derived previously. Experiments were performed to investigate whether homozygosity at the Mov loci would be compatible with normal development. Animals heterozygous at an Mov locus were mated, and the genotype of the offspring was analyzed. From parents heterozygous at the loci Mov1 to Mov12, respectively, homozygous offspring were obtained with the expected Mendelian frequency. In contrast, no homozygous offspring or embryos older than day 15 of gestation were obtained from parents heterozygous at the Mov13 locus. When pregnant Mov13 females at day 13 and day 14 of gestation were analyzed, approximately 25% of the embryos were degenerated. Genotyping revealed that these degenerated embryos were invariably homozygous and the normal appearing embryos were either heterozygous or negative for M-MuLV. These results suggest that integration of M-MuLV at the Mov13 locus leads to insertion mutagenesis, resulting in embryonic arrest between day 12 and day 13 of gestation. It is possible that the Mov13 locus represents a gene or gene complex involved in the early embryonic development of the mouse.     

HLA-linked rheumatoid arthritis.

Twenty-eight pedigrees were ascertained through pairs of first-degree relatives diagnosed with rheumatoid arthritis (RA). RA was confirmed in 77 pedigree members including probands; the absence of disease was verified in an additional 261 pedigree members. Pedigree members were serologically typed for HLA. We used likelihood analysis to statistically characterize the HLA-linked RA susceptibility locus. The genetic model assumed tight linkage to HLA. The analysis supported the existence of an HLA-linked RA susceptibility locus, estimated the susceptibility allele frequency as 2.16%, and estimated the lifetime penetrance as 41% in male homozygotes and as 48% in female homozygotes. Inheritance was recessive in males and was nearly recessive in females. In addition, the analysis attributed 78% of the variance within genotypes to genetic or environmental effects shared by siblings. The genetic model inferred in this analysis is consistent with previous association, linkage, and familial aggregation studies of RA. The inferred HLA-linked RA susceptibility locus accounts for approximately one-half of familial RA, although it accounts for only approximately one-fifth of the RA in the population. Although other genes may account for the remaining familial RA, a large portion of RA cases may occur sporadically.     

Association studies of QTL for multi-allele markers by mixed models

In this paper, we extend association study methods of both Fan et al. [Hum Hered 2002;53:130-145], in which a quantitative trait locus (QTL) and a multi-allele marker are considered for trio families, and Fan and Xiong [Biostatistics 2003, in press], in which a QTL and a bi-allelic marker are considered for nuclear families. The objective is to build mixed models for association study between a QTL and a multi-allelic marker for nuclear families with any number of offspring. Two types of nuclear family data are considered: the first is genetic data of offspring from at least one heterozygous parents, and the second is genetic data of offspring of nuclear family. (1) For the data of offspring from at least one heterozygous parents, we assume that at least one parent is heterozygous at the marker locus, and we may infer clearly the transmission of parental marker alleles to the offspring. We show that it can be used in association study in the presence of linkage. The theoretical basis is the difference between the conditional mean of trait value given an allele is transmitted and the conditional mean of trait value given the allele is not transmitted from a heterozygous parent. To build valid models, we calculate the variance covariance structure of trait values of offspring. Besides, the reduction of the number of parameters is discussed under an assumption of tight linkage between the trait locus and the marker. (2) For the data of offspring of nuclear family, we show that it can be used in general association study. In this case, the theoretical basis is the difference between the conditional mean of trait values given an allele is transmitted from a parent and the population mean. Then, we calculate variance-covariance structure of trait values of offspring. (3) Based on the theoretical analysis, mixed models are built for each type of the data, and related test statistics are proposed for association study. By power calculation and comparison, we show that, in some instances, the proposed test statistics have higher power than that by collapsing alleles to be new ones. The proposed models are used to analyze chromosomes 4 and chromosome 16 data of the Oxford asthma data, Genetic Analysis Workshop 12.     

Fine mapping of the diabetes-susceptibility locus, IDDM4, on chromosome 11q13.

Genomewide linkage studies of type 1 diabetes (or insulin-dependent diabetes mellitus [IDDM]) indicate that several unlinked susceptibility loci can explain the clustering of the disease in families. One such locus has been mapped to chromosome 11q13 (IDDM4). In the present report we have analyzed 707 affected sib pairs, obtaining a peak multipoint maximum LOD score (MLS) of 2.7 (lambda(s)=1.09) with linkage (MLS>=0.7) extending over a 15-cM region. The problem is, therefore, to fine map the locus to permit structural analysis of positional candidate genes. In a two-stage approach, we first scanned the 15-cM linked region for increased or decreased transmission, from heterozygous parents to affected siblings in 340 families, of the three most common alleles of each of 12 microsatellite loci. One of the 36 alleles showed decreased transmission (50% expected, 45.1% observed [P=.02, corrected P=.72]) at marker D11S1917. Analysis of an additional 1,702 families provided further support for negative transmission (48%) of D11S1917 allele 3 to affected offspring and positive transmission (55%) to unaffected siblings (test of heterogeneity P=3x10-4, corrected P=. 01]). A second polymorphic marker, H0570polyA, was isolated from a cosmid clone containing D11S1917, and genotyping of 2,042 families revealed strong linkage disequilibrium between the two markers (15 kb apart), with a specific haplotype, D11S1917*03-H0570polyA*02, showing decreased transmission (46.4%) to affected offspring and increased transmission (56.6%) to unaffected siblings (test of heterogeneity P=1.5x10-6, corrected P=4.3x10-4). These results not only provide sufficient justification for analysis of the gene content of the D11S1917 region for positional candidates but also show that, in the mapping of genes for common multifactorial diseases, analysis of both affected and unaffected siblings is of value and that both predisposing and nonpredisposing alleles should be anticipated.     

Mapping and Validation of the Major Sex-Determining Region in Nile Tilapia (Oreochromis niloticus L.) Using RAD Sequencing

Sex in Oreochromis niloticus (Nile tilapia) is principally determined by an XX/XY locus but other genetic and environmental factors also influence sex ratio. Restriction Associated DNA (RAD) sequencing was used in two families derived from crossing XY males with females from an isogenic clonal line, in order to identify Single Nucleotide Polymorphisms (SNPs) and map the sex-determining region(s). We constructed a linkage map with 3,802 SNPs, which corresponded to 3,280 informative markers, and identified a major sex-determining region on linkage group 1, explaining nearly 96% of the phenotypic variance. This sex-determining region was mapped in a 2 cM interval, corresponding to approximately 1.2 Mb in the O. niloticus draft genome. In order to validate this, a diverse family (4 families; 96 individuals in total) and population (40 broodstock individuals) test panel were genotyped for five of the SNPs showing the highest association with phenotypic sex. From the expanded data set, SNPs Oni23063 and Oni28137 showed the highest association, which persisted both in the case of family and population data. Across the entire dataset all females were found to be homozygous for these two SNPs. Males were heterozygous, with the exception of five individuals in the population and two in the family dataset. These fish possessed the homozygous genotype expected of females. Progeny sex ratios (over 95% females) from two of the males with the “female” genotype indicated that they were neomales (XX males). Sex reversal induced by elevated temperature during sexual differentiation also resulted in phenotypic males with the “female” genotype. This study narrows down the region containing the main sex-determining locus, and provides genetic markers tightly linked to this locus, with an association that persisted across the population. These markers will be of use in refining the production of genetically male O. niloticus for aquaculture.     

A sib-pair strategy for the use of restriction fragment length polymorphisms to study the mode of transmission of type II diabetes.

Three models for explaining the joint distribution of DNA insertional elements and type II diabetes in sibships are given as potential candidates for resolving the mode of transmission for this common non-Mendelian disorder. While the distributions predicted from all models are subject to the same rigid marginal constraints, only two models can be considered genetic (a combined major locus multifactorial background model and a two-locus model). A sequential probability ratio test is proposed to distinguish between these two models. Once the probands have been ascertained and using realistic parameter estimates, it is shown that, on the average, fewer than 100 affected sib-pairs are required to reach a decision at the 1% significance level with 99% power.     

Identifying quantitative trait locus by genetic background interactions in association studies

Association studies are designed to identify main effects of alleles across a potentially wide range of genetic backgrounds. To control for spurious associations, effects of the genetic background itself are often incorporated into the linear model, either in the form of subpopulation effects in the case of structure or in the form of genetic relationship matrices in the case of complex pedigrees. In this context epistatic interactions between loci can be captured as an interaction effect between the associated locus and the genetic background. In this study I developed genetic and statistical models to tie the locus by genetic background interaction idea back to more standard concepts of epistasis when genetic background is modeled using an additive relationship matrix. I also simulated epistatic interactions in four-generation randomly mating pedigrees and evaluated the ability of the statistical models to identify when a biallelic associated locus was epistatic to other loci. Under additive-by-additive epistasis, when interaction effects of the associated locus were quite large (explaining 20% of the phenotypic variance), epistasis was detected in 79% of pedigrees containing 320 individuals. The epistatic model also predicted the genotypic value of progeny better than a standard additive model in 78% of simulations. When interaction effects were smaller (although still fairly large, explaining 5% of the phenotypic variance), epistasis was detected in only 9% of pedigrees containing 320 individuals and the epistatic and additive models were equally effective at predicting the genotypic values of progeny. Epistasis was detected with the same power whether the overall epistatic effect was the result of a single pairwise interaction or the sum of nine pairwise interactions, each generating one ninth of the epistatic variance. The power to detect epistasis was highest (94%) at low QTL minor allele frequency, fell to a minimum (60%) at minor allele frequency of about 0.2, and then plateaued at about 80% as alleles reached intermediate frequencies. The power to detect epistasis declined when the linkage disequilibrium between the DNA marker and the functional polymorphism was not complete.     

Associative overdominance: Evaluating the effects of inbreeding and linkage disequilibrium

Expressions are obtained for the expected phenotypic values of homozygous and heterozygous genotypes for a neutral marker locus linked to a locus segregating for a recessive deleterious gene. The phenotypic values are functions of the allele frequencies at the marker locus, the inbreeding coefficient and the degree of association of the deleterious gene with the marker alleles. The analysis is extended to more than two alleles at the marker locus. Either linkage disequilibrium or inbreeding alone can produce an apparent superiority of heterozygotes for the marker locus (unless specified otherwise, the terms ‘homozygote’ and ‘heterozygote’ will refer to the marker locus). The effect of linkage disequilibrium on the difference between the heterozygote and homozygote values can be positive (associative overdominance) or negative (associative underdominance), depending on the frequencies of the marker alleles and the degree of their association with the deleterious gene. Inbreeding has always a positive effect. In general, the expected value of a homozygote is a positive function of its allele frequency. When the various homozygous genotypes are combined into one class and the various heterozygous genotypes into another, the phenotypic difference of the two classes is a function of the evenness of the allelic frequency distribution. Inbreeding is a more likely explanation of associative overdominance if the frequency of the deleterious gene is low, but its effect on the character high. Conversely, linkage disequilibrium is more likely if the frequency is high and the effect low. The degrees of association between marker alleles and the deleterious gene can, in principle, be estimated from the observed phenotypic scores and used to calculate expected multi-locus genotype scores. This could provide the basis for statistical tests of the associative overdominance hypothesis as an explanation of observed correlations between multi-locus heterozygosity and phenotypic traits.