The cluster variation method for efficient linkage analysis on extended pedigrees.

BACKGROUND: Computing exact multipoint LOD scores for extended pedigrees rapidly becomes infeasible as the number of markers and untyped individuals increase. When markers are excluded from the computation, significant power may be lost. Therefore accurate approximate methods which take into account all markers are desirable. METHODS: We present a novel method for efficient estimation of LOD scores on extended pedigrees. Our approach is based on the Cluster Variation Method, which deterministically estimates likelihoods by performing exact computations on tractable subsets of variables (clusters) of a Bayesian network. First a distribution over inheritances on the marker loci is approximated with the Cluster Variation Method. Then this distribution is used to estimate the LOD score for each location of the trait locus. RESULTS: First we demonstrate that significant power may be lost if markers are ignored in the multi-point analysis. On a set of pedigrees where exact computation is possible we compare the estimates of the LOD scores obtained with our method to the exact LOD scores. Secondly, we compare our method to a state of the art MCMC sampler. When both methods are given equal computation time, our method is more efficient. Finally, we show that CVM scales to large problem instances. CONCLUSION: We conclude that the Cluster Variation Method is as accurate as MCMC and generally is more efficient. Our method is a promising alternative to approaches based on MCMC sampling.

Diagnostic genome profiling in mental retardation.

Mental retardation (MR) occurs in 2%-3% of the general population. Conventional karyotyping has a resolution of 5-10 million bases and detects chromosomal alterations in approximately 5% of individuals with unexplained MR. The frequency of smaller submicroscopic chromosomal alterations in these patients is unknown. Novel molecular karyotyping methods, such as array-based comparative genomic hybridization (array CGH), can detect submicroscopic chromosome alterations at a resolution of 100 kb. In this study, 100 patients with unexplained MR were analyzed using array CGH for DNA copy-number changes by use of a novel tiling-resolution genomewide microarray containing 32,447 bacterial artificial clones. Alterations were validated by fluorescence in situ hybridization and/or multiplex ligation-dependent probe amplification, and parents were tested to determine de novo occurrence. Reproducible DNA copy-number changes were present in 97% of patients. The majority of these alterations were inherited from phenotypically normal parents, which reflects normal large-scale copy-number variation. In 10% of the patients, de novo alterations considered to be clinically relevant were found: seven deletions and three duplications. These alterations varied in size from 540 kb to 12 Mb and were scattered throughout the genome. Our results indicate that the diagnostic yield of this approach in the general population of patients with MR is at least twice as high as that of standard GTG-banded karyotyping.

Molecular characterisation of patients with subtelomeric 22q abnormalities using chromosome specific array-based comparative genomic hybridisation.

The 22q13 deletion syndrome is associated with global developmental delay, absent or delayed speech, and generalised hypotonia. In this study, the size and nature of 22q13 deletions (n=9) were studied in detail by high-resolution chromosome specific array-based comparative genomic hybridisation (array CGH). The deletion sizes varied considerably between the different patients, that is, the largest deletion spanning 8.4 Mb with the breakpoint mapping to 22q13.2 and the smallest deletion spanning 3.3 Mb with the breakpoint mapping to 22q13.31. In one case, a unique subtelomeric 3.9 Mb deletion associated with a 2.0 Mb duplication of 22q13 was observed, adding to a growing number of similar cases identified for other chromosome ends. Remarkably, this patient had signs suggestive of retinitis pigmentosa, which has never been reported before in the 22q13 deletion syndrome. The identification of two pairs of recurrent proximal breakpoints on 22q13 suggests that these specific regions may be prone to recombination, due to yet unknown genome architectural features. In addition to the copy number changes on 22q13, a duplication of approximately 330 kb on 22q11.1 was observed and shown to be a genetic large-scale copy number variation without clinical consequences. The current study failed to reveal relationships between the clinical features and the deletion sizes. Global developmental delay and absent or severely delayed speech were observed in all patients, whereas hypotonia was present in 89% of the cases (8/9). This study underscores the utility of array CGH for characterising the size and nature of subtelomeric deletions, such as monosomy 22q13, and underlines the considerable variability in deletion size in the 22q13 deletion syndrome regardless of the clinical phenotype.