The effect of genetic drift in a young genetically isolated population

A cross-omics integrative study of metabolic signatures of chronic obstructive pulmonary disease

Background

Chronic obstructive pulmonary disease (COPD) is a common lung disorder characterized by persistent and progressive airflow limitation as well as systemic changes. Metabolic changes in blood may help detect COPD in an earlier stage and predict prognosis.

Methods

We conducted a comprehensive study of circulating metabolites, measured by proton Nuclear Magnetic Resonance Spectroscopy, in relation with COPD and lung function. The discovery sample consisted of 5557 individuals from two large population-based studies in the Netherlands, the Rotterdam Study and the Erasmus Rucphen Family study. Significant findings were replicated in 12,205 individuals from the Lifelines-DEEP study, FINRISK and the Prospective Investigation of the Vasculature in Uppsala Seniors (PIVUS) studies. For replicated metabolites further investigation of causality was performed, utilizing genetics in the Mendelian randomization approach.

Results

There were 602 cases of COPD and 4955 controls used in the discovery meta-analysis. Our logistic regression results showed that higher levels of plasma Glycoprotein acetyls (GlycA) are significantly associated with COPD (OR = 1.16, P = 5.6 × 10^− 4 in the discovery and OR = 1.30, P = 1.8 × 10^− 6 in the replication sample). A bi-directional two-sample Mendelian randomization analysis suggested that circulating blood GlycA is not causally related to COPD, but that COPD causally increases GlycA levels. Using the prospective data of the same sample of Rotterdam Study in Cox-regression, we show that the circulating GlycA level is a predictive biomarker of COPD incidence (HR = 1.99, 95%CI 1.52–2.60, comparing those in the highest and lowest quartile of GlycA) but is not significantly associated with mortality in COPD patients (HR = 1.07, 95%CI 0.94–1.20).

Conclusions

Our study shows that circulating blood GlycA is a biomarker of early COPD pathology.

Exome Sequencing Analysis Identifies Rare Variants in ATM and RPL8 That Are Associated With Shorter Telomere Length

Telomeres are important for maintaining genomic stability. Telomere length has been associated with aging, disease, and mortality and is highly heritable (∼82%). In this study, we aimed to identify rare genetic variants associated with telomere length using whole-exome sequence data. We studied 1,303 participants of the Erasmus Rucphen Family (ERF) study, 1,259 of the Rotterdam Study (RS), and 674 of the British Heart Foundation Family Heart Study (BHF-FHS). We conducted two analyses, first we analyzed the family-based ERF study and used the RS and BHF-FHS for replication. Second, we combined the summary data of the three studies in a meta-analysis. Telomere length was measured by quantitative polymerase chain reaction in blood. We identified nine rare variants significantly associated with telomere length (p-value < 1.42 × 10^–7, minor allele frequency of 0.2–0.5%) in the ERF study. Eight of these variants (in C11orf65, ACAT1, NPAT, ATM, KDELC2, and EXPH5) were located on chromosome 11q22.3 that contains ATM, a gene involved in telomere maintenance. Although we were unable to replicate the variants in the RS and BHF-FHS (p-value ≥ 0.21), segregation analysis showed that all variants segregate with shorter telomere length in a family. In the meta-analysis of all studies, a nominally significant association with LTL was observed with a rare variant in RPL8 (p-value = 1.48 × 10⁻⁶), which has previously been associated with age. Additionally, a novel rare variant in the known RTEL1 locus showed suggestive evidence for association (p-value = 1.18 × 10^–4) with LTL. To conclude, we identified novel rare variants associated with telomere length. Larger samples size are needed to confirm these findings and to identify additional variants.

Defining the genetic control of human blood plasma N-glycome using genome-wide association study

Glycosylation is a common post-translational modification of proteins. Glycosylation is associated with a number of human diseases. Defining genetic factors altering glycosylation may provide a basis for novel approaches to diagnostic and pharmaceutical applications. Here we report a genome-wide association study of the human blood plasma N-glycome composition in up to 3811 people measured by Ultra Performance Liquid Chromatography (UPLC) technology. Starting with the 36 original traits measured by UPLC, we computed an additional 77 derived traits leading to a total of 113 glycan traits. We studied associations between these traits and genetic polymorphisms located on human autosomes. We discovered and replicated 12 loci. This allowed us to demonstrate an overlap in genetic control between total plasma protein and IgG glycosylation. The majority of revealed loci contained genes that encode enzymes directly involved in glycosylation (FUT3/FUT6, FUT8, B3GAT1, ST6GAL1, B4GALT1, ST3GAL4, MGAT3 and MGAT5) and a known regulator of plasma protein fucosylation (HNF1A). However, we also found loci that could possibly reflect other more complex aspects of glycosylation process. Functional genomic annotation suggested the role of several genes including DERL3, CHCHD10, TMEM121, IGH and IKZF1. The hypotheses we generated may serve as a starting point for further functional studies in this research area.

Allele frequency dynamics in a pedigreed natural population

A central goal of population genetics is to understand how genetic drift, natural selection, and gene flow shape allele frequencies through time. However, the actual processes underlying these changes-variation in individual survival, reproductive success, and movement-are often difficult to quantify. Fully understanding these processes requires the population pedigree, the set of relationships among all individuals in the population through time. Here, we use extensive pedigree and genomic information from a long-studied natural population of Florida Scrub-Jays (Aphelocoma coerulescens) to directly characterize the relative roles of different evolutionary processes in shaping patterns of genetic variation through time. We performed gene dropping simulations to estimate individual genetic contributions to the population and model drift on the known pedigree. We found that observed allele frequency changes are generally well predicted by accounting for the different genetic contributions of founders. Our results show that the genetic contribution of recent immigrants is substantial, with some large allele frequency shifts that otherwise may have been attributed to selection actually due to gene flow. We identified a few SNPs under directional short-term selection after appropriately accounting for gene flow. Using models that account for changes in population size, we partitioned the proportion of variance in allele frequency change through time. Observed allele frequency changes are primarily due to variation in survival and reproductive success, with gene flow making a smaller contribution. This study provides one of the most complete descriptions of short-term evolutionary change in allele frequencies in a natural population to date.

Meta-analysis of genome-wide association studies identifies 8 novel loci involved in shape variation of human head hair

Shape variation of human head hair shows striking variation within and between human populations, while its genetic basis is far from being understood. We performed a series of genome-wide association studies (GWASs) and replication studies in a total of 28 964 subjects from 9 cohorts from multiple geographic origins. A meta-analysis of three European GWASs identified 8 novel loci (1p36.23 ERRFI1/SLC45A1, 1p36.22 PEX14, 1p36.13 PADI3, 2p13.3 TGFA, 11p14.1 LGR4, 12q13.13 HOXC13, 17q21.2 KRTAP, and 20q13.33 PTK6), and confirmed 4 previously known ones (1q21.3 TCHH/TCHHL1/LCE3E, 2q35 WNT10A, 4q21.21 FRAS1, and 10p14 LINC00708/GATA3), all showing genome-wide significant association with hair shape (P < 5e-8). All except one (1p36.22 PEX14) were replicated with nominal significance in at least one of the 6 additional cohorts of European, Native American and East Asian origins. Three additional previously known genes (EDAR, OFCC1, and PRSS53) were confirmed at the nominal significance level. A multivariable regression model revealed that 14 SNPs from different genes significantly and independently contribute to hair shape variation, reaching a cross-validated AUC value of 0.66 (95% CI: 0.62–0.70) and an AUC value of 0.64 in an independent validation cohort, providing an improved accuracy compared with a previous model. Prediction outcomes of 2504 individuals from a multiethnic sample were largely consistent with general knowledge on the global distribution of hair shape variation. Our study thus delivers target genes and DNA variants for future functional studies to further evaluate the molecular basis of hair shape in humans.

Whole-Genome Linkage Scan Combined With Exome Sequencing Identifies Novel Candidate Genes for Carotid Intima-Media Thickness

Carotid intima-media thickness (cIMT) is an established heritable marker for subclinical atherosclerosis. In this study, we aim to identify rare variants with large effects driving differences in cIMT by performing genome-wide linkage analysis of individuals in the extremes of cIMT trait distribution (>90th percentile) in a large family-based study from a genetically isolated population in the Netherlands. Linked regions were subsequently explored by fine-mapping using exome sequencing. We observed significant evidence of linkage on chromosomes 2p16.3 [rs1017418, heterogeneity LOD (HLOD) = 3.35], 19q13.43 (rs3499, HLOD = 9.09), 20p13 (rs1434789, HLOD = 4.10), and 21q22.12 (rs2834949, HLOD = 3.59). Fine-mapping using exome sequencing data identified a non-coding variant (rs62165235) in PNPT1 gene under the linkage peak at chromosome 2 that is likely to have a regulatory function. The variant was associated with quantitative cIMT in the family-based study population (effect = 0.27, p-value = 0.013). Furthermore, we identified several genes under the linkage peak at chromosome 21 highly expressed in tissues relevant for atherosclerosis. To conclude, our linkage analysis identified four genomic regions significantly linked to cIMT. Further analyses are needed to demonstrate involvement of identified candidate genes in development of atherosclerosis.

