Polygenic risk score portability for common diseases across genetically diverse populations.

BACKGROUND: Polygenic risk scores (PRS) derived from European individuals have reduced portability across global populations, limiting their clinical implementation at worldwide scale. Here, we investigate the performance of a wide range of PRS models across four ancestry groups (Africans, Europeans, East Asians, and South Asians) for 14 conditions of high-medical interest. METHODS: To select the best-performing model per trait, we first compared PRS performances for publicly available scores, and constructed new models using different methods (LDpred2, PRS-CSx and SNPnet). We used 285 K European individuals from the UK Biobank (UKBB) for training and 18 K, including diverse ancestries, for testing. We then evaluated PRS portability for the best models in Europeans and compared their accuracies with respect to the best PRS per ancestry. Finally, we validated the selected PRS models using an independent set of 8,417 individuals from Biobank of the Americas-Genomelink (BbofA-GL); and performed a PRS-Phewas. RESULTS: We confirmed a decay in PRS performances relative to Europeans when the evaluation was conducted using the best-PRS model for Europeans (51.3% for South Asians, 46.6% for East Asians and 39.4% for Africans). We observed an improvement in the PRS performances when specifically selecting ancestry specific PRS models (phenotype variance increase: 1.62 for Africans, 1.40 for South Asians and 0.96 for East Asians). Additionally, when we selected the optimal model conditional on ancestry for CAD, HDL-C and LDL-C, hypertension, hypothyroidism and T2D, PRS performance for studied populations was more comparable to what was observed in Europeans. Finally, we were able to independently validate tested models for Europeans, and conducted a PRS-Phewas, identifying cross-trait interplay between cardiometabolic conditions, and between immune-mediated components. CONCLUSION: Our work comprehensively evaluated PRS accuracy across a wide range of phenotypes, reducing the uncertainty with respect to which PRS model to choose and in which ancestry group. This evaluation has let us identify specific conditions where implementing risk-prioritization strategies could have practical utility across diverse ancestral groups, contributing to democratizing the implementation of PRS.

Predictive Utility of Polygenic Risk Scores for Coronary Heart Disease in Three Major Racial and Ethnic Groups.

Because polygenic risk scores (PRSs) for coronary heart disease (CHD) are derived from mainly European ancestry (EA) cohorts, their validity in African ancestry (AA) and Hispanic ethnicity (HE) individuals is unclear. We investigated associations of "restricted" and genome-wide PRSs with CHD in three major racial and ethnic groups in the U.S. The eMERGE cohort (mean age 48 ± 14 years, 58% female) included 45,645 EA, 7,597 AA, and 2,493 HE individuals. We assessed two restricted PRSs (PRSTikkanen and PRSTada; 28 and 50 variants, respectively) and two genome-wide PRSs (PRSmetaGRS and PRSLDPred; 1.7 M and 6.6 M variants, respectively) derived from EA cohorts. Over a median follow-up of 11.1 years, 2,652 incident CHD events occurred. Hazard and odds ratios for the association of PRSs with CHD were similar in EA and HE cohorts but lower in AA cohorts. Genome-wide PRSs were more strongly associated with CHD than restricted PRSs were. PRSmetaGRS, the best performing PRS, was associated with CHD in all three cohorts; hazard ratios (95% CI) per 1 SD increase were 1.53 (1.46-1.60), 1.53 (1.23-1.90), and 1.27 (1.13-1.43) for incident CHD in EA, HE, and AA individuals, respectively. The hazard ratios were comparable in the EA and HE cohorts (pinteraction = 0.77) but were significantly attenuated in AA individuals (pinteraction= 2.9 × 10-3). These results highlight the potential clinical utility of PRSs for CHD as well as the need to assemble diverse cohorts to generate ancestry- and ethnicity PRSs.

Validating and automating learning of cardiometabolic polygenic risk scores from direct-to-consumer genetic and phenotypic data: implications for scaling precision health research.

