Genome-wide association studies in oesophageal adenocarcinoma and Barrett's oesophagus: a large-scale meta-analysis

Background: Esophageal adenocarcinoma (EA) is one of the fastest rising cancers in western countries. Barrett’s Esophagus (BE) is the premalignant precursor of EA. However, only a subset of BE patients develop EA, which complicates the clinical management in the absence of valid predictors. Genetic risk factors for BE and EA are incompletely understood. This study aimed to identify novel genetic risk factors for BE and EA.Methods: Within an international consortium of groups involved in the genetics of BE/EA, we performed the first meta-analysis of all genome-wide association studies (GWAS) available, involving 6,167 BE patients, 4,112 EA patients, and 17,159 representative controls, all of European ancestry, genotyped on Illumina high-density SNP-arrays, collected from four separate studies within North America, Europe, and Australia. Meta-analysis was conducted using the fixed-effects inverse variance-weighting approach. We used the standard genome-wide significant threshold of 5×10-8 for this study. We also conducted an association analysis following reweighting of loci using an approach that investigates annotation enrichment among the genome-wide significant loci. The entire GWAS-data set was also analyzed using bioinformatics approaches including functional annotation databases as well as gene-based and pathway-based methods in order to identify pathophysiologically relevant cellular pathways.Findings: We identified eight new associated risk loci for BE and EA, within or near the CFTR (rs17451754, P=4·8×10-10), MSRA (rs17749155, P=5·2×10-10), BLK (rs10108511, P=2·1×10-9), KHDRBS2 (rs62423175, P=3·0×10-9), TPPP/CEP72 (rs9918259, P=3·2×10-9), TMOD1 (rs7852462, P=1·5×10-8), SATB2 (rs139606545, P=2·0×10-8), and HTR3C/ABCC5 genes (rs9823696, P=1·6×10-8). A further novel risk locus at LPA (rs12207195, posteriori probability=0·925) was identified after re-weighting using significantly enriched annotations. This study thereby doubled the number of known risk loci. The strongest disease pathways identified (P<10-6) belong to muscle cell differentiation and to mesenchyme development/differentiation, which fit with current pathophysiological BE/EA concepts. To our knowledge, this study identified for the first time an EA-specific association (rs9823696, P=1·6×10-8) near HTR3C/ABCC5 which is independent of BE development (P=0·45).Interpretation: The identified disease loci and pathways reveal new insights into the etiology of BE and EA. Furthermore, the EA-specific association at HTR3C/ABCC5 may constitute a novel genetic marker for the prediction of transition from BE to EA. Mutations in CFTR, one of the new risk loci identified in this study, cause cystic fibrosis (CF), the most common recessive disorder in Europeans. Gastroesophageal reflux (GER) belongs to the phenotypic CF-spectrum and represents the main risk factor for BE/EA. Thus, the CFTR locus may trigger a common GER-mediated pathophysiology.

Pleiotropic Analysis of Lung Cancer and Blood Triglycerides

Epidemiologically related traits may share genetic risk factors, and pleiotropic analysis could identify individual loci associated with these traits. Because of their shared epidemiological associations, we conducted pleiotropic analysis of genome-wide association studies of lung cancer (12 160 lung cancer case patients and 16 838 control subjects) and cardiovascular disease risk factors (blood lipids from 188 577 subjects, type 2 diabetes from 148 821 subjects, body mass index from 123 865 subjects, and smoking phenotypes from 74 053 subjects). We found that 6p22.1 (rs6904596, ZNF184) was associated with both lung cancer (P = 5.50x10(-6)) and blood triglycerides (P = 1.39x10(-5)). We replicated the association in 6097 lung cancer case patients and 204 657 control subjects (P = 2.40 × 10(-4)) and in 71 113 subjects with triglycerides data (P = .01). rs6904596 reached genome-wide significance in lung cancer meta-analysis (odds ratio = 1.15, 95% confidence interval = 1.10 to 1.21 ,: Pcombined = 5.20x10(-9)). The large sample size provided by the lipid GWAS data and the shared genetic risk factors between the two traits contributed to the uncovering of a hitherto unidentified genetic locus for lung cancer.

