Renal cell carcinoma (RCC) is the most lethal urologic cancer. Only two common susceptibility loci for RCC have been confirmed to date. To identify additional RCC common susceptibility loci, we conducted an independent genome-wide association study (GWAS). We analyzed 533 191 single nucleotide polymorphisms (SNPs) for association with RCC in 894 cases and 1516 controls of European descent recruited from MD Anderson Cancer Center in the primary scan, and validated the top 500 SNPs in silico in 3772 cases and 8505 controls of European descent involved in the only published GWAS of RCC. We identified two common variants in linkage disequilibrium, rs718314 and rs1049380 (r2 = 0.64, D ′ = 0.84), in the inositol 1,4,5-triphosphate receptor, type 2 (ITPR2) gene on 12p11.23 as novel susceptibility loci for RCC (P = 8.89 × 10−10 and P = 6.07 × 10−9, respectively, in meta-analysis) with an allelic odds ratio of 1.19 [95% confidence interval (CI): 1.13–1.26] for rs718314 and 1.18 (95% CI: 1.12–1.25) for rs1049380. It has been recently identified that rs718314 in ITPR2 is associated with waist–hip ratio (WHR) phenotype. To our knowledge, this is the first genetic locus associated with both cancer risk and WHR.

We conducted a genome-wide association study to search for risk alleles associated with Tetralogy of Fallot (TOF), using a northern European discovery set of 835 cases and 5159 controls. A region on chromosome 12q24 was associated (P = 1.4 × 10−7) and replicated convincingly (P = 3.9 × 10−5) in 798 cases and 2931 controls [per allele odds ratio (OR) = 1.27 in replication cohort, P = 7.7 × 10−11 in combined populations]. Single nucleotide polymorphisms in the glypican 5 gene on chromosome 13q32 were also associated (P = 1.7 × 10−7) and replicated convincingly (P = 1.2 × 10−5) in 789 cases and 2927 controls (per allele OR = 1.31 in replication cohort, P = 3.03 × 10−11 in combined populations). Four additional regions on chromosomes 10, 15 and 16 showed suggestive association accompanied by nominal replication. This study, the first genome-wide association study of a congenital heart malformation phenotype, provides evidence that common genetic variation influences the risk of TOF.

Among inherited retinal dystrophies, autosomal recessive retinitis pigmentosa (arRP) is the most genetically heterogenous condition with 32 genes currently known that account for ∼60 % of patients. Molecular diagnosis thus requires the tedious systematic sequencing of 506 exons. To rapidly identify the causative mutations, we devised a strategy that combines gene mapping and phenotype assessment in small non-consanguineous families. Two unrelated sibships with arRP had whole-genome scan using SNP microchips. Chromosomal regions were selected by calculating a score based on SNP coverage and genotype identity of affected patients. Candidate genes from the regions with the highest scores were then selected based on phenotype concordance of affected patients with previously described phenotype for each candidate gene. For families RP127 and RP1459, 33 and 40 chromosomal regions showed possible linkage, respectively. By comparing the scores with the phenotypes, we ended with one best candidate gene for each family, namely tubby-like protein 1 (TULP1) and C2ORF71 for RP127 and RP1459, respectively. We found that RP127 patients were compound heterozygous for two novel TULP1 mutations, p.Arg311Gln and p.Arg342Gln, and that RP1459 patients were compound heterozygous for two novel C2ORF71 mutations, p.Leu777PhefsX34 and p.Leu777AsnfsX28. Phenotype assessment showed that TULP1 patients had severe early onset arRP and that C2ORF71 patients had a cone rod dystrophy type of arRP. Only two affected individuals in each sibship were sufficient to lead to mutation identification by screening the best candidate gene selected by a combination of gene mapping and phenotype characterization.

