Characterization of CYP2B6 and CYP2A6 Pharmacogenetic Variation in Sub-Saharan African Populations

Genetic variation in CYP2B6 and CYP2A6 is known to impact interindividual response to antiretrovirals, nicotine, and bupropion, among other drugs. However, the full catalogue of clinically relevant pharmacogenetic variants in these genes is yet to be established, especially across African populations. This study therefore aimed to characterize the star allele (haplotype) distribution in CYP2B6 and CYP2A6 across diverse and understudied sub-Saharan African (SSA) populations. We called star alleles from 961 high-depth full genomes using StellarPGx, Aldy, and PyPGx. In addition, we performed CYP2B6 and CYP2A6 star allele frequency comparisons between SSA and other global biogeographical groups represented in the new 1000 Genomes Project high-coverage dataset (n = 2,000). This study presents frequency information for star alleles in CYP2B6 (e.g., *6 and *18; frequency of 21-47% and 2-19%, respectively) and CYP2A6 (e.g., *4, *9, and *17; frequency of 0-6%, 3-10%, and 6-20%, respectively), and predicted phenotypes (for CYP2B6), across various African populations. In addition, 50 potentially novel African-ancestry star alleles were computationally predicted by StellarPGx in CYP2B6 and CYP2A6 combined. For each of these genes, over 4% of the study participants had predicted novel star alleles. Three novel star alleles in CYP2A6 (*54, *55, and *56) and CYP2B6 apiece, and several suballeles were further validated via targeted Single-Molecule Real-Time resequencing. Our findings are important for informing the design of comprehensive pharmacogenetic testing platforms, and are highly relevant for personalized medicine strategies, especially relating to antiretroviral medication and smoking cessation treatment in Africa and the African diaspora. More broadly, this study highlights the importance of sampling diverse African ethnolinguistic groups for accurate characterization of the pharmacogene variation landscape across the continent.

CIAO1 and MMS19 deficiency: A lethal neurodegenerative phenotype caused by cytosolic Fe-S cluster protein assembly disorders

PURPOSE

The functionality of many cellular proteins depends on cofactors; yet, they have only been implicated in a minority of Mendelian diseases. Here, we describe the first 2 inherited disorders of the cytosolic iron-sulfur protein assembly system.

METHODS

Genetic testing via genome sequencing was applied to identify the underlying disease cause in 3 patients with microcephaly, congenital brain malformations, progressive developmental and neurologic impairments, recurrent infections, and a fatal outcome. Studies in patient-derived skin fibroblasts and zebrafish models were performed to investigate the biochemical and cellular consequences.

RESULTS

Metabolic analysis showed elevated uracil and thymine levels in body fluids but no pathogenic variants in DPYD, encoding dihydropyrimidine dehydrogenase. Genome sequencing identified compound heterozygosity in 2 patients for missense variants in CIAO1, encoding cytosolic iron-sulfur assembly component 1, and homozygosity for an in-frame 3-nucleotide deletion in MMS19, encoding the MMS19 homolog, cytosolic iron-sulfur assembly component, in the third patient. Profound alterations in the proteome, metabolome, and lipidome were observed in patient-derived fibroblasts. We confirmed the detrimental effect of deficiencies in CIAO1 and MMS19 in zebrafish models.

CONCLUSION

A general failure of cytosolic and nuclear iron-sulfur protein maturation caused pleiotropic effects. The critical function of the cytosolic iron-sulfur protein assembly machinery for antiviral host defense may well explain the recurrent severe infections occurring in our patients.

Characterization of CYP2D6 Pharmacogenetic Variation in Sub-Saharan African Populations

Cytochrome P450 2D6 (CYP2D6) is a key enzyme in drug response owing to its involvement in the metabolism of ~ 25% of clinically prescribed medications. The encoding CYP2D6 gene is highly polymorphic, and many pharmacogenetics studies have been performed worldwide to investigate the distribution of CYP2D6 star alleles (haplotypes); however, African populations have been relatively understudied to date. In this study, the distributions of CYP2D6 star alleles and predicted drug metabolizer phenotypes-derived from activity scores-were examined across multiple sub-Saharan African populations based on bioinformatics analysis of 961 high-depth whole genome sequences. This was followed by characterization of novel star alleles and suballeles in a subset of the participants via targeted high-fidelity Single-Molecule Real-Time resequencing (Pacific Biosciences). This study revealed varying frequencies of known CYP2D6 alleles and predicted phenotypes across different African ethnolinguistic groups. Twenty-seven novel CYP2D6 star alleles were predicted computationally and two of them were further validated. This study highlights the importance of studying variation in key pharmacogenes such as CYP2D6 in the African context to better understand population-specific allele frequencies. This will aid in the development of better genotyping panels and star allele detection approaches with a view toward supporting effective implementation of precision medicine strategies in Africa and across the African diaspora.

