Search

Table S1 provides a summary of the genetic association results in the European only cohort

Figure S4 provides examples of audiogram configurations

Purpose: Adverse drug reactions such as ototoxicity, which occurs in approximately one-fifth of adult patients who receive cisplatin treatment, can incur large socioeconomic burdens on patients with testicular cancer who develop this cancer during early adulthood. Recent genome-wide association studies have identified genetic variants in ACYP2 and WFS1 that are associated with cisplatin-induced ototoxicity. We sought to explore the role of these genetic susceptibility factors to cisplatin-induced ototoxicity in patients with testicular cancer.Experimental Design: Extensive clinical and demographic data were collected for 229 patients with testicular cancer treated with cisplatin. Patients were genotyped for two variants, ACYP2 rs1872328 and WFS1 rs62283056, that have previously been associated with hearing loss in cisplatin-treated patients. Analyses were performed to investigate the association of these variants with ototoxicity in this cohort of adult patients with testicular cancer.Results: Pharmacogenomic analyses revealed that ACYP2 rs1872328 was significantly associated with cisplatin-induced ototoxicity [P = 2.83 × 10−3, OR (95% CI):14.7 (2.6–84.2)]. WFS1 rs62283056 was not significantly associated with ototoxicity caused by cisplatin (P = 0.39); however, this variant was associated with hearing loss attributable to any cause [P = 5.67 × 10−3, OR (95% CI): 3.2 (1.4–7.7)].Conclusions: This study has provided the first evidence for the role of ACYP2 rs1872328 in cisplatin-induced ototoxicity in patients with testicular cancer. These results support the use of this information to guide the development of strategies to prevent cisplatin-induced ototoxicity across cancers. Further, this study has highlighted the importance of phenotypic differences in replication studies and has provided further evidence for the role of WFS1 rs62283056 in susceptibility to hearing loss, which may be worsened by cisplatin treatment. Clin Cancer Res; 24(8); 1866–71. ©2018 AACR.

Huntington disease (HD) is a neurodegenerative disorder that is caused by a CAG repeat expansion in HTT. The length of this repeat, however, only explains a proportion of the variability in age of onset in patients. Genome-wide association studies have identified modifiers that contribute toward a proportion of the observed variance. By incorporating tissue-specific transcriptomic information with these results, additional modifiers can be identified. We performed a transcriptome-wide association study assessing heritable differences in genetically determined expression in diverse tissues, with genome-wide data from over 4000 patients. Functional validation of prioritized genes was undertaken in isogenic HD stem cells and patient brains. Enrichment analyses were performed with biologically relevant gene sets to identify the core pathways. HD-associated gene coexpression modules were assessed for associations with neurological phenotypes in an independent cohort and to guide drug repurposing analyses. Transcriptomic analyses identified genes that were associated with age of HD onset and displayed colocalization with gene expression signals in brain tissue (FAN1, GPR161, PMS2, SUMF2), with supporting evidence from functional experiments. This included genes involved in DNA repair, as well as novel-candidate modifier genes that have been associated with other neurological conditions. Further, cortical coexpression modules were also associated with cognitive decline and HD-related traits in a longitudinal cohort. In summary, the combination of population-scale gene expression information with HD patient genomic data identified novel modifier genes for the disorder. Further, these analyses expanded the pathways potentially involved in modifying HD onset and prioritized candidate therapeutics for future study., Graphical Abstract Graphical Abstract Schematic representation of methodology employed to identify transcriptomic modifiers of Huntington disease (HD). Transcriptome-wide association analyses facilitated (A) the identification of top transcriptomic HD onset modifier genes, (B) relevant biologically diverse gene sets that are enriched for this trait and (C) human cortical coexpression modules that are relevant for HD onset and related phenotypes in an aging cohort, along with signature matching of HD signals for drug repurposing.

Purpose In some Huntington disease (HD) patients, the “loss of interruption” (LOI) variant eliminates an interrupting codon in the HTT CAG-repeat tract, which causes earlier age of onset (AOO). The magnitude of this effect is uncertain, since previous studies included few LOI carriers, and the variant also causes CAG size misestimation. We developed a rapid LOI detection screen, enabling unbiased frequency estimation among manifest HD patients. Additionally, we combined published data with clinical data from newly identified patients to accurately characterize the LOI’s effect on AOO. Methods We developed a LOI detection polymerase chain reaction (PCR) assay, and screened patients to estimate the frequency of the LOI variant and its effect on AOO. Results Mean onset for LOI carriers (n = 49) is 20.4 years earlier than expected based on diagnosed CAG size. After correcting for CAG size underestimation, the variant is still associated with onset 9.5 years earlier. The LOI is present in 1.02% of symptomatic HD patients, and in 32.2% of symptomatic reduced penetrance (RP) range patients (36–39 CAGs). Conclusion The LOI causes significantly earlier onset, greater than expected by CAG length, particularly in persons with 36–39 CAG repeats. Detection of this variant has implications for HD families, especially for those in the RP range.

Background: Ototoxicity (hearing loss, tinnitus and/or vertigo) is a serious adverse event of cisplatin treatment in children with cancer. The heterogeneity in ototoxicity occurrence after similar treatment suggests a role for genetic susceptibility. This study investigated the association between carriership of novel single nucleotide polymorphisms (SNPs) and cisplatin-induced hearing loss (CIHL) in childhood cancer patients.Material and methods: The discovery cohort included cisplatin treated, non-cranial irradiated pediatric cancer patients within the European PanCareLIFE (PCL) study (N=390). CIHL at end of cancer treatment was defined as Muenster grade >=2b, assessed by pure tone audiometry. DNA was genotyped using the Infinium© Global Screening Array. Logistic regression models were applied including age at diagnosis, sex, cisplatin total cumulative dose and principal components 1-4, assuming an additive effect of the minor allele. Replication of the findings was performed in two independent, similarly treated cohorts (N=192 and N=188). Functional validation experiments in cultured human HeLa cell lines were performed to determine the effect of knockdown of the SNPs nearest identified gene on cisplatin-induced toxicity.Results: In the PCL discovery cohort, 8 SNPs reached a suggestive significance of P

