Mutations generate sequence diversity and provide a substrate for selection. The rate of de novo mutations is therefore of major importance to evolution. We conducted a study of genomewide mutation rate by sequencing the entire genomes of 78 Icelandic parent-offspring trios at high coverage. Here we show that in our samples, with an average father’s age of 29.7, the average de novo mutation rate is 1.20×10−8 per nucleotide per generation. Most strikingly, the diversity in mutation rate of single-nucleotide polymorphism (SNP) is dominated by the age of the father at conception of the child. The effect is an increase of about 2 mutations per year. After accounting for random Poisson variation, father’s age is estimated to explain nearly all of the remaining variation in the de novo mutation counts. These observations shed light on the importance of the father’s age on the risk of diseases such as schizophrenia and autism.

Objective

Rare copy number variants have been implicated in different neurodevelopmental disorders, with the same copy number variants often increasing risk of more than one of these phenotypes. In a discovery sample of 22 schizophrenia patients with an early onset of illness (10–15 years of age), the authors observed in one patient a maternally derived 15q11-q13 duplication overlapping the Prader-Willi/Angelman syndrome critical region. This prompted investigation of the role of 15q11-q13 duplications in psychotic illness.

Method

The authors scanned 7,582 patients with schizophrenia or schizoaffective disorder and 41,370 comparison subjects without known psychiatric illness for copy number variants at 15q11-q13 and determined the parental origin of duplications using methylation-sensitive Southern hybridization analysis.

Results

Duplications were found in four case patients and five comparison subjects. All four case patients had maternally derived duplications (0.05%), while only three of the five comparison duplications were maternally derived (0.007%), resulting in a significant excess of maternally derived duplications in case patients (odds ratio=7.3). This excess is compatible with earlier observations that risk for psychosis in people with Prader-Willi syndrome caused by maternal uniparental disomy is much higher than in those caused by deletion of the paternal chromosome.

Conclusions

These findings suggest that the presence of two maternal copies of a fragment of chromosome 15q11.2-q13.1 that overlaps with the Prader-Willi/Angelman syndrome critical region may be a rare risk factor for schizophrenia and other psychoses. Given that maternal duplications of this region are among the most consistent cytogenetic observations in autism, the findings provide further support for a shared genetic etiology between autism and psychosis.

IgA nephropathy (IgAN), major cause of kidney failure worldwide, is common in Asians, moderately prevalent in Europeans, and rare in Africans. It is not known if these differences represent variation in genes, environment, or ascertainment. In a recent GWAS, we localized five IgAN susceptibility loci on Chr.6p21 (HLA-DQB1/DRB1, PSMB9/TAP1, and DPA1/DPB2 loci), Chr.1q32 (CFHR3/R1 locus), and Chr.22q12 (HORMAD2 locus). These IgAN loci are associated with risk of other immune-mediated disorders such as type I diabetes, multiple sclerosis, or inflammatory bowel disease. We tested association of these loci in eight new independent cohorts of Asian, European, and African-American ancestry (N = 4,789), followed by meta-analysis with risk-score modeling in 12 cohorts (N = 10,755) and geospatial analysis in 85 world populations. Four susceptibility loci robustly replicated and all five loci were genome-wide significant in the combined cohort (P = 5×10−32–3×10−10), with heterogeneity detected only at the PSMB9/TAP1 locus (I2 = 0.60). Conditional analyses identified two new independent risk alleles within the HLA-DQB1/DRB1 locus, defining multiple risk and protective haplotypes within this interval. We also detected a significant genetic interaction, whereby the odds ratio for the HORMAD2 protective allele was reversed in homozygotes for a CFHR3/R1 deletion (P = 2.5×10−4). A seven–SNP genetic risk score, which explained 4.7% of overall IgAN risk, increased sharply with Eastward and Northward distance from Africa (r = 0.30, P = 3×10−128). This model paralleled the known East–West gradient in disease risk. Moreover, the prediction of a South–North axis was confirmed by registry data showing that the prevalence of IgAN–attributable kidney failure is increased in Northern Europe, similar to multiple sclerosis and type I diabetes. Variation at IgAN susceptibility loci correlates with differences in disease prevalence among world populations. These findings inform genetic, biological, and epidemiological investigations of IgAN and permit cross-comparison with other complex traits that share genetic risk loci and geographic patterns with IgAN.

Author Summary

