Lost in Transition – The Island of Susak (1951–2001)

The isolated population of the Island of Susak was thoroughly studied by a multidisciplinary team of the Croatian Academy of Sciences and Arts in early 1950’s. Recently, a team of scientists revisited the island. This paper describes the main characteristics of the transition process during which a massive exodus occurred with 90% of the island’s population migrating to New Jersey, USA. We summarise the differences in lifestyle, economy, social structure and sense of identity between the historic (1950’s) and contemporary (2001) Susak population. We applied contemporary methods (analysis of microsatellite DNA polymorphisms) to investigate local myths about extreme levels of inbreeding and genetic homogeneity among the Susak islanders. Analysis of short-tandemrepeat (STR) loci showed that Susak displayed characteristics of a small homogeneous breeding isolate. The average heterozygosity was found to be low compared to outbred populations. The signature of a recent severe bottleneck could be detected. Analysis of 8 markers located on Xq13–21 in 71 individuals suggested extensive level of linkage disequilibrium (LD). A migrant study was designed to investigate the effects of large environmental changes (Susak vs. USA) and inbreeding (Susak vs. Croatian general population) on some biologically important quantitative traits, such as blood pressure and serum lipids. The results confirmed the positive correlation between inbreeding level and blood pressure that has been reported in the literature on several occasions. The last remnants of this traditional island community will soon be lost forever

Uncovering Networks from Genome-Wide Association Studies via Circular Genomic Permutation

Genome-wide association studies (GWAS) aim to detect single nucleotide polymorphisms (SNP) associated with trait variation. However, due to the large number of tests, standard analysis techniques impose highly stringent significance thresholds, leaving potentially associated SNPs undetected, and much of the trait genetic variation unexplained. Pathway- and network-based methodologies applied to GWAS aim to detect associations missed by standard single-marker approaches. The complex and non-random architecture of the genome makes it a challenge to derive an appropriate testing framework for such methodologies. We developed a rapid and simple permutation approach that uses GWAS SNP association results to establish the significance of pathway associations while accounting for the linkage disequilibrium structure of SNPs and the clustering of functionally related elements in the genome. All SNPs used in the GWAS are placed in a “circular genome” according to their location. Then the complete set of SNP association P values are permuted by rotation with respect to the genomic locations of the SNPs. Once these “simulated” P values are assigned, the joint gene P values are calculated using Fisher’s combination test, and the association of pathways is tested using the hypergeometric test. The circular genomic permutation approach was applied to a human genome-wide association dataset. The data consists of 719 individuals from the ORCADES study genotyped for ∼300,000 SNPs and measured for 51 traits ranging from physical to biochemical measurements. KEGG pathways (n = 225) were used as the sets of pathways to be tested. Our results demonstrate that the circular genomic permutations provide robust association P values. The non-permuted hypergeometric analysis generates ∼1400 pathway-trait combination results with an association P value more significant than P ≤ 0.05, whereas applying circular genomic permutation reduces the number of significant results to a more credible 40% of that value. The circular permutation software (“genomicper”) is available as an R package at http://cran.r-project.org/

Comparative assessment of methods for estimating individual genome-wide homozygosity-by-descent from human genomic data

<p>Abstract</p> <p>Background</p> <p>Genome-wide homozygosity estimation from genomic data is becoming an increasingly interesting research topic. The aim of this study was to compare different methods for estimating individual homozygosity-by-descent based on the information from human genome-wide scans rather than genealogies. We considered the four most commonly used methods and investigated their applicability to single-nucleotide polymorphism (SNP) data in both a simulation study and by using the human genotyped data. A total of 986 inhabitants from the isolated Island of Vis, Croatia (where inbreeding is present, but no pedigree-based inbreeding was observed at the level of F > 0.0625) were included in this study. All individuals were genotyped with the Illumina HumanHap300 array with 317,503 SNP markers.</p> <p>Results</p> <p>Simulation data suggested that multi-point FEstim is the method most strongly correlated to true homozygosity-by-descent. Correlation coefficients between the homozygosity-by-descent estimates were high but only for inbred individuals, with nearly absolute correlation between single-point measures.</p> <p>Conclusions</p> <p>Deciding who is really inbred is a methodological challenge where multi-point approaches can be very helpful once the set of SNP markers is filtered to remove linkage disequilibrium. The use of several different methodological approaches and hence different homozygosity measures can help to distinguish between homozygosity-by-state and homozygosity-by-descent in studies investigating the effects of genomic autozygosity on human health.</p

