No abstract

A more powerful robust test for linkage is developed from the methodology of Haseman and Elston [Behav Genet 2(1):3-19, 1972]. This new robust test uses weighted least-squares (WLS) methods to detect linkage between a quantitative trait and a polymorphic marker. For comparison, the characteristics of a test for linakge that uses known trait genotypes for the parents are also studied. Sample sizes needed to detect linkage, calculated using asymptotic results, are compared for 1) the usual Haseman-Elston method, 2) the WLS method, and 3) the method that uses parental trait genotype data. The WLS method needs at most twice the number of sib pairs as does the method that uses information on the trait genotypes of the parents. The small sample properties of the Haseman-Elston (H-E) and WLS tests are investigated by simulation. The power calculations for the H-E method are found to be accurate. The power of the WLS method is overestimated when fewer than 300 sib pairs are studied, but the WLS method is nonetheless more powerful than the usual H-E method. In samples of fewer than 300 sib pairs, the WLS test tends to be anticonservative. Treating all sib pairs from sibships of size 3 or 5 as independent does not increase the significance of the tests.

The robust method for detecting linkage developed by Haseman and Elston [The investigation of linkage between a quantitative trait and a marker locus. Behav Genet 2:3-19, 1972] for data from sib pairs is extended to any type of noninbred relative pair. The regression of the squared relative-pair trait difference on the estimated proportion of genes identical by descent (i.b.d.) at a marker locus is shown to depend upon the recombination fraction between the two loci; the regression coefficient is negative if the trait and marker loci are linked. A test for linkage based on data from any informative type of relative pair can thus be obtained by testing that this regression coefficient is less than zero. Formulae for the asymptotic power of such tests for linkage based upon independent relative pairs are developed. Results are also given for the special case in which the proportion of genes shared i.b.d. for relative pairs is known. Finally, a general algorithm is described that will incorporate all available pedigree data to calculate an estimate of the proportion of genes that a relative pair shares i.b.d. at a marker locus.

We fine mapped two primary biliary cirrhosis (PBC) risk loci, CLEC16A (C-type lectin domain family 16 member A)–suppressor of cytokine signaling 1 (SOCS1) and Spi-B protein (SPIB) and sequenced a locus, sialic acid acetylesterase (SIAE), proposed to harbor autoimmunity-associated mutations. In all, 1450 PBC cases and 2957 healthy controls were genotyped for 84 single-nucleotide polymorphisms (SNPs) across the CLEC16A-SOCS1 and SPIB loci. All 10 exons of the SIAE gene were resequenced in 381 cases and point substitutions of unknown significance assayed for activity and secretion. Fine mapping identified 26 SNPs across the CLEC16A-SOCS1 and 11 SNPs across the SPIB locus with significant association to PBC, the strongest signals at the CLEC16A-SOCS1 locus emanating from a SOCS1 intergenic SNP (rs243325; P = 9.91 × 10−9) and at the SPIB locus from a SPIB intronic SNP (rs34944112; P = 3.65 × 10−9). Among the associated SNPs at the CLEC16A-SOCS1 locus, two within the CLEC16A gene as well as one SOCS1 SNP (rs243325) remained significant after conditional logistic regression and contributed independently to risk. Sequencing of the SIAE gene and functional assays of newly identified variants revealed six patients with functional non-synonymous SIAE mutations (Fisher’s P = 9 × 10−4 vs controls) We demonstrate independent effects on risk of PBC for CLEC16A, SOCS1 and SPIB variants, while identifying functionally defective SIAE variants as potential factors in risk for PBC.

Susceptibility to primary biliary cirrhosis (PBC) is strongly associated with HLA region polymorphisms. To determine if associations can be explained by classical HLA determinants we studied Italian 676 cases and 1440 controls with genotyped with dense single nucleotide polymorphisms (SNPs) for which classical HLA alleles and amino acids were imputed. Although previous genome-wide association studies and our results show stronger SNP associations near DQB1, we demonstrate that the HLA signals can be attributed to classical DRB1 and DPB1 genes. Strong support for the predominant role of DRB1 is provided by our conditional analyses. We also demonstrate an independent association of DPB1. Specific HLA-DRB1 genes (*08, *11 and *14) account for most of the DRB1 association signal. Consistent with previous studies, DRB1*08 (p = 1.59 × 10−11) was the strongest predisposing allele where as DRB1*11 (p = 1.42 × 10−10) was protective. Additionally DRB1*14 and the DPB1 association (DPB1*03:01) (p = 9.18 × 10−7) were predisposing risk alleles. No signal was observed in the HLA class 1 or class 3 regions. These findings better define the association of PBC with HLA and specifically support the role of classical HLA-DRB1 and DPB1 genes and alleles in susceptibility to PBC.