Mari Nelis

Systematic identification of trans eQTLs as putative drivers of known disease associations

Identifying the downstream effects of disease-associated SNPs is challenging. To help overcome this problem, we performed expression quantitative trait locus (eQTL) meta-analysis in non-transformed peripheral blood samples from 5,311 individuals with replication in 2,775 individuals. We identified and replicated trans eQTLs for 233 SNPs (reflecting 103 independent loci) that were previously associated with complex traits at genome-wide significance. Some of these SNPs affect multiple genes in trans that are known to be altered in individuals with disease: rs4917014, previously associated with systemic lupus erythematosus (SLE), altered gene expression of C1QB and five type I interferon response genes, both hallmarks of SLE. DeepSAGE RNA sequencing showed that rs4917014 strongly alters the 3′ UTR levels of IKZF1 in cis, and chromatin immunoprecipitation and sequencing analysis of the trans-regulated genes implicated IKZF1 as the causal gene. Variants associated with cholesterol metabolism and type 1 diabetes showed similar phenomena, indicating that large-scale eQTL mapping provides insight into the downstream effects of many trait-associated variants.

A Genome-wide Association Study of Lung Cancer Identifies a Region of Chromosome 5p15 Associated with Risk for Adenocarcinoma

Three genetic loci for lung cancer risk have been identified by genome-wide association studies (GWAS), but inherited susceptibility to specific histologic types of lung cancer is not well established. We conducted a GWAS of lung cancer and its major histologic types, genotyping 515,922 single-nucleotide polymorphisms (SNPs) in 5739 lung cancer cases and 5848 controls from one population-based case-control study and three cohort studies. Results were combined with summary data from ten additional studies, for a total of 13,300 cases and 19,666 controls of European descent. Four studies also provided histology data for replication, resulting in 3333 adenocarcinomas (AD), 2589 squamous cell carcinomas (SQ), and 1418 small cell carcinomas (SC). In analyses by histology, rs2736100 (TERT), on chromosome 5p15.33, was associated with risk of adenocarcinoma (odds ratio [OR]=1.23, 95% confidence interval [CI]=1.13-1.33, p=3.02x10(-7)), but not with other histologic types (OR=1.01, p=0.84 and OR=1.00, p=0.93 for SQ and SC, respectively). This finding was confirmed in each replication study and overall meta-analysis (OR=1.24, 95% CI=1.17-1.31, p=3.74x10(-14) for AD; OR=0.99, p=0.69 and OR=0.97, p=0.48 for SQ and SC, respectively). Other previously reported association signals on 15q25 and 6p21 were also refined, but no additional loci reached genome-wide significance. In conclusion, a lung cancer GWAS identified a distinct hereditary contribution to adenocarcinoma.

A new highly penetrant form of obesity due to deletions on chromosome 16p11.2

Obesity has become a major worldwide challenge to public health, owing to an interaction between the Western ‘obesogenic’ environment and a strong genetic contribution. Recent extensive genome-wide association studies (GWASs) have identified numerous single nucleotide polymorphisms associated with obesity, but these loci together account for only a small fraction of the known heritable component. Thus, the ‘common disease, common variant’ hypothesis is increasingly coming under challenge. Here we report a highly penetrant form of obesity, initially observed in 31 subjects who were heterozygous for deletions of at least 593 kilobases at 16p11.2 and whose ascertainment included cognitive deficits. Nineteen similar deletions were identified from GWAS data in 16,053 individuals from eight European cohorts. These deletions were absent from healthy non-obese controls and accounted for 0.7% of our morbid obesity cases (body mass index (BMI) ≥ 40 kg m-2 or BMI standard deviation score ≥ 4; P = 6.4 × 10-8, odds ratio 43.0), demonstrating the potential importance in common disease of rare variants with strong effects. This highlights a promising strategy for identifying missing heritability in obesity and other complex traits: cohorts with extreme phenotypes are likely to be enriched for rare variants, thereby improving power for their discovery. Subsequent analysis of the loci so identified may well reveal additional rare variants that further contribute to the missing heritability, as recently reported for SIM1 (ref. 3). The most productive approach may therefore be to combine the ‘power of the extreme’ in small, well-phenotyped cohorts, with targeted follow-up in case-control and population cohorts.

Common variants in KCNN3 are associated with lone atrial fibrillation