Polygenic scores (PGSs) have emerged as promising tools for complex trait risk prediction. The application of these scores to pharmacogenomics provides new opportunities to improve the prediction of treatment outcomes. To gain insight into this area of research, we conducted a systematic review and accompanying analysis. This review uncovered 51 papers examining the use of PGSs for drug-related outcomes, with the majority of these papers focusing on the treatment of psychiatric disorders (n = 30). Due to difficulties in collecting large cohorts of uniformly treated patients, the majority of pharmacogenomic PGSs were derived from large-scale genome-wide association studies of disease phenotypes that were related to the pharmacogenomic phenotypes under investigation (e.g., schizophrenia-derived PGSs for antipsychotic response prediction). Examination of the research participants included in these studies revealed that the majority of cohort participants were of European descent (78.4%). These biases were also reflected in research affiliations, which were heavily weighted towards institutions located in Europe and North America, with no first or last authors originating from institutions in Africa or South Asia. There was also substantial variability in the methods used to develop PGSs, with between 3 and 6.6 million variants included in the PGSs. Finally, we observed significant inconsistencies in the reporting of PGS analyses and results, particularly in terms of risk model development and application, coupled with a lack of data transparency and availability, with only three pharmacogenomics PGSs deposited on the Polygenic Score Catalog. These findings highlight current gaps and key areas for future pharmacogenomic PGS research.

We report an inborn error of metabolism caused by an expansion of a GCA-repeat tract in the 5′ untranslated region of the gene encoding glutaminase (GLS) that was identified through detailed clinical and biochemical phenotyping, combined with whole-genome sequencing. The expansion was observed in three unrelated patients who presented with an early-onset delay in overall development, progressive ataxia, and elevated levels of glutamine. In addition to ataxia, one patient also showed cerebellar atrophy. The expansion was associated with a relative deficiency of GLS messenger RNA transcribed from the expanded allele, which probably resulted from repeat-mediated chromatin changes upstream of the GLS repeat. Our discovery underscores the importance of careful examination of regions of the genome that are typically excluded from or poorly captured by exome sequencing.

Aim: To improve the identification and interpretation of pharmacogenetic variants through the integration of disease and drug-related traits. Materials & methods: We hypothesized that integrating genome-wide disease and pharmacogenomic data may drive new insights into drug toxicity and response by identifying shared genetic architecture. Pleiotropic variants were identified using a methodological framework incorporating colocalization analysis. Results: Using genome-wide association studies summary statistics from the UK Biobank, European Bioinformatics Institute genome-wide association studies catalog and the Pharmacogenomics Research Network, we validated pleiotropy at the ABCG2 locus between allopurinol response and gout and identified novel pleiotropy between antihypertensive-induced new-onset diabetes, Crohn’s disease and inflammatory bowel disease at the IL18RAP/SLC9A4 locus. Conclusion: New mechanistic insights and genetic loci can be uncovered by identifying pleiotropy between disease and drug-related traits.

The authors regret not all contributing authors correctly acknowledged funding. Jonathan Peter’s IMARI-Africa project (AFRISCAR) is part of the EDCTP2 programme supported by the European Union (grant number TMA2017SF-1981 ). The authors would like to apologise for any inconvenience caused.

Summary Background Huntington's disease is a fatal neurodegenerative disorder that is caused by CAG-CAA repeat expansion, encoding polyglutamine, in the huntingtin (HTT) gene. Current age-of-clinical-onset prediction models for Huntington's disease are based on polyglutamine length and explain only a proportion of the variability in age of onset observed between patients. These length-based assays do not interrogate the underlying genetic variation, because known genetic variants in this region do not alter the protein coding sequence. Given that individuals with identical repeat lengths can present with Huntington's disease decades apart, the search for genetic modifiers of clinical age of onset has become an active area of research. Recent developments Results from three independent genetic studies of Huntington's disease have shown that glutamine-encoding CAA variants that interrupt DNA CAG repeat tracts, but do not alter polyglutamine length or polyglutamine homogeneity, are associated with substantial differences in age of onset of Huntington's disease in carriers. A variant that results in the loss of CAA interruption is associated with early onset and is particularly relevant to individuals that carry alleles in the reduced penetrance range (ie, CAG 36–39). Approximately a third of clinically manifesting carriers of reduced penetrance alleles, defined by current diagnostics, carry this variant. Somatic repeat instability, modified by interrupted CAG tracts, is the most probable cause mediating this effect. This relationship is supported by genome-wide screens for disease modifiers, which have revealed the importance of DNA-repair genes in Huntington's disease (ie, FAN1, LIG1, MLH1, MSH3, PMS1, and PMS2). Where next? Focus needs to be placed on refining our understanding of the effect of the loss-of-interruption and duplication-of-interruption variants and other interrupting sequence variants on age of onset, and assessing their effect in disease-relevant brain tissues, as well as in diverse population groups, such as individuals from Africa and Asia. Diagnostic tests should be augmented or updated, since current tests do not assess the underlying DNA sequence variation, especially when assessing individuals that carry alleles in the reduced penetrance range. Future studies should explore somatic repeat instability and DNA repair as new therapeutic targets to modify age of onset in Huntington's disease and in other repeat-mediated disorders. Disease-modifying therapies could potentially be developed by therapeutically targeting these processes. Promising approaches include therapeutically targeting the expanded repeat or directly perturbing key DNA-repair genes (eg, with antisense oligonucleotides or small molecules). Targeting the CAG repeat directly with naphthyridine-azaquinolone, a compound that induces contractions, and altering the expression of MSH3, represent two viable therapeutic strategies. However, as a first step, the capability of such novel therapeutic approaches to delay clinical onset in animal models should be assessed.

Adverse drug reactions (ADRs) are a major problem in modern medicine, representing up to the fourth-highest cause of mortality. Pharmacogenomic tests are 1 of the most promising methods to tackle the challenge of ADRs. The objective of this study was to analyze the clinical and demographic information of the pan-Canadian active surveillance network, Canadian Pharmacogenomics Network for Drug Safety (CPNDS). Information entered into the database by trained active surveillors between May 15, 2005 and May 9, 2017 was collected and analyzed. Specific data included for analysis were number of ADR reports, reports of drug use without ADRs, date of onset of ADR, suspected drugs, concomitant drugs, and fatal ADR cases. The CPNDS database consisted of 93,974 reports of medication use, including 10,475 reports of ADRs, of which 72.6% occurred in pediatric patients (≤21 years old). Self-reported ancestries were predominantly Europe (38.2%), Canada (9.6%), and East Asia (4.9%). The 5 most frequent ADRs were cutaneous ADRs, peripheral neuropathy, cardiotoxicity, central nervous system toxicity, and ototoxicity. The 5 drugs most commonly suspected to cause ADRs were methotrexate, vincristine, doxorubicin, cisplatin, and L-asparaginase. The CPNDS database is a valuable resource to identify clinical and genomic predictors of ADRs. The database also highlights our candidate ADRs for pharmacogenomic discovery research to identify additional ADR biomarkers. Additionally, the database provides information that can be used for developing strategies to prevent ADRs and raises awareness of ADRs among Canadian healthcare professionals.

