Leveraging in Vitro Models for Clinically Relevant Rare CYP2D6 Variants in Pharmacogenomics.

Cytochrome P450 2D6 (CYP2D6) is responsible for the metabolism of up to 20% of small-molecule drugs and therefore, may impact the safety and efficacy of medicines in broad therapeutic areas. CYP2D6 is highly polymorphic, and the frequency of variants can differ across racial and ethnic populations, significantly affecting enzymatic function and drug metabolism. However, rare variants of CYP2D6 present a unique challenge for academia, industry, and regulatory agencies alike due to the lack of feasibility of characterizing their clinical relevance in clinical trials, particularly in variants that exhibit population-specific frequencies in racial and ethnic groups that are poorly represented in clinical trials. Despite significant advancement in pharmacogenomics, the substrate specificity and related clinical relevance of these CYP2D6 rare variants remain largely unclear, and further efforts are warranted to characterize the burden of these variants on adverse drug reactions and drug efficacy. Thus, cell-based in vitro systems can be used to inform substrate-specific effects and the overall relevance of a rare variant. Liver microsomes, cell-based expression systems, ex vivo primary samples, and purified variant protein have all been used with various substrates to potentially predict the clinical impact of new substrates. In this review, we identify rare variants of CYP2D6 that demonstrate differences across races in prevalence and thus are often unassessed in clinical trials. Accordingly, we examine current pharmacogenomic in vitro models used to analyze the functional impact of these rare variants in a substrate-specific manner. SIGNIFICANCE STATEMENT: Variants of CYP2D6 play a clinically relevant role in drug metabolism, leading to potential safety and efficacy concerns. Although the influence of prevalent variants is often well characterized, rare variants are traditionally not included in clinical trials. This review captures the clinical relevance of rare variants in CYP2D6 by highlighting in vitro models that analyze their impact on the metabolism of CYP2D6 substrates.

MEGSA: A Powerful and Flexible Framework for Analyzing Mutual Exclusivity of Tumor Mutations.

The central challenges in tumor sequencing studies is to identify driver genes and pathways, investigate their functional relationships, and nominate drug targets. The efficiency of these analyses, particularly for infrequently mutated genes, is compromised when subjects carry different combinations of driver mutations. Mutual exclusivity analysis helps address these challenges. To identify mutually exclusive gene sets (MEGS), we developed a powerful and flexible analytic framework based on a likelihood ratio test and a model selection procedure. Extensive simulations demonstrated that our method outperformed existing methods for both statistical power and the capability of identifying the exact MEGS, particularly for highly imbalanced MEGS. Our method can be used for de novo discovery, for pathway-guided searches, or for expanding established small MEGS. We applied our method to the whole-exome sequencing data for 13 cancer types from The Cancer Genome Atlas (TCGA). We identified multiple previously unreported non-pairwise MEGS in multiple cancer types. For acute myeloid leukemia, we identified a MEGS with five genes (FLT3, IDH2, NRAS, KIT, and TP53) and a MEGS (NPM1, TP53, and RUNX1) whose mutation status was strongly associated with survival (p = 6.7 × 10(-4)). For breast cancer, we identified a significant MEGS consisting of TP53 and four infrequently mutated genes (ARID1A, AKT1, MED23, and TBL1XR1), providing support for their role as cancer drivers.

A Powerful Procedure for Pathway-Based Meta-analysis Using Summary Statistics Identifies 43 Pathways Associated with Type II Diabetes in European Populations.