Detection of Significant Association Between Variants in Cannabinoid Receptor 1 Gene (CNR1) and Personality in African-American Population

Background: Several studies have revealed significant associations between single nucleotide polymorphisms (SNPs) in the cannabinoid receptor 1 (CNR1) gene and a broad spectrum of psychiatric disorders such as major depressive disorder (MDD), attention deficit hyperactivity disorder (ADHD), and schizophrenia. Personality traits that are highly related to susceptibility to these conditions have been associated with the CNR1 variants in subjects of Caucasian origin. However, there are no reported studies regarding the effects of CNR1 polymorphisms on personality traits in the African-American (AA) population. Methods: We performed an imputation-based association analysis for 26 CNR1 variants with five dimensions of personality in 3,046 AAs. Results: SNPs rs806372 and rs2180619 showed a significant association with extraversion after Bonferroni correction for multiple testing (p < 0.0019). Further, several extraversion-associated SNPs were significantly associated with conscientiousness, agreeableness, and openness. SNP priority score analysis indicated that SNPs rs806368, rs806371, and rs2180619 play a role in the modulation of personality and psychiatric conditions. Conclusion:CNR1 is important in determining personality traits in the AA population.

A rare missense variant in RCL1 segregates with depression in extended families

Depression is the most prevalent psychiatric disorder with a complex and elusive etiology that is moderately heritable. Identification of genes would greatly facilitate the elucidation of the biological mechanisms underlying depression, however, its complex etiology has proved to be a major bottleneck in the identification of its genetic risk factors, especially in genome-wide association-like studies. In this study, we exploit the properties of a genetic isolate and its family-based structure to explore whether relatively rare exonic variants influence the burden of depressive symptoms in families. Using a multistep approach involving linkage and haplotype analyses followed by exome sequencing in the Erasmus Rucphen Family (ERF) study, we identified a rare (minor allele frequency (MAF)=1%) missense c.1114C>T mutation (rs115482041) in the RCL1 gene segregating with depression across multiple generations. Rs115482041 showed significant association with depressive symptoms (N=2393, β_T-allele=2.33, P-value=1 × 10^-4) and explained 2.9% of the estimated genetic variance of depressive symptoms (22%) in ERF. Despite being twice as rare (MAF<0.5%), c.1114C>T showed similar effect and significant association with depressive symptoms in samples from the independent population-based Rotterdam study (N=1604, β_T-allele=3.60, P-value=3 × 10^-2). A comparison of RCL1 expression in human and mouse brain revealed a striking co-localization of RCL1 with the layer 1 interlaminar subclass of astrocytes found exclusively in higher-order primates. Our findings identify RCL1 as a novel candidate gene for depression and offer insights into mechanisms through which RCL1 may be relevant for depression.

A Genome-Wide Linkage Study for Chronic Obstructive Pulmonary Disease in a Dutch Genetic Isolate Identifies Novel Rare Candidate Variants

Chronic obstructive pulmonary disease (COPD) is a complex and heritable disease, associated with multiple genetic variants. Specific familial types of COPD may be explained by rare variants, which have not been widely studied. We aimed to discover rare genetic variants underlying COPD through a genome-wide linkage scan. Affected-only analysis was performed using the 6K Illumina Linkage IV Panel in 142 cases clustered in 27 families from a genetic isolate, the Erasmus Rucphen Family (ERF) study. Potential causal variants were identified by searching for shared rare variants in the exome-sequence data of the affected members of the families contributing most to the linkage peak. The identified rare variants were then tested for association with COPD in a large meta-analysis of several cohorts. Significant evidence for linkage was observed on chromosomes 15q14-15q25 [logarithm of the odds (LOD) score = 5.52], 11p15.4-11q14.1 (LOD = 3.71) and 5q14.3-5q33.2 (LOD = 3.49). In the chromosome 15 peak, that harbors the known COPD locus for nicotinic receptors, and in the chromosome 5 peak we could not identify shared variants. In the chromosome 11 locus, we identified four rare (minor allele frequency (MAF) <0.02), predicted pathogenic, missense variants. These were shared among the affected family members. The identified variants localize to genes including neuroblast differentiation-associated protein (AHNAK), previously associated with blood biomarkers in COPD, phospholipase C Beta 3 (PLCB3), shown to increase airway hyper-responsiveness, solute carrier family 22-A11 (SLC22A11), involved in amino acid metabolism and ion transport, and metallothionein-like protein 5 (MTL5), involved in nicotinate and nicotinamide metabolism. Association of SLC22A11 and MTL5 variants were confirmed in the meta-analysis of 9,888 cases and 27,060 controls. In conclusion, we have identified novel rare variants in plausible genes related to COPD. Further studies utilizing large sample whole-genome sequencing should further confirm the associations at chromosome 11 and investigate the chromosome 15 and 5 linked regions.

A combined linkage, microarray and exome analysis suggests MAP3K11 as a candidate gene for left ventricular hypertrophy

BACKGROUND

Electrocardiographic measures of left ventricular hypertrophy (LVH) are used as predictors of cardiovascular risk. We combined linkage and association analyses to discover novel rare genetic variants involved in three such measures and two principal components derived from them.

METHODS

The study was conducted among participants from the Erasmus Rucphen Family Study (ERF), a Dutch family-based sample from the southwestern Netherlands. Variance components linkage analyses were performed using Merlin. Regions of interest (LOD > 1.9) were fine-mapped using microarray and exome sequence data.

RESULTS

We observed one significant LOD score for the second principal component on chromosome 15 (LOD score = 3.01) and 12 suggestive LOD scores. Several loci contained variants identified in GWAS for these traits; however, these did not explain the linkage peaks, nor did other common variants. Exome sequence data identified two associated variants after multiple testing corrections were applied.

CONCLUSIONS

We did not find common SNPs explaining these linkage signals. Exome sequencing uncovered a relatively rare variant in MAPK3K11 on chromosome 11 (MAF = 0.01) that helped account for the suggestive linkage peak observed for the first principal component. Conditional analysis revealed a drop in LOD from 2.01 to 0.88 for MAP3K11, suggesting that this variant may partially explain the linkage signal at this chromosomal location. MAP3K11 is related to the JNK pathway and is a pro-apoptotic kinase that plays an important role in the induction of cardiomyocyte apoptosis in various pathologies, including LVH.

Sea-level driven glacial-age refugia and post-glacial mixing on subtropical coasts, a palaeohabitat and genetic study

Using a novel combination of palaeohabitat modelling and genetic mixture analyses, we identify and assess a sea-level-driven recolonization process following the Last Glacial Maximum (LGM). Our palaeohabitat modelling reveals dramatic changes in estuarine habitat distribution along the coast of California (USA) and Baja California (Mexico). At the LGM (approx. 20 kya), when sea level was approximately 130 m lower, the palaeo-shoreline was too steep for tidal estuarine habitat formation, eliminating this habitat type from regions where it is currently most abundant, and limiting such estuaries to a northern and a southern refugium separated by 1000 km. We assess the recolonization of estuaries formed during post-LGM sea-level rise through examination of refugium-associated alleles and approximate Bayesian computation in three species of estuarine fishes. Results reveal sourcing of modern populations from both refugia, which admix in the newly formed habitat between the refuges. We infer a dramatic peak in habitat area between 15 and 10 kya with subsequent decline. Overall, this approach revealed a previously undocumented dynamic and integrated relationship between sea-level change, coastal processes and population genetics. These results extend glacial refugial dynamics to unglaciated subtropical coasts and have significant implications for biotic response to predicted sea-level rise.

Exome-Wide Meta-Analysis Identifies Rare 3'-UTR Variant in ERCC1/CD3EAP Associated with Symptoms of Sleep Apnea