Purpose: Adverse drug reactions such as ototoxicity, which occurs in approximately one-fifth of adult patients who receive cisplatin treatment, can incur large socioeconomic burdens on patients with testicular cancer who develop this cancer during early adulthood. Recent genome-wide association studies have identified genetic variants in ACYP2 and WFS1 that are associated with cisplatin-induced ototoxicity. We sought to explore the role of these genetic susceptibility factors to cisplatin-induced ototoxicity in patients with testicular cancer. Experimental Design: Extensive clinical and demographic data were collected for 229 patients with testicular cancer treated with cisplatin. Patients were genotyped for two variants, ACYP2 rs1872328 and WFS1 rs62283056, that have previously been associated with hearing loss in cisplatin-treated patients. Analyses were performed to investigate the association of these variants with ototoxicity in this cohort of adult patients with testicular cancer. Results: Pharmacogenomic analyses revealed that ACYP2 rs1872328 was significantly associated with cisplatin-induced ototoxicity [P = 2.83 × 10−3, OR (95% CI):14.7 (2.6–84.2)]. WFS1 rs62283056 was not significantly associated with ototoxicity caused by cisplatin (P = 0.39); however, this variant was associated with hearing loss attributable to any cause [P = 5.67 × 10−3, OR (95% CI): 3.2 (1.4–7.7)]. Conclusions: This study has provided the first evidence for the role of ACYP2 rs1872328 in cisplatin-induced ototoxicity in patients with testicular cancer. These results support the use of this information to guide the development of strategies to prevent cisplatin-induced ototoxicity across cancers. Further, this study has highlighted the importance of phenotypic differences in replication studies and has provided further evidence for the role of WFS1 rs62283056 in susceptibility to hearing loss, which may be worsened by cisplatin treatment. Clin Cancer Res; 24(8); 1866–71. ©2018 AACR.

Huntington disease (HD) is the most common monogenic neurodegenerative disorder in populations of European ancestry, but occurs at lower prevalence in populations of East Asian or black African descent. New mutations for HD result from CAG repeat expansions of intermediate alleles (IAs), usually of paternal origin. The differing prevalence of HD may be related to the rate of new mutations in a population, but no comparative estimates of IA frequency or the HD new mutation rate are available. In this study, we characterize IA frequency and the CAG repeat distribution in fifteen populations of diverse ethnic origin. We estimate the HD new mutation rate in a series of populations using molecular IA expansion rates. The frequency of IAs was highest in Hispanic Americans and Northern Europeans, and lowest in black Africans and East Asians. The prevalence of HD correlated with the frequency of IAs by population and with the proportion of IAs found on the HD-associated A1 haplotype. The HD new mutation rate was estimated to be highest in populations with the highest frequency of IAs. In European ancestry populations, one in 5,372 individuals from the general population and 7.1% of individuals with an expanded CAG repeat in the HD range are estimated to have a molecular new mutation. Our data suggest that the new mutation rate for HD varies substantially between populations, and that IA frequency and haplotype are closely linked to observed epidemiological differences in the prevalence of HD across major ancestry groups in different countries.

Huntington disease (HD) is a dominant neurodegenerative disorder caused by a CAG repeat expansion in the Huntingtin (HTT) gene. HD occurs worldwide, but the causative mutation is found on different HTT haplotypes in distinct ethnic groups. In Latin America, HD is thought to have European origins, but indigenous Amerindian ancestry has not been investigated. Here, we report dense HTT haplotypes in 62 mestizo Peruvian HD families, 17 HD families from across Latin America, and 42 controls of defined Peruvian Amerindian ethnicity to determine the origin of HD in populations of admixed Amerindian and European descent. HD in Peru occurs most frequently on the A1 HTT haplotype (73%), as in Europe, but on an unexpected indigenous variant also found in Amerindian controls. This Amerindian A1 HTT haplotype predominates over the European A1 variant among geographically disparate Latin American controls and in HD families from across Latin America, supporting an indigenous origin of the HD mutation in mestizo American populations. We also show that a proportion of HD mutations in Peru occur on a C1 HTT haplotype of putative Amerindian origin (14%). The majority of HD mutations in Latin America may therefore occur on haplotypes of Amerindian ancestry rather than on haplotypes resulting from European admixture. Despite the distinct ethnic ancestry of Amerindian and European A1 HTT, alleles on the parent A1 HTT haplotype allow for development of identical antisense molecules to selectively silence the HD mutation in the greatest proportion of patients in both Latin American and European populations.

Stevens-Johnson syndrome and toxic epidermal necrolysis (SJS/TEN) are potentially life-threatening, immune-mediated adverse reactions characterized by widespread erythema, epidermal necrosis, and detachment of skin and mucosa. Efforts to grow and develop functional international collaborations and a multidisciplinary interactive network focusing on SJS/TEN as an uncommon but high burden disease will be necessary to improve efforts in prevention, early diagnosis and improved acute and long-term management. SJS/TEN 2019: From Science to Translation was a 1.5-day scientific program held April 26-27, 2019, in Vancouver, Canada. The meeting successfully engaged clinicians, researchers, and patients and conducted many productive discussions on research and patient care needs.

Anthracyclines are highly effective chemotherapeutic agents; however, their clinical utility is limited by severe anthracycline-induced cardiotoxicity (ACT). Genome-wide association studies (GWAS) have uncovered several genetic variants associated with ACT, but the impact of these findings requires further elucidation. We conducted a transcriptome-wide association study (TWAS) using our previous GWAS summary statistics (n = 280 patients) to identify gene expression-related associations with ACT. We identified a genetic association between decreased expression of GDF5 and ACT (Z-score = −4.30, P = 1.70 × 10−5), which was replicated in an independent cohort (n = 845 patients, P = 3.54 × 10−3). Additionally, cell viability of GDF5-silenced human cardiac myocytes was significantly decreased in response to anthracycline treatment. Subsequent gene set enrichment and pathway analyses of the TWAS data revealed that genes essential for survival, cardioprotection and response to anthracyclines, as well as genes involved in ribosomal, spliceosomal and cardiomyopathy pathways are important for the development of ACT.