Enrichment of pathogenic alleles in the brittle cornea gene, ZNF469, in keratoconus

Keratoconus, a common inherited ocular disorder resulting in progressive corneal thinning, is the leading indication for corneal transplantation in the developed world. Genome-wide association studies have identified common SNPs 100 kb upstream of ZNF469 strongly associated with corneal thickness. Homozygous mutations in ZNF469 and PR domain-containing protein 5 (PRDM5) genes result in brittle cornea syndrome (BCS) Types 1 and 2, respectively. BCS is an autosomal recessive generalized connective tissue disorder associated with extreme corneal thinning and a high risk of corneal rupture. Some individuals with heterozygous PRDM5 mutations demonstrate a carrier ocular phenotype, which includes a mildly reduced corneal thickness, keratoconus and blue sclera. We hypothesized that heterozygous variants in PRDM5 and ZNF469 predispose to the development of isolated keratoconus. We found a significant enrichment of potentially pathologic heterozygous alleles in ZNF469 associated with the development of keratoconus (P = 0.00102) resulting in a relative risk of 12.0. This enrichment of rare potentially pathogenic alleles in ZNF469 in 12.5% of keratoconus patients represents a significant mutational load and highlights ZNF469 as the most significant genetic factor responsible for keratoconus identified to date

Enrichment of pathogenic alleles in the brittle cornea gene, ZNF469, in keratoconus

Keratoconus, a common inherited ocular disorder resulting in progressive corneal thinning, is the leading indication for corneal transplantation in the developed world. Genome-wide association studies have identified common SNPs 100 kb upstream of ZNF469 strongly associated with corneal thickness. Homozygous mutations in ZNF469 and PR domain-containing protein 5 (PRDM5) genes result in brittle cornea syndrome (BCS) Types 1 and 2, respectively. BCS is an autosomal recessive generalized connective tissue disorder associated with extreme corneal thinning and a high risk of corneal rupture. Some individuals with heterozygous PRDM5 mutations demonstrate a carrier ocular phenotype, which includes a mildly reduced corneal thickness, keratoconus and blue sclera. We hypothesized that heterozygous variants in PRDM5 and ZNF469 predispose to the development of isolated keratoconus. We found a significant enrichment of potentially pathologic heterozygous alleles in ZNF469 associated with the development of keratoconus (P = 0.00102) resulting in a relative risk of 12.0. This enrichment of rare potentially pathogenic alleles in ZNF469 in 12.5% of keratoconus patients represents a significant mutational load and highlights ZNF469 as the most significant genetic factor responsible for keratoconus identified to dat

Lost in Transition – The Island of Susak (1951–2001)

The isolated population of the Island of Susak was thoroughly studied by a multidisciplinary team of the Croatian Academy of Sciences and Arts in early 1950’s. Recently, a team of scientists revisited the island. This paper describes the main characteristics of the transition process during which a massive exodus occurred with 90% of the island’s population migrating to New Jersey, USA. We summarise the differences in lifestyle, economy, social structure and sense of identity between the historic (1950’s) and contemporary (2001) Susak population. We applied contemporary methods (analysis of microsatellite DNA polymorphisms) to investigate local myths about extreme levels of inbreeding and genetic homogeneity among the Susak islanders. Analysis of short-tandemrepeat (STR) loci showed that Susak displayed characteristics of a small homogeneous breeding isolate. The average heterozygosity was found to be low compared to outbred populations. The signature of a recent severe bottleneck could be detected. Analysis of 8 markers located on Xq13–21 in 71 individuals suggested extensive level of linkage disequilibrium (LD). A migrant study was designed to investigate the effects of large environmental changes (Susak vs. USA) and inbreeding (Susak vs. Croatian general population) on some biologically important quantitative traits, such as blood pressure and serum lipids. The results confirmed the positive correlation between inbreeding level and blood pressure that has been reported in the literature on several occasions. The last remnants of this traditional island community will soon be lost forever