INTRODUCTION: A major challenge to enabling precision health at a global scale is the bias between those who enroll in state sponsored genomic research and those suffering from chronic disease. More than 30 million people have been genotyped by direct-to-consumer (DTC) companies such as 23andMe, Ancestry DNA, and MyHeritage, providing a potential mechanism for democratizing access to medical interventions and thus catalyzing improvements in patient outcomes as the cost of data acquisition drops. However, much of these data are sequestered in the initial provider network, without the ability for the scientific community to either access or validate. Here, we present a novel geno-pheno platform that integrates heterogeneous data sources and applies learnings to common chronic disease conditions including Type 2 diabetes (T2D) and hypertension. METHODS: We collected genotyped data from a novel DTC platform where participants upload their genotype data files and were invited to answer general health questionnaires regarding cardiometabolic traits over a period of 6 months. Quality control, imputation, and genome-wide association studies were performed on this dataset, and polygenic risk scores were built in a case-control setting using the BASIL algorithm. RESULTS: We collected data on N = 4,550 (389 cases / 4,161 controls) who reported being affected or previously affected for T2D and N = 4,528 (1,027 cases / 3,501 controls) for hypertension. We identified 164 out of 272 variants showing identical effect direction to previously reported genome-significant findings in Europeans. Performance metric of the PRS models was AUC = 0.68, which is comparable to previously published PRS models obtained with larger datasets including clinical biomarkers. DISCUSSION: DTC platforms have the potential of inverting research models of genome sequencing and phenotypic data acquisition. Quality control (QC) mechanisms proved to successfully enable traditional GWAS and PRS analyses. The direct participation of individuals has shown the potential to generate rich datasets enabling the creation of PRS cardiometabolic models. More importantly, federated learning of PRS from reuse of DTC data provides a mechanism for scaling precision health care delivery beyond the small number of countries who can afford to finance these efforts directly. CONCLUSIONS: The genetics of T2D and hypertension have been studied extensively in controlled datasets, and various polygenic risk scores (PRS) have been developed. We developed predictive tools for both phenotypes trained with heterogeneous genotypic and phenotypic data generated outside of the clinical environment and show that our methods can recapitulate prior findings with fidelity. From these observations, we conclude that it is possible to leverage DTC genetic repositories to identify individuals at risk of debilitating diseases based on their unique genetic landscape so that informed, timely clinical interventions can be incorporated.

A Functional Polymorphism Downstream of Vitamin A Regulator Gene CYP26B1 Is Associated with Hand Osteoarthritis.

While genetic analyses have revealed ~100 risk loci associated with osteoarthritis (OA), only eight have been linked to hand OA. Besides, these studies were performed in predominantly European and Caucasian ancestries. Here, we conducted a genome-wide association study in the Han Chinese population to identify genetic variations associated with the disease. We recruited a total of 1136 individuals (n = 420 hand OA-affected; n = 716 unaffected control subjects) of Han Chinese ancestry. We carried out genotyping using Axiom Asia Precisi on Medicine Research Array, and we employed the RegulomeDB database and RoadMap DNase I Hypersensitivity Sites annotations to further narrow down our potential candidate variants. Genetic variants identified were tested in the Geisinger's hand OA cohort selected from the Geisinger MyCode community health initiative (MyCode®). We also performed a luciferase reporter assay to confirm the potential impact of top candidate single-nucleotide polymorphisms (SNPs) on hand OA. We identified six associated SNPs (p-value = 6.76 × 10-7-7.31 × 10-6) clustered at 2p13.2 downstream of the CYP26B1 gene. The strongest association signal identified was rs883313 (p-value = 6.76 × 10-7, odds ratio (OR) = 1.76), followed by rs12713768 (p-value = 1.36 × 10-6, OR = 1.74), near or within the enhancer region closest to the CYP26B1 gene. Our findings showed that the major risk-conferring CC haplotype of SNPs rs12713768 and rs10208040 [strong linkage disequilibrium (LD); D' = 1, r2 = 0.651] drives 18.9% of enhancer expression activity. Our findings highlight that the SNP rs12713768 is associated with susceptibility to and severity of hand OA in the Han Chinese population and that the suggested retinoic acid signaling pathway may play an important role in its pathogenesis.

Uterine fibroid polygenic risk score (PRS) associates and predicts risk for uterine fibroid.

