Transcriptomic Profiles Associated with Experimental Placebo Effects in Chronic Pain.

Gene expression networks associated with placebo effects are understudied; in this study, we identified transcriptomic profiles associated with placebo responsivity. Participants suffering from chronic pain underwent a verbal suggestion and conditioning paradigm with individually tailored thermal painful stimulations to elicit conditioned placebo effects. Participants reported pain intensity on a visual analog scale (VAS) anchored from zero = no pain to 100 = maximum imaginable pain. RNA was extracted from venous blood and RNA sequencing and validation tests were performed to identify differentially expressed genes (DEGs) associated with placebo effects, controlling for sex and level of pain. Unbiased enrichment analyses were performed to identify biological processes associated with placebo effects. Of the 10,700 protein-coding genes that passed quality control filters, 667 were found to be associated with placebo effects (FDR <0.05). Most genes (97%) upregulated were associated with larger placebo effects. The 17 top transcriptome-wide significant genes were further validated via RT-qPCR in an independent cohort of chronic pain participants. Six of them (CCDC85B, FBXL15, HAGH, PI3, SELENOM, and TNFRSF4) showed positive and significant (P < 0.05) correlation with placebo effects in the cohort. The overall DEGs were highly enriched in regulation of expression of SLITs and ROBOs (R-HSA-9010553, FDR = 1.26e-33), metabolism of RNA (R-HSA-8953854, FDR = 1.34e-30), Huntington's disease (hsa05016, FDR = 9.84e-31), and ribosome biogenesis (GO:0042254, FDR = 2.67e-15); alternations in these pathways might jeopardize the proneness to elicit placebo effects. Future studies are needed to replicate this finding and better understand the unique molecular dynamics of people who are more or less affected by pain and placebo.

Rapid effects of valproic acid on the fetal brain transcriptome: Implications for brain development and autism.

There is an increased incidence of autism among the children of women who take the anti-epileptic, mood-stabilizing drug, valproic acid (VPA) during pregnancy; moreover, exposure to VPA in utero causes autistic-like symptoms in rodents and non-human primates. Analysis of RNA-seq data obtained from E12.5 fetal mouse brains 3 hours after VPA administration to the pregnant dam revealed that VPA rapidly and significantly increased or decreased the expression of approximately 7,300 genes. No significant sex differences in VPA-induced gene expression were observed. Expression of 399 autism risk genes was significantly altered by VPA as was expression of 255 genes that have been reported to play fundamental roles in fetal brain development but are not otherwise linked to autism. Expression of genes associated with intracellular signaling pathways, neurogenesis, and excitation-inhibition balance as well as synaptogenesis, neuronal fate determination, axon and dendritic development, neuroinflammation, circadian rhythms, and epigenetic modulation of gene expression was dysregulated by VPA. The goal of this study was to identify mouse genes that are: (a) significantly up- or down-regulated by VPA in the fetal brain and (b) known to be associated with autism and/or to play a role in embryonic neurodevelopmental processes, perturbation of which has the potential to alter brain connectivity and, consequently behavior, in the adult. The set of genes meeting these criteria provides potential targets for future hypothesis-driven studies to elucidate the proximal causes of errors in brain connectivity underlying neurodevelopmental disorders such as autism.

Nociceptive and Transcriptomic Responses in a Swine Diabetic Wound Model Treated With a Topical Angiotensin 1 Receptor Antagonist.

