Modeling.

Modeling is coming to the fore as it is now widely accepted and indeed expected during drug discovery and development. Modeling integrates knowledge, increases understanding and provides the ability to predict an outcome either before it occurs or when it is not possible to measure. This makes modeling an attractive option for inhaled drugs as it is not possible to routinely measure what is occurring to the drug (pharmacokinetics) and what effect the drug is having (pharmacodynamics) at local microscopic sites of such a diverse and complex organ as the lung. Many pieces of information (data and knowledge) exist like the pieces of a jigsaw puzzle and modeling brings the pieces together in a scientific and mechanistically coherent manner to increase understanding of both the efficacy and safety of inhaled drugs.

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.

Physiologically Based Pharmacokinetic Modelling of Inhaled Nemiralisib: Mechanistic Components for Pulmonary Absorption, Systemic Distribution, and Oral Absorption.

BACKGROUND AND OBJECTIVES: Physiologically based pharmacokinetic (PBPK) modelling has evolved to accommodate different routes of drug administration and enables prediction of drug concentrations in tissues as well as plasma. The inhalation route of administration has proven successful in treating respiratory diseases but can also be used for rapid systemic delivery, holding great promise for treatment of diseases requiring systemic exposure. The objective of this work was to develop a PBPK model that predicts plasma and tissue concentrations following inhalation administration of the PI3Kδ inhibitor nemiralisib. METHODS: A PBPK model was built in GastroPlus® that includes a complete mechanistic description of pulmonary absorption, systemic distribution and oral absorption following inhalation administration of nemiralisib. The availability of clinical data obtained after intravenous, oral and inhalation administration enabled validation of the model with observed data and accurate assessment of pulmonary drug absorption. The PBPK model described in this study incorporates novel use of key parameters such as lung systemic absorption rate constants derived from human physiological lung blood flows, and implementation of the specific permeability-surface area product per millilitre of tissue cell volume (SpecPStc) to predict tissue distribution. RESULTS: The inhaled PBPK model was verified using plasma and bronchoalveolar lavage fluid concentration data obtained in human subjects. Prediction of tissue concentrations using the permeability-limited systemic disposition tissue model was further validated using tissue concentration data obtained in the rat following intravenous infusion administration to steady state. CONCLUSIONS: Fully mechanistic inhaled PBPK models such as the model described herein could be applied for cross molecule assessments with respect to lung retention and systemic exposure, both in terms of pharmacology and toxicology, and may facilitate clinical indication selection.

PharmVar GeneFocus: CYP2B6.

The Pharmacogene Variation Consortium (PharmVar) catalogs star (*) allele nomenclature for the polymorphic human CYP2B6 gene. Genetic variation within the CYP2B6 gene locus impacts the metabolism or bioactivation of clinically important drugs. Of particular importance are efficacy and safety concerns regarding: efavirenz, which is used for the treatment of HIV type-1 infection; methadone, a mainstay in the treatment of opioid use disorder and as an analgesic; ketamine, used as an antidepressant and analgesic; and bupropion, which is prescribed to treat depression and for smoking cessation. This GeneFocus provides a comprehensive overview and summary of CYP2B6 and describes how haplotype information catalogued by PharmVar is utilized by the Pharmacogenomics Knowledgebase (PharmGKB) and the Clinical Pharmacogenetics Implementation Consortium (CPIC).

ITGB2 (Integrin ß2) Immunomodulatory Gene Variants in Premature Infants With Necrotizing Enterocolitis.