Obstructive sleep apnea (OSA) is a common sleep breathing disorder associated with an increased risk of cardiovascular and cerebrovascular diseases and mortality. Although OSA is fairly heritable (~40%), there have been only few studies looking into the genetics of OSA. In the present study, we aimed to identify genetic variants associated with symptoms of sleep apnea by performing a whole-exome sequence meta-analysis of symptoms of sleep apnea in 1,475 individuals of European descent. We identified 17 rare genetic variants with at least suggestive evidence of significance. Replication in an independent dataset confirmed the association of a rare genetic variant (rs2229918; minor allele frequency = 0.3%) with symptoms of sleep apnea (p-value_meta = 6.98 × 10⁻⁹, β_meta = 0.99). Rs2229918 overlaps with the 3′ untranslated regions of ERCC1 and CD3EAP genes on chromosome 19q13. Both genes are expressed in tissues in the neck area, such as the tongue, muscles, cartilage and the trachea. Further, CD3EAP is localized in the nucleus and mitochondria and involved in the tumor necrosis factor-alpha/nuclear factor kappa B signaling pathway. Our results and biological functions of CD3EAP/ERCC1 genes suggest that the 19q13 locus is interesting for further OSA research.

Variants in TTC25 affect autistic trait in patients with autism spectrum disorder and general population

Autism spectrum disorder (ASD) is a highly heritable neurodevelopmental disorder with a complex genetic architecture. To identify genetic variants underlying ASD, we performed single-variant and gene-based genome-wide association studies using a dense genotyping array containing over 2.3 million single-nucleotide variants in a discovery sample of 160 families with at least one child affected with non-syndromic ASD using a binary (ASD yes/no) phenotype and a quantitative autistic trait. Replication of the top findings was performed in Psychiatric Genomics Consortium and Erasmus Rucphen Family (ERF) cohort study. Significant association of quantitative autistic trait was observed with the TTC25 gene at 17q21.2 (effect size=10.2, P-value=3.4 × 10^-7) in the gene-based analysis. The gene also showed nominally significant association in the cohort-based ERF study (effect=1.75, P-value=0.05). Meta-analysis of discovery and replication improved the association signal (P-value_meta=1.5 × 10^-8). No genome-wide significant signal was observed in the single-variant analysis of either the binary ASD phenotype or the quantitative autistic trait. Our study has identified a novel gene TTC25 to be associated with quantitative autistic trait in patients with ASD. The replication of association in a cohort-based study and the effect estimate suggest that variants in TTC25 may also be relevant for broader ASD phenotype in the general population. TTC25 is overexpressed in frontal cortex and testis and is known to be involved in cilium movement and thus an interesting candidate gene for autistic trait.

Exome-sequencing in a large population-based study reveals a rare Asn396Ser variant in the LIPG gene associated with depressive symptoms

Despite a substantial genetic component, efforts to identify common genetic variation underlying depression have largely been unsuccessful. In the current study we aimed to identify rare genetic variants that might have large effects on depression in the general population. Using high-coverage exome-sequencing, we studied the exonic variants in 1265 individuals from the Rotterdam study (RS), who were assessed for depressive symptoms. We identified a missense Asn396Ser mutation (rs77960347) in the endothelial lipase (LIPG) gene, occurring with an allele frequency of 1% in the general population, which was significantly associated with depressive symptoms (P-value=5.2 × 10^-08, β=7.2). Replication in three independent data sets (N=3612) confirmed the association of Asn396Ser (P-value=7.1 × 10^-03, β=2.55) with depressive symptoms. LIPG is predicted to have enzymatic function in steroid biosynthesis, cholesterol biosynthesis and thyroid hormone metabolic processes. The Asn396Ser variant is predicted to have a damaging effect on the function of LIPG. Within the discovery population, carriers also showed an increased burden of white matter lesions (P-value=3.3 × 10^-02) and a higher risk of Alzheimer's disease (odds ratio=2.01; P-value=2.8 × 10^-02) compared with the non-carriers. Together, these findings implicate the Asn396Ser variant of LIPG in the pathogenesis of depressive symptoms in the general population.

Discovery and functional prioritization of Parkinson's disease candidate genes from large-scale whole exome sequencing

BACKGROUND

Whole-exome sequencing (WES) has been successful in identifying genes that cause familial Parkinson's disease (PD). However, until now this approach has not been deployed to study large cohorts of unrelated participants. To discover rare PD susceptibility variants, we performed WES in 1148 unrelated cases and 503 control participants. Candidate genes were subsequently validated for functions relevant to PD based on parallel RNA-interference (RNAi) screens in human cell culture and Drosophila and C. elegans models.

RESULTS

Assuming autosomal recessive inheritance, we identify 27 genes that have homozygous or compound heterozygous loss-of-function variants in PD cases. Definitive replication and confirmation of these findings were hindered by potential heterogeneity and by the rarity of the implicated alleles. We therefore looked for potential genetic interactions with established PD mechanisms. Following RNAi-mediated knockdown, 15 of the genes modulated mitochondrial dynamics in human neuronal cultures and four candidates enhanced α-synuclein-induced neurodegeneration in Drosophila. Based on complementary analyses in independent human datasets, five functionally validated genes-GPATCH2L, UHRF1BP1L, PTPRH, ARSB, and VPS13C-also showed evidence consistent with genetic replication.

CONCLUSIONS

By integrating human genetic and functional evidence, we identify several PD susceptibility gene candidates for further investigation. Our approach highlights a powerful experimental strategy with broad applicability for future studies of disorders with complex genetic etiologies.

A Combined Linkage and Exome Sequencing Analysis for Electrocardiogram Parameters in the Erasmus Rucphen Family Study

Electrocardiogram (ECG) measurements play a key role in the diagnosis and prediction of cardiac arrhythmias and sudden cardiac death. ECG parameters, such as the PR, QRS, and QT intervals, are known to be heritable and genome-wide association studies of these phenotypes have been successful in identifying common variants; however, a large proportion of the genetic variability of these traits remains to be elucidated. The aim of this study was to discover loci potentially harboring rare variants utilizing variance component linkage analysis in 1547 individuals from a large family-based study, the Erasmus Rucphen Family Study (ERF). Linked regions were further explored using exome sequencing. Five suggestive linkage peaks were identified: two for QT interval (1q24, LOD = 2.63; 2q34, LOD = 2.05), one for QRS interval (1p35, LOD = 2.52) and two for PR interval (9p22, LOD = 2.20; 14q11, LOD = 2.29). Fine-mapping using exome sequence data identified a C > G missense variant (c.713C > G, p.Ser238Cys) in the FCRL2 gene associated with QT (rs74608430; P = 2.8 × 10^-4, minor allele frequency = 0.019). Heritability analysis demonstrated that the SNP explained 2.42% of the trait’s genetic variability in ERF (P = 0.02). Pathway analysis suggested that the gene is involved in cytosolic Ca²⁺ levels (P = 3.3 × 10^-3) and AMPK stimulated fatty acid oxidation in muscle (P = 4.1 × 10^-3). Look-ups in bioinformatics resources showed that expression of FCRL2 is associated with ARHGAP24 and SETBP1 expression. This finding was not replicated in the Rotterdam study. Combining the bioinformatics information with the association and linkage analyses, FCRL2 emerges as a strong candidate gene for QT interval.

Associations with intraocular pressure across Europe: The European Eye Epidemiology (E3) Consortium

Raised intraocular pressure (IOP) is the most important risk factor for developing glaucoma, the second commonest cause of blindness globally. Understanding associations with IOP and variations in IOP between countries may teach us about mechanisms underlying glaucoma. We examined cross-sectional associations with IOP in 43,500 European adults from 12 cohort studies belonging to the European Eye Epidemiology (E³) consortium. Each study conducted multivariable linear regression with IOP as the outcome variable and results were pooled using random effects meta-analysis. The association of standardized study IOP with latitude was tested using meta-regression. Higher IOP was observed in men (0.18 mmHg; 95 % CI 0.06, 0.31; P = 0.004) and with higher body mass index (0.21 mmHg per 5 kg/m²; 95 % CI 0.14, 0.28; P < 0.001), shorter height (−0.17 mmHg per 10 cm; 95 % CI –0.25, −0.08; P < 0.001), higher systolic blood pressure (0.17 mmHg per 10 mmHg; 95 % CI 0.12, 0.22; P < 0.001) and more myopic refraction (0.06 mmHg per Dioptre; 95 % CI 0.03, 0.09; P < 0.001). An inverted U-shaped trend was observed between age and IOP, with IOP increasing up to the age of 60 and decreasing in participants older than 70 years. We found no significant association between standardized IOP and study location latitude (P = 0.76). Novel findings of our study include the association of lower IOP in taller people and an inverted-U shaped association of IOP with age. We found no evidence of significant variation in IOP across Europe. Despite the limited range of latitude amongst included studies, this finding is in favour of collaborative pooling of data from studies examining environmental and genetic determinants of IOP in Europeans.

Candidate SNP Markers of Gender-Biased Autoimmune Complications of Monogenic Diseases Are Predicted by a Significant Change in the Affinity of TATA-Binding Protein for Human Gene Promoters

