Genome-wide association studies (GWAS) have now identified at least 2,000 common variants that appear associated with common diseases or related traits (http://www.genome.gov/gwastudies), hundreds of which have been convincingly replicated. It is generally thought that the associated markers reflect the effect of a nearby common (minor allele frequency >0.05) causal site, which is associated with the marker, leading to extensive resequencing efforts to find causal sites. We propose as an alternative explanation that variants much less common than the associated one may create "synthetic associations" by occurring, stochastically, more often in association with one of the alleles at the common site versus the other allele. Although synthetic associations are an obvious theoretical possibility, they have never been systematically explored as a possible explanation for GWAS findings. Here, we use simple computer simulations to show the conditions under which such synthetic associations will arise and how they may be recognized. We show that they are not only possible, but inevitable, and that under simple but reasonable genetic models, they are likely to account for or contribute to many of the recently identified signals reported in genome-wide association studies. We also illustrate the behavior of synthetic associations in real datasets by showing that rare causal mutations responsible for both hearing loss and sickle cell anemia create genome-wide significant synthetic associations, in the latter case extending over a 2.5-Mb interval encompassing scores of "blocks" of associated variants. In conclusion, uncommon or rare genetic variants can easily create synthetic associations that are credited to common variants, and this possibility requires careful consideration in the interpretation and follow up of GWAS signals.

A wealth of genetic associations for cardiovascular and metabolic phenotypes in humans has been accumulating over the last decade, in particular a large number of loci derived from recent genome wide association studies (GWAS). True complex disease-associated loci often exert modest effects, so their delineation currently requires integration of diverse phenotypic data from large studies to ensure robust meta-analyses. We have designed a gene-centric 50 K single nucleotide polymorphism (SNP) array to assess potentially relevant loci across a range of cardiovascular, metabolic and inflammatory syndromes. The array utilizes a "cosmopolitan" tagging approach to capture the genetic diversity across approximately 2,000 loci in populations represented in the HapMap and SeattleSNPs projects. The array content is informed by GWAS of vascular and inflammatory disease, expression quantitative trait loci implicated in atherosclerosis, pathway based approaches and comprehensive literature searching. The custom flexibility of the array platform facilitated interrogation of loci at differing stringencies, according to a gene prioritization strategy that allows saturation of high priority loci with a greater density of markers than the existing GWAS tools, particularly in African HapMap samples. We also demonstrate that the IBC array can be used to complement GWAS, increasing coverage in high priority CVD-related loci across all major HapMap populations. DNA from over 200,000 extensively phenotyped individuals will be genotyped with this array with a significant portion of the generated data being released into the academic domain facilitating in silico replication attempts, analyses of rare variants and cross-cohort meta-analyses in diverse populations. These datasets will also facilitate more robust secondary analyses, such as explorations with alternative genetic models, epistasis and gene-environment interactions.

High-throughput sequencing platforms are generating massive amounts of genetic variation data for diverse genomes, but it remains a challenge to pinpoint a small subset of functionally important variants. To fill these unmet needs, we developed the ANNOVAR tool to annotate single nucleotide variants (SNVs) and insertions/deletions, such as examining their functional consequence on genes, inferring cytogenetic bands, reporting functional importance scores, finding variants in conserved regions, or identifying variants reported in the 1000 Genomes Project and dbSNP. ANNOVAR can utilize annotation databases from the UCSC Genome Browser or any annotation data set conforming to Generic Feature Format version 3 (GFF3). We also illustrate a 'variants reduction' protocol on 4.7 million SNVs and indels from a human genome, including two causal mutations for Miller syndrome, a rare recessive disease. Through a stepwise procedure, we excluded variants that are unlikely to be causal, and identified 20 candidate genes including the causal gene. Using a desktop computer, ANNOVAR requires ∼4 min to perform gene-based annotation and ∼15 min to perform variants reduction on 4.7 million variants, making it practical to handle hundreds of human genomes in a day. ANNOVAR is freely available at http://www.openbioinformatics.org/annovar/.

This study aimed to identify the genetics underlying dominant forms of inherited retinal dystrophies using whole exome sequencing (WES) in six families extensively screened for known mutations or genes. Thirty-eight individuals were subjected to WES. Causative variants were searched among single nucleotide variants (SNVs) and insertion/deletion variants (indels) and whenever no potential candidate emerged, copy number variant (CNV) analysis was performed. Variants or regions harboring a candidate variant were prioritized and segregation of the variant with the disease was further assessed using Sanger sequencing in case of SNVs and indels, and quantitative PCR (qPCR) for CNVs. SNV and indel analysis led to the identification of a previously reported mutation in PRPH2. Two additional mutations linked to different forms of retinal dystrophies were identified in two families: a known frameshift deletion in RPGR, a gene responsible for X-linked retinitis pigmentosa and p.Ser163Arg in C1QTNF5 associated with Late-Onset Retinal Degeneration. A novel heterozygous deletion spanning the entire region of PRPF31 was also identified in the affected members of a fourth family, which was confirmed with qPCR. This study allowed the identification of the genetic cause of the retinal dystrophy and the establishment of a correct diagnosis in four families, including a large heterozygous deletion in PRPF31, typically considered one of the pitfalls of this method. Since all findings in this study are restricted to known genes, we propose that targeted sequencing using gene-panel is an optimal first approach for the genetic screening and that once known genetic causes are ruled out, WES might be used to uncover new genes involved in inherited retinal dystrophies.

Calcium, magnesium, potassium, sodium, chloride and phosphorus are the major dietary minerals involved in various biological functions and are commonly measured in the blood serum. Sufficient mineral intake is especially important for children due to their rapid growth. Currently, the genetic mechanisms influencing serum mineral levels are poorly understood, especially for children. We carried out a genome-wide association (GWA) study on 5,602 European-American children and 4,706 African-American children who had mineral measures available in their electronic medical records (EMR). While no locus met the criteria for genome-wide significant association, our results demonstrated a nominal association of total serum calcium levels with a missense variant in the calcium -sensing receptor (CASR) gene on 3q13 (rs1801725, P = 1.96 × 10(-3)) in the African-American pediatric cohort, a locus previously reported in Caucasians. We also confirmed the association result in our pediatric European-American cohort (P = 1.38 × 10(-4)). We further replicated two other loci associated with serum calcium levels in the European-American cohort (rs780094, GCKR, P = 4.26 × 10(-3); rs10491003, GATA3, P = 0.02). In addition, we replicated a previously reported locus on 1q21, demonstrating association of serum magnesium levels with MUC1 (rs4072037, P = 2.04 × 10(-6)). Moreover, in an extended gene-based association analysis we uncovered evidence for association of calcium levels with the previously reported gene locus DGKD in both European-American children and African-American children. Taken together, our results support a role for CASR and DGKD mediated calcium regulation in both African-American and European-American children, and corroborate the association of calcium levels with GCKR and GATA3, and the association of magnesium levels with MUC1 in the European-American children.

Eosinophilic esophagitis (EoE) is an allergic disorder characterized by infiltration of the oesophagus with eosinophils. We had previously reported association of the TSLP/WDR36 locus with EoE. Here we report genome-wide significant associations at four additional loci; c11orf30 and STAT6, which have been previously associated with both atopic and autoimmune diseases, and two EoE-specific loci, ANKRD27 that regulates the trafficking of melanogenic enzymes to epidermal melanocytes and CAPN14, that encodes a calpain whose expression is highly enriched in the oesophagus. The identification of five EoE loci, not only expands our aetiological understanding of the disease but may also represent new therapeutic targets to treat the most debilitating aspect of EoE, oesophageal inflammation and remodelling.

Candidate gene studies have revealed limited genetic bases for opioid analgesic response variability. Genome-wide association studies facilitate impartial queries of common genetic variants, allowing identification of novel genetic contributions to drug effect. Illumina (Illumina Inc, San Diego, CA, USA) single nucleotide polymorphism (SNP) arrays were used to investigate SNP associations with total morphine requirement as a quantitative trait locus and with postoperative pain in a retrospective population of opioid-naïve children ages 4-18years who had undergone day surgery tonsillectomy and adenoidectomy. In an independent replication cohort, significant genome-wide association studies-identified SNPs were assayed using TaqMan probes. Among 617 comprehensively phenotyped children, the 277 subjects of European Caucasian (EC) ancestry demonstrated nominal association between morphine dose and a series of novel SNPs (top rs795484, P=1.01 × 10(-6) and rs1277441, P=2.77 × 10(-6)) at the TAOK3 locus. Age, body mass index, and physical status were included covariates. Morphine requirement averaged 132.4 μg/kg (SD 40.9). Each minor allele at rs795484 (guanine [G]>adenine [A]) contributed +17.6 μg/kg (95% confidence interval [CI] 10.7-24.4) to dose. Effect direction and magnitude were replicated in an independent cohort of 75 EC children (P<0.05). No association with morphine dose was detected in African Americans (AA) (n=241). Postoperative pain scores ≥ 7/10 were associated with rs795484 (G>A) in the EC cohort (odds ratio 2.35, 95% CI 1.56-3.52, P<0.00005) and this association replicated in AA children (odds ratio 1.76, 95% CI 1.14-2.71, P<0.01). Variants in TAOK3 encoding the serine/threonine-protein kinase, TAO3, are associated with increased morphine requirement in children of EC ancestry and with increased acute postoperative pain in both EC and AA subjects.

Rare copy-number variation (CNV) is an important source of risk for autism spectrum disorders (ASDs). We analyzed 2,446 ASD-affected families and confirmed an excess of genic deletions and duplications in affected versus control groups (1.41-fold, p = 1.0 × 10(-5)) and an increase in affected subjects carrying exonic pathogenic CNVs overlapping known loci associated with dominant or X-linked ASD and intellectual disability (odds ratio = 12.62, p = 2.7 × 10(-15), ∼3% of ASD subjects). Pathogenic CNVs, often showing variable expressivity, included rare de novo and inherited events at 36 loci, implicating ASD-associated genes (CHD2, HDAC4, and GDI1) previously linked to other neurodevelopmental disorders, as well as other genes such as SETD5, MIR137, and HDAC9. Consistent with hypothesized gender-specific modulators, females with ASD were more likely to have highly penetrant CNVs (p = 0.017) and were also overrepresented among subjects with fragile X syndrome protein targets (p = 0.02). Genes affected by de novo CNVs and/or loss-of-function single-nucleotide variants converged on networks related to neuronal signaling and development, synapse function, and chromatin regulation.

The goal of this paper is to review recent research on copy number variations (CNVs) and their association with complex and rare diseases. In the latter part of this paper, we focus on how large biorepositories such as the electronic medical record and genomics (eMERGE) consortium may be best leveraged to systematically mine for potentially pathogenic CNVs, and we end with a discussion of how such variants might be reported back for inclusion in electronic medical records as part of medical history.

Estimates from large scale genome sequencing studies indicate that each human carries up to 20 genetic variants that are predicted to results in loss of function (LOF) of protein-coding genes. While some are known disease-causing variants or common, tolerated, LOFs in non-essential genes, the majority remain of unknown consequence. We explore the possibility of using imputed GWAS data from large biorepositories such as the electronic medical record and genomics (eMERGE) consortium to determine the effects of rare LOFs. Here, we show that two hypocholesterolemia-associated LOF mutations in the PCSK9 gene can be accurately imputed into large-scale GWAS datasets which raises the possibility of assessing LOFs through genomics-linked medical records.

Acute otitis media (AOM) is among the most common pediatric diseases, and the most frequent reason for antibiotic treatment in children. Risk of AOM is dependent on environmental and host factors, as well as a significant genetic component. We identify genome-wide significance at a locus on 6q25.3 (rs2932989, Pmeta=2.15 × 10-09), and show that the associated variants are correlated with the methylation status of the FNDC1 gene (cg05678571, P=1.43 × 10-06), and further show it is an eQTL for FNDC1 (P=9.3 × 10-05). The mouse homologue, Fndc1, is expressed in middle ear tissue and its expression is upregulated upon lipopolysaccharide treatment. In this first GWAS of AOM and the largest OM genetic study to date, we identify the first genome-wide significant locus associated with AOM.

Statin treatment and variants in the gene encoding HMG-CoA reductase are associated with reductions in both the concentration of LDL cholesterol and the risk of coronary heart disease, but also with modest hyperglycaemia, increased bodyweight, and modestly increased risk of type 2 diabetes, which in no way offsets their substantial benefits. We sought to investigate the associations of LDL cholesterol-lowering PCSK9 variants with type 2 diabetes and related biomarkers to gauge the likely effects of PCSK9 inhibitors on diabetes risk.

Hepatoblastoma (HBL), the most common childhood liver cancer is cured with surgical resection after chemotherapy or with liver transplantation if local invasion and multifocality preclude resection. However, variable survival rates of 60-80% and debilitating chemotherapy sequelae argue for more informed treatment selection, which is not possible by grading the Wnt-β-catenin over activity present in most HBL tumors. A hypothesis-generating whole transcriptome analysis shows that HBL tumors removed at transplantation are enriched most for cancer signaling pathways which depend predominantly on epidermal growth factor (EGF) signaling, and to a lesser extent, on aberrant Wnt-β-catenin signaling. We therefore evaluated whether EGFR, ASAP1, ERBB2 and ERBB4, which signal downstream after ligation of EGF, and which show aberrant expression in several other invasive cancers, would also predict HBL tumor invasiveness. Immunohistochemistry of HBL tumors (n = 60), which are histologically heterogeneous, shows that compared with well-differentiated fetal cells, less differentiated embryonal and undifferentiated small cells (SCU) progressively lose EGFR and ASAP1 expression. This trend is exaggerated in unresectable, locally invasive or metastatic tumors, in which embryonal tumor cells are EGFR-negative, while SCU cells are EGFR-negative and ASAP1-negative. Loss of EGFR-ASAP1 signaling characterizes undifferentiated and invasive HBL. EGFR-expressing HBL tumors present novel therapeutic targeting opportunities.

To facilitate the clinical implementation of genomic medicine by next-generation sequencing, it will be critically important to obtain accurate and consistent variant calls on personal genomes. Multiple software tools for variant calling are available, but it is unclear how comparable these tools are or what their relative merits in real-world scenarios might be.

Social communication difficulties represent an autistic trait that is highly heritable and persistent during the course of development. However, little is known about the underlying genetic architecture of this phenotype.

We performed risk assessment for Crohn's disease (CD) and ulcerative colitis (UC), the two common forms of inflammatory bowel disease (IBD), by using data from the International IBD Genetics Consortium's Immunochip project. This data set contains ~17,000 CD cases, ~13,000 UC cases, and ~22,000 controls from 15 European countries typed on the Immunochip. This custom chip provides a more comprehensive catalog of the most promising candidate variants by picking up the remaining common variants and certain rare variants that were missed in the first generation of GWAS. Given this unprecedented large sample size and wide variant spectrum, we employed the most recent machine-learning techniques to build optimal predictive models. Our final predictive models achieved areas under the curve (AUCs) of 0.86 and 0.83 for CD and UC, respectively, in an independent evaluation. To our knowledge, this is the best prediction performance ever reported for CD and UC to date.

In addition to HLA genetic incompatibility, non-HLA difference between donor and recipients of transplantation leading to allograft rejection are now becoming evident. We aimed to create a unique genome-wide platform to facilitate genomic research studies in transplant-related studies. We designed a genome-wide genotyping tool based on the most recent human genomic reference datasets, and included customization for known and potentially relevant metabolic and pharmacological loci relevant to transplantation.

The standard expression quantitative trait loci (eQTL) detects polymorphisms associated with gene expression without revealing causality. We introduce a coupled Bayesian regression approach--eQTeL, which leverages epigenetic data to estimate regulatory and gene interaction potential, and identifies combination of regulatory single-nucleotide polymorphisms (SNPs) that explain the gene expression variance. On human heart data, eQTeL not only explains a significantly greater proportion of expression variance but also predicts gene expression more accurately than other methods. Based on realistic simulated data, we demonstrate that eQTeL accurately detects causal regulatory SNPs, including those with small effect sizes. Using various functional data, we show that SNPs detected by eQTeL are enriched for allele-specific protein binding and histone modifications, which potentially disrupt binding of core cardiac transcription factors and are spatially proximal to their target. eQTeL SNPs capture a substantial proportion of genetic determinants of expression variance and we estimate that 58% of these SNPs are putatively causal.

Absence of the anterior (ACL) or posterior cruciate ligament (PCL) are rare congenital malformations that result in knee joint instability, with a prevalence of 1.7 per 100,000 live births and can be associated with other lower-limb abnormalities such as ACL agnesia and absence of the menisci of the knee. While a few cases of absence of ACL/PCL are reported in the literature, a number of large familial case series of related conditions such as ACL agnesia suggest a potential underlying monogenic etiology. We performed whole exome sequencing of a family with two individuals affected by ACL/PCL.

Genetic association studies have identified 21 loci associated with atopic dermatitis risk predominantly in populations of European ancestry. To identify further susceptibility loci for this common, complex skin disease, we performed a meta-analysis of >15 million genetic variants in 21,399 cases and 95,464 controls from populations of European, African, Japanese and Latino ancestry, followed by replication in 32,059 cases and 228,628 controls from 18 studies. We identified ten new risk loci, bringing the total number of known atopic dermatitis risk loci to 31 (with new secondary signals at four of these loci). Notably, the new loci include candidate genes with roles in the regulation of innate host defenses and T cell function, underscoring the important contribution of (auto)immune mechanisms to atopic dermatitis pathogenesis.