Genome-wide Association Studies in Ancestrally Diverse Populations: Opportunities, Methods, Pitfalls, and Recommendations.

Genome-wide association studies (GWASs) have focused primarily on populations of European descent, but it is essential that diverse populations become better represented. Increasing diversity among study participants will advance our understanding of genetic architecture in all populations and ensure that genetic research is broadly applicable. To facilitate and promote research in multi-ancestry and admixed cohorts, we outline key methodological considerations and highlight opportunities, challenges, solutions, and areas in need of development. Despite the perception that analyzing genetic data from diverse populations is difficult, it is scientifically and ethically imperative, and there is an expanding analytical toolbox to do it well.

Mendelian randomization analyses implicate biogenesis of translation machinery in human aging.

Reduced provision of protein translation machinery promotes healthy aging in a number of animal models. In humans, however, inborn impairments in translation machinery are a known cause of several developmental disorders, collectively termed ribosomopathies. Here, we use casual inference approaches in genetic epidemiology to investigate whether adult, tissue-specific biogenesis of translation machinery drives human aging. We assess naturally occurring variation in the expression of genes encoding subunits specific to the two RNA polymerases (Pols) that transcribe ribosomal and transfer RNAs, namely Pol I and III, and the variation in expression of ribosomal protein (RP) genes, using Mendelian randomization. We find each causally associated with human longevity (ß = -0.15 ± 0.047, P = 9.6 × 10-4, q = 0.015; ß = -0.13 ± 0.040, P = 1.4 × 10-3, q = 0.023; ß = -0.048 ± 0.016, P = 3.5 × 10-3, q = 0.056, respectively), and this does not appear to be mediated by altered susceptibility to a single disease. We find that reduced expression of Pol III, RPs, or Pol I promotes longevity from different organs, namely visceral adipose, liver, and skeletal muscle, echoing the tissue specificity of ribosomopathies. Our study shows the utility of leveraging genetic variation in expression to elucidate how essential cellular processes impact human aging. The findings extend the evolutionary conservation of protein synthesis as a critical process that drives animal aging to include humans.

Prevalence and correlates of common mental disorders among participants of the Uganda Genome Resource: Opportunities for psychiatric genetics research.

Genetics research has potential to alleviate the burden of mental disorders in low- and middle-income-countries through identification of new mechanistic pathways which can lead to efficacious drugs or new drug targets. However, there is currently limited genetics data from Africa. The Uganda Genome Resource provides opportunity for psychiatric genetics research among underrepresented people from Africa. We aimed at determining the prevalence and correlates of major depressive disorder (MDD), suicidality, post-traumatic stress disorder (PTSD), alcohol abuse, generalised anxiety disorder (GAD) and probable attention-deficit hyperactivity disorder (ADHD) among participants of the Uganda Genome Resource. Standardised tools assessed for each mental disorder. Prevalence of each disorder was calculated with 95% confidence intervals. Multivariate logistic regression models evaluated the association between each mental disorder and associated demographic and clinical factors. Among 985 participants, prevalence of the disorders were: current MDD 19.3%, life-time MDD 23.3%, suicidality 10.6%, PTSD 3.1%, alcohol abuse 5.7%, GAD 12.9% and probable ADHD 9.2%. This is the first study to determine the prevalence of probable ADHD among adult Ugandans from a general population. We found significant association between sex and alcohol abuse (adjusted odds ratio [AOR] = 0.26 [0.14,0.45], p < 0.001) and GAD (AOR = 1.78 [1.09,2.49], p = 0.019) respectively. We also found significant association between body mass index and suicidality (AOR = 0.85 [0.73,0.99], p = 0.041), alcohol abuse (AOR = 0.86 [0.78,0.94], p = 0.003) and GAD (AOR = 0.93 [0.87,0.98], p = 0.008) respectively. We also found a significant association between high blood pressure and life-time MDD (AOR = 2.87 [1.08,7.66], p = 0.035) and probable ADHD (AOR = 1.99 [1.00,3.97], p = 0.050) respectively. We also found a statistically significant association between tobacco smoking and alcohol abuse (AOR = 3.2 [1.56,6.67], p = 0.002). We also found ever been married to be a risk factor for probable ADHD (AOR = 2.12 [0.88,5.14], p = 0.049). The Uganda Genome Resource presents opportunity for psychiatric genetics research among underrepresented people from Africa.

Widespread genomic influences on phenotype in Dravet syndrome, a 'monogenic' condition.