HNF1B variants associate with promoter methylation and regulate gene networks activated in prostate and ovarian cancer

Two independent regions within HNF1B are consistently identified in prostate and ovarian cancer genome-wide association studies (GWAS); their functional roles are unclear. We link prostate cancer (PC) risk SNPs rs11649743 and rs3760511 with elevated HNF1B gene expression and allele-specific epigenetic silencing, and outline a mechanism by which common risk variants could effect functional changes that increase disease risk: functional assays suggest that HNF1B is a pro-differentiation factor that suppresses epithelial-to-mesenchymal transition (EMT) in unmethylated, healthy tissues. This tumor-suppressor activity is lost when HNF1B is silenced by promoter methylation in the progression to PC. Epigenetic inactivation of HNF1B in ovarian cancer also associates with known risk SNPs, with a similar impact on EMT. This represents one of the first comprehensive studies into the pleiotropic role of a GWAS-associated transcription factor across distinct cancer types, and is the first to describe a conserved role for a multi-cancer genetic risk factor.

Evidence for three genetic loci involved in both anorexia nervosa risk and variation of body mass index

The maintenance of normal body weight is disrupted in patients with anorexia nervosa (AN) for prolonged periods of time. Prior to the onset of AN, premorbid body mass index (BMI) spans the entire range from underweight to obese. After recovery, patients have reduced rates of overweight and obesity. As such, loci involved in body weight regulation may also be relevant for AN and vice versa. Our primary analysis comprised a cross-trait analysis of the 1000 single-nucleotide polymorphisms (SNPs) with the lowest P-values in a genome-wide association meta-analysis (GWAMA) of AN (GCAN) for evidence of association in the largest published GWAMA for BMI (GIANT). Subsequently we performed sex-stratified analyses for these 1000 SNPs. Functional ex vivo studies on four genes ensued. Lastly, a look-up of GWAMA-derived BMI-related loci was performed in the AN GWAMA. We detected significant associations (P-values <5 × 10(-5), Bonferroni-corrected P<0.05) for nine SNP alleles at three independent loci. Interestingly, all AN susceptibility alleles were consistently associated with increased BMI. None of the genes (chr. 10: CTBP2, chr. 19: CCNE1, chr. 2: CARF and NBEAL1; the latter is a region with high linkage disequilibrium) nearest to these SNPs has previously been associated with AN or obesity. Sex-stratified analyses revealed that the strongest BMI signal originated predominantly from females (chr. 10 rs1561589; Poverall: 2.47 × 10(-06)/Pfemales: 3.45 × 10(-07)/Pmales: 0.043). Functional ex vivo studies in mice revealed reduced hypothalamic expression of Ctbp2 and Nbeal1 after fasting. Hypothalamic expression of Ctbp2 was increased in diet-induced obese (DIO) mice as compared with age-matched lean controls. We observed no evidence for associations for the look-up of BMI-related loci in the AN GWAMA. A cross-trait analysis of AN and BMI loci revealed variants at three chromosomal loci with potential joint impact. The chromosome 10 locus is particularly promising given that the association with obesity was primarily driven by females. In addition, the detected altered hypothalamic expression patterns of Ctbp2 and Nbeal1 as a result of fasting and DIO implicate these genes in weight regulation.

The Genetic Landscape of Renal Complications in Type 1 Diabetes