Primary ciliary dyskinesia (PCD) is an inherited disorder characterized by recurrent infections of the upper and lower respiratory tract, reduced fertility in males and situs inversus in about 50% of affected individuals (Kartagener syndrome). It is caused by motility defects in the respiratory cilia that are responsible for airway clearance, the flagella that propel sperm cells and the nodal monocilia that determine left-right asymmetry1. Recessive mutations that cause PCD have been identified in genes encoding components of the outer dynein arms, radial spokes and cytoplasmic pre-assembly factors of axonemal dyneins, but these mutations account for only about 50% of cases of PCD. We exploited the unique properties of dog populations to positionally clone a new PCD gene, CCDC39. We found that loss-of-function mutations in the human ortholog underlie a substantial fraction of PCD cases with axonemal disorganization and abnormal ciliary beating. Functional analyses indicated that CCDC39 localizes to ciliary axonemes and is essential for assembly of inner dynein arms and the dynein regulatory complex.

Background

The release of the porcine genome sequence offers great perspectives for Pig genetics and genomics, and more generally will contribute to the understanding of mammalian genome biology and evolution. The process of producing a complete genome sequence of high quality, while facilitated by high-throughput sequencing technologies, remains a difficult task. The porcine genome was sequenced using a combination of a hierarchical shotgun strategy and data generated with whole genome shotgun. In addition to the BAC contig map used for the clone-by-clone approach, genomic mapping resources for the pig include two radiation hybrid (RH) panels at two different resolutions. These two panels have been used extensively for the physical mapping of pig genes and markers prior to the availability of the pig genome sequence.

Results

In order to contribute to the assembly of the pig genome, we genotyped the two radiation hybrid (RH) panels with a SNP array (the Illumina porcineSNP60 array) and produced high density physical RH maps for each pig autosome. We first present the methods developed to obtain high density RH maps with 38,379 SNPs from the SNP array genotyping. We then show how they were useful to identify problems in a draft of the pig genome assembly, and how the RH maps enabled the problems to be corrected in the porcine genome sequence. Finally, we used the RH maps to predict the position of 2,703 SNPs and 1,328 scaffolds currently unplaced on the porcine genome assembly.

Conclusions

A complete process, from genotyping of a high density SNP array on RH panels, to the construction of genome-wide high density RH maps, and finally their exploitation for validating and improving a genome assembly is presented here. The study includes the cross-validation of RH based findings with independent information from genetic data and comparative mapping with the Human genome. Several additional resources are also provided, in particular the predicted genomic location of currently unplaced SNPs and associated scaffolds summing up to a total of 72 megabases, that can be useful for the exploitation of the pig genome assembly.

While gliomas are the most common primary brain tumors, their etiology is largely unknown. To identify novel risk loci for glioma, we conducted genome-wide association (GWA) analysis of two case–control series from France and Germany (2269 cases and 2500 controls). Pooling these data with previously reported UK and US GWA studies provided data on 4147 glioma cases and 7435 controls genotyped for 424 460 common tagging single-nucleotide polymorphisms. Using these data, we demonstrate two statistically independent associations between glioma and rs11979158 and rs2252586, at 7p11.2 which encompasses the EGFR gene (population-corrected statistics, Pc = 7.72 × 10−8 and 2.09 × 10−8, respectively). Both associations were independent of tumor subtype, and were independent of EGFR amplification, p16INK4a deletion and IDH1 mutation status in tumors; compatible with driver effects of the variants on glioma development. These findings show that variation in 7p11.2 is a determinant of inherited glioma risk.

Variants associated with meconium ileus in cystic fibrosis (CF) were identified in 3,763 patients by GWAS. Five SNPs at two loci near SLC6A14 (min P=1.28×10−12 at rs3788766), chr Xq23-24 and SLC26A9 (min P=9.88×10−9 at rs4077468), chr 1q32.1 accounted for ~5% of the phenotypic variability, and were replicated in an independent patient collection (n=2,372; P=0.001 and 0.0001 respectively). By incorporating that disease-causing mutations in CFTR alter electrolyte and fluid flux across epithelia into an hypothesis-driven genome-wide analysis (GWAS-HD), we identified the same SLC6A14 and SLC26A9 associated SNPs, while establishing evidence for the involvement of SNPs in a third solute carrier gene, SLC9A3. In addition, GWAS-HD provided evidence of association between meconium ileus and multiple constituents of the apical plasma membrane where CFTR resides (P=0.0002, testing 155 apical genes jointly and replicated, P=0.022). These findings suggest that modulating activities of apical membrane constituents could complement current therapeutic paradigms for cystic fibrosis.

