Augustine Kong

Large recurrent microdeletions associated with schizophrenia.

Reduced fecundity, associated with severe mental disorders, places negative selection pressure on risk alleles and may explain, in part, why common variants have not been found that confer risk of disorders such as autism, schizophrenia and mental retardation. Thus, rare variants may account for a larger fraction of the overall genetic risk than previously assumed. In contrast to rare single nucleotide mutations, rare copy number variations (CNVs) can be detected using genome-wide single nucleotide polymorphism arrays. This has led to the identification of CNVs associated with mental retardation and autism. In a genome-wide search for CNVs associating with schizophrenia, we used a population-based sample to identify de novo CNVs by analysing 9,878 transmissions from parents to offspring. The 66 de novo CNVs identified were tested for association in a sample of 1,433 schizophrenia cases and 33,250 controls. Three deletions at 1q21.1, 15q11.2 and 15q13.3 showing nominal association with schizophrenia in the first sample (phase I) were followed up in a second sample of 3,285 cases and 7,951 controls (phase II). All three deletions significantly associate with schizophrenia and related psychoses in the combined sample. The identification of these rare, recurrent risk variants, having occurred independently in multiple founders and being subject to negative selection, is important in itself. CNV analysis may also point the way to the identification of additional and more prevalent risk variants in genes and pathways involved in schizophrenia.

Allele-Sharing Models: LOD Scores and Accurate Linkage Tests

Starting with a test statistic for linkage analysis based on allele sharing, we propose an associated one-parameter model. Under general missing-data patterns, this model allows exact calculation of likelihood ratios and LOD scores and has been implemented by a simple modification of existing software. Most important, accurate linkage tests can be performed. Using an example, we show that some previously suggested approaches to handling less than perfectly informative data can be unacceptably conservative. Situations in which this model may not perform well are discussed, and an alternative model that requires additional computations is suggested.

Sequential Imputations and Bayesian Missing Data Problems

Abstract For missing data problems, Tanner and Wong have described a data augmentation procedure that approximates the actual posterior distribution of the parameter vector by a mixture of complete data posteriors. Their method of constructing the complete data sets is closely related to the Gibbs sampler. Both required iterations, and, similar to the EM algorithm, convergence can be slow. We introduce in this article an alternative procedure that involves imputing the missing data sequentially and computing appropriate importance sampling weights. In many applications this new procedure works very well without the need for iterations. Sensitivity analysis, influence analysis, and updating with new data can be performed cheaply. Bayesian prediction and model selection can also be incorporated. Examples taken from a wide range of applications are used for illustration.

Refining the impact of TCF7L2 gene variants on type 2 diabetes and adaptive evolution

We recently described an association between risk of type 2diabetes and variants in the transcription factor 7-like 2 gene (TCF7L2; formerly TCF4), with a population attributable risk (PAR) of 17%–28% in three populations of European ancestry. Here, we refine the definition of the TCF7L2 type 2diabetes risk variant, HapBT2D, to the ancestral T allele of a SNP, rs7903146, through replication in West African and Danish type 2 diabetes case-control studies and an expanded Icelandic study. We also identify another variant of the same gene, HapA, that shows evidence of positive selection in East Asian, European and West African populations. Notably, HapA shows a suggestive association with body mass index and altered concentrations of the hunger-satiety hormones ghrelin and leptin in males, indicating that the selective advantage of HapA may have been mediated through effects on energy metabolism.

Loci on chromosomes 2 (NIDDM1) and 15 interact to increase susceptibility to diabetes in Mexican Americans

Complex disorders such as diabetes, cardiovascular disease, asthma, hypertension and psychiatric illnesses account for a large and disproportionate share of health care costs, but remain poorly characterized with respect to aetiology. The transmission of such disorders is complex, reflecting the actions and interactions of multiple genetic and environmental factors. Genetic analyses that allow for the simultaneous consideration of susceptibility from multiple regions may improve the ability to map genes for complex disorders, but such analyses are currently computationally intensive and narrowly focused. We describe here an approach to assessing the evidence for statistical interactions between unlinked regions that allows multipoint allele–sharing analysis to take the evidence for linkage at one region into account in assessing the evidence for linkage over the rest of the genome. Using this method, we show that the interaction of genes on chromosomes 2 (NIDDM1) and 15 (near CYP19) makes a contribution to susceptibility to type 2 diabetes in Mexican Americans from Starr County, Texas.

Identification of low-frequency and rare sequence variants associated with elevated or reduced risk of type 2 diabetes.