IgA nephropathy (IgAN) is the most common cause of kidney failure in Asia, has lower prevalence in Europe, and is very infrequent among populations of African ancestry. A long-standing question in the field is whether these differences represent variation in genes, environment, or ascertainment. In a recent genome-wide association study of 5,966 individuals, we identified five susceptibility loci for this trait. In this paper, we study the largest IgAN case-control cohort reported to date, composed of 10,775 individuals of European, Asian, and African-American ancestry. We confirm that all five loci are significant contributors to disease risk across this multi-ethnic cohort. In addition, we identify two novel independent susceptibility alleles within the HLA-DQB1/DRB1 locus and a new genetic interaction between loci on Chr.1p36 and Chr.22q22. We develop a seven–SNP genetic risk score that explains nearly 5% of variation in disease risk. In geospatial analysis of 85 world populations, the genetic risk score closely parallels worldwide patterns of disease prevalence. The genetic risk score also predicts an unsuspected Northward risk gradient in Europe. This genetic prediction is verified by examination of registry data demonstrating, similarly to other immune-mediated diseases such as multiple sclerosis and type I diabetes, a previously unrecognized increase in IgAN–attributable kidney failure in Northern European countries.

Androgenetic alopecia (AGA) is a highly heritable condition and the most common form of hair loss in humans. Susceptibility loci have been described on the X chromosome and chromosome 20, but these loci explain a minority of its heritable variance. We conducted a large-scale meta-analysis of seven genome-wide association studies for early-onset AGA in 12,806 individuals of European ancestry. While replicating the two AGA loci on the X chromosome and chromosome 20, six novel susceptibility loci reached genome-wide significance (p = 2.62×10−9–1.01×10−12). Unexpectedly, we identified a risk allele at 17q21.31 that was recently associated with Parkinson's disease (PD) at a genome-wide significant level. We then tested the association between early-onset AGA and the risk of PD in a cross-sectional analysis of 568 PD cases and 7,664 controls. Early-onset AGA cases had significantly increased odds of subsequent PD (OR = 1.28, 95% confidence interval: 1.06–1.55, p = 8.9×10−3). Further, the AGA susceptibility alleles at the 17q21.31 locus are on the H1 haplotype, which is under negative selection in Europeans and has been linked to decreased fertility. Combining the risk alleles of six novel and two established susceptibility loci, we created a genotype risk score and tested its association with AGA in an additional sample. Individuals in the highest risk quartile of a genotype score had an approximately six-fold increased risk of early-onset AGA [odds ratio (OR) = 5.78, p = 1.4×10−88]. Our results highlight unexpected associations between early-onset AGA, Parkinson's disease, and decreased fertility, providing important insights into the pathophysiology of these conditions.

Author Summary

While most genome-wide association studies (GWAS) focus on the identification of susceptibility loci for a specific disease, this hypothesis-free approach also enables the identification of unexpected associations between different diseases by taking advantage of the previously published GWAS associations. Androgenetic Alopecia (AGA, also known as male pattern baldness) is the most common type of hair loss in humans. Parkinson's disease is reported to occur more commonly in men than in women; however, there are no studies investigating the link between AGA and Parkinson's disease. Here, we show that a specific genetic locus, chromosome 17q21.31, which is associated with Parkinson's disease, is also a susceptibility locus for early-onset AGA. We further investigate the association between early-onset AGA and Parkinson's disease, irrespective of genotype, directly in a large-scale web-based study. We find that men with early-onset AGA have 28% higher risk of developing Parkinson's disease. The early-onset AGA locus on chromosome 17q21.31 has also been linked to decreased fertility previously. Future studies of this locus may implicate novel biological pathways affecting these three conditions.

To identify genes involved in osteoarthritis (OA), the most prevalent form of joint disease, we performed a genome-wide association study (GWAS) in which we tested 500,510 Single Nucelotide Polymorphisms (SNPs) in 1341 OA cases and 3496 Dutch Caucasian controls. SNPs associated with at least two OA-phenotypes were analysed in 14,938 OA cases and approximately 39,000 controls. The C-allele of rs3815148 on chromosome 7q22 (MAF 23%, 172 kb upstream of the GPR22 gene) was consistently associated with a 1.14-fold increased risk (95%CI: 1.09–1.19) for knee- and/or hand-OA (p=8×10−8), and also with a 30% increased risk for knee-OA progression (95%CI: 1.03–1.64, p=0.03). This SNP is in almost complete linkage disequilibrium with rs3757713 (located 68 kb upstream of GPR22) which is associated with GPR22 expression levels in lymphoblast cell lines (p=4×10−12). GPR22 encodes an G-protein coupled receptor with unkown ligand (orphan receptor). Immunohistochemistry experiments showed absence of GPR22 in normal mouse articular cartilage or synovium. However, GPR22 positive chondrocytes were found in the upper layers of the articular cartilage of mouse knee joints that were challenged by in vivo papain treatment or in the presence of interleukin-1 driven inflammation. GRP22 positive chondrocyte-like cells were also found in osteophytes in instability-induced OA. In addition, GPR22 is also present in areas of the brain involved in locomotor function. Our findings reveal a novel common variant on chromosome 7q22 to influence susceptibility for prevalence and progression of OA.