Meta-analysis of multiple genome-wide association studies (GWAS) has become an effective approach for detecting single nucleotide polymorphism (SNP) associations with complex traits. However, it is difficult to integrate the readily accessible SNP-level summary statistics from a meta-analysis into more powerful multi-marker testing procedures, which generally require individual-level genetic data. We developed a general procedure called Summary based Adaptive Rank Truncated Product (sARTP) for conducting gene and pathway meta-analysis that uses only SNP-level summary statistics in combination with genotype correlation estimated from a panel of individual-level genetic data. We demonstrated the validity and power advantage of sARTP through empirical and simulated data. We conducted a comprehensive pathway-based meta-analysis with sARTP on type 2 diabetes (T2D) by integrating SNP-level summary statistics from two large studies consisting of 19,809 T2D cases and 111,181 controls with European ancestry. Among 4,713 candidate pathways from which genes in neighborhoods of 170 GWAS established T2D loci were excluded, we detected 43 T2D globally significant pathways (with Bonferroni corrected p-values < 0.05), which included the insulin signaling pathway and T2D pathway defined by KEGG, as well as the pathways defined according to specific gene expression patterns on pancreatic adenocarcinoma, hepatocellular carcinoma, and bladder carcinoma. Using summary data from 8 eastern Asian T2D GWAS with 6,952 cases and 11,865 controls, we showed 7 out of the 43 pathways identified in European populations remained to be significant in eastern Asians at the false discovery rate of 0.1. We created an R package and a web-based tool for sARTP with the capability to analyze pathways with thousands of genes and tens of thousands of SNPs.

Association of high-evidence gastric cancer susceptibility loci and somatic gene expression levels with survival.

Eleven high-evidence single-nucleotide polymorphisms (SNPs) at nine loci for gastric cancer (GC) risk were reported, but their associations with survival remain unknown. In this study, we examined associations between SNP and GC survival by anatomic location and histology among 1147 incident cases from the Shanxi Upper Gastrointestinal Genetics Project. We further examined whether SNPs were expression quantitative trait loci in normal and tumor gastric tissues, and whether tumor versus normal somatic mRNA differences in 126 cases were associated with survival. No SNPs were associated with GC survival overall. However, subtype-specific associations were observed for gastric cardia adenocarcinomas at MUC1/TRIM46/1q22 rs2070803 [HRAA versus GA+GG = 2.16; 95% confidence interval (CI) = 1.24-3.78; P = 0.0068] and LTA/TNF/6p21.33 rs1799724 (HRTT+CT versus CC = 1.30; 95% CI = 1.07-1.57; P = 0.0077), and for diffuse-type GC at PSCA/8q24.3 rs2294008 (HRTT versus CT+CC = 1.99; 95% CI = 1.33-2.97; P = 7.8E-04). Rs2294008T was a cis-expression quantitative trait loci for PSCA, upregulating mRNA in normal gastric (ß = 0.60; P = 5.7E-21) and GC (ß = 0.30; P = 0.0089) tissues. Cases in the highest quartile (the smallest downregulation of tumor PSCA) had shortest survival than cases with the most downregulated PSCA (median survival of 0.47 years in the highest quartile versus 3.73 years in the lowest quartile; hazard ratio = 9.70; 95% CI = 2.46-38.4; P = 0.0012). Less striking effects for mRNA levels were observed for MTX1 at 1q22 in gastric cardia adenocarcinoma and for JRK at 8q24.3 in diffuse GC. Our results suggest three high-evidence GC risk loci have prognostic importance in GC subtypes. Future studies in well-characterized independent populations are warranted to validate our findings and further investigate the clinical utility of these variants in predicting GC prognosis.

Inherited variation in circadian rhythm genes and risks of prostate cancer and three other cancer sites in combined cancer consortia.

