SLCO1B1 Genetic Variation Influence on Atorvastatin Systemic Exposure in Pediatric Hypercholesterolemia.

This clinical study examined the influence of SLCO1B1 c.521T>C (rs4149056) on plasma atorvastatin concentrations in pediatric hypercholesterolemia. The participants (8-21 years), including heterozygous (c.521T/C, n = 13), homozygous (c.521C/C, n = 2) and controls (c.521T/T, n = 13), completed a single-oral-dose pharmacokinetic study. Similar to in adults, the atorvastatin (AVA) area-under-concentration-time curve from 0 to 24 h (AUC0-24) was 1.7-fold and 2.8-fold higher in participants with c.521T/C and c.521C/C compared to the c.521T/T participants, respectively. The inter-individual variability in AVA exposure within these genotype groups ranged from 2.3 to 4.8-fold, indicating that additional factors contribute to the inter-individual variability in the AVA dose-exposure relationship. A multivariate model reinforced the SLCO1B1 c.521T>C variant as the central factor contributing to AVA systemic exposure in this pediatric cohort, accounting for ~65% of the variability in AVA AUC0-24. Furthermore, lower AVA lactone concentrations in participants with increased body mass index contributed to higher exposure within the c.521T/T and c.521T/C genotype groups. Collectively, these factors contributing to higher systemic exposure could increase the risk of toxicity and should be accounted for when individualizing the dosing of atorvastatin in eligible pediatric patients.

Updated DPYD HapB3 haplotype structure and implications for pharmacogenomic testing.

The DPYD gene encodes dihydropyrimidine dehydrogenase, the rate-limiting enzyme for the metabolism of fluoropyrimidines 5-fluorouracil and capecitabine. Genetic variants in DPYD have been associated with altered enzyme activity, therefore accurate detection and interpretation is critical to predict metabolizer status for individualized fluoropyrimidine therapy. The most commonly observed deleterious variation is the causal variant linked to the previously described HapB3 haplotype, c.1129-5923C>G (rs75017182) in intron 10, which introduces a cryptic splice site. A benign synonymous variant in exon 11, c.1236G>A (rs56038477) is also linked to HapB3 and is commonly used for testing. Previously, these single-nucleotide polymorphisms (SNPs) have been reported to be in perfect linkage disequilibrium (LD); therefore, c.1236G>A is often utilized as a proxy for the function-altering intronic variant. Clinical genotyping of DPYD identified a patient who had c.1236G>A, but not c.1129-5923C>G, suggesting that these two SNPs may not be in perfect LD, as previously assumed. Additional individuals with c.1236G>A, but not c.1129-5923C>G, were identified in the Children's Mercy Data Warehouse and the All of Us Research Program version 7 cohort substantiating incomplete SNP linkage. Consequently, testing only c.1236G>A can generate false-positive results in some cases and lead to suboptimal dosing that may negatively impact patient therapy and prospect of survival. Our data show that DPYD genotyping should include the functional variant c.1129-5923C>G, and not the c.1236G>A proxy, to accurately predict DPD activity.

PharmVar GeneFocus: SLCO1B1.

The Pharmacogene Variation Consortium (PharmVar) is now providing star (*) allele nomenclature for the highly polymorphic human SLCO1B1 gene encoding the organic anion transporting polypeptide 1B1 (OATP1B1) drug transporter. Genetic variation within the SLCO1B1 gene locus impacts drug transport, which can lead to altered pharmacokinetic profiles of several commonly prescribed drugs. Variable OATP1B1 function is of particular importance regarding hepatic uptake of statins and the risk of statin-associated musculoskeletal symptoms. To introduce this important drug transporter gene into the PharmVar database and serve as a unified reference of haplotype variation moving forward, an international group of gene experts has performed an extensive review of all published SLCO1B1 star alleles. Previously published star alleles were self-assigned by authors and only loosely followed the star nomenclature system that was first developed for cytochrome P450 genes. This nomenclature system has been standardized by PharmVar and is now applied to other important pharmacogenes such as SLCO1B1. In addition, data from the 1000 Genomes Project and investigator-submitted data were utilized to confirm existing haplotypes, fill knowledge gaps, and/or define novel star alleles. The PharmVar-developed SLCO1B1 nomenclature has been incorporated by the Clinical Pharmacogenetics Implementation Consortium (CPIC) 2022 guideline on statin-associated musculoskeletal symptoms.