BACKGROUND: Painful, treatment-resistant wounds are prevalent among diabetic patients and significantly affect health-related quality of life (HRQOL). Topical treatments may help alleviate pain without risk of dependence or side effects. However, there is a lack of topical wound compounds targeting pain-specific receptors. One possible target is proinflammatory angiotensin 1 receptor (AT1R), which is upregulated in diabetic skin and has been implicated in nociception. OBJECTIVES: We investigated the effects of topical valsartan, an AT1R antagonist, on pain (nociceptive thresholds) and gene expression changes (transcriptomics) in a swine model of diabetic wounds. METHODS: Eight wounds were surgically induced in diabetic, hyperglycemic Yucatan miniature swine ( n = 4). Topical AT1R antagonist was applied to wounds on one side and vehicle on the other side. Nocifensive testing was conducted at baseline and then weekly, beginning 7 days after wound induction. Mechanical and thermal stimuli were applied to the wound margins until a nocifensive reaction was elicited or a predetermined cutoff was reached. After 7 weeks of testing, tissue from the dorsal horn, dorsal root ganglion, and wounds were sequenced and analyzed with DESeq2. Unbiased pathway analyses using Metascape were conducted on differentially expressed genes. RESULTS: There was no significant difference in mechanical tolerance threshold between AT1R antagonist-treated and vehicle-treated wounds ( p = .106). Thermal tolerance was significantly higher in AT1R antagonist-treated wounds compared to vehicle-treated ( p = .015). Analysis of differentially expressed genes revealed enriched pathways of interest: interleukin-18 signaling in dorsal horn laminae IV-V and sensory perception of mechanical stimulus in wound tissue. DISCUSSION: In this study, wounds modeling diabetic ulcers were created in hyperglycemic swine and treated with a topical AT1R antagonist. AT1R-antagonist-treated wounds had a higher tolerance threshold than vehicle-treated wounds for thermal hyperalgesia, but not mechanical allodynia. Pathway analyses of differentially expressed genes revealed several pathways of interest for future pain research. Although further studies are needed to confirm the findings, this study can improve nursing care by providing information about a potential future treatment that may be used to decrease pain and improve HRQOL in patients with diabetic wounds.

Differential Gene Expression in Pain-Related Genes are not Affected by the Presence of Dementia.

BACKGROUND: Prior work has demonstrated differences in the transcriptome between those with and without chronic musculoskeletal pain. AIMS: The aim of this study was to explore whether pain-related gene expression is similar between individuals with and without dementia. DESIGN: This was a descriptive study using a one-time assessment. SETTINGS: PARTICIPANTS/SUBJECTS: A total of 20 older adults living in a continuing care retirement community, 50% of whom had dementia were inlcuded in this study. All were female and the mean age of participants was 89 (SD = 6). METHODS: Pain was evaluated based on the PROMIS Pain Intensity Short Form 3a. Whole blood was collected by venipuncture into Tempus vacutainer tubes (3 ml) and the RNA was extracted at the Translational Genomics Laboratory at the University of Maryland Baltimore. Analyses included a differential expression analysis, a weighted gene co-expression network analysis, and a pathway enrichment analysis. RESULTS: Eighty-three genes were differentially expressed between individuals with and without pain (p <.05). After normalizing gene counts and removing the low expressed genes, 18,028 genes were left in the final analysis. There was no clustering of the samples related to study variables of pain or dementia. CONCLUSION: The findings from this study provided some preliminary support that pain-related gene expression is similar between individuals with and without dementia.

A single-cell genomic atlas for maturation of the human cerebellum during early childhood.

Inflammation early in life is a clinically established risk factor for autism spectrum disorders and schizophrenia, yet the impact of inflammation on human brain development is poorly understood. The cerebellum undergoes protracted postnatal maturation, making it especially susceptible to perturbations contributing to the risk of developing neurodevelopmental disorders. Here, using single-cell genomics of postmortem cerebellar brain samples, we characterized the postnatal development of cerebellar neurons and glia in 1- to 5-year-old children, comparing individuals who had died while experiencing inflammation with those who had died as a result of an accident. Our analyses revealed that inflammation and postnatal cerebellar maturation are associated with extensive, overlapping transcriptional changes primarily in two subtypes of inhibitory neurons: Purkinje neurons and Golgi neurons. Immunohistochemical analysis of a subset of these postmortem cerebellar samples revealed no change to Purkinje neuron soma size but evidence for increased activation of microglia in those children who had experienced inflammation. Maturation-associated and inflammation-associated gene expression changes included genes implicated in neurodevelopmental disorders. A gene regulatory network model integrating cell type-specific gene expression and chromatin accessibility identified seven temporally specific gene networks in Purkinje neurons and suggested that inflammation may be associated with the premature down-regulation of developmental gene expression programs.

Genome wide association joint analysis reveals 99 risk loci for pain susceptibility and pleiotropic relationships with psychiatric, metabolic, and immunological traits.