Huntington disease (HD) is a neurodegenerative disorder that is caused by a CAG repeat expansion in the HTT gene. In an attempt to identify genomic modifiers that contribute towards the age of onset of HD, we performed a transcriptome wide association study assessing heritable differences in genetically determined expression in diverse tissues, employing genome wide data from over 4,000 patients. This identified genes that showed evidence for colocalization and replication, with downstream functional validation being performed in isogenic HD stem cells and patient brains. Enrichment analyses detected associations with various biologically-relevant gene sets and striatal coexpression modules that are mediated by CAG length. Further, cortical coexpression modules that are relevant for HD onset were also associated with cognitive decline and HD-related traits in a longitudinal cohort. In summary, the combination of population-scale gene expression information with HD patient genomic data identified novel modifier genes for the disorder.

Huntington disease (HD) is caused by a CAG repeat expansion in the huntingtin (HTT) gene. Although the length of this repeat is inversely correlated with age of onset (AOO), it does not fully explain the variability in AOO. We assessed the sequence downstream of the CAG repeat in HTT [reference: (CAG)n-CAA-CAG], since variants within this region have been previously described, but no study of AOO has been performed. These analyses identified a variant that results in complete loss of interrupting (LOI) adenine nucleotides in this region [(CAG)n-CAG-CAG]. Analysis of multiple HD pedigrees showed that this LOI variant is associated with dramatically earlier AOO (average of 25 years) despite the same polyglutamine length as in individuals with the interrupting penultimate CAA codon. This LOI allele is particularly frequent in persons with reduced penetrance alleles who manifest with HD and increases the likelihood of presenting clinically with HD with a CAG of 36-39 repeats. Further, we show that the LOI variant is associated with increased somatic repeat instability, highlighting this as a significant driver of this effect. These findings indicate that the number of uninterrupted CAG repeats, which is lengthened by the LOI, is the most significant contributor to AOO of HD and is more significant than polyglutamine length, which is not altered in these individuals. In addition, we identified another variant in this region, where the CAA-CAG sequence is duplicated, which was associated with later AOO. Identification of these cis-acting modifiers have potentially important implications for genetic counselling in HD-affected families.

Although previous research identified candidate genetic polymorphisms associated with cisplatin nephrotoxicity, varying outcome definitions potentially contributed to the variability in the effect size and direction of this relationship. We selected genetic variants that have been significantly associated with cisplatin-induced nephrotoxicity in more than one published study (SLC22A2 rs316019, ERCC1 rs11615 and rs3212986, ERCC2 rs1799793 and rs13181) and performed a replication analysis to confirm associations between these genetic polymorphisms and cisplatin nephrotoxicity using various outcome definitions. We included 282 germ cell testicular cancer patients treated with cisplatin from 2009&ndash, 2014, aged &gt, 17 years recruited by the Canadian Pharmacogenomics Network for Drug Safety. Nephrotoxicity was defined using four grading tools: (1) Common Terminology Criteria for Adverse Events (CTCAE) v4.03 for acute kidney injury (AKI) or CTCAE-AKI, (2) adjusted cisplatin-induced AKI, (3) elevation of serum creatinine, and (4) reduction in the estimated glomerular filtration rate (eGFR). Significant associations were only found when using the CTCAE v4.03 definition: genotype CA of the ERCC1 rs3212986 was associated with decreased risk of cisplatin nephrotoxicity (ORadj = 0.24, 95% CI:0.08&ndash, 0.70, p = 0.009) compared to genotype CC. In contrast, addition of allele A at SLC22A2 rs316019 was associated with increased risk (ORadj = 4.41, 95% CI:1.96&ndash, 9.88, p &lt, 0.001) while genotype AC was associated with a higher risk of cisplatin nephrotoxicity (ORadj = 5.06, 95% CI:1.69&ndash, 15.16, p = 0.004) compared to genotype CC. Our study showed that different case definitions led to variability in the genetic risk ascertainment of cisplatin nephrotoxicity. Therefore, consensus on a set of clinically relevant outcome definitions that all such studies should follow is needed.