Transcriptome-wide association analyses identify an association between ARL14EP and polycystic ovary syndrome

Polycystic ovary syndrome (PCOS) is a common endocrine disorder, which is accompanied by a variety of comorbidities including metabolic, reproductive, and psychiatric disorders. Genome-wide association studies have identified several genetic variants that are associated with PCOS. However, these variants often occur outside of coding regions and require further investigation to understand their contribution to PCOS. A transcriptome-wide association study (TWAS) was performed to uncover heritable gene expression profiles that are associated with PCOS in two independent cohorts. Causal gene prioritization was subsequently performed and expression of genes prioritized through these analyses was examined in 49 PCOS patients and 30 controls. TWAS analyses revealed that increased expression of ARL14EP was significantly associated with PCOS risk in the discovery (P = 1.6 × 10^-6) and replication cohorts (P = 2.0 × 10^-13). Gene prioritization pipelines provided further evidence that ARL14EP is the most likely causal gene at this locus. ARL14EP gene expression was shown to be significantly different between PCOS cases and controls, after adjusting for body mass index, age and testosterone levels (P = 1.2 × 10^-13). This study has provided evidence for the role of ARL14EP in PCOS. Given that ARL14EP has been reported to play an important role in chromatin remodeling, variants affecting the expression of ARL14EP may also affect the expression of other genes that contribute to PCOS pathogenesis.

HLA‐DPB1 is associated with ANCA ‐associated vasculitis in children

OBJECTIVE

Antineutrophil cytoplasmic antibody (ANCA)-associated vasculitis (AAV) is a rare, life-threatening inflammation of blood vessels that can affect both adults and children. Compared to adult-onset disease, AAV is especially rare in children, with an annual prevalence of 0.5-6.4 cases per million children. The etiology of AAV remains largely unknown, and both environmental and genetic factors are likely involved. The present study was undertaken to explore the genetic susceptibility factors recently identified in adult patients, including HLA-DP and HLA-DQ, in pediatric patients.

METHODS

We performed a genome-wide association study of pediatric AAV in patients of European ancestry (n = 63 AAV cases, n = 315 population-matched controls).

RESULTS

We identified a significant genetic association between pediatric AAV and the HLA-DPB1*04:01 allele (P = 1.5 × 10^-8 , odds ratio [OR] 3.5), with a stronger association observed in children with proteinase 3-ANCA positivity than in children with myeloperoxidase-ANCA positivity. Among the HLA alleles, the HLA-DPB1*04:01 allele was the most highly associated with AAV, although not significantly, in a follow-up adult AAV cohort (P = 2.6 × 10^-4 , OR 0.4). T cell receptor and interferon signaling pathways were also shown to be enriched in the pediatric AAV cohort.

CONCLUSION

The HLA-DPB1 locus showed an association with pediatric AAV, as similarly shown previously in adult AAV. Despite the difference in the age of onset, these findings suggest that childhood- and adult-onset vasculitis share a common genetic predisposition. The identification of genetic variants contributing to AAV is an important step to improved classification tools and treatment strategies.

PharmVar GeneFocus: CYP3A5

The Pharmacogene Variation Consortium (PharmVar) catalogs star (*) allele nomenclature for the polymorphic human CYP3A5 gene. Genetic variation within the CYP3A5 gene locus impacts the metabolism of several clinically important drugs, including the immunosuppressants tacrolimus, sirolimus, cyclosporine, and the benzodiazepine midazolam. Variable CYP3A5 activity is of clinical importance regarding tacrolimus metabolism. This GeneFocus provides a CYP3A5 gene summary with a focus on aspects regarding standardized nomenclature. In addition, this review also summarizes recent changes and updates, including the retirement of several allelic variants and provides an overview of how PharmVar CYP3A5 star allele nomenclature is utilized by the Pharmacogenomics Knowledgebase (PharmGKB) and the Clinical Pharmacogenetics Implementation Consortium (CPIC).

Patient-specific genetic factors predict treatment failure in sofosbuvir-treated patients with chronic hepatitis C

BACKGROUND & AIMS

According to pivotal clinical trials, cure rates for sofosbuvir-based antiviral therapy exceed 96%. Treatment failure is usually assumed to be because of virological resistance-associated substitutions or clinical risk factors, yet the role of patient-specific genetic factors has not been well explored. We determined if patient-specific genetic factors help predict patients likely to fail sofosbuvir treatment in real-world treatment situations.

METHODS