Diabetes is the leading cause of end stage renal disease. Despite evidence for a substantial heritability of diabetic kidney disease, efforts to identify genetic susceptibility variants have had limited success. We extended previous efforts in three dimensions, examining a more comprehensive set of genetic variants in larger numbers of subjects with type 1 diabetes characterized for a wider range of cross-sectional diabetic kidney disease phenotypes. In 2,843 subjects, we estimated that the heritability of diabetic kidney disease was 35% ( p=6x10-3 ). Genome-wide association analysis and replication in 12,540 individuals identified no single variants reaching stringent levels of significance and, despite excellent power, provided little independent confirmation of previously published associated variants. Whole exome sequencing in 997 subjects failed to identify any large-effect coding alleles of lower frequency influencing the risk of diabetic kidney disease. However, sets of alleles increasing body mass index ( p=2.2×10-5) and the risk of type 2 diabetes (p=6.1x10-4 ) were associated with the risk of diabetic kidney disease. We also found genome-wide genetic correlation between diabetic kidney disease and failure at smoking cessation ( p=1.1×10-4 ). Pathway analysis implicated ascorbate and aldarate metabolism ( p=9×10-6), and pentose and glucuronate interconversions ( p=3×10-6) in pathogenesis of diabetic kidney disease. These data provide further evidence for the role of genetic factors influencing diabetic kidney disease in those with type 1 diabetes and highlight some key pathways that may be responsible. Altogether these results reveal important biology behind the major cause of kidney disease.  

Risk factors and biomarkers of age-related macular degeneration

A biomarker can be a substance or structure measured in body parts, fluids or products that can affect or predict disease incidence. As age-related macular degeneration (AMD) is the leading cause of blindness in the developed world, much research and effort has been invested in the identification of different biomarkers to predict disease incidence, identify at risk individuals, elucidate causative pathophysiological etiologies, guide screening, monitoring and treatment parameters, and predict disease outcomes. To date, a host of genetic, environmental, proteomic, and cellular targets have been identified as both risk factors and potential biomarkers for AMD. Despite this, their use has been confined to research settings and has not yet crossed into the clinical arena. A greater understanding of these factors and their use as potential biomarkers for AMD can guide future research and clinical practice. This article will discuss known risk factors and novel, potential biomarkers of AMD in addition to their application in both academic and clinical settings.

Novel Genetic Determinants of Diabetic Kidney Disease

Diabetic kidney disease (DKD) is a devastating diabetes complication, with known heritability not fully revealed by previous genetics studies. We performed the largest genome-wide association study of type 1 DKD to date, in a 13-cohort consortium of 15,590 individuals of European ancestry genotyped on the Illumina HumanCoreExome Beadchip, which allows exploration of coding variation in addition to genomic markers.

As prior work has shown that different characterizations of the DKD phenotype highlight distinct genetic associations, we investigated a spectrum of DKD definitions based on proteinuria and renal function criteria. Controls were DKD-free after a minimum of 15 years diabetes duration; cases had diabetes for at least 10 years prior to DKD diagnosis. We also performed a quantitative trait analysis of estimated glomerular filtration rate in all participants.