BACKGROUND & AIMS: We performed genetic analyses of a multiethnic cohort of patients with idiosyncratic drug-induced liver injury (DILI) to identify variants associated with susceptibility. METHODS: We performed a genome-wide association study of 2048 individuals with DILI (cases) and 12,429 individuals without (controls). Our analysis included subjects of European (1806 cases and 10,397 controls), African American (133 cases and 1,314 controls), and Hispanic (109 cases and 718 controls) ancestry. We analyzed DNA from 113 Icelandic cases and 239,304 controls to validate our findings. RESULTS: We associated idiosyncratic DILI with rs2476601, a nonsynonymous polymorphism that encodes a substitution of tryptophan with arginine in the protein tyrosine phosphatase, non-receptor type 22 gene (PTPN22) (odds ratio [OR], 1.44; 95% CI, 1.28–1.62; P=1.2×10(−9) and replicated the finding in the validation set (OR, 1.48; 95% CI, 1.09–1.99; P=.01). The minor allele frequency showed the same effect size (OR > 1) among ethnic groups. The strongest association was with amoxicillin and clavulanate-associated DILI in persons of European ancestry (OR, 1.62; 95% CI, 1.32–1.98; P=4.0×10(−6); allele frequency=13.3%), but the polymorphism was associated with DILI of other causes (OR, 1.37; 95% CI, 1.21–1.56; P= 1.5×10(−6); allele frequency=11.5%). Among amoxicillin- and clavulanate-associated cases of European ancestry, rs2476601 doubled the risk for DILI among those with the HLA risk alleles A*02:01 and DRB1*15:01. CONCLUSIONS: In a genome-wide association study, we identified rs2476601 in PTPN22 as a non-HLA variant that associates with risk of liver injury caused by multiple drugs and validated our finding in a separate cohort. This variant has been associated with increased risk of autoimmune diseases, providing support for the concept that alterations in immune regulation contribute to idiosyncratic DILI.

Huntington disease (HD) is an autosomal dominant neurological disorder that is caused by a CAG repeat expansion, translated into polyglutamine, in the huntingtin (HTT) gene. Although the length of this repeat polymorphism is inversely correlated with age of onset (AOO), it does not fully explain the variability in AOO. Genomic studies have provided evidence for the involvement of DNA repair in modifying this trait, potentially through somatic repeat instability. We therefore assessed genetic variants within the 12bp interrupting sequence between the pathogenic CAG repeat and the polymorphic proline (CCG) tract in the HTT gene and identified variants that result in complete loss of interruption (LOI) between the adjacent CAG/CCG repeats. Analysis of multiple HD pedigrees showed that this variant is associated with dramatically earlier AOO and is particularly relevant to HD patients with reduced penetrance alleles. On average AOO of HD is hastened by an average of 25 years in LOI carriers. This finding indicates that the number of uninterrupted CAG repeats is the most significant contributor to AOO of HD and is more impactful than polyglutamine length, which is not altered in these patients. We show that the LOI variant is associated with increases in both somatic and germline repeat instability, demonstrating a potential mechanism for this effect. Screening individuals from the general population (n=2,674 alleles) suggests that the variant occurs only in expanded CAG repeat alleles. Identification of this modifier has important clinical implications for disease management of HD families, especially for those in the reduced penetrance ranges.

Long-term exposure to efavirenz (EFV) measured in hair samples may predict response to antiretroviral treatment (ART). Polymorphisms in CYP2B6 are known to alter EFV levels. The aim of this study was to assess the relationship between CYP2B6 genotype, EFV levels measured in hair, and virological outcomes on ART in a real-world setting. We measured EFV levels in hair from HIV-positive South African females who had been receiving EFV-based treatment for at least 3 months from the South African Black (SAB) (n = 81) and Cape Mixed Ancestry (CMA) (n = 53) populations. Common genetic variation in CYP2B6 was determined in 15 individuals from each population using bidirectional Sanger sequencing. Prioritized variants (n = 16) were subsequently genotyped in the entire patient cohort (n = 134). The predictive value of EFV levels in hair and selected variants in CYP2B6 on virological treatment outcomes was assessed. Previously described alleles (CYP2B6*2, CYP2B6*5, CYP2B6*6, CYP2B6*17, and CYP2B6*18), as well as two novel alleles (CYP2B6*31 and CYP2B6*32), were detected in this study. Compared to noncarriers, individuals homozygous for CYP2B6*6 had ∼109% increased EFV levels in hair (p = .016) and CYP2B6*18 heterozygotes demonstrated 82% higher EFV hair levels (p = .0006). This study confirmed that alleles affecting CYP2B6 metabolism and subsequent EFV exposure are present at significant frequencies in both the SAB and CMA populations. Furthermore, this study demonstrated that the use of hair samples for testing EFV concentrations may be a useful tool in determining long-term drug exposure in resource-limited countries.

Importance Cisplatin-induced ototoxic effects are an important complication that affects testicular cancer survivors as a consequence of treatment. The identification of genetic variants associated with this adverse drug reaction will further our mechanistic understanding of its development and potentially lead to strategies to prevent ototoxic effects. Objective To identify the genetic variants associated with cisplatin-induced ototoxic effects in adult testicular cancer patients. Design, Setting, and Participants This retrospective study was performed by the Canadian Pharmacogenomics Network for Drug Safety using patients recruited from 5 adult oncology treatment centers across Canada. Male patients who were 17 years or older, diagnosed with germ cell testicular cancer, and previously treated with cisplatin-based chemotherapy were recruited from July 2009 to April 2013 using active surveillance methodology. Cisplatin-induced ototoxic effects were independently diagnosed by 2 audiologists. Patients were genotyped for 7907 variants using a custom pharmacogenomic array. Logistic regression was used to identify genetic variants that were significantly associated with ototoxic effects. The validity of these findings was confirmed through independent replication and cell-based functional assays. Exposures Cisplatin-based chemotherapy. Main Outcomes and Measures Cisplatin-induced ototoxic effects. Results After exclusions, 188 patients (median [interquartile range] age, 31 [24-39] years) were enrolled in this study to form the discovery and replication cohorts. Association and fine-mapping analyses identified a protein-coding variant,rs4788863in SLC16A5 , that was associated with protection against cisplatin-induced ototoxic effects in 2 independent cohorts (combined cohort: odds ratio, 0.06; 95% CI, 0.02-0.22; P = 2.17 × 10 −7 ). Functional validation of this transporter gene revealed that in vitro SLC16A5 –silencing altered cellular responses to cisplatin treatment, supporting a role for SLC16A5 in the development of cisplatin-induced ototoxic effects. These results were further supported by the literature, which provided confirmatory evidence for the role that SLC16A5 plays in hearing. Conclusions and Relevance This study has identified a novel association between protein-coding variation in SLC16A5 and cisplatin-induced ototoxic effects. These findings have provided insight into the molecular mechanisms of this adverse drug reaction in adult patients with germ cell testicular cancer. Given that previous studies have shown that cimetidine, an SLC16A5-inhibitor, prevents murine cisplatin-induced ototoxic effects, the findings from this study have important implications for otoprotectant strategies in humans.

Mutations in the GJB2 gene, encoding connexin 26, could account for 50% of congenital, nonsyndromic, recessive deafness cases in some Caucasian/Asian populations. There is a scarcity of published data in sub-Saharan Africans. We Sanger sequenced the coding region of the GJB2 gene in 205 Cameroonian and Xhosa South Africans with congenital, nonsyndromic deafness; and performed bioinformatic analysis of variations in the GJB2 gene, incorporating data from the 1000 Genomes Project. Amongst Cameroonian patients, 26.1% were familial. The majority of patients (70%) suffered from sensorineural hearing loss. Ten GJB2 genetic variants were detected by sequencing. A previously reported pathogenic mutation, g.3741_3743delTTC (p.F142del), and a putative pathogenic mutation, g.3816G>A (p.V167M), were identified in single heterozygous samples. Amongst eight the remaining variants, two novel variants, g.3318-41G>A and g.3332G>A, were reported. There were no statistically significant differences in allele frequencies between cases and controls. Principal Components Analyses differentiated between Africans, Asians, and Europeans, but only explained 40% of the variation. The present study is the first to compare African GJB2 sequences with the data from the 1000 Genomes Project and have revealed the low variation between population groups. This finding has emphasized the hypothesis that the prevalence of mutations in GJB2 in nonsyndromic deafness amongst European and Asian populations is due to founder effects arising after these individuals migrated out of Africa, and not to a putative “protective” variant in the genomic structure of GJB2 in Africans. Our results confirm that mutations in GJB2 are not associated with nonsyndromic deafness in Africans.

Multiple sclerosis (MS) is a disease of the central nervous system treated with disease-modifying therapies, including the biologic, interferon-β (IFN-β). Up to 60% of IFN-β-exposed MS patients develop abnormal biochemical liver test results1,2, and 1 in 50 experiences drug-induced liver injury3. Since genomic variation contributes to other forms of drug-induced liver injury4,5, we aimed to identify biomarkers of IFN-β-induced liver injury using a two-stage genome-wide association study. The rs2205986 variant, previously linked to differential expression of IRF6, surpassed genome-wide significance in the combined two-stage analysis (P = 2.3 × 10–8, odds ratio = 8.3, 95% confidence interval = 3.6–19.2). Analysis of an independent cohort of IFN-β-treated MS patients identified via electronic medical records showed that rs2205986 was also associated with increased peak levels of aspartate aminotransferase (P = 7.6 × 10–5) and alkaline phosphatase (P = 4.9 × 10-4). We show that these findings may be applicable to predicting IFN-β-induced liver injury, offering insight into its safer use. A genome-wide study identifies common variation near IRF6 associated with liver injury induced by IFN-β in individuals treated for multiple sclerosis. The risk locus is associated with differential expression of IRF6 and may help guide safer use of this biologic.

Identifying rare genetic variants that drive the onset of disease is challenging, even before considering the additional genetic and environmental influences that likely exist in complex diseases. We recently published a study proposing a rare variant in the NR1H3 gene (p.R415Q, rs61731956) as responsible for the onset of multiple sclerosis (MS) in two multi-incident families (Wang et al., 2016). This publication has generated much discussion, and fortunately the possibility to validate a finding or prove it spurious can occur rapidly in genetic studies. All novel discoveries must be replicated, and best efforts should be made to ensure that these replications use the appropriate samples and approach, and provide the correct interpretation of the results. This Matters Arising Response paper addresses the Minikel and MacArthur (2016) and The International Multiple Sclerosis Genetics Consortium (2016) Matters Arising papers, published concurrently in Neuron.

Differences in response to medications have a strong genetic component. By leveraging publically available data, the spectrum of such genomic variation can be investigated extensively. Pharmacogenomic variation was extracted from the 1000 Genomes Project Phase 3 data (2504 individuals, 26 global populations). A total of 12 084 genetic variants were found in 120 pharmacogenes, with the majority (90.0%) classified as rare variants (global minor allele frequency 0.5) for one or more super-population comparison. A median of three clinical variants (PharmGKB level 1A/B) was found per individual, and 55.4% of individuals carried loss-of-function variants, varying by super-population (East Asian 60.9%>African 60.1%>South Asian 60.3%>European 49.3%>Admixed 39.2%). Genome sequencing can therefore identify clinical pharmacogenomic variation, and future studies need to consider rare variation to understand the spectrum of genetic diversity contributing to drug response.

BACKGROUND Although antipsychotics are integral to the treatment of schizophrenia, drug efficacy varies between patients. Although it has been shown that antipsychotic treatment response outcomes are heritable, our understanding of the genetic factors that are involved remains incomplete. Therefore, this study aims to use an unbiased scan of the genome to identify the genetic variants contributing toward antipsychotic treatment response outcomes. MATERIALS AND METHODS This study utilized whole-exome sequencing of patients on extreme ends of the treatment response spectrum (n=11) in combination with results from previous antipsychotic studies to design a panel of variants that were genotyped in two well-characterized first-episode schizophrenia cohorts (n=103 and 87). Association analyses were carried out to determine whether these variants were significantly associated with antipsychotic treatment response outcomes. RESULTS Association analyses in the discovery cohort identified two nonsynonymous variants that were significantly associated with antipsychotic treatment response outcomes (P

Multiple sclerosis (MS) is an inflammatory disease characterized by myelin loss and neuronal dysfunction. Despite the aggregation observed in some families, pathogenic mutations have remained elusive. In this study, we describe the identification of NR1H3 p.Arg415Gln in seven MS patients from two multi-incident families presenting severe and progressive disease, with an average age at onset of 34 years. Additionally, association analysis of common variants in NR1H3 identified rs2279238 conferring a 1.35-fold increased risk of developing progressive MS. The p.Arg415Gln position is highly conserved in orthologs and paralogs, and disrupts NR1H3 heterodimerization and transcriptional activation of target genes. Protein expression analysis revealed that mutant NR1H3 (LXRA) alters gene expression profiles, suggesting a disruption in transcriptional regulation as one of the mechanisms underlying MS pathogenesis. Our study indicates that pharmacological activation of LXRA or its targets may lead to effective treatments for the highly debilitating and currently untreatable progressive phase of MS.

South Africa, like many other developing countries, stands to benefit from novel diagnostics and drugs developed by pharmacogenomics guidance due to high prevalence of disease burden in the region. This includes both communicable (e.g., HIV/AIDS and tuberculosis) and non-communicable (e.g., diabetes and cardiovascular) diseases. For example, although only 0.7% of the world's population lives in South Africa, the country carries 17% of the global HIV/AIDS burden and 5% of the global tuberculosis burden. Nobel Peace Prize Laureate Archbishop Emeritus Desmond Tutu has coined the term Rainbow Nation, referring to a land of wealth in its many diverse peoples and cultures. It is now timely and necessary to reflect on how best to approach new genomics biotechnologies in a manner that carefully considers the public health needs and extant disease burden in the region. The aim of this paper is to document and review the advances in pharmacogenomics in South Africa and importantly, to evaluate the direction that future research should take. Previous research has shown that the populations in South Africa exhibit unique allele frequencies and novel genetic variation in pharmacogenetically relevant genes, often differing from other African and global populations. The high level of genetic diversity, low linkage disequilibrium and the presence of rare variants in these populations question the feasibility of the use of current commercially available genotyping platforms, and may partially account for genotype- phenotype discordance observed in past studies. However, the employment of high throughput technologies for genomic research, within the context of large clinical trials, combined with interdisciplinary studies and appropriate regulatory guidelines, should aid in acceleration of pharmacogenomic discoveries in high priority therapeutic areas in South Africa. Finally, we suggest that projects such as the H3Africa Initiative, the SAHGP and PGENI should play an integral role in the coordination of genomic research in South Africa, but also other African countries, by providing infrastructure and capital to local researchers, as well as providing aid in addressing the computational and statistical bottlenecks encountered at present.

Summary Genetic variation of the CYP2D6 gene has been associated with altered drug metabolism; however, limited studies have investigated CYP2D6 sequence diversity in African populations. We devised a CYP2D6 genotyping strategy to analyse the South African Xhosa population and genotype a Xhosa schizophrenia cohort, as CYP2D6 metabolises many antipsychotics and antidepressants. The entire CYP2D6 gene locus was sequenced in 15 Xhosa control individuals and the data generated were used to design a comprehensive genotyping strategy. Over 25 CYP2D6 alleles were genotyped in Xhosa controls and Xhosa schizophrenia patients using long-range PCR, DNA sequencing and single nucleotide primer extension analysis. Bioinformatic algorithms were used to predict the functional consequences of relevant mutations and samples were assigned CYP2D6 activity scores. A unique allele distribution was revealed and two rare novel alleles, CYP2D6*73 and CYP2D6*74, were identified. No significant differences in allele frequencies were detected between Xhosa controls and schizophrenia patients. This study provides i) comprehensive data on a poorly characterised population, ii) a valuable CYP2D6 genotyping strategy and iii) due to their unique genetic profile, provides the basis for pharmacogenetic intervention for Xhosa individuals.

