Search

Variants identified in earlier genome-wide association studies (GWAS) on differentiated thyroid carcinoma (DTC) explain about 10% of the overall estimated genetic contribution and could not provide complete insights into biological mechanisms involved in DTC susceptibility. Integrating systems biology information from model organisms, genome-wide expression data from tumor and matched normal tissue and GWAS data could help identifying DTC-associated genes, and pathways or functional networks in which they are involved. We performed data mining of GWAS data of the EPITHYR consortium (1551 cases and 1957 controls) using various pathways and protein-protein interaction (PPI) annotation databases and gene expression data from The Cancer Genome Atlas. We identified eight DTC-associated genes at known loci 2q35 (DIRC3), 8p12 (NRG1), 9q22 (FOXE1, TRMO, HEMGN, ANP32B, NANS) and 14q13 (MBIP). Using the EW&#95;dmGWAS approach we found that gene networks related to glycogenolysis, glycogen metabolism, insulin metabolism and signal transduction pathways associated with muscle contraction were overrepresented with association signals (false discovery rate adjusted p-value < 0.05). Additionally, suggestive association of 21 KEGG and 75 REACTOME pathways with DTC indicate a link between DTC susceptibility and functions related to metabolism of cholesterol, amino sugar and nucleotide sugar metabolism, steroid biosynthesis, and downregulation of ERBB2 signaling pathways. Together, our results provide novel insights into biological mechanisms contributing to DTC risk.

Single-nucleotide polymorphisms (SNPs) in over 180 loci have been associated with breast cancer (BC) through genome-wide association studies involving mostly unselected population-based case-control series. Some of them modify BC risk of women carrying a BRCA1 or BRCA2 (BRCA1/2) mutation and may also explain BC risk variability in BC-prone families with no BRCA1/2 mutation. Here, we assessed the contribution of SNPs of the iCOGS array in GENESIS consisting of BC cases with no BRCA1/2 mutation and a sister with BC, and population controls. Genotyping data were available for 1281 index cases, 731 sisters with BC, 457 unaffected sisters and 1272 controls. In addition to the standard SNP-level analysis using index cases and controls, we performed pedigree-based association tests to capture transmission information in the sibships. We also performed gene- and pathway-level analyses to maximize the power to detect associations with lower-frequency SNPs or those with modest effect sizes. While SNP-level analyses identified 18 loci, gene-level analyses identified 112 genes. Furthermore, 31 Kyoto Encyclopedia of Genes and Genomes and 7 Atlas of Cancer Signaling Network pathways were highlighted (false discovery rate of 5%). Using results from the "index case-control" analysis, we built pathway-derived polygenic risk scores (PRS) and assessed their performance in the population-based CECILE study and in a data set composed of GENESIS-affected sisters and CECILE controls. Although these PRS had poor predictive value in the general population, they performed better than a PRS built using our SNP-level findings, and we found that the joint effect of family history and PRS needs to be considered in risk prediction models., (© 2020 The Authors. International Journal of Cancer published by John Wiley & Sons Ltd on behalf of UICC.)

Genome-wide association studies (GWAS) explore the genetic causes of complex diseases. However, classical approaches ignore the biological context of the genetic variants and genes under study. To address this shortcoming, one can use biological networks, which model functional relationships, to search for functionally related susceptibility loci. Many such network methods exist, each arising from different mathematical frameworks, pre-processing steps, and assumptions about the network properties of the susceptibility mechanism. Unsurprisingly, this results in disparate solutions. To explore how to exploit these heterogeneous approaches, we selected six network methods and applied them to GENESIS, a nationwide French study on familial breast cancer. First, we verified that network methods recovered more interpretable results than a standard GWAS. We addressed the heterogeneity of their solutions by studying their overlap, computing what we called the consensus. The key gene in this consensus solution was COPS5, a gene related to multiple cancer hallmarks. Another issue we observed was that network methods were unstable, selecting very different genes on different subsamples of GENESIS. Therefore, we proposed a stable consensus solution formed by the 68 genes most consistently selected across multiple subsamples. This solution was also enriched in genes known to be associated with breast cancer susceptibility (BLM, CASP8, CASP10, DNAJC1, FGFR2, MRPS30, and SLC4A7, P-value = 3 × 10-4). The most connected gene was CUL3, a regulator of several genes linked to cancer progression. Lastly, we evaluated the biases of each method and the impact of their parameters on the outcome. In general, network methods preferred highly connected genes, even after random rewirings that stripped the connections of any biological meaning. In conclusion, we present the advantages of network-guided GWAS, characterize their shortcomings, and provide strategies to address them. To compute the consensus networks, implementations of all six methods are available at https://github.com/hclimente/gwas-tools., Competing Interests: The authors have declared that no competing interests exist.