Atrial fibrillation (AF) is the most common sustained arrhythmia. Previous studies have identified several genetic loci associated with typical AF. We sought to identify common genetic variants underlying lone AF. This condition affects a subset of individuals without overt heart disease and with an increased heritability of AF. We report a meta-analysis of genome-wide association studies conducted using 1,335 individuals with lone AF (cases) and 12,844 unaffected individuals (referents). Cases were obtained from the German AF Network, Heart and Vascular Health Study, the Atherosclerosis Risk in Communities Study, the Cleveland Clinic and Massachusetts General Hospital. We identified an association on chromosome 1q21 to lone AF (rs13376333, adjusted odds ratio = 1.56; P = 6.3 × 10−12), and we replicated this association in two independent cohorts with lone AF (overall combined odds ratio = 1.52, 95% CI 1.40–1.64; P = 1.83 × 10−21). rs13376333 is intronic to KCNN3, which encodes a potassium channel protein involved in atrial repolarization.

Genetic Structure of Europeans: A View from the North–East

Using principal component (PC) analysis, we studied the genetic constitution of 3,112 individuals from Europe as portrayed by more than 270,000 single nucleotide polymorphisms (SNPs) genotyped with the Illumina Infinium platform. In cohorts where the sample size was >100, one hundred randomly chosen samples were used for analysis to minimize the sample size effect, resulting in a total of 1,564 samples. This analysis revealed that the genetic structure of the European population correlates closely with geography. The first two PCs highlight the genetic diversity corresponding to the northwest to southeast gradient and position the populations according to their approximate geographic origin. The resulting genetic map forms a triangular structure with a) Finland, b) the Baltic region, Poland and Western Russia, and c) Italy as its vertexes, and with d) Central- and Western Europe in its centre. Inter- and intra- population genetic differences were quantified by the inflation factor lambda (λ) (ranging from 1.00 to 4.21), fixation index (Fst) (ranging from 0.000 to 0.023), and by the number of markers exhibiting significant allele frequency differences in pair-wise population comparisons. The estimated lambda was used to assess the real diminishing impact to association statistics when two distinct populations are merged directly in an analysis. When the PC analysis was confined to the 1,019 Estonian individuals (0.1% of the Estonian population), a fine structure emerged that correlated with the geography of individual counties. With at least two cohorts available from several countries, genetic substructures were investigated in Czech, Finnish, German, Estonian and Italian populations. Together with previously published data, our results allow the creation of a comprehensive European genetic map that will greatly facilitate inter-population genetic studies including genome wide association studies (GWAS).

Rare variants of large effect in BRCA2 and CHEK2 affect risk of lung cancer

We conducted imputation to the 1000 Genomes Project of four genome-wide association studies of lung cancer in populations of European ancestry (11,348 cases and 15,861 controls) and genotyped an additional 10,246 cases and 38,295 controls for follow-up. We identified large-effect genome-wide associations for squamous lung cancer with the rare variants BRCA2 p.Lys3326X (rs11571833, odds ratio (OR) = 2.47, P = 4.74 × 10−20) and CHEK2 p.Ile157Thr (rs17879961, OR = 0.38, P = 1.27 × 10−13). We also showed an association between common variation at 3q28 (TP63, rs13314271, OR = 1.13, P = 7.22 × 10−10) and lung adenocarcinoma that had been previously reported only in Asians. These findings provide further evidence for inherited genetic susceptibility to lung cancer and its biological basis. Additionally, our analysis demonstrates that imputation can identify rare disease-causing variants with substantive effects on cancer risk from preexisting genome-wide association study data.

Influence of common genetic variation on lung cancer risk: meta-analysis of 14 900 cases and 29 485 controls

Recent genome-wide association studies (GWASs) have identified common genetic variants at 5p15.33, 6p21–6p22 and 15q25.1 associated with lung cancer risk. Several other genetic regions including variants of CHEK2 (22q12), TP53BP1 (15q15) and RAD52 (12p13) have been demonstrated to influence lung cancer risk in candidate- or pathway-based analyses. To identify novel risk variants for lung cancer, we performed a meta-analysis of 16 GWASs, totaling 14 900 cases and 29 485 controls of European descent. Our data provided increased support for previously identified risk loci at 5p15 (P = 7.2 × 10−16), 6p21 (P = 2.3 × 10−14) and 15q25 (P = 2.2 × 10−63). Furthermore, we demonstrated histology-specific effects for 5p15, 6p21 and 12p13 loci but not for the 15q25 region. Subgroup analysis also identified a novel disease locus for squamous cell carcinoma at 9p21 (CDKN2A/p16INK4A/p14ARF/CDKN2B/p15INK4B/ANRIL; rs1333040, P = 3.0 × 10−7) which was replicated in a series of 5415 Han Chinese (P = 0.03; combined analysis, P = 2.3 × 10−8). This large analysis provides additional evidence for the role of inherited genetic susceptibility to lung cancer and insight into biological differences in the development of the different histological types of lung cancer.

Shared and Unique Components of Human Population Structure and Genome-Wide Signals of Positive Selection in South Asia

