Colon cancer arises from the accumulation of multiple genetic and epigenetic alterations to normal colonic tissue. microRNAs (miRNAs) are small, non-coding regulatory RNAs that post-transcriptionally regulate gene expression. Differential miRNA expression in cancer versus normal tissue is a common event and may be pivotal for tumor onset and progression.

We conducted a large genome-wide association study (GWAS) of the immune responses to primary smallpox vaccination in a combined cohort of 1,653 subjects. We did not observe any polymorphisms associated with standard vaccine response outcomes (e.g., neutralizing antibody, T cell ELISPOT response, or T cell cytokine production); however, we did identify a cluster of SNPs on chromosome 5 (5q31.2) that were significantly associated (p-value: 1.3 x 10-12 - 1.5x10-36) with IFNα response to in vitro poxvirus stimulation. Examination of these SNPs led to the functional testing of rs1131769, a non-synonymous SNP in TMEM173 causing an Arg-to-His change at position 232 in the STING protein-a major regulator of innate immune responses to viral infections. Our findings demonstrate differences in the ability of the two STING variants to phosphorylate the downstream intermediates TBK1 and IRF3 in response to multiple STING ligands. Further downstream in the STING pathway, we observed significantly reduced expression of type I IFNs (including IFNα) and IFN-response genes in cells carrying the H232 variant. Subsequent molecular modeling of both alleles predicted altered ligand binding characteristics between the two variants, providing a potential mechanism underlying differences in inter-individual responses to poxvirus infection. Our data indicate that possession of the H232 variant may impair STING-mediated innate immunity to poxviruses. These results clarify prior studies evaluating functional effects of genetic variants in TMEM173 and provide novel data regarding genetic control of poxvirus immunity.

As patient derived xenograft (PDX) models are increasingly used for preclinical drug development, strategies to account for the nonhuman component of PDX RNA expression data are critical to its interpretation. A bioinformatics pipeline to separate donor tumor and mouse stroma transcriptome profiles was devised and tested. To examine the molecular fidelity of PDX versus donor tumors, we compared mRNA differences between paired PDX-donor tumors from nine ovarian cancer patients. 1,935 differentially expressed genes were identified between PDX and donor tumors. Over 90% (n = 1767) of these genes were down-regulated in PDX models and enriched in stroma-specific functions. Several protein kinases were also differentially expressed in PDX tumors, e.g. PDGFRA, PDGFRB and CSF1R. Upon in silico removal of these PDX-donor tumor differentially expressed genes, a stronger transcriptional resemblance between PDX-donor tumor pairs was seen (average correlation coefficient increases from 0.91 to 0.95). We devised and validated an effective bioinformatics strategy to separate mouse stroma expression from human tumor expression for PDX RNAseq. In addition, we showed most of the PDX-donor differentially expressed genes were implicated in stromal components. The molecular similarities and differences between PDX and donor tumors have implications in future therapeutic trial designs and treatment response evaluations using PDX models.

Long-term treatment outcomes for patients with high grade ovarian cancers have not changed despite innovations in therapies. There is no recommended assay for predicting patient response to second-line therapy, thus clinicians must make treatment decisions based on each individual patient. Patient-derived xenograft (PDX) tumors have been shown to predict drug sensitivity in ovarian cancer patients, but the time frame for intraperitoneal (IP) tumor generation, expansion, and drug screening is beyond that for tumor recurrence and platinum resistance to occur, thus results do not have clinical utility. We describe a drug sensitivity screening assay using a drug delivery microdevice implanted for 24 h in subcutaneous (SQ) ovarian PDX tumors to predict treatment outcomes in matched IP PDX tumors in a clinically relevant time frame. The SQ tumor response to local microdose drug exposure was found to be predictive of the growth of matched IP tumors after multi-week systemic therapy using significantly fewer animals (10 SQ vs 206 IP). Multiplexed immunofluorescence image analysis of phenotypic tumor response combined with a machine learning classifier could predict IP treatment outcomes against three second-line cytotoxic therapies with an average AUC of 0.91.