Some variations of human genome [for example, single nucleotide polymorphisms (SNPs)] are markers of hereditary diseases and drug responses. Analysis of them can help to improve treatment. Computer-based analysis of millions of SNPs in the 1000 Genomes project makes a search for SNP markers more targeted. Here, we combined two computer-based approaches: DNA sequence analysis and keyword search in databases. In the binding sites for TATA-binding protein (TBP) in human gene promoters, we found candidate SNP markers of gender-biased autoimmune diseases, including rs1143627 [cachexia in rheumatoid arthritis (double prevalence among women)]; rs11557611 [demyelinating diseases (thrice more prevalent among young white women than among non-white individuals)]; rs17231520 and rs569033466 [both: atherosclerosis comorbid with related diseases (double prevalence among women)]; rs563763767 [Hughes syndrome-related thrombosis (lethal during pregnancy)]; rs2814778 [autoimmune diseases (excluding multiple sclerosis and rheumatoid arthritis) underlying hypergammaglobulinemia in women]; rs72661131 and rs562962093 (both: preterm delivery in pregnant diabetic women); and rs35518301, rs34166473, rs34500389, rs33981098, rs33980857, rs397509430, rs34598529, rs33931746, rs281864525, and rs63750953 (all: autoimmune diseases underlying hypergammaglobulinemia in women). Validation of these predicted candidate SNP markers using the clinical standards may advance personalized medicine.

Rare Exome Sequence Variants in CLCN6 Reduce Blood Pressure Levels and Hypertension Risk

BACKGROUND

Rare genetic variants influence blood pressure (BP).

METHODS AND RESULTS

Whole-exome sequencing was performed on DNA samples from 17 956 individuals of European ancestry and African ancestry (14 497, first-stage discovery and 3459, second-stage discovery) to examine the effect of rare variants on hypertension and 4 BP traits: systolic BP, diastolic BP, pulse pressure, and mean arterial pressure. Tests of ≈170 000 common variants (minor allele frequency, ≥1%; statistical significance, P≤2.9×10(-7)) and gene-based tests of rare variants (minor allele frequency, <1%; ≈17 000 genes; statistical significance, P≤1.5×10(-6)) were evaluated for each trait and ancestry, followed by multiethnic meta-analyses. In the first-stage discovery, rare coding variants (splicing, stop-gain, stop-loss, nonsynonymous variants, or indels) in CLCN6 were associated with lower diastolic BP (cumulative minor allele frequency, 1.3%; β=-3.20; P=4.1×10(-6)) and were independent of a nearby common variant (rs17367504) previously associated with BP. CLCN6 rare variants were also associated with lower systolic BP (β=-4.11; P=2.8×10(-4)), mean arterial pressure (β=-3.50; P=8.9×10(-6)), and reduced hypertension risk (odds ratio, 0.72; P=0.017). Meta-analysis of the 2-stage discovery samples showed that CLCN6 was associated with lower diastolic BP at exome-wide significance (cumulative minor allele frequency, 1.1%; β=-3.30; P=5.0×10(-7)).

CONCLUSIONS

These findings implicate the effect of rare coding variants in CLCN6 in BP variation and offer new insights into BP regulation.

Prediction of male-pattern baldness from genotypes

The global demand for products that effectively prevent the development of male-pattern baldness (MPB) has drastically increased. However, there is currently no established genetic model for the estimation of MPB risk. We conducted a prediction analysis using single-nucleotide polymorphisms (SNPs) identified from previous GWASs of MPB in a total of 2725 German and Dutch males. A logistic regression model considering the genotypes of 25 SNPs from 12 genomic loci demonstrates that early-onset MPB risk is predictable at an accuracy level of 0.74 when 14 SNPs were included in the model, and measured using the area under the receiver-operating characteristic curves (AUC). Considering age as an additional predictor, the model can predict normal MPB status in middle-aged and elderly individuals at a slightly lower accuracy (AUC 0.69-0.71) when 6-11 SNPs were used. A variance partitioning analysis suggests that 55.8% of early-onset MPB genetic liability can be explained by common autosomal SNPs and 23.3% by X-chromosome SNPs. For normal MPB status in elderly individuals, the proportion of explainable variance is lower (42.4% for autosomal and 9.8% for X-chromosome SNPs). The gap between GWAS findings and the variance partitioning results could be explained by a large body of common DNA variants with small effects that will likely be identified in GWAS of increased sample sizes. Although the accuracy obtained here has not reached a clinically desired level, our model was highly informative for up to 19% of Europeans, thus may assist decision making on early MPB intervention actions and in forensic investigations.

Heritabilities, proportions of heritabilities explained by GWAS findings, and implications of cross-phenotype effects on PR interval

Electrocardiogram (ECG) measurements are a powerful tool for evaluating cardiac function and are widely used for the diagnosis and prediction of a variety of conditions, including myocardial infarction, cardiac arrhythmias, and sudden cardiac death. Recently, genome-wide association studies (GWASs) identified a large number of genes related to ECG parameter variability, specifically for the QT, QRS, and PR intervals. The aims of this study were to establish the heritability of ECG traits, including indices of left ventricular hypertrophy, and to directly assess the proportion of those heritabilities explained by GWAS variants. These analyses were conducted in a large, Dutch family-based cohort study, the Erasmus Rucphen Family study using variance component methods implemented in the SOLAR (Sequential Oligogenic Linkage Analysis Routines) software package. Heritability estimates ranged from 34 % for QRS and Cornell voltage product to 49 % for 12-lead sum. Trait-specific GWAS findings for each trait explained a fraction of their heritability (17 % for QRS, 4 % for QT, 2 % for PR, 3 % for Sokolow–Lyon index, and 4 % for 12-lead sum). The inclusion of all ECG-associated single nucleotide polymorphisms explained an additional 6 % of the heritability of PR. In conclusion, this study shows that, although GWAS explain a portion of ECG trait variability, a large amount of heritability remains to be explained. In addition, larger GWAS for PR are likely to detect loci already identified, particularly those observed for QRS and 12-lead sum.

Genome-wide association analysis on five isolated populations identifies variants of the HLA-DOA gene associated with white wine liking

Wine is the most popular alcoholic beverage around the world and because of its importance in society has been widely studied. Understanding what drives its flavor has been a quest for decades but much is still unknown and will be determined at least in part by individual taste preferences. Recently studies in the genetics of taste have uncovered the role of different genes in the determination of food preferences giving new insight on its physiology. In this context we have performed a genome-wide association study on red and white wine liking using three isolated populations collected in Italy, and replicated our results on two additional populations coming from the Netherland and Central Asia for a total of 3885 samples. We have found a significant association (P=2.1 × 10(-8)) between white wine liking and rs9276975:C>T a polymorphism in the HLA-DOA gene encoding a non-canonical MHC II molecule, which regulates other MHC II molecules. The same association was also found with red wine liking (P=8.3 × 10(-6)). Sex-separated analysis have also revealed that the effect of HLA-DOA is twice as large in women as compared to men suggesting an interaction between this polymorphism and gender. Our results are one of the first examples of genome-wide association between liking of a commonly consumed food and gene variants. Moreover, our results suggest a role of the MHC system in the determination of food preferences opening new insight in this field in general.

Genome of The Netherlands population-specific imputations identify an ABCA6 variant associated with cholesterol levels

Variants associated with blood lipid levels may be population-specific. To identify low-frequency variants associated with this phenotype, population-specific reference panels may be used. Here we impute nine large Dutch biobanks (~35,000 samples) with the population-specific reference panel created by the Genome of the Netherlands Project and perform association testing with blood lipid levels. We report the discovery of five novel associations at four loci (P value <6.61 × 10⁻⁴), including a rare missense variant in ABCA6 (rs77542162, p.Cys1359Arg, frequency 0.034), which is predicted to be deleterious. The frequency of this ABCA6 variant is 3.65-fold increased in the Dutch and its effect (β_LDL-C=0.135, β_TC=0.140) is estimated to be very similar to those observed for single variants in well-known lipid genes, such as LDLR.

Frequencies of rare variants fluctuate over populations, hampering gene discovery. Here the authors use a population-specific reference panel, the Genome of the Netherlands, to discover four novel loci involved in lipid metabolism, including an exonic variant in ABCA6.

Genome-wide association study for refractive astigmatism reveals genetic co-determination with spherical equivalent refractive error: the CREAM consortium

