Partitioning the Heritability of Tourette Syndrome and Obsessive Compulsive Disorder Reveals Differences in Genetic Architecture

The direct estimation of heritability from genome-wide common variant data as implemented in the program Genome-wide Complex Trait Analysis (GCTA) has provided a means to quantify heritability attributable to all interrogated variants. We have quantified the variance in liability to disease explained by all SNPs for two phenotypically-related neurobehavioral disorders, obsessive-compulsive disorder (OCD) and Tourette Syndrome (TS), using GCTA. Our analysis yielded a heritability point estimate of 0.58 (se = 0.09, p = 5.64e-12) for TS, and 0.37 (se = 0.07, p = 1.5e-07) for OCD. In addition, we conducted multiple genomic partitioning analyses to identify genomic elements that concentrate this heritability. We examined genomic architectures of TS and OCD by chromosome, MAF bin, and functional annotations. In addition, we assessed heritability for early onset and adult onset OCD. Among other notable results, we found that SNPs with a minor allele frequency of less than 5% accounted for 21% of the TS heritability and 0% of the OCD heritability. Additionally, we identified a significant contribution to TS and OCD heritability by variants significantly associated with gene expression in two regions of the brain (parietal cortex and cerebellum) for which we had available expression quantitative trait loci (eQTLs). Finally we analyzed the genetic correlation between TS and OCD, revealing a genetic correlation of 0.41 (se = 0.15, p = 0.002). These results are very close to previous heritability estimates for TS and OCD based on twin and family studies, suggesting that very little, if any, heritability is truly missing (i.e., unassayed) from TS and OCD GWAS studies of common variation. The results also indicate that there is some genetic overlap between these two phenotypically-related neuropsychiatric disorders, but suggest that the two disorders have distinct genetic architectures.     

Assumptions and Properties of Limiting Pathway Models for Analysis of Epistasis in Complex Traits

For most complex traits, results from genome-wide association studies show that the proportion of the phenotypic variance attributable to the additive effects of individual SNPs, that is, the heritability explained by the SNPs, is substantially less than the estimate of heritability obtained by standard methods using correlations between relatives. This difference has been called the “missing heritability”. One explanation is that heritability estimates from family (including twin) studies are biased upwards. Zuk et al. revisited overestimation of narrow sense heritability from twin studies as a result of confounding with non-additive genetic variance. They propose a limiting pathway (LP) model that generates significant epistatic variation and its simple parametrization provides a convenient way to explore implications of epistasis. They conclude that over-estimation of narrow sense heritability from family data (‘phantom heritability’) may explain an important proportion of missing heritability. We show that for highly heritable quantitative traits large phantom heritability estimates from twin studies are possible only if a large contribution of common environment is assumed. The LP model is underpinned by strong assumptions that are unlikely to hold, including that all contributing pathways have the same mean and variance and are uncorrelated. Here, we relax the assumptions that underlie the LP model to be more biologically plausible. Together with theoretical, empirical, and pragmatic arguments we conclude that in outbred populations the contribution of additive genetic variance is likely to be much more important than the contribution of non-additive variance.     

Global Genetic Variations Predict Brain Response to Faces

Face expressions are a rich source of social signals. Here we estimated the proportion of phenotypic variance in the brain response to facial expressions explained by common genetic variance captured by ∼500,000 single nucleotide polymorphisms. Using genomic-relationship-matrix restricted maximum likelihood (GREML), we related this global genetic variance to that in the brain response to facial expressions, as assessed with functional magnetic resonance imaging (fMRI) in a community-based sample of adolescents (n = 1,620). Brain response to facial expressions was measured in 25 regions constituting a face network, as defined previously. In 9 out of these 25 regions, common genetic variance explained a significant proportion of phenotypic variance (40–50%) in their response to ambiguous facial expressions; this was not the case for angry facial expressions. Across the network, the strength of the genotype-phenotype relationship varied as a function of the inter-individual variability in the number of functional connections possessed by a given region (R² = 0.38, p<0.001). Furthermore, this variability showed an inverted U relationship with both the number of observed connections (R² = 0.48, p<0.001) and the magnitude of brain response (R² = 0.32, p<0.001). Thus, a significant proportion of the brain response to facial expressions is predicted by common genetic variance in a subset of regions constituting the face network. These regions show the highest inter-individual variability in the number of connections with other network nodes, suggesting that the genetic model captures variations across the adolescent brains in co-opting these regions into the face network.     