Coffee is the most commonly used stimulant and caffeine is its main psychoactive ingredient. The heritability of coffee consumption has been estimated at around 50%. We performed a meta-analysis of four genome-wide association studies of coffee consumption among coffee drinkers from Iceland (n = 2680), the Netherlands (n = 2791), the Sorbs Slavonic population isolate in Germany (n = 771) and the USA (n = 369) using both directly genotyped and imputed single nucleotide polymorphisms (SNPs) (2.5 million SNPs). SNPs at the two most significant loci were also genotyped in a sample set from Iceland (n = 2430) and a Danish sample set consisting of pregnant women (n = 1620). Combining all data, two sequence variants significantly associated with increased coffee consumption: rs2472297-T located between CYP1A1 and CYP1A2 at 15q24 (P = 5.4 · 10−14) and rs6968865-T near aryl hydrocarbon receptor (AHR) at 7p21 (P = 2.3 · 10−11). An effect of ∼0.2 cups a day per allele was observed for both SNPs. CYP1A2 is the main caffeine metabolizing enzyme and is also involved in drug metabolism. AHR detects xenobiotics, such as polycyclic aryl hydrocarbons found in roasted coffee, and induces transcription of CYP1A1 and CYP1A2. The association of these SNPs with coffee consumption was present in both smokers and non-smokers.

DNA repair genes are important for maintaining genomic stability and limiting carcinogenesis. We analyzed all single nucleotide polymorphisms (SNPs) of 125 DNA repair genes covered by the Illumina HumanHap300 (v1.1) BeadChips in a previously conducted genome-wide association study (GWAS) of 1,154 lung cancer cases and 1,137 controls and replicated the top-hits of XRCC4 SNPs in an independent set of 597 cases and 611 controls in Texas populations. We found that six of 20 XRCC4 SNPs were associated with a decreased risk of lung cancer with a P value of 0.01 or lower in the discovery dataset, of which the most significant SNP was rs10040363 (P for allelic test = 4.89 ×10−4). Moreover, the data in this region allowed us to impute a potentially functional SNP rs2075685 (imputed P for allelic test = 1.3 ×10−3). A luciferase reporter assay demonstrated that the rs2075685G>T change in the XRCC4 promoter increased expression of the gene. In the replication study of rs10040363, rs1478486, rs9293329, and rs2075685, however, only rs10040363 achieved a borderline association with a decreased risk of lung cancer in a dominant model (adjusted OR = 0.80, 95% CI = 0.62–1.03, P = 0.079). In the final combined analysis of both the Texas GWAS discovery and replication datasets, the strength of the association was increased for rs10040363 (adjusted OR = 0.77, 95% CI = 0.66–0.89, Pdominant = 5×10−4 and P for trend = 5×10−4) and rs1478486 (adjusted OR = 0.82, 95% CI = 0.71 −0.94, Pdominant = 6×10−3 and P for trend = 3.5×10−3). Finally, we conducted a meta-analysis of these XRCC4 SNPs with available data from published GWA studies of lung cancer with a total of 12,312 cases and 47,921 controls, in which none of these XRCC4 SNPs was associated with lung cancer risk. It appeared that rs2075685, although associated with increased expression of a reporter gene and lung cancer risk in the Texas populations, did not have an effect on lung cancer risk in other populations. This study underscores the importance of replication using published data in larger populations.

Published genome-wide association studies (GWASs) have identified few variants in the known biological pathways involved in lung cancer etiology. To mine the possibly hidden causal single nucleotide polymorphisms (SNPs), we explored all SNPs in the extrinsic apoptosis pathway from our published GWAS dataset for 1154 lung cancer cases and 1137 cancer-free controls. In an initial association analysis of 611 tagSNPs in 41 apoptosis-related genes, we identified only 10 tagSNPs associated with lung cancer risk with a P value <10−2, including four tagSNPs in DAPK1 and three tagSNPs in TNFSF8. Unlike DAPK1 SNPs, TNFSF8 rs2181033 tagged other four predicted functional but untyped SNPs (rs776576, rs776577, rs31813148 and rs2075533) in the promoter region. Therefore, we further tested binding affinity of these four SNPs by performing the electrophoretic mobility shift assay. We found that only rs2075533T allele modified levels of nuclear proteins bound to DNA, leading to significantly decreased expression of luciferase reporter constructs by 5- to –10-fold in H1299, HeLa and HCT116 cell lines compared with the C allele. We also performed a replication study of the untyped rs2075533 in an independent Texas population but did not confirm the protective effect. We further performed a mini meta-analysis for SNPs of TNFSF8 obtained from other four published lung cancer GWASs with 12  214 cases and 47  721 controls, and we found that only rs3181366 (r2 = 0.69 with the untyped rs2075533) was associated to lung cancer risk (P = 0.008). Our findings suggest a possible role of novel TNFSF8 variants in susceptibility to lung cancer.