To identify genetic variants associated with refractive astigmatism in the general population, meta-analyses of genome-wide association studies were performed for: White Europeans aged at least 25 years (20 cohorts, N = 31,968); Asian subjects aged at least 25 years (7 cohorts, N = 9,295); White Europeans aged <25 years (4 cohorts, N = 5,640); and all independent individuals from the above three samples combined with a sample of Chinese subjects aged <25 years (N = 45,931). Participants were classified as cases with refractive astigmatism if the average cylinder power in their two eyes was at least 1.00 diopter and as controls otherwise. Genome-wide association analysis was carried out for each cohort separately using logistic regression. Meta-analysis was conducted using a fixed effects model. In the older European group the most strongly associated marker was downstream of the neurexin-1 (NRXN1) gene (rs1401327, P = 3.92E-8). No other region reached genome-wide significance, and association signals were lower for the younger European group and Asian group. In the meta-analysis of all cohorts, no marker reached genome-wide significance: The most strongly associated regions were, NRXN1 (rs1401327, P = 2.93E-07), TOX (rs7823467, P = 3.47E-07) and LINC00340 (rs12212674, P = 1.49E-06). For 34 markers identified in prior GWAS for spherical equivalent refractive error, the beta coefficients for genotype versus spherical equivalent, and genotype versus refractive astigmatism, were highly correlated (r = -0.59, P = 2.10E-04). This work revealed no consistent or strong genetic signals for refractive astigmatism; however, the TOX gene region previously identified in GWAS for spherical equivalent refractive error was the second most strongly associated region. Analysis of additional markers provided evidence supporting widespread genetic co-susceptibility for spherical and astigmatic refractive errors.

The challenges of genome-wide interaction studies: lessons to learn from the analysis of HDL blood levels

Genome-wide association studies (GWAS) have revealed 74 single nucleotide polymorphisms (SNPs) associated with high-density lipoprotein cholesterol (HDL) blood levels. This study is, to our knowledge, the first genome-wide interaction study (GWIS) to identify SNP×SNP interactions associated with HDL levels. We performed a GWIS in the Rotterdam Study (RS) cohort I (RS-I) using the GLIDE tool which leverages the massively parallel computing power of Graphics Processing Units (GPUs) to perform linear regression on all genome-wide pairs of SNPs. By performing a meta-analysis together with Rotterdam Study cohorts II and III (RS-II and RS-III), we were able to filter 181 interaction terms with a p-value<1 · 10⁻⁸ that replicated in the two independent cohorts. We were not able to replicate any of these interaction term in the AGES, ARIC, CHS, ERF, FHS and NFBC-66 cohorts (N_total = 30,011) when adjusting for multiple testing. Our GWIS resulted in the consistent finding of a possible interaction between rs774801 in ARMC8 (ENSG00000114098) and rs12442098 in SPATA8 (ENSG00000185594) being associated with HDL levels. However, p-values do not reach the preset Bonferroni correction of the p-values. Our study suggest that even for highly genetically determined traits such as HDL the sample sizes needed to detect SNP×SNP interactions are large and the 2-step filtering approaches do not yield a solution. Here we present our analysis plan and our reservations concerning GWIS.

The dystrophin gene and cognitive function in the general population

The aim of our study is to investigate whether single-nucleotide dystrophin gene (DMD) variants associate with variability in cognitive functions in healthy populations. The study included 1240 participants from the Erasmus Rucphen family (ERF) study and 1464 individuals from the Rotterdam Study (RS). The participants whose exomes were sequenced and who were assessed for various cognitive traits were included in the analysis. To determine the association between DMD variants and cognitive ability, linear (mixed) modeling with adjustment for age, sex and education was used. Moreover, Sequence Kernel Association Test (SKAT) was used to test the overall association of the rare genetic variants present in the DMD with cognitive traits. Although no DMD variant surpassed the prespecified significance threshold (P<1 × 10(-4)), rs147546024:A>G showed strong association (β = 1.786, P-value = 2.56 × 10(-4)) with block-design test in the ERF study, while another variant rs1800273:G>A showed suggestive association (β = -0.465, P-value = 0.002) with Mini-Mental State Examination test in the RS. Both variants are highly conserved, although rs147546024:A>G is an intronic variant, whereas rs1800273:G>A is a missense variant in the DMD which has a predicted damaging effect on the protein. Further gene-based analysis of DMD revealed suggestive association (P-values = 0.087 and 0.074) with general cognitive ability in both cohorts. In conclusion, both single variant and gene-based analyses suggest the existence of variants in the DMD which may affect cognitive functioning in the general populations.

Genome-wide analysis of multi-ancestry cohorts identifies new loci influencing intraocular pressure and susceptibility to glaucoma

Elevated intraocular pressure (IOP) is an important risk factor in developing glaucoma and IOP variability may herald glaucomatous development or progression. We report the results of a genome-wide association study meta-analysis of 18 population cohorts from the International Glaucoma Genetics Consortium (IGGC), comprising 35,296 multiethnic participants for IOP. We confirm genetic association of known loci for IOP and primary open angle glaucoma (POAG) and identify four new IOP loci located on chromosome 3q25.31 within the FNDC3B gene (p=4.19×10⁻⁰⁸ for rs6445055), two on chromosome 9 (p=2.80×10⁻¹¹ for rs2472493 near ABCA1 and p=6.39×10⁻¹¹ for rs8176693 within ABO) and one on chromosome 11p11.2 (best p=1.04×10⁻¹¹ for rs747782). Separate meta-analyses of four independent POAG cohorts, totaling 4,284 cases and 95,560 controls, show that three of these IOP loci are also associated with POAG.

Genome-wide analyses of borderline personality features

The heritability of borderline personality (BP) features has been established in multiple twin and family studies. Using data from the borderline subscale of the Personality Assessment Inventory Borderline Features Scale (PAI-BOR) collected in two Dutch cohorts (N=7125), the Netherlands Twin Register and The Netherlands Study of Depression and Anxiety, we show that heritability of the PAI-BOR total score using genome-wide single-nucleotide polymorphism (SNPs) is estimated at 23%, and that the genetic variance is substantially higher in affect instability items compared with the other three subscales of the PAI-BOR (42.7% vs non-significant estimates for self-harm, negative relations and identity problems). We present results from a first genome-wide association study of BP features, which shows a promising signal on chromosome 5 corresponding to SERINC5, a protein involved in myelination. Reduced myelination has been suggested as possibly having a role in the development of psychiatric disorders characterized by lack of social interaction. The signal was confirmed in a third independent Dutch cohort drawn from the Erasmus Rucphen Family study (N=1301). Our analyses were complemented by investigating the heterogeneity that was implied by the differences in genetic variance components in the four subscales of the PAI-BOR. These analyses show that the association of SNPs tagging SERINC5 differs substantially across the 24 items of the PAI-BOR. Further, using reverse regression we showed that the effects were present only in subjects with higher scores on the PAI-BOR. Taken together, these results suggest that future genome-wide analyses can benefit substantially by taking into account the phenotypic and genetic heterogeneity of BP features.

Association analysis of bitter receptor genes in five isolated populations identifies a significant correlation between TAS2R43 variants and coffee liking

Coffee, one of the most popular beverages in the world, contains many different physiologically active compounds with a potential impact on people’s health. Despite the recent attention given to the genetic basis of its consumption, very little has been done in understanding genes influencing coffee preference among different individuals. Given its markedly bitter taste, we decided to verify if bitter receptor genes (TAS2Rs) variants affect coffee liking. In this light, 4066 people from different parts of Europe and Central Asia filled in a field questionnaire on coffee liking. They have been consequently recruited and included in the study. Eighty-eight SNPs covering the 25 TAS2R genes were selected from the available imputed ones and used to run association analysis for coffee liking. A significant association was detected with three SNP: one synonymous and two functional variants (W35S and H212R) on the TAS2R43 gene. Both variants have been shown to greatly reduce in vitro protein activity. Surprisingly the wild type allele, which corresponds to the functional form of the protein, is associated to higher liking of coffee. Since the hTAS2R43 receptor is sensible to caffeine, we verified if the detected variants produced differences in caffeine bitter perception on a subsample of people coming from the FVG cohort. We found a significant association between differences in caffeine perception and the H212R variant but not with the W35S, which suggests that the effect of the TAS2R43 gene on coffee liking is mediated by caffeine and in particular by the H212R variant.

No other significant association was found with other TAS2R genes. In conclusion, the present study opens new perspectives in the understanding of coffee liking. Further studies are needed to clarify the role of the TAS2R43 gene in coffee hedonics and to identify which other genes and pathways are involved in its genetics.

Heritability of dietary food intake patterns