ABSTRACT: Aberrant toll-like receptor (TLR) activation is central to necrotizing enterocolitis (NEC) pathogenesis. ß2 integrins regulate TLR signaling, and integrin ß2 (ITGB2) deficiency causes TLR hyperresponsiveness. To test the hypothesis that ITGB2 genetic variants modulate NEC susceptibility, we sequenced the exonic ITGB2 locus to compare the prevalence of deleterious variants among 221 preterm infants with and without NEC. ITGB2 variants were not associated with NEC in our entire cohort (NEC [9/56] versus controls [16/165], P = 0.19) or in extremely low birthweight infants (ELBW, controls [7.9%] versus NEC [18.2%]; P = 0.11) but were increased compared to the populace (4.5%, gnomad.broadinstitute.org). Combined annotation-dependent depletion -predicted deleterious ITGB2 variants increased proportionately with increasing NEC severity in ELBW infants (controls [6.7%] versus medical NEC [16.7%] versus surgical NEC [19%] (P = 0.03). Although ITGB2 variants were not associated with NEC in our preterm cohort, subgroup analysis showed a trend towards enrichment with NEC severity in ELBW infants.

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.

Physiologically Based Pharmacokinetic Modelling for First-In-Human Predictions: An Updated Model Building Strategy Illustrated with Challenging Industry Case Studies.

Physiologically based pharmacokinetic modelling is well established in the pharmaceutical industry and is accepted by regulatory agencies for the prediction of drug-drug interactions. However, physiologically based pharmacokinetic modelling is valuable to address a much wider range of pharmaceutical applications, and new regulatory impact is expected as its full power is leveraged. As one example, physiologically based pharmacokinetic modelling is already routinely used during drug discovery for in-vitro to in-vivo translation and pharmacokinetic modelling in preclinical species, and this leads to the application of verified models for first-in-human pharmacokinetic predictions. A consistent cross-industry strategy in this application area would increase confidence in the approach and facilitate further learning. With this in mind, this article aims to enhance a previously published first-in-human physiologically based pharmacokinetic model-building strategy. Based on the experience of scientists from multiple companies participating in the GastroPlus™ User Group Steering Committee, new Absorption, Distribution, Metabolism and Excretion knowledge is integrated and decision trees proposed for each essential component of a first-in-human prediction. We have reviewed many relevant scientific publications to identify new findings and highlight gaps that need to be addressed. Finally, four industry case studies for more challenging compounds illustrate and highlight key components of the strategy.

Correction: Whole exome sequencing reveals HSPA1L as a genetic risk factor for spontaneous preterm birth.

[This corrects the article DOI: 10.1371/journal.pgen.1007394.].

Whole exome sequencing reveals HSPA1L as a genetic risk factor for spontaneous preterm birth.

Preterm birth is a leading cause of morbidity and mortality in infants. Genetic and environmental factors play a role in the susceptibility to preterm birth, but despite many investigations, the genetic basis for preterm birth remain largely unknown. Our objective was to identify rare, possibly damaging, nucleotide variants in mothers from families with recurrent spontaneous preterm births (SPTB). DNA samples from 17 Finnish mothers who delivered at least one infant preterm were subjected to whole exome sequencing. All mothers were of northern Finnish origin and were from seven multiplex families. Additional replication samples of European origin consisted of 93 Danish sister pairs (and two sister triads), all with a history of a preterm delivery. Rare exonic variants (frequency <1%) were analyzed to identify genes and pathways likely to affect SPTB susceptibility. We identified rare, possibly damaging, variants in genes that were common to multiple affected individuals. The glucocorticoid receptor signaling pathway was the most significant (p<1.7e-8) with genes containing these variants in a subgroup of ten Finnish mothers, each having had 2-4 SPTBs. This pathway was replicated among the Danish sister pairs. A gene in this pathway, heat shock protein family A (Hsp70) member 1 like (HSPA1L), contains two likely damaging missense alleles that were found in four different Finnish families. One of the variants (rs34620296) had a higher frequency in cases compared to controls (0.0025 vs. 0.0010, p = 0.002) in a large preterm birth genome-wide association study (GWAS) consisting of mothers of general European ancestry. Sister pairs in replication samples also shared rare, likely damaging HSPA1L variants. Furthermore, in silico analysis predicted an additional phosphorylation site generated by rs34620296 that could potentially affect chaperone activity or HSPA1L protein stability. Finally, in vitro functional experiment showed a link between HSPA1L activity and decidualization. In conclusion, rare, likely damaging, variants in HSPA1L were observed in multiple families with recurrent SPTB.

The NSIGHT1-randomized controlled trial: rapid whole-genome sequencing for accelerated etiologic diagnosis in critically ill infants.

Genetic disorders are a leading cause of morbidity and mortality in infants in neonatal and pediatric intensive care units (NICU/PICU). While genomic sequencing is useful for genetic disease diagnosis, results are usually reported too late to guide inpatient management. We performed an investigator-initiated, partially blinded, pragmatic, randomized, controlled trial to test the hypothesis that rapid whole-genome sequencing (rWGS) increased the proportion of NICU/PICU infants receiving a genetic diagnosis within 28 days. The participants were families with infants aged <4 months in a regional NICU and PICU, with illnesses of unknown etiology. The intervention was trio rWGS. Enrollment from October 2014 to June 2016, and follow-up until November 2016. Of all, 26 female infants, 37 male infants, and 2 infants of undetermined sex were randomized to receive rWGS plus standard genetic tests (n = 32, cases) or standard genetic tests alone (n = 33, controls). The study was terminated early due to loss of equipoise: 73% (24) controls received genomic sequencing as standard tests, and 15% (five) controls underwent compassionate cross-over to receive rWGS. Nevertheless, intention to treat analysis showed the rate of genetic diagnosis within 28 days of enrollment (the primary end-point) to be higher in cases (31%, 10 of 32) than controls (3%, 1 of 33; difference, 28% [95% CI, 10-46%]; p = 0.003). Among infants enrolled in the first 25 days of life, the rate of neonatal diagnosis was higher in cases (32%, 7 of 22) than controls (0%, 0 of 23; difference, 32% [95% CI, 11-53%];p = 0.004). Median age at diagnosis (25 days [range 14-90] in cases vs. 130 days [range 37-451] in controls) and median time to diagnosis (13 days [range 1-84] in cases, vs. 107 days [range 21-429] in controls) were significantly less in cases than controls (p = 0.04). In conclusion, rWGS increased the proportion of NICU/PICU infants who received timely diagnoses of genetic diseases.

The Pharmacogene Variation (PharmVar) Consortium: Incorporation of the Human Cytochrome P450 (CYP) Allele Nomenclature Database.

The Human Cytochrome P450 (CYP) Allele Nomenclature Database, a critical resource for the pharmacogenetics and genomics communities, has transitioned to the Pharmacogene Variation (PharmVar) Consortium. In this report we provide a summary of the current database, provide an overview of the PharmVar consortium, and highlight the PharmVar database which will serve as the new home for pharmacogene nomenclature.

Erratum: Constellation: a tool for rapid, automated phenotype assignment of a highly polymorphic pharmacogene, CYP2D6, from whole-genome sequences.

[This corrects the article DOI: 10.1038/npjgenmed.2015.7.].

In vivo characterization of CYP2D6*12, *29 and *84 using dextromethorphan as a probe drug: a case report.

CYP2D6*84 was first described in a Black South African subject, however, its function remains unknown. Astrolabe, a probabilistic scoring tool developed in our laboratory to call genotypes from whole genome sequence, identified CYP2D6*84 in a trio. The father presented with intermediate metabolism when challenged with the CYP2D6 probe drug dextromethorphan (DM/dextrorphan [DX] = 0.0839). Since his second allele, CYP2D6*12, is nonfunctional, the observed activity is derived by CYP2D6*84. This finding suggests that the allele's hallmark P267H causes decreased activity toward DM and that this allele should receive a value of 0.5 for Activity Score calculations. The mother's DM/DX of 0.0543 was consistent with the decreased activity classification of CYP2D6*29. The child, a critically ill neonate, was not phenotyped, but predicted to be a normal metabolizer.

Constellation: a tool for rapid, automated phenotype assignment of a highly polymorphic pharmacogene, CYP2D6, from whole-genome sequences.

An important component of precision medicine-the use of whole-genome sequencing (WGS) to guide lifelong healthcare-is electronic decision support to inform drug choice and dosing. To achieve this, automated identification of genetic variation in genes involved in drug absorption, distribution, metabolism, excretion and response (ADMER) is required. CYP2D6 is a major enzyme for drug bioactivation and elimination. CYP2D6 activity is predominantly governed by genetic variation; however, it is technically arduous to haplotype. Not only is the nucleotide sequence of CYP2D6 highly polymorphic, but the locus also features diverse structural variations, including gene deletion, duplication, multiplication events and rearrangements with the nonfunctional, neighbouring CYP2D7 and CYP2D8 genes. We developed Constellation, a probabilistic scoring system, enabling automated ascertainment of CYP2D6 activity scores from 2×100 paired-end WGS. The consensus reference method included TaqMan genotyping assays, quantitative copy-number variation determination and Sanger sequencing. When compared with the consensus reference Constellation had an analytic sensitivity of 97% (59 of 61 diplotypes) and analytic specificity of 95% (116 of 122 haplotypes). All extreme phenotypes, i.e., poor and ultrarapid metabolisers were accurately identified by Constellation. Constellation is anticipated to be extensible to functional variation in all ADMER genes, and to be performed at marginal incremental financial and computational costs in the setting of diagnostic WGS.

Clinical detection of deletion structural variants in whole-genome sequences.

Optimal management of acutely ill infants with monogenetic diseases requires rapid identification of causative haplotypes. Whole-genome sequencing (WGS) has been shown to identify pathogenic nucleotide variants in such infants. Deletion structural variants (DSVs, >50 nt) are implicated in many genetic diseases, and tools have been designed to identify DSVs using short-read WGS. Optimisation and integration of these tools into a WGS pipeline could improve diagnostic sensitivity and specificity of WGS. In addition, it may improve turnaround time when compared with current CNV assays, enhancing utility in acute settings. Here we describe DSV detection methods for use in WGS for rapid diagnosis in acutely ill infants: SKALD (Screening Konsensus and Annotation of Large Deletions) combines calls from two tools (Breakdancer and GenomeStrip) with calibrated filters and clinical interpretation rules. In four WGS runs, the average analytic precision (positive predictive value) of SKALD was 78%, and recall (sensitivity) was 27%, when compared with validated reference DSV calls. When retrospectively applied to a cohort of 36 families with acutely ill infants SKALD identified causative DSVs in two. The first was heterozygous deletion of exons 1-3 of MMP21 in trans with a heterozygous frame-shift deletion in two siblings with transposition of the great arteries and heterotaxy. In a newborn female with dysmorphic features, ventricular septal defect and persistent pulmonary hypertension, SKALD identified the breakpoints of a heterozygous, de novo 1p36.32p36.13 deletion. In summary, consensus DSV calling, implemented in an 8-h computational pipeline with parameterised filtering, has the potential to increase the diagnostic yield of WGS in acutely ill neonates and discover novel disease genes.

A 26-hour system of highly sensitive whole genome sequencing for emergency management of genetic diseases.

While the cost of whole genome sequencing (WGS) is approaching the realm of routine medical tests, it remains too tardy to help guide the management of many acute medical conditions. Rapid WGS is imperative in light of growing evidence of its utility in acute care, such as in diagnosis of genetic diseases in very ill infants, and genotype-guided choice of chemotherapy at cancer relapse. In such situations, delayed, empiric, or phenotype-based clinical decisions may meet with substantial morbidity or mortality. We previously described a rapid WGS method, STATseq, with a sensitivity of >96 % for nucleotide variants that allowed a provisional diagnosis of a genetic disease in 50 h. Here improvements in sequencing run time, read alignment, and variant calling are described that enable 26-h time to provisional molecular diagnosis with >99.5 % sensitivity and specificity of genotypes. STATseq appears to be an appropriate strategy for acutely ill patients with potentially actionable genetic diseases.