A general approach for haplotype phasing across the full spectrum of relatedness

Genome-wide association study identifies distinct genetic contributions to prognosis and susceptibility in Crohn's disease

For most immune-mediated diseases, the main determinant of patient well-being is not the diagnosis itself but instead the course that the disease takes over time (prognosis). Prognosis may vary substantially between patients for reasons that are poorly understood. Familial studies support a genetic contribution to prognosis, but little evidence has been found for a proposed association between prognosis and the burden of susceptibility variants. To better characterize how genetic variation influences disease prognosis, we performed a within-cases genome-wide association study in two cohorts of patients with Crohn's disease. We identified four genome-wide significant loci, none of which showed any association with disease susceptibility. Conversely, the aggregated effect of all 170 disease susceptibility loci was not associated with disease prognosis. Together, these data suggest that the genetic contribution to prognosis in Crohn's disease is largely independent of the contribution to disease susceptibility and point to a biology of prognosis that could provide new therapeutic opportunities.

Recombinant Haplotypes Narrow the ARMS2/HTRA1 Association Signal for Age-Related Macular Degeneration

Age-related macular degeneration (AMD) is the leading cause of blindness in ageing societies, triggered by both environmental and genetic factors. The strongest genetic signal for AMD with odds ratios of up to 2.8 per adverse allele was found previously over a chromosomal region in 10q26 harboring two genes, ARMS2 and HTRA1, although with little knowledge as to which gene or genetic variation is functionally relevant to AMD pathology. In this study, we analyzed rare recombinant haplotypes in 16,144 AMD cases and 17,832 controls from the International AMD Genomics Consortium and identified variants in ARMS2 but not HTRA1 to exclusively carry the AMD risk with P-values between 1.0 × 10-773 and 6.7 × 10-5 This now allows prioritization of the gene of interest for subsequent functional studies.

Genome-Wide Association Study of Polymorphisms Predisposing to Bronchiolitis

Bronchiolitis is a major cause of hospitalization among infants. Severe bronchiolitis is associated with later asthma, suggesting a common genetic predisposition. Genetic background of bronchiolitis is not well characterized. To identify polymorphisms associated with bronchiolitis, we conducted a genome-wide association study (GWAS) in which 5,300,000 single nucleotide polymorphisms (SNPs) were tested for association in a Finnish–Swedish population of 217 children hospitalized for bronchiolitis and 778 controls. The most promising SNPs (n = 77) were genotyped in a Dutch replication population of 416 cases and 432 controls. Finally, we used a set of 202 Finnish bronchiolitis cases to further investigate candidate SNPs. We did not detect genome-wide significant associations, but several suggestive association signals (p < 10⁻⁵) were observed in the GWAS. In the replication population, three SNPs were nominally associated (p < 0.05). Of them, rs269094 was an expression quantitative trait locus (eQTL) for KCND3, previously shown to be associated with occupational asthma. In the additional set of Finnish cases, the association for another SNP (rs9591920) within a noncoding RNA locus was further strengthened. Our results provide a first genome-wide examination of the genetics underlying bronchiolitis. These preliminary findings require further validation in a larger sample size.

MixFit: Methodology for Computing Ancestry-Related Genetic Scores at the Individual Level and Its Application to the Estonian and Finnish Population Studies

Ancestry information at the individual level can be a valuable resource for personalized medicine, medical, demographical and history research, as well as for tracing back personal history. We report a new method for quantitatively determining personal genetic ancestry based on genome-wide data. Numerical ancestry component scores are assigned to individuals based on comparisons with reference populations. These comparisons are conducted with an existing analytical pipeline making use of genotype phasing, similarity matrix computation and our addition-multidimensional best fitting by MixFit. The method is demonstrated by studying Estonian and Finnish populations in geographical context. We show the main differences in the genetic composition of these otherwise close European populations and how they have influenced each other. The components of our analytical pipeline are freely available computer programs and scripts one of which was developed in house (available at: www.geenivaramu.ee/en/tools/mixfit).

A survey of sub-Saharan gene flow into the Mediterranean at risk loci for coronary artery disease

This study tries to find detectable signals of gene flow of Sub-Saharan origin into the Mediterranean in four genomic regions previously associated with coronary artery disease. A total of 366 single-nucleotide polymorphisms were genotyped in 772 individuals from 10 Mediterranean countries. Population structure analyses were performed, in which a noticeable Sub-Saharan component was found in the studied samples. The overall percentage of this Sub-Saharan component presents differences between the two Mediterranean coasts. D-statistics suggest possible Sub-Saharan introgression into one of the studied genomic regions (10q11). We also found differences in linkage disequilibrium patterns between the two Mediterranean coasts, possibly attributable to differential Sub-Saharan admixture. Our results confirm the potentially important role of human demographic history when performing epidemiological studies.

A genome wide association study suggests the association of muskelin with early onset bipolar disorder: Implications for a GABAergic epileptogenic neurogenesis model

BACKGROUND

Although multiple genes have been implicated in bipolar disorder (BD), they explain only a small proportion of its heritability. Identifying additional BD risk variants may be impaired by phenotypic heterogeneity, which is usually not taken into account in genome-wide association studies (GWAS). BD with early age at onset is a more homogeneous familial form of the disorder associated with greater symptom severity.

METHODS

We conducted a GWAS of early-onset BD (onset of mania/hypomania ≤19 years old) in a discovery sample of 419 cases and 1034 controls and a replication sample of 181 cases and 777 controls. These two samples were meta-analyzed, followed by replication of one signal in a third independent sample of 141 cases and 746 controls.

RESULTS

No single nucleotide polymorphism (SNP) associations were genome-wide significant in the discovery sample. Of the top 15 SNPs in the discovery analysis, rs114034759 in the muskelin (MKLN1) gene was nominally significant in the replication analysis, and was among the top associations in the meta-analysis (p=2.63E-06, OR=1.9). In the third sample, this SNP was again associated with early-onset BD (p=0.036, OR=1.6). Gene expression analysis showed that the rs114034759 risk allele is associated with decreased hippocampal MKLN1 expression.

LIMITATIONS

The sample sizes of the early-onset BD subgroups were relatively small.

CONCLUSIONS

Our results suggest MKLN1 is associated with early-onset BD. MKLN1 regulates cellular trafficking of GABA-A receptors, which is involved in synaptic transmission and plasticity, and is implicated in the mechanism of action of a group of antiepileptic mood stabilizers. These results therefore indicate that GABAergic neurotransmission may be implicated in early-onset BD. We propose that an increase in GABA-A receptors in the hippocampus in BD patients due to lower MKLN1 expression might increase the excitability during the GABA-excited early phase of young neurons, leading to an increased risk of developing a manic/hypomanic episode. Further studies are needed to test this model.

SCOPA and META-SCOPA: software for the analysis and aggregation of genome-wide association studies of multiple correlated phenotypes

Background

Genome-wide association studies (GWAS) of single nucleotide polymorphisms (SNPs) have been successful in identifying loci contributing genetic effects to a wide range of complex human diseases and quantitative traits. The traditional approach to GWAS analysis is to consider each phenotype separately, despite the fact that many diseases and quantitative traits are correlated with each other, and often measured in the same sample of individuals. Multivariate analyses of correlated phenotypes have been demonstrated, by simulation, to increase power to detect association with SNPs, and thus may enable improved detection of novel loci contributing to diseases and quantitative traits.

Results

We have developed the SCOPA software to enable GWAS analysis of multiple correlated phenotypes. The software implements “reverse regression” methodology, which treats the genotype of an individual at a SNP as the outcome and the phenotypes as predictors in a general linear model. SCOPA can be applied to quantitative traits and categorical phenotypes, and can accommodate imputed genotypes under a dosage model. The accompanying META-SCOPA software enables meta-analysis of association summary statistics from SCOPA across GWAS. Application of SCOPA to two GWAS of high-and low-density lipoprotein cholesterol, triglycerides and body mass index, and subsequent meta-analysis with META-SCOPA, highlighted stronger association signals than univariate phenotype analysis at established lipid and obesity loci. The META-SCOPA meta-analysis also revealed a novel signal of association at genome-wide significance for triglycerides mapping to GPC5 (lead SNP rs71427535, p = 1.1x10⁻⁸), which has not been reported in previous large-scale GWAS of lipid traits.

Conclusions

The SCOPA and META-SCOPA software enable discovery and dissection of multiple phenotype association signals through implementation of a powerful reverse regression approach.