Dravet syndrome is an archetypal rare severe epilepsy, considered 'monogenic', typically caused by loss-of-function SCN1A variants. Despite a recognizable core phenotype, its marked phenotypic heterogeneity is incompletely explained by differences in the causal SCN1A variant or clinical factors. In 34 adults with SCN1A-related Dravet syndrome, we show additional genomic variation beyond SCN1A contributes to phenotype and its diversity, with an excess of rare variants in epilepsy-related genes as a set and examples of blended phenotypes, including one individual with an ultra-rare DEPDC5 variant and focal cortical dysplasia. The polygenic risk score for intelligence was lower, and for longevity, higher, in Dravet syndrome than in epilepsy controls. The causal, major-effect, SCN1A variant may need to act against a broadly compromised genomic background to generate the full Dravet syndrome phenotype, whilst genomic resilience may help to ameliorate the risk of premature mortality in adult Dravet syndrome survivors.

Mendelian randomisation study of body composition and depression in people of East Asian ancestry highlights potential setting-specific causality.

BACKGROUND: Extensive evidence links higher body mass index (BMI) to higher odds of depression in people of European ancestry. However, our understanding of the relationship across different settings and ancestries is limited. Here, we test the relationship between body composition and depression in people of East Asian ancestry. METHODS: Multiple Mendelian randomisation (MR) methods were used to test the relationship between (a) BMI and (b) waist-hip ratio (WHR) with depression. Firstly, we performed two-sample MR using genetic summary statistics from a recent genome-wide association study (GWAS) of depression (with 15,771 cases and 178,777 controls) in people of East Asian ancestry. We selected 838 single nucleotide polymorphisms (SNPs) correlated with BMI and 263 SNPs correlated with WHR as genetic instrumental variables to estimate the causal effect of BMI and WHR on depression using the inverse-variance weighted (IVW) method. We repeated these analyses stratifying by home location status: China versus UK or USA. Secondly, we performed one-sample MR in the China Kadoorie Biobank (CKB) in 100,377 participants. This allowed us to test the relationship separately in (a) males and females and (b) urban and rural dwellers. We also examined (c) the linearity of the BMI-depression relationship. RESULTS: Both MR analyses provided evidence that higher BMI was associated with lower odds of depression. For example, a genetically-instrumented 1-SD higher BMI in the CKB was associated with lower odds of depressive symptoms [OR: 0.77, 95% CI: 0.63, 0.95]. There was evidence of differences according to place of residence. Using the IVW method, higher BMI was associated with lower odds of depression in people of East Asian ancestry living in China but there was no evidence for an association in people of East Asian ancestry living in the USA or UK. Furthermore, higher genetic BMI was associated with differential effects in urban and rural dwellers within China. CONCLUSIONS: This study provides the first MR evidence for an inverse relationship between BMI and depression in people of East Asian ancestry. This contrasts with previous findings in European populations and therefore the public health response to obesity and depression is likely to need to differ based on sociocultural factors for example, ancestry and place of residence. This highlights the importance of setting-specific causality when using genetic causal inference approaches and data from diverse populations to test hypotheses. This is especially important when the relationship tested is not purely biological and may involve sociocultural factors.

Is auditory processing measured by the N100 an endophenotype for psychosis? A family study and a meta-analysis.

BACKGROUND: The N100, an early auditory event-related potential, has been found to be altered in patients with psychosis. However, it is unclear if the N100 is a psychosis endophenotype that is also altered in the relatives of patients. METHODS: We conducted a family study using the auditory oddball paradigm to compare the N100 amplitude and latency across 243 patients with psychosis, 86 unaffected relatives, and 194 controls. We then conducted a systematic review and a random-effects meta-analysis pooling our results and 14 previously published family studies. We compared data from a total of 999 patients, 1192 relatives, and 1253 controls in order to investigate the evidence and degree of N100 differences. RESULTS: In our family study, patients showed reduced N100 amplitudes and prolonged N100 latencies compared to controls, but no significant differences were found between unaffected relatives and controls. The meta-analysis revealed a significant reduction of the N100 amplitude and delay of the N100 latency in both patients with psychosis (standardized mean difference [s.m.d.] = -0.48 for N100 amplitude and s.m.d. = 0.43 for N100 latency) and their relatives (s.m.d. = - 0.19 for N100 amplitude and s.m.d. = 0.33 for N100 latency). However, only the N100 latency changes in relatives remained significant when excluding studies with affected relatives. CONCLUSIONS: N100 changes, especially prolonged N100 latencies, are present in both patients with psychosis and their relatives, making the N100 a promising endophenotype for psychosis. Such changes in the N100 may reflect changes in early auditory processing underlying the etiology of psychosis.

Type 2 diabetes risks and determinants in second-generation migrants and mixed ethnicity people of South Asian and African Caribbean descent in the UK.