Genome-wide association studies (GWAS) have identified three genomic regions, at 15q24-25.1, 5p15.33 and 6p21.33, which associate with risk of lung cancer. Large meta-analyses of GWA data have failed to find additional associations of genome-wide significance. In this study, we sought to confirm 7 variants with suggestive association to lung cancer (P<10−5) in a recently published meta-analysis. In a GWA dataset of 1,447 lung cancer cases and 36,256 controls in Iceland, three correlated variants on 15q15.2 (rs504417, rs11853991 and rs748404) showed a significant association with lung cancer whereas rs4254535 on 2p14, rs1530057 on 3p24.1, rs6438347 on 3q13.31 and rs1926203 on 10q23.31 did not. The most significant variant, rs748404, was genotyped in additional 1,299 lung cancer cases and 4,102 controls from the Netherlands, Spain and the USA and the results combined with published GWAS data. In this analysis, the T allele of rs748404 reached genome-wide significance (OR=1.15, P=1.1×10−9). Another variant at the same locus, rs12050604, showed association with lung cancer (OR=1.09, 3.6×10−6) and remained significant after adjustment for rs748404 and vice versa. rs748404 is located 140 kb centromeric of the TP53BP1 gene that has been implicated in lung cancer risk. Two fully correlated, non-synonymous coding variants in TP53BP1, rs2602141 (Q1136K) and rs560191 (E353D), showed association with lung cancer in our sample set; however, this association did not remain significant after adjustment for rs748404. Our data show that one or more lung cancer risk variants of genome-wide significance and distinct from the coding variants in TP53BP1 are located at 15q15.2.

Background

Germline CDKN2A mutations have been observed in 20-40% of high-risk melanoma-prone families, however little is known about their prevalence in population-based series of melanoma cases and controls.

Methods

We resequenced the CDKN2A gene, including the p14ARF variant and promoter regions, in approximately 703 registry-ascertained melanoma cases and 691 population-based controls from Iceland, a country in which the incidence of melanoma has increased rapidly.

Results

We identified a novel germline variant, G89D that was strongly associated with increased melanoma risk and appeared to be an Icelandic founder mutation. The G89D variant was present in about 2% of Icelandic invasive cutaneous malignant melanoma cases. Relatives of affected G89D carriers were at significantly increased risk of melanoma, head & neck cancers, and pancreatic carcinoma compared to relatives of other melanoma patients. Nineteen other germline variants were identified, but none conferred an unequivocal risk of melanoma.

Conclusions

This population-based study of Icelandic melanoma cases and controls showed a frequency of disease-related CDKN2A mutant alleles ranging from 0.7% to 1.0%, thus expanding our knowledge about the frequency of CDKN2A mutations in different populations. In contrast to North America and Australia where a broad spectrum of mutations was observed at a similar frequency, in Iceland, functional CDKN2A mutations consists of only one or two different variants. Additional genetic and/or environmental factors are likely critical for explaining the high incidence rates for melanoma in Iceland. This study adds to the geographic regions for which population-based estimates of CDKN2A mutation frequencies are available.

Objective

To address the need for standardization of osteoarthritis (OA) phenotypes by examining the effect of heterogeneity among symptomatic (SOA) and radiographic osteoarthritis (ROA) phenotypes.

Methods

Descriptions of OA phenotypes of the 28 studies involved in the TREAT-OA consortium were collected. To investigate whether different OA definitions result in different association results, we created hip OA definitions used within the consortium in the Rotterdam Study-I and tested the association of hip OA with gender, age and BMI using one-way ANOVA. For radiographic OA, we standardized the hip, knee and hand ROA definitions and calculated prevalence's of ROA before and after standardization in 9 cohort studies. This procedure could only be performed in cohort studies and standardization of SOA definitions was not feasible at this moment.

Results

In this consortium, all studies with symptomatic OA phenotypes (knee, hip and hand) used a different definition and/or assessment of OA status. For knee, hip and hand radiographic OA 5, 4 and 7 different definitions were used, respectively. Different hip OA definitions do lead to different association results. For example, we showed in the Rotterdam Study-I that hip OA defined as “at least definite JSN and one definite osteophyte” was not associated with gender (p=0.22), but defined as “at least one definite osteophyte” was significantly associated with gender (p=3×10−9). Therefore, a standardization process was undertaken for radiographic OA definitions. Before standardization a wide range of ROA prevalence's was observed in the 9 cohorts studied. After standardization the range in prevalence of knee and hip ROA was small. Standardization of SOA phenotypes was not possible due to the case-control design of the studies.