Analysis of Whole Exome Sequencing with Cardiometabolic Traits Using Family-Based Linkage and Association in the IRAS Family Study

Family-based methods are a potentially powerful tool to identify trait-defining genetic variants in extended families, particularly when used to complement conventional association analysis. We utilized two-point linkage analysis and single variant association analysis to evaluate whole exome sequencing (WES) data from 1205 Hispanic Americans (78 families) from the Insulin Resistance Atherosclerosis Family Study. WES identified 211,612 variants above the minor allele frequency threshold of ≥0.005. These variants were tested for linkage and/or association with 50 cardiometabolic traits after quality control checks. Two-point linkage analysis yielded 10,580,600 logarithm of the odds (LOD) scores with 1148 LOD scores ≥3, 183 LOD scores ≥4, and 29 LOD scores ≥5. The maximal novel LOD score was 5.50 for rs2289043:T>C, in UNC5C with subcutaneous adipose tissue volume. Association analysis identified 13 variants attaining genome-wide significance (P < 5 × 10^-08 ), with the strongest association between rs651821:C>T in APOA5 and triglyceride levels (P  =  3.67 × 10^-10 ). Overall, there was a 5.2-fold increase in the number of informative variants detected by WES compared to exome chip analysis in this population, nearly 30% of which were novel variants relative to the Database of Single Nucleotide Polymorphisms (dbSNP) build 138. Thus, integration of results from two-point linkage and single-variant association analysis from WES data enabled identification of novel signals potentially contributing to cardiometabolic traits.

ERAP1 association with ankylosing spondylitis is attributable to common genotypes rather than rare haplotype combinations

We investigated the proposal that ankylosing spondylitis (AS) is associated with unusual ERAP1 genotypes. ERAP1 haplotypes were constructed for 213 AS cases and 46 rheumatoid arthritis controls using family data. Haplotypes were generated from five common ERAP1 single nucleotide polymorphisms (SNPs)-rs2287987 (M349V), rs30187 (K528R), rs10050860 (D575N), rs17482078 (R725Q), and rs27044 (Q730E). Haplotype frequencies were compared using Fisher's exact test. ERAP1 haplotypes imputed from the International Genetics of AS Consortium (IGAS) Immunochip study were also studied. In the family study, we identified only four common ERAP1 haplotypes ("VRNQE," "MKDRQ," "MRDRE," and "MKDRE") in both AS cases and controls apart from two rare (<0.5%) previously unreported haplotypes. There were no examples of the unusual ERAP1 haplotype combination ("*001/*005") previously reported by others in 53% of AS cases. As expected, K528-bearing haplotypes were increased in the AS family study (AS 43% vs. control 35%), due particularly to an increase in the MKDRQ haplotype (AS 35% vs. control 25%, P = 0.01). This trend was replicated in the imputed Immunochip data for the two K528-bearing haplotypes MKDRQ (AS 33% vs. controls 27%, P = 1.2 × 10^-24) and MKDRE (AS 8% vs. controls 7%, P = 0.004). The ERAP1 association with AS is therefore predominantly attributable to common ERAP1 haplotypes and haplotype combinations.

Genetic Variation in the SLC8A1 Calcium Signaling Pathway Is Associated With Susceptibility to Kawasaki Disease and Coronary Artery Abnormalities

Background—

Kawasaki disease (KD) is an acute pediatric vasculitis in which host genetics influence both susceptibility to KD and the formation of coronary artery aneurysms. Variants discovered by genome-wide association studies and linkage studies only partially explain the influence of genetics on KD susceptibility.

Methods and Results—

To search for additional functional genetic variation, we performed pathway and gene stability analysis on a genome-wide association study data set. Pathway analysis using European genome-wide association study data identified 100 significantly associated pathways ( P <5×10 − 4 ). Gene stability selection identified 116 single nucleotide polymorphisms in 26 genes that were responsible for driving the pathway associations, and gene ontology analysis demonstrated enrichment for calcium transport ( P =1.05×10 − 4 ). Three single nucleotide polymorphisms in solute carrier family 8, member 1 ( SLC8A1 ), a sodium/calcium exchanger encoding NCX1, were validated in an independent Japanese genome-wide association study data set (meta-analysis P =0.0001). Patients homozygous for the A (risk) allele of rs13017968 had higher rates of coronary artery abnormalities ( P =0.029). NCX1, the protein encoded by SLC8A1 , was expressed in spindle-shaped and inflammatory cells in the aneurysm wall. Increased intracellular calcium mobilization was observed in B cell lines from healthy controls carrying the risk allele.

Conclusions—

Pathway-based association analysis followed by gene stability selection proved to be a valuable tool for identifying risk alleles in a rare disease with complex genetics. The role of SLC8A1 polymorphisms in altering calcium flux in cells that mediate coronary artery damage in KD suggests that this pathway may be a therapeutic target and supports the study of calcineurin inhibitors in acute KD.

rehh 2.0: a reimplementation of the R package rehh to detect positive selection from haplotype structure