Uterine fibroids (UF) are common pelvic tumors in women, heritable, and genome-wide association studies (GWAS) have identified ~ 30 loci associated with increased risk in UF. Using summary statistics from a previously published UF GWAS performed in a non-Hispanic European Ancestry (NHW) female subset from the Electronic Medical Records and Genomics (eMERGE) Network, we constructed a polygenic risk score (PRS) for UF. UF-PRS was developed using PRSice and optimized in the separate clinical population of BioVU. PRS was validated using parallel methods of 10-fold cross-validation logistic regression and phenome-wide association study (PheWAS) in a seperate subset of eMERGE NHW females (validation set), excluding samples used in GWAS. PRSice determined pt < 0.001 and after linkage disequilibrium pruning (r2 < 0.2), 4458 variants were in the PRS which was significant (pseudo-R2 = 0.0018, p = 0.041). 10-fold cross-validation logistic regression modeling of validation set revealed the model had an area under the curve (AUC) value of 0.60 (95% confidence interval [CI] 0.58-0.62) when plotted in a receiver operator curve (ROC). PheWAS identified six phecodes associated with the PRS with the most significant phenotypes being 218 'benign neoplasm of uterus' and 218.1 'uterine leiomyoma' (p = 1.94 × 10-23, OR 1.31 [95% CI 1.26-1.37] and p = 3.50 × 10-23, OR 1.32 [95% CI 1.26-1.37]). We have developed and validated the first PRS for UF. We find our PRS has predictive ability for UF and captures genetic architecture of increased risk for UF that can be used in further studies.

Human-Disease Phenotype Map Derived from PheWAS across 38,682 Individuals.

Phenome-wide association studies (PheWASs) have been a useful tool for testing associations between genetic variations and multiple complex traits or diagnoses. Linking PheWAS-based associations between phenotypes and a variant or a genomic region into a network provides a new way to investigate cross-phenotype associations, and it might broaden the understanding of genetic architecture that exists between diagnoses, genes, and pleiotropy. We created a network of associations from one of the largest PheWASs on electronic health record (EHR)-derived phenotypes across 38,682 unrelated samples from the Geisinger's biobank; the samples were genotyped through the DiscovEHR project. We computed associations between 632,574 common variants and 541 diagnosis codes. Using these associations, we constructed a "disease-disease" network (DDN) wherein pairs of diseases were connected on the basis of shared associations with a given genetic variant. The DDN provides a landscape of intra-connections within the same disease classes, as well as inter-connections across disease classes. We identified clusters of diseases with known biological connections, such as autoimmune disorders (type 1 diabetes, rheumatoid arthritis, and multiple sclerosis) and cardiovascular disorders. Previously unreported relationships between multiple diseases were identified on the basis of genetic associations as well. The network approach applied in this study can be used to uncover interactions between diseases as a result of their shared, potentially pleiotropic SNPs. Additionally, this approach might advance clinical research and even clinical practice by accelerating our understanding of disease mechanisms on the basis of similar underlying genetic associations.

Genetic Architecture of Abdominal Aortic Aneurysm in the Million Veteran Program.

BACKGROUND: Abdominal aortic aneurysm (AAA) is an important cause of cardiovascular mortality; however, its genetic determinants remain incompletely defined. In total, 10 previously identified risk loci explain a small fraction of AAA heritability. METHODS: We performed a genome-wide association study in the Million Veteran Program testing ≈18 million DNA sequence variants with AAA (7642 cases and 172 172 controls) in veterans of European ancestry with independent replication in up to 4972 cases and 99 858 controls. We then used mendelian randomization to examine the causal effects of blood pressure on AAA. We examined the association of AAA risk variants with aneurysms in the lower extremity, cerebral, and iliac arterial beds, and derived a genome-wide polygenic risk score (PRS) to identify a subset of the population at greater risk for disease. RESULTS: Through a genome-wide association study, we identified 14 novel loci, bringing the total number of known significant AAA loci to 24. In our mendelian randomization analysis, we demonstrate that a genetic increase of 10 mm Hg in diastolic blood pressure (odds ratio, 1.43 [95% CI, 1.24-1.66]; P=1.6×10-6), as opposed to systolic blood pressure (odds ratio, 1.06 [95% CI, 0.97-1.15]; P=0.2), likely has a causal relationship with AAA development. We observed that 19 of 24 AAA risk variants associate with aneurysms in at least 1 other vascular territory. A 29-variant PRS was strongly associated with AAA (odds ratioPRS, 1.26 [95% CI, 1.18-1.36]; PPRS=2.7×10-11 per SD increase in PRS), independent of family history and smoking risk factors (odds ratioPRS+family history+smoking, 1.24 [95% CI, 1.14-1.35]; PPRS=1.27×10-6). Using this PRS, we identified a subset of the population with AAA prevalence greater than that observed in screening trials informing current guidelines. CONCLUSIONS: We identify novel AAA genetic associations with therapeutic implications and identify a subset of the population at significantly increased genetic risk of AAA independent of family history. Our data suggest that extending current screening guidelines to include testing to identify those with high polygenic AAA risk, once the cost of genotyping becomes comparable with that of screening ultrasound, would significantly increase the yield of current screening at reasonable cost.