Genome Wide Meta-analysis Highlights the Role of Genetic Variation in RARRES2 in the Regulation of Circulating Serum Chemerin

Chemerin is an adipokine proposed to link obesity and chronic inflammation of adipose tissue. Genetic factors determining chemerin release from adipose tissue are yet unknown. We conducted a meta-analysis of genome-wide association studies (GWAS) for serum chemerin in three independent cohorts from Europe: Sorbs and KORA from Germany and PPP-Botnia from Finland (total N = 2,791). In addition, we measured mRNA expression of genes within the associated loci in peripheral mononuclear cells by micro-arrays, and within adipose tissue by quantitative RT-PCR and performed mRNA expression quantitative trait and expression-chemerin association studies to functionally substantiate our loci. Heritability estimate of circulating chemerin levels was 16.2% in the Sorbs cohort. Thirty single nucleotide polymorphisms (SNPs) at chromosome 7 within the retinoic acid receptor responder 2 (RARRES2)/Leucine Rich Repeat Containing (LRRC61) locus reached genome-wide significance (p<5.0×10⁻⁸) in the meta-analysis (the strongest evidence for association at rs7806429 with p = 7.8×10⁻¹⁴, beta = −0.067, explained variance 2.0%). All other SNPs within the cluster were in linkage disequilibrium with rs7806429 (minimum r² = 0.43 in the Sorbs cohort). The results of the subgroup analyses of males and females were consistent with the results found in the total cohort. No significant SNP-sex interaction was observed. rs7806429 was associated with mRNA expression of RARRES2 in visceral adipose tissue in women (p<0.05 after adjusting for age and body mass index). In conclusion, the present meta-GWAS combined with mRNA expression studies highlights the role of genetic variation in the RARRES2 locus in the regulation of circulating chemerin concentrations.     

A Variant In the Abo Gene Explains the Variation in Soluble E-Selectin Levels—Results from Dense Genotyping in Two Independent Populations

Background

Elevated soluble (s) E-selectin levels have been associated with various cardiovascular diseases. Recently, genetic variants in the ABO blood group have been related to E-selectin levels in a small cohort of patients with type 1 diabetes. We evaluated whether this association is reproducible in two large samples of Caucasians.

Methodology/ Principal Findings

Data of the present study was drawn from the population-based MONICA/KORA Augsburg study (n = 1,482) and the patients-based LURIC study (n = 1,546). A high-density genotyping array (50K IBC Chip) containing single-nucleotide polymorphisms (SNPs) from E-selectin candidate genes selected on known biology of E-selectin metabolism, mouse genetic studies, and human genetic association studies, was used for genotyping. Linear regression analyses with adjustment for age and sex (and survey in KORA) were applied to assess associations between gene variants and sE-selectin concentrations. A number of 12 SNPs (in KORA) and 13 SNPs (in LURIC), all from the ABO blood group gene, were significantly associated with the log-transformed concentration of E-selectin. The strongest association was observed for rs651007 with a change of log-transformed sE-selectin per one copy of the minor allele of −0.37 ng/ml (p = 1.87×10⁻¹⁰³) in KORA and −0.35 ng/ml (p = 5.11×10⁻⁸⁴) in LURIC. Inclusion of rs651007 increased the explained sE-selectin variance by 0.256 in KORA and 0.213 in LURIC. All SNPs had minor allele frequencies above 20% showing a substantial gene variation.

Conclusions/ Significance

Our findings in two independent samples indicate that the genetic variants at the ABO locus affect sE-selectin levels. Since distinct genome-wide association studies linked the ABO gene with myocardial infarction (MI) in the presence of coronary atherosclerosis and with coronary artery disease, these findings may not only enhance our understanding of adhesion molecule biology, but may also provide a focus for several novel research avenues.     

The Contribution of Diet and Genotype to Iron Status in Women: A Classical Twin Study

This is the first published report examining the combined effect of diet and genotype on body iron content using a classical twin study design. The aim of this study was to determine the relative contribution of genetic and environmental factors in determining iron status. The population was comprised of 200 BMI- and age-matched pairs of MZ and DZ healthy twins, characterised for habitual diet and 15 iron-related candidate genetic markers. Variance components analysis demonstrated that the heritability of serum ferritin (SF) and soluble transferrin receptor was 44% and 54% respectively. Measured single nucleotide polymorphisms explained 5% and selected dietary factors 6% of the variance in iron status; there was a negative association between calcium intake and body iron (p = 0.02) and SF (p = 0.04).     