Identifying genomic regions with unusually high local haplotype homozygosity represents a powerful strategy to characterize candidate genes responding to natural or artificial positive selection. To that end, statistics measuring the extent of haplotype homozygosity within (e.g. EHH, iHS) and between (Rsb or XP-EHH) populations have been proposed in the literature. The rehh package for r was previously developed to facilitate genome-wide scans of selection, based on the analysis of long-range haplotypes. However, its performance was not sufficient to cope with the growing size of available data sets. Here, we propose a major upgrade of the rehh package, which includes an improved processing of the input files, a faster algorithm to enumerate haplotypes, as well as multithreading. As illustrated with the analysis of large human haplotype data sets, these improvements decrease the computation time by more than one order of magnitude. This new version of rehh will thus allow performing iHS-, Rsb- or XP-EHH-based scans on large data sets. The package rehh 2.0 is available from the CRAN repository (http://cran.r-project.org/web/packages/rehh/index.html) together with help files and a detailed manual.

A Pedigree-Based Map of Recombination in the Domestic Dog Genome

Meiotic recombination in mammals has been shown to largely cluster into hotspots, which are targeted by the chromatin modifier PRDM9. The canid family, including wolves and dogs, has undergone a series of disrupting mutations in this gene, rendering PRDM9 inactive. Given the importance of PRDM9, it is of great interest to learn how its absence in the dog genome affects patterns of recombination placement. We have used genotypes from domestic dog pedigrees to generate sex-specific genetic maps of recombination in this species. On a broad scale, we find that placement of recombination events in dogs is consistent with that in mice and apes, in that the majority of recombination occurs toward the telomeres in males, while female crossing over is more frequent and evenly spread along chromosomes. It has been previously suggested that dog recombination is more uniform in distribution than that of humans; however, we found that recombination in dogs is less uniform than in humans. We examined the distribution of recombination within the genome, and found that recombination is elevated immediately upstream of the transcription start site and around CpG islands, in agreement with previous studies, but that this effect is stronger in male dogs. We also found evidence for positive crossover interference influencing the spacing between recombination events in dogs, as has been observed in other species including humans and mice. Overall our data suggests that dogs have similar broad scale properties of recombination to humans, while fine scale recombination is similar to other species lacking PRDM9.

Gene expression elucidates functional impact of polygenic risk for schizophrenia

Over 100 genetic loci harbor schizophrenia-associated variants, yet how these variants confer liability is uncertain. The CommonMind Consortium sequenced RNA from dorsolateral prefrontal cortex of people with schizophrenia (N = 258) and control subjects (N = 279), creating a resource of gene expression and its genetic regulation. Using this resource, ∼20% of schizophrenia loci have variants that could contribute to altered gene expression and liability. In five loci, only a single gene was involved: FURIN, TSNARE1, CNTN4, CLCN3 or SNAP91. Altering expression of FURIN, TSNARE1 or CNTN4 changed neurodevelopment in zebrafish; knockdown of FURIN in human neural progenitor cells yielded abnormal migration. Of 693 genes showing significant case-versus-control differential expression, their fold changes were ≤ 1.33, and an independent cohort yielded similar results. Gene co-expression implicates a network relevant for schizophrenia. Our findings show that schizophrenia is polygenic and highlight the utility of this resource for mechanistic interpretations of genetic liability for brain diseases.

Pleiotropic Mechanisms Indicated for Sex Differences in Autism

Sexual dimorphism in common disease is pervasive, including a dramatic male preponderance in autism spectrum disorders (ASDs). Potential genetic explanations include a liability threshold model requiring increased polymorphism risk in females, sex-limited X-chromosome contribution, gene-environment interaction driven by differences in hormonal milieu, risk influenced by genes sex-differentially expressed in early brain development, or contribution from general mechanisms of sexual dimorphism shared with secondary sex characteristics. Utilizing a large single nucleotide polymorphism (SNP) dataset, we identify distinct sex-specific genome-wide significant loci. We investigate genetic hypotheses and find no evidence for increased genetic risk load in females, but evidence for sex heterogeneity on the X chromosome, and contribution of sex-heterogeneous SNPs for anthropometric traits to ASD risk. Thus, our results support pleiotropy between secondary sex characteristic determination and ASDs, providing a biological basis for sex differences in ASDs and implicating non brain-limited mechanisms.

Imputing rare variants in families using a two-stage approach

BACKGROUND

Recent focus on studying rare variants makes imputation accuracy of rare variants an important issue. Many approaches have been proposed to increase imputation accuracy among rare variants, from reference panel selection to combinations of existing methods to multistage analyses. We aimed to bring the strengths of these new approaches together with our proposed two-stage imputation for family data.

METHODS

Our imputation methods were tested on the region from 46.75Mb to 49.25Mb on chromosome 3. We did quality control based on the proportion of missing genotypes per variant and individual, leaving 495 individuals with 761 genome-wide association studies (GWAS) variants only, 45 with 14,077 sequence variants only, and 419 with both GWAS and sequencing data. All data were prephased using SHAPEIT2 with a duo hidden Markov model algorithm prior to performing imputation. Imputations were performed 100 times, each time masking the sequence data for 1 individual and imputing it from the GWAS data. We used well-imputed genotypes, defined as a probability of greater than 0.9, above 2 different minor allele frequency cutoffs-0.01 and 0.05-from Impute2 as input for Merlin, and compared these results to Impute2 and Merlin separately. The imputed results were evaluated using correlation measurement and the imputation quality score.

RESULTS

Our method improved imputation accuracy, measured by imputation quality score, for variants with minor allele frequency between 0.01 and 0.40, but failed to improve accuracy for variants with minor allele frequency less than 0.01 when we used a minor allele frequency cutoff of 0.01 for the Impute2 results. In contrast, our 2-stage approach with a minor allele frequency cutoff of 0.05 performed the worst of all methods for variants with minor allele frequency between 0.01 and 0.40.

CONCLUSIONS

This method gave promising results, but may be further improved by changing the inclusion criteria of Impute2 variants. More analyses are needed on a larger region with different inclusion thresholds to assess the accuracy of this approach.

Increasing Generality and Power of Rare-Variant Tests by Utilizing Extended Pedigrees

Recently, multiple studies have performed whole-exome or whole-genome sequencing to identify groups of rare variants associated with complex traits and diseases. They have primarily utilized case-control study designs that often require thousands of individuals to reach acceptable statistical power. Family-based studies can be more powerful because a rare variant can be enriched in an extended pedigree and segregate with the phenotype. Although many methods have been proposed for using family data to discover rare variants involved in a disease, a majority of them focus on a specific pedigree structure and are designed to analyze either binary or continuously measured outcomes. In this article, we propose RareIBD, a general and powerful approach to identifying rare variants involved in disease susceptibility. Our method can be applied to large extended families of arbitrary structure, including pedigrees with only affected individuals. The method accommodates both binary and quantitative traits. A series of simulation experiments suggest that RareIBD is a powerful test that outperforms existing approaches. In addition, our method accounts for individuals in top generations, which are not usually genotyped in extended families. In contrast to available statistical tests, RareIBD generates accurate p values even when genetic data from these individuals are missing. We applied RareIBD, as well as other methods, to two extended family datasets generated by different genotyping technologies and representing different ethnicities. The analysis of real data confirmed that RareIBD is the only method that properly controls type I error.

Single nucleotide polymorphisms in an intergenic chromosome 2q region associated with tissue factor pathway inhibitor plasma levels and venous thromboembolism

Essentials Tissue factor pathway inhibitor (TFPI) regulates the blood coagulation cascade. We replicated previously reported linkage of TFPI plasma levels to the chromosome 2q region. The putative causal locus, rs62187992, was associated with TFPI plasma levels and thrombosis. rs62187992 was marginally associated with TFPI expression in human aortic endothelial cells. Click to hear Ann Gil's presentation on new insights into thrombin activatable fibrinolysis inhibitor SUMMARY: Background Tissue factor pathway inhibitor (TFPI) regulates fibrin clot formation, and low TFPI plasma levels increase the risk of arterial thromboembolism and venous thromboembolism (VTE). TFPI plasma levels are also heritable, and a previous linkage scan implicated the chromosome 2q region, but no specific genes. Objectives To replicate the finding of the linkage region in an independent sample, and to identify the causal locus. Methods We first performed a linkage analysis of microsatellite markers and TFPI plasma levels in 251 individuals from the F5L Family Study, and replicated the finding of the linkage peak on chromosome 2q (LOD = 3.06). We next defined a follow-up region that included 112 603 single nucleotide polymorphisms (SNPs) under the linkage peak, and meta-analyzed associations between these SNPs and TFPI plasma levels across the F5L Family Study and the Marseille Thrombosis Association (MARTHA) Study, a study of 1033 unrelated VTE patients. SNPs with false discovery rate q-values of < 0.10 were tested for association with TFPI plasma levels in 892 patients with coronary artery disease in the AtheroGene Study. Results and Conclusions One SNP, rs62187992, was associated with TFPI plasma levels in all three samples (β = + 0.14 and P = 4.23 × 10-6 combined; β = + 0.16 and P = 0.02 in the F5L Family Study; β = + 0.13 and P = 6.3 × 10-4 in the MARTHA Study; β = + 0.17 and P = 0.03 in the AtheroGene Study), and contributed to the linkage peak in the F5L Family Study. rs62187992 was also associated with clinical VTE (odds ratio 0.90, P = 0.03) in the INVENT Consortium of > 7000 cases and their controls, and was marginally associated with TFPI expression (β = + 0.19, P = 0.08) in human aortic endothelial cells, a primary site of TFPI synthesis. The biological mechanisms underlying these associations remain to be elucidated.

Genetic Complexity of Crohn's Disease in Two Large Ashkenazi Jewish Families

BACKGROUND & AIMS

Crohn's disease (CD) is a highly heritable disease that is particularly common in the Ashkenazi Jewish population. We studied 2 large Ashkenazi Jewish families with a high prevalence of CD in an attempt to identify novel genetic risk variants.

METHODS

Ashkenazi Jewish patients with CD and a positive family history were recruited from the University College London Hospital. We used genome-wide, single-nucleotide polymorphism data to assess the burden of common CD-associated risk variants and for linkage analysis. Exome sequencing was performed and rare variants that were predicted to be deleterious and were observed at a high frequency in cases were prioritized. We undertook within-family association analysis after imputation and assessed candidate variants for evidence of association with CD in an independent cohort of Ashkenazi Jewish individuals. We examined the effects of a variant in DUOX2 on hydrogen peroxide production in HEK293 cells.

RESULTS

We identified 2 families (1 with >800 members and 1 with >200 members) containing 54 and 26 cases of CD or colitis, respectively. Both families had a significant enrichment of previously described common CD-associated risk variants. No genome-wide significant linkage was observed. Exome sequencing identified candidate variants, including a missense mutation in DUOX2 that impaired its function and a frameshift mutation in CSF2RB that was associated with CD in an independent cohort of Ashkenazi Jewish individuals.

CONCLUSIONS

In a study of 2 large Ashkenazi Jewish with multiple cases of CD, we found the genetic basis of the disease to be complex, with a role for common and rare genetic variants. We identified a frameshift mutation in CSF2RB that was replicated in an independent cohort. These findings show the value of family studies and the importance of the innate immune system in the pathogenesis of CD.

Identification of type 2 diabetes subgroups through topological analysis of patient similarity

Type 2 diabetes (T2D) is a heterogeneous complex disease affecting more than 29 million Americans alone with a rising prevalence trending toward steady increases in the coming decades. Thus, there is a pressing clinical need to improve early prevention and clinical management of T2D and its complications. Clinicians have understood that patients who carry the T2D diagnosis have a variety of phenotypes and susceptibilities to diabetes-related complications. We used a precision medicine approach to characterize the complexity of T2D patient populations based on high-dimensional electronic medical records (EMRs) and genotype data from 11,210 individuals. We successfully identified three distinct subgroups of T2D from topology-based patient-patient networks. Subtype 1 was characterized by T2D complications diabetic nephropathy and diabetic retinopathy; subtype 2 was enriched for cancer malignancy and cardiovascular diseases; and subtype 3 was associated most strongly with cardiovascular diseases, neurological diseases, allergies, and HIV infections. We performed a genetic association analysis of the emergent T2D subtypes to identify subtype-specific genetic markers and identified 1279, 1227, and 1338 single-nucleotide polymorphisms (SNPs) that mapped to 425, 322, and 437 unique genes specific to subtypes 1, 2, and 3, respectively. By assessing the human disease-SNP association for each subtype, the enriched phenotypes and biological functions at the gene level for each subtype matched with the disease comorbidities and clinical differences that we identified through EMRs. Our approach demonstrates the utility of applying the precision medicine paradigm in T2D and the promise of extending the approach to the study of other complex, multifactorial diseases.

Ancestral Origins and Genetic History of Tibetan Highlanders

The origin of Tibetans remains one of the most contentious puzzles in history, anthropology, and genetics. Analyses of deeply sequenced (30×-60×) genomes of 38 Tibetan highlanders and 39 Han Chinese lowlanders, together with available data on archaic and modern humans, allow us to comprehensively characterize the ancestral makeup of Tibetans and uncover their origins. Non-modern human sequences compose ∼6% of the Tibetan gene pool and form unique haplotypes in some genomic regions, where Denisovan-like, Neanderthal-like, ancient-Siberian-like, and unknown ancestries are entangled and elevated. The shared ancestry of Tibetan-enriched sequences dates back to ∼62,000-38,000 years ago, predating the Last Glacial Maximum (LGM) and representing early colonization of the plateau. Nonetheless, most of the Tibetan gene pool is of modern human origin and diverged from that of Han Chinese ∼15,000 to ∼9,000 years ago, which can be largely attributed to post-LGM arrivals. Analysis of ∼200 contemporary populations showed that Tibetans share ancestry with populations from East Asia (∼82%), Central Asia and Siberia (∼11%), South Asia (∼6%), and western Eurasia and Oceania (∼1%). Our results support that Tibetans arose from a mixture of multiple ancestral gene pools but that their origins are much more complicated and ancient than previously suspected. We provide compelling evidence of the co-existence of Paleolithic and Neolithic ancestries in the Tibetan gene pool, indicating a genetic continuity between pre-historical highland-foragers and present-day Tibetans. In particular, highly differentiated sequences harbored in highlanders' genomes were most likely inherited from pre-LGM settlers of multiple ancestral origins (SUNDer) and maintained in high frequency by natural selection.

Genome-wide scans reveal variants at EDAR predominantly affecting hair straightness in Han Chinese and Uyghur populations

Hair straightness/curliness is one of the most conspicuous features of human variation and is particularly diverse among populations. A recent genome-wide scan found common variants in the Trichohyalin (TCHH) gene that are associated with hair straightness in Europeans, but different genes might affect this phenotype in other populations. By sampling 2899 Han Chinese, we performed the first genome-wide scan of hair straightness in East Asians, and found EDAR (rs3827760) as the predominant gene (P = 4.67 × 10^-16), accounting for 3.66 % of the total variance. The candidate gene approach did not find further significant associations, suggesting that hair straightness may be affected by a large number of genes with subtle effects. Notably, genetic variants associated with hair straightness in Europeans are generally low in frequency in Han Chinese, and vice versa. To evaluate the relative contribution of these variants, we performed a second genome-wide scan in 709 samples from the Uyghur, an admixed population with both eastern and western Eurasian ancestries. In Uyghurs, both EDAR (rs3827760: P = 1.92 × 10^-12) and TCHH (rs11803731: P = 1.46 × 10^-3) are associated with hair straightness, but EDAR (OR 0.415) has a greater effect than TCHH (OR 0.575). We found no significant interaction between EDAR and TCHH (P = 0.645), suggesting that these two genes affect hair straightness through different mechanisms. Furthermore, haplotype analysis indicates that TCHH is not subject to selection. While EDAR is under strong selection in East Asia, it does not appear to be subject to selection after the admixture in Uyghurs. These suggest that hair straightness is unlikely a trait under selection.

Genome-wide association meta-analysis in Chinese and European individuals identifies ten new loci associated with systemic lupus erythematosus

Systemic lupus erythematosus (SLE; OMIM 152700) is a genetically complex autoimmune disease. Genome-wide association studies (GWASs) have identified more than 50 loci as robustly associated with the disease in single ancestries, but genome-wide transancestral studies have not been conducted. We combined three GWAS data sets from Chinese (1,659 cases and 3,398 controls) and European (4,036 cases and 6,959 controls) populations. A meta-analysis of these studies showed that over half of the published SLE genetic associations are present in both populations. A replication study in Chinese (3,043 cases and 5,074 controls) and European (2,643 cases and 9,032 controls) subjects found ten previously unreported SLE loci. Our study provides further evidence that the majority of genetic risk polymorphisms for SLE are contained within the same regions across both populations. Furthermore, a comparison of risk allele frequencies and genetic risk scores suggested that the increased prevalence of SLE in non-Europeans (including Asians) has a genetic basis.

A thrifty variant in CREBRF strongly influences body mass index in Samoans

Samoans are a unique founder population with a high prevalence of obesity, making them well suited for identifying new genetic contributors to obesity. We conducted a genome-wide association study (GWAS) in 3,072 Samoans, discovered a variant, rs12513649, strongly associated with body mass index (BMI) (P = 5.3 × 10(-14)), and replicated the association in 2,102 additional Samoans (P = 1.2 × 10(-9)). Targeted sequencing identified a strongly associated missense variant, rs373863828 (p.Arg457Gln), in CREBRF (meta P = 1.4 × 10(-20)). Although this variant is extremely rare in other populations, it is common in Samoans (frequency of 0.259), with an effect size much larger than that of any other known common BMI risk variant (1.36-1.45 kg/m(2) per copy of the risk-associated allele). In comparison to wild-type CREBRF, the Arg457Gln variant when overexpressed selectively decreased energy use and increased fat storage in an adipocyte cell model. These data, in combination with evidence of positive selection of the allele encoding p.Arg457Gln, support a 'thrifty' variant hypothesis as a factor in human obesity.

High Level of Fasting Plasma Proenkephalin-A Predicts Deterioration of Kidney Function and Incidence of CKD

High levels of proenkephalin-A (pro-ENK) have been associated with decreased eGFR in an acute setting. Here, we examined whether pro-ENK levels predict CKD and decline of renal function in a prospective cohort of 2568 participants without CKD (eGFR>60 ml/min per 1.73 m²) at baseline. During a mean follow-up of 16.6 years, 31.7% of participants developed CKD. Participants with baseline pro-ENK levels in the highest tertile had significantly greater yearly mean decline of eGFR (P_trend<0.001) and rise of cystatin C (P_trend=0.01) and creatinine (P_trend<0.001) levels. Furthermore, compared with participants in the lowest tertile, participants in the highest tertile of baseline pro-ENK concentration had increased CKD incidence (odds ratio, 1.51; 95% confidence interval, 1.18 to 1.94) when adjusted for multiple factors. Adding pro-ENK to a model of conventional risk factors in net reclassification improvement analysis resulted in reclassification of 14.14% of participants. Genome-wide association analysis in 4150 participants of the same cohort revealed the strongest association of pro-ENK levels with rs1012178 near the PENK gene, where the minor T-allele associated with a 0.057 pmol/L higher pro-ENK level per allele (P=4.67x10^-21). Furthermore, the T-allele associated with a 19% increased risk of CKD per allele (P=0.03) and a significant decrease in the instrumental variable estimator for eGFR (P<0.01) in a Mendelian randomization analysis. In conclusion, circulating plasma pro-ENK level predicts incident CKD and may aid in identifying subjects in need of primary preventive regimens. Additionally, the Mendelian randomization analysis suggests a causal relationship between pro-ENK level and deterioration of kidney function over time.

Biallelic Mutations in TMEM126B Cause Severe Complex I Deficiency with a Variable Clinical Phenotype

Complex I deficiency is the most common biochemical phenotype observed in individuals with mitochondrial disease. With 44 structural subunits and over 10 assembly factors, it is unsurprising that complex I deficiency is associated with clinical and genetic heterogeneity. Massively parallel sequencing (MPS) technologies including custom, targeted gene panels or unbiased whole-exome sequencing (WES) are hugely powerful in identifying the underlying genetic defect in a clinical diagnostic setting, yet many individuals remain without a genetic diagnosis. These individuals might harbor mutations in poorly understood or uncharacterized genes, and their diagnosis relies upon characterization of these orphan genes. Complexome profiling recently identified TMEM126B as a component of the mitochondrial complex I assembly complex alongside proteins ACAD9, ECSIT, NDUFAF1, and TIMMDC1. Here, we describe the clinical, biochemical, and molecular findings in six cases of mitochondrial disease from four unrelated families affected by biallelic (c.635G>T [p.Gly212Val] and/or c.401delA [p.Asn134Ilefs^∗2]) TMEM126B variants. We provide functional evidence to support the pathogenicity of these TMEM126B variants, including evidence of founder effects for both variants, and establish defects within this gene as a cause of complex I deficiency in association with either pure myopathy in adulthood or, in one individual, a severe multisystem presentation (chronic renal failure and cardiomyopathy) in infancy. Functional experimentation including viral rescue and complexome profiling of subject cell lines has confirmed TMEM126B as the tenth complex I assembly factor associated with human disease and validates the importance of both genome-wide sequencing and proteomic approaches in characterizing disease-associated genes whose physiological roles have been previously undetermined.

PADRE: Pedigree-Aware Distant-Relationship Estimation

Accurate estimation of shared ancestry is an important component of many genetic studies; current prediction tools accurately estimate pairwise genetic relationships up to the ninth degree. Pedigree-aware distant-relationship estimation (PADRE) combines relationship likelihoods generated by estimation of recent shared ancestry (ERSA) with likelihoods from family networks reconstructed by pedigree reconstruction and identification of a maximum unrelated set (PRIMUS), improving the power to detect distant relationships between pedigrees. Using PADRE, we estimated relationships from simulated pedigrees and three extended pedigrees, correctly predicting 20% more fourth- through ninth-degree simulated relationships than when using ERSA alone. By leveraging pedigree information, PADRE can even identify genealogical relationships between individuals who are genetically unrelated. For example, although 95% of 13(th)-degree relatives are genetically unrelated, in simulations, PADRE correctly predicted 50% of 13(th)-degree relationships to within one degree of relatedness. The improvement in prediction accuracy was consistent between simulated and actual pedigrees. We also applied PADRE to the HapMap3 CEU samples and report new cryptic relationships and validation of previously described relationships between families. PADRE greatly expands the range of relationships that can be estimated by using genetic data in pedigrees.

Fast and accurate long-range phasing in a UK Biobank cohort

Recent work has leveraged the extensive genotyping of the Icelandic population to perform long-range phasing (LRP), enabling accurate imputation and association analysis of rare variants in target samples typed on genotyping arrays. Here we develop a fast and accurate LRP method, Eagle, that extends this paradigm to populations with much smaller proportions of genotyped samples by harnessing long (>4-cM) identical-by-descent (IBD) tracts shared among distantly related individuals. We applied Eagle to N ≈ 150,000 samples (0.2% of the British population) from the UK Biobank, and we determined that it is 1-2 orders of magnitude faster than existing methods while achieving similar or better phasing accuracy (switch error rate ≈ 0.3%, corresponding to perfect phase in a majority of 10-Mb segments). We also observed that, when used within an imputation pipeline, Eagle prephasing improved downstream imputation accuracy in comparison to prephasing in batches using existing methods, as necessary to achieve comparable computational cost.

Haplotype estimation for biobank-scale data sets

The UK Biobank (UKB) has recently released genotypes on 152,328 individuals together with extensive phenotypic and lifestyle information. We present a new phasing method SHAPEIT3 that can handle such biobank scale datasets and results in switch error rates as low as ~0.3%. The method exhibits O(NlogN) scaling in sample size (N), enabling fast and accurate phasing of even larger cohorts.

A genome-wide association scan implicates DCHS2, RUNX2, GLI3, PAX1 and EDAR in human facial variation

We report a genome-wide association scan for facial features in ∼6,000 Latin Americans. We evaluated 14 traits on an ordinal scale and found significant association (P values<5 × 10⁻⁸) at single-nucleotide polymorphisms (SNPs) in four genomic regions for three nose-related traits: columella inclination (4q31), nose bridge breadth (6p21) and nose wing breadth (7p13 and 20p11). In a subsample of ∼3,000 individuals we obtained quantitative traits related to 9 of the ordinal phenotypes and, also, a measure of nasion position. Quantitative analyses confirmed the ordinal-based associations, identified SNPs in 2q12 associated to chin protrusion, and replicated the reported association of nasion position with SNPs in PAX3. Strongest association in 2q12, 4q31, 6p21 and 7p13 was observed for SNPs in the EDAR, DCHS2, RUNX2 and GLI3 genes, respectively. Associated SNPs in 20p11 extend to PAX1. Consistent with the effect of EDAR on chin protrusion, we documented alterations of mandible length in mice with modified Edar funtion.

Humans show great diversity in facial appearance and this variation is highly heritable. Here, Andres Ruiz-Linares and colleagues examined facial features in admixed Latin Americans and identify genome-wide associations for 14 facial traits, including four gene loci (RUNX2, GLI3, DCHS2 and PAX1) influencing nose morphology.

Multi-ethnic genome-wide association study identifies novel locus for type 2 diabetes susceptibility

Genome-wide association studies (GWAS) have traditionally been undertaken in homogeneous populations from the same ancestry group. However, with the increasing availability of GWAS in large-scale multi-ethnic cohorts, we have evaluated a framework for detecting association of genetic variants with complex traits, allowing for population structure, and developed a powerful test of heterogeneity in allelic effects between ancestry groups. We have applied the methodology to identify and characterise loci associated with susceptibility to type 2 diabetes (T2D) using GWAS data from the Resource for Genetic Epidemiology on Adult Health and Aging, a large multi-ethnic population-based cohort, created for investigating the genetic and environmental basis of age-related diseases. We identified a novel locus for T2D susceptibility at genome-wide significance (P<5 × 10⁻⁸) that maps to TOMM40-APOE, a region previously implicated in lipid metabolism and Alzheimer's disease. We have also confirmed previous reports that single-nucleotide polymorphisms at the TCF7L2 locus demonstrate the greatest extent of heterogeneity in allelic effects between ethnic groups, with the lowest risk observed in populations of East Asian ancestry.

Multiallelic copy number variation in the complement component 4A (C4A) gene is associated with late-stage age-related macular degeneration (AMD)

BACKGROUND

Age-related macular degeneration (AMD) is the leading cause of vision loss in Western societies with a strong genetic component. Candidate gene studies as well as genome-wide association studies strongly implicated genetic variations in complement genes to be involved in disease risk. So far, no association of AMD with complement component 4 (C4) was reported probably due to the complex nature of the C4 locus on chromosome 6.

METHODS

We used multiplex ligation-dependent probe amplification (MLPA) to determine the copy number of the C4 gene as well as of both relevant isoforms, C4A and C4B, and assessed their association with AMD using logistic regression models.

RESULTS

Here, we report on the analysis of 2645 individuals (1536 probands and 1109 unaffected controls), across three different centers, for multiallelic copy number variation (CNV) at the C4 locus. We find strong statistical significance for association of increased copy number of C4A (OR 0.81 (0.73; 0.89);P = 4.4 × 10(-5)), with the effect most pronounced in individuals over 78 years (OR 0.67 (0.55; 0.81)) and females (OR 0.77 (0.68; 0.87)). Furthermore, this association is independent of known AMD-associated risk variants in the nearby CFB/C2 locus, particularly in females and in individuals over 78 years.

CONCLUSIONS

Our data strengthen the notion that complement dysregulation plays a crucial role in AMD etiology, an important finding for early intervention strategies and future therapeutics. In addition, for the first time, we provide evidence that multiallelic CNVs are associated with AMD pathology.

Transcript Isoform Variation Associated with Cytosine Modification in Human Lymphoblastoid Cell Lines

Cytosine modification on DNA is variable among individuals, which could correlate with gene expression variation. The effect of cytosine modification on interindividual transcript isoform variation (TIV), however, remains unclear. In this study, we assessed the extent of cytosine modification-specific TIV in lymphoblastoid cell lines (LCLs) derived from unrelated individuals of European and African descent. Our study detected cytosine modification-specific TIVs for 17% of the analyzed genes at a 5% false discovery rate. Forty-five percent of the TIV-associated cytosine modifications correlated with the overall gene expression levels as well, with the corresponding CpG sites overrepresented in transcript initiation sites, transcription factor binding sites, and distinct histone modification peaks, suggesting that alternative isoform transcription underlies the TIVs. Our analysis also revealed 33% of the TIV-associated cytosine modifications that affected specific exons, with the corresponding CpG sites overrepresented in exon/intron junctions, splicing branching points, and transcript termination sites, implying that the TIVs are attributable to alternative splicing or transcription termination. Genetic and epigenetic regulation of TIV shared target preference but exerted independent effects on 61% of the common exon targets. Cytosine modification-specific TIVs detected from LCLs were differentially enriched in those detected from various tissues in The Cancer Genome Atlas, indicating their developmental dependency. Genes containing cytosine modification-specific TIVs were enriched in pathways of cancers and metabolic disorders. Our study demonstrated a prominent effect of cytosine modification variation on the transcript isoform spectrum over gross transcript abundance and revealed epigenetic contributions to diseases that were mediated through cytosine modification-specific TIV.

High rates of phasing errors in highly polymorphic species with low levels of linkage disequilibrium

Short read sequencing of diploid individuals does not permit the direct inference of the sequence on each of the two homologous chromosomes. Although various phasing software packages exist, they were primarily tailored for and tested on human data, which differ from other species in factors that influence phasing, such as SNP density, amounts of linkage disequilibrium (LD) and sample sizes. Despite becoming increasingly popular for other species, the reliability of phasing in non-human data has not been evaluated to a sufficient extent. We scrutinized the phasing accuracy for Drosophila melanogaster, a species with high polymorphism levels and reduced LD relative to humans. We phased two D. melanogaster populations and compared the results to the known haplotypes. The performance increased with size of the reference panel and was highest when the reference panel and phased individuals were from the same population. Full genomic SNP data and inclusion of sequence read information also improved phasing. Despite humans and Drosophila having similar switch error rates between polymorphic sites, the distances between switch errors were much shorter in Drosophila with only fragments <300-1500 bp being correctly phased with ≥95% confidence. This suggests that the higher SNP density cannot compensate for the higher recombination rate in D. melanogaster. Furthermore, we show that populations that have gone through demographic events such as bottlenecks can be phased with higher accuracy. Our results highlight that statistically phased data are particularly error prone in species with large population sizes or populations lacking suitable reference panels.

Inferring Population Size History from Large Samples of Genome-Wide Molecular Data - An Approximate Bayesian Computation Approach

Inferring the ancestral dynamics of effective population size is a long-standing question in population genetics, which can now be tackled much more accurately thanks to the massive genomic data available in many species. Several promising methods that take advantage of whole-genome sequences have been recently developed in this context. However, they can only be applied to rather small samples, which limits their ability to estimate recent population size history. Besides, they can be very sensitive to sequencing or phasing errors. Here we introduce a new approximate Bayesian computation approach named PopSizeABC that allows estimating the evolution of the effective population size through time, using a large sample of complete genomes. This sample is summarized using the folded allele frequency spectrum and the average zygotic linkage disequilibrium at different bins of physical distance, two classes of statistics that are widely used in population genetics and can be easily computed from unphased and unpolarized SNP data. Our approach provides accurate estimations of past population sizes, from the very first generations before present back to the expected time to the most recent common ancestor of the sample, as shown by simulations under a wide range of demographic scenarios. When applied to samples of 15 or 25 complete genomes in four cattle breeds (Angus, Fleckvieh, Holstein and Jersey), PopSizeABC revealed a series of population declines, related to historical events such as domestication or modern breed creation. We further highlight that our approach is robust to sequencing errors, provided summary statistics are computed from SNPs with common alleles.

Molecular data sampled from extant individuals contains considerable information about their demographic history. In particular, one classical question in population genetics is to reconstruct past population size changes from such data. Relating these changes to various climatic, geological or anthropogenic events allows characterizing the main factors driving genetic diversity and can have major outcomes for conservation. Until recently, mostly very simple histories, including one or two population size changes, could be estimated from genetic data. This has changed with the sequencing of entire genomes in many species, and several methods allow now inferring complex histories consisting of several tens of population size changes. However, analyzing entire genomes, while accounting for recombination, remains a statistical and numerical challenge. These methods, therefore, can only be applied to small samples with a few diploid genomes. We overcome this limitation by using an approximate estimation approach, where observed genomes are summarized using a small number of statistics related to allele frequencies and linkage disequilibrium. In contrast to previous approaches, we show that our method allows us to reconstruct also the most recent part (the last 100 generations) of the population size history. As an illustration, we apply it to large samples of whole-genome sequences in four cattle breeds.

The phenotypic legacy of admixture between modern humans and Neandertals

Many modern human genomes retain DNA inherited from interbreeding with archaic hominins, such as Neandertals, yet the influence of this admixture on human traits is largely unknown. We analyzed the contribution of common Neandertal variants to over 1000 electronic health record (EHR)-derived phenotypes in ~28,000 adults of European ancestry. We discovered and replicated associations of Neandertal alleles with neurological, psychiatric, immunological, and dermatological phenotypes. Neandertal alleles together explained a significant fraction of the variation in risk for depression and skin lesions resulting from sun exposure (actinic keratosis), and individual Neandertal alleles were significantly associated with specific human phenotypes, including hypercoagulation and tobacco use. Our results establish that archaic admixture influences disease risk in modern humans, provide hypotheses about the effects of hundreds of Neandertal haplotypes, and demonstrate the utility of EHR data in evolutionary analyses.

A genome-wide association scan in admixed Latin Americans identifies loci influencing facial and scalp hair features

We report a genome-wide association scan in over 6,000 Latin Americans for features of scalp hair (shape, colour, greying, balding) and facial hair (beard thickness, monobrow, eyebrow thickness). We found 18 signals of association reaching genome-wide significance (P values 5 × 10(-8) to 3 × 10(-119)), including 10 novel associations. These include novel loci for scalp hair shape and balding, and the first reported loci for hair greying, monobrow, eyebrow and beard thickness. A newly identified locus influencing hair shape includes a Q30R substitution in the Protease Serine S1 family member 53 (PRSS53). We demonstrate that this enzyme is highly expressed in the hair follicle, especially the inner root sheath, and that the Q30R substitution affects enzyme processing and secretion. The genome regions associated with hair features are enriched for signals of selection, consistent with proposals regarding the evolution of human hair.

Polymorphisms in HLA Class II Genes Are Associated With Susceptibility to Staphylococcus aureus Infection in a White Population

BACKGROUND

Staphylococcus aureus can cause life-threatening infections. Human susceptibility to S. aureus infection may be influenced by host genetic variation.

METHODS

A genome-wide association study (GWAS) in a large health plan-based cohort included biologic specimens from 4701 culture-confirmed S. aureus cases and 45 344 matched controls; 584 535 single-nucleotide polymorphisms (SNPs) were genotyped on an array specific to individuals of European ancestry. Coverage was increased by imputation of >25 million common SNPs, using the 1000 Genomes Reference panel. In addition, human leukocyte antigen (HLA) serotypes were also imputed.

RESULTS

Logistic regression analysis, performed under the assumption of an additive genetic model, revealed several imputed SNPs (eg, rs115231074: odds ratio [OR], 1.22 [P = 1.3 × 10(-10)]; rs35079132: OR, 1.24 [P = 3.8 × 10(-8)]) achieving genome-wide significance on chromosome 6 in the HLA class II region. One adjacent genotyped SNP was nearly genome-wide significant (rs4321864: OR, 1.13; P = 8.8 × 10(-8)). These polymorphisms are located near the genes encoding HLA-DRA and HLA-DRB1. Results of further logistic regression analysis, in which the most significant GWAS SNPs were conditioned on HLA-DRB1*04 serotype, showed additional support for the strength of association between HLA class II genetic variants and S. aureus infection.

CONCLUSIONS

Our study results are the first reported evidence of human genetic susceptibility to S. aureus infection.

Genetic variants near MLST8 and DHX57 affect the epigenetic age of the cerebellum

DNA methylation (DNAm) levels lend themselves for defining an epigenetic biomarker of aging known as the ‘epigenetic clock'. Our genome-wide association study (GWAS) of cerebellar epigenetic age acceleration identifies five significant (P<5.0 × 10⁻⁸) SNPs in two loci: 2p22.1 (inside gene DHX57) and 16p13.3 near gene MLST8 (a subunit of mTOR complex 1 and 2). We find that the SNP in 16p13.3 has a cis-acting effect on the expression levels of MLST8 (P=6.9 × 10⁻¹⁸) in most brain regions. In cerebellar samples, the SNP in 2p22.1 has a cis-effect on DHX57 (P=4.4 × 10⁻⁵). Gene sets found by our GWAS analysis of cerebellar age acceleration exhibit significant overlap with those of Alzheimer's disease (P=4.4 × 10⁻¹⁵), age-related macular degeneration (P=6.4 × 10⁻⁶), and Parkinson's disease (P=2.6 × 10⁻⁴). Overall, our results demonstrate the utility of a new paradigm for understanding aging and age-related diseases: it will be fruitful to use epigenetic tissue age as endophenotype in GWAS.

This genome-wide association study identifies five significant SNPs in two loci which are associated with the epigenetic age of post-mortem cerebellar tissue according to a DNA methylation based biomarker of human aging.

Epigenetic age of the pre-frontal cortex is associated with neuritic plaques, amyloid load, and Alzheimer's disease related cognitive functioning

There is an urgent need to develop molecular biomarkers of brain age in order to advance our understanding of age related neurodegeneration. Recently, we developed a highly accurate epigenetic biomarker of tissue age (known as epigenetic clock) which is based on DNA methylation levels. Here we use n=700 dorsolateral prefrontal cortex (DLPFC) samples from Caucasian subjects of the Religious Order Study and the Rush Memory and Aging Project to examine the association between epigenetic age and Alzheimer's disease (AD) related cognitive decline, and AD related neuropathological markers. Epigenetic age acceleration of DLPFC is correlated with several neuropathological measurements including diffuse plaques (r=0.12, p=0.0015), neuritic plaques (r=0.11, p=0.0036), and amyloid load (r=0.091, p=0.016). Further, it is associated with a decline in global cognitive functioning (β=-0.500, p=0.009), episodic memory (β=-0.411, p=0.009) and working memory (β=-0.405, p=0.011) among individuals with AD. The neuropathological markers may mediate the association between epigenetic age and cognitive decline. Genetic complex trait analysis (GCTA) revealed that epigenetic age acceleration is heritable (h2=0.41) and has significant genetic correlations with diffuse plaques (r=0.24, p=0.010) and possibly working memory (r=-0.35, p=0.065). Overall, these results suggest that the epigenetic clock may lend itself as a molecular biomarker of brain age.

Oxidative phosphorylation and lacunar stroke: Genome-wide enrichment analysis of common variants

OBJECTIVE

We investigated whether oxidative phosphorylation (OXPHOS) abnormalities were associated with lacunar stroke, hypothesizing that these would be more strongly associated in patients with multiple lacunar infarcts and leukoaraiosis (LA).

METHODS

In 1,012 MRI-confirmed lacunar stroke cases and 964 age-matched controls recruited from general practice surgeries, we investigated associations between common genetic variants within the OXPHOS pathway and lacunar stroke using a permutation-based enrichment approach. Cases were phenotyped using MRI into those with multiple infarcts or LA (MLI/LA) and those with isolated lacunar infarcts (ILI) based on the number of subcortical infarcts and degree of LA, using the Fazekas grading. Using gene-level association statistics, we tested for enrichment of genes in the OXPHOS pathway with all lacunar stroke and the 2 subtypes.

RESULTS

There was a specific association with strong evidence of enrichment in the top 1% of genes in the MLI/LA (subtype p = 0.0017) but not in the ILI subtype (p = 1). Genes in the top percentile for the all lacunar stroke analysis were not significantly enriched (p = 0.07).

CONCLUSIONS

Our results implicate the OXPHOS pathway in the pathogenesis of lacunar stroke, and show the association is specific to patients with the MLI/LA subtype. They show that MRI-based subtyping of lacunar stroke can provide insights into disease pathophysiology, and imply that different radiologic subtypes of lacunar stroke subtypes have distinct underlying pathophysiologic processes.

Serum bilirubin concentration is modified by UGT1A1 haplotypes and influences risk of type-2 diabetes in the Norfolk Island genetic isolate

Background

Located in the Pacific Ocean between Australia and New Zealand, the unique population isolate of Norfolk Island has been shown to exhibit increased prevalence of metabolic disorders (type-2 diabetes, cardiovascular disease) compared to mainland Australia. We investigated this well-established genetic isolate, utilising its unique genomic structure to increase the ability to detect related genetic markers. A pedigree-based genome-wide association study of 16 routinely collected blood-based clinical traits in 382 Norfolk Island individuals was performed.

Results

A striking association peak was located at chromosome 2q37.1 for both total bilirubin and direct bilirubin, with 29 SNPs reaching statistical significance (P < 1.84 × 10⁻⁷). Strong linkage disequilibrium was observed across a 200 kb region spanning the UDP-glucuronosyltransferase family, including UGT1A1, an enzyme known to metabolise bilirubin. Given the epidemiological literature suggesting negative association between CVD-risk and serum bilirubin we further explored potential associations using stepwise multivariate regression, revealing significant association between direct bilirubin concentration and type-2 diabetes risk. In the Norfolk Island cohort increased direct bilirubin was associated with a 28 % reduction in type-2 diabetes risk (OR: 0.72, 95 % CI: 0.57-0.91, P = 0.005). When adjusted for genotypic effects the overall model was validated, with the adjusted model predicting a 30 % reduction in type-2 diabetes risk with increasing direct bilirubin concentrations (OR: 0.70, 95 % CI: 0.53-0.89, P = 0.0001).

Conclusions

In summary, a pedigree-based GWAS of blood-based clinical traits in the Norfolk Island population has identified variants within the UDPGT family directly associated with serum bilirubin levels, which is in turn implicated with reduced risk of developing type-2 diabetes within this population.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0291-z) contains supplementary material, which is available to authorized users.