We recruited sofosbuvir-treated patients with chronic hepatitis C from five Canadian treatment sites, and performed a case-control pharmacogenomics study assessing both previously published and novel genetic polymorphisms. Specifically studied were variants predicted to impair CES1-dependent production of sofosbuvir's active metabolite, interferon-λ signalling variants expected to impact a patient's immune response to the virus and an HLA variant associated with increased spontaneous and treatment-induced viral clearance.

RESULTS

Three hundred and fifty-nine sofosbuvir-treated patients were available for analyses after exclusions, with 34 (9.5%) failing treatment. We identified CES1 variants as novel predictors for treatment failure in European patients (rs115629050 or rs4513095; odds ratio (OR): 5.43; 95% confidence interval (CI): 1.64-18.01; P = .0057), replicated associations with IFNL4 variants predicted to increase interferon-λ signalling (eg rs12979860; OR: 2.25; 95% CI: 1.25-4.06; P = .0071) and discovered a novel association with a coding variant predicted to enhance the activity of IFNL4's receptor (rs2834167 in IL10RB; OR: 1.81; 95% CI: 1.01-3.24; P = .047).

CONCLUSIONS

Ultimately, this work demonstrates that patient-specific genetic factors could be used as a tool to identify patients at higher risk of treatment failure and allow for these patients to receive effective therapy sooner.

A pharmacogenomic investigation of the cardiac safety profile of ondansetron in children and pregnant women

BACKGROUND

Ondansetron is a highly effective antiemetic for the treatment of nausea and vomiting. However, this medication has also been associated with QT prolongation. Pharmacogenomic information on therapeutic response to ondansetron exists, but no investigation has been performed on genetic factors that influence the cardiac safety of this medication.

METHODS

Three patient groups receiving ondansetron were recruited and followed prospectively (pediatric post-surgical patients n = 101; pediatric oncology patients n = 98; pregnant women n = 62). Electrocardiograms were conducted at baseline, and 5- and 30-min post-ondansetron administration, to determine the effect of ondansetron treatment on QT interval. Pharmacogenomic associations were assessed via analyses of comprehensive CYP2D6 genotyping and genome-wide association study data.

RESULTS

In the entire cohort, 62 patients (24.1%) met the criteria for prolonged QT, with 1.2% of the cohort exhibiting unsafe QT prolongation. The most significant shift from baseline occurred at five minutes post-ondansetron administration (P = 9.8 × 10^-4). CYP2D6 activity score was not associated with prolonged QT. Genome-wide analyses identified novel associations with a missense variant in TLR3 (rs3775291; P = 2.00 × 10^-7) and a variant linked to the expression of SLC36A1 (rs34124313; P = 1.97 × 10^-7).

CONCLUSIONS

This study has provided insight into the genomic basis of ondansetron-induced cardiac changes and has emphasized the importance of genes that have been implicated in serotonin-related traits. These biologically-relevant findings represent the first step towards understanding this adverse event with the overall goal to improve the safety of this commonly used antiemetic medication.

A systematic review and analysis of the use of polygenic scores in pharmacogenomics

Polygenic scores (PGS) 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 PGS 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 PGS 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 PGS 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 PGS, with between 3 and 6.6 million variants included in the PGS. 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 PGS deposited on the PGS Catalog. These findings highlight current gaps and key areas for future pharmacogenomic PGS research.

Novel variant in glycophorin c gene protects against ribavirin-induced anemia during chronic hepatitis C treatment

BACKGROUND

The current use of ribavirin in difficult-to-cure chronic hepatitis C patients (HCV) and patients with severe respiratory infections is constrained by the issue of ribavirin-induced hemolytic anemia that affects 30% of treated patients, requiring dosage modification or discontinuation. Though some genetic variants have been identified predicting this adverse effect, known clinical and genetic factors do not entirely explain the risk of ribavirin-induced anemia.

METHODS

We assessed the associations of previously identified variants in inosine triphosphatase (ITPA), solute carrier 28A2 (SLC28A2) and vitamin D receptor (VDR) genes with ribavirin-induced anemia defined as hemoglobin decline of ≥30 g/L on treatment, followed by a staged discovery (n = 114), replication (n = 74), and combined (n = 188) genome-wide association study to uncover potential new predictive variants.

RESULTS

We identified a novel association in the gene coding glycophorin C (rs6741425; OR:0.12, 95%CI:0.04-0.34, P = 2.94 × 10^-6) that predicts protection against ribavirin-induced anemia. We also replicated the associations of ITPA and VDR genetic variants with the development of ribavirin-induced anemia (rs1127354; OR:0.13, 95%CI:0.04-0.41, P = 8.66 ×10^-5; and rs1544410; OR:1.65, 95%CI:1.01-2.70, P = 0.0437).