Chronic pain is at epidemic proportions in the United States, represents a significant burden on our public health system, and is coincident with a growing opioid crisis. While numerous genome-wide association studies have been reported for specific pain-related traits, many of these studies were underpowered, and the genetic relationship among these traits remains poorly understood. Here, we conducted a joint analysis of genome-wide association study summary statistics from seventeen pain susceptibility traits in the UK Biobank. This analysis revealed 99 genome-wide significant risk loci, 65 of which overlap loci identified in earlier studies. The remaining 34 loci are novel. We applied leave-one-trait-out meta-analyses to evaluate the influence of each trait on the joint analysis, which suggested that loci fall into four categories: loci associated with nearly all pain-related traits; loci primarily associated with a single trait; loci associated with multiple forms of skeletomuscular pain; and loci associated with headache-related pain. Overall, 664 genes were mapped to the 99 loci by genomic proximity, eQTLs, and chromatin interaction and ~15% of these genes showed differential expression in individuals with acute or chronic pain compared to healthy controls. Risk loci were enriched for genes involved in neurological and inflammatory pathways. Genetic correlation and two-sample Mendelian randomization indicated that psychiatric, metabolic, and immunological traits mediate some of these effects.

Diagnostic Performance of a Tumor Marker Gene Test to Personalize Serum CA19-9 Reference Ranges.

PURPOSE: CA19-9 synthesis is influenced by common variants in the fucosyltransferase (FUT) enzymes FUT3 and FUT2. We developed a clinical test to detect FUT variants, and evaluated its diagnostic performance for pancreatic ductal adenocarcinoma (PDAC). EXPERIMENTAL DESIGN: A representative set of controls from the Cancer of the Pancreas Screening study was identified for each FUT functional group. Diagnostic sensitivity was determined first in a testing set of 234 PDAC cases, followed by a 134-case validation set, all of whom had undergone resection with curative intent without neoadjuvant therapy. Tumor marker gene testing was performed in the Johns Hopkins Molecular Diagnostics Laboratory. CA19-9 levels were measured in the Hopkins Clinical Chemistry lab. Receiver operating characteristic (ROC) curve analysis was used to evaluate the discriminative ability of CA19-9 alone versus with the gene test. RESULTS: Applying the CA19-9 standard cutoff (<36 U/mL) to all 716 subjects yielded a 68.8% sensitivity in the test set of cases, 67.2% in the validation set, at 91.4% specificity. Applying 99th percentile cutoffs according to each individual's FUT group (3, 34.9, 41.8, and 89.2, for the FUT3-null, FUT-low, FUT-intermediate, and FUT-high groups, respectively) yielded a diagnostic sensitivity for CA19-9 in the first set of cases of 66.7%, 65.7% in the validation set, at 98.9% specificity. ROC analysis for CA19-9 alone yielded an AUC of 0.84; with the tumor marker gene test, AUC improved to 0.92 (P < 0.001). CONCLUSIONS: Using a tumor marker gene test to personalize an individual's CA19-9 reference range significantly improves diagnostic accuracy.

Rapid effects of valproic acid on the fetal brain transcriptome: Implications for brain development and autism.

There is an increased incidence of autism among the children of women who take the anti-epileptic, mood stabilizing drug, valproic acid (VPA) during pregnancy; moreover, exposure to VPA in utero causes autistic-like symptoms in rodents and non-human primates. Analysis of RNAseq data ob-tained from E12.5 fetal mouse brains 3 hours after VPA administration revealed that VPA significant-ly increased or decreased the expression of approximately 7,300 genes. No significant sex differ-ences in VPA-induced gene expression were observed. Expression of genes associated with neu-rodevelopmental disorders (NDDs) such as autism as well as neurogenesis, axon growth and syn-aptogenesis, GABAergic, glutaminergic and dopaminergic synaptic transmission, perineuronal nets, and circadian rhythms was dysregulated by VPA. Moreover, expression of 399 autism risk genes was significantly altered by VPA as was expression of 252 genes that have been reported to play fundamental roles in the development of the nervous system but are not otherwise linked to autism. The goal of this study was to identify mouse genes that are: (a) significantly up- or down-regulated by VPA in the fetal brain and (b) known to be associated with autism and/or to play a role in embryonic neurodevelopmental processes, perturbation of which has the potential to alter brain connectivity in the postnatal and adult brain. The set of genes meeting these criteria pro-vides potential targets for future hypothesis-driven approaches to elucidating the proximal underly-ing causes of defective brain connectivity in NDDs such as autism.