Resistance to platinum compounds is a major determinant of patient survival in high-grade serous ovarian cancer (HGSOC). To understand mechanisms of platinum resistance and identify potential therapeutic targets in resistant HGSOC, we generated a data resource composed of dynamic (±carboplatin) protein, post-translational modification, and RNA sequencing (RNA-seq) profiles from intra-patient cell line pairs derived from 3 HGSOC patients before and after acquiring platinum resistance. These profiles reveal extensive responses to carboplatin that differ between sensitive and resistant cells. Higher fatty acid oxidation (FAO) pathway expression is associated with platinum resistance, and both pharmacologic inhibition and CRISPR knockout of carnitine palmitoyltransferase 1A (CPT1A), which represents a rate limiting step of FAO, sensitize HGSOC cells to platinum. The results are further validated in patient-derived xenograft models, indicating that CPT1A is a candidate therapeutic target to overcome platinum resistance. All multiomic data can be queried via an intuitive gene-query user interface (https://sites.google.com/view/ptrc-cell-line).

Despite high levels of MMR-II usage in the US, mumps outbreaks continue to occur. Evidence suggests that mumps vaccine-induced humoral immunity wanes over time. Relatively few studies have examined cell-mediated immunity or reported on sex-based differences. To better understand sex-based differences in the immune response to mumps vaccine, we measured neutralizing antibody titers and mumps-specific cytokine/chemokine responses in a cohort of 748 adolescents and young adults after two doses of MMR vaccine. We observed significantly higher neutralizing antibody titers in females than in males (120.8 IU/mL, 98.7 IU/mL, p = 0.038) but significantly higher secretion levels of MIP-1α, MIP-1β, TNFα, IL-6, IFNγ, and IL-1β in males compared to females. These data demonstrate that sex influences mumps-specific humoral and cell-mediated immune response outcomes, a phenomenon that should be considered during efforts to improve vaccines and prevent future outbreaks.

One of the greatest barriers to curative treatment of neuroblastoma is its frequent metastatic outgrowth prior to diagnosis, especially in cases driven by amplification of the MYCN oncogene. However, only a limited number of regulatory proteins that contribute to this complex MYCN-mediated process have been elucidated. Here we show that the growth arrest-specific 7 (GAS7) gene, located at chromosome band 17p13.1, is preferentially deleted in high-risk MYCN-driven neuroblastoma. GAS7 expression was also suppressed in MYCN-amplified neuroblastoma lacking 17p deletion. GAS7 deficiency led to accelerated metastasis in both zebrafish and mammalian models of neuroblastoma with overexpression or amplification of MYCN. Analysis of expression profiles and the ultrastructure of zebrafish neuroblastoma tumors with MYCN overexpression identified that GAS7 deficiency led to (i) downregulation of genes involved in cell-cell interaction, (ii) loss of contact among tumor cells as critical determinants of accelerated metastasis, and (iii) increased levels of MYCN protein. These results provide the first genetic evidence that GAS7 depletion is a critical early step in the cascade of events culminating in neuroblastoma metastasis in the context of MYCN overexpression. SIGNIFICANCE: Heterozygous deletion or MYCN-mediated repression of GAS7 in neuroblastoma releases an important brake on tumor cell dispersion and migration to distant sites, providing a novel mechanism underlying tumor metastasis in MYCN-driven neuroblastoma.See related commentary by Menard, p. 2815.

Obesity and diabetes are major modifiable risk factors for pancreatic cancer. Interactions between genetic variants and diabetes/obesity have not previously been comprehensively investigated in pancreatic cancer at the genome-wide level.

Focal adhesion kinase (FAK) is an intensely studied non-receptor tyrosine kinase with roles in cancer and other common human diseases. Despite the large interest in FAK, the in vivo contribution of FAK auto-phosphorylation site tyrosine (Y) 397 to FAK function is incompletely understood. To study FAK Y397 in vivo we analyzed mice with 'non-phosphorylatable' Y-to-phenylalanine (F) and 'phospho-mimicking' Y-to-glutamate (E) mutations in the germline. We found that FAK Y397F mice die early during embryogenesis with abnormal angiogenesis like FAK kinase-dead mice. When Y397 is mutated to a glutamate mice survive beyond mid-gestation like mice where Y397 is lost by deletion of FAK exon 15. In culture, defects in proliferation, invasion and gene expression were more severe with the FAK Y397F than with the FAK Y397E mutation despite the inability of FAK Y397E to bind SRC. Conditional expression of FAK Y397F or Y397E in unchallenged avascular epidermis, however, resulted in no appreciable phenotype. We conclude that FAK Y397 is required for the highly dynamic tissue remodeling during development but dispensable for normal homeostasis of avascular epidermis. In contrast to the Y397F mutation, FAK Y397E retains sufficient biological activity to allow for development beyond mid-gestation.

The lack of standardization of the currently used commercial anti-rubella IgG antibody assays leads to frequent misinterpretation of results for samples with low/equivocal antibody concentration. The use of alternative approaches in rubella serology could add new information leading to a fuller understanding of rubella protective immunity and neutralizing antibody response after vaccination.

Although influenza causes significant morbidity and mortality in the elderly, the factors underlying the reduced vaccine immunogenicity and efficacy in this age group are not completely understood. Age and immunosenescence factors, and their impact on humoral immunity after influenza vaccination, are of growing interest for the development of better vaccines for the elderly.

Pancreatic cancer is the fourth leading cause of cancer death in the developed world. Both inherited high-penetrance mutations in BRCA2 (ref. 2), ATM, PALB2 (ref. 4), BRCA1 (ref. 5), STK11 (ref. 6), CDKN2A and mismatch-repair genes and low-penetrance loci are associated with increased risk. To identify new risk loci, we performed a genome-wide association study on 9,925 pancreatic cancer cases and 11,569 controls, including 4,164 newly genotyped cases and 3,792 controls in 9 studies from North America, Central Europe and Australia. We identified three newly associated regions: 17q25.1 (LINC00673, rs11655237, odds ratio (OR) = 1.26, 95% confidence interval (CI) = 1.19-1.34, P = 1.42 × 10(-14)), 7p13 (SUGCT, rs17688601, OR = 0.88, 95% CI = 0.84-0.92, P = 1.41 × 10(-8)) and 3q29 (TP63, rs9854771, OR = 0.89, 95% CI = 0.85-0.93, P = 2.35 × 10(-8)). We detected significant association at 2p13.3 (ETAA1, rs1486134, OR = 1.14, 95% CI = 1.09-1.19, P = 3.36 × 10(-9)), a region with previous suggestive evidence in Han Chinese. We replicated previously reported associations at 9q34.2 (ABO), 13q22.1 (KLF5), 5p15.33 (TERT and CLPTM1), 13q12.2 (PDX1), 1q32.1 (NR5A2), 7q32.3 (LINC-PINT), 16q23.1 (BCAR1) and 22q12.1 (ZNRF3). Our study identifies new loci associated with pancreatic cancer risk.

Gene set analysis (GSA) has been used for analysis of microarray data to aid the interpretation and to increase statistical power. With the advent of next-generation sequencing, the use of GSA is even more relevant, as studies are often conducted on a small number of samples. We propose the use of soft truncation thresholding and the Gamma Method (GM) to determine significant gene set (GS), where a generalized linear model is used to assess per-gene significance. The approach was compared to other methods using an extensive simulation study and RNA-seq data from smallpox vaccine study. The GM was found to outperform other proposed methods. Application of the GM to the smallpox vaccine study found the GSs to be moderately associated with response, including focal adhesion (p = 0.04) and extracellular matrix receptor interaction (p = 0.05). The application of GSA to RNA-seq data will provide new insights into the genomic basis of complex traits.

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.

Immune responses to current rubella vaccines demonstrate significant inter-individual variability. We performed mRNA-Seq profiling on PBMCs from high and low antibody responders to rubella vaccination to delineate transcriptional differences upon viral stimulation. Generalized linear models were used to assess the per gene fold change (FC) for stimulated versus unstimulated samples or the interaction between outcome and stimulation. Model results were evaluated by both FC and p-value. Pathway analysis and self-contained gene set tests were performed for assessment of gene group effects. Of 17,566 detected genes, we identified 1,080 highly significant differentially expressed genes upon viral stimulation (p<1.00E(-15), FDR<1.00E(-14)), including various immune function and inflammation-related genes, genes involved in cell signaling, cell regulation and transcription, and genes with unknown function. Analysis by immune outcome and stimulation status identified 27 genes (p≤0.0006 and FDR≤0.30) that responded differently to viral stimulation in high vs. low antibody responders, including major histocompatibility complex (MHC) class I genes (HLA-A, HLA-B and B2M with p = 0.0001, p = 0.0005 and p = 0.0002, respectively), and two genes related to innate immunity and inflammation (EMR3 and MEFV with p = 1.46E(-08) and p = 0.0004, respectively). Pathway and gene set analysis also revealed transcriptional differences in antigen presentation and innate/inflammatory gene sets and pathways between high and low responders. Using mRNA-Seq genome-wide transcriptional profiling, we identified antigen presentation and innate/inflammatory genes that may assist in explaining rubella vaccine-induced immune response variations. Such information may provide new scientific insights into vaccine-induced immunity useful in rational vaccine development and immune response monitoring.

Interferon-induced antiviral genes are crucial players in innate antiviral defense and potential determinants of immune response heterogeneity. We selected 114 candidate single-nucleotide polymorphisms (SNPs) from 12 antiviral genes using an LD tagSNP selection approach and genotyped them in a cohort of 738 school children immunized with two doses of rubella vaccine. Associations between SNPs/haplotypes and rubella virus-specific immune measures were assessed using linear regression methodologies. We identified 23 significant associations (p < 0.05) between polymorphisms within the 2'-5'-oligoadenylate synthetase (OAS) gene cluster, and rubella virus-specific IL-2, IL-10, IL-6 secretion, and antibody levels. The minor allele variants of three OAS1 SNPs (rs3741981/Ser162Gly, rs1051042/Thr361Arg, rs2660), located in a linkage disequilibrium block of functional importance, were significantly associated with an increase in rubella virus-specific IL-2/T(h)1 response (p

The goal of annual influenza vaccination is to reduce mortality and morbidity associated with this disease through the generation of protective immune responses. The objective of the current study was to examine markers of immunosenescence and identify immunosenescence-related differences in gene expression, gene regulation, cytokine secretion, and immunologic changes in an older study population receiving seasonal influenza A/H1N1 vaccination. Surprisingly, prior studies in this cohort revealed weak correlations between immunosenescence markers and humoral immune response to vaccination. In this report, we further examined the relationship of each immunosenescence marker (age, T cell receptor excision circle frequency, telomerase expression, percentage of CD28- CD4+ T cells, percentage of CD28- CD8+ T cells, and the CD4/CD8 T cell ratio) with additional markers of immune response (serum cytokine and chemokine expression) and measures of gene expression and/or regulation. Many of the immunosenescence markers indeed correlated with distinct sets of individual DNA methylation sites, miRNA expression levels, mRNA expression levels, serum cytokines, and leukocyte subsets. However, when the individual immunosenescence markers were grouped by pathways or functional terms, several shared biological functions were identified: antigen processing and presentation pathways, MAPK, mTOR, TCR, BCR, and calcium signaling pathways, as well as key cellular metabolic, proliferation and survival activities. Furthermore, the percent of CD4+ and/or CD8+ T cells lacking CD28 expression also correlated with miRNAs regulating clusters of genes known to be involved in viral infection. Integrated (DNA methylation, mRNA, miRNA, and protein levels) network biology analysis of immunosenescence-related pathways and genesets identified both known pathways (e.g., chemokine signaling, CTL, and NK cell activity), as well as a gene expression module not previously annotated with a known function. These results may improve our ability to predict immune responses to influenza and aid in new vaccine development, and highlight the need for additional studies to better define and characterize immunosenescence.

Human antibody response to measles vaccine is highly variable in the population. Host genes contribute to inter-individual antibody response variation. The killer cell immunoglobulin-like receptors (KIR) are recognized to interact with HLA molecules and possibly influence humoral immune response to viral antigens. To expand on and improve our previous work with HLA genes, and to explore the genetic contribution of KIR genes to the inter-individual variability in measles vaccine-induced antibody responses, we performed a large population-based study in 2,506 healthy immunized subjects (ages 11 to 41 years) to identify HLA and KIR associations with measles vaccine-induced neutralizing antibodies. After correcting for the large number of statistical tests of allele effects on measles-specific neutralizing antibody titers, no statistically significant associations were found for either HLA or KIR loci. However, suggestive associations worthy of follow-up in other cohorts include B*57:01, DQB1*06:02, and DRB1*15:05 alleles. Specifically, the B*57:01 allele (1,040 mIU/mL; p = 0.0002) was suggestive of an association with lower measles antibody titer. In contrast, the DQB1*06:02 (1,349 mIU/mL; p = 0.0004) and DRB1*15:05 (2,547 mIU/mL; p = 0.0004) alleles were suggestive of an association with higher measles antibodies. Notably, the associations with KIR genotypes were strongly nonsignificant, suggesting that KIR loci in terms of copy number and haplotypes are not likely to play a major role in antibody response to measles vaccination. These findings refine our knowledge of the role of HLA and KIR alleles in measles vaccine-induced immunity.

MicroRNAs are important mediators of post-transcriptional regulation of gene expression through RNA degradation and translational repression, and are emerging biomarkers of immune system activation/response after vaccination.

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.