Conjunctival fibrosis and the innate barriers to Chlamydia trachomatis intracellular infection: a genome wide association study

Chlamydia trachomatis causes both trachoma and sexually transmitted infections. These diseases have similar pathology and potentially similar genetic predisposing factors. We aimed to identify polymorphisms and pathways associated with pathological sequelae of ocular Chlamydia trachomatis infections in The Gambia. We report a discovery phase genome-wide association study (GWAS) of scarring trachoma (1090 cases, 1531 controls) that identified 27 SNPs with strong, but not genome-wide significant, association with disease (5 × 10⁻⁶ > P > 5 × 10⁻⁸). The most strongly associated SNP (rs111513399, P = 5.38 × 10⁻⁷) fell within a gene (PREX2) with homology to factors known to facilitate chlamydial entry to the host cell. Pathway analysis of GWAS data was significantly enriched for mitotic cell cycle processes (P = 0.001), the immune response (P = 0.00001) and for multiple cell surface receptor signalling pathways. New analyses of published transcriptome data sets from Gambia, Tanzania and Ethiopia also revealed that the same cell cycle and immune response pathways were enriched at the transcriptional level in various disease states. Although unconfirmed, the data suggest that genetic associations with chlamydial scarring disease may be focussed on processes relating to the immune response, the host cell cycle and cell surface receptor signalling.