Circadian disruption has been linked to carcinogenesis in animal models, but the evidence in humans is inconclusive. Genetic variation in circadian rhythm genes provides a tool to investigate such associations. We examined associations of genetic variation in nine core circadian rhythm genes and six melatonin pathway genes with risk of colorectal, lung, ovarian and prostate cancers using data from the Genetic Associations and Mechanisms in Oncology (GAME-ON) network. The major results for prostate cancer were replicated in the Prostate, Lung, Colorectal and Ovarian (PLCO) cancer screening trial, and for colorectal cancer in the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO). The total number of cancer cases and controls was 15,838/18,159 for colorectal, 14,818/14,227 for prostate, 12,537/17,285 for lung and 4,369/9,123 for ovary. For each cancer site, we conducted gene-based and pathway-based analyses by applying the summary-based Adaptive Rank Truncated Product method (sARTP) on the summary association statistics for each SNP within the candidate gene regions. Aggregate genetic variation in circadian rhythm and melatonin pathways were significantly associated with the risk of prostate cancer in data combining GAME-ON and PLCO, after Bonferroni correction (ppathway < 0.00625). The two most significant genes were NPAS2 (pgene = 0.0062) and AANAT (pgene = 0.00078); the latter being significant after Bonferroni correction. For colorectal cancer, we observed a suggestive association with the circadian rhythm pathway in GAME-ON (ppathway = 0.021); this association was not confirmed in GECCO (ppathway = 0.76) or the combined data (ppathway = 0.17). No significant association was observed for ovarian and lung cancer. These findings support a potential role for circadian rhythm and melatonin pathways in prostate carcinogenesis. Further functional studies are needed to better understand the underlying biologic mechanisms.

Genetic modifiers of neurofibromatosis type 1-associated café-au-lait macule count identified using multi-platform analysis.

Neurofibromatosis type 1 (NF1) is an autosomal dominant, monogenic disorder of dysregulated neurocutaneous tissue growth. Pleiotropy, variable expressivity and few NF1 genotype-phenotype correlates limit clinical prognostication in NF1. Phenotype complexity in NF1 is hypothesized to derive in part from genetic modifiers unlinked to the NF1 locus. In this study, we hypothesized that normal variation in germline gene expression confers risk for certain phenotypes in NF1. In a set of 79 individuals with NF1, we examined the association between gene expression in lymphoblastoid cell lines with NF1-associated phenotypes and sequenced select genes with significant phenotype/expression correlations. In a discovery cohort of 89 self-reported European-Americans with NF1 we examined the association between germline sequence variants of these genes with café-au-lait macule (CALM) count, a tractable, tumor-like phenotype in NF1. Two correlated, common SNPs (rs4660761 and rs7161) between DPH2 and ATP6V0B were significantly associated with the CALM count. Analysis with tiled regression also identified SNP rs4660761 as significantly associated with CALM count. SNP rs1800934 and 12 rare variants in the mismatch repair gene MSH6 were also associated with CALM count. Both SNPs rs7161 and rs4660761 (DPH2 and ATP6V0B) were highly significant in a mega-analysis in a combined cohort of 180 self-reported European-Americans; SNP rs1800934 (MSH6) was near-significant in a meta-analysis assuming dominant effect of the minor allele. SNP rs4660761 is predicted to regulate ATP6V0B, a gene associated with melanosome biology. Individuals with homozygous mutations in MSH6 can develop an NF1-like phenotype, including multiple CALMs. Through a multi-platform approach, we identified variants that influence NF1 CALM count.

Whole exome sequencing in families at high risk for Hodgkin lymphoma: identification of a predisposing mutation in the KDR gene.

Hodgkin lymphoma shows strong familial aggregation but no major susceptibility genes have been identified to date. The goal of this study was to identify high-penetrance variants using whole exome sequencing in 17 Hodgkin lymphoma prone families with three or more affected cases or obligate carriers (69 individuals), followed by targeted sequencing in an additional 48 smaller HL families (80 individuals). Alignment and variant calling were performed using standard methods. Dominantly segregating, rare, coding or potentially functional variants were further prioritized based on predicted deleteriousness, conservation, and potential importance in lymphoid malignancy pathways. We selected 23 genes for targeted sequencing. Only the p.A1065T variant in KDR (kinase insert domain receptor) also known as VEGFR2 (vascular endothelial growth factor receptor 2) was replicated in two independent Hodgkin lymphoma families. KDR is a type III receptor tyrosine kinase, the main mediator of vascular endothelial growth factor induced proliferation, survival, and migration. Its activity is associated with several diseases including lymphoma. Functional experiments have shown that p.A1065T, located in the activation loop, can promote constitutive autophosphorylation on tyrosine in the absence of vascular endothelial growth factor and that the kinase activity was abrogated after exposure to kinase inhibitors. A few other promising mutations were identified but appear to be "private". In conclusion, in the largest sequenced cohort of Hodgkin lymphoma families to date, we identified a causal mutation in the KDR gene. While independent validation is needed, this mutation may increase downstream tumor cell proliferation activity and might be a candidate for targeted therapy.