Genetics of Microstructure of the Corpus Callosum in Older Adults

The current study sought to examine the relative influence of genetic and environmental factors on corpus callosum (CC) microstructure in a community sample of older adult twins. Analyses were undertaken in 284 healthy older twins (66% female; 79 MZ and 63 DZ pairs) from the Older Australian Twins Study. The average age of the sample was 69.82 (SD = 4.76) years. Brain imaging scans were collected and DTI measures were estimated for the whole CC as well as its five subregions. Parcellation of the CC was performed using Analyze. In addition, white matter lesion (WMLs) burden was estimated. Heritability and genetic correlation analyses were undertaken using the SOLAR software package. Age, sex, scanner, handedness and blood pressure were considered as covariates. Heritability (h²) analysis for the DTI metrics of whole CC, indicated significant h² for fractional anisotropy (FA) (h² = 0.56; p = 2.89×10⁻¹⁰), mean diffusivity (MD) (h² = 0.52; p = 0.30×10⁻⁶), radial diffusivity (RD) (h² = 0.49; p = 0.2×10⁻⁶) and axial diffusivity (AD) (h² = 0.37; p = 8.15×10⁻⁵). We also performed bivariate genetic correlation analyses between (i) whole CC DTI measures and (ii) whole CC DTI measures with total brain WML burden. Across the DTI measures for the whole CC, MD and RD shared 84% of the common genetic variance, followed by MD- AD (77%), FA - RD (52%), RD - AD (37%) and FA – MD (11%). For total WMLs, significant genetic correlations indicated that there was 19% shared common genetic variance with whole CC MD, followed by CC RD (17%), CC AD (16%) and CC FA (5%). Our findings suggest that the CC microstructure is under moderate genetic control. There was also evidence of shared genetic factors between the CC DTI measures. In contrast, there was less shared genetic variance between WMLs and the CC DTI metrics, suggesting fewer common genetic variants.     

Genetic Variation in the Human Brain Dopamine System Influences Motor Learning and Its Modulation by L-Dopa

Dopamine is important to learning and plasticity. Dopaminergic drugs are the focus of many therapies targeting the motor system, where high inter-individual differences in response are common. The current study examined the hypothesis that genetic variation in the dopamine system is associated with significant differences in motor learning, brain plasticity, and the effects of the dopamine precursor L-Dopa. Skilled motor learning and motor cortex plasticity were assessed using a randomized, double-blind, placebo-controlled, crossover design in 50 healthy adults during two study weeks, one with placebo and one with L-Dopa. The influence of five polymorphisms with established effects on dopamine neurotransmission was summed using a gene score, with higher scores corresponding to higher dopaminergic neurotransmission. Secondary hypotheses examined each polymorphism individually. While training on placebo, higher gene scores were associated with greater motor learning (p = .03). The effect of L-Dopa on learning varied with the gene score (gene score*drug interaction, p = .008): participants with lower gene scores, and thus lower endogenous dopaminergic neurotransmission, showed the largest learning improvement with L-Dopa relative to placebo (p<.0001), while L-Dopa had a detrimental effect in participants with higher gene scores (p = .01). Motor cortex plasticity, assessed via transcranial magnetic stimulation (TMS), also showed a gene score*drug interaction (p = .02). Individually, DRD2/ANKK1 genotype was significantly associated with motor learning (p = .02) and its modulation by L-Dopa (p<.0001), but not with any TMS measures. However, none of the individual polymorphisms explained the full constellation of findings associated with the gene score. These results suggest that genetic variation in the dopamine system influences learning and its modulation by L-Dopa. A polygene score explains differences in L-Dopa effects on learning and plasticity most robustly, thus identifying distinct biological phenotypes with respect to L-Dopa effects on learning and plasticity. These findings may have clinical applications in post-stroke rehabilitation or the treatment of Parkinson''s disease.     

Alzheimer’s Disease: Analyzing the Missing Heritability