The manifestation of coronary artery disease (CAD) follows a well-choreographed series of events that includes damage of arterial endothelial cells and deposition of lipids in the sub-endothelial layers. Genome-wide association studies (GWAS) of multiple populations with distinctive genetic and lifestyle backgrounds are a crucial step in understanding global CAD pathophysiology. In this study, we report a GWAS on the genetic basis of arterial stenosis as measured by cardiac catheterization in a Lebanese population. The locus of the phosphatase and actin regulator 1 gene (PHACTR1) showed association with coronary stenosis in a discovery experiment with genome wide data in 1,949 individuals (rs9349379, OR = 1.37, p = 1.57×10−5). The association was replicated in an additional 2,547 individuals (OR = 1.31, p = 8.85×10−6), leading to genome-wide significant association in a combined analysis (OR = 1.34, p = 8.02×10−10). Results from this GWAS support a central role of PHACTR1 in CAD susceptibility irrespective of lifestyle and ethnic divergences. This association provides a plausible component for understanding molecular mechanisms involved in the formation of stenosis in cardiac vessels and a potential drug target against CAD.

By applying an imputation strategy based on the 1000 Genomes project to two genome-wide association studies (GWAS), we detected a susceptibility locus for venous thrombosis on chromosome 11p11.2 that was missed by previous GWAS analyses that had been conducted on the same datasets. A comprehensive linkage disequilibrium and haplotype analysis of the whole locus where twelve SNPs exhibited association p-values lower than 2.23 10−11 and the use of independent case-control samples demonstrated that the culprit variant was a rare variant located ∼1 Mb away from the original hits, not tagged by current genome-wide genotyping arrays and even not well imputed in the original GWAS samples. This variant was in fact the rs1799963, also known as the FII G20210A prothrombin mutation. This work may be of major interest not only for its scientific impact but also for its methodological findings.

Our specific aims were to evaluate the power of bivariate analysis and to compare its performance with traditional univariate analysis in samples of unrelated subjects under varying sampling selection designs. Bivariate association analysis was based on the seemingly unrelated regression (SUR) model that allows different genetic models for different traits. We conducted extensive simulations for the case of two correlated quantitative phenotypes, with the quantitative trait locus making equal or unequal contributions to each phenotype. Our simulation results confirmed that the power of bivariate analysis is affected by the size, direction and source of the phenotypic correlations between traits. They also showed that the optimal sampling scheme depends on the size and direction of the induced genetic correlation. In addition, we demonstrated the efficacy of SUR-based bivariate test by applying it to a real Genome-Wide Association Study (GWAS) of Bone Mineral Density (BMD) values measured at the lumbar spine (LS) and at the femoral neck (FN) in a sample of unrelated males with low BMD (LS Z-scores ≤−2) and with high BMD (LS and FN Z-scores >0.5). A substantial amount of top hits in bivariate analysis did not reach nominal significance in any of the two single-trait analyses. Altogether, our studies suggest that bivariate analysis is of practical significance for GWAS of correlated phenotypes.

Background

Genome-wide association studies (GWAS) require large sample sizes to obtain adequate statistical power, but it may be possible to increase the power by incorporating complementary data. In this study we investigated the feasibility of automatically retrieving information from the medical literature and leveraging this information in GWAS.

Methods

We developed a method that searches through PubMed abstracts for pre-assigned keywords and key concepts, and uses this information to assign prior probabilities of association for each single nucleotide polymorphism (SNP) with the phenotype of interest - the Adjusting Association Priors with Text (AdAPT) method. Association results from a GWAS can subsequently be ranked in the context of these priors using the Bayes False Discovery Probability (BFDP) framework. We initially tested AdAPT by comparing rankings of known susceptibility alleles in a previous lung cancer GWAS, and subsequently applied it in a two-phase GWAS of oral cancer.