South Asia harbors one of the highest levels genetic diversity in Eurasia, which could be interpreted as a result of its long-term large effective population size and of admixture during its complex demographic history. In contrast to Pakistani populations, populations of Indian origin have been underrepresented in previous genomic scans of positive selection and population structure. Here we report data for more than 600,000 SNP markers genotyped in 142 samples from 30 ethnic groups in India. Combining our results with other available genome-wide data, we show that Indian populations are characterized by two major ancestry components, one of which is spread at comparable frequency and haplotype diversity in populations of South and West Asia and the Caucasus. The second component is more restricted to South Asia and accounts for more than 50% of the ancestry in Indian populations. Haplotype diversity associated with these South Asian ancestry components is significantly higher than that of the components dominating the West Eurasian ancestry palette. Modeling of the observed haplotype diversities suggests that both Indian ancestry components are older than the purported Indo-Aryan invasion 3,500 YBP. Consistent with the results of pairwise genetic distances among world regions, Indians share more ancestry signals with West than with East Eurasians. However, compared to Pakistani populations, a higher proportion of their genes show regionally specific signals of high haplotype homozygosity. Among such candidates of positive selection in India are MSTN and DOK5, both of which have potential implications in lipid metabolism and the etiology of type 2 diabetes.

Genetic associations of 115 polymorphisms with cancers of the upper aerodigestive tract across 10 european countries: The ARCAGE project

Cancers of the upper aerodigestive tract (UADT) include malignant tumors of the oral cavity, pharynx, larynx, and esophagus and account for 6.4% of all new cancers in Europe. In the context of a multicenter case-control study conducted in 14 centers within 10 European countries and comprising 1,511 cases and 1,457 controls (ARCAGE study), 115 single nucleotide polymorphisms (SNP) from 62 a priori-selected genes were studied in relation to UADT cancer. We found 11 SNPs that were statistically associated with UADT cancers overall (5.75 expected). Considering the possibility of false-positive results, we focused on SNPs in CYP2A6, MDM2, tumor necrosis factor (TNF), and gene amplified in squamous cell carcinoma 1 (GASC1), for which low P values for trend (P trend<0.01) were observed in the main effects analyses of UADT cancer overall or by subsite. The rare variant of CYP2A6 -47A>C (rs28399433), a phase I metabolism gene, was associated with reduced UADT cancer risk (P trend=0.01). Three SNPs in the MDM2 gene, involved in cell cycle control, were associated with UADT cancer. MDM2 IVS5+1285A>G (rs3730536) showed a strong codominant effect (P trend=0.007). The rare variants of two SNPs in the TNF gene were associated with a decreased risk; for TNF IVS1+123G>A (rs1800610), the P trend was 0.007. Variants in two SNPs of GASC1 were found to be strongly associated with increased UADT cancer risk (for both, P trend=0.008). This study is the largest genetic epidemiologic study on UADT cancers in Europe. Our analysis points to potentially relevant genes in various pathways.

Genome-wide analysis of BMI in adolescents and young adults reveals additional insight into the effects of genetic loci over the life course

Genetic loci for body mass index (BMI) in adolescence and young adulthood, a period of high risk for weight gain, are understudied, yet may yield important insight into the etiology of obesity and early intervention. To identify novel genetic loci and examine the influence of known loci on BMI during this critical time period in late adolescence and early adulthood, we performed a two-stage meta-analysis using 14 genome-wide association studies in populations of European ancestry with data on BMI between ages 16 and 25 in up to 29 880 individuals. We identified seven independent loci (P < 5.0 × 10⁻⁸) near FTO (P = 3.72 × 10⁻²³), TMEM18 (P = 3.24 × 10⁻¹⁷), MC4R (P = 4.41 × 10⁻¹⁷), TNNI3K (P = 4.32 × 10⁻¹¹), SEC16B (P = 6.24 × 10⁻⁹), GNPDA2 (P = 1.11 × 10⁻⁸) and POMC (P = 4.94 × 10⁻⁸) as well as a potential secondary signal at the POMC locus (rs2118404, P = 2.4 × 10⁻⁵ after conditioning on the established single-nucleotide polymorphism at this locus) in adolescents and young adults. To evaluate the impact of the established genetic loci on BMI at these young ages, we examined differences between the effect sizes of 32 published BMI loci in European adult populations (aged 18-90) and those observed in our adolescent and young adult meta-analysis. Four loci (near PRKD1, TNNI3K, SEC16B and CADM2) had larger effects and one locus (near SH2B1) had a smaller effect on BMI during adolescence and young adulthood compared with older adults (P < 0.05). These results suggest that genetic loci for BMI can vary in their effects across the life course, underlying the importance of evaluating BMI at different ages.