Conclusion

Phenotype definitions influence the prevalence of OA and association with clinical variables. ROA phenotypes within the TREAT-OA consortium were standardized to reduce heterogeneity and improve power in future genetics studies.

Aims

We tested the hypothesis that the 9p21 risk locus promotes atherosclerosis by examining the association between rs10757278 and coronary artery disease (CAD) severity and progression determined by semi-quantitative angiographic scores.

Methods and results

The rs10757278 single nucleotide polymorphism (SNP) was genotyped as the marker for the 9p21 locus in 2334 Caucasian patients undergoing cardiac catheterization (mean age 63, male 67%). Angiographic CAD was assessed using two semi-quantitative scoring systems with one estimating severity (Gensini) and the other extent (Sullivan). A subset of 308 patients who underwent two or more coronary angiograms at least 6 months apart were examined for net change in Gensini and Sullivan scores over time to determine the rate of CAD progression by genotype and were further classified as ‘progressors’ or ‘non-progressors’ based on absolute change per year in angiographic severity score. We replicated the association between the rs10757278 SNP and myocardial infarction and binary (presence/absence) angiographic classifications of CAD. Furthermore, we observed a significant additive association with this SNP, and both severity and extent of CAD using angiographic scores, after adjustment for age, gender, body mass index, traditional cardiovascular risk factors, myocardial infarction, and statin use (Gensini P = 0.016, Sullivan P = 0.005). In addition, there was a significant linear association with CAD progression before and after adjustment for covariates (Gensini P = 0.023, Sullivan P = 0.003) with homozygotes for the risk variant having three-fold greater odds of CAD progression compared with the referent group.

Conclusion

The 9p21 risk locus is associated with angiographically defined severity, extent, and progression of CAD, suggesting a role for this locus in influencing atherosclerosis and its progression.

Genome-wide association studies have identified many genetic variants associated with complex traits. However, at only a minority of loci have the molecular mechanisms mediating these associations been characterized. In parallel, whilst cis-regulatory patterns of gene expression have been extensively explored, the identification of trans-regulatory effects in humans has attracted less attention. We demonstrate that the Type 2 diabetes and HDL-cholesterol associated cis-acting eQTL of the maternally-expressed transcription factor KLF14 acts as a master trans-regulator of adipose gene expression. Expression levels of genes regulated by this trans-eQTL are highly-correlated with concurrently-measured metabolic traits, and a subset of the trans-genes harbor variants directly-associated with metabolic phenotypes. This trans-eQTL network provides a mechanistic understanding of the effect of the KLF14 locus on metabolic disease risk, providing a potential model for other complex traits.

Osteoarthritis (OA) is the most common form of arthritis and a major cause of disability. This study evaluates the association in Caucasian populations of two single nucleotide polymorphisms (SNPs) mapping to the Human Leukocyte Antigen (HLA) region and deriving from a genome wide association scan (GWAS) of knee OA in Japanese populations. The frequencies for rs10947262 were compared in 36,408 controls and 5,749 knee OA cases from European-descent populations. rs7775228 was tested in 32,823 controls and 1,837 knee OA cases of European descent. The risk (major) allele at rs10947262 in Caucasian samples was not significantly associated with an odds ratio (OR)  = 1.07 (95%CI 0.94 -1.21; p = 0.28). For rs7775228 the meta-analysis resulted in OR = 0.94 (95%CI 0.81-1.09; p = 0.42) for the allele associated with risk in the Japanese GWAS. In Japanese individuals these two SNPs are in strong linkage disequilibrium (LD) (r2 = 0.86) with the HLA class II haplotype DRB1*1502 DQA1*0103 DQB1*0601 (frequency 8%). In Caucasian and Chinese samples, using imputed data, these SNPs appear not to be in LD with that haplotype (r2<0.07). The rs10947262 and rs7775228 variants are not associated with risk of knee OA in European descent populations and they do not appear tag the same HLA class II haplotype as they do in Japanese individuals.