A genome-wide association study of polycystic ovary syndrome identified from electronic health records.

BACKGROUND: Polycystic ovary syndrome is the most common endocrine disorder affecting women of reproductive age. A number of criteria have been developed for clinical diagnosis of polycystic ovary syndrome, with the Rotterdam criteria being the most inclusive. Evidence suggests that polycystic ovary syndrome is significantly heritable, and previous studies have identified genetic variants associated with polycystic ovary syndrome diagnosed using different criteria. The widely adopted electronic health record system provides an opportunity to identify patients with polycystic ovary syndrome using the Rotterdam criteria for genetic studies. OBJECTIVE: To identify novel associated genetic variants under the same phenotype definition, we extracted polycystic ovary syndrome cases and unaffected controls based on the Rotterdam criteria from the electronic health records and performed a discovery-validation genome-wide association study. STUDY DESIGN: We developed a polycystic ovary syndrome phenotyping algorithm on the basis of the Rotterdam criteria and applied it to 3 electronic health record-linked biobanks to identify cases and controls for genetic study. In the discovery phase, we performed an individual genome-wide association study using the Geisinger MyCode and the Electronic Medical Records and Genomics cohorts, which were then meta-analyzed. We attempted validation of the significant association loci (P<1×10-6) in the BioVU cohort. All association analyses used logistic regression, assuming an additive genetic model, and adjusted for principal components to control for population stratification. An inverse-variance fixed-effect model was adopted for meta-analysis. In addition, we examined the top variants to evaluate their associations with each criterion in the phenotyping algorithm. We used the STRING database to characterize protein-protein interaction network. RESULTS: Using the same algorithm based on the Rotterdam criteria, we identified 2995 patients with polycystic ovary syndrome and 53,599 population controls in total (2742 cases and 51,438 controls from the discovery phase; 253 cases and 2161 controls in the validation phase). We identified 1 novel genome-wide significant variant rs17186366 (odds ratio [OR]=1.37 [1.23, 1.54], P=2.8×10-8) located near SOD2. In addition, 2 loci with suggestive association were also identified: rs113168128 (OR=1.72 [1.42, 2.10], P=5.2×10-8), an intronic variant of ERBB4 that is independent from the previously published variants, and rs144248326 (OR=2.13 [1.52, 2.86], P=8.45×10-7), a novel intronic variant in WWTR1. In the further association tests of the top 3 single-nucleotide polymorphisms with each criterion in the polycystic ovary syndrome algorithm, we found that rs17186366 (SOD2) was associated with polycystic ovaries and hyperandrogenism, whereas rs11316812 (ERBB4) and rs144248326 (WWTR1) were mainly associated with oligomenorrhea or infertility. We also validated the previously reported association with DENND1A1. Using the STRING database to characterize protein-protein interactions, we found both ERBB4 and WWTR1 can interact with YAP1, which has been previously associated with polycystic ovary syndrome. CONCLUSION: Through a discovery-validation genome-wide association study on polycystic ovary syndrome identified from electronic health records using an algorithm based on Rotterdam criteria, we identified and validated a novel genome-wide significant association with a variant near SOD2. We also identified a novel independent variant within ERBB4 and a suggestive association with WWTR1. With previously identified polycystic ovary syndrome gene YAP1, the ERBB4-YAP1-WWTR1 network suggests involvement of the epidermal growth factor receptor and the Hippo pathway in the multifactorial etiology of polycystic ovary syndrome.