Rod-cone dystrophy (RCD), also called retinitis pigmentosa, is characterized by rod followed by cone photoreceptor degeneration, leading to gradual visual loss. Mutations in over 65 genes have been associated with non-syndromic RCD explaining 60% to 70% of cases, with novel gene defects possibly accounting for the unsolved cases. Homozygosity mapping and whole-exome sequencing applied to a case of autosomal recessive non-syndromic RCD from a consanguineous union identified a homozygous variant in WDR34. Mutations in WDR34 have been previously associated with severe ciliopathy syndromes possibly associated with a retinal dystrophy. This is the first report of a homozygous mutation in WDR34 associated with non-syndromic RCD., (© 2020 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.)

Pathogenic variants in BRCA1 and BRCA2 only explain the underlying genetic cause of about 10% of hereditary breast and ovarian cancer families. Because of cost-effectiveness, multigene panel testing is often performed even if the clinical utility of testing most of the genes remains questionable. The purpose of our study was to assess the contribution of rare, deleterious-predicted variants in DNA repair genes in familial breast cancer (BC) in a well-characterized and homogeneous population. We analyzed 113 DNA repair genes selected from either an exome sequencing or a candidate gene approach in the GENESIS study, which includes familial BC cases with no BRCA1 or BRCA2 mutation and having a sister with BC (N = 1,207), and general population controls (N = 1,199). Sequencing data were filtered for rare loss-of-function variants (LoF) and likely deleterious missense variants (MV). We confirmed associations between LoF and MV in PALB2, ATM and CHEK2 and BC occurrence. We also identified for the first time associations between FANCI, MAST1, POLH and RTEL1 and BC susceptibility. Unlike other associated genes, carriers of an ATM LoF had a significantly higher risk of developing BC than carriers of an ATM MV (OR LoF = 17.4 vs. OR MV = 1.6; p Het = 0.002). Hence, our approach allowed us to specify BC relative risks associated with deleterious-predicted variants in PALB2, ATM and CHEK2 and to add MAST1, POLH, RTEL1 and FANCI to the list of DNA repair genes possibly involved in BC susceptibility. We also highlight that different types of variants within the same gene can lead to different risk estimates., (© 2018 The Authors. International Journal of Cancer published by John Wiley & Sons Ltd on behalf of UICC.)

Background: Genetic factors may influence an individual's sensitivity to ionising radiation and therefore modify his/her risk of developing papillary thyroid carcinoma (PTC). Previously, we reported that common single nucleotide polymorphisms (SNPs) within the DNA damage recognition gene ATM contribute to PTC risk in Belarusian children exposed to fallout from the Chernobyl power plant accident. Here we explored in the same population the contribution of a panel of DNA repair-related SNPs in genes acting downstream of ATM., Methods: The association of 141 SNPs located in 43 DNA repair genes was examined in 75 PTC cases and 254 controls from the Gomel region in Belarus. All subjects were younger than 15 years at the time of the Chernobyl accident. Conditional logistic regressions accounting for radiation dose were performed with PLINK using the additive allelic inheritance model, and a linkage disequilibrium (LD)-based Bonferroni correction was used for correction for multiple testing., Results: The intronic SNP rs2296675 in MGMT was associated with an increased PTC risk [per minor allele odds ratio (OR) 2.54 95% CI 1.50, 4.30, P per allele  = 0.0006, P corr.= 0.05], and gene-wide association testing highlighted a possible role for ERCC5 (P Gene  = 0.01) and PCNA (P Gene  = 0.05) in addition to MGMT (P Gene  = 0.008)., Conclusions: These findings indicate that several genes acting in distinct DNA repair mechanisms contribute to PTC risk. Further investigation is needed to decipher the functional properties of the methyltransferase encoded by MGMT and to understand how alteration of such functions may lead to the development of the most common type of thyroid cancer.