A sequence variant (rs7216389-T) near the ORMDL3 gene on chromosome 17q21 was recently found to be associated with childhood asthma. We sought to evaluate the effect of rs7216389-T on asthma subphenotypes and its correlation with expression levels of neighboring genes. The association of rs7216389-T with asthma was replicated in six European and one Asian study cohort (N=4917 cases N=34 589 controls). In addition, we found that the association of rs7216389-T was confined to cases with early onset of asthma, particularly in early childhood (age: 0–5 years OR=1.51, P=6.89·10−9) and adolescence (age: 14–17 years OR=1.71, P=5.47·10−9). A weaker association was observed for onset between 6 and 13 years of age (OR=1.17, P=0.035), but none for adult-onset asthma (OR=1.07, P=0.12). Cases were further stratified by sex, asthma severity and atopy status. An association with greater asthma severity was observed among early-onset asthma cases (P=0.0012), but no association with sex or atopy status was observed among the asthma cases. An association between sequence variants and the expression of genes in the 17q21 region was assessed in white blood cell RNA samples collected from Icelandic individuals (n=743). rs7216389 associated with the expression of GSDMB and ORMDL3 genes. However, other sequence variants showing a weaker association with asthma compared with that of rs7216389 were more strongly associated with the expression of both genes. Thus, the contribution of rs7216389-T to the development of asthma is unlikely to operate only through an impact on the expression of ORMDL3 or GSDMB genes.

Background

Genome-wide association studies have identified three chromosomal regions at 15q25, 5p15, and 6p21 as being associated with the risk of lung cancer. To confirm these associations in independent studies and investigate heterogeneity of these associations within specific subgroups, we conducted a coordinated genotyping study within the International Lung Cancer Consortium based on independent studies that were not included in previous genome-wide association studies.

Methods

Genotype data for single-nucleotide polymorphisms at chromosomes 15q25 (rs16969968, rs8034191), 5p15 (rs2736100, rs402710), and 6p21 (rs2256543, rs4324798) from 21 case–control studies for 11 645 lung cancer case patients and 14 954 control subjects, of whom 85% were white and 15% were Asian, were pooled. Associations between the variants and the risk of lung cancer were estimated by logistic regression models. All statistical tests were two-sided.

Results

Associations between 15q25 and the risk of lung cancer were replicated in white ever-smokers (rs16969968: odds ratio [OR] = 1.26, 95% confidence interval [CI] = 1.21 to 1.32, Ptrend = 2 × 10−26), and this association was stronger for those diagnosed at younger ages. There was no association in never-smokers or in Asians between either of the 15q25 variants and the risk of lung cancer. For the chromosome 5p15 region, we confirmed statistically significant associations in whites for both rs2736100 (OR = 1.15, 95% CI = 1.10 to 1.20, Ptrend = 1 × 10−10) and rs402710 (OR = 1.14, 95% CI = 1.09 to 1.19, Ptrend = 5 × 10−8) and identified similar associations in Asians (rs2736100: OR = 1.23, 95% CI = 1.12 to 1.35, Ptrend = 2 × 10−5; rs402710: OR = 1.15, 95% CI = 1.04 to 1.27, Ptrend = .007). The associations between the 5p15 variants and lung cancer differed by histology; odds ratios for rs2736100 were highest in adenocarcinoma and for rs402710 were highest in adenocarcinoma and squamous cell carcinomas. This pattern was observed in both ethnic groups. Neither of the two variants on chromosome 6p21 was associated with the risk of lung cancer.

Conclusions

In this international genetic association study of lung cancer, previous associations found in white populations were replicated and new associations were identified in Asian populations. Future genetic studies of lung cancer should include detailed stratification by histology.

Despite significant progress in the identification of genetic loci for age-related macular degeneration (AMD), not all of the heritability has been explained. To identify variants which contribute to the remaining genetic susceptibility, we performed the largest meta-analysis of genome-wide association studies to date for advanced AMD. We imputed 6 036 699 single-nucleotide polymorphisms with the 1000 Genomes Project reference genotypes on 2594 cases and 4134 controls with follow-up replication of top signals in 5640 cases and 52 174 controls. We identified two new common susceptibility alleles, rs1999930 on 6q21-q22.3 near FRK/COL10A1 [odds ratio (OR) 0.87; P = 1.1 × 10−8] and rs4711751 on 6p12 near VEGFA (OR 1.15; P = 8.7 × 10−9). In addition to the two novel loci, 10 previously reported loci in ARMS2/HTRA1 (rs10490924), CFH (rs1061170, and rs1410996), CFB (rs641153), C3 (rs2230199), C2 (rs9332739), CFI (rs10033900), LIPC (rs10468017), TIMP3 (rs9621532) and CETP (rs3764261) were confirmed with genome-wide significant signals in this large study. Loci in the recently reported genes ABCA1 and COL8A1 were also detected with suggestive evidence of association with advanced AMD. The novel variants identified in this study suggest that angiogenesis (VEGFA) and extracellular collagen matrix (FRK/COL10A1) pathways contribute to the development of advanced AMD.