LPA Variants Are Associated With Residual Cardiovascular Risk in Patients Receiving Statins.

BACKGROUND: Coronary heart disease (CHD) is a leading cause of death globally. Although therapy with statins decreases circulating levels of low-density lipoprotein cholesterol and the incidence of CHD, additional events occur despite statin therapy in some individuals. The genetic determinants of this residual cardiovascular risk remain unknown. METHODS: We performed a 2-stage genome-wide association study of CHD events during statin therapy. We first identified 3099 cases who experienced CHD events (defined as acute myocardial infarction or the need for coronary revascularization) during statin therapy and 7681 controls without CHD events during comparable intensity and duration of statin therapy from 4 sites in the Electronic Medical Records and Genomics Network. We then sought replication of candidate variants in another 160 cases and 1112 controls from a fifth Electronic Medical Records and Genomics site, which joined the network after the initial genome-wide association study. Finally, we performed a phenome-wide association study for other traits linked to the most significant locus. RESULTS: The meta-analysis identified 7 single nucleotide polymorphisms at a genome-wide level of significance within the LPA/PLG locus associated with CHD events on statin treatment. The most significant association was for an intronic single nucleotide polymorphism within LPA/PLG (rs10455872; minor allele frequency, 0.069; odds ratio, 1.58; 95% confidence interval, 1.35-1.86; P=2.6×10-10). In the replication cohort, rs10455872 was also associated with CHD events (odds ratio, 1.71; 95% confidence interval, 1.14-2.57; P=0.009). The association of this single nucleotide polymorphism with CHD events was independent of statin-induced change in low-density lipoprotein cholesterol (odds ratio, 1.62; 95% confidence interval, 1.17-2.24; P=0.004) and persisted in individuals with low-density lipoprotein cholesterol ≤70 mg/dL (odds ratio, 2.43; 95% confidence interval, 1.18-4.75; P=0.015). A phenome-wide association study supported the effect of this region on coronary heart disease and did not identify noncardiovascular phenotypes. CONCLUSIONS: Genetic variations at the LPA locus are associated with CHD events during statin therapy independently of the extent of low-density lipoprotein cholesterol lowering. This finding provides support for exploring strategies targeting circulating concentrations of lipoprotein(a) to reduce CHD events in patients receiving statins.

Dissecting genetic factors affecting phenylephrine infusion rates during anesthesia: a genome-wide association study employing EHR data.

BACKGROUND: The alpha-adrenergic agonist phenylephrine is often used to treat hypotension during anesthesia. In clinical situations, low blood pressure may require prompt intervention by intravenous bolus or infusion. Differences in responsiveness to phenylephrine treatment are commonly observed in clinical practice. Candidate gene studies indicate genetic variants may contribute to this variable response. METHODS: Pharmacological and physiological data were retrospectively extracted from routine clinical anesthetic records. Response to phenylephrine boluses could not be reliably assessed, so infusion rates were used for analysis. Unsupervised k-means clustering was conducted on clean data containing 4130 patients based on phenylephrine infusion rate and blood pressure parameters, to identify potential phenotypic subtypes. Genome-wide association studies (GWAS) were performed against average infusion rates in two cohorts: phase I (n = 1205) and phase II (n = 329). Top genetic variants identified from the meta-analysis were further examined to see if they could differentiate subgroups identified by k-means clustering. RESULTS: Three subgroups of patients with different response to phenylephrine were clustered and characterized: resistant (high infusion rate yet low mean systolic blood pressure (SBP)), intermediate (low infusion rate and low SBP), and sensitive (low infusion rate with high SBP). Differences among clusters were tabulated to assess for possible confounding influences. Comorbidity hierarchical clustering showed the resistant group had a higher prevalence of confounding factors than the intermediate and sensitive groups although overall prevalence is below 6%. Three loci with P < 1 × 10-6 were associated with phenylephrine infusion rate. Only rs11572377 with P = 6.09 × 10-7, a 3'UTR variant of EDN2, encoding a secretory vasoconstricting peptide, could significantly differentiate resistant from sensitive groups (P = 0.015 and 0.018 for phase I and phase II) or resistant from pooled sensitive and intermediate groups (P = 0.047 and 0.018). CONCLUSIONS: Retrospective analysis of electronic anesthetic records data coupled with the genetic data identified genetic variants contributing to variable sensitivity to phenylephrine infusion during anesthesia. Although the identified top gene, EDN2, has robust biological relevance to vasoconstriction by binding to endothelin type A (ETA) receptors on arterial smooth muscle cells, further functional as well as replication studies are necessary to confirm this association.