The quality and quantity of food intake affect body weight, but little is known about the genetics of such human dietary intake patterns in relation to the genetics of BMI. We aimed to estimate the heritability of dietary intake patterns and genetic correlation with BMI in participants of the Erasmus Rucphen Family study. The study included 1,690 individuals (42 % men; age range, 19-92), of whom 41.4 % were overweight and 15.9 % were obese. Self-report questionnaires were used to assess the number of days (0-7) on which participants consumed vegetables, fruit, fruit juice, fish, unhealthy snacks, fastfood, and soft drinks. Principal component analysis was applied to examine the correlations between the questionnaire items and to generate dietary intake pattern scores. Heritability and the shared genetic and shared non-genetic (environmental) correlations were estimated using the family structure of the cohort. Principal component analysis suggested that the questionnaire items could be grouped in a healthy and unhealthy dietary intake pattern, explaining 22 and 18 % of the phenotypic variance, respectively. The dietary intake patterns had a heritability of 0.32 for the healthy and 0.27 for the unhealthy pattern. Genetic correlations between the dietary intake patterns and BMI were not significant, but we found a significant environmental correlation between the unhealthy dietary intake pattern and BMI. Specific dietary intake patterns are associated with the risk of obesity and are heritable traits. The genetic factors that determine specific dietary intake patterns do not significantly overlap with the genetic factors that determine BMI.

Nine loci for ocular axial length identified through genome-wide association studies, including shared loci with refractive error

Refractive errors are common eye disorders of public health importance worldwide. Ocular axial length (AL) is the major determinant of refraction and thus of myopia and hyperopia. We conducted a meta-analysis of genome-wide association studies for AL, combining 12,531 Europeans and 8,216 Asians. We identified eight genome-wide significant loci for AL (RSPO1, C3orf26, LAMA2, GJD2, ZNRF3, CD55, MIP, and ALPPL2) and confirmed one previously reported AL locus (ZC3H11B). Of the nine loci, five (LAMA2, GJD2, CD55, ALPPL2, and ZC3H11B) were associated with refraction in 18 independent cohorts (n = 23,591). Differential gene expression was observed for these loci in minus-lens-induced myopia mouse experiments and human ocular tissues. Two of the AL genes, RSPO1 and ZNRF3, are involved in Wnt signaling, a pathway playing a major role in the regulation of eyeball size. This study provides evidence of shared genes between AL and refraction, but importantly also suggests that these traits may have unique pathways.

Region-based association analysis of human quantitative traits in related individuals

Regional-based association analysis instead of individual testing of each SNP was introduced in genome-wide association studies to increase the power of gene mapping, especially for rare genetic variants. For regional association tests, the kernel machine-based regression approach was recently proposed as a more powerful alternative to collapsing-based methods. However, the vast majority of existing algorithms and software for the kernel machine-based regression are applicable only to unrelated samples. In this paper, we present a new method for the kernel machine-based regression association analysis of quantitative traits in samples of related individuals. The method is based on the GRAMMAR+ transformation of phenotypes of related individuals, followed by use of existing kernel machine-based regression software for unrelated samples. We compared the performance of kernel-based association analysis on the material of the Genetic Analysis Workshop 17 family sample and real human data by using our transformation, the original untransformed trait, and environmental residuals. We demonstrated that only the GRAMMAR+ transformation produced type I errors close to the nominal value and that this method had the highest empirical power. The new method can be applied to analysis of related samples by using existing software for kernel-based association analysis developed for unrelated samples.

Meta-analysis of telomere length in 19,713 subjects reveals high heritability, stronger maternal inheritance and a paternal age effect

Telomere length (TL) has been associated with aging and mortality, but individual differences are also influenced by genetic factors, with previous studies reporting heritability estimates ranging from 34 to 82%. Here we investigate the heritability, mode of inheritance and the influence of parental age at birth on TL in six large, independent cohort studies with a total of 19,713 participants. The meta-analysis estimate of TL heritability was 0.70 (95% CI 0.64-0.76) and is based on a pattern of results that is highly similar for twins and other family members. We observed a stronger mother-offspring (r=0.42; P-value=3.60 × 10(-61)) than father-offspring correlation (r=0.33; P-value=7.01 × 10(-5)), and a significant positive association with paternal age at offspring birth (β=0.005; P-value=7.01 × 10(-5)). Interestingly, a significant and quite substantial correlation in TL between spouses (r=0.25; P-value=2.82 × 10(-30)) was seen, which appeared stronger in older spouse pairs (mean age ≥55 years; r=0.31; P-value=4.27 × 10(-23)) than in younger pairs (mean age<55 years; r=0.20; P-value=3.24 × 10(-10)). In summary, we find a high and very consistent heritability estimate for TL, evidence for a maternal inheritance component and a positive association with paternal age.

Refining genome-wide linkage intervals using a meta-analysis of genome-wide association studies identifies loci influencing personality dimensions

Personality traits are complex phenotypes related to psychosomatic health. Individually, various gene finding methods have not achieved much success in finding genetic variants associated with personality traits. We performed a meta-analysis of four genome-wide linkage scans (N=6149 subjects) of five basic personality traits assessed with the NEO Five-Factor Inventory. We compared the significant regions from the meta-analysis of linkage scans with the results of a meta-analysis of genome-wide association studies (GWAS) (N∼17 000). We found significant evidence of linkage of neuroticism to chromosome 3p14 (rs1490265, LOD=4.67) and to chromosome 19q13 (rs628604, LOD=3.55); of extraversion to 14q32 (ATGG002, LOD=3.3); and of agreeableness to 3p25 (rs709160, LOD=3.67) and to two adjacent regions on chromosome 15, including 15q13 (rs970408, LOD=4.07) and 15q14 (rs1055356, LOD=3.52) in the individual scans. In the meta-analysis, we found strong evidence of linkage of extraversion to 4q34, 9q34, 10q24 and 11q22, openness to 2p25, 3q26, 9p21, 11q24, 15q26 and 19q13 and agreeableness to 4q34 and 19p13. Significant evidence of association in the GWAS was detected between openness and rs677035 at 11q24 (P-value=2.6 × 10(-06), KCNJ1). The findings of our linkage meta-analysis and those of the GWAS suggest that 11q24 is a susceptible locus for openness, with KCNJ1 as the possible candidate gene.

Rapid variance components-based method for whole-genome association analysis

The variance component tests used in genome-wide association studies (GWAS) including large sample sizes become computationally exhaustive when the number of genetic markers is over a few hundred thousand. We present an extremely fast variance components-based two-step method, GRAMMAR-Gamma, developed as an analytical approximation within a framework of the score test approach. Using simulated and real human GWAS data sets, we show that this method provides unbiased estimates of the SNP effect and has a power close to that of the likelihood ratio test-based method. The computational complexity of our method is close to its theoretical minimum, that is, to the complexity of the analysis that ignores genetic structure. The running time of our method linearly depends on sample size, whereas this dependency is quadratic for other existing methods. Simulations suggest that GRAMMAR-Gamma may be used for association testing in whole-genome resequencing studies of large human cohorts.

Meta-analysis of genome-wide association studies for personality

Personality can be thought of as a set of characteristics that influence people's thoughts, feelings and behavior across a variety of settings. Variation in personality is predictive of many outcomes in life, including mental health. Here we report on a meta-analysis of genome-wide association (GWA) data for personality in 10 discovery samples (17,375 adults) and five in silico replication samples (3294 adults). All participants were of European ancestry. Personality scores for Neuroticism, Extraversion, Openness to Experience, Agreeableness and Conscientiousness were based on the NEO Five-Factor Inventory. Genotype data of ≈ 2.4M single-nucleotide polymorphisms (SNPs; directly typed and imputed using HapMap data) were available. In the discovery samples, classical association analyses were performed under an additive model followed by meta-analysis using the weighted inverse variance method. Results showed genome-wide significance for Openness to Experience near the RASA1 gene on 5q14.3 (rs1477268 and rs2032794, P=2.8 × 10(-8) and 3.1 × 10(-8)) and for Conscientiousness in the brain-expressed KATNAL2 gene on 18q21.1 (rs2576037, P=4.9 × 10(-8)). We further conducted a gene-based test that confirmed the association of KATNAL2 to Conscientiousness. In silico replication did not, however, show significant associations of the top SNPs with Openness and Conscientiousness, although the direction of effect of the KATNAL2 SNP on Conscientiousness was consistent in all replication samples. Larger scale GWA studies and alternative approaches are required for confirmation of KATNAL2 as a novel gene affecting Conscientiousness.

Dissecting the genetic make-up of North-East Sardinia using a large set of haploid and autosomal markers

Sardinia has been used for genetic studies because of its historical isolation, genetic homogeneity and increased prevalence of certain rare diseases. Controversy remains concerning the genetic substructure and the extent of genetic homogeneity, which has implications for the design of genome-wide association studies (GWAS). We revisited this issue by examining the genetic make-up of a sample from North-East Sardinia using a dense set of autosomal, Y chromosome and mitochondrial markers to assess the potential of the sample for GWAS and fine mapping studies. We genotyped individuals for 500K single-nucleotide polymorphisms, Y chromosome markers and sequenced the mitochondrial hypervariable (HVI-HVII) regions. We identified major haplogroups and compared these with other populations. We estimated linkage disequilibrium (LD) and haplotype diversity across autosomal markers, and compared these with other populations. Our results show that within Sardinia there is no major population substructure and thus it can be considered a genetically homogenous population. We did not find substantial differences in the extent of LD in Sardinians compared with other populations. However, we showed that at least 9% of genomic regions in Sardinians differed in LD structure, which is helpful for identifying functional variants using fine mapping. We concluded that Sardinia is a powerful setting for genetic studies including GWAS and other mapping approaches.