Background

Recent genome-wide association studies (GWAS) of myocardial infarction (MI) and other forms of coronary artery disease (CAD) have led to the discovery of at least 13 genetic loci. In addition to the effect size, power to detect associations is largely driven by sample size. Therefore, to maximize the chance of finding novel susceptibility loci for CAD and MI, the Coronary ARtery DIsease Genome-wide Replication And Meta-analysis (CARDIoGRAM) consortium was formed.

Methods and Results

CARDIoGRAM combines data from all published and several unpublished GWAS in individuals with European ancestry; includes >22 000 cases with CAD, MI, or both and >60 000 controls; and unifies samples from the Atherosclerotic Disease VAscular functioN and genetiC Epidemiology study, CADomics, Cohorts for Heart and Aging Research in Genomic Epidemiology, deCODE, the German Myocardial Infarction Family Studies I, II, and III, Ludwigshafen Risk and Cardiovascular Heath Study/AtheroRemo, MedStar, Myocardial Infarction Genetics Consortium, Ottawa Heart Genomics Study, PennCath, and the Wellcome Trust Case Control Consortium. Genotyping was carried out on Affymetrix or Illumina platforms followed by imputation of genotypes in most studies. On average, 2.2 million single nucleotide polymorphisms were generated per study. The results from each study are combined using meta-analysis. As proof of principle, we meta-analyzed risk variants at 9p21 and found that rs1333049 confers a 29% increase in risk for MI per copy (P=2×10−20).

Conclusion

CARDIoGRAM is poised to contribute to our understanding of the role of common genetic variation on risk for CAD and MI.

Through complementary application of SNP genotyping, whole-genome sequencing and imputation in 38,384 Icelanders, we have discovered a previously unidentified sick sinus syndrome susceptibility gene, MYH6, encoding the alpha heavy chain subunit of cardiac myosin. A missense variant in this gene, c.2161C>T, results in the conceptual amino acid substitution p.Arg721Trp, has an allelic frequency of 0.38% in Icelanders and associates with sick sinus syndrome with an odds ratio = 1 2.53 and P = 1.5 × 10−29. We show that the lifetime risk of being diagnosed with sick sinus syndrome is around 6% for non-carriers of c.2161C>T but is approximately 50% for carriers of the c.2161C>T variant.

Family studies of individual tissues have shown that gene expression traits are genetically heritable. Here, we investigate cis and trans components of heritability both within and across tissues by applying variance-components methods to 722 Icelanders from family cohorts, using identity-by-descent (IBD) estimates from long-range phased genome-wide SNP data and gene expression measurements for ∼19,000 genes in blood and adipose tissue. We estimate the proportion of gene expression heritability attributable to cis regulation as 37% in blood and 24% in adipose tissue. Our results indicate that the correlation in gene expression measurements across these tissues is primarily due to heritability at cis loci, whereas there is little sharing of trans regulation across tissues. One implication of this finding is that heritability in tissues composed of heterogeneous cell types is expected to be more dominated by cis regulation than in tissues composed of more homogeneous cell types, consistent with our blood versus adipose results as well as results of previous studies in lymphoblastoid cell lines. Finally, we obtained similar estimates of the cis components of heritability using IBD between unrelated individuals, indicating that transgenerational epigenetic inheritance does not contribute substantially to the “missing heritability” of gene expression in these tissue types.

Author Summary

An important goal in biology is to understand how genotype affects gene expression. Because gene expression varies across tissues, the relationship between genotype and gene expression may be tissue-specific. In this study, we used heritability approaches to study the regulation of gene expression in two tissue types, blood and adipose tissue, as well as the regulation of gene expression that is shared across these tissues. Heritability can be partitioned into cis and trans effects by assessing identity-by-descent (IBD) at the genomic location close to the expressed gene or genome-wide, respectively, and applying variance-components methods to partition the heritability of each gene. We estimated the proportion of gene expression heritability explained by cis regulation as 37% in blood and 24% in adipose tissue. Notably, the heritability shared across tissue types was primarily due to cis regulation. Thus, the relative contribution of cis versus trans regulation is expected to increase with the number of cell types present in the tissue being assayed, just as observed in our study and in a comparison to previous work on lymphoblastoid cell lines (LCL). We specifically ruled out a substantial contribution of transgenerational epigenetic inheritance to heritability of gene expression in these cohorts by repeating our heritability analyses using segments shared IBD in distantly related Icelanders.

Summary

Background