The role of common genetic variation in educational attainment and income: evidence from the National Child Development Study

We investigated the role of common genetic variation in educational attainment and household income. We used data from 5,458 participants of the National Child Development Study to estimate: 1) the associations of rs9320913, rs11584700 and rs4851266 and socioeconomic position and educational phenotypes; and 2) the univariate chip-heritability of each phenotype, and the genetic correlation between each phenotype and educational attainment at age 16. The three SNPs were associated with most measures of educational attainment. Common genetic variation contributed to 6 of 14 socioeconomic background phenotypes, and 17 of 29 educational phenotypes. We found evidence of genetic correlations between educational attainment at age 16 and 4 of 14 social background and 8 of 28 educational phenotypes. This suggests common genetic variation contributes both to differences in educational attainment and its relationship with other phenotypes. However, we remain cautious that cryptic population structure, assortative mating, and dynastic effects may influence these associations.

Genome-Wide Association Study of Late-Onset Myasthenia Gravis: Confirmation of TNFRSF11A and Identification of ZBTB10 and Three Distinct HLA Associations

To investigate the genetics of late-onset myasthenia gravis (LOMG), we conducted a genome-wide association study imputation of>6 million single nucleotide polymorphisms (SNPs) in 532 LOMG cases (anti-acetylcholine receptor [AChR] antibody positive; onset age≥50 years) and 2,128 controls matched for sex and population substructure. The data confirm reported TNFRSF11A associations (rs4574025, P = 3.9 × 10^-7, odds ratio [OR] 1.42) and identify a novel candidate gene, ZBTB10, achieving genome-wide significance (rs6998967, P = 8.9 × 10^-10, OR 0.53). Several other SNPs showed suggestive significance including rs2476601 (P = 6.5 × 10^-6, OR 1.62) encoding the PTPN22 R620W variant noted in early-onset myasthenia gravis (EOMG) and other autoimmune diseases. In contrast, EOMG-associated SNPs in TNIP1 showed no association in LOMG, nor did other loci suggested for EOMG. Many SNPs within the major histocompatibility complex (MHC) region showed strong associations in LOMG, but with smaller effect sizes than in EOMG (highest OR ~2 versus ~6 in EOMG). Moreover, the strongest associations were in opposite directions from EOMG, including an OR of 0.54 for DQA1*05:01 in LOMG (P = 5.9 × 10^-12) versus 2.82 in EOMG (P = 3.86 × 10^-45). Association and conditioning studies for the MHC region showed three distinct and largely independent association peaks for LOMG corresponding to (a) MHC class II (highest attenuation when conditioning on DQA1), (b) HLA-A and (c) MHC class III SNPs. Conditioning studies of human leukocyte antigen (HLA) amino acid residues also suggest potential functional correlates. Together, these findings emphasize the value of subgrouping myasthenia gravis patients for clinical and basic investigations and imply distinct predisposing mechanisms in LOMG.