The Eysenck personality factors: psychometric structure, reliability, heritability and phenotypic and genetic correlations with psychological distress in an isolated Croatian population

We report the psychometric structure of a Croatian translation of the Eysenck Personality Questionnaire-Revised (short-form), its correlations with psychological distress (General Health Questionnaire-30), its heritability, and personality–psychological distress genetic correlations. The setting is a large (≈1000), family-based sample of men and women from an isolated Croatian island. The neuroticism and extraversion traits and the lie scale showed good psychometric characteristics. The translated psychoticism scale was unsatisfactory in this sample. It had a very low internal consistency, probably due in part to heavily biased item responses. There were significant additive genetic contributions to variation in neuroticism, extraversion, and psychological distress. Psychological distress had a very high genetic correlation with neuroticism, and a moderate genetic correlation with extraversion

Bayesian methods for instrumental variable analysis with genetic instruments (‘Mendelian randomization’): example with urate transporter SLC2A9 as an instrumental variable for effect of urate levels on metabolic syndrome

The ‘Mendelian randomization’ approach uses genotype as an instrumental variable to distinguish between causal and non-causal explanations of biomarker–disease associations. Classical methods for instrumental variable analysis are limited to linear or probit models without latent variables or missing data, rely on asymptotic approximations that are not valid for weak instruments and focus on estimation rather than hypothesis testing. We describe a Bayesian approach that overcomes these limitations, using the JAGS program to compute the log-likelihood ratio (lod score) between causal and non-causal explanations of a biomarker–disease association. To demonstrate the approach, we examined the relationship of plasma urate levels to metabolic syndrome in the ORCADES study of a Scottish population isolate, using genotype at six single-nucleotide polymorphisms in the urate transporter gene SLC2A9 as an instrumental variable. In models that allow for intra-individual variability in urate levels, the lod score favouring a non-causal over a causal explanation was 2.34. In models that do not allow for intra-individual variability, the weight of evidence against a causal explanation was weaker (lod score 1.38). We demonstrate the ability to test one of the key assumptions of instrumental variable analysis—that the effects of the instrument on outcome are mediated only through the intermediate variable—by constructing a test for residual effects of genotype on outcome, similar to the tests of ‘overidentifying restrictions’ developed for classical instrumental variable analysis. The Bayesian approach described here is flexible enough to deal with any instrumental variable problem, and does not rely on asymptotic approximations that may not be valid for weak instruments. The approach can easily be extended to combine information from different study designs. Statistical power calculations show that instrumental variable analysis with genetic instruments will typically require combining information from moderately large cohort and cross-sectional studies of biomarkers with information from very large genetic case–control studies

Characterisation of Genome-Wide Association Epistasis Signals for Serum Uric Acid in Human Population Isolates

Genome-wide association (GWA) studies have identified a number of loci underlying variation in human serum uric acid (SUA) levels with the SLC2A9 gene having the largest effect identified so far. Gene-gene interactions (epistasis) are largely unexplored in these GWA studies. We performed a full pair-wise genome scan in the Italian MICROS population (n = 1201) to characterise epistasis signals in SUA levels. In the resultant epistasis profile, no SNP pairs reached the Bonferroni adjusted threshold for the pair-wise genome-wide significance. However, SLC2A9 was found interacting with multiple loci across the genome, with NFIA - SLC2A9 and SLC2A9 - ESRRAP2 being significant based on a threshold derived for interactions between GWA significant SNPs and the genome and jointly explaining 8.0% of the phenotypic variance in SUA levels (3.4% by interaction components). Epistasis signal replication in a CROATIAN population (n = 1772) was limited at the SNP level but improved dramatically at the gene ontology level. In addition, gene ontology terms enriched by the epistasis signals in each population support links between SUA levels and neurological disorders. We conclude that GWA epistasis analysis is useful despite relatively low power in small isolated populations