Aim: To align predicted and measured CYP2C19 phenotype in a South African cohort. Materials & methods: Genotyping of CYP2C19*2, *3, *9, *15, *17, *27 and *28 was performed using PCR-RFLP, and an activity score (AS) system was used to predict phenotype. True phenotype was measured using plasma concentrations of omeprazole and its metabolite 5′-hydroxyomperazole. Results: Partial genotype-phenotype discrepancies were reported, and an adapted AS system was developed, which showed a marked improvement in phenotype prediction. Results highlight the need for a more comprehensive CYP2C19 genotyping approach to improve prediction of omeprazole metabolism. Conclusion: Evidence for the utility of a CYP2C19 AS system is provided, for which the accuracy can be further improved by means of comprehensive genotyping and substrate-specific modification.

Multiple sclerosis (MS) is an inflammatory disease characterized by myelin loss and neuronal dysfunction. Despite the aggregation observed in some families, pathogenic mutations have remained elusive. In this study we describe the identification of NR1H3 p.Arg415Gln in seven MS patients from two multi-incident families presenting severe and progressive disease, with an average age at onset of 34 years. Additionally, association analysis of common variants in NR1H3 identified rs2279238 conferring a 1.35-fold increased risk of developing progressive MS. The p.Arg415Gln position is highly conserved in orthologs and paralogs, and disrupts NR1H3 heterodimerization and transcriptional activation of target genes. Protein expression analysis revealed that mutant NR1H3 (LXRA) alters gene expression profiles, suggesting a disruption in transcriptional regulation as one of the mechanisms underlying MS pathogenesis. Our study indicates that pharmacological activation of LXRA or its targets may lead to effective treatments for the highly debilitating and currently untreatable progressive phase of MS.