Alzheimer’s disease (AD) is a complex disorder influenced by environmental and genetic factors. Recent work has identified 11 AD markers in 10 loci. We used Genome-wide Complex Trait Analysis to analyze >2 million SNPs for 10,922 individuals from the Alzheimer’s Disease Genetics Consortium to assess the phenotypic variance explained first by known late-onset AD loci, and then by all SNPs in the Alzheimer’s Disease Genetics Consortium dataset. In all, 33% of total phenotypic variance is explained by all common SNPs. APOE alone explained 6% and other known markers 2%, meaning more than 25% of phenotypic variance remains unexplained by known markers, but is tagged by common SNPs included on genotyping arrays or imputed with HapMap genotypes. Novel AD markers that explain large amounts of phenotypic variance are likely to be rare and unidentifiable using genome-wide association studies. Based on our findings and the current direction of human genetics research, we suggest specific study designs for future studies to identify the remaining heritability of Alzheimer’s disease.     

Ubiquitous Polygenicity of Human Complex Traits: Genome-Wide Analysis of 49 Traits in Koreans

Recent studies in population of European ancestry have shown that 30%∼50% of heritability for human complex traits such as height and body mass index, and common diseases such as schizophrenia and rheumatoid arthritis, can be captured by common SNPs and that genetic variation attributed to chromosomes are in proportion to their length. Using genome-wide estimation and partitioning approaches, we analysed 49 human quantitative traits, many of which are relevant to human diseases, in 7,170 unrelated Korean individuals genotyped on 326,262 SNPs. For 43 of the 49 traits, we estimated a nominally significant (P<0.05) proportion of variance explained by all SNPs on the Affymetrix 5.0 genotyping array (). On average across 47 of the 49 traits for which the estimate of is non-zero, common SNPs explain approximately one-third (range of 7.8% to 76.8%) of narrow sense heritability.The estimate of is highly correlated with the proportion of SNPs with association P<0.031 (r ² = 0.92). Longer genomic segments tend to explain more phenotypic variation, with a correlation of 0.78 between the estimate of variance explained by individual chromosomes and their physical length, and 1% of the genome explains approximately 1% of the genetic variance. Despite the fact that there are a few SNPs with large effects for some traits, these results suggest that polygenicity is ubiquitous for most human complex traits and that a substantial proportion of the “missing heritability” is captured by common SNPs.     

Genetic Variation in DNA Repair Pathways and Risk of Non-Hodgkin's Lymphoma

Molecular and genetic evidence suggests that DNA repair pathways may contribute to lymphoma susceptibility. Several studies have examined the association of DNA repair genes with lymphoma risk, but the findings from these reports have been inconsistent. Here we provide the results of a focused analysis of genetic variation in DNA repair genes and their association with the risk of non-Hodgkin''s lymphoma (NHL). With a population of 1,297 NHL cases and 1,946 controls, we have performed a two-stage case/control association analysis of 446 single nucleotide polymorphisms (SNPs) tagging the genetic variation in 81 DNA repair genes. We found the most significant association with NHL risk in the ATM locus for rs227060 (OR = 1.27, 95% CI: 1.13–1.43, p = 6.77×10⁻⁵), which remained significant after adjustment for multiple testing. In a subtype-specific analysis, associations were also observed for the ATM locus among both diffuse large B-cell lymphomas (DLBCL) and small lymphocytic lymphomas (SLL), however there was no association observed among follicular lymphomas (FL). In addition, our study provides suggestive evidence of an interaction between SNPs in MRE11A and NBS1 associated with NHL risk (OR = 0.51, 95% CI: 0.34–0.77, p = 0.0002). Finally, an imputation analysis using the 1,000 Genomes Project data combined with a functional prediction analysis revealed the presence of biologically relevant variants that correlate with the observed association signals. While the findings generated here warrant independent validation, the results of our large study suggest that ATM may be a novel locus associated with the risk of multiple subtypes of NHL.     

Complex Variation in Measures of General Intelligence and Cognitive Change