AIMS/HYPOTHESIS: Excess risks of type 2 diabetes in UK South Asians (SA) and African Caribbeans (AC) compared with Europeans remain unexplained. We studied risks and determinants of type 2 diabetes in first- and second-generation (born in the UK) migrants, and in those of mixed ethnicity. METHODS: Data from the UK Biobank, a population-based cohort of ~500,000 participants aged 40-69 at recruitment, were used. Type 2 diabetes was assigned using self-report and HbA1c. Ethnicity was both self-reported and genetically assigned using admixture level scores. European, mixed European/South Asian (MixESA), mixed European/African Caribbean (MixEAC), SA and AC groups were analysed, matched for age and sex to enable comparison. In the frames of this cross-sectional study, we compared type 2 diabetes in second- vs first-generation migrants, and mixed ethnicity vs non-mixed groups. Risks and explanations were analysed using logistic regression and mediation analysis, respectively. RESULTS: Type 2 diabetes prevalence was markedly elevated in SA (599/3317 = 18%) and AC (534/4180 = 13%) compared with Europeans (140/3324 = 4%). Prevalence was lower in second- vs first-generation SA (124/1115 = 11% vs 155/1115 = 14%) and AC (163/2200 = 7% vs 227/2200 = 10%). Favourable adiposity (i.e. lower waist/hip ratio or BMI) contributed to lower risk in second-generation migrants. Type 2 diabetes in mixed populations (MixESA: 52/831 = 6%, MixEAC: 70/1045 = 7%) was lower than in comparator ethnic groups (SA: 18%, AC: 13%) and higher than in Europeans (4%). Greater socioeconomic deprivation accounted for 17% and 42% of the excess type 2 diabetes risk in MixESA and MixEAC compared with Europeans, respectively. Replacing self-reported with genetically assigned ethnicity corroborated the mixed ethnicity analysis. CONCLUSIONS/INTERPRETATION: Type 2 diabetes risks in second-generation SA and AC migrants are a fifth lower than in first-generation migrants. Mixed ethnicity risks were markedly lower than SA and AC groups, though remaining higher than in Europeans. Distribution of environmental risk factors, largely obesity and socioeconomic status, appears to play a key role in accounting for ethnic differences in type 2 diabetes risk.

Integrating polygenic risk scores in the prediction of type 2 diabetes risk and subtypes in British Pakistanis and Bangladeshis: A population-based cohort study.

BACKGROUND: Type 2 diabetes (T2D) is highly prevalent in British South Asians, yet they are underrepresented in research. Genes & Health (G&H) is a large, population study of British Pakistanis and Bangladeshis (BPB) comprising genomic and routine health data. We assessed the extent to which genetic risk for T2D is shared between BPB and European populations (EUR). We then investigated whether the integration of a polygenic risk score (PRS) for T2D with an existing risk tool (QDiabetes) could improve prediction of incident disease and the characterisation of disease subtypes. METHODS AND FINDINGS: In this observational cohort study, we assessed whether common genetic loci associated with T2D in EUR individuals were replicated in 22,490 BPB individuals in G&H. We replicated fewer loci in G&H (n = 76/338, 22%) than would be expected given power if all EUR-ascertained loci were transferable (n = 101, 30%; p = 0.001). Of the 27 transferable loci that were powered to interrogate this, only 9 showed evidence of shared causal variants. We constructed a T2D PRS and combined it with a clinical risk instrument (QDiabetes) in a novel, integrated risk tool (IRT) to assess risk of incident diabetes. To assess model performance, we compared categorical net reclassification index (NRI) versus QDiabetes alone. In 13,648 patients free from T2D followed up for 10 years, NRI was 3.2% for IRT versus QDiabetes (95% confidence interval (CI): 2.0% to 4.4%). IRT performed best in reclassification of individuals aged less than 40 years deemed low risk by QDiabetes alone (NRI 5.6%, 95% CI 3.6% to 7.6%), who tended to be free from comorbidities and slim. After adjustment for QDiabetes score, PRS was independently associated with progression to T2D after gestational diabetes (hazard ratio (HR) per SD of PRS 1.23, 95% CI 1.05 to 1.42, p = 0.028). Using cluster analysis of clinical features at diabetes diagnosis, we replicated previously reported disease subgroups, including Mild Age-Related, Mild Obesity-related, and Insulin-Resistant Diabetes, and showed that PRS distribution differs between subgroups (p = 0.002). Integrating PRS in this cluster analysis revealed a Probable Severe Insulin Deficient Diabetes (pSIDD) subgroup, despite the absence of clinical measures of insulin secretion or resistance. We also observed differences in rates of progression to micro- and macrovascular complications between subgroups after adjustment for confounders. Study limitations include the absence of an external replication cohort and the potential biases arising from missing or incorrect routine health data. CONCLUSIONS: Our analysis of the transferability of T2D loci between EUR and BPB indicates the need for larger, multiancestry studies to better characterise the genetic contribution to disease and its varied aetiology. We show that a T2D PRS optimised for this high-risk BPB population has potential clinical application in BPB, improving the identification of T2D risk (especially in the young) on top of an established clinical risk algorithm and aiding identification of subgroups at diagnosis, which may help future efforts to stratify care and treatment of the disease.