Genome-wide significant risk loci for mood disorders in the Old Order Amish founder population.

Genome-wide association studies (GWAS) of mood disorders in large case-control cohorts have identified numerous risk loci, yet pathophysiological mechanisms remain elusive, primarily due to the very small effects of common variants. We sought to discover risk variants with larger effects by conducting a genome-wide association study of mood disorders in a founder population, the Old Order Amish (OOA, n = 1,672). Our analysis revealed four genome-wide significant risk loci, all of which were associated with >2-fold relative risk. Quantitative behavioral and neurocognitive assessments (n = 314) revealed effects of risk variants on sub-clinical depressive symptoms and information processing speed. Network analysis suggested that OOA-specific risk loci harbor novel risk-associated genes that interact with known neuropsychiatry-associated genes via gene interaction networks. Annotation of the variants at these risk loci revealed population-enriched, non-synonymous variants in two genes encoding neurodevelopmental transcription factors, CUX1 and CNOT1. Our findings provide insight into the genetic architecture of mood disorders and a substrate for mechanistic and clinical studies.

Differential Gene Expression Among Patients With Heart Failure Experiencing Pain.

BACKGROUND: Chronic pain is frequently experienced by patients with heart failure (HF) and is associated with higher mortality, higher symptom burden, and worsened health-related quality of life. However, the genomic mechanisms underlying chronic pain in HF are understudied. Building an understanding of the mechanistic underpinnings of pain may inform novel interventions. OBJECTIVE: The objective was to identify genes associated with pain from messenger RNA sequence data collected from patients with HF with and without pain. METHODS: The current study analyzed data from 40 patients with HF previously enrolled in a clinical trial. Pain presence was measured using the Health Utilities Index Mark-3. Genes were tested for differential expression using DESeq2, and differentially expressed genes were analyzed for protein-protein interaction (PPI) and relevant ontological pathways using Metascape. Genes located within the core of the PPI network were considered key in disease-relevant biological pathways. Differentially expressed genes within this PPI network were reviewed in existing literature to narrow down candidate genes of interest. These target genes of interest were reanalyzed in a second sample of 24 patients with HF using validation quantitative polymerase chain reaction. RESULTS: A total of 334 genes (279 upregulated, 55 downregulated) were differentially expressed between patients with and without pain in the primary sample of 40. These genes were largely aligned with neutrophil degranulation pathways. Seven genes of interest were identified from a core network of 15 co-expressed genes in the PPI network and existing literature. Three of these seven genes, matrix metallopeptidase 8 ( MMP8 ), proprotein convertase subtilisin/kexin type 9 ( PCSK9 ), and neutrophil defensin 3 ( DEFA3 ), were upregulated in patients with pain versus without pain in both the primary and validation samples. All seven genes of interest are involved in immune, inflammatory, and atherosclerotic processes. DISCUSSION: These results identify potential genes that may play a mechanistic role in chronic pain in HF. Further research is needed to evaluate these potential genes among clearly delineated pain phenotypes.

Evaluating depressive symptoms, BDNF Val66Met, and APOE-ε4 as moderators of response to computerized cognitive training in heart failure.