CONCLUSIONS

GYPC variation affecting erythrocyte membrane strength is important in predicting risk for developing ribavirin-induced anemia. ITPA and VDR genetic variants are also important predictors of this adverse reaction.

TCERG1L allelic variation is associated with cisplatin-induced hearing loss in childhood cancer, a PanCareLIFE study

In children with cancer, the heterogeneity in ototoxicity occurrence after similar treatment suggests a role for genetic susceptibility. Using a genome-wide association study (GWAS) approach, we identified a genetic variant in TCERG1L (rs893507) to be associated with hearing loss in 390 non-cranial irradiated, cisplatin-treated children with cancer. These results were replicated in two independent, similarly treated cohorts (n = 192 and 188, respectively) (combined cohort: P = 5.3 × 10-10, OR 3.11, 95% CI 2.2-4.5). Modulating TCERG1L expression in cultured human cells revealed significantly altered cellular responses to cisplatin-induced cytokine secretion and toxicity. These results contribute to insights into the genetic and pathophysiological basis of cisplatin-induced ototoxicity.

Transcriptome-wide association study uncovers the role of essential genes in anthracycline-induced cardiotoxicity

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.

The Extent and Impact of Variation in ADME Genes in Sub-Saharan African Populations

Introduction: Investigating variation in genes involved in the absorption, distribution, metabolism, and excretion (ADME) of drugs are key to characterizing pharmacogenomic (PGx) relationships. ADME gene variation is relatively well characterized in European and Asian populations, but data from African populations are under-studied-which has implications for drug safety and effective use in Africa. Results: We identified significant ADME gene variation in African populations using data from 458 high-coverage whole genome sequences, 412 of which are novel, and from previously available African sequences from the 1,000 Genomes Project. ADME variation was not uniform across African populations, particularly within high impact coding variation. Copy number variation was detected in 116 ADME genes, with equal ratios of duplications/deletions. We identified 930 potential high impact coding variants, of which most are discrete to a single African population cluster. Large frequency differences (i.e., >10%) were seen in common high impact variants between clusters. Several novel variants are predicted to have a significant impact on protein structure, but additional functional work is needed to confirm the outcome of these for PGx use. Most variants of known clinical outcome are rare in Africa compared to European populations, potentially reflecting a clinical PGx research bias to European populations. Discussion: The genetic diversity of ADME genes across sub-Saharan African populations is large. The Southern African population cluster is most distinct from that of far West Africa. PGx strategies based on European variants will be of limited use in African populations. Although established variants are important, PGx must take into account the full range of African variation. This work urges further characterization of variants in African populations including in vitro and in silico studies, and to consider the unique African ADME landscape when developing precision medicine guidelines and tools for African populations.

Integrating disease and drug-related phenotypes for improved identification of pharmacogenomic variants

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.

StellarPGx: A Nextflow Pipeline for Calling Star Alleles in Cytochrome P450 Genes

Bioinformatics pipelines for calling star alleles (haplotypes) in cytochrome P450 (CYP) genes are important for the implementation of precision medicine. Genotyping CYP genes using high throughput sequencing data is complicated, e.g., by being highly polymorphic, not to mention the structural variations especially in CYP2D6, CYP2A6, and CYP2B6. Genome graph-based variant detection approaches have been shown to be reliable for genotyping HLA alleles. However, their application to enhancing star allele calling in CYP genes has not been extensively explored. We present StellarPGx, a Nextflow pipeline for accurately genotyping CYP genes by combining genome graph-based variant detection, read coverage information from the original reference-based alignments, and combinatorial diplotype assignments. The implementation of StellarPGx using Nextflow facilitates its portability, reproducibility, and scalability on various user platforms. StellarPGx is currently able to genotype 12 important pharmacogenes belonging to the CYP1, 2, and 3 families. For purposes of validation, we use CYP2D6 as a model gene owing to its high degree of polymorphisms (over 130 star alleles defined to date, including complex structural variants) and clinical importance. We applied StellarPGx and three existing callers to 109 whole genome sequenced samples for which the Genetic Testing Reference Material Coordination Program (GeT-RM) has recently provided consensus truth CYP2D6 diplotypes. StellarPGx had the highest CYP2D6 diplotype concordance (99%) with GeT-RM compared with Cyrius (98%), Aldy (82%), and Stargazer (84%). This exemplifies the high accuracy of StellarPGx and highlights its importance for both research and clinical pharmacogenomics applications. The StellarPGx pipeline is open-source and available from https://github.com/SBIMB/StellarPGx.

Bi-allelic GOT2 Mutations Cause a Treatable Malate-Aspartate Shuttle-Related Encephalopathy