The genetic association between personality and major depression or bipolar disorder. A polygenic score analysis using genome-wide association data

The relationship between major depressive disorder (MDD) and bipolar disorder (BD) remains controversial. Previous research has reported differences and similarities in risk factors for MDD and BD, such as predisposing personality traits. For example, high neuroticism is related to both disorders, whereas openness to experience is specific for BD. This study examined the genetic association between personality and MDD and BD by applying polygenic scores for neuroticism, extraversion, openness to experience, agreeableness and conscientiousness to both disorders. Polygenic scores reflect the weighted sum of multiple single-nucleotide polymorphism alleles associated with the trait for an individual and were based on a meta-analysis of genome-wide association studies for personality traits including 13,835 subjects. Polygenic scores were tested for MDD in the combined Genetic Association Information Network (GAIN-MDD) and MDD2000+ samples (N=8921) and for BD in the combined Systematic Treatment Enhancement Program for Bipolar Disorder and Wellcome Trust Case-Control Consortium samples (N=6329) using logistic regression analyses. At the phenotypic level, personality dimensions were associated with MDD and BD. Polygenic neuroticism scores were significantly positively associated with MDD, whereas polygenic extraversion scores were significantly positively associated with BD. The explained variance of MDD and BD, ∼0.1%, was highly comparable to the variance explained by the polygenic personality scores in the corresponding personality traits themselves (between 0.1 and 0.4%). This indicates that the proportions of variance explained in mood disorders are at the upper limit of what could have been expected. This study suggests shared genetic risk factors for neuroticism and MDD on the one hand and for extraversion and BD on the other.

Insulin-resistance and metabolic syndrome are related to executive function in women in a large family-based study

While type 2 diabetes is well-known to be associated with poorer cognitive performance, few studies have reported on the association of metabolic syndrome (MetS) and contributing factors, such as insulin-resistance (HOMA-IR), low adiponectin-, and high C-reactive protein (CRP)- levels. We studied whether these factors are related to cognitive function and which of the MetS components are independently associated. The study was embedded in an ongoing family-based cohort study in a Dutch population. All participants underwent physical examinations, biomedical measurements, and neuropsychological testing. Linear regression models were used to determine the association between MetS, HOMA-IR, adiponectin levels, CRP, and cognitive test scores. Cross-sectional analyses were performed in 1,898 subjects (mean age 48 years, 43% men). People with MetS had significantly higher HOMA-IR scores, lower adiponectin levels, and higher CRP levels. MetS and high HOMA-IR were associated with poorer executive function in women (P = 0.03 and P = 0.009). MetS and HOMA-IR are associated with poorer executive function in women.

A genome-wide association study of optic disc parameters

The optic nerve head is involved in many ophthalmic disorders, including common diseases such as myopia and open-angle glaucoma. Two of the most important parameters are the size of the optic disc area and the vertical cup-disc ratio (VCDR). Both are highly heritable but genetically largely undetermined. We performed a meta-analysis of genome-wide association (GWA) data to identify genetic variants associated with optic disc area and VCDR. The gene discovery included 7,360 unrelated individuals from the population-based Rotterdam Study I and Rotterdam Study II cohorts. These cohorts revealed two genome-wide significant loci for optic disc area, rs1192415 on chromosome 1p22 (p = 6.72×10⁻¹⁹) within 117 kb of the CDC7 gene and rs1900004 on chromosome 10q21.3-q22.1 (p = 2.67×10⁻³³) within 10 kb of the ATOH7 gene. They revealed two genome-wide significant loci for VCDR, rs1063192 on chromosome 9p21 (p = 6.15×10⁻¹¹) in the CDKN2B gene and rs10483727 on chromosome 14q22.3-q23 (p = 2.93×10⁻¹⁰) within 40 kbp of the SIX1 gene. Findings were replicated in two independent Dutch cohorts (Rotterdam Study III and Erasmus Rucphen Family study; N = 3,612), and the TwinsUK cohort (N = 843). Meta-analysis with the replication cohorts confirmed the four loci and revealed a third locus at 16q12.1 associated with optic disc area, and four other loci at 11q13, 13q13, 17q23 (borderline significant), and 22q12.1 for VCDR. ATOH7 was also associated with VCDR independent of optic disc area. Three of the loci were marginally associated with open-angle glaucoma. The protein pathways in which the loci of optic disc area are involved overlap with those identified for VCDR, suggesting a common genetic origin.

Morphologic characteristics of the optic nerve head are involved in many ophthalmic diseases. Its size, called the optic disc area, is an important measure and has been associated with e.g. myopia and open-angle glaucoma (OAG). Another important and clinical parameter of the optic disc is the vertical cup-disc ratio (VCDR). Although studies have shown a high heritability of optic disc area and VCDR, its genetic determinants are still undetermined. We therefore conducted a genome-wide association (GWA) study on these quantitative traits, using data of over 11,000 Caucasian participants, and related the findings to myopia and OAG. We found evidence for association of three loci with optic disc area: CDC7/TGFBR3 region, ATOH7, and SALL1; and six with VCDR: CDKN2B, SIX1, SCYL1, CHEK2, ATOH7, and DCLK1; and additionally one borderline significant locus: BCAS3. None of the loci could be related to myopia. There was marginal evidence for association of ATOH7, CDKN2B, and SIX1 with OAG, which remains to be confirmed. The present study reveals new insights into the physiological development of the optic nerve and may shed light on the pathophysiological protein pathways leading to (neuro-) ophthalmologic diseases such as OAG.

PedHunter 2.0 and its usage to characterize the founder structure of the Old Order Amish of Lancaster County

BACKGROUND

Because they are a closed founder population, the Old Order Amish (OOA) of Lancaster County have been the subject of many medical genetics studies. We constructed four versions of Anabaptist Genealogy Database (AGDB) using three sources of genealogies and multiple updates. In addition, we developed PedHunter, a suite of query software that can solve pedigree-related problems automatically and systematically.

METHODS

We report on how we have used new features in PedHunter to quantify the number and expected genetic contribution of founders to the OOA. The queries and utility of PedHunter programs are illustrated by examples using AGDB in this paper. For example, we calculated the number of founders expected to be contributing genetic material to the present-day living OOA and estimated the mean relative founder representation for each founder. New features in PedHunter also include pedigree trimming and pedigree renumbering, which should prove useful for studying large pedigrees.

RESULTS

With PedHunter version 2.0 querying AGDB version 4.0, we identified 34,160 presumed living OOA individuals and connected them into a 14-generation pedigree descending from 554 founders (332 females and 222 males) after trimming. From the analysis of cumulative mean relative founder representation, 128 founders (78 females and 50 males) accounted for over 95% of the mean relative founder contribution among living OOA descendants.

DISCUSSION/CONCLUSIONS

The OOA are a closed founder population in which a modest number of founders account for the genetic variation present in the current OOA population. Improvements to the PedHunter software will be useful in future studies of both the OOA and other populations with large and computerized genealogies.

Genetic architecture of plasma adiponectin overlaps with the genetics of metabolic syndrome-related traits

OBJECTIVE

Adiponectin, a hormone secreted by adipose tissue, is of particular interest in metabolic syndrome, because it is inversely correlated with obesity and insulin sensitivity. However, it is not known to what extent the genetics of plasma adiponectin and the genetics of obesity and insulin sensitivity are interrelated. We aimed to evaluate the heritability of plasma adiponectin and its genetic correlation with the metabolic syndrome and metabolic syndrome-related traits and the association between these traits and 10 ADIPOQ single nucleotide polymorphisms (SNPs).

RESEARCH DESIGN AND METHODS

We made use of a family-based population, the Erasmus Rucphen Family study (1,258 women and 967 men). Heritability analysis was performed using a polygenic model. Genetic correlations were estimated using bivariate heritability analyses. Genetic association analysis was performed using a mixed model.

RESULTS

Plasma adiponectin showed a heritability of 55.1%. Genetic correlations between plasma adiponectin HDL cholesterol and plasma insulin ranged from 15 to 24% but were not significant for fasting glucose, triglycerides, blood pressure, homeostasis model assessment of insulin resistance (HOMA-IR), and C-reactive protein. A significant association with plasma adiponectin was found for ADIPOQ variants rs17300539 and rs182052. A nominally significant association was found with plasma insulin and HOMA-IR and ADIPOQ variant rs17300539 after adjustment for plasma adiponectin.