Association between circulating levels of sex steroid hormones and Barrett's esophagus in men: a case-control analysis.

BACKGROUND & AIMS: Esophageal adenocarcinoma is believed to result from the progression of gastroesophageal reflux disease to erosive esophagitis and re-epithelialization of the esophagus with a columnar cell population termed Barrett's esophagus (BE). Men develop BE and esophageal adenocarcinoma more frequently than women, yet little is known about the mechanisms of this difference. We assessed whether sex steroid hormones were associated with BE in a male population. METHODS: We analyzed data from the Barrett's Esophagus Early Detection Case Control Study, based at the Walter Reed National Military Medical Center. Blood samples were collected from 174 men with BE and 213 men without BE (controls, based on endoscopic analysis); 13 sex steroid hormones were measured by mass spectrometry and sex hormone binding globulin was measured by enzyme-linked immunosorbent assay. We also calculated free estradiol, free testosterone, and free dihydrotestosterone (DHT). We used multivariable logistic regression to estimate odds ratios (ORs) and 95% confidence intervals (CIs) adjusted for age, race, smoking status, alcohol consumption, body mass index, heartburn, regurgitation, and gastroesophageal symptom score (excluding heartburn and regurgitation). RESULTS: Levels of free testosterone and free DHT were associated positively with BE risk; patients in the highest quartile for these hormones were most likely to have BE (free testosterone: OR, 5.36; 95% CI, 2.21-13.03; P = .0002; free DHT: OR, 4.25; 95% CI, 1.87-9.66; P = .001). Level of estrone sulfate was associated inversely with BE risk (P for trend = .02). No other hormone was associated with BE risk. Relationships were not modified by age or BMI. CONCLUSIONS: In an analysis of men, levels of free testosterone and free DHT were significantly associated with BE.

Genetic polymorphisms in the 9p21 region associated with risk of multiple cancers.