Results

Known lung cancer susceptibility SNPs were consistently ranked higher by AdAPT BFDPs than by p-values. In the oral cancer GWAS, we sought to replicate the top five SNPs as ranked by AdAPT BFDPs, of which rs991316, located in the ADH gene region of 4q23, displayed a statistically significant association with oral cancer risk in the replication phase (per-rare-allele log additive p-value [ptrend] = 2.5×10−3). The combined OR for having one additional rare allele was 0.83 (95% CI: 0.76–0.90), and this association was independent of previously identified susceptibility SNPs that are associated with overall UADT cancer in this gene region. We also investigated if rs991316 was associated with other cancers of the upper aerodigestive tract (UADT), but no additional association signal was found.

Conclusion

This study highlights the potential utility of systematically incorporating prior knowledge from the medical literature in genome-wide analyses using the AdAPT methodology. AdAPT is available online (url: http://services.gate.ac.uk/lld/gwas/service/config).

The insulin (INS) region is the second most important locus associated with Type 1 Diabetes (T1D). The study of the DNA methylation pattern of the 7 CpGs proximal to the TSS in the INS gene promoter revealed that T1D patients have a lower level of methylation of CpG -19, -135 and -234 (p = 2.10−16) and a higher methylation of CpG -180 than controls, while methylation was comparable for CpG -69, -102, -206. The magnitude of the hypomethylation relative to a control population was 8–15% of the corresponding levels in controls and was correlated in CpGs -19 and -135 (r = 0.77) and CpG -135 and -234 (r = 0.65). 70/485 (14%) of T1D patients had a simultaneous decrease in methylation of CpG -19, -135, -234 versus none in 317 controls. CpG methylation did not correlate with glycated hemoglobin or with T1D duration. The methylation of CpG -69, -102, -180, -206, but not CpG -19, -135, -234 was strongly influenced by the cis-genotype at rs689, a SNP known to show a strong association with T1D. We hypothesize that part of this genetic association could in fact be mediated at the statistical and functional level by the underlying changes in neighboring CpG methylation. Our observation of a CpG-specific, locus-specific methylation pattern, although it can provide an epigenetic biomarker of a multifactorial disease, does not indicate whether the reported epigenetic pattern preexists or follows the establishment of T1D. To explore the effect of chronic hyperglycemia on CpG methylation, we studied non obese patients with type 2 diabetes (T2D) who were found to have decreased CpG-19 methylation versus age-matched controls, similar to T1D (p = 2.10−6) but increased CpG-234 methylation (p = 5.10−8), the opposite of T1D. The causality and natural history of the different epigenetic changes associated with T1D or T2D remain to be determined.

We conducted a two-stage genome-wide association study of renal cell carcinoma (RCC) in 3,772 cases and 8,505 controls of European background from 11 studies, and followed up 6 SNPs in three replication studies of 2,198 cases and 4,918 controls. Two loci on the regions of 2p21 and 11q13.3 were associated with RCC susceptibility below genome-wide significance. Two correlated variants (r2 = 0.99 in controls), rs11894252 (P = 1.8×10−8) and rs7579899 (P = 2.3×10−9), map to EPAS1 on 2p21, which encodes hypoxia-inducible- factor-2 alpha, a transcription factor previously implicated in RCC. The second locus, rs7105934, at 11q13, contains no characterized genes (P = 7.8×10−14). In addition, we observed a promising association on 12q24.31 for rs4765623 which maps to the scavenger receptor class B, member 1 (SCARB1) gene (P = 2.6×10−8). Our study reports novel genomic regions associated with RCC risk that may lead to new etiological insights.

Background

Venous Thrombosis (VT) is a common multifactorial disease associated with a major public health burden. Genetics factors are known to contribute to the susceptibility of the disease but how many genes are involved and their contribution to VT risk still remain obscure. We aimed to identify genetic variants associated with VT risk.

Methodology/Principal Findings

We conducted a genome-wide association study (GWAS) based on 551,141 SNPs genotyped in 1,542 cases and 1,110 controls. Twelve SNPs reached the genome-wide significance level of 2.0×10−8 and encompassed four known VT-associated loci, ABO, F5, F11 and FGG. By means of haplotype analyses, we also provided novel arguments in favor of a role of HIVEP1, PROCR and STAB2, three loci recently hypothesized to participate in the susceptibility to VT. However, no novel VT-associated loci came out of our GWAS. Using a recently proposed statistical methodology, we also showed that common variants could explain about 35% of the genetic variance underlying VT susceptibility among which 3% could be attributable to the main identified VT loci. This analysis additionally suggested that the common variants left to be identified are not uniformly distributed across the genome and that chromosome 20, itself, could contribute to ∼7% of the total genetic variance.

Conclusions/Significance

This study might also provide a valuable source of information to expand our understanding of biological mechanisms regulating quantitative biomarkers for VT.

Background

Thirty-two common variants associated with body mass index (BMI) have been identified in genome-wide association studies, explaining ∼1.45% of BMI variation in general population cohorts. We performed a genome-wide association study in a sample of young adults enriched for extremely overweight individuals. We aimed to identify new loci associated with BMI and to ascertain whether using an extreme sampling design would identify the variants known to be associated with BMI in general populations.

Methodology/Principal Findings

From two large Danish cohorts we selected all extremely overweight young men and women (n = 2,633), and equal numbers of population-based controls (n = 2,740, drawn randomly from the same populations as the extremes, representing ∼212,000 individuals). We followed up novel (at the time of the study) association signals (p<0.001) from the discovery cohort in a genome-wide study of 5,846 Europeans, before attempting to replicate the most strongly associated 28 SNPs in an independent sample of Danish individuals (n = 20,917) and a population-based cohort of 15-year-old British adolescents (n = 2,418). Our discovery analysis identified SNPs at three loci known to be associated with BMI with genome-wide confidence (P<5×10−8; FTO, MC4R and FAIM2). We also found strong evidence of association at the known TMEM18, GNPDA2, SEC16B, TFAP2B, SH2B1 and KCTD15 loci (p<0.001), and nominal association (p<0.05) at a further 8 loci known to be associated with BMI. However, meta-analyses of our discovery and replication cohorts identified no novel associations.

Significance

Our results indicate that the detectable genetic variation associated with extreme overweight is very similar to that previously found for general BMI. This suggests that population-based study designs with enriched sampling of individuals with the extreme phenotype may be an efficient method for identifying common variants that influence quantitative traits and a valid alternative to genotyping all individuals in large population-based studies, which may require tens of thousands of subjects to achieve similar power.

Background

Single nucleotide polymorphisms (SNPs) are the most abundant source of genetic variation among individuals of a species. New genotyping technologies allow examining hundreds to thousands of SNPs in a single reaction for a wide range of applications such as genetic diversity analysis, linkage mapping, fine QTL mapping, association studies, marker-assisted or genome-wide selection. In this paper, we evaluated the potential of highly-multiplexed SNP genotyping for genetic mapping in maritime pine (Pinus pinaster Ait.), the main conifer used for commercial plantation in southwestern Europe.

Results

We designed a custom GoldenGate assay for 1,536 SNPs detected through the resequencing of gene fragments (707 in vitro SNPs/Indels) and from Sanger-derived Expressed Sequenced Tags assembled into a unigene set (829 in silico SNPs/Indels). Offspring from three-generation outbred (G2) and inbred (F2) pedigrees were genotyped. The success rate of the assay was 63.6% and 74.8% for in silico and in vitro SNPs, respectively. A genotyping error rate of 0.4% was further estimated from segregating data of SNPs belonging to the same gene. Overall, 394 SNPs were available for mapping. A total of 287 SNPs were integrated with previously mapped markers in the G2 parental maps, while 179 SNPs were localized on the map generated from the analysis of the F2 progeny. Based on 98 markers segregating in both pedigrees, we were able to generate a consensus map comprising 357 SNPs from 292 different loci. Finally, the analysis of sequence homology between mapped markers and their orthologs in a Pinus taeda linkage map, made it possible to align the 12 linkage groups of both species.

Conclusions

Our results show that the GoldenGate assay can be used successfully for high-throughput SNP genotyping in maritime pine, a conifer species that has a genome seven times the size of the human genome. This SNP-array will be extended thanks to recent sequencing effort using new generation sequencing technologies and will include SNPs from comparative orthologous sequences that were identified in the present study, providing a wider collection of anchor points for comparative genomics among the conifers.

Background

Alcohol has been shown to critically modulate cyclic adenosine-3′,5′ monophosphate (cAMP) signaling. A number of downstream effectors that respond to the cAMP signals (e.g., protein kinase A, cAMP response element binding protein) have, in turn, been examined in relation to alcohol consumption. These studies did not, however, delineate the point at which the actions of alcohol on the cAMP cascade might translate into differences in drinking behavior. To further understand the role of cAMP synthesis in alcohol drinking and dependence, we investigated a specific adenylyl cyclase isoform, adenylyl cyclase (AC) Type 7, whose activity is selectively enhanced by ethanol.

Methods

We measured alcohol consumption and preference in mice in which one copy of the Adcy7 gene was disrupted (Adcy7+/−). To demonstrate relevance of this gene for alcohol dependence in humans, we tested the association of polymorphisms in the ADCY7 gene with alcohol dependence in a sample of 1703 alcohol-dependent individuals and 1347 control subjects.

Results

We show that Adcy7+/− female mice have higher preference for alcohol than wild-type mice, whereas there is little difference in alcohol consumption or preference between Adcy7+/− male mice and wild-type control subjects. In the human sample, we found that single nucleotide polymorphisms in ADCY7 associate with alcohol dependence in women, and these markers are also associated with ADCY7 expression (messenger RNA) levels.

Conclusions

These findings implicate adenylyl cyclase Type 7 as a critical component of the molecular pathways contributing to alcohol drinking and the development of alcohol dependence.

Background

Usher syndrome (USH) combines sensorineural deafness with blindness. It is inherited in an autosomal recessive mode. Early diagnosis is critical for adapted educational and patient management choices, and for genetic counseling. To date, nine causative genes have been identified for the three clinical subtypes (USH1, USH2 and USH3). Current diagnostic strategies make use of a genotyping microarray that is based on the previously reported mutations. The purpose of this study was to design a more accurate molecular diagnosis tool.

Methods

We sequenced the 366 coding exons and flanking regions of the nine known USH genes, in 54 USH patients (27 USH1, 21 USH2 and 6 USH3).

Results

Biallelic mutations were detected in 39 patients (72%) and monoallelic mutations in an additional 10 patients (18.5%). In addition to biallelic mutations in one of the USH genes, presumably pathogenic mutations in another USH gene were detected in seven patients (13%), and another patient carried monoallelic mutations in three different USH genes. Notably, none of the USH3 patients carried detectable mutations in the only known USH3 gene, whereas they all carried mutations in USH2 genes. Most importantly, the currently used microarray would have detected only 30 of the 81 different mutations that we found, of which 39 (48%) were novel.

Conclusions

Based on these results, complete exon sequencing of the currently known USH genes stands as a definite improvement for molecular diagnosis of this disease, which is of utmost importance in the perspective of gene therapy.

Background

Rod-cone dystrophy, also known as retinitis pigmentosa (RP), and cone-rod dystrophy (CRD) are degenerative retinal dystrophies leading to blindness. To identify new genes responsible for these diseases, we have studied one large non consanguineous French family with autosomal dominant (ad) CRD.

Methods

Family members underwent detailed ophthalmological examination. Linkage analysis using microsatellite markers and a whole-genome SNP analysis with the use of Affymetrix 250 K SNP chips were performed. Five candidate genes within the candidate region were screened for mutations by direct sequencing.

Results

We first excluded the involvement of known adRP and adCRD genes in the family by genotyping and linkage analysis. Then, we undertook a whole-genome scan on 22 individuals in the family. The analysis revealed a 41.3-Mb locus on position 2q24.2-2q33.1. This locus was confirmed by linkage analysis with specific markers of this region. The maximum LOD score was 2.86 at θ = 0 for this locus. Five candidate genes, CERKL, BBS5, KLHL23, NEUROD1, and SF3B1 within this locus, were not mutated.

Conclusion

A novel locus for adCRD, named CORD12, has been mapped to chromosome 2q24.2-2q33.1 in a non consanguineous French family.

The only established genetic determinant of non-Mendelian forms of Alzheimer’s disease (AD) is the ε4 allele of the apolipoprotein E gene (APOE). Recently, it has been reported that the P86L polymorphism of the calcium homeostasis modulator 1 gene (CALHM1) is associated with the risk of developing AD. In order to independently assess this association, we performed a meta-analysis of 7,873 AD cases and 13,274 controls of Caucasian origin (from a total of 24 centres in Belgium, Finland, France, Italy, Spain, Sweden, the UK and the USA). Our results indicate that the CALHM1 P86L polymorphism is likely not a genetic determinant of AD but may modulate age at onset by interacting with the effect of the ε4 allele of the APOE gene.

Blood low-density lipoprotein (LDL) cholesterol, high-density lipoprotein (HDL) cholesterol and triglyceride levels are risk factors for cardiovascular disease. To dissect the polygenic basis of these traits, we conducted genome-wide association screens in 19,840 individuals and replication in up to 20,623 individuals. We identified 30 distinct loci associated with lipoprotein concentrations (each with P < 5 × 10-8), including 11 loci that reached genome-wide significance for the first time. The 11 newly defined loci include common variants associated with LDL cholesterol near ABCG8, MAFB, HNF1A and TIMD4; with HDL cholesterol near ANGPTL4, FADS1-FADS2-FADS3, HNF4A, LCAT, PLTP and TTC39B; and with triglycerides near AMAC1L2, FADS1-FADS2-FADS3 and PLTP. The proportion of individuals exceeding clinical cut points for high LDL cholesterol, low HDL cholesterol and high triglycerides varied according to an allelic dosage score (P < 10-15 for each trend). These results suggest that the cumulative effect of multiple common variants contributes to polygenic dyslipidemia.

Summary

Background

Information is scarce about the combined effects on breast cancer incidence of low-penetrance genetic susceptibility polymorphisms and environmental factors (reproductive, behavioural, and anthropometric risk factors for breast cancer). To test for evidence of gene–environment interactions, we compared genotypic relative risks for breast cancer across the other risk factors in a large UK prospective study.

Methods

We tested gene–environment interactions in 7610 women who developed breast cancer and 10 196 controls without the disease, studying the effects of 12 polymorphisms (FGFR2-rs2981582, TNRC9-rs3803662, 2q35-rs13387042, MAP3K1-rs889312, 8q24-rs13281615, 2p-rs4666451, 5p12-rs981782, CASP8-rs1045485, LSP1-rs3817198, 5q-rs30099, TGFB1-rs1982073, and ATM-rs1800054) in relation to prospectively collected information about ten established environmental risk factors (age at menarche, parity, age at first birth, breastfeeding, menopausal status, age at menopause, use of hormone replacement therapy, body-mass index, height, and alcohol consumption).

Findings

After allowance for multiple testing none of the 120 comparisons yielded significant evidence of a gene–environment interaction. By contrast with previous suggestions, there was little evidence that the genotypic relative risks were affected by use of hormone replacement therapy, either overall or for oestrogen-receptor-positive disease. Only one of the 12 polymorphisms was correlated with any of the ten other risk factors: carriers of the high-risk C allele of MAP3K1-rs889312 were significantly shorter than non-carriers (mean height 162·4 cm [95% CI 162·1–162·7] vs 163·1 cm [162·9–163·2]; p=0·01 after allowance for multiple testing).

Interpretation

Risks of breast cancer associated with low-penetrance susceptibility polymorphisms do not vary significantly with these ten established environmental risk factors.

Funding

Cancer Research UK and the UK Medical Research Council.

We carried out a genome-wide association study of lung cancer (3,259 cases and 4,159 controls), followed by replication in 2,899 cases and 5,573 controls. Two uncorrelated disease markers at 5p15.33, rs402710 and rs2736100 were detected by the genome-wide data (P = 2 × 10-7 and P = 4 × 10-6) and replicated by the independent study series (P = 7 × 10-5 and P = 0.016). The susceptibility region contains two genes, TERT and CLPTM1L, suggesting that one or both may have a role in lung cancer etiology.

Spondyloarthritis (SpA) is a chronic inflammatory disorder with a strong genetic predisposition dominated by the role of HLA-B27. However, the contribution of other genes to the disease susceptibility has been clearly demonstrated. We previously reported significant evidence of linkage of SpA to chromosome 9q31–34. The current study aimed to characterize this locus, named SPA2. First, we performed a fine linkage mapping of SPA2 (24 cM) with 28 microsatellite markers in 149 multiplex families, which allowed us to reduce the area of investigation to an 18 cM (13 Mb) locus delimited by the markers D9S279 and D9S112. Second, we constructed a linkage disequilibrium (LD) map of this region with 1,536 tag single-nucleotide polymorphisms (SNPs) in 136 families (263 patients). The association was assessed using a transmission disequilibrium test. One tag SNP, rs4979459, yielded a significant P-value (4.9×10−5). Third, we performed an extension association study with rs4979459 and 30 surrounding SNPs in LD with it, in 287 families (668 patients), and in a sample of 139 cases and 163 controls. Strong association was observed in both familial and case/control datasets for several SNPs. In the replication study, carried with 8 SNPs in an independent sample of 232 cases and 149 controls, one SNP, rs6478105, yielded a nominal P-value<3×10−2. Pooled case/control study (371 cases and 312 controls) as well as combined analysis of extension and replication data showed very significant association (P<5×10−4) for 6 of the 8 latter markers (rs7849556, rs10817669, rs10759734, rs6478105, rs10982396, and rs10733612). Finally, haplotype association investigations identified a strongly associated haplotype (P<8.8×10−5) consisting of these 6 SNPs and located in the direct vicinity of the TNFSF15 gene. In conclusion, we have identified within the SPA2 locus a haplotype strongly associated with predisposition to SpA which is located near to TNFSF15, one of the major candidate genes in this region.

Author Summary

Spondyloarthritis (SpA) is a common variety of articular inflammatory disorder characterized by axial and/or peripheral arthritis, frequently associated with extra-articular manifestations such as psoriasis, uveitis, and inflammatory bowel diseases (ulcerative colitis or Crohn's disease (CD)). SpA is a complex disorder with high heritability. The MHC class I HLA-B27 allele is a very strong risk factor for its development, but other genetic factors located outside the MHC also play a role in disease susceptibility. By a previous whole-genome linkage investigation, we have demonstrated that a region located on the chromosome 9q31–34 was involved in SpA susceptibility. The present study aimed to further characterize this locus. Using a stepwise linkage and association approach, we identified a haplotype spanning 6 single-nucleotide polymorphisms strongly associated with SpA and located in a genomic region paralogous to the MHC, near to the TNFSF15 gene. Interestingly, polymorphisms of this gene have previously been shown to be associated with CD. This original finding offers a new research track for the understanding of SpA pathophysiology, which is still poorly understood, as well as new hope for diagnostic and therapeutic innovation.

Several new risk factors for Crohn's disease have been identified in recent genome-wide association studies. To advance gene discovery further we have combined the data from three studies (a total of 3,230 cases and 4,829 controls) and performed replication in 3,664 independent cases with a mixture of population-based and family-based controls. The results strongly confirm 11 previously reported loci and provide genome-wide significant evidence for 21 new loci, including the regions containing STAT3, JAK2, ICOSLG, CDKAL1, and ITLN1. The expanded molecular understanding of the basis of disease offers promise for informed therapeutic development.