Auton A, Brooks LD, Durbin RM, Garrison EP, Kang HM, Korbel JO, Marchini JL, McCarthy S, McVean GA, Abecasis GR. A global reference for human genetic variation

The 1000 Genomes Project set out to provide a comprehensive description of common human genetic variation by applying whole-genome sequencing to a diverse set of individuals from multiple populations. Here we report completion of the project, having reconstructed the genomes of 2,504 individuals from 26 populations using a combination of low-coverage whole-genome sequencing, deep exome sequencing, and dense microarray genotyping. We characterized a broad spectrum of genetic variation, in total over 88 million variants (84.7 million single nucleotide polymorphisms (SNPs), 3.6 million short insertions/deletions (indels), and 60,000 structural variants), all phased onto high-quality haplotypes. This resource includes >99% of SNP variants with a frequency of >1% for a variety of ancestries. We describe the distribution of genetic variation across the global sample, and discuss the implications for common disease studies.

Results for the final phase of the 1000 Genomes Project are presented including whole-genome sequencing, targeted exome sequencing, and genotyping on high-density SNP arrays for 2,504 individuals across 26 populations, providing a global reference data set to support biomedical genetics.

The 1000 Genomes Project has sought to comprehensively catalogue human genetic variation across populations, providing a valuable public genomic resource. The data obtained so far have found applications ranging from association studies and fine mapping studies to the filtering of likely neutral variants in rare-disease cohorts. The authors now report on the final phase of the project, phase 3, which covers previously uncharacterized areas of human genetic diversity in terms of the populations sampled and categories of characterized variation. The sample now includes more than 2,500 individuals from 26 global populations, with low coverage whole-genome and deep exome sequencing, as well as dense microarray genotyping. They find that while most common variants are shared across populations, rarer variants are often restricted to closely related populations. The authors also demonstrate the use of the phase 3 dataset as a reference panel for imputation to improve the resolution in genetic association studies.