Combining information from multiple SNPs may capture a greater amount of genetic variation than from the sum of individual SNP effects and help identifying missing heritability. Regions may capture variation from multiple common variants of small effect, multiple rare variants or a combination of both. We describe regional heritability mapping of human cognition. Measures of crystallised (g_c) and fluid intelligence (g_f) in late adulthood (64–79 years) were available for 1806 individuals genotyped for 549,692 autosomal single nucleotide polymorphisms (SNPs). The same individuals were tested at age 11, enabling us the rare opportunity to measure cognitive change across most of their lifespan. 547,750 SNPs ranked by position are divided into 10, 908 overlapping regions of 101 SNPs to estimate the genetic variance each region explains, an approach that resembles classical linkage methods. We also estimate the genetic variation explained by individual autosomes and by SNPs within genes. Empirical significance thresholds are estimated separately for each trait from whole genome scans of 500 permutated data sets. The 5% significance threshold for the likelihood ratio test of a single region ranged from 17–17.5 for the three traits. This is the equivalent to nominal significance under the expectation of a chi-squared distribution (between 1df and 0) of P<1.44×10⁻⁵. These thresholds indicate that the distribution of the likelihood ratio test from this type of variance component analysis should be estimated empirically. Furthermore, we show that estimates of variation explained by these regions can be grossly overestimated. After applying permutation thresholds, a region for g_f on chromosome 5 spanning the PRRC1 gene is significant at a genome-wide 10% empirical threshold. Analysis of gene methylation on the temporal cortex provides support for the association of PRRC1 and fluid intelligence (P = 0.004), and provides a prime candidate gene for high throughput sequencing of these uniquely informative cohorts.     

Genetic Variation in the Base Excision Repair Pathway,Environmental Risk Factors,and Colorectal Adenoma Risk

Cigarette smoking, high alcohol intake, and low dietary folate levels are risk factors for colorectal adenomas. Oxidative damage caused by these three factors can be repaired through the base excision repair pathway (BER). We hypothesized that genetic variation in BER might modify colorectal adenoma risk. In a sigmoidoscopy-based study, we examined associations between 182 haplotype tagging SNPs in 14 BER genes, and colorectal adenoma risk, and examined their potential role as modifiers of the effect cigarette smoking, alcohol intake, and dietary folate levels. Among all individuals, no statistically significant associations between BER SNPs and adenoma risk persisted after correction for multiple comparisons. However, among Asian-Pacific Islanders we observed two SNPs in FEN1 and one in NTHL1, and among African-Americans one SNP in APEX1 that were associated with colorectal adenoma risk. Significant associations were also observed between SNPs in the NEIL2 gene and rectal adenoma risk. Three SNPS modified the effect of smoking (MUTYH interaction p = 0.002; OGG1 interaction p = 0.013); FEN1 interaction p = 0.013)), one SNP in LIG3 modified the effect of alcohol consumption (interaction p = 0.024) and two SNPs in LIG3 modified the effect of dietary folate (interaction p = 0.001 and p = 0.08) on colorectal adenoma risk. These findings support a role for genetic variants in the BER pathway as potential modifiers of colorectal adenoma risk. Our findings strengthen the role of oxidative damage induced by key lifestyle and dietary risk factors in colorectal adenoma formation.     

Polymorphisms within the Novel Type 2 Diabetes Risk Locus MTNR1B Determine β-Cell Function

Background

Very recently, a novel type 2 diabetes risk gene, i.e., MTNR1B, was identified and reported to affect fasting glycemia. Using our thoroughly phenotyped cohort of subjects at an increased risk for type 2 diabetes, we assessed the association of common genetic variation within the MTNR1B locus with obesity and prediabetes traits, namely impaired insulin secretion and insulin resistance.

Methodology/Principal Findings

We genotyped 1,578 non-diabetic subjects, metabolically characterized by oral glucose tolerance test, for five tagging single nucleotide polymorphisms (SNPs) covering 100% of common genetic variation (minor allele frequency >0.05) within the MTNR1B locus (rs10830962, rs4753426, rs12804291, rs10830963, rs3781638). In a subgroup (N = 513), insulin sensitivity was assessed by hyperinsulinemic-euglycemic clamp, and in a further subgroup (N = 301), glucose-stimulated insulin secretion was determined by intravenous glucose tolerance test. After appropriate adjustment for confounding variables and Bonferroni correction for multiple comparisons, none of the tagging SNPs was reliably associated with measures of adiposity. SNPs rs10830962, rs4753426, and rs10830963 were significantly associated with higher fasting plasma glucose concentrations (p<0.0001) and reduced OGTT- and IVGTT-induced insulin release (p≤0.0007 and p≤0.01, respectively). By contrast, SNP rs3781638 displayed significant association with lower fasting plasma glucose levels and increased OGTT-induced insulin release (p<0.0001 and p≤0.0002, respectively). Moreover, SNP rs3781638 revealed significant association with elevated fasting- and OGTT-derived insulin sensitivity (p≤0.0021). None of the MTNR1B tagging SNPs altered proinsulin-to-insulin conversion.