BACKGROUND: Depressive symptoms, brain-derived neurotrophic factor (BDNF) Val66Met, and apolipoprotein (APOE)-ε4 may moderate response to computerized cognitive training (CCT) interventions among patients with heart failure (HF). OBJECTIVES: The purpose of this study was to examine moderators of intervention response to CCT over 8 months among patients with HF enrolled in a 3-arm randomized controlled trial. Outcomes were memory, serum BDNF, working memory, instrumental activities of daily living (IADLs), and health-related quality of life (HRQL). METHODS: 256 patients with HF were randomized to CCT, computerized crossword puzzles active control, and usual care control groups for 8 weeks. Data were collected at enrollment, baseline, 10 weeks, and 4 and 8 months. Mixed effects models were computed to evaluate moderators. RESULTS: As previously reported, there were no statistically significant group by time effects in outcomes among the 3 groups over 8 months. Tests of moderation indicated that depressive symptoms and presence of BDNF Val66Met and APOE-ε4 were not statistically significant moderators of intervention response in outcomes of delayed recall memory, serum BDNF, working memory, IADLs, and HRQL. In post hoc analysis evaluating baseline global cognitive function, gender, age, and HF severity as moderators, no significant effects were found. HF severity was imbalanced among groups (P = .049) which may have influenced results. CONCLUSIONS: Studies are needed to elucidate biological mechanisms of cognitive dysfunction in HF and test novel interventions to improve memory, serum BDNF, working memory, IADLs and HRQL. Patients may need to be stratified or randomized by HF severity within intervention trials.

Early and Late Transcriptional Changes in Blood, Neural, and Colon Tissues in Rat Models of Stress-Induced and Comorbid Pain Hypersensitivity Reveal Regulatory Roles in Neurological Disease.

Background: Irritable bowel syndrome (IBS) and temporomandibular disorder (TMD) are two chronic pain conditions that frequently overlap in the same individual, more commonly in women. Stress is a significant risk factor, exacerbating or triggering one or both conditions. However, the mechanisms underlying IBS-TMD co-morbidity are mostly unknown. Aim: To detect both specific and common stress-induced visceral hypersensitivity (SIH) and comorbid TMD-IBS pain hypersensitivity (CPH) genetic signatures over time. Method: Twenty-four female rats were randomly assigned to one of three experimental groups: naïve, SIH, and CPH (orofacial pain plus stress). RNA was extracted from blood, colon, spinal cord, and dorsal root ganglion 1 or 7 weeks after the stress paradigm. We combined differential gene expression and co-expression network analyses to define both SIH and CPH expression profiles across tissues and time. Results: The transcriptomic profile in blood and colon showed increased expression of genes enriched in inflammatory and neurological biological processes in CPH compared to SIH rats, both at 1 and 7 weeks after stress. In lumbosacral spinal tissue, both SIH and CPH rats compared to naïve revealed decreased expression of genes related to synaptic activity and increased expression of genes enriched in "angiogenesis," "Neurotrophin," and "PI3K-Akt" pathways. Compared to SIH, CPH rats showed increased expression of angiogenesis-related genes 1 week after exposure to stress, while 7 weeks post-stress the expression of these genes was higher in SIH rats. In dorsal root ganglia (DRG), CPH rats showed decreased expression of immune response genes at week 1 and inhibition of nerve myelination genes at 7 weeks compared to naïve. For all tissues, we observed higher expression of genes involved in ATP production in SIH compared to CPH at 1 week and this was reversed 7 weeks after the induction of stress. Conclusion: Our study highlights an increased inflammatory response in CPH compared to SIH rats in the blood and colon. DRG and spinal transcriptomic profiles of both CPH and SIH rats showed inhibition of synaptic activity along with activation of angiogenesis. Targeting these biological processes may lead to a more profound understanding of the mechanisms underlying IBS-TMD comorbidities and new diagnostic and therapeutic strategies.

Smoking Modifies Pancreatic Cancer Risk Loci on 2q21.3.