Polygenic risk scores indicate extreme ages at onset of breast cancer in female BRCA1/2 pathogenic variant carriers.

BACKGROUND: Clinical management of women carrying a germline pathogenic variant (PV) in the BRCA1/2 genes demands for accurate age-dependent estimators of breast cancer (BC) risks, which were found to be affected by a variety of intrinsic and extrinsic factors. Here we assess the contribution of polygenic risk scores (PRSs) to the occurrence of extreme phenotypes with respect to age at onset, namely, primary BC diagnosis before the age of 35 years (early diagnosis, ED) and cancer-free survival until the age of 60 years (late/no diagnosis, LD) in female BRCA1/2 PV carriers. METHODS: Overall, estrogen receptor (ER)-positive, and ER-negative BC PRSs as developed by Kuchenbaecker et al. for BC risk discrimination in female BRCA1/2 PV carriers were employed for PRS computation in a curated sample of 295 women of European descent carrying PVs in the BRCA1 (n=183) or the BRCA2 gene (n=112), and did either fulfill the ED criteria (n=162, mean age at diagnosis: 28.3 years, range: 20 to 34 years) or the LD criteria (n=133). Binomial logistic regression was applied to assess the association of standardized PRSs with either ED or LD under adjustment for patient recruitment criteria for germline testing and localization of BRCA1/2 PVs in the corresponding BC or ovarian cancer (OC) cluster regions. RESULTS: For BRCA1 PV carriers, the standardized overall BC PRS displayed the strongest association with ED (odds ratio (OR) = 1.62; 95% confidence interval (CI): 1.16-2.31, p<0.01). Additionally, statistically significant associations of selection for the patient recruitment criteria for germline testing and localization of pathogenic PVs outside the BRCA1 OC cluster region with ED were observed. For BRCA2 PV carriers, the standardized PRS for ER-negative BC displayed the strongest association (OR = 2.27, 95% CI: 1.45-3.78, p<0.001). CONCLUSIONS: PRSs contribute to the development of extreme phenotypes of female BRCA1/2 PV carriers with respect to age at primary BC diagnosis. Construction of optimized PRS SNP sets for BC risk stratification in BRCA1/2 PV carriers should be the task of future studies with larger, well-defined study samples. Furthermore, our results provide further evidence, that localization of PVs in BC/OC cluster regions might be considered in BC risk calculations for unaffected BRCA1/2 PV carriers.

Genetic copy number variants, cognition and psychosis: a meta-analysis and a family study.

The burden of large and rare copy number genetic variants (CNVs) as well as certain specific CNVs increase the risk of developing schizophrenia. Several cognitive measures are purported schizophrenia endophenotypes and may represent an intermediate point between genetics and the illness. This paper investigates the influence of CNVs on cognition. We conducted a systematic review and meta-analysis of the literature exploring the effect of CNV burden on general intelligence. We included ten primary studies with a total of 18,847 participants and found no evidence of association. In a new psychosis family study, we investigated the effects of CNVs on specific cognitive abilities. We examined the burden of large and rare CNVs (>200 kb, <1% MAF) as well as known schizophrenia-associated CNVs in patients with psychotic disorders, their unaffected relatives and controls (N = 3428) from the Psychosis Endophenotypes International Consortium (PEIC). The carriers of specific schizophrenia-associated CNVs showed poorer performance than non-carriers in immediate (P = 0.0036) and delayed (P = 0.0115) verbal recall. We found suggestive evidence that carriers of schizophrenia-associated CNVs had poorer block design performance (P = 0.0307). We do not find any association between CNV burden and cognition. Our findings show that the known high-risk CNVs are not only associated with schizophrenia and other neurodevelopmental disorders, but are also a contributing factor to impairment in cognitive domains such as memory and perceptual reasoning, and act as intermediate biomarkers of disease risk.

Very low-depth whole-genome sequencing in complex trait association studies.