Coronary risk in relation to genetic variation in MEOX2 and TCF15 in a Flemish population

Background

In mice MEOX2/TCF15 heterodimers are highly expressed in heart endothelial cells and are involved in the transcriptional regulation of lipid transport. In a general population, we investigated whether genetic variation in these genes predicted coronary heart disease (CHD).

Results

In 2027 participants randomly recruited from a Flemish population (51.0 % women; mean age 43.6 years), we genotyped six SNPs in MEOX2 and four in TCF15. Over 15.2 years (median), CHD, myocardial infarction, coronary revascularisation and ischaemic cardiomyopathy occurred in 106, 53, 78 and 22 participants. For SNPs, we contrasted CHD risk in minor-allele heterozygotes and homozygotes (variant) vs. major-allele homozygotes (reference) and for haplotypes carriers (variant) vs. non-carriers. In multivariable-adjusted analyses with correction for multiple testing, CHD risk was associated with MEOX2 SNPs (P ≤ 0.049), but not with TCF15 SNPs (P ≥ 0.29). The MEOX2 GTCCGC haplotype (frequency 16.5 %) was associated with the sex- and age-standardised CHD incidence (5.26 vs. 3.03 events per 1000 person-years; P = 0.036); the multivariable-adjusted hazard ratio [HR] of CHD was 1.78 (95 % confidence interval, 1.25–2.56; P = 0.0054). For myocardial infarction, coronary revascularisation, and ischaemic cardiomyopathy, the corresponding HRs were 1.96 (1.16–3.31), 1.87 (1.20–2.91) and 3.16 (1.41–7.09), respectively. The MEOX2 GTCCGC haplotype significantly improved the prediction of CHD over and beyond traditional risk factors and was associated with similar population-attributable risk as smoking (18.7 % vs. 16.2 %).