CONCLUSIONS

The significant genetic correlation between plasma adiponectin and HDL cholesterol and plasma insulin should be taken into account in the interpretation of genome-wide association studies. Association of ADIPOQ SNPs with plasma adiponectin was replicated, and we showed association between one ADIPOQ SNP and plasma insulin and HOMA-IR.

ProbABEL package for genome-wide association analysis of imputed data

Background

Over the last few years, genome-wide association (GWA) studies became a tool of choice for the identification of loci associated with complex traits. Currently, imputed single nucleotide polymorphisms (SNP) data are frequently used in GWA analyzes. Correct analysis of imputed data calls for the implementation of specific methods which take genotype imputation uncertainty into account.

Results

We developed the ProbABEL software package for the analysis of genome-wide imputed SNP data and quantitative, binary, and time-till-event outcomes under linear, logistic, and Cox proportional hazards models, respectively. For quantitative traits, the package also implements a fast two-step mixed model-based score test for association in samples with differential relationships, facilitating analysis in family-based studies, studies performed in human genetically isolated populations and outbred animal populations.

Conclusions

ProbABEL package provides fast efficient way to analyze imputed data in genome-wide context and will facilitate future identification of complex trait loci.

A meta-analysis of genome-wide data from five European isolates reveals an association of COL22A1, SYT1, and GABRR2 with serum creatinine level

BACKGROUND

Serum creatinine (S CR) is the most important biomarker for a quick and non-invasive assessment of kidney function in population-based surveys. A substantial proportion of the inter-individual variability in S CR level is explicable by genetic factors.

METHODS

We performed a meta-analysis of genome-wide association studies of S CR undertaken in five population isolates ('discovery cohorts'), all of which are part of the European Special Population Network (EUROSPAN) project. Genes showing the strongest evidence for an association with SCR (candidate loci) were replicated in two additional population-based samples ('replication cohorts').

RESULTS

After the discovery meta-analysis, 29 loci were selected for replication. Association between SCR level and polymorphisms in the collagen type XXII alpha 1 (COL22A1) gene, on chromosome 8, and in the synaptotagmin-1 (SYT1) gene, on chromosome 12, were successfully replicated in the replication cohorts (p value = 1.0 x 10(-6) and 1.7 x 10(-4), respectively). Evidence of association was also found for polymorphisms in a locus including the gamma-aminobutyric acid receptor rho-2 (GABRR2) gene and the ubiquitin-conjugating enzyme E2-J1 (UBE2J1) gene (replication p value = 3.6 x 10(-3)). Previously reported findings, associating glomerular filtration rate with SNPs in the uromodulin (UMOD) gene and in the schroom family member 3 (SCHROOM3) gene were also replicated.

CONCLUSIONS

While confirming earlier results, our study provides new insights in the understanding of the genetic basis of serum creatinine regulatory processes. In particular, the association with the genes SYT1 and GABRR2 corroborate previous findings that highlighted a possible role of the neurotransmitters GABAA receptors in the regulation of the glomerular basement membrane and a possible interaction between GABAA receptors and synaptotagmin-I at the podocyte level.

Comparison of participant information and informed consent forms of five European studies in genetic isolated populations

Family-based research in genetically isolated populations is an effective approach for identifying loci influencing variation in disease traits. In common with all studies in humans, those in genetically isolated populations need ethical approval; however, existing ethical frameworks may be inadequate to protect participant privacy and confidentiality and to address participants' information needs in such populations. Using the ethical-legal guidelines of the Council for International Organizations of Medical Sciences (CIOMS) as a template, we compared the participant information leaflets and consent forms of studies in five European genetically isolated populations to identify additional information that should be incorporated into information leaflets and consent forms to guarantee satisfactorily informed consent. We highlight the additional information that participants require on the research purpose and the reasons why their population was chosen; on the potential risks and benefits of participation; on the opportunities for benefit sharing; on privacy; on the withdrawal of consent and on the disclosure of genetic data. This research raises some important issues that should be addressed properly and identifies relevant types of information that should be incorporated into information leaflets for this type of study.

A powerful genome-wide feasible approach to detect parent-of-origin effects in studies of quantitative traits

There is currently a lot of interest in the role of genomic imprinting in mammalian development. Many human diseases, such as cancer, obesity, diabetes and behavioral traits, may be related to imprinted genes. When searching for genes related to complex disorders, the power of genome-wide association analysis can be improved by introducing parent-of-origin effects into the analyses. For quantitative traits, family-based TDT analysis has successfully implemented such an approach. Although attractive for several reasons, TDT-based tests are known to be less powerful than methods based on measured genotype approaches. In this study, we describe a fast, powerful method for detecting parent-of-origin effects in studies of quantitative traits using a measured genotype framework. First, for each locus studied, we estimate the probabilities of an allele's parental origin using multipoint haplotype reconstruction. Next, we introduce the parental origin of these alleles as a covariate in regression models during the second step of GRAMMAR, a fast approximation to the measured genotype approach. We show that, compared with a TDT-based analysis, our method has a higher power to detect a locus exhibiting a parent-of-origin effect. Moreover, our method is applicable to a wider range of data, including pedigree structures that are not very informative for TDT. The method gives no false positives in the absence of parent-of-origin effects, under both additive and dominant models. As this method is an extension of the rapid GRAMMAR analysis, it is fast enough to be suitable for genome-wide association scans.

Association between type 2 diabetes loci and measures of fatness

Background

Type 2 diabetes (T2D) is a metabolic disorder characterized by disturbances of carbohydrate, fat and protein metabolism and insulin resistance. The majority of T2D patients are obese and obesity by itself may be a cause of insulin resistance. Our aim was to evaluate whether the recently identified T2D risk alleles are associated with human measures of fatness as characterized with Dual Energy X-ray Absorptiometry (DEXA).

Methodology/Principal Findings

Genotypes and phenotypes of approximately 3,000 participants from cross-sectional ERF study were analyzed. Nine single nucleotide polymorphisms (SNPs) in CDKN2AB, CDKAL1, FTO, HHEX, IGF2BP2, KCNJ11, PPARG, SLC30A8 and TCF7L2 were genotyped. We used linear regression to study association between individual SNPs and the combined allelic risk score with body mass index (BMI), fat mass index (FMI), fat percentage (FAT), waist circumference (WC) and waist to hip ratio (WHR). Significant association was observed between rs8050136 (FTO) and BMI (p = 0.003), FMI (p = 0.007) and WC (p = 0.03); fat percentage was borderline significant (p = 0.053). No other SNPs alone or combined in a risk score demonstrated significant association to the measures of fatness.

Conclusions/Significance

From the recently identified T2D risk variants only the risk variant of the FTO gene (rs8050136) showed statistically significant association with BMI, FMI, and WC.

Heterozygous NTF4 mutations impairing neurotrophin-4 signaling in patients with primary open-angle glaucoma

Glaucoma, a main cause of blindness in the developed world, is characterized by progressive degeneration of retinal ganglion cells (RGCs), resulting in irreversible loss of vision. Although members of the neurotrophin gene family in various species are known to support the survival of numerous neuronal populations, including RGCs, it is less clear whether they are also required for survival and maintenance of adult neurons in humans. Here, we report seven different heterozygous mutations in the Neurotrophin-4 (NTF4) gene accounting for about 1.7% of primary open-angle glaucoma patients of European origin. Molecular modeling predicted a decreased affinity of neurotrophin 4 protein (NT-4) mutants with its specific tyrosine kinase receptor B (TrkB). Expression of recombinant NT-4 carrying the most frequent mutation was demonstrated to lead to decreased activation of TrkB. These findings suggest a pathway in the pathophysiology of glaucoma through loss of neurotrophic function and may eventually open the possibility of using ligands activating TrkB to prevent the progression of the disease.

A multiple splitting approach to linkage analysis in large pedigrees identifies a linkage to asthma on chromosome 12

Large genealogies are potentially very informative for linkage analysis. However, the software available for exact non-parametric multipoint linkage analysis is limited with respect to the complexity of the families it can handle. A solution is to split the large pedigrees into sub-families meeting complexity constraints. Different methods have been proposed to "best" split large genealogies. Here, we propose a new procedure in which linkage is performed on several carefully chosen sub-pedigree sets from the genealogy instead of using just a single sub-pedigree set. Our multiple splitting procedure capitalizes on the sensitivity of linkage results to family structure and has been designed to control computational feasibility and global type I error. We describe and apply this procedure to the extreme case of the highly complex Hutterite pedigree and use it to perform a genome-wide linkage analysis on asthma. The detection of a genome-wide significant linkage for asthma on chromosome 12q21 illustrates the potential of this multiple splitting approach.