Frequent somatic transfer of mitochondrial DNA into the nuclear genome of human cancer cells.

Mitochondrial genomes are separated from the nuclear genome for most of the cell cycle by the nuclear double membrane, intervening cytoplasm, and the mitochondrial double membrane. Despite these physical barriers, we show that somatically acquired mitochondrial-nuclear genome fusion sequences are present in cancer cells. Most occur in conjunction with intranuclear genomic rearrangements, and the features of the fusion fragments indicate that nonhomologous end joining and/or replication-dependent DNA double-strand break repair are the dominant mechanisms involved. Remarkably, mitochondrial-nuclear genome fusions occur at a similar rate per base pair of DNA as interchromosomal nuclear rearrangements, indicating the presence of a high frequency of contact between mitochondrial and nuclear DNA in some somatic cells. Transmission of mitochondrial DNA to the nuclear genome occurs in neoplastically transformed cells, but we do not exclude the possibility that some mitochondrial-nuclear DNA fusions observed in cancer occurred years earlier in normal somatic cells.

Genome-wide and candidate gene approaches of clopidogrel efficacy using pharmacodynamic and clinical end points-Rationale and design of the International Clopidogrel Pharmacogenomics Consortium (ICPC).

RATIONALE: The P2Y12 receptor inhibitor clopidogrel is widely used in patients with acute coronary syndrome, percutaneous coronary intervention, or ischemic stroke. Platelet inhibition by clopidogrel shows wide interpatient variability, and high on-treatment platelet reactivity is a risk factor for atherothrombotic events, particularly in high-risk populations. CYP2C19 polymorphism plays an important role in this variability, but heritability estimates suggest that additional genetic variants remain unidentified. The aim of the International Clopidogrel Pharmacogenomics Consortium (ICPC) is to identify genetic determinants of clopidogrel pharmacodynamics and clinical response. STUDY DESIGN: Based on the data published on www.clinicaltrials.gov, clopidogrel intervention studies containing genetic and platelet function data were identified for participation. Lead investigators were invited to share DNA samples, platelet function test results, patient characteristics, and cardiovascular outcomes to perform candidate gene and genome-wide studies. RESULTS: In total, 17 study sites from 13 countries participate in the ICPC, contributing individual patient data from 8,829 patients. Available adenosine diphosphate-stimulated platelet function tests included vasodilator-stimulated phosphoprotein assay, light transmittance aggregometry, and the VerifyNow P2Y12 assay. A proof-of-principle analysis based on genotype data provided by each group showed a strong and consistent association between CYP2C19*2 and platelet reactivity (P value=5.1 × 10-40). CONCLUSION: The ICPC aims to identify new loci influencing clopidogrel efficacy by using state-of-the-art genetic approaches in a large cohort of clopidogrel-treated patients to better understand the genetic basis of on-treatment response variability.

Polymorphic analysis of CYP2C9 gene in Vietnamese population.