Large-scale whole genome and exome resequencing studies have revealed that humans have a high level of deleterious rare variation, which has important implications for the design of future pharmacogenetics studies.Current pharmacogenetic guidelines focus on the implementation of common variation into dosing guidelines. However, it is becoming apparent that rare variation may also play an important role in differential drug response. Current sequencing technologies offer the opportunity to examine rare variation, but there are many challenges associated with such analyses. Nonetheless, if a comprehensive picture of the role that genetic variants play in treatment outcomes is to be obtained, it will be necessary to include the entire spectrum of variation, including rare variants, into pharmacogenetic research.In order to implement pharmacogenetics in the clinic, patients should be genotyped for clinically actionable pharmacogenetic variants and patients responding unfavourably to treatment after pharmacogenetics-based dosing should be identified and resequenced to identify additional functionally relevant variants, including rare variants. All derived information should be added to a central database to allow for the updating of existing dosing guidelines. By routinely implementing such strategies, pharmacogenetics-based treatment guidelines will continue to improve.

Aim: Many antipsychotic pharmacogenetics studies have been performed examining candidate genes or known variation; however, our understanding of the genetic factors involved in antipsychotic pharmacogenetic traits remains limited. Materials & methods: A well-characterized cohort of first-episode schizophrenia (FES) patients was used to identify a subset of nonresponders and responders to antipsychotic treatment for exome sequencing (n = 11). The variation observed in the responders and nonresponders was subsequently compared and a prioritization strategy was employed to identify variants for genotyping in the entire FES cohort (n = 103) as well as an additional Xhosa schizophrenia cohort (n = 222). Results: Examination of coding variation revealed a potential role for rare loss-of-function variants in treatment response outcomes. One variant, rs11368509, was found to be weakly associated with better treatment outcomes in the FES cohort (p = 0.057) and the Xhosa schizophrenia cohort (p = 0.016). In addition, the majority of the loss-of-function variation that was considered likely to be involved in antipsychotic treatment response was either novel or rare in Asian and European populations. Conclusion: This pilot study has highlighted the importance of exome sequencing for antipsychotic pharmacogenomics studies, particularly in African individuals. Furthermore, the results emphasize once again the complexity of antipsychotic pharmacogenomics and the need for future research. Original submitted 8 July 2013; Revision submitted 24 October 2013

Genetics and the clinical response to warfarin and edoxaban: findings from the randomized, double-blind ENGAGE AF-TIMI 48 trial Mega et al. (2015) Lancet

Because of the unmet needs of current pharmacotherapy for schizophrenia, antipsychotic pharmacogenetic research is of utmost importance. However, to date, few clinically applicable antipsychotic pharmacogenomic alleles have been identified. Nonetheless, next-generation sequencing technologies are expected to aid in the identification of clinically significant variants for this complex phenotype. The aim of this study was therefore to critically examine the ability of next-generation sequencing technologies to reliably detect variation present in pharmacogenes.Candidate antipsychotic pharmacogenes and very important pharmacogenes were identified from the literature and the Pharmacogenomics Knowledgebase. Thereafter, the percentage sequence similarity observed between these genes and their corresponding pseudogenes and paralogues, as well as the percentage low-complexity sequence and GC content of each gene, was calculated. These sequence attributes were subsequently compared with the 'inaccessible' regions of these genes as described by the 1000 Genomes Project.It was found that the percentage 'inaccessible genome' correlated well with GC content (P=9.96×10), low-complexity sequence (P=0.0002) and the presence of pseudogenes/paralogues (P=8.02×10). In addition, it was found that many of the pharmacogenes were not ideally suited to next-generation sequencing because of these genomic complexities. These included the CYP and HLA genes, both of which are of importance to many fields of pharmacogenetics.Current short read sequencing technologies are unable to comprehensively capture the variation in all pharmacogenes. Therefore, until high-throughput sequencing technologies advance further, it may be necessary to combine next-generation sequencing with other genotyping strategies.

