Risk and Resilience Variants in the Retinoic Acid Metabolic and Developmental Pathways Associated with Risk of FASD Outcomes.

Fetal Alcohol Spectrum Disorder (FASD) is a common neurodevelopmental disorder that affects an estimated 2-5% of North Americans. FASD is induced by prenatal alcohol exposure (PAE) during pregnancy and while there is a clear genetic contribution, few genetic factors are currently identified or understood. In this study, using a candidate gene approach, we performed a genetic variant analysis of retinoic acid (RA) metabolic and developmental signaling pathway genes on whole exome sequencing data of 23 FASD-diagnosed individuals. We found risk and resilience alleles in ADH and ALDH genes known to normally be involved in alcohol detoxification at the expense of RA production, causing RA deficiency, following PAE. Risk and resilience variants were also identified in RA-regulated developmental pathway genes, especially in SHH and WNT pathways. Notably, we also identified significant variants in the causative genes of rare neurodevelopmental disorders sharing comorbidities with FASD, including STRA6 (Matthew-Wood), SOX9 (Campomelic Dysplasia), FDG1 (Aarskog), and 22q11.2 deletion syndrome (TBX1). Although this is a small exploratory study, the findings support PAE-induced RA deficiency as a major etiology underlying FASD and suggest risk and resilience variants may be suitable biomarkers to determine the risk of FASD outcomes following PAE.

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 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.

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.

Association Between HLA-DPB1 and Antineutrophil Cytoplasmic Antibody-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.

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.

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).

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.

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 Systematic Review and Analysis of the Use of Polygenic Scores in Pharmacogenomics.

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.

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.

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.

Secondary biogenic amine deficiencies: genetic etiology, therapeutic interventions, and clinical effects.

Monoamine neurotransmitter disorders present predominantly with neurologic features, including dystonic or dyskinetic cerebral palsy and movement disorders. Genetic conditions that lead to secondary defects in the synthesis, catabolism, transport, and metabolism of biogenic amines can lead to neurotransmitter abnormalities, which can present with similar features. Eleven patients with secondary neurotransmitter abnormalities were enrolled between 2011 and 2015. All patients underwent research-based whole exome and/or whole genome sequencing (WES/WGS). A trial of treatment with levodopa/carbidopa and 5-hydroxytryptophan was initiated. In six families with abnormal neurotransmitter profiles and neurological phenotypes, variants in known disease-causing genes (KCNJ6, SCN2A, CSTB in 2 siblings, NRNX1, KIF1A and PAK3) were identified, while one patient had a variant of uncertain significance in a candidate gene (DLG4) that may explain her phenotype. In 3 patients, no compelling candidate genes were identified. A trial of neurotransmitter replacement therapy led to improvement in motor and behavioral symptoms in all but two patients. The patient with KCNJ6 variant did not respond to L-dopa therapy, but rather experienced increased dyskinetic movements even at low dose of medication. The patient's symptoms harboring the NRNX1 deletion remained unaltered. This study demonstrates the utility of genome-wide sequencing in further understanding the etiology and pathophysiology of neurometabolic conditions, and the potential of secondary neurotransmitter deficiencies to serve as novel therapeutic targets. As there was a largely favorable response to therapy in our case series, a careful trial of neurotransmitter replacement therapy should be considered in patients with cerebrospinal fluid (CSF) monoamines below reference range.

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.

Adult GAMT deficiency: A literature review and report of two siblings.

Guanidinoacetate methyltransferase (GAMT) deficiency is a creatine deficiency disorder and an inborn error of metabolism presenting with progressive intellectual and neurological deterioration. As most cases are identified and treated in early childhood, adult phenotypes that can help in understanding the natural history of the disorder are rare. We describe two adult cases of GAMT deficiency from a consanguineous family in Pakistan that presented with a history of global developmental delay, cognitive impairments, excessive drooling, behavioral abnormalities, contractures and apparent bone deformities initially presumed to be the reason for abnormal gait. Exome sequencing identified a homozygous nonsense variant in GAMT: NM_000156.5:c.134G>A (p.Trp45*). We also performed a literature review and compiled the genetic and clinical characteristics of all adult cases of GAMT deficiency reported to date. When compared to the adult cases previously reported, the musculoskeletal phenotype and the rapidly progressive nature of neurological and motor decline seen in our patients is striking. This study presents an opportunity to gain insights into the adult presentation of GAMT deficiency and highlights the need for in-depth evaluation and reporting of clinical features to expand our understanding of the phenotypic spectrum.

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.

Lay abstract Disease-focused genome-wide association studies (GWAS) have identified a plethora of actionable genetic variants over the last 20 years. International collaboration and technological breakthroughs have enabled rapid genotyping and data collection, which has correspondingly increased sample size and power. Contrastingly, recruitment of well-characterized cohorts of patients for pharmacogenomics research has proven challenging. Given the greater number of associated genetic variants and larger cohort sizes in common disease GWAS, we hypothesized that integrating genome-wide disease and pharmacogenomic data may drive new insights into drug toxicity and drug efficacy phenotypes, beyond the standard scope of current pharmacogenetic analyses. Using GWAS summary statistics from the UK Biobank, European Bioinformatics Institute GWAS catalog, and the Pharmacogenomics Research Network, and a methodological framework incorporating colocalization analysis, we validated pleiotropy at the ABCG2 locus between allopurinol response, gout, and serum urate and identified novel pleiotropy between antihypertensive-induced new-onset diabetes, Crohn's disease and inflammatory bowel disease at the IL18RAP/SLC9A4 locus.

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.

A systematic comparison of pharmacogene star allele calling bioinformatics algorithms: a focus on CYP2D6 genotyping.

Genetic variation in genes encoding cytochrome P450 enzymes has important clinical implications for drug metabolism. Bioinformatics algorithms for genotyping these highly polymorphic genes using high-throughput sequence data and automating phenotype prediction have recently been developed. The CYP2D6 gene is often used as a model during the validation of these algorithms due to its clinical importance, high polymorphism, and structural variations. However, the validation process is often limited to common star alleles due to scarcity of reference datasets. In addition, there has been no comprehensive benchmark of these algorithms to date. We performed a systematic comparison of three star allele calling algorithms using 4618 simulations as well as 75 whole-genome sequence samples from the GeT-RM project. Overall, we found that Aldy and Astrolabe are better suited to call both common and rare diplotypes compared to Stargazer, which is affected by population structure. Aldy was the best performing algorithm in calling CYP2D6 structural variants followed by Stargazer, whereas Astrolabe had limitations especially in calling hybrid rearrangements. We found that ensemble genotyping, characterised by taking a consensus of genotypes called by all three algorithms, has higher haplotype concordance but it is prone to ambiguities whenever complete discrepancies between the tools arise. Further, we evaluated the effects of sequencing coverage and indel misalignment on genotyping accuracy. Our account of the strengths and limitations of these algorithms is extremely important to clinicians and researchers in the pharmacogenomics and precision medicine communities looking to haplotype CYP2D6 and other pharmacogenes using high-throughput sequencing data.

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.