Through whole-genome sequencing of 2,630 Icelanders and imputation into 11,114 Icelandic cases and 267,140 controls followed by testing in Danish and Iranian samples, we discovered 4 previously unreported variants affecting risk of type 2 diabetes (T2D). A low-frequency (1.47%) variant in intron 1 of CCND2, rs76895963[G], reduces risk of T2D by half (odds ratio (OR) = 0.53, P = 5.0 × 10−21) and is correlated with increased CCND2 expression. Notably, this variant is also associated with both greater height and higher body mass index (1.17 cm per allele, P = 5.5 × 10−12 and 0.56 kg/m2 per allele, P = 6.5 × 10−7, respectively). In addition, two missense variants in PAM, encoding p.Asp563Gly (frequency of 4.98%) and p.Ser539Trp (frequency of 0.65%), confer moderately higher risk of T2D (OR = 1.23, P = 3.9 × 10−10 and OR = 1.47, P = 1.7 × 10−5, respectively), and a rare (0.20%) frameshift variant in PDX1, encoding p.Gly218Alafs*12, associates with high risk of T2D (OR = 2.27, P = 7.3 × 10−7).

Graphical templates for model registration

A new method of model registration is proposed using graphical templates. A graph of landmarks is chosen in the template image. All possible candidates for these landmarks are found in the data image using local operators. A dynamic programming algorithm on decomposable subgraphs of the template graph finds the optimal match to a subset of the candidate points in polynomial time. This combination of local operators to describe points of interest/landmarks and a graph to describe their geometric orientation in the plane, yields fast and precise matches of the model to the data, with no initialization required.

A genome-wide scan for preeclampsia in the Netherlands

Preeclampsia, hallmarked by de novo hypertension and proteinuria in pregnancy, has a familial tendency. Recently, a large Icelandic genome-wide scan provided evidence for a maternal susceptibility locus for preeclampsia on chromosome 2p13 which was confirmed by a genome scan from Australia and New Zealand (NZ). The current study reports on a genome-wide scan of Dutch affected sib-pair families. In total 67 Dutch affected sib-pair families, comprising at least two siblings with proteinuric preeclampsia, eclampsia or HELLP-syndrome, were typed for 293 polymorphic markers throughout the genome and linkage analysis was performed. The highest allele sharing lod score of 1.99 was seen on chromosome 12q at 109.5 cM. Two peaks overlapped in the same regions between the Dutch and Icelandic genome-wide scan at chromosome 3p and chromosome 15q. No overlap was seen on 2p. Re-analysis in 38 families without HELLP-syndrome (preeclampsia families) and 34 families with at least one sibling with HELLP syndrome (HELLP families), revealed two peaks with suggestive evidence for linkage in the non-HELLP families on chromosome 10q (lod score 2.38, D10S1432, 93.9 cM) and 22q (lod score 2.41, D22S685, 32.4 cM). The peak on 12q appeared to be associated with HELLP syndrome; it increased to a lod score of 2.1 in the HELLP families and almost disappeared in the preeclampsia families. A nominal peak on chromosome 11 in the preeclampsia families showed overlap with the second highest peak in the Australian/NZ study. Results from our Dutch genome-wide scan indicate that HELLP syndrome might have a different genetic background than preeclampsia.

Blocking Gibbs sampling in very large probabilistic expert systems

Abstract We introduce a methodology for performing approximate computations in very complex probabilistic systems (e.g. huge pedigrees). Our approach, called blocking Gibbs, combines exact local computations with Gibbs sampling in a way that complements the strengths of both. The methodology is illustrated on a real-world problem involving a heavily inbred pedigreee containing 20 000 individuals. We present results showing that blocking-Gibbs sampling converges much faster than plain Gibbs sampling for very complex problems.

A rare nonsynonymous sequence variant in C3 is associated with high risk of age-related macular degeneration

Through whole-genome sequencing of 2,230 Icelanders, we detected a rare nonsynonymous SNP (minor allele frequency = 0.55%) in the C3 gene encoding a p.Lys155Gln substitution in complement factor 3, which, following imputation into a set of Icelandic cases with age-related macular degeneration (AMD) and controls, associated with disease (odds ratio (OR) = 3.45; P = 1.1 × 10−7). This signal is independent of the previously reported common SNPs in C3 encoding p.Pro314Leu and p.Arg102Gly that associate with AMD. The association of p.Lys155Gln was replicated in AMD case-control samples of European ancestry with OR = 4.22 and P = 1.6 × 10−10, resulting in OR = 3.65 and P = 8.8 × 10−16 for all studies combined. In vitro studies have suggested that the p.Lys155Gln substitution reduces C3b binding to complement factor H, potentially creating resistance to inhibition by this factor. This resistance to inhibition in turn is predicted to result in enhanced complement activation.