Conclusions/Significance

In conclusion, common genetic variation within MTNR1B determines glucose-stimulated insulin secretion and plasma glucose concentrations. Their impact on β-cell function might represent the prevailing pathomechanism how MTNR1B variants increase the type 2 diabetes risk.     

First Genome-Wide Association Study on Anxiety-Related Behaviours in Childhood

Background

Twin studies have shown that anxiety in a general population sample of children involves both domain-general and trait-specific genetic effects. For this reason, in an attempt to identify genes responsible for these effects, we investigated domain-general and trait-specific genetic associations in the first genome-wide association (GWA) study on anxiety-related behaviours (ARBs) in childhood.

Methods

The sample included 2810 7-year-olds drawn from the Twins Early Development Study (TEDS) with data available for parent-rated anxiety and genome-wide DNA markers. The measure was the Anxiety-Related Behaviours Questionnaire (ARBQ), which assesses four anxiety traits and also yields a general anxiety composite. Affymetrix GeneChip 6.0 DNA arrays were used to genotype nearly 700,000 single-nucleotide polymorphisms (SNPs), and IMPUTE v2 was used to impute more than 1 million SNPs. Several GWA associations from this discovery sample were followed up in another TEDS sample of 4804 children. In addition, Genome-wide Complex Trait Analysis (GCTA) was used on the discovery sample, to estimate the total amount of variance in ARBs that can be accounted for by SNPs on the array.

Results

No SNP associations met the demanding criterion of genome-wide significance that corrects for multiple testing across the genome (p<5×10⁻⁸). Attempts to replicate the top associations did not yield significant results. In contrast to the substantial twin study estimates of heritability which ranged from 0.50 (0.03) to 0.61 (0.01), the GCTA estimates of phenotypic variance accounted for by the SNPs were much lower 0.01 (0.11) to 0.19 (0.12).

Conclusions

Taken together, these GWAS and GCTA results suggest that anxiety – similar to height, weight and intelligence − is affected by many genetic variants of small effect, but unlike these other prototypical polygenic traits, genetic influence on anxiety is not well tagged by common SNPs.     

Genetic Variants in the Fat and Obesity Associated (FTO) Gene and Risk of Alzheimer's Disease

Background

Recent studies showed that polymorphisms in the Fat and Obesity-Associated (FTO) gene have robust effects on obesity, obesity-related traits and endophenotypes associated with Alzheimer''s disease (AD).

Methods

We used 1,877 Caucasian cases and controls from the NIA-LOAD study and 1,093 Caribbean Hispanics to further explore the association of FTO with AD. Using logistic regression, we assessed 42 SNPs in introns 1 and 2, the region previously reported to be associated with AD endophenotypes, which had been derived by genome-wide screenings. In addition, we performed gene expression analyses of neuropathologically confirmed AD cases and controls of two independent datasets (19 AD cases, 10 controls; 176 AD cases, 188 controls) using within- and between-group factors ANOVA of log₁₀ transformed rank invariant normalized expression data.

Results

In the NIALOAD study, one SNP was significantly associated with AD and three additional markers were close to significance (rs6499640, rs10852521, rs16945088, rs8044769, FDR p-value: 0.05<p<0.09). Two of the SNPs are in strong LD (D′>0.9) with the previously reported SNPs. In the Caribbean Hispanic dataset, we identified three SNPs (rs17219084, rs11075996, rs11075997, FDR p-value: 0.009<p<0.01) that were associated with AD. These results were confirmed by haplotype analyses and in a metaanalysis in which we included the ADNI dataset. FTO had a significantly lower expresssion in AD cases compared to controls in two independent datasets derived from human cortex and amygdala tissue, respectively (p = 2.18×10−5 and p<0.0001).

Conclusions

Our data support the notion that genetic variation in Introns 1 and 2 of the FTO gene may contribute to AD risk.     

Extent of height variability explained by known height-associated genetic variants in an isolated population of the Adriatic coast of Croatia

Background

Human height is a classical example of a polygenic quantitative trait. Recent large-scale genome-wide association studies (GWAS) have identified more than 200 height-associated loci, though these variants explain only 2∼10% of overall variability of normal height. The objective of this study was to investigate the variance explained by these loci in a relatively isolated population of European descent with limited admixture and homogeneous genetic background from the Adriatic coast of Croatia.