Impact of CYP2C:TG Haplotype on CYP2C19 Substrates Clearance In Vivo, Protein Content, and In Vitro Activity.

A novel haplotype composed of two non-coding variants, CYP2C18 NM_000772.3:c.*31T (rs2860840) and NM_000772.2:c.819+2182G (rs11188059), referred to as "CYP2C:TG," was recently associated with ultrarapid metabolism of various CYP2C19 substrates. As the underlying mechanism and clinical relevance of this effect remain uncertain, we analyzed existing in vivo and in vitro data to determine the magnitude of the CYP2C:TG haplotype effect. We assessed variability in pharmacokinetics of CYP2C19 substrates, including citalopram, sertraline, voriconazole, omeprazole, pantoprazole, and rabeprazole in 222 healthy volunteers receiving one of these six drugs. We also determined its impact on CYP2C8, CYP2C9, CYP2C18, and CYP2C19 protein abundance in 135 human liver tissue samples, and on CYP2C18/CYP2C19 activity in vitro using N-desmethyl atomoxetine formation. No effects were observed according to CYP2C:TG haplotype or to CYP2C19*1+TG alleles (i.e., CYP2C19 alleles containing the CYP2C:TG haplotype). In contrast, CYP2C19 intermediate (e.g., CYP2C19*1/*2) and poor metabolizers (e.g., CYP2C19*2/*2) showed significantly higher exposure in vivo, lower CYP2C19 protein abundance in human liver microsomes, and lower activity in vitro compared with normal, rapid (i.e., CYP2C19*1/*17), and ultrarapid metabolizers (i.e., CYP2C19*17/*17). Moreover, a tendency toward lower exposure was observed in ultrarapid metabolizers compared with rapid metabolizers and normal metabolizers. Furthermore, when the CYP2C19*17 allele was present, CYP2C18 protein abundance was increased suggesting that genetic variation in CYP2C19 may be relevant to the overall metabolism of certain drugs by regulating not only its expression levels, but also those of CYP2C18. Considering all available data, we conclude that there is insufficient evidence supporting clinical CYP2C:TG testing to inform drug therapy.

A DEAD-box RNA helicase mediates meiotic silencing by unpaired DNA.

During the sexual phase of Neurospora crassa, unpaired genes are subject to a silencing mechanism known as meiotic silencing by unpaired DNA (MSUD). MSUD targets the transcripts of an unpaired gene and utilizes typical RNA interference factors for its process. Using a reverse genetic screen, we have identified a meiotic silencing gene called sad-9, which encodes a DEAD-box RNA helicase. While not essential for vegetative growth, SAD-9 plays a crucial role in both sexual development and MSUD. Our results suggest that SAD-9, with the help of the SAD-2 scaffold protein, recruits the SMS-2 Argonaute to the perinuclear region, the center of MSUD activity.

Characterization of Reference Materials for CYP3A4 and CYP3A5: A (GeT-RM) Collaborative Project.

Pharmacogenetic testing for CYP3A4 is increasingly provided by clinical and research laboratories; however, only a limited number of quality control and reference materials are currently available for many of the CYP3A4 variants included in clinical tests. To address this need, the Division of Laboratory Systems, CDC-based Genetic Testing Reference Material Coordination Program (GeT-RM), in collaboration with members of the pharmacogenetic testing and research communities and the Coriell Institute for Medical Research, has characterized 30 DNA samples derived from Coriell cell lines for CYP3A4. Samples were distributed to five volunteer laboratories for genotyping using a variety of commercially available and laboratory-developed tests. Sanger and next-generation sequencing were also utilized by some of the laboratories. Whole-genome sequencing data from the 1000 Genomes Projects were utilized to inform genotype. Twenty CYP3A4 alleles were identified in the 30 samples characterized for CYP3A4: CYP3A4∗4, ∗5, ∗6, ∗7, ∗8, ∗9, ∗10, ∗11, ∗12, ∗15, ∗16, ∗18, ∗19, ∗20, ∗21, ∗22, ∗23, ∗24, ∗35, and a novel allele, CYP3A4∗38. Nineteen additional samples with preexisting data for CYP3A4 or CYP3A5 were re-analyzed to generate comprehensive reference material panels for these genes. These publicly available and well-characterized materials can be used to support the quality assurance and quality control programs of clinical laboratories performing clinical pharmacogenetic testing.

The Identification of Novel CYP2D6 Variants in US Hmong: Results From Genome Sequencing and Clinical Genotyping.

Objective: Hmong individuals represent a unique East Asian subpopulation in whom limited information concerning pharmacogenetic variation exists. The objectives of this study were to comprehensively characterize the highly polymorphic CYP2D6 gene in Hmong, estimate allele and phenotype frequencies and to compare results between two testing platforms. Methods: DNA from 48 self-identified Hmong participants were sequenced using a targeted next-generation sequencing (NGS) panel. Star allele calls were made using Astrolabe, manual inspection of NGS variant calls and confirmatory Sanger sequencing. Structural variation was determined by long-range (XL)-PCR and digital droplet PCR (ddPCR). The consensus diplotypes were subsequently translated into phenotype utilizing the activity score system. Clinical grade pharmacogenetic testing was obtained for 12 of the 48 samples enabling an assessment of concordance between the consensus calls and those determined by clinical testing platforms. Results: A total of 13 CYP2D6 alleles were identified. The most common alleles were CYP2D6*10 and its structural arrangements (37.5%, 36/96) and the *5 gene deletion (13.5%, 13/96). Three novel suballeles (*10.007, *36.004, and *75.002) were also identified. Phenotype frequencies were as follows: ultrarapid metabolizers (4.2%, 2/48), normal metabolizers (41.7%, 20/48) and intermediate metabolizers (52.1%, 25/48); none of the 48 participants were predicted to be poor metabolizers. Concordance of diplotype and phenotype calls between the consensus and clinical testing were 66.7 and 50%, respectively. Conclusion: Our study to explore CYP2D6 genotypes in the Hmong population suggests that this subpopulation is unique regarding CYP2D6 allelic variants; also, a higher portion of Hmong participants (50%) are predicted to have an intermediate metabolizer phenotype for CYP2D6 compared to other East Asians which range between 27 and 44%. Results from different testing methods varied considerably. These preliminary findings underscore the importance of thoroughly interrogating unique subpopulations to accurately predict a patient's CYP2D6 metabolizer status.

The Impact of the CYP2D6 "Enhancer" Single Nucleotide Polymorphism on CYP2D6 Activity.

rs5758550 has been associated with enhanced transcription and suggested to be a useful marker of CYP2D6 activity. As there are limited and inconsistent data regarding the utility of this distant "enhancer" single nucleotide polymorphism (SNP), our goal was to further assess the impact of rs5758550 on CYP2D6 activity toward two probe substrates, atomoxetine (ATX) and dextromethorphan (DM), using in vivo urinary metabolite (DM; n = 188) and pharmacokinetic (ATX; n = 70) and in vitro metabolite formation (ATX and DM; n = 166) data. All subjects and tissues were extensively genotyped, the "enhancer" SNP phased with established CYP2D6 haplotypes either computationally or experimentally, and the impact on CYP2D6 activity investigated using several linear models of varying complexity to determine the proportion of variability in CYP2D6 activity captured by each model. For all datasets and models, the "enhancer" SNP had no or only a modest impact on CYP2D6 activity prediction. An increased effect, when present, was more pronounced for ATX than DM suggesting potential substate-dependency. In addition, CYP2D6*2 alleles with the "enhancer" SNP were associated with modestly higher metabolite formation rates in vitro, but not in vivo; no effect was detected for CYP2D6*1 alleles with "enhancer" SNP. In summary, it remains inconclusive whether the small effects detected in this investigation are indeed caused by the "enhancer" SNP or are rather due to its incomplete linkage with other variants within the gene. Taken together, there does not appear to be sufficient evidence to warrant the "enhancer" SNP be included in clinical CYP2D6 pharmacogenetic testing.

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

CYP2C8, CYP2C9, and CYP2C19 Characterization Using Next-Generation Sequencing and Haplotype Analysis: A GeT-RM Collaborative Project.

Pharmacogenetic tests typically target selected sequence variants to identify haplotypes that are often defined by star (∗) allele nomenclature. Due to their design, these targeted genotyping assays are unable to detect novel variants that may change the function of the gene product and thereby affect phenotype prediction and patient care. In the current study, 137 DNA samples that were previously characterized by the Genetic Testing Reference Material (GeT-RM) program using a variety of targeted genotyping methods were recharacterized using targeted and whole genome sequencing analysis. Sequence data were analyzed using three genotype calling tools to identify star allele diplotypes for CYP2C8, CYP2C9, and CYP2C19. The genotype calls from next-generation sequencing (NGS) correlated well to those previously reported, except when novel alleles were present in a sample. Six novel alleles and 38 novel suballeles were identified in the three genes due to identification of variants not covered by targeted genotyping assays. In addition, several ambiguous genotype calls from a previous study were resolved using the NGS and/or long-read NGS data. Diplotype calls were mostly consistent between the calling algorithms, although several discrepancies were noted. This study highlights the utility of NGS for pharmacogenetic testing and demonstrates that there are many novel alleles that are yet to be discovered, even in highly characterized genes such as CYP2C9 and CYP2C19.

Characterization of Reference Materials for TPMT and NUDT15: A GeT-RM Collaborative Project.

Pharmacogenetic testing is increasingly provided by clinical and research laboratories; however, only a limited number of quality control and reference materials are currently available for many of the TPMT and NUDT15 variants included in clinical tests. To address this need, the Division of Laboratory Systems, Centers for Disease Control and Prevention-based Genetic Testing Reference Material (GeT-RM) coordination program, in collaboration with members of the pharmacogenetic testing and research communities and the Coriell Institute for Medical Research, has characterized 19 DNA samples derived from Coriell cell lines. DNA samples were distributed to four volunteer testing laboratories for genotyping using a variety of commercially available and laboratory developed tests and/or Sanger sequencing. Of the 12 samples characterized for TPMT, newly identified variants include TPMT∗2, ∗6, ∗12, ∗16, ∗21, ∗24, ∗32, ∗33, and ∗40; for the 7 NUDT15 reference material samples, newly identified variants are NUDT15∗2, ∗3, ∗4, ∗5, ∗6, and ∗9. In addition, a novel haplotype, TPMT∗46, was identified in this study. Preexisting data on an additional 11 Coriell samples, as well as some supplemental testing, were used to create comprehensive reference material panels for TPMT and NUDT15. These publicly available and well-characterized materials can be used to support the quality assurance and quality control programs of clinical laboratories performing clinical pharmacogenetic testing.

Fluorescent and bimolecular-fluorescent protein tagging of genes at their native loci in Neurospora crassa using specialized double-joint PCR plasmids.

The double-joint polymerase chain reaction (DJ-PCR) is a technique that can be used to construct vectors for targeted genome integration without laborious subcloning steps. Here we report the availability of plasmids that facilitate DJ-PCR-based construction of Neurospora crassa tagging vectors. These plasmids allow the creation of green or red fluorescent protein (GFP or RFP) tagging vectors for protein localization studies, as well as split-yellow fluorescent protein (YFP) tagging vectors for bimolecular fluorescence complementation (BiFC) analyses. We have demonstrated the utility of each plasmid with the tagging of known meiotic silencing proteins. Microscopic analysis of the tagged strains indicates that SMS-2 and QIP form macromolecular complexes in the perinuclear region during meiosis.

Identification of CYP2D6 Haplotypes that Interfere with Commonly Used Assays for Copy Number Variation Characterization.

Cytochrome P450 2D6 (CYP2D6) copy number (CN) variation affects the metabolism of numerous prescribed drugs. Sequence variation within primer or probe target regions of hydrolysis probe CN assays can generate false-positive calls for CN loss. Furthermore, CYP2D6-CYP2D7 hybrids and gene conversions can cause difficult to interpret discordant CN calls. The identification of haplotypes with CN variations and structural arrangements is important to predict phenotype accurately. During clinical testing with hydrolysis probe assays targeting three CYP2D6 regions (intron 2, intron 6, and exon 9), samples with haplotypes causing inconsistent CN calls were identified. To resolve these cases, next-generation sequencing and allele-specific Sanger sequencing was performed. Sequence analysis of 16 samples, all but one from subjects of African descent, identified six novel suballeles containing single-nucleotide polymorphisms, which cause false-positive calls for CN loss in introns 2 and 6. Five samples with an exon 9 CN loss contained CYP2D6/CYP2D7 hybrids (∗13 or ∗36) and one sample was found to have a novel haplotype, CYP2D6∗141. Interestingly, CYP2D6∗141 contains a CYP2D7-derived exon 9 conversion and core single-nucleotide polymorphisms that are otherwise found in CYP2D6∗17 and ∗27. Although these variants are rare, they can cause inconsistent CN calls that typically are reported as no calls or indeterminant, and thus may deprive patients, particularly those of African descent, from taking full benefit of pharmacogenetic testing.

Involvement of RNA granule proteins in meiotic silencing by unpaired DNA.

In Neurospora crassa, expression from an unpaired gene is suppressed by a mechanism known as meiotic silencing by unpaired DNA (MSUD). MSUD utilizes common RNA interference (RNAi) factors to silence target mRNAs. Here, we report that Neurospora CAR-1 and CGH-1, homologs of two Caenorhabditis elegans RNA granule components, are involved in MSUD. These fungal proteins are found in the perinuclear region and P-bodies, much like their worm counterparts. They interact with components of the meiotic silencing complex (MSC), including the SMS-2 Argonaute. This is the first time MSUD has been linked to RNA granule proteins.

Long-Distance Phasing of a Tentative "Enhancer" Single-Nucleotide Polymorphism With CYP2D6 Star Allele Definitions.

BACKGROUND: The CYP2D6 gene locus has been extensively studied over decades, yet a portion of variability in CYP2D6 activity cannot be explained by known sequence variations within the gene, copy number variation, or structural rearrangements. It was proposed that rs5758550, located 116 kb downstream of the CYP2D6 gene locus, increases gene expression and thus contributes to variability in CYP2D6 activity. This finding has, however, not been validated. The purpose of the study was to address a major technological barrier, i.e., experimentally linking rs5758550, also referred to as the "enhancer" single-nucleotide polymorphism (SNP), to CYP2D6 haplotypes >100 kb away. To overcome this challenge is essential to ultimately determine the contribution of the "enhancer" SNP to interindividual variability in CYP2D6 activity. METHODS: A large ethnically mixed population sample (n=3,162) was computationally phased to determine linkage between the "enhancer" SNP and CYP2D6 haplotypes (or star alleles). To experimentally validate predicted linkages, DropPhase2D6, a digital droplet PCR (ddPCR)-based method was developed. 10X Genomics Linked-Reads were utilized as a proof of concept. RESULTS: Phasing predicted that the "enhancer" SNP can occur on numerous CYP2D6 haplotypes including CYP2D6*1, *2, *5, and *41 and suggested that linkage is incomplete, i.e., a portion of these alleles do not have the "enhancer" SNP. Phasing also revealed differences among the European and African ancestry data sets regarding the proportion of alleles with and without the "enhancer" SNP. DropPhase2D6 was utilized to confirm or refute the predicted "enhancer" SNP location for individual samples, e.g., of n=3 samples genotyped as *1/*41, rs5758550 was on the *41 allele of two samples and on the *1 allele of one sample. Our findings highlight that the location of the "enhancer" SNP must not be assigned by "default." Furthermore, linkage between the "enhancer" SNP and CYP2D6 star allele haplotypes was confirmed with 10X Genomics technology. CONCLUSIONS: Since the "enhancer" SNP can be present on a portion of normal, decreased, or no function alleles, the phase of the "enhancer" SNP must be considered when investigating the impact of the "enhancer" SNP on CYP2D6 activity.

An NCBP3-Domain Protein Mediates Meiotic Silencing by Unpaired DNA.

In the filamentous fungus Neurospora crassa, genes unpaired during meiosis are silenced by a process known as meiotic silencing by unpaired DNA (MSUD). MSUD utilizes common RNA interference (RNAi) proteins, such as Dicer and Argonaute, to target homologous mRNAs for silencing. Previously, we demonstrated that nuclear cap-binding proteins NCBP1 and NCBP2 are involved in MSUD. We report here that SAD-8, a protein similar to human NCBP3, also mediates silencing. Although SAD-8 is not essential for either vegetative or sexual development, it is required for MSUD. SAD-8 localizes predominantly in the nucleus and interacts with both NCBP1 and NCBP2. Similar to NCBP1 and NCBP2, SAD-8 interacts with a component (Argonaute) of the perinuclear meiotic silencing complex (MSC), further implicating the involvement of cap-binding proteins in silencing.

Characterization of Reference Materials for Genetic Testing of CYP2D6 Alleles: A GeT-RM Collaborative Project.

Pharmacogenetic testing increasingly is available from clinical and research laboratories. However, only a limited number of quality control and other reference materials currently are available for the complex rearrangements and rare variants that occur in the CYP2D6 gene. To address this need, the Division of Laboratory Systems, CDC-based Genetic Testing Reference Material Coordination Program, in collaboration with members of the pharmacogenetic testing and research communities and the Coriell Cell Repositories (Camden, NJ), has characterized 179 DNA samples derived from Coriell cell lines. Testing included the recharacterization of 137 genomic DNAs that were genotyped in previous Genetic Testing Reference Material Coordination Program studies and 42 additional samples that had not been characterized previously. DNA samples were distributed to volunteer testing laboratories for genotyping using a variety of commercially available and laboratory-developed tests. These publicly available samples will support the quality-assurance and quality-control programs of clinical laboratories performing CYP2D6 testing.

The Nuclear Cap-Binding Complex Mediates Meiotic Silencing by Unpaired DNA.

In the filamentous fungus Neurospora crassa, cross walls between individual cells are normally incomplete, making the entire fungal network vulnerable to attack by viruses and selfish DNAs. Accordingly, several genome surveillance mechanisms are maintained to help the fungus combat these repetitive elements. One of these defense mechanisms is called meiotic silencing by unpaired DNA (MSUD), which identifies and silences unpaired genes during meiosis. Utilizing common RNA interference (RNAi) proteins, such as Dicer and Argonaute, MSUD targets mRNAs homologous to the unpaired sequence to achieve silencing. In this study, we have identified an additional silencing component, namely the cap-binding complex (CBC). Made up of cap-binding proteins CBP20 and CBP80, CBC associates with the 5' cap of mRNA transcripts in eukaryotes. The loss of CBC leads to a deficiency in MSUD activity, suggesting its role in mediating silencing. As confirmed in this study, CBC is predominantly nuclear, although it is known to travel in and out of the nucleus to facilitate RNA transport. As seen in animals but not in plants, CBP20's robust nuclear import depends on CBP80 in Neurospora CBC interacts with a component (Argonaute) of the perinuclear meiotic silencing complex (MSC), directly linking the two cellular factors.

Complex formation of RNA silencing proteins in the perinuclear region of Neurospora crassa.

In Neurospora, genes not paired during meiosis are targeted by meiotic silencing by unpaired DNA (MSUD). Here, our bimolecular fluorescence complementation (BiFC) study suggests that RNA-directed RNA polymerase, Dicer, Argonaute, and others form a silencing complex in the perinuclear region, with intimate interactions among the majority of them. We have also shown that SAD-2 is likely the anchor for this assembly.