Genetic variations in CYP2C9 are associated to inter-individual variability of drugs metabolism and response. The only report has been done previously mainly focusing on the common variant alleles of CYP2C9 in Vietnamese Kinh subjects. However, little is known about the complete spectrum of this gene polymorphism in different ethnic groups of Vietnam. We sequenced the promoter region and all exons of CYP2C9 in 100 healthy unrelated Vietnamese Kinh subjects. Additionally, common CYP2C9 variants, *2 and *3, were also analyzed by RFLP-PCR in extra 194 Kinh subjects and 279 of other four ethnic groups in Vietnam. The results of these common variants observed from five ethnic groups were compared with other populations in the world. Seven previously reported alleles and two genotypes were determined in Kinh subjects. The percentage of CYP2C9*1 and CYP2C9*3 alleles are 96.5 and 3.5%, respectively. We found one novel non-synonymous variant in exon 7 leading to amino acid change at 363 position from proline to histidine. Functional analysis by SIFT and Polyphen-2 indicated that this mutation is intolerant and probably damaging. Prevalence of CYP2C9*2 observed in Vietnamese population was significantly lower compared with that of other populations in the South and West of Asia as well as in Europe. This study provides information of genetic distribution pattern of CYP2C9 in Vietnamese, which would be useful for optimizing drug therapies in Vietnam.

The landscape of cancer genes and mutational processes in breast cancer.

All cancers carry somatic mutations in their genomes. A subset, known as driver mutations, confer clonal selective advantage on cancer cells and are causally implicated in oncogenesis, and the remainder are passenger mutations. The driver mutations and mutational processes operative in breast cancer have not yet been comprehensively explored. Here we examine the genomes of 100 tumours for somatic copy number changes and mutations in the coding exons of protein-coding genes. The number of somatic mutations varied markedly between individual tumours. We found strong correlations between mutation number, age at which cancer was diagnosed and cancer histological grade, and observed multiple mutational signatures, including one present in about ten per cent of tumours characterized by numerous mutations of cytosine at TpC dinucleotides. Driver mutations were identified in several new cancer genes including AKT2, ARID1B, CASP8, CDKN1B, MAP3K1, MAP3K13, NCOR1, SMARCD1 and TBX3. Among the 100 tumours, we found driver mutations in at least 40 cancer genes and 73 different combinations of mutated cancer genes. The results highlight the substantial genetic diversity underlying this common disease.

Variant GADL1 and response to lithium therapy in bipolar I disorder.

BACKGROUND: Lithium has been a first-line choice for maintenance treatment of bipolar disorders to prevent relapse of mania and depression, but many patients do not have a response to lithium treatment. METHODS: We selected subgroups from a sample of 1761 patients of Han Chinese descent with bipolar I disorder who were recruited by the Taiwan Bipolar Consortium. We assessed their response to lithium treatment using the Alda scale and performed a genomewide association study on samples from one subgroup of 294 patients with bipolar I disorder who were receiving lithium treatment. We then tested the single-nucleotide polymorphisms (SNPs) that showed the strongest association with a response to lithium for association in a replication sample of 100 patients and tested them further in a follow-up sample of 24 patients. We sequenced the exons, exon-intron boundaries, and part of the promoter of the gene encoding glutamate decarboxylase-like protein 1 (GADL1) in 94 patients who had a response to lithium and in 94 patients who did not have a response in the genomewide association sample. RESULTS: Two SNPs in high linkage disequilibrium, rs17026688 and rs17026651, that are located in the introns of GADL1 showed the strongest associations in the genomewide association study (P=5.50×10(-37) and P=2.52×10(-37), respectively) and in the replication sample of 100 patients (P=9.19×10(-15) for each SNP). These two SNPs had a sensitivity of 93% for predicting a response to lithium and differentiated between patients with a good response and those with a poor response in the follow-up cohort. Resequencing of GADL1 revealed a novel variant, IVS8+48delG, which lies in intron 8 of the gene, is in complete linkage disequilibrium with rs17026688 and is predicted to affect splicing. CONCLUSIONS: Genetic variations in GADL1 are associated with the response to lithium maintenance treatment for bipolar I disorder in patients of Han Chinese descent. (Funded by Academia Sinica and others.).

Genome-wide association and replication study of anti-tuberculosis drugs-induced liver toxicity.