Methodology/Principal Findings

In a sample of 1304 individuals from the island population of Hvar, Croatia, we performed genome-wide SNP typing and assessed the variance explained by genetic scores constructed from different panels of height-associated SNPs extracted from five published studies. The combined information of the 180 SNPs reported by Lango Allen el al. explained 7.94% of phenotypic variation in our sample. Genetic scores based on 20∼50 SNPs reported by the remaining individual GWA studies explained 3∼5% of height variance. These percentages of variance explained were within ranges comparable to the original studies and heterogeneity tests did not detect significant differences in effect size estimates between our study and the original reports, if the estimates were obtained from populations of European descent.

Conclusions/Significance

We have evaluated the portability of height-associated loci and the overall fitting of estimated effect sizes reported in large cohorts to an isolated population. We found proportions of explained height variability were comparable to multiple reference GWAS in cohorts of European descent. These results indicate similar genetic architecture and comparable effect sizes of height loci among populations of European descent.     

IDO1 and IDO2 Non-Synonymous Gene Variants: Correlation with Crohn's Disease Risk and Clinical Phenotype

BackgroundCrohn''s disease (CD) is a chronic inflammatory disease of the gastrointestinal tract. Genetic polymorphisms can confer CD risk and influence disease phenotype. Indoleamine 2,3 dioxygenase-1 (IDO1) is one of the most over-expressed genes in CD and mediates potent anti-inflammatory effects via tryptophan metabolism along the kynurenine pathway. We aimed to determine whether non-synonymous polymorphisms in IDO1 or IDO2 (a gene paralog) are important either as CD risk alleles or as modifiers of CD phenotype.MethodsUtilizing a prospectively collected database, clinically phenotyped CD patients (n = 734) and non-IBD controls (n = 354) were genotyped for established IDO1 and IDO2 non-synonymous single nucleotide polymorphisms (SNPs) and novel genetic variants elucidated in the literature. Allelic frequencies between CD and non-IBD controls were compared. Genotype-phenotype analysis was conducted. IDO1 enzyme activity was assessed by calculating the serum kynurenine to tryptophan ratio (K/T).ResultsIDO1 SNPs were rare (1.7% non-IBD vs 1.1% CD; p = NS) and not linked to Crohn''s disease diagnosis in this population. IDO1 SNPs did however associate with a severe clinical course, presence of perianal disease, extraintestinal manifestations and a reduced serum K/T ratio during active disease suggesting lower IDO1 function. IDO2 minor allele variants were common and one of them, rs45003083, associated with reduced risk of Crohn''s disease (p = 0.025). No IDO2 SNPs associated with a particular Crohn''s disease clinical phenotype.ConclusionsThis work highlights the functional importance of IDO enzymes in human Crohn''s disease and establishes relative rates of IDO genetic variants in a US population.     

Genome-Wide Search for Gene-Gene Interactions in Colorectal Cancer

Genome-wide association studies (GWAS) have successfully identified a number of single-nucleotide polymorphisms (SNPs) associated with colorectal cancer (CRC) risk. However, these susceptibility loci known today explain only a small fraction of the genetic risk. Gene-gene interaction (GxG) is considered to be one source of the missing heritability. To address this, we performed a genome-wide search for pair-wise GxG associated with CRC risk using 8,380 cases and 10,558 controls in the discovery phase and 2,527 cases and 2,658 controls in the replication phase. We developed a simple, but powerful method for testing interaction, which we term the Average Risk Due to Interaction (ARDI). With this method, we conducted a genome-wide search to identify SNPs showing evidence for GxG with previously identified CRC susceptibility loci from 14 independent regions. We also conducted a genome-wide search for GxG using the marginal association screening and examining interaction among SNPs that pass the screening threshold (p<10⁻⁴). For the known locus rs10795668 (10p14), we found an interacting SNP rs367615 (5q21) with replication p = 0.01 and combined p = 4.19×10⁻⁸. Among the top marginal SNPs after LD pruning (n = 163), we identified an interaction between rs1571218 (20p12.3) and rs10879357 (12q21.1) (nominal combined p = 2.51×10⁻⁶; Bonferroni adjusted p = 0.03). Our study represents the first comprehensive search for GxG in CRC, and our results may provide new insight into the genetic etiology of CRC.