MOTIVATION: Very low-depth sequencing has been proposed as a cost-effective approach to capture low-frequency and rare variation in complex trait association studies. However, a full characterization of the genotype quality and association power for very low-depth sequencing designs is still lacking. RESULTS: We perform cohort-wide whole-genome sequencing (WGS) at low depth in 1239 individuals (990 at 1× depth and 249 at 4× depth) from an isolated population, and establish a robust pipeline for calling and imputing very low-depth WGS genotypes from standard bioinformatics tools. Using genotyping chip, whole-exome sequencing (75× depth) and high-depth (22×) WGS data in the same samples, we examine in detail the sensitivity of this approach, and show that imputed 1× WGS recapitulates 95.2% of variants found by imputed GWAS with an average minor allele concordance of 97% for common and low-frequency variants. In our study, 1× further allowed the discovery of 140 844 true low-frequency variants with 73% genotype concordance when compared to high-depth WGS data. Finally, using association results for 57 quantitative traits, we show that very low-depth WGS is an efficient alternative to imputed GWAS chip designs, allowing the discovery of up to twice as many true association signals than the classical imputed GWAS design. AVAILABILITY AND IMPLEMENTATION: The HELIC genotype and WGS datasets have been deposited to the European Genome-phenome Archive (https://www.ebi.ac.uk/ega/home): EGAD00010000518; EGAD00010000522; EGAD00010000610; EGAD00001001636, EGAD00001001637. The peakplotter software is available at https://github.com/wtsi-team144/peakplotter, the transformPhenotype app can be downloaded at https://github.com/wtsi-team144/transformPhenotype. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Combination therapy as a potential risk factor for the development of type 2 diabetes in patients with schizophrenia: the GOMAP study.

BACKGROUND: Schizophrenia (SCZ) is associated with increased risk of type 2 diabetes (T2D). The potential diabetogenic effect of concomitant application of psychotropic treatment classes in patients with SCZ has not yet been evaluated. The overarching goal of the Genetic Overlap between Metabolic and Psychiatric disease (GOMAP) study is to assess the effect of pharmacological, anthropometric, lifestyle and clinical measurements, helping elucidate the mechanisms underlying the aetiology of T2D. METHODS: The GOMAP case-control study (Genetic Overlap between Metabolic and Psychiatric disease) includes hospitalized patients with SCZ, some of whom have T2D. We enrolled 1653 patients with SCZ; 611 with T2D and 1042 patients without T2D. This is the first study of SCZ and T2D comorbidity at this scale in the Greek population. We retrieved detailed information on first- and second-generation antipsychotics (FGA, SGA), antidepressants and mood stabilizers, applied as monotherapy, 2-drug combination, or as 3- or more drug combination. We assessed the effects of psychotropic medication, body mass index, duration of schizophrenia, number of hospitalizations and physical activity on risk of T2D. Using logistic regression, we calculated crude and adjusted odds ratios (OR) to identify associations between demographic factors and the psychiatric medications. RESULTS: Patients with SCZ on a combination of at least three different classes of psychiatric drugs had a higher risk of T2D [OR 1.81 (95% CI 1.22-2.69); p = 0.003] compared to FGA alone therapy, after adjustment for age, BMI, sex, duration of SCZ and number of hospitalizations. We did not find evidence for an association of SGA use or the combination of drugs belonging to two different classes of psychiatric medications with increased risk of T2D [1.27 (0.84-1.93), p = 0.259 and 0.98 (0.71-1.35), p = 0.885, respectively] compared to FGA use. CONCLUSIONS: We find an increased risk of T2D in patients with SCZ who take a combination of at least three different psychotropic medication classes compared to patients whose medication consists only of one or two classes of drugs.

Cell-type-specific enrichment of risk-associated regulatory elements at ovarian cancer susceptibility loci.

Understanding the regulatory landscape of the human genome is a central question in complex trait genetics. Most single-nucleotide polymorphisms (SNPs) associated with cancer risk lie in non-protein-coding regions, implicating regulatory DNA elements as functional targets of susceptibility variants. Here, we describe genome-wide annotation of regions of open chromatin and histone modification in fallopian tube and ovarian surface epithelial cells (FTSECs, OSECs), the debated cellular origins of high-grade serous ovarian cancers (HGSOCs) and in endometriosis epithelial cells (EECs), the likely precursor of clear cell ovarian carcinomas (CCOCs). The regulatory architecture of these cell types was compared with normal human mammary epithelial cells and LNCaP prostate cancer cells. We observed similar positional patterns of global enhancer signatures across the three different ovarian cancer precursor cell types, and evidence of tissue-specific regulatory signatures compared to non-gynecological cell types. We found significant enrichment for risk-associated SNPs intersecting regulatory biofeatures at 17 known HGSOC susceptibility loci in FTSECs (P = 3.8 × 10(-30)), OSECs (P = 2.4 × 10(-23)) and HMECs (P = 6.7 × 10(-15)) but not for EECs (P = 0.45) or LNCaP cells (P = 0.88). Hierarchical clustering of risk SNPs conditioned on the six different cell types indicates FTSECs and OSECs are highly related (96% of samples using multi-scale bootstrapping) suggesting both cell types may be precursors of HGSOC. These data represent the first description of regulatory catalogues of normal precursor cells for different ovarian cancer subtypes, and provide unique insights into the tissue specific regulatory variation with respect to the likely functional targets of germline genetic susceptibility variants for ovarian cancer.

DNA glycosylases involved in base excision repair may be associated with cancer risk in BRCA1 and BRCA2 mutation carriers.

Single Nucleotide Polymorphisms (SNPs) in genes involved in the DNA Base Excision Repair (BER) pathway could be associated with cancer risk in carriers of mutations in the high-penetrance susceptibility genes BRCA1 and BRCA2, given the relation of synthetic lethality that exists between one of the components of the BER pathway, PARP1 (poly ADP ribose polymerase), and both BRCA1 and BRCA2. In the present study, we have performed a comprehensive analysis of 18 genes involved in BER using a tagging SNP approach in a large series of BRCA1 and BRCA2 mutation carriers. 144 SNPs were analyzed in a two stage study involving 23,463 carriers from the CIMBA consortium (the Consortium of Investigators of Modifiers of BRCA1 and BRCA2). Eleven SNPs showed evidence of association with breast and/or ovarian cancer at p<0.05 in the combined analysis. Four of the five genes for which strongest evidence of association was observed were DNA glycosylases. The strongest evidence was for rs1466785 in the NEIL2 (endonuclease VIII-like 2) gene (HR: 1.09, 95% CI (1.03-1.16), p = 2.7 × 10(-3)) for association with breast cancer risk in BRCA2 mutation carriers, and rs2304277 in the OGG1 (8-guanine DNA glycosylase) gene, with ovarian cancer risk in BRCA1 mutation carriers (HR: 1.12 95%CI: 1.03-1.21, p = 4.8 × 10(-3)). DNA glycosylases involved in the first steps of the BER pathway may be associated with cancer risk in BRCA1/2 mutation carriers and should be more comprehensively studied.

Risks of Breast, Ovarian, and Contralateral Breast Cancer for BRCA1 and BRCA2 Mutation Carriers.

Importance: The clinical management of BRCA1 and BRCA2 mutation carriers requires accurate, prospective cancer risk estimates. Objectives: To estimate age-specific risks of breast, ovarian, and contralateral breast cancer for mutation carriers and to evaluate risk modification by family cancer history and mutation location. Design, Setting, and Participants: Prospective cohort study of 6036 BRCA1 and 3820 BRCA2 female carriers (5046 unaffected and 4810 with breast or ovarian cancer or both at baseline) recruited in 1997-2011 through the International BRCA1/2 Carrier Cohort Study, the Breast Cancer Family Registry and the Kathleen Cuningham Foundation Consortium for Research into Familial Breast Cancer, with ascertainment through family clinics (94%) and population-based studies (6%). The majority were from large national studies in the United Kingdom (EMBRACE), the Netherlands (HEBON), and France (GENEPSO). Follow-up ended December 2013; median follow-up was 5 years. Exposures: BRCA1/2 mutations, family cancer history, and mutation location. Main Outcomes and Measures: Annual incidences, standardized incidence ratios, and cumulative risks of breast, ovarian, and contralateral breast cancer. Results: Among 3886 women (median age, 38 years; interquartile range [IQR], 30-46 years) eligible for the breast cancer analysis, 5066 women (median age, 38 years; IQR, 31-47 years) eligible for the ovarian cancer analysis, and 2213 women (median age, 47 years; IQR, 40-55 years) eligible for the contralateral breast cancer analysis, 426 were diagnosed with breast cancer, 109 with ovarian cancer, and 245 with contralateral breast cancer during follow-up. The cumulative breast cancer risk to age 80 years was 72% (95% CI, 65%-79%) for BRCA1 and 69% (95% CI, 61%-77%) for BRCA2 carriers. Breast cancer incidences increased rapidly in early adulthood until ages 30 to 40 years for BRCA1 and until ages 40 to 50 years for BRCA2 carriers, then remained at a similar, constant incidence (20-30 per 1000 person-years) until age 80 years. The cumulative ovarian cancer risk to age 80 years was 44% (95% CI, 36%-53%) for BRCA1 and 17% (95% CI, 11%-25%) for BRCA2 carriers. For contralateral breast cancer, the cumulative risk 20 years after breast cancer diagnosis was 40% (95% CI, 35%-45%) for BRCA1 and 26% (95% CI, 20%-33%) for BRCA2 carriers (hazard ratio [HR] for comparing BRCA2 vs BRCA1, 0.62; 95% CI, 0.47-0.82; P=.001 for difference). Breast cancer risk increased with increasing number of first- and second-degree relatives diagnosed as having breast cancer for both BRCA1 (HR for ≥2 vs 0 affected relatives, 1.99; 95% CI, 1.41-2.82; P<.001 for trend) and BRCA2 carriers (HR, 1.91; 95% CI, 1.08-3.37; P=.02 for trend). Breast cancer risk was higher if mutations were located outside vs within the regions bounded by positions c.2282-c.4071 in BRCA1 (HR, 1.46; 95% CI, 1.11-1.93; P=.007) and c.2831-c.6401 in BRCA2 (HR, 1.93; 95% CI, 1.36-2.74; P<.001). Conclusions and Relevance: These findings provide estimates of cancer risk based on BRCA1 and BRCA2 mutation carrier status using prospective data collection and demonstrate the potential importance of family history and mutation location in risk assessment.

Identification of independent association signals and putative functional variants for breast cancer risk through fine-scale mapping of the 12p11 locus.

BACKGROUND: Multiple recent genome-wide association studies (GWAS) have identified a single nucleotide polymorphism (SNP), rs10771399, at 12p11 that is associated with breast cancer risk. METHOD: We performed a fine-scale mapping study of a 700 kb region including 441 genotyped and more than 1300 imputed genetic variants in 48,155 cases and 43,612 controls of European descent, 6269 cases and 6624 controls of East Asian descent and 1116 cases and 932 controls of African descent in the Breast Cancer Association Consortium (BCAC; http://bcac.ccge.medschl.cam.ac.uk/ ), and in 15,252 BRCA1 mutation carriers in the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA). Stepwise regression analyses were performed to identify independent association signals. Data from the Encyclopedia of DNA Elements project (ENCODE) and the Cancer Genome Atlas (TCGA) were used for functional annotation. RESULTS: Analysis of data from European descendants found evidence for four independent association signals at 12p11, represented by rs7297051 (odds ratio (OR) = 1.09, 95 % confidence interval (CI) = 1.06-1.12; P = 3 × 10(-9)), rs805510 (OR = 1.08, 95 % CI = 1.04-1.12, P = 2 × 10(-5)), and rs1871152 (OR = 1.04, 95 % CI = 1.02-1.06; P = 2 × 10(-4)) identified in the general populations, and rs113824616 (P = 7 × 10(-5)) identified in the meta-analysis of BCAC ER-negative cases and BRCA1 mutation carriers. SNPs rs7297051, rs805510 and rs113824616 were also associated with breast cancer risk at P < 0.05 in East Asians, but none of the associations were statistically significant in African descendants. Multiple candidate functional variants are located in putative enhancer sequences. Chromatin interaction data suggested that PTHLH was the likely target gene of these enhancers. Of the six variants with the strongest evidence of potential functionality, rs11049453 was statistically significantly associated with the expression of PTHLH and its nearby gene CCDC91 at P < 0.05. CONCLUSION: This study identified four independent association signals at 12p11 and revealed potentially functional variants, providing additional insights into the underlying biological mechanism(s) for the association observed between variants at 12p11 and breast cancer risk.

HMG-coenzyme A reductase inhibition, type 2 diabetes, and bodyweight: evidence from genetic analysis and randomised trials.

BACKGROUND: Statins increase the risk of new-onset type 2 diabetes mellitus. We aimed to assess whether this increase in risk is a consequence of inhibition of 3-hydroxy-3-methylglutaryl-CoA reductase (HMGCR), the intended drug target. METHODS: We used single nucleotide polymorphisms in the HMGCR gene, rs17238484 (for the main analysis) and rs12916 (for a subsidiary analysis) as proxies for HMGCR inhibition by statins. We examined associations of these variants with plasma lipid, glucose, and insulin concentrations; bodyweight; waist circumference; and prevalent and incident type 2 diabetes. Study-specific effect estimates per copy of each LDL-lowering allele were pooled by meta-analysis. These findings were compared with a meta-analysis of new-onset type 2 diabetes and bodyweight change data from randomised trials of statin drugs. The effects of statins in each randomised trial were assessed using meta-analysis. FINDINGS: Data were available for up to 223 463 individuals from 43 genetic studies. Each additional rs17238484-G allele was associated with a mean 0·06 mmol/L (95% CI 0·05-0·07) lower LDL cholesterol and higher body weight (0·30 kg, 0·18-0·43), waist circumference (0·32 cm, 0·16-0·47), plasma insulin concentration (1·62%, 0·53-2·72), and plasma glucose concentration (0·23%, 0·02-0·44). The rs12916 SNP had similar effects on LDL cholesterol, bodyweight, and waist circumference. The rs17238484-G allele seemed to be associated with higher risk of type 2 diabetes (odds ratio [OR] per allele 1·02, 95% CI 1·00-1·05); the rs12916-T allele association was consistent (1·06, 1·03-1·09). In 129 170 individuals in randomised trials, statins lowered LDL cholesterol by 0·92 mmol/L (95% CI 0·18-1·67) at 1-year of follow-up, increased bodyweight by 0·24 kg (95% CI 0·10-0·38 in all trials; 0·33 kg, 95% CI 0·24-0·42 in placebo or standard care controlled trials and -0·15 kg, 95% CI -0·39 to 0·08 in intensive-dose vs moderate-dose trials) at a mean of 4·2 years (range 1·9-6·7) of follow-up, and increased the odds of new-onset type 2 diabetes (OR 1·12, 95% CI 1·06-1·18 in all trials; 1·11, 95% CI 1·03-1·20 in placebo or standard care controlled trials and 1·12, 95% CI 1·04-1·22 in intensive-dose vs moderate dose trials). INTERPRETATION: The increased risk of type 2 diabetes noted with statins is at least partially explained by HMGCR inhibition. FUNDING: The funding sources are cited at the end of the paper.

Identification of a BRCA2-specific modifier locus at 6p24 related to breast cancer risk.

Common genetic variants contribute to the observed variation in breast cancer risk for BRCA2 mutation carriers; those known to date have all been found through population-based genome-wide association studies (GWAS). To comprehensively identify breast cancer risk modifying loci for BRCA2 mutation carriers, we conducted a deep replication of an ongoing GWAS discovery study. Using the ranked P-values of the breast cancer associations with the imputed genotype of 1.4 M SNPs, 19,029 SNPs were selected and designed for inclusion on a custom Illumina array that included a total of 211,155 SNPs as part of a multi-consortial project. DNA samples from 3,881 breast cancer affected and 4,330 unaffected BRCA2 mutation carriers from 47 studies belonging to the Consortium of Investigators of Modifiers of BRCA1/2 were genotyped and available for analysis. We replicated previously reported breast cancer susceptibility alleles in these BRCA2 mutation carriers and for several regions (including FGFR2, MAP3K1, CDKN2A/B, and PTHLH) identified SNPs that have stronger evidence of association than those previously published. We also identified a novel susceptibility allele at 6p24 that was inversely associated with risk in BRCA2 mutation carriers (rs9348512; per allele HR = 0.85, 95% CI 0.80-0.90, P = 3.9 × 10(-8)). This SNP was not associated with breast cancer risk either in the general population or in BRCA1 mutation carriers. The locus lies within a region containing TFAP2A, which encodes a transcriptional activation protein that interacts with several tumor suppressor genes. This report identifies the first breast cancer risk locus specific to a BRCA2 mutation background. This comprehensive update of novel and previously reported breast cancer susceptibility loci contributes to the establishment of a panel of SNPs that modify breast cancer risk in BRCA2 mutation carriers. This panel may have clinical utility for women with BRCA2 mutations weighing options for medical prevention of breast cancer.

Genome-wide association study in BRCA1 mutation carriers identifies novel loci associated with breast and ovarian cancer risk.

BRCA1-associated breast and ovarian cancer risks can be modified by common genetic variants. To identify further cancer risk-modifying loci, we performed a multi-stage GWAS of 11,705 BRCA1 carriers (of whom 5,920 were diagnosed with breast and 1,839 were diagnosed with ovarian cancer), with a further replication in an additional sample of 2,646 BRCA1 carriers. We identified a novel breast cancer risk modifier locus at 1q32 for BRCA1 carriers (rs2290854, P = 2.7 × 10(-8), HR = 1.14, 95% CI: 1.09-1.20). In addition, we identified two novel ovarian cancer risk modifier loci: 17q21.31 (rs17631303, P = 1.4 × 10(-8), HR = 1.27, 95% CI: 1.17-1.38) and 4q32.3 (rs4691139, P = 3.4 × 10(-8), HR = 1.20, 95% CI: 1.17-1.38). The 4q32.3 locus was not associated with ovarian cancer risk in the general population or BRCA2 carriers, suggesting a BRCA1-specific association. The 17q21.31 locus was also associated with ovarian cancer risk in 8,211 BRCA2 carriers (P = 2×10(-4)). These loci may lead to an improved understanding of the etiology of breast and ovarian tumors in BRCA1 carriers. Based on the joint distribution of the known BRCA1 breast cancer risk-modifying loci, we estimated that the breast cancer lifetime risks for the 5% of BRCA1 carriers at lowest risk are 28%-50% compared to 81%-100% for the 5% at highest risk. Similarly, based on the known ovarian cancer risk-modifying loci, the 5% of BRCA1 carriers at lowest risk have an estimated lifetime risk of developing ovarian cancer of 28% or lower, whereas the 5% at highest risk will have a risk of 63% or higher. Such differences in risk may have important implications for risk prediction and clinical management for BRCA1 carriers.