Large, rare chromosomal deletions and duplications known as copy number variants (CNVs) have been implicated in neurodevelopmental disorders similar to attention-deficit hyperactivity disorder (ADHD). We aimed to establish whether burden of CNVs was increased in ADHD, and to investigate whether identified CNVs were enriched for loci previously identified in autism and schizophrenia.

Methods

We undertook a genome-wide analysis of CNVs in 410 children with ADHD and 1156 unrelated ethnically matched controls from the 1958 British Birth Cohort. Children of white UK origin, aged 5–17 years, who met diagnostic criteria for ADHD or hyperkinetic disorder, but not schizophrenia and autism, were recruited from community child psychiatry and paediatric outpatient clinics. Single nucleotide polymorphisms (SNPs) were genotyped in the ADHD and control groups with two arrays; CNV analysis was limited to SNPs common to both arrays and included only samples with high-quality data. CNVs in the ADHD group were validated with comparative genomic hybridisation. We assessed the genome-wide burden of large (>500 kb), rare (<1% population frequency) CNVs according to the average number of CNVs per sample, with significance assessed via permutation. Locus-specific tests of association were undertaken for test regions defined for all identified CNVs and for 20 loci implicated in autism or schizophrenia. Findings were replicated in 825 Icelandic patients with ADHD and 35 243 Icelandic controls.

Findings

Data for full analyses were available for 366 children with ADHD and 1047 controls. 57 large, rare CNVs were identified in children with ADHD and 78 in controls, showing a significantly increased rate of CNVs in ADHD (0·156 vs 0·075; p=8·9×10−5). This increased rate of CNVs was particularly high in those with intellectual disability (0·424; p=2·0×10−6), although there was also a significant excess in cases with no such disability (0·125, p=0·0077). An excess of chromosome 16p13.11 duplications was noted in the ADHD group (p=0·0008 after correction for multiple testing), a finding that was replicated in the Icelandic sample (p=0·031). CNVs identified in our ADHD cohort were significantly enriched for loci previously reported in both autism (p=0·0095) and schizophrenia (p=0·010).

Interpretation

Our findings provide genetic evidence of an increased rate of large CNVs in individuals with ADHD and suggest that ADHD is not purely a social construct.

Funding

Action Research; Baily Thomas Charitable Trust; Wellcome Trust; UK Medical Research Council; European Union.

Most genome-wide association (GWA) studies have focused on populations of European ancestry with limited assessment of the influence of the sequence variants on populations of other ethnicities. To determine whether markers that we have recently shown to associate with Bone Mineral Density (BMD) in Europeans also associate with BMD in East-Asians we analysed 50 markers from 23 genomic loci in samples from Korea (n = 1,397) and two Chinese Hong Kong sample sets (n = 3,869 and n = 785). Through this effort we identified fourteen loci that associated with BMD in East-Asian samples using a false discovery rate (FDR) of 0.05; 1p36 (ZBTB40, P = 4.3×10−9), 1p31 (GPR177, P = 0.00012), 3p22 (CTNNB1, P = 0.00013), 4q22 (MEPE, P = 0.0026), 5q14 (MEF2C, P = 1.3×10−5), 6q25 (ESR1, P = 0.0011), 7p14 (STARD3NL, P = 0.00025), 7q21 (FLJ42280, P = 0.00017), 8q24 (TNFRSF11B, P = 3.4×10−5), 11p15 (SOX6, P = 0.00033), 11q13 (LRP5, P = 0.0033), 13q14 (TNFSF11, P = 7.5×10−5), 16q24 (FOXL1, P = 0.0010) and 17q21 (SOST, P = 0.015). Our study marks an early effort towards the challenge of cataloguing bone density variants shared by many ethnicities by testing BMD variants that have been established in Europeans, in East-Asians.

The increasing availability of personal genomic tests has led to discussions about the validity and utility of such tests and the balance of benefits and harms. A multidisciplinary workshop was convened by the National Institutes of Health and the Centers for Disease Control and Prevention to review the scientific foundation for using personal genomics in risk assessment and disease prevention and to develop recommendations for targeted research. The clinical validity and utility of personal genomics is a moving target with rapidly developing discoveries but little translation research to close the gap between discoveries and health impact. Workshop participants made recommendations in five domains: (1) developing and applying scientific standards for assessing personal genomic tests; (2) developing and applying a multidisciplinary research agenda, including observational studies and clinical trials to fill knowledge gaps in clinical validity and utility; (3) enhancing credible knowledge synthesis and information dissemination to clinicians and consumers; (4) linking scientific findings to evidence-based recommendations for use of personal genomics; and (5) assessing how the concept of personal utility can affect health benefits, costs, and risks by developing appropriate metrics for evaluation. To fulfill the promise of personal genomics, a rigorous multidisciplinary research agenda is needed.