Early-infantile encephalopathies with epilepsy are devastating conditions mandating an accurate diagnosis to guide proper management. Whole-exome sequencing was used to investigate the disease etiology in four children from independent families with intellectual disability and epilepsy, revealing bi-allelic GOT2 mutations. In-depth metabolic studies in individual 1 showed low plasma serine, hypercitrullinemia, hyperlactatemia, and hyperammonemia. The epilepsy was serine and pyridoxine responsive. Functional consequences of observed mutations were tested by measuring enzyme activity and by cell and animal models. Zebrafish and mouse models were used to validate brain developmental and functional defects and to test therapeutic strategies. GOT2 encodes the mitochondrial glutamate oxaloacetate transaminase. GOT2 enzyme activity was deficient in fibroblasts with bi-allelic mutations. GOT2, a member of the malate-aspartate shuttle, plays an essential role in the intracellular NAD(H) redox balance. De novo serine biosynthesis was impaired in fibroblasts with GOT2 mutations and GOT2-knockout HEK293 cells. Correcting the highly oxidized cytosolic NAD-redox state by pyruvate supplementation restored serine biosynthesis in GOT2-deficient cells. Knockdown of got2a in zebrafish resulted in a brain developmental defect associated with seizure-like electroencephalography spikes, which could be rescued by supplying pyridoxine in embryo water. Both pyridoxine and serine synergistically rescued embryonic developmental defects in zebrafish got2a morphants. The two treated individuals reacted favorably to their treatment. Our data provide a mechanistic basis for the biochemical abnormalities in GOT2 deficiency that may also hold for other MAS defects.

Length of Uninterrupted CAG, Independent of Polyglutamine Size, Results in Increased Somatic Instability, Hastening Onset of Huntington Disease

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.