Conclusions

Genetic variation in MEOX2, but not TCF15, is a strong predictor of CHD. Further experimental studies should elucidate the underlying molecular mechanisms.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0272-2) contains supplementary material, which is available to authorized users.

Imputation of KIR Types from SNP Variation Data

Large population studies of immune system genes are essential for characterizing their role in diseases, including autoimmune conditions. Of key interest are a group of genes encoding the killer cell immunoglobulin-like receptors (KIRs), which have known and hypothesized roles in autoimmune diseases, resistance to viruses, reproductive conditions, and cancer. These genes are highly polymorphic, which makes typing expensive and time consuming. Consequently, despite their importance, KIRs have been little studied in large cohorts. Statistical imputation methods developed for other complex loci (e.g., human leukocyte antigen [HLA]) on the basis of SNP data provide an inexpensive high-throughput alternative to direct laboratory typing of these loci and have enabled important findings and insights for many diseases. We present KIR^∗IMP, a method for imputation of KIR copy number. We show that KIR^∗IMP is highly accurate and thus allows the study of KIRs in large cohorts and enables detailed investigation of the role of KIRs in human disease.

Molgenis-impute: imputation pipeline in a box

Background

Genotype imputation is an important procedure in current genomic analysis such as genome-wide association studies, meta-analyses and fine mapping. Although high quality tools are available that perform the steps of this process, considerable effort and expertise is required to set up and run a best practice imputation pipeline, particularly for larger genotype datasets, where imputation has to scale out in parallel on computer clusters.

Results

Here we present MOLGENIS-impute, an ‘imputation in a box’ solution that seamlessly and transparently automates the set up and running of all the steps of the imputation process. These steps include genome build liftover (liftovering), genotype phasing with SHAPEIT2, quality control, sample and chromosomal chunking/merging, and imputation with IMPUTE2. MOLGENIS-impute builds on MOLGENIS-compute, a simple pipeline management platform for submission and monitoring of bioinformatics tasks in High Performance Computing (HPC) environments like local/cloud servers, clusters and grids. All the required tools, data and scripts are downloaded and installed in a single step. Researchers with diverse backgrounds and expertise have tested MOLGENIS-impute on different locations and imputed over 30,000 samples so far using the 1,000 Genomes Project and new Genome of the Netherlands data as the imputation reference. The tests have been performed on PBS/SGE clusters, cloud VMs and in a grid HPC environment.

Conclusions

MOLGENIS-impute gives priority to the ease of setting up, configuring and running an imputation. It has minimal dependencies and wraps the pipeline in a simple command line interface, without sacrificing flexibility to adapt or limiting the options of underlying imputation tools. It does not require knowledge of a workflow system or programming, and is targeted at researchers who just want to apply best practices in imputation via simple commands. It is built on the MOLGENIS compute workflow framework to enable customization with additional computational steps or it can be included in other bioinformatics pipelines. It is available as open source from: https://github.com/molgenis/molgenis-imputation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13104-015-1309-3) contains supplementary material, which is available to authorized users.