Germline variation and smoking are independently associated with pancreatic ductal adenocarcinoma (PDAC). We conducted genome-wide smoking interaction analysis of PDAC using genotype data from four previous genome-wide association studies in individuals of European ancestry (7,937 cases and 11,774 controls). Examination of expression quantitative trait loci data from the Genotype-Tissue Expression Project followed by colocalization analysis was conducted to determine whether there was support for common SNP(s) underlying the observed associations. Statistical tests were two sided and P < 5 × 10-8 was considered statistically significant. Genome-wide significant evidence of qualitative interaction was identified on chr2q21.3 in intron 5 of the transmembrane protein 163 (TMEM163) and upstream of the cyclin T2 (CCNT2). The most significant SNP using the Empirical Bayes method, in this region that included 45 significantly associated SNPs, was rs1818613 [per allele OR in never smokers 0.87, 95% confidence interval (CI), 0.82-0.93; former smokers 1.00, 95% CI, 0.91-1.07; current smokers 1.25, 95% CI 1.12-1.40, P interaction = 3.08 × 10-9). Examination of the Genotype-Tissue Expression Project data demonstrated an expression quantitative trait locus in this region for TMEM163 and CCNT2 in several tissue types. Colocalization analysis supported a shared SNP, rs842357, in high linkage disequilibrium with rs1818613 (r 2 = 0. 94) driving both the observed interaction and the expression quantitative trait loci signals. Future studies are needed to confirm and understand the differential biologic mechanisms by smoking status that contribute to our PDAC findings. SIGNIFICANCE: This large genome-wide interaction study identifies a susceptibility locus on 2q21.3 that significantly modified PDAC risk by smoking status, providing insight into smoking-associated PDAC, with implications for prevention.

A Transcriptome-Wide Association Study Identifies Novel Candidate Susceptibility Genes for Pancreatic Cancer.

BACKGROUND: Although 20 pancreatic cancer susceptibility loci have been identified through genome-wide association studies in individuals of European ancestry, much of its heritability remains unexplained and the genes responsible largely unknown. METHODS: To discover novel pancreatic cancer risk loci and possible causal genes, we performed a pancreatic cancer transcriptome-wide association study in Europeans using three approaches: FUSION, MetaXcan, and Summary-MulTiXcan. We integrated genome-wide association studies summary statistics from 9040 pancreatic cancer cases and 12 496 controls, with gene expression prediction models built using transcriptome data from histologically normal pancreatic tissue samples (NCI Laboratory of Translational Genomics [n = 95] and Genotype-Tissue Expression v7 [n = 174] datasets) and data from 48 different tissues (Genotype-Tissue Expression v7, n = 74-421 samples). RESULTS: We identified 25 genes whose genetically predicted expression was statistically significantly associated with pancreatic cancer risk (false discovery rate < .05), including 14 candidate genes at 11 novel loci (1p36.12: CELA3B; 9q31.1: SMC2, SMC2-AS1; 10q23.31: RP11-80H5.9; 12q13.13: SMUG1; 14q32.33: BTBD6; 15q23: HEXA; 15q26.1: RCCD1; 17q12: PNMT, CDK12, PGAP3; 17q22: SUPT4H1; 18q11.22: RP11-888D10.3; and 19p13.11: PGPEP1) and 11 at six known risk loci (5p15.33: TERT, CLPTM1L, ZDHHC11B; 7p14.1: INHBA; 9q34.2: ABO; 13q12.2: PDX1; 13q22.1: KLF5; and 16q23.1: WDR59, CFDP1, BCAR1, TMEM170A). The association for 12 of these genes (CELA3B, SMC2, and PNMT at novel risk loci and TERT, CLPTM1L, INHBA, ABO, PDX1, KLF5, WDR59, CFDP1, and BCAR1 at known loci) remained statistically significant after Bonferroni correction. CONCLUSIONS: By integrating gene expression and genotype data, we identified novel pancreatic cancer risk loci and candidate functional genes that warrant further investigation.

Analysis of Heritability and Genetic Architecture of Pancreatic Cancer: A PanC4 Study.