Our top finding was a missense mutation in COL4A3, rs55703767 (Asp326Tyr); the minor allele is common in Europeans (20%) and East Asians (13%) but not Africans (2%). This SNP had a genome-wide significant association with traditionally defined DKD (macroalbuminuria or end-stage renal disease [ESRD], (OR= 0.79, P=1.9×10-9), and a suggestive association with macroalbuminuria (OR= 0.79, P=1.6×10-6) and ESRD (OR= 0.79, P=4.5×10-5) individually. Though its PolyPhen score is 0.3 (benign), this SNP has been implicated as a splice site disruptor.

The COL4A3 gene encodes the alpha 3 subunit of Type IV collagen, the major structural component of basement membranes. Pathogenic mutations in COL4A3 have been identified in thin basement membrane nephropathy, familial focal segmental glomerulosclerosis, and Alport syndrome. A proxy (r2=0.6) for rs55703767 had no significant associations in the CKDGen consortium, suggesting its pathogenicity occurs solely in the setting of hyperglycemia.

By significantly increasing sample size we have discovered a novel locus underlying DKD risk, paving the way for better understanding of pathology, prevention, and treatment.

A Rare Loss-of-Function Variant in Scavenger Receptor Class B Type I (SCARB1) Raises HDL Cholesterol and Increases Risk of Coronary Disease

Scavenger receptor BI (SR-BI) is the major receptor for high-density lipoprotein (HDL)
cholesterol (HDL-C). In humans, high amounts of HDL-C in plasma are associated with a
lower risk of coronary heart disease (CHD). Mice that have depleted Scarb1 (SR-BI
knockout mice) have markedly elevated HDL-C levels but, paradoxically, increased
atherosclerosis. The impact of SR-BI on HDL metabolism and CHD risk in humans remains
unclear. Through targeted sequencing of coding regions of lipid-modifying genes in 328
individuals with extremely high plasma HDL-C levels, we identified a homozygote for a lossof-function
variant, in which leucine replaces proline 376 (P376L), in SCARB1, the gene
encoding SR-BI. The P376L variant impairs posttranslational processing of SR-BI and
abrogates selective HDL cholesterol uptake in transfected cells, in hepatocyte-like cells
derived from induced pluripotent stem cells from the homozygous subject, and in mice.
Large population-based studies revealed that subjects who are heterozygous carriers of
the P376L variant have significantly increased levels of plasma HDL-C. P376L carriers have
a profound HDL-related phenotype and an increased risk of CHD (odds ratio = 1.79, which is
statistically significant).

Prevalence and risk of Down syndrome in monozygotic and dizygotic multiple pregnancies in Europe: implications for prenatal screening

OBJECTIVE: To determine risk of Down syndrome (DS) in multiple relative to singleton pregnancies, and compare prenatal diagnosis rates and pregnancy outcome.

DESIGN: Population-based prevalence study based on EUROCAT congenital anomaly registries.

SETTING: Eight European countries.

POPULATION: 14.8 million births 1990-2009; 2.89% multiple births.

METHODS: DS cases included livebirths, fetal deaths from 20 weeks, and terminations of pregnancy for fetal anomaly (TOPFA). Zygosity is inferred from like/unlike sex for birth denominators, and from concordance for DS cases.

MAIN OUTCOME MEASURES: Relative risk (RR) of DS per fetus/baby from multiple versus singleton pregnancies and per pregnancy in monozygotic/dizygotic versus singleton pregnancies. Proportion of prenatally diagnosed and pregnancy outcome.

STATISTICAL ANALYSIS: Poisson and logistic regression stratified for maternal age, country and time.

RESULTS: Overall, the adjusted (adj) RR of DS for fetus/babies from multiple versus singleton pregnancies was 0.58 (95% CI 0.53-0.62), similar for all maternal ages except for mothers over 44, for whom it was considerably lower. In 8.7% of twin pairs affected by DS, both co-twins were diagnosed with the condition. The adjRR of DS for monozygotic versus singleton pregnancies was 0.34 (95% CI 0.25-0.44) and for dizygotic versus singleton pregnancies 1.34 (95% CI 1.23-1.46). DS fetuses from multiple births were less likely to be prenatally diagnosed than singletons (adjOR 0.62 [95% CI 0.50-0.78]) and following diagnosis less likely to be TOPFA (adjOR 0.40 [95% CI 0.27-0.59]).

CONCLUSIONS: The risk of DS per fetus/baby is lower in multiple than singleton pregnancies. These estimates can be used for genetic counselling and prenatal screening.

The genetic basis of psychological traits in infancy: Implications for understanding the causes of developmental psychopathology

This chapter reviews genetic studies that have aimed to identify genes influencing psychological traits in infancy (from birth to age 12 months), and considers how this research informs us about the causes of developmental psychopathology. Specifically, this chapter systematically reviews findings from studies that associated common genetic variants with individual variation in infants’ attention, temperament and behaviour, and attachment disorganisation. DRD4 and 5-HTTLPR genes were the most frequently studied candidate genes. Possibly the most coherent set of results relates to the L-DRD4 genotype, which is significantly associated with infant attention, temperament, and attachment style. Research in infant genetics has been strengthened by a careful focus on uniform age ranges within studies, by several longitudinal studies, and by exploration of gene-environment interactions between genes and maternal characteristics. However there is also considerable inconsistency in results in this field and possible reasons for this are discussed. The chapter outlines the main genetic methods that have been used and what new genetic approaches such as polygenic risk scoring could offer infant genetics. Recent findings suggest that some traits during infancy predict individual differences in developmental psychopathology in childhood. It is argued that infant genetic research has considerable potential for the identification of populations at risk for psychopathology in later life, and this remains an area open for future research.

Rare coding variants in the phospholipase D3 gene confer risk for Alzheimer's disease.

Genome-wide association studies (GWAS) have identified several risk variants for late-onset Alzheimer's disease (LOAD)1, 2. These common variants have replicable but small effects on LOAD risk and generally do not have obvious functional effects. Low-frequency coding variants, not detected by GWAS, are predicted to include functional variants with larger effects on risk. To identify low-frequency coding variants with large effects on LOAD risk, we carried out whole-exome sequencing (WES) in 14 large LOAD families and follow-up analyses of the candidate variants in several large LOAD case–control data sets. A rare variant in PLD3 (phospholipase D3; Val232Met) segregated with disease status in two independent families and doubled risk for Alzheimer’s disease in seven independent case–control series with a total of more than 11,000 cases and controls of European descent. Gene-based burden analyses in 4,387 cases and controls of European descent and 302 African American cases and controls, with complete sequence data for PLD3, reveal that several variants in this gene increase risk for Alzheimer’s disease in both populations. PLD3 is highly expressed in brain regions that are vulnerable to Alzheimer’s disease pathology, including hippocampus and cortex, and is expressed at significantly lower levels in neurons from Alzheimer’s disease brains compared to control brains. Overexpression of PLD3 leads to a significant decrease in intracellular amyloid-β precursor protein (APP) and extracellular Aβ42 and Aβ40 (the 42- and 40-residue isoforms of the amyloid-β peptide), and knockdown of PLD3 leads to a significant increase in extracellular Aβ42 and Aβ40. Together, our genetic and functional data indicate that carriers of PLD3 coding variants have a twofold increased risk for LOAD and that PLD3 influences APP processing. This study provides an example of how densely affected families may help to identify rare variants with large effects on risk for disease or other complex traits.

Alleles that increase risk for type 2 diabetes mellitus are not associated with increased risk for Alzheimer's disease

Although epidemiological studies suggest that type 2 diabetes mellitus (T2DM) increases the risk of late-onset Alzheimer's disease (LOAD), the biological basis of this relationship is not well understood. The aim of this study was to examine the genetic comorbidity between the 2 disorders and to investigate whether genetic liability to T2DM, estimated by a genotype risk scores based on T2DM associated loci, is associated with increased risk of LOAD. This study was performed in 2 stages. In stage 1, we combined genotypes for the top 15 T2DM-associated polymorphisms drawn from approximately 3000 individuals (1349 cases and 1351 control subjects) with extracted and/or imputed data from 6 genome-wide studies (>10,000 individuals; 4507 cases, 2183 controls, 4989 population controls) to form a genotype risk score and examined if this was associated with increased LOAD risk in a combined meta-analysis. In stage 2, we investigated the association of LOAD with an expanded T2DM score made of 45 well-established variants drawn from the 6 genome-wide studies. Results were combined in a meta-analysis. Both stage 1 and stage 2 T2DM risk scores were not associated with LOAD risk (odds ratio = 0.988; 95% confidence interval, 0.972-1.004; p = 0.144 and odds ratio = 0.993; 95% confidence interval, 0.983-1.003; p = 0.149 per allele, respectively). Contrary to expectation, genotype risk scores based on established T2DM candidates were not associated with increased risk of LOAD. The observed epidemiological associations between T2DM and LOAD could therefore be a consequence of secondary disease processes, pleiotropic mechanisms, and/or common environmental risk factors. Future work should focus on well-characterized longitudinal cohorts with extensive phenotypic and genetic data relevant to both LOAD and T2DM.