Outcome Definition Influences the Relationship Between Genetic Polymorphisms of ERCC1, ERCC2, SLC22A2 and Cisplatin Nephrotoxicity in Adult Testicular Cancer Patients

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-2014, aged >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-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-9.88; p < 0.001) while genotype AC was associated with a higher risk of cisplatin nephrotoxicity (ORadj = 5.06; 95% CI:1.69-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.

Glutaminase Deficiency Caused by Short Tandem Repeat Expansion in GLS

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.

Analyses of Adverse Drug Reactions-Nationwide Active Surveillance Network: Canadian Pharmacogenomics Network for Drug Safety Database

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.

Pharmacogenomics of Vincristine-Induced Peripheral Neuropathy Implicates Pharmacokinetic and Inherited Neuropathy Genes

Vincristine is an effective chemotherapeutic drug for various cancers, including acute lymphoblastic leukemia (ALL). Unfortunately, clinical utility is restricted by dose‐limiting vincristine‐induced peripheral neuropathies (VIPN). We sought to determine the association of VIPN with a recently identified risk variant, CEP72 rs924607, and drug absorption, distribution, metabolism, and excretion (ADME) gene variants in pediatric ALL. This was followed by a meta‐analysis of pharmacogenomic data from over 500 patients. CEP72 rs924607 was significantly associated with VIPN (P = 0.02; odds ratio (OR) = 3.4). ADME analyses identified associations between VIPN and ABCC1 rs3784867 (P = 5.34 × 10⁻⁵; OR = 4.9), and SLC5A7 rs1013940 (P = 9.00 × 10⁻⁴; OR= 8.6); genes involved in vincristine transport and inherited neuropathies, respectively. Meta‐analysis identified an association with a variant related to TTPA (rs10504361: P = 6.85 × 10⁻⁴; OR = 2.0), a heritable neuropathy‐related gene. This study provides essential corroboratory evidence for CEP72 rs924607 and highlights the importance of drug transporter and inherited neuropathy genes in VIPN.

Gain-of-function KCNJ6 Mutation in a Severe Hyperkinetic Movement Disorder Phenotype

Here, we describe a fourth case of a human with a de novo KCNJ6 (GIRK2) mutation, who presented with clinical findings of severe hyperkinetic movement disorder and developmental delay, similar to the Keppen-Lubinsky syndrome but without lipodystrophy. Whole-exome sequencing of the patient's DNA revealed a heterozygous de novo variant in the KCNJ6 (c.512T>G, p.Leu171Arg). We conducted in vitro functional studies to determine if this Leu-to-Arg mutation alters the function of GIRK2 channels. Heterologous expression of the mutant GIRK2 channel alone produced an aberrant basal inward current that lacked G protein activation, lost K⁺ selectivity and gained Ca²⁺ permeability. Notably, the inward current was inhibited by the Na⁺ channel blocker QX-314, similar to the previously reported weaver mutation in murine GIRK2. Expression of a tandem dimer containing GIRK1 and GIRK2(p.Leu171Arg) did not lead to any currents, suggesting heterotetramers are not functional. In neurons expressing p.Leu171Arg GIRK2 channels, these changes in channel properties would be expected to generate a sustained depolarization, instead of the normal G protein-gated inhibitory response, which could be mitigated by expression of other GIRK subunits. The identification of the p.Leu171Arg GIRK2 mutation potentially expands the Keppen-Lubinsky syndrome phenotype to include severe dystonia and ballismus. Our study suggests screening for dominant KCNJ6 mutations in the evaluation of patients with severe movement disorders, which could provide evidence to support a causal role of KCNJ6 in neurological channelopathies.

Further Investigation of the Role of ACYP2 and WFS1 Pharmacogenomic Variants in the Development of Cisplatin-Induced Ototoxicity in Testicular Cancer Patients

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.

Association Between SLC16A5 Genetic Variation and Cisplatin-Induced Ototoxic Effects in Adult Patients With Testicular Cancer

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, rs4788863 in 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.

A de novo mosaic mutation in SPAST with two novel alternative alleles and chromosomal copy number variant in a boy with spastic paraplegia and autism spectrum disorder

Here we report a 12 year old male with an extreme presentation of spastic paraplegia along with autism and dysmorphisms. Whole exome sequencing identified a predicted pathogenic pair of missense variants in SPAST at the same chromosomal location, each with a different alternative allele, while a chromosome microarray identified a 1.73 Mb paternally inherited copy gain of 1q21.1q21.2 resulting in a blended phenotype of both Spastic paraplegia 4 and 1q21.1 microduplication syndrome. We believe that the extreme phenotype observed is likely caused by the presence of cells which contain only mutant SPAST, but that the viability of the patient is possible due mosaicism of mutant alleles observed in different proportions across tissues.

Identification of a large intronic transposal insertion in SLC17A5 causing sialic acid storage disease

BACKGROUND

Sialic acid storage diseases are neurodegenerative disorders characterized by accumulation of sialic acid in the lysosome. These disorders are caused by mutations in SLC17A5, the gene encoding sialin, a sialic acid transporter located in the lysosomal membrane. The most common form of sialic acid storage disease is the slowly progressive Salla disease, presenting with hypotonia, ataxia, epilepsy, nystagmus and findings of cerebral and cerebellar atrophy. Hypomyelination and corpus callosum hypoplasia are typical as well. We report a 16 year-old boy with an atypically mild clinical phenotype of sialic acid storage disease characterized by psychomotor retardation and a mixture of spasticity and rigidity but no ataxia, and only weak features of hypomyelination and thinning of corpus callosum on MRI of the brain.

RESULTS

The thiobarbituric acid method showed elevated levels of free sialic acid in urine and fibroblasts, indicating sialic acid storage disease. Initial Sanger sequencing of SLC17A5 coding regions did not show any pathogenic variants, although exon 9 could not be sequenced. Whole exome sequencing followed by RNA and genomic DNA analysis identified a homozygous 6040 bp insertion in intron 9 of SLC17A5 corresponding to a long interspersed element-1 retrotransposon (KF425758.1). This insertion adds two splice sites, both resulting in a frameshift which in turn creates a premature stop codon 4 bp into intron 9.

CONCLUSIONS

This study describes a novel pathogenic variant in SLC17A5, namely an intronic transposal insertion, in a patient with mild biochemical and clinical phenotypes. The presence of a small fraction of normal transcript may explain the mild phenotype. This case illustrates the importance of including lysosomal sialic acid storage disease in the differential diagnosis of developmental delay with postnatal onset and hypomyelination, as well as intronic regions in the genetic investigation of inborn errors of metabolism.

Fine-mapping of antipsychotic response genome-wide association studies reveals novel regulatory mechanisms

Aim: Noncoding variation has demonstrated regulatory effects on disease treatment outcomes. This study investigated the potential functionality of previously implicated noncoding variants on schizophrenia treatment response. Materials & methods: Predicted regulatory potential of variation identified from antipsychotic response genome-wide association studies was determined. Prioritized variants were assessed for association(s) with treatment outcomes in a South African first episode schizophrenia cohort (n = 103). Results: Bioinformatic and association results implicated a relationship between regulatory variants, expression of MANBA, COL9A2 and NFKB1, and treatment response. Three SNPs were associated with poor outcomes (rs230493: p = 1.88 × 10-6; rs3774959: p = 1.75 × 10-5; and rs230504: p = 1.48 × 10-4). Conclusion: This study has thoroughly investigated previous GWAS to pinpoint variants that may play a causal role in poor schizophrenia treatment outcomes, and provides potential candidate genes for further study in the field of antipsychotic response.

Exome Sequencing and the Management of Neurometabolic Disorders

BACKGROUND

Whole-exome sequencing has transformed gene discovery and diagnosis in rare diseases. Translation into disease-modifying treatments is challenging, particularly for intellectual developmental disorder. However, the exception is inborn errors of metabolism, since many of these disorders are responsive to therapy that targets pathophysiological features at the molecular or cellular level.

METHODS

To uncover the genetic basis of potentially treatable inborn errors of metabolism, we combined deep clinical phenotyping (the comprehensive characterization of the discrete components of a patient's clinical and biochemical phenotype) with whole-exome sequencing analysis through a semiautomated bioinformatics pipeline in consecutively enrolled patients with intellectual developmental disorder and unexplained metabolic phenotypes.

RESULTS

We performed whole-exome sequencing on samples obtained from 47 probands. Of these patients, 6 were excluded, including 1 who withdrew from the study. The remaining 41 probands had been born to predominantly nonconsanguineous parents of European descent. In 37 probands, we identified variants in 2 genes newly implicated in disease, 9 candidate genes, 22 known genes with newly identified phenotypes, and 9 genes with expected phenotypes; in most of the genes, the variants were classified as either pathogenic or probably pathogenic. Complex phenotypes of patients in five families were explained by coexisting monogenic conditions. We obtained a diagnosis in 28 of 41 probands (68%) who were evaluated. A test of a targeted intervention was performed in 18 patients (44%).

CONCLUSIONS

Deep phenotyping and whole-exome sequencing in 41 probands with intellectual developmental disorder and unexplained metabolic abnormalities led to a diagnosis in 68%, the identification of 11 candidate genes newly implicated in neurometabolic disease, and a change in treatment beyond genetic counseling in 44%. (Funded by BC Children's Hospital Foundation and others.).

CYP2B6*6 and CYP2B6*18 Predict Long-Term Efavirenz Exposure Measured in Hair Samples in HIV-Positive South African Women

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.

Evaluation of predictive CYP2C19 genotyping assays relative to measured phenotype in a South African cohort

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.

Considerations for rare variants in drug metabolism genes and the clinical implications

INTRODUCTION

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.

AREAS COVERED

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.

EXPERT OPINION

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.

Introduction of the AmpliChip CYP450 Test to a South African cohort: a platform comparative prospective cohort study

Background

Adverse drug reactions and lack of therapeutic efficacy associated with currently prescribed pharmacotherapeutics may be attributed, in part, to inter-individual variability in drug metabolism. Studies on the pharmacogenetics of Cytochrome P450 (CYP) enzymes offer insight into this variability. The objective of this study was to compare the AmpliChip CYP450 Test^® (AmpliChip) to alternative genotyping platforms for phenotype prediction of CYP2C19 and CYP2D6 in a representative cohort of the South African population.

Methods

AmpliChip was used to screen for thirty-three CYP2D6 and three CYP2C19 alleles in two different cohorts. As a comparison cohort 2 was then genotyped using a CYP2D6 specific long range PCR with sequencing (CYP2D6 XL-PCR + Sequencing) platform and a PCR-RFLP platform for seven CYP2C19 alleles.

Results

Even though there was a low success rate for the AmpliChip, allele frequencies for both CYP2D6 and CYP2C19 were very similar between the two different cohorts. The CYP2D6 XL-PCR + Sequencing platform detected CYP2D6*5 more reliably and could correctly distinguish between CYP2D6*2 and *41 in the Black African individuals. Alleles not covered by the AmpliChip were identified and four novel CYP2D6 alleles were also detected. CYP2C19 PCR-RFLP identified CYP2C19*9,*15, *17 and *27 in the Black African individuals, with *2, *17 and *27 being relatively frequent in the cohort. Eliminating mismatches and identifying additional alleles will contribute to improving phenotype prediction for both enzymes. Phenotype prediction differed between platforms for both genes.

Conclusion

Comprehensive genotyping of CYP2D6 and CYP2C19 with the platforms used in this study, would be more appropriate than AmpliChip for phenotypic prediction in the South African population. Pharmacogenetically important novel alleles may remain undiscovered when using assays that are designed according to Caucasian specific variation, unless alternate strategies are utilised.

Association of MB-COMT polymorphisms with schizophrenia-susceptibility and symptom severity in an African cohort

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.

Extended haplotype studies in South African and Dutch variegate porphyria families carrying the recurrent p.R59W mutation confirm a common ancestry

BACKGROUND

Variegate porphyria (VP) is due to a partial deficiency of protoporphyrinogen oxidase (PPOX), the seventh enzyme in the haem biosynthetic pathway. Clinically, VP is characterized by photosensitivity and acute neurovisceral attacks that can manifest separately or together in affected individuals. The disease is inherited in an autosomal dominant fashion with incomplete penetrance and PPOX gene mutations associated with VP are usually unique to patients and their families. In South Africa, however, VP is highly prevalent as the result of a founder mutation, designated p.R59W. Previous genealogical and haplotype studies showed a link between South African and Dutch carriers of p.R59W and it was suggested that this mutation was introduced to South Africa by Dutch settlers at the end of the 17th century.

OBJECTIVES

To perform extended haplotype analysis in six South African and Dutch VP families with the p.R59W mutation.

METHODS

Haplotyping of 13 microsatellite markers flanking the PPOX gene on chromosome 1q22-23 and five informative single nucleotide polymorphisms within and around the gene.

RESULTS

A core haplotype cosegregated in all families studied.

CONCLUSIONS

Our data deliver further confirmation that the South African and Dutch VP families carrying mutation p.R59W shared a common ancestor.

Whole-genome resequencing in pharmacogenomics: moving away from past disparities to globally representative applications

Africa suffers from a high burden of disease; nonetheless, it has been one of the most under-represented continents with regard to genomic research. It can be argued that this disproportionate research is related to the fact that the genome architecture of African individuals is poorly suited to SNP-based genome-wide association studies, given existing genotyping platforms. However, this argument is no longer plausible with the arrival of next-generation sequencing technologies, which allow for the analysis of entire genomes. Using pharmacogenes to critically examine the merit of next-generation sequencing technologies in pharmacogenomics, we found a substantial amount of novel/uncharacterized variation, which was predicted to alter protein function. This variation was predominantly observed in African individuals, emphasizing the benefit of next-generation sequencing technologies specifically for these individuals. We also observed an improvement in the reliability of sequencing technologies in a relatively short time. Therefore, as sequencing technologies develop and decrease in cost, the ability to reliably detect variation will improve and these technologies will begin to replace other less comprehensive genotyping assays.

Psychiatric genetics in South Africa: cutting a rough diamond

Psychiatric disorders place a considerable healthcare burden on South African society. Incorporating genetic technologies into future treatment plans offers a potential mechanism to reduce this burden. This review focuses on psychiatric genetic research that has been performed in South African populations with regards to obsessive-compulsive disorder, schizophrenia and bipolar disorder. Preliminary findings from these studies suggest that data obtained in developed countries cannot necessarily be extrapolated to South African population groups. Psychiatric genetic studies in South Africa seem to involve relatively low-cost methodologies and only a limited number of large national collaborative studies. Future research in South Africa should therefore aim to incorporate high-throughput technologies into large scale psychiatric studies through the development of collaborations. On a global level, the vast majority of psychiatric genetic studies have been performed in non-African populations. South Africa, as the leading contributor to scientific research in Africa, may provide a foundation for addressing this disparity and strengthening psychiatric genetic research on the continent. Although the elucidation of the genetic architecture of psychiatric disorders has proved challenging, examining the unique genetic profiles found in South African populations could provide valuable insight into the genetics of psychiatric disorders.

Pharmacogenomic Research in South Africa: Lessons Learned and Future Opportunities in the Rainbow Nation

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.

Characterization of the genetic profile of CYP2C19 in two South African populations

AIMS

This study was aimed at elucidating the common sequence variation present in the CYP2C19 gene within the South African Xhosa population and comparing it with the Cape Mixed Ancestry (CMA) population for possible future pharmacogenetic applications.

MATERIALS & METHODS

Common sequence variation was identified through the resequencing of 15 Xhosa individuals. The detected variants were prioritized for genotyping in an additional 85 Xhosa and 75 CMA individuals, while 5 -upstream variants were analyzed using dual luciferase reporter assays.

RESULTS

Resequencing of the Xhosa population revealed 30 variants, including the novel CYP2C19*27 and CYP2C19*28 alleles. CYP2C19*27, characterized by -1041G>A, caused a twofold decrease in luciferase activity, while CYP2C19*28 is characterized by the nonsynonymous V374I variant. In addition, the previously characterized variants, CYP2C19*2, CYP2C19*9 and CYP2C19*17, were present in both populations, while CYP2C19*3 was only observed in the CMA population.

CONCLUSION

Our data demonstrate that both the Xhosa and CMA populations exhibit unique genetic profiles that could influence the outcome of drug therapy in these populations.

Elucidation of CYP2D6 genetic diversity in a unique African population: implications for the future application of pharmacogenetics in the Xhosa population

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.