BACKGROUND: Pancreatic cancer is the fourth-leading cause of cancer death in both men and women in the United States. The currently identified common susceptibility loci account for a small fraction of estimated heritability. We sought to estimate overall heritability of pancreatic cancer and partition the heritability by variant frequencies and functional annotations. METHODS: Analysis using the genome-based restricted maximum likelihood method (GREML) was conducted on Pancreatic Cancer Case-Control Consortium (PanC4) genome-wide association study (GWAS) data from 3,568 pancreatic cancer cases and 3,363 controls of European Ancestry. RESULTS: Applying linkage disequilibrium- and minor allele frequency-stratified GREML (GREML-LDMS) method to imputed GWAS data, we estimated the overall heritability of pancreatic cancer to be 21.2% (SE = 4.8%). Across the functional groups (intronic, intergenic, coding, and regulatory variants), intronic variants account for most of the estimated heritability (12.4%). Previously identified GWAS loci explained 4.1% of the total phenotypic variation of pancreatic cancer. Mutations in hereditary pancreatic cancer susceptibility genes are present in 4% to 10% of patients with pancreatic cancer, yet our GREML-LDMS results suggested these regions explain only 0.4% of total phenotypic variance for pancreatic cancer. CONCLUSIONS: Although higher than previous studies, our estimated 21.2% overall heritability may still be downwardly biased due to the inherent limitation that the contribution of rare variants in genes with a substantive overall impact on disease are not captured when applying these commonly used methods to imputed GWAS data. IMPACT: Our work demonstrated the importance of rare and common variants in pancreatic cancer risk.

A New Fast Phasing Method Based On Haplotype Subtraction.

We developed a novel phasing approach, based solely on molecules and genotype frequency, that does not rely on inference of new alleles. We initiated the project because of errors that were detected in the phased 1000 Genomes Project data. The algorithm first combined identical genotypes into clusters and ranked them by descending frequency. Using alleles defined in homozygotes, it combined them to produce expected genotypes that were dismissed and subtracted them from remaining genotypes to define additional new putative alleles. Putative alleles had to be confirmed by identifying them in independent genotypes, and the process was iterated until all alleles were identified. The new approach was validated using single-molecule sequencing of eight loci, 145 (8 to 35 per locus) alleles were identified, and an average 98.2% (range, 95.0% to 99.9%) of 1000 genome individuals at these loci were explained. The accuracy of the new method was compared with that from PHASE and SHAPEIT2 to the experimentally determined genotypes based on single-molecule sequencing. Our method was comparable to PHASE and SHAPEIT2 in accuracy but was, on average, 14.6- and 10.8-fold faster, respectively.

Agnostic Pathway/Gene Set Analysis of Genome-Wide Association Data Identifies Associations for Pancreatic Cancer.

BACKGROUND: Genome-wide association studies (GWAS) identify associations of individual single-nucleotide polymorphisms (SNPs) with cancer risk but usually only explain a fraction of the inherited variability. Pathway analysis of genetic variants is a powerful tool to identify networks of susceptibility genes. METHODS: We conducted a large agnostic pathway-based meta-analysis of GWAS data using the summary-based adaptive rank truncated product method to identify gene sets and pathways associated with pancreatic ductal adenocarcinoma (PDAC) in 9040 cases and 12 496 controls. We performed expression quantitative trait loci (eQTL) analysis and functional annotation of the top SNPs in genes contributing to the top associated pathways and gene sets. All statistical tests were two-sided. RESULTS: We identified 14 pathways and gene sets associated with PDAC at a false discovery rate of less than 0.05. After Bonferroni correction (P ≤ 1.3 × 10-5), the strongest associations were detected in five pathways and gene sets, including maturity-onset diabetes of the young, regulation of beta-cell development, role of epidermal growth factor (EGF) receptor transactivation by G protein-coupled receptors in cardiac hypertrophy pathways, and the Nikolsky breast cancer chr17q11-q21 amplicon and Pujana ATM Pearson correlation coefficient (PCC) network gene sets. We identified and validated rs876493 and three correlating SNPs (PGAP3) and rs3124737 (CASP7) from the Pujana ATM PCC gene set as eQTLs in two normal derived pancreas tissue datasets. CONCLUSION: Our agnostic pathway and gene set analysis integrated with functional annotation and eQTL analysis provides insight into genes and pathways that may be biologically relevant for risk of PDAC, including those not previously identified.

Genome-wide meta-analysis identifies five new susceptibility loci for pancreatic cancer.