Primary cilia are sensory organelles present on most mammalian cells. The assembly and maintenance of primary cilia are facilitated by intraflagellar transport (IFT), a bidirectional protein trafficking along the cilium. Mutations in genes coding for IFT components have been associated with a group of diseases called ciliopathies. These genetic disorders can affect a variety of organs including the retina. Using whole exome sequencing in three families, we identified mutations in Intraflagellar Transport 172 Homolog [IFT172 (Chlamydomonas)] that underlie an isolated retinal degeneration and Bardet-Biedl syndrome. Extensive functional analyses of the identified mutations in cell culture, rat retina and in zebrafish demonstrated their hypomorphic or null nature. It has recently been reported that mutations in IFT172 cause a severe ciliopathy syndrome involving skeletal, renal, hepatic and retinal abnormalities (Jeune and Mainzer-Saldino syndromes). Here, we report for the first time that mutations in this gene can also lead to an isolated form of retinal degeneration. The functional data for the mutations can partially explain milder phenotypes; however, the involvement of modifying alleles in the IFT172-associated phenotypes cannot be excluded. These findings expand the spectrum of disease associated with mutations in IFT172 and suggest that mutations in genes originally reported to be associated with syndromic ciliopathies should also be considered in subjects with non-syndromic retinal dystrophy., (© The Author 2014. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Pulmonary veno-occlusive disease (PVOD) is a rare and devastating cause of pulmonary hypertension that is characterized histologically by widespread fibrous intimal proliferation of septal veins and preseptal venules and is frequently associated with pulmonary capillary dilatation and proliferation. PVOD is categorized into a separate pulmonary arterial hypertension-related group in the current classification of pulmonary hypertension. PVOD presents either sporadically or as familial cases with a seemingly recessive mode of transmission. Using whole-exome sequencing, we detected recessive mutations in EIF2AK4 (also called GCN2) that cosegregated with PVOD in all 13 families studied. We also found biallelic EIF2AK4 mutations in 5 of 20 histologically confirmed sporadic cases of PVOD. All mutations, either in a homozygous or compound-heterozygous state, disrupted the function of the gene. These findings point to EIF2AK4 as the major gene that is linked to PVOD development and contribute toward an understanding of the complex genetic architecture of pulmonary hypertension.

The high prevalence of type 2 diabetes and its uneven distribution among human populations is both a major public health concern and a puzzle in evolutionary biology. Why is this deleterious disease so common, while the associated genetic variants should be removed by natural selection? The 'thrifty genotype' hypothesis proposed that the causal genetic variants were advantageous and selected for during the majority of human evolution. It remains, however, unclear whether genetic data support this scenario. In this study, we characterized patterns of selection at 10 variants associated with type 2 diabetes, contrasting one herder and one farmer population from Central Asia. We aimed at identifying which alleles (risk or protective) are under selection, dating the timing of selective events, and investigating the effect of lifestyle on selective patterns. We did not find any evidence of selection on risk variants, as predicted by the thrifty genotype hypothesis. Instead, we identified clear signatures of selection on protective variants, in both populations, dating from the beginning of the Neolithic, which suggests that this major transition was accompanied by a selective advantage for non-thrifty variants. Combining our results with worldwide data further suggests that East Asia was particularly prone to such recent selection of protective haplotypes. As much effort has been devoted so far to searching for thrifty variants, we argue that more attention should be paid to the evolution of non-thrifty variants.

Mutations in the CRB1 gene are associated with variable phenotypes of severe retinal dystrophies, ranging from leber congenital amaurosis (LCA) to rod-cone dystrophy, also called retinitis pigmentosa (RP). Moreover, retinal dystrophies resulting from CRB1 mutations may be accompanied by specific fundus features: preservation of the para-arteriolar retinal pigment epithelium (PPRPE) and retinal telangiectasia with exudation (also referred to as Coats-like vasculopathy). In this publication, we report seven novel mutations and classify over 150 reported CRB1 sequence variants that were found in more that 240 patients. The data from previous reports were used to analyze a potential correlation between CRB1 variants and the clinical features of respective patients. This meta-analysis suggests that the differential phenotype of patients with CRB1 mutations is due to additional modifying factors rather than particular mutant allele combination., (© 2011 Wiley Periodicals, Inc.)

Purpose: To identify the genetic defect of a consanguineous Portuguese family with rod-cone dystrophy and varying degrees of decreased audition., Methods: A detailed ophthalmic and auditory examination was performed on a Portuguese patient with severe autosomal recessive rod-cone dystrophy. Known genetic defects were excluded by performing autosomal recessive retinitis pigmentosa (arRP) genotyping microarray analysis and by Sanger sequencing of the coding exons and flanking intronic regions of eyes shut homolog-drosophila (EYS) and chromosome 2 open reading frame 71 (C2orf71). Subsequently, genome-wide homozygosity mapping was performed in DNA samples from available family members using a 700K single nucleotide polymorphism (SNP) microarray. Candidate genes present in the significantly large homozygous regions were screened for mutations using Sanger sequencing., Results: The largest homozygous region (~11 Mb) in the affected family members was mapped to chromosome 9, which harbors deafness, autosomal recessive 31 (DFNB31; a gene previously associated with Usher syndrome). Mutation analysis of DFNB31 in the index patient identified a novel one-base-pair deletion (c.737delC), which is predicted to lead to a truncated protein (p.Pro246HisfsX13) and co-segregated with the disease in the family. Ophthalmic examination of the index patient and the affected siblings showed severe rod-cone dystrophy. Pure tone audiometry revealed a moderate hearing loss in the index patient, whereas the affected siblings were reported with more profound and early onset hearing impairment., Conclusions: We report a novel truncating mutation in DFNB31 associated with severe rod-cone dystrophy and varying degrees of hearing impairment in a consanguineous family of Portuguese origin. This is the second report of DFNB31 implication in Usher type 2.

HLA typing, demographic and immunological risk factors for nevirapine and efavirenz reactions were studied in a French cohort of HIV patients. Cases with isolated rash were significantly associated with HLA-DRB101 allele. No liver toxicity was observed and no association was detected with the percentage of CD4 T-cells. This study suggests that HLA-DRB101 allele plays an important role in susceptibility to cutaneous reactions associated with nevirapine and efavirenz in HIV patients.

Background: Stevens-Johnson syndrome (SJS) and its severe form, toxic epidermal necrolysis (TEN), are rare but life-threatening cutaneous adverse reactions to drugs, especially to allopurinol, carbamazepine, lamotrigine, phenobarbital, phenytoine, sulfamethoxazole, oxicam and nevirapine. Recently, a strong association between carbamazepine and allopurinol induced SJS or TEN has been described with respectively, HLA-B*1502 and HLA-B*5801 in a Han Chinese population from Taiwan and other Asian countries., Objective: The objective is to further investigate the relationship between SJS/TEN and HLA-B in a large number of patients in a European population., Methods: HLA-B genotyping was performed on 150 patients included in a European study (RegiSCAR) of SJS and TEN. We focused on patients related to 'high-risk' drugs including: 31 cases related to allopurinol, 28 to sulfamethoxazole, 19 to lamotrigine and 14 to oxicam., Results: Sixty-one percent of 31 allopurinol-induced SJS/TEN patients carried the HLA-B*5801 allele and the figure was 55% for 27 patients of European ancestry [odds ratio=80 (34-187)], (P<10(-6)) as previously observed in Han Chinese. For other drugs, two rare alleles showed a weaker association with SJS/TEN in a limited number of patients: B*38 for sulfamethoxazole or lamotrigine-related patients, and B*73 for oxicam., Conclusion: At variance with prior results in Asia, this study shows that even when HLA-B alleles behave as strong risk factors, as for allopurinol, they are neither sufficient nor necessary to explain the disease. Further investigations are necessary to delineate the exact role of the HLA region in SJS/TEN, and to look for other associations in other regions of the genome.

Stevens-Johnson syndrome (SJS) and toxic epidermal necrolysis (TEN) are rare but severe cutaneous adverse drug reactions, which can be caused by a certain number of specific drugs among which is carbamazepine, an antiepileptic agent. A very strong association of carbamazepine-induced SJS with HLA-B*1502 has recently been described in the Han Chinese population. Here in, we report preliminary results from a European study (RegiSCAR) of 12 carbamazepine-induced SJS/TEN cases (nine French and three German). Among these only four had a HLA-B*1502 allele. Remarkably, these four patients had an Asian ancestry, whereas the others did not as far as we have ascertained. This shows that although the HLA region may contain important genes for SJS, the HLA-B*1502 allele is not a universal marker for this disease and that ethnicity matters.

Whereas the human linkage map appears on limited evidence to be constant over populations, maps of linkage disequilibrium (LD) vary among populations that differ in gene history. The greatest difference is between populations of sub-Saharan origin and populations remotely derived from Africa after a major bottleneck that reduced their heterozygosity and altered their Malecot parameters, increasing the intercept M that reflects association in founders and decreasing the exponential decline epsilon. Variation among populations within this ethnic dichotomy is much smaller. These observations validate use of a cosmopolitan LD map based on a sizeable sample representing a large population reliably typed for markers at high density. Then an LD map for a region or isolate within an ethnic group may be created by fitting the sample LD to the cosmopolitan map, estimating Malecot parameters simultaneously. The cosmopolitan map scaled by epsilon recovers 95% of the information that a local map at the same density gives and therefore more than the information in a low-resolution local map. Relative to a Eurasian cosmopolitan map the scaling factors are estimated to be 0.82 for isolates of European descent, 1.53 for Yorubans, and 1.74 for African Americans. These observations are consistent with a common bottleneck (perhaps but not necessarily speciation) approximately 173,500 years ago, if the bottleneck associated with migration out of Africa was 100,000 years ago. Eurasian populations (especially isolates with numerous cases) are efficient for genome scans, and populations of recent African origin (such as African Americans) are efficient for identification of causal polymorphisms within a candidate sequence.

Linkage tests to localize oligogenes have been extended during the past year. Using simulated data and multiplex selection we find that several tests on affected sib pairs have comparable power and type I error. Three variants of SIBPAL2 are favoured when substantial numbers of normal sibs are included, but performance relative to the BETA benchmark degrades rapidly as normal sibs are depleted by selective sampling or typing. Neglect of this fact may explain recent failure of retrospective collaboration to confirm asthma candidates in the 5q cytokine region that are supported by other studies. A fully quantitative trait favours variance components under complete ascertainment and two options in SIBPAL2 under multiplex selection, with substantial gain in power from covariance analysis if the covariate is independent of the candidate locus. A dichotomy and liability threshold give virtually identical results in the SOLAR variance components program. Comparison with single-marker parametric analysis suggests that extension to multiple markers would be competitive with nonparametric methods in power, and superior in depth of genetic analysis. The simulated examples illustrate common problems encountered with linkage scans for oligogenes. They show that nonparametric methods provide no panacea for analytical problems posed by different phenotypes and methods of ascertainment.

Atopy and asthma are complex genetic diseases resulting from the interactions of a number of genetic and environmental factors. We had previously reported allelic association between the IL9 marker on chromosome 5q31-33 and atopy. In order to further investigate the role of susceptibility genes on 5q31-33 in the development of atopy and asthma we have studied 240 UK families comprising 131 families selected at random, 60 multiplex families with affected sib pairs, and 49 single proband nuclear families. Polymorphic markers on 5q31-33 were genotyped and both single and multipoint linkage analysis was undertaken using the BETA program. We have used both affection status and quantitative scores for atopy and asthma for phenotypic variables, combining data into scores for asthma and atopy. The strongest suggestion of linkage using multipoint analysis was centred around D5S410 with a maximum Lod of 1.946 at location 171.3 cM and a standard error of 3.3 for the asthma quantitative score. There was no evidence of linkage with atopy, the atopy quantitative score or total serum IgE.

Previous studies have suggested that atopy is linked to the beta chain of the high affinity IgE receptor (Fcepsilon R1-beta) on chromosome 11q13. Fcepsilon R1-beta polymorphisms, I181L, V183L, and E237G, are reported to be associated with asthma and atopy. The aim of this study was to investigate linkage to Fcepsilon R1-beta in a UK population and to assess the frequency of the polymorphisms and their association with asthma and atopy. A sample of 131 families was recruited at random with a sample of 109 families ascertained via an asthmatic proband. Each subject completed a written and video-assisted questionnaire and underwent bronchial challenge and skin prick testing. Serum total and specific IgE levels were measured. Quantitative scores were derived for asthma and atopy using principal component analysis. Four microsatellite markers were genotyped, including Fcepsilon R1-beta. The frequency of the I181L and V183L polymorphisms were determined by sequencing, and the E237G polymorphism was determined using the amplification refractory mutation system. We found no evidence for linkage to Fcepsilon R1-beta and only weak evidence for linkage to the less informative marker E237G. We found no examples of the I181L/V183L polymorphism in our population sample. Our study has failed to strengthen the evidence for a candidate gene on chromosome 11q13.

The central problem of complex inheritance is to map oligogenes for disease susceptibility, integrating linkage and association over samples that differ in several ways. Combination of evidence over multiple samples with 1,037 families supports loci contributing to asthma susceptibility in the cytokine region on 5q [maximum logarithm of odds (lod) = 2.61 near IL-4], but no evidence for atopy. The principal problems with retrospective collaboration on linkage appear to have been solved, providing far more information than a single study. A multipoint lod table evaluated at commonly agreed reference loci is required for both collaboration and metaanalysis, but variations in ascertainment, pedigree structure, phenotype definition, and marker selection are tolerated. These methods are invariant with statistical methods that increase the power of lods and are applicable to all diseases, motivating collaboration rather than competition. In contrast to linkage, positional cloning by allelic association has yet to be extended to multiple samples, a prerequisite for efficient combination with linkage and the greatest current challenge to genetic epidemiology.

On the causal hypothesis, most genetic determinants of disease are single-nucleotide polymorphisms (SNPs) that are likely to be selected as markers for positional cloning. On the proximity hypothesis, most disease determinants will not be included among markers but may be detected through linkage disequilibrium with other SNPs. In that event, allelic association among SNPs is an essential factor in positional cloning. Recent simulation based on monotonic population expansion suggests that useful association does not usually extend beyond 3 kb. This is contradicted by significant disequilibrium at much greater distances, with corresponding reduction in the number of SNPs required for a cost-effective genome scan. A plausible explanation is that cyclical expansions follow population bottlenecks that establish new disequilibria. Data on more than 1,000 locus pairs indicate that most disequilibria trace to the Neolithic, with no apparent difference between haplotypes that are random or selected through a major disease gene. Short duration may be characteristic of alleles contributing to disease susceptibility and haplotypes characteristic of particular ethnic groups. Alleles that are highly polymorphic in all ethnic groups may be older, neutral, or advantageous, in weak disequilibrium with nearby markers, and therefore less useful for positional cloning of disease genes. Significant disequilibrium at large distance makes the number of suitably chosen SNPs required for genome screening as small as 30,000, or 1 per 100 kb, with greater density (including less common SNPs) reserved for candidate regions.

Combination of evidence over samples, each of which is too small to be conclusive, is the central problem in complex inheritance. There are three approaches: meta-analysis, prospective collaboration, and retrospective collaboration. Our experience with the first and last, which are the most cost-effective, is discussed.

An efficient test of deviation from Hardy-Weinberg frequencies with one degree of freedom was applied to 44 marker loci in a genome scan, and 7 loci had a significant excess of apparent homozygotes (chi2 (1) > 6) suggestive of typing error. In this example evidence for linkage did not increase when outliers were censored. Statistical quality control is an essential part of genotyping, and the effect of mistyping and map error should be considered in evaluating any genome scan.

Allelic association has proven useful to refine the location of major genes prior to positional cloning, but it is of uncertain value for genome scans in complex inheritance. We have extended kinship theory to give information content for linkage and allelic association. Application to pairs of closely linked markers as a surrogate for marker x oligogene pairs indicates that association is largely determined by regional founders, with little effect of subsequent demography. Sub-Saharan Africa has the least allelic association, consistent with settlement of other regions by small numbers of founders. Recent speculation about substantial advantages of isolates over large populations, of constant size over expansion, and of F1 hybrids over incrosses is not supported by theory or data. On the contrary, fewer affected cases, less opportunity for replication, and more stochastic variation tend to make isolates less informative for allelic association, as they are for linkage.

Recently a graphical study of linkage disequilibrium around the CAPN3 locus failed to refine the 1.3-Mb interval suggested by haplotype sharing. On the contrary, the Malecot model as implemented in the ALLASS program maps CAPN3 within 3 kb of its true location (23 kb from the locus midpoint), overcoming identified problems with small samples, interrelated sibships, and short duration.

Recombination acts on the genetic map, not on the physical map. On the other hand, the physical map is usually more accurate. Choice of the genetic or physical map for positional cloning by allelic association depends on the goodness of fit of data to each map under an established model. Huntington disease illustrates the usual case in which the greater reliability of physical data outweighs recombinational heterogeneity. Hemochromatosis represents an exceptional case in which unrecognized recombinational heterogeneity retarded positional cloning for a decade. The Malecot model performs well for major genes, but no approach assuming either equilibrium or disequilibrium has been validated for oligogenes contributing to common disease. In this case of greatest interest, the power of allelic association relative to linkage is less clear than for major genes.

We have undertaken a study of the haplotypes among French potential bone marrow donors in order to define the geographical regions of France with the maximum of polymorphism and also to develop a strategy for optimal donor recruitment. A maximum likelihood estimator was used to calculate haplotype frequencies and their support limits for each region and for the whole of France. The observed differences between the regions were statistically significant. For each region, the minimum number of haplotypes necessary to explain 50% of the total frequency was calculated and compared with the equivalent values, and confidence intervals, obtained by repeated random samplings from the overall file. This approach shows that some regions (e.g., Provence) appear to be richer in terms of the numbers of haplotypes observed, and others (e.g., Bretagne) poorer. In the latter case, however, the frequencies of the most common haplotypes are greater. The haplotype frequencies of the whole sample were used to calculate the probability of finding a match for the next potential recipient for given sizes of the donor file, assuming random selection of donors. They were also used to calculate expected numbers of the major phenotypes, assuming Hardy-Weinberg equilibrium, and these were compared with those observed in the real data file. In this way, a large number of under-represented and nonrepresented phenotypes were identified. For each of these phenotypes, the most probable haplotypes and the regions in which these have the greatest frequencies have been identified. A search for donors with such particular phenotypes would be much more fruitful if directed towards these regions.

A likelihood approach to HLA serology has been developed in which the aim is not to define a recognition set for a serum but to describe the serum's ability to react with each and every antigen in the test cells, this ability being quantified in terms of the probability of a positive reaction. For a given set of probabilities, one for each antigen, it is possible to derive the probability of the observed set of reactions (the likelihood of the set of probabilities). The maximum possible value of the likelihood for any possible combination of the probability set can then be sought, but this requires a maximization of likelihood with respect to 60-100 independent parameters. Theoretical considerations of the shape of the likelihood surface prove that, in this particular case, this is a feasible proposition. This approach allows the recognition of three groups of antigens: those for which there is considerable evidence of a specificity, those for which there is either no specificity or a very weak specificity, and those for which there is insufficient evidence on which to base a conclusion. The existence of a specificity can be tested using a log likelihood ratio as a statistic, but the usual assumption of a chi 2 distribution of this statistic cannot automatically be made in this situation. Therefore, the distribution is estimated by simulation. A serologist using this approach would receive considerably more information as to the serum's reaction patterns and valid statistics for the existence, or not, of a specificity.