The CYP3A4 enzyme is the most abundant human cytochrome P450 (CYP) and is regarded as the most important enzyme involved in drug metabolism. Inter-individual and inter-population variability in gene expression and enzyme activity are thought to be influenced, in part, by genetic variation. Although Southern African individuals have been shown to exhibit the highest levels of genetic diversity, they have been under-represented in pharmacogenetic research to date. Therefore, the aim of this study was to identify genetic variation within CYP3A4 in three South African population groups comprising of 29 Khoisan, 65 Xhosa and 65 Mixed Ancestry (MA) individuals. To identify known and novel CYP3A4 variants, 15 individuals were randomly selected from each of the population groups for bi-directional Sanger sequencing of ~600 bp of the 5'-upstream region and all thirteen exons including flanking intronic regions. Genetic variants detected were genotyped in the rest of the cohort. In total, 24 SNPs were detected, including CYP3A4(*)12, CYP3A4(*)15, and the reportedly functional CYP3A4(*)1B promoter polymorphism, as well as two novel non-synonymous variants. These putatively functional variants, p.R162W and p.Q200H, were present in two of the three populations and all three populations, respectively, and in silico analysis predicted that the former would damage the protein product. Furthermore, the three populations were shown to exhibit distinct genetic profiles. These results confirm that South African populations show unique patterns of variation in the genes encoding xenobiotic metabolizing enzymes. This research suggests that population-specific genetic profiles for CYP3A4 and other drug metabolizing genes would be essential to make full use of pharmacogenetics in Southern Africa. Further investigation is needed to determine if the identified genetic variants influence CYP3A4 metabolism phenotype in these populations.

Background Huntington’s disease (HD) is diagnosed in 1 in 7300 individuals in Western populations but the frequency and penetrance of the causative CAG repeat expansion is unknown. Effects of population ageing, which may increase the rate of late-onset HD, remain unclear. Aims To directly estimate the frequency and penetrance of CAG repeat alleles associated with HD, and model changes in prevalence resulting from increased ascertainment of late-onset cases. Methods CAG repeat length was evaluated in 7315 individuals from three population-based cohorts in British Columbia, the United States, and Scotland. The frequency of CAG 36–38 repeat genotypes was compared to the prevalence of HD patients with genetically confirmed CAG 36–38 in a multisource clinical ascertainment in British Columbia, Canada. Penetrance of 36–38 CAG repeat alleles for HD was directly estimated for individuals ≥65 years of age. Age-specific prevalence rates were used to model change in prevalence as a function of population ageing and increased ascertainment of patients ≥65 years of age. Results 18 of 7315 individuals had ≥36 CAG, revealing that approximately 1 in 400 individuals in the general population have an expanded CAG repeat associated with HD (0.246%). Individuals with CAG 36–37 genotypes are the most common (36, 0.096%; 37, 0.082%; 38, 0.027%; 39, 0.000%; ≥40, 0.041%). The prevalence of HD is expected to increase as a result of both population ageing and increased ascertainment of late-onset cases. Conclusions The relatively infrequent diagnosis of HD at 36–38 CAG repeats suggests low penetrance in this range. Another contributing factor may be reduced ascertainment of HD in those of older age. Our data imply that population ageing will lead to higher prevalence rates of HD. Improved ascertainment of late-onset HD, particularly in the reduced penetrance range, may lead to further increases.

The catechol-O-methyltransferase (COMT) gene is an attractive schizophrenia candidate gene, encoding a catabolic dopamine enzyme. The enzyme exists as two distinct isoforms, with the membrane bound enzyme (i.e. MB-COMT) being predominantly expressed in the brain. Since African populations remain underrepresented in genetic/genomic research, we performed an association study to determine whether MB-COMT genetic variants are associated with schizophrenia-susceptibility and symptom severity in the South African Xhosa population. Fourteen candidate polymorphisms were selected by means of a literature search and in silico analyses and were subsequently genotyped in a cohort of 238 Xhosa schizophrenia patients and 240 healthy Xhosa controls. Genetic association was tested with schizophrenia-susceptibility as well as symptom severity within the patient group. Polymorphisms of interest were also analysed using functional assays. Two SNPs, rs2020917 (OR=0.54, 95% CI 0.37-0.79; P=0.0011) and rs737865 (OR=0.52, 95% CI 0.36-0.74; P=0.0002), in the P2 promoter region were significantly associated with schizophrenia as well as an increase (increase=11.2%, 95% CI 3.7%-19.2%; P=0.0031) in reporter gene expression. The minor alleles of these SNPs were underrepresented in the schizophrenia cohort, indicating a possible protective effect. The P2 region also formed part of a haplotype found to be associated with the severity of the negative symptoms of the disorder. The data generated by this study indicate that genetic variation of MB-COMT could be associated with schizophrenia and negative symptom severity in the Xhosa population and may therefore be one of the genomic loci contributing towards the disorder in the South African community. Future large-scale studies in other African schizophrenia populations are required to further elucidate the significance of these findings.

Vural Ozdemir1,*, Erik Fisher2, Edward S. Dove1, Hilary Burton3, Galen E. B. Wright4, Mario Masellis5,*, and Louise Warnich4,* 1Centre of Genomics and Policy, Department of Human Genetics, Faculty of Medicine, McGill University, Montreal, QC, Canada 2School of Politics and Global Studies, Consortium for Science, Policy & Outcomes, Center for Nanotechnology in Society, College of Liberal Arts and Sciences, Arizona State University, Tempe, AZ, USA 3Foundation for Genomics and Population Health, 2 Worts Causeway, Cambridge, UK 4Department of Genetics, Stellenbosch University, Stellenbosch, South Africa 5L.C. Campbell Cognitive Neurology Research Unit, Sunnybrook Research Institute, Department of Medicine (Neurology), Sunnybrook Health Sciences Centre, University of Toronto, Toronto, ON, Canada