In 2020, 146,063 deaths due to pancreatic cancer are estimated to occur in Europe and the United States combined. To identify common susceptibility alleles, we performed the largest pancreatic cancer GWAS to date, including 9040 patients and 12,496 controls of European ancestry from the Pancreatic Cancer Cohort Consortium (PanScan) and the Pancreatic Cancer Case-Control Consortium (PanC4). Here, we find significant evidence of a novel association at rs78417682 (7p12/TNS3, P = 4.35 × 10-8). Replication of 10 promising signals in up to 2737 patients and 4752 controls from the PANcreatic Disease ReseArch (PANDoRA) consortium yields new genome-wide significant loci: rs13303010 at 1p36.33 (NOC2L, P = 8.36 × 10-14), rs2941471 at 8q21.11 (HNF4G, P = 6.60 × 10-10), rs4795218 at 17q12 (HNF1B, P = 1.32 × 10-8), and rs1517037 at 18q21.32 (GRP, P = 3.28 × 10-8). rs78417682 is not statistically significantly associated with pancreatic cancer in PANDoRA. Expression quantitative trait locus analysis in three independent pancreatic data sets provides molecular support of NOC2L as a pancreatic cancer susceptibility gene.

Haplotype Counting for Sensitive Chimerism Testing: Potential for Early Leukemia Relapse Detection.

Fields of forensics, transplantation, and paternity rely on human identity testing. Currently, this is accomplished through amplification of microsatellites followed by capillary electrophoresis. An alternative and theoretically better approach uses multiple single-nucleotide polymorphisms located within a small region of DNA, a method we initially developed using HLA-A and called haplotype counting. Herein, we validated seven additional polymorphic loci, sequenced a total of 45 individuals from three of the 1000 Genomes populations (15 from each), and determined the number of haplotypes, heterozygosity, and polymorphic information content for each locus. In addition, we developed a multiplex PCR that amplifies five of these loci simultaneously. Using this strategy with a small cohort of leukemic patients who underwent allogeneic bone marrow transplantation, we first attempted to define a threshold (0.26% recipient) by examining seven patients who tested all donor and did not relapse. Although this initial threshold will need to be confirmed in a larger cohort, we detected increased recipient DNA above this threshold 90 to 145 days earlier than microsatellite positivity, and 127 to 142 days before clinical relapse in four of eight patients (50%). Haplotype counting using these novel loci may be useful for ultrasensitive detection in fields such as bone marrow transplantation, solid organ transplant rejection, patient identification, and forensics.

Three new pancreatic cancer susceptibility signals identified on chromosomes 1q32.1, 5p15.33 and 8q24.21.

Genome-wide association studies (GWAS) have identified common pancreatic cancer susceptibility variants at 13 chromosomal loci in individuals of European descent. To identify new susceptibility variants, we performed imputation based on 1000 Genomes (1000G) Project data and association analysis using 5,107 case and 8,845 control subjects from 27 cohort and case-control studies that participated in the PanScan I-III GWAS. This analysis, in combination with a two-staged replication in an additional 6,076 case and 7,555 control subjects from the PANcreatic Disease ReseArch (PANDoRA) and Pancreatic Cancer Case-Control (PanC4) Consortia uncovered 3 new pancreatic cancer risk signals marked by single nucleotide polymorphisms (SNPs) rs2816938 at chromosome 1q32.1 (per allele odds ratio (OR) = 1.20, P = 4.88x10 -15), rs10094872 at 8q24.21 (OR = 1.15, P = 3.22x10 -9) and rs35226131 at 5p15.33 (OR = 0.71, P = 1.70x10 -8). These SNPs represent independent risk variants at previously identified pancreatic cancer risk loci on chr1q32.1 ( NR5A2), chr8q24.21 ( MYC) and chr5p15.33 ( CLPTM1L- TERT) as per analyses conditioned on previously reported susceptibility variants. We assessed expression of candidate genes at the three risk loci in histologically normal ( n = 10) and tumor ( n = 8) derived pancreatic tissue samples and observed a marked reduction of NR5A2 expression (chr1q32.1) in the tumors (fold change -7.6, P = 5.7x10 -8). This finding was validated in a second set of paired ( n = 20) histologically normal and tumor derived pancreatic tissue samples (average fold change for three NR5A2 isoforms -31.3 to -95.7, P = 7.5x10 -4-2.0x10 -3). Our study has identified new susceptibility variants independently conferring pancreatic cancer risk that merit functional follow-up to identify target genes and explain the underlying biology.