BACKGROUND: Drug-induced liver injury (DILI) is a well-recognized adverse event of anti tuberculosis drugs (ATD) possibly associated with genetic variations. The objective of this study was to perform genome-wide association study (GWAS) to identify genetic variants associated with the risk for ATD induced liver toxicity in Ethiopian patients. RESULT: Treatment-naïve newly diagnosed tuberculosis patients (n = 646) were enrolled prospectively and treated with rifampicin based short course anti-tuberculosis therapy. Whole genome genotyping was done using Illumina Omni Express Exome Bead Chip genotyping array with 951,117 single nucleotide polymorphisms (SNPs) on 48 DILI cases and 354 ATD tolerants. Replication study was carried out for 50 SNPs with the lowest P-values (top SNPs) using an independent cohort consisting of 27 DILI cases and 217 ATD tolerants. In the combined analysis, the top SNP identified was rs10946737 (P = 4.4 × 10-6, OR = 3.4, 95 % confidence interval = 2.2-5.3) in the intron of FAM65B in chromosome 6. In addition, we identified a cluster of SNPs with suggestive genome-wide significance in the intron of ATP/GTP binding protein-like 4 (AGBL4). CONCLUSION: We identified genetic variants that are potentially associated with ATD induced liver toxicity. Further studies with larger sample sizes are essential to confirm the findings.

Identification of susceptibility gene associated with female primary Sjögren's syndrome in Han Chinese by genome-wide association study.

Primary Sjögren's syndrome (PSS) is an autoimmune disease targeting exocrine glands. It ten times more dominantly affects women than men with an onset peak at menopause. The genetic factor predisposing women to PSS remains unclear. Therefore, we aimed to identify susceptibility loci for PSS in women. We performed genome-wide association study (GWAS) using 242 female PSS patients and 1444 female control in Han Chinese population residing in Taiwan. Replication was conducted in an independent cohort of 178 female PSS and 14,432 control subjects. We identified rs117026326 on GTF2I with GWAS significance (P = 1.10 × 10-15) and rs13079920 on RBMS3 with suggestive significance (P = 2.90 × 10-5) associating with PSS in women. The association of RBMS3 was further evidenced by imputation in which rs13072846 (P = 4.89 × 10-5) was identified and confirmed as female PSS associating SNP within the same LD with rs13079920. PSS pathogenesis involves both immune (effector) and exocrine (target) system. We suggested that while GTF2I is a previously reported associating gene which may function in immune system, RBMS3 is a novel susceptibility gene that predisposes women to PSS potentially through modulating acinar apoptosis and TGF-ß signaling in target exocrine system.

Success stories in genomic medicine from resource-limited countries.

In recent years, the translation of genomic discoveries into mainstream medical practice and public health has gained momentum, facilitated by the advent of new technologies. However, there are often major discrepancies in the pace of implementation of genomic medicine between developed and developing/resource-limited countries. The main reason does not only lie in the limitation of resources but also in the slow pace of adoption of the new findings and the poor understanding of the potential that this new discipline offers to rationalize medical diagnosis and treatment. Here, we present and critically discuss examples from the successful implementation of genomic medicine in resource-limited countries, focusing on pharmacogenomics, genome informatics, and public health genomics, emphasizing in the latter case genomic education, stakeholder analysis, and economics in pharmacogenomics. These examples can be considered as model cases and be readily replicated for the wide implementation of pharmacogenomics and genomic medicine in other resource-limited environments.

Pharmacogenomics of adverse drug reactions: implementing personalized medicine.

Pharmacogenomics aims to investigate the genetic basis of inter-individual differences in drug responses, such as efficacy, dose requirements and adverse events. Research in pharmacogenomics has grown over the past decade, evolving from a candidate-gene approach to genome-wide association studies (GWASs). Genetic variants in genes coding for drug metabolism, drug transport and more recently human-leukocyte antigens (HLAs) have been linked to inter-individual differences in the risk of adverse drug reactions (ADRs). The tight association of specific HLA alleles with Stevens-Johnson syndrome, toxic epidermal necrolysis, drug hypersensitivity syndrome and drug-induced liver injury underscore the importance of HLA in the pathogenesis of these idiosyncratic drug hypersensitivity reactions. However, as with the search for the genetic basis for common diseases, pharmacogenomic research, including GWAS, has so far been a disappointment in discovering major gene variants responsible for the efficacy of drugs used to treat common diseases. This review focuses on the pharmacogenomics of ADRs, the underlying mechanisms and the potential use of genomic biomarkers in clinical practice for dose adjustment and the avoidance of drug toxicity. We also discuss obstacles to the implementation of pharmacogenomics and the direction of future translational research.