The chromosome 9p21 region has been implicated in the pathogenesis of multiple cancers. We analyzed 9p21 single nucleotide polymorphisms (SNPs) from eight genome-wide association studies (GWAS) with data deposited in dbGaP, including studies of esophageal squamous cell carcinoma (ESCC), gastric cancer (GC), pancreatic cancer, renal cell carcinoma (RCC), lung cancer (LC), breast cancer (BrC), bladder cancer (BC) and prostate cancer (PrC). The number of subjects ranged from 2252 (PrC) to 7619 (LC). SNP-level analyses for each cancer were conducted by logistic regression or random-effects meta-analysis. A subset-based statistical approach (ASSET) was performed to combine SNP-level P values across multiple cancers. We calculated gene-level P values using the adaptive rank truncated product method. We identified that rs1063192 and rs2157719 in the CDKN2A/2B region were significantly associated with ESCC and rs2764736 (3' of TUSC1) was associated with BC (P ≤ 2.59 × 10(-6)). ASSET analyses identified four SNPs significantly associated with multiple cancers: rs3731239 (CDKN2A intronic) with ESCC, GC and BC (P = 3.96 × 10(-) (4)); rs10811474 (3' of IFNW1) with RCC and BrC (P = 0.001); rs12683422 (LINGO2 intronic) with RCC and BC (P = 5.93 × 10(-) (4)) and rs10511729 (3' of ELAVL2) with LC and BrC (P = 8.63 × 10(-) (4)). At gene level, CDKN2B, CDKN2A and CDKN2B-AS1 were significantly associated with ESCC (P ≤ 4.70 × 10(-) (5)). Rs10511729 and rs10811474 were associated with cis-expression of 9p21 genes in corresponding cancer tissues in the expression quantitative trait loci analysis. In conclusion, we identified several genetic variants in the 9p21 region associated with the risk of multiple cancers, suggesting that this region may contribute to a shared susceptibility across different cancer types.

Genetic variants in fas signaling pathway genes and risk of gastric cancer.

Populations in north central China are at high risk for gastric cancers (GC), and altered FAS-mediated cell signaling and/or apoptosis may contribute to this risk. We examined the association of 554 single nucleotide polymorphisms (SNPs) in 53 Fas signaling-related genes using a pathway-based approach in 1758 GC cases (1126 gastric cardia adenocarcinomas (GCA) and 632 gastric noncardia adenocarcinomas (GNCA)), and 2111 controls from a genome-wide association study (GWAS) of GC in ethnic Chinese. SNP associations with risk of overall GC, GCA and GNCA were evaluated using unconditional logistic regressions controlling for age, sex and study. Gene- and pathway-based associations were tested using the adaptive rank-truncated product (ARTP) method. Statistical significance was evaluated empirically by permutation. Significant pathway-based associations were observed for Fas signaling with risk of overall GC (p = 5.5E-04) and GCA (p = 6.3E-03), but not GNCA (p= 8.1E-02). Among examined genes in the Fas signaling pathway, MAP2K4, FAF1, MAPK8, CASP10, CASP8, CFLAR, MAP2K1, CAP8AP2, PAK2 and IKBKB were associated with risk of GC (nominal p < 0.05), and FAF1 and MAPK8 were significantly associated with risk of both GCA and GNCA (nominal p< 0.05). Our examination of genetic variation in the Fas signaling pathway is consistent with an association of altered Fas signaling and/or apoptosis with risk of GC. As one of the first attempts to investigate a pathway-level association, our results suggest that these genes and the Fas signaling pathway warrant further evaluation in relation to GC risk in other populations.

Characterization of T gene sequence variants and germline duplications in familial and sporadic chordoma.

Chordoma is a rare bone cancer that is believed to originate from notochordal remnants. We previously identified germline T duplication as a major susceptibility mechanism in several chordoma families. Recently, a common genetic variant in T (rs2305089) was significantly associated with the risk of sporadic chordoma. We sequenced all T exons in 24 familial cases and 54 unaffected family members from eight chordoma families (three with T duplications), 103 sporadic cases, and 160 unrelated controls. We also measured T copy number variation in all sporadic cases. We confirmed the association between the previously reported variant rs2305089 and risk of familial [odds ratio (OR) = 2.6, 95% confidence interval (CI) = 0.93, 7.25, P = 0.067] and sporadic chordoma (OR = 2.85, 95% CI = 1.89, 4.29, P < 0.0001). We also identified a second common variant, rs1056048, that was strongly associated with chordoma in families (OR = 4.14, 95% CI = 1.43, 11.92, P = 0.0086). Among sporadic cases, another common variant (rs3816300) was significantly associated with risk when jointly analyzed with rs2305089. The association with rs3816300 was significantly stronger in cases with early age onset. In addition, we identified three rare variants that were only observed among sporadic chordoma cases, all of which have potential functional relevance based on in silico predictions. Finally, we did not observe T duplication in any sporadic chordoma case. Our findings further highlight the importance of the T gene in the pathogenesis of both familial and sporadic chordoma and suggest a complex susceptibility related to T.

Genetic variants reflecting higher vitamin e status in men are associated with reduced risk of prostate cancer.

Vitamin E (α-tocopherol) plays a key role in the regulation of cell growth and differentiation and has been studied as a potential chemopreventive agent for prostate cancer. The association of serum vitamin E concentrations with cancer risk may be modified by genetic variations in vitamin E-related genes. We examined whether variants in vitamin E-related genes were associated with risk of prostate cancer in a nested case-control study using 483 prostate cancer cases and 542 matched controls of European ancestry from a large U.S. multicenter trial that had available measurements of serum vitamin E concentrations and genotyping of 3 genome-wide association study meta-analysis-identified single-nucleotide polymorphisms (SNPs) associated with circulating vitamin E. ORs and 95% CIs were calculated using unconditional logistic regression adjusted for age, family history of prostate cancer, and serum total cholesterol. Findings suggest lower prostate cancer risk for men whose genotypes reflect higher vitamin E (i.e., α-tocopherol) status. An SNP (rs964184) near budding-site selection protein 13 (yeast) (BUD13), zinc finger protein 259 (ZNF259), and apolipoprotein A5 (APOA5) on 11q23.3 was significantly associated with prostate cancer risk (per-allele OR = 0.75; 95% CI: 0.58, 0.98; P-trend = 0.03). The association between rs964184 and prostate cancer risk was stronger among homozygous carriers of the minor allele (OR = 0.27; 95% CI: 0.09, 0.83). Another variant, rs11057830 in scavenger receptor class-B member 1 (SCARB1) on 12p24.31, approached statistical significance (OR = 0.32; 95% CI: 0.10, 1.01, P = 0.05; 2 minor allele copies). This study suggests that polymorphisms near BUD13/ZNF259/APOA5, involved in vitamin E transport and metabolism, may be associated with lower risk of prostate cancer. This trial was registered at clinicaltrials.gov as NCT00002540.

Genetic variants in sex hormone metabolic pathway genes and risk of esophageal squamous cell carcinoma.

In China, esophageal cancer is the fourth leading cause of cancer death where essentially all cases are histologically esophageal squamous cell carcinoma (ESCC), in contrast to esophageal adenocarcinoma in the West. Globally, ESCC is 2.4 times more common among men than women and recently it has been suggested that sex hormones may be associated with the risk of ESCC. We examined the association between genetic variants in sex hormone metabolic genes and ESCC risk in a population from north central China with high-incidence rates. A total of 1026 ESCC cases and 1452 controls were genotyped for 797 unique tag single-nucleotide polymorphisms (SNPs) in 51 sex hormone metabolic genes. SNP-, gene- and pathway-based associations with ESCC risk were evaluated using unconditional logistic regression adjusted for age, sex and geographical location and the adaptive rank truncated product (ARTP) method. Statistical significance was determined through use of permutation for pathway- and gene-based associations. No associations were observed for the overall sex hormone metabolic pathway (P = 0.14) or subpathways (androgen synthesis: P = 0.30, estrogen synthesis: P = 0.15 and estrogen removal: P = 0.19) with risk of ESCC. However, six individual genes (including SULT2B1, CYP1B1, CYP3A7, CYP3A5, SHBG and CYP11A1) were significantly associated with ESCC risk (P < 0.05). Our examination of genetic variation in the sex hormone metabolic pathway is consistent with a potential association with risk of ESCC. These positive findings warrant further evaluation in relation to ESCC risk and replication in other populations.

Genetic variants in DNA repair pathway genes and risk of esophageal squamous cell carcinoma and gastric adenocarcinoma in a Chinese population.

The DNA repair pathways help to maintain genomic integrity and therefore genetic variation in the pathways could affect the propensity to develop cancer. Selected germline single nucleotide polymorphisms (SNPs) in the pathways have been associated with esophageal cancer and gastric cancer (GC) but few studies have comprehensively examined the pathway genes. We aimed to investigate associations between DNA repair pathway genes and risk of esophageal squamous cell carcinoma (ESCC) and GC, using data from a genome-wide association study in a Han Chinese population where ESCC and GC are the predominant cancers. In sum, 1942 ESCC cases, 1758 GC cases and 2111 controls from the Shanxi Upper Gastrointestinal Cancer Genetics Project (discovery set) and the Linxian Nutrition Intervention Trials (replication set) were genotyped for 1675 SNPs in 170 DNA repair-related genes. Logistic regression models were applied to evaluate SNP-level associations. Gene- and pathway-level associations were determined using the resampling-based adaptive rank-truncated product approach. The DNA repair pathways overall were significantly associated with risk of ESCC (P = 6.37 × 10(-4)), but not with GC (P = 0.20). The most significant gene in ESCC was CHEK2 (P = 2.00 × 10(-6)) and in GC was CLK2 (P = 3.02 × 10(-4)). We observed several other genes significantly associated with either ESCC (SMUG1, TDG, TP53, GTF2H3, FEN1, POLQ, HEL308, RAD54B, MPG, FANCE and BRCA1) or GC risk (MRE11A, RAD54L and POLE) (P < 0.05). We provide evidence for an association between specific genes in the DNA repair pathways and the risk of ESCC and GC. Further studies are warranted to validate these associations and to investigate underlying mechanisms.

Considerations for Use of Pharmacodynamic Biomarkers to Support Biosimilar Development - (III) A Randomized Trial with Interferon Beta-1a Products.

The US Food and Drug Administration (FDA) has taken steps to bring efficiency to the development of biosimilars, including establishing guidance for the use of pharmacokinetic and pharmacodynamic (PD) similarity study data without a comparative clinical study with efficacy end point(s). To better understand the potential role for PD biomarkers in biosimilar development and inform best practices for biomarker selection and analysis, we conducted a randomized, double-blinded, placebo-controlled, single-dose, parallel-arm clinical study in healthy participants. Eighty-four healthy participants (n = 12 per dose arm) received either placebo or one of three doses of either interferon ß-1a (7.5-30 µg) or pegylated interferon ß-1a (31.25-125 µg) to evaluate the maximum change from baseline and the baseline-adjusted area under the effect curve for the biomarkers neopterin in serum and myxovirus resistance protein 1 in blood. Both PD biomarkers increased following product administration with clear separation from baseline (neopterin: 3.4-fold and 3.9-fold increase for interferon ß-1a and pegylated interferon ß-1a, respectively; myxovirus resistance protein 1: 19.0-fold and 47.2-fold increase for interferon ß-1a and pegylated interferon ß-1a, respectively). The dose-response curves support that therapeutic doses were adequately sensitive to detect differences in both PD biomarkers for consideration in a PD similarity study design. Because baseline levels of both biomarkers are low compared with on-treatment values, there was little difference in using PD measures adjusted to baseline compared with the results without baseline adjustment. This study illustrates potential methodologies for evaluating PD biomarkers and an approach to address information gaps when limited information is publicly available for one or more PD biomarkers.

Evaluating the Utility of Proteomics for the Identification of Circulating Pharmacodynamic Biomarkers of IFNß-1a Biologics.

Proteomics has the potential to identify pharmacodynamic (PD) biomarkers for similarity assessment of proposed biosimilars without relying on clinical efficacy end points. In this study, with 36 healthy participants randomized to therapeutic doses of interferon-beta 1a products (IFNß-1a) or pegylated-IFNß-1a (pegIFNß-1a) approved to treat multiple sclerosis or placebo, we evaluated the utility of a proteomic assay that profiles > 7,000 plasma proteins. IFNß-1a and pegIFNß-1a resulted in 248 and 528 differentially expressed protein analytes, respectively, between treatment and placebo groups over the time course. Thirty-one proteins were prioritized based on a maximal fold change ≥ 2 from baseline, baseline adjusted area under the effect curve (AUEC) and overlap between the 2 products. Of these, the majority had a significant AUEC compared with placebo in response to either product; 8 proteins showed > 4-fold maximal change from baseline. We identified previously reported candidates, beta-2microglobulin and interferon-induced GTP-binding protein (Mx1) with ~ 50% coefficient of variation (CV) for AUEC, and many new candidates (including I-TAC, C1QC, and IP-10) with CVs ranging from 26%-129%. Upstream regulator analysis of differentially expressed proteins predicted activation of IFNß1 signaling as well as other cytokine, enzyme, and transcription signaling networks by both products. Although independent replication is required to confirm present results, our study demonstrates the utility of proteomics for the identification of individual and composite candidate PD biomarkers that may be leveraged to support clinical pharmacology studies for biosimilar approvals, especially when biologics have complex mechanisms of action or do not have previously characterized PD biomarkers.

Evidence for alteration of EZH2, BMI1, and KDM6A and epigenetic reprogramming in human papillomavirus type 16 E6/E7-expressing keratinocytes.

A number of epigenetic alterations occur in both the virus and host cellular genomes during human papillomavirus (HPV)-associated carcinogenesis, and investigations of such alterations, including changes in chromatin proteins and histone modifications, have the potential to lead to therapeutic epigenetic reversion. We report here that transformed HPV16 E6/E7-expressing primary human foreskin keratinocytes (HFKs) (E6/E7 cells) demonstrate increased expression of the PRC2 methyltransferase EZH2 at both the mRNA and protein levels but do not exhibit the expected increase in trimethylated H3K27 (H3K27me3) compared to normal keratinocytes. In contrast, these cells show a reduction in global H3K27me3 levels in vitro, as well as upregulation of the KDM6A demethylase. We further show for the first time that transformation with the HPV16 E6 and E7 oncogenes also results in an increase in phosphorylated EZH2 serine 21 (P-EZH2-Ser21), mediated by active Akt, and in a downregulation of the PRC1 protein BMI1 in these cells. High-grade squamous cervical intraepithelial lesions also showed a loss of H3K27me3 in the presence of increased expression of EZH2. Correlating with the loss of H3K27me3, E6/E7 cells exhibited derepression of specific EZH2-, KMD6A-, and BMI1-targeted HOX genes. These results suggest that the observed reduction in H3K27me3 may be due to a combination of reduced activities/levels of specific polycomb proteins and increases in demethylases. The dysregulation of multiple chromatin proteins resulting in the loss of global H3K27me3 and the transcriptional reprogramming in HPV16 E6/E7-infected cells could provide an epigenetic signature associated with risk and/or progression of HPV16-associated cancers, as well as the potential for epigenetic reversion in the future.

Expanding Approved Patient Populations for Rare Disease Treatment Using In Vitro Data.

In vitro cell-based data can be used to support the extension of pharmaceutical approval to patient subsets with unique genetic variants. A set of conditions should be satisfied to support the extension of approval. The disease mechanism should be well described, and the impact of variants on protein function should be reasonably understood. The incidence data should show that clinical trials for the variants in question are not practical. The overall safety and efficacy of the drug should be clear in adequate and well-controlled clinical trials. The clinical trial should include patients found to be responders and nonresponders so that both positive and negative predictive power of the in vitro assay may be measured. The mechanism of action of the drug should be clearly defined and should be consistent with the disease mechanism. The assay system should be qualified, including the following points: (i) each variant construct should be confirmed by bidirectional sequencing; (ii) the in vitro assay should directly measure the variant protein function in comparison with the reference protein; (iii) the assay should be formally validated to the extent possible, clearly demonstrating precision, reproducibility, and sensitivity used to support the efficacy claim; and (iv) the primary data should be available for inspection and analytical validation. The overall goal is a robust and validated cell-based system that can be shown to predict the outcome of targeted therapy.

New science, drug regulation, and emergent public health issues: The work of FDA's division of applied regulatory science.

The U.S. Food and Drug Administration (FDA) Division of Applied Regulatory Science (DARS) moves new science into the drug review process and addresses emergent regulatory and public health questions for the Agency. By forming interdisciplinary teams, DARS conducts mission-critical research to provide answers to scientific questions and solutions to regulatory challenges. Staffed by experts across the translational research spectrum, DARS forms synergies by pulling together scientists and experts from diverse backgrounds to collaborate in tackling some of the most complex challenges facing FDA. This includes (but is not limited to) assessing the systemic absorption of sunscreens, evaluating whether certain drugs can convert to carcinogens in people, studying drug interactions with opioids, optimizing opioid antagonist dosing in community settings, removing barriers to biosimilar and generic drug development, and advancing therapeutic development for rare diseases. FDA tasks DARS with wide ranging issues that encompass regulatory science; DARS, in turn, helps the Agency solve these challenges. The impact of DARS research is felt by patients, the pharmaceutical industry, and fellow regulators. This article reviews applied research projects and initiatives led by DARS and conducts a deeper dive into select examples illustrating the impactful work of the Division.