Comparing clinical and genomic features based on the tumor location in patients with resected pancreatic cancer.

BACKGROUND: Pancreatic cancer is anatomically divided into pancreatic head and body/tail cancers, and some studies have reported differences in prognosis. However, whether this discrepancy is induced from the difference of tumor biology is hotly debated. Therefore, we aimed to evaluate the differences in clinical outcomes and tumor biology depending on the tumor location. METHODS: In this retrospective cohort study, we identified 800 patients with pancreatic ductal adenocarcinoma who had undergone upfront curative-intent surgery. Cox regression analysis was performed to explore the prognostic impact of the tumor location. Among them, 153 patients with sufficient tumor tissue and blood samples who provided informed consent for next-generation sequencing were selected as the cohort for genomic analysis. RESULTS: Out of the 800 patients, 500 (62.5%) had pancreatic head cancer, and 300 (37.5%) had body/tail cancer. Tumor location in the body/tail of the pancreas was not identified as a significant predictor of survival outcomes compared to that in the head in multivariate analysis (hazard ratio, 0.94; 95% confidence interval, 0.77-1.14; P = 0.511). Additionally, in the genomic analyses of 153 patients, there were no significant differences in mutational landscapes, distribution of subtypes based on transcriptomic profiling, and estimated infiltration levels of various immune cells between pancreatic head and body/tail cancers. CONCLUSIONS: We could not find differences in prognosis and tumor biology depending on tumor location in pancreatic ductal adenocarcinoma. Discrepancies in prognosis may represent a combination of lead time, selection bias, and clinical differences, including the surgical burden between tumor sites.

Context-dependent genomic locus effects on antibody production in recombinant Chinese hamster ovary cells generated through random integration.

High-yield production of therapeutic protein using Chinese hamster ovary (CHO) cells requires stable cell line development (CLD). CLD typically uses random integration of transgenes; however, this results in clonal variation and subsequent laborious clone screening. Therefore, site-specific integration of a protein expression cassette into a desired chromosomal locus showing high transcriptional activity and stability, referred to as a hot spot, is emerging. Although positional effects are important for therapeutic protein expression, the sequence-specific mechanisms by which hotspots work are not well understood. In this study, we performed whole-genome sequencing (WGS) to locate randomly inserted vectors in the genome of recombinant CHO cells expressing high levels of monoclonal antibodies (mAbs) and experimentally validated these locations and vector compositions. The integration site was characterized by active histone marks and potential enhancer activities, and clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) mediated indel mutations in the region upstream of the integration site led to a significant reduction in specific antibody productivity by up to 30%. Notably, the integration site and its core region did not function equivalently outside the native genomic context, showing a minimal effect on the increase in exogenous protein expression in the host cell line. We also observed a superior production capacity of the mAb expressing cell line compared to that of the host cell line. Collectively, this study demonstrates that developing recombinant CHO cell lines to produce therapeutic proteins at high levels requires a balance of factors including transgene configuration, genomic locus landscape, and host cell properties.

Genetic assessment of pathogenic germline alterations in lysosomal genes among Asian patients with pancreatic ductal adenocarcinoma.

BACKGROUND: Lysosomes are closely linked to autophagic activity, which plays a vital role in pancreatic ductal adenocarcinoma (PDAC) biology. The survival of PDAC patients is still poor, and the identification of novel genetic factors for prognosis and treatment is highly required to prevent PDAC-related deaths. This study investigated the germline variants related to lysosomal dysfunction in patients with PDAC and to analyze whether they contribute to the development of PDAC. METHODS: The germline putative pathogenic variants (PPV) in genes involved in lysosomal storage disease (LSD) was compared between patients with PDAC (n = 418) and healthy controls (n = 845) using targeted panel and whole-exome sequencing. Furthermore, pancreatic organoids from wild-type and KrasG12D mice were used to evaluate the effect of lysosomal dysfunction on PDAC development. RNA sequencing (RNA-seq) analysis was performed with established PDAC patient-derived organoids (PDOs) according to the PPV status. RESULTS: The PPV in LSD-related genes was higher in patients with PDAC than in healthy controls (8.13 vs. 4.26%, Log2 OR = 1.65, P = 3.08 × 10-3). The PPV carriers of LSD-related genes with PDAC were significantly younger than the non-carriers (mean age 61.5 vs. 65.3 years, P = 0.031). We further studied a variant of the lysosomal enzyme, galactosylceramidase (GALC), which was the most frequently detected LSD variant in our cohort. Autophagolysosomal activity was hampered when GALC was downregulated, which was accompanied by paradoxically elevated autophagic flux. Furthermore, the number of proliferating Ki-67+ cells increased significantly in pancreatic organoids derived from Galc knockout KrasG12D mice. Moreover, GALC PPV carriers tended to show drug resistance in both PDAC cell line and PDAC PDO, and RNA-seq analysis revealed that various metabolism and gene repair pathways were upregulated in PDAC PDOs harboring a GALC variant. CONCLUSIONS: Genetically defined lysosomal dysfunction is frequently observed in patients with young-onset PDAC. This might contribute to PDAC development by altering metabolism and impairing autophagolysosomal activity, which could be potentially implicated in therapeutic applications for PDAC.

Assessment of attenuation of varicella-zoster virus vaccines based on genomic comparison.

Live attenuated varicella-zoster virus (VZV) vaccines are used to prevent chickenpox and shingles. Single nucleotide polymorphisms (SNPs) that occur during the attenuation of parental strains are critical indicators of vaccine safety. To assess the attenuation of commercial VZV vaccines, genetic variants were comprehensively examined through high-throughput sequencing of viral DNA isolated from four VZV vaccines (Barycela, VarilRix, VariVax, and SKY Varicella). Whole-genome comparison of the four vaccines with the wild-type strain (Dumas) revealed that the sequences are highly conserved on a genome-wide scale. Among the 196 common variants across the four vaccines, 195 were already present in the genome of the parental strain (pOka), indicating that the variants occurred during the generation of the parental strain from the Dumas strain. Compared to the pOka genome, the vaccines exhibited distinct variant frequencies on a genome-wide and within an attenuation-related open reading frame. In particular, attenuation-associated 42 SNPs showed that Barycela, VarilRix, VariVax, and SKY Varicella are in ascending order regarding similarity with pOka-like genotypes, which in turn, might provide genomic evidence for the levels of attenuation. Finally, the phylogenetic network analysis demonstrated that genetic distances from the parental strain correlated with the attenuation levels of the vaccines.

Proteogenomic landscape of human pancreatic ductal adenocarcinoma in an Asian population reveals tumor cell-enriched and immune-rich subtypes.

We report a proteogenomic analysis of pancreatic ductal adenocarcinoma (PDAC). Mutation-phosphorylation correlations identified signaling pathways associated with somatic mutations in significantly mutated genes. Messenger RNA-protein abundance correlations revealed potential prognostic biomarkers correlated with patient survival. Integrated clustering of mRNA, protein and phosphorylation data identified six PDAC subtypes. Cellular pathways represented by mRNA and protein signatures, defining the subtypes and compositions of cell types in the subtypes, characterized them as classical progenitor (TS1), squamous (TS2-4), immunogenic progenitor (IS1) and exocrine-like (IS2) subtypes. Compared with the mRNA data, protein and phosphorylation data further classified the squamous subtypes into activated stroma-enriched (TS2), invasive (TS3) and invasive-proliferative (TS4) squamous subtypes. Orthotopic mouse PDAC models revealed a higher number of pro-tumorigenic immune cells in TS4, inhibiting T cell proliferation. Our proteogenomic analysis provides significantly mutated genes/biomarkers, cellular pathways and cell types as potential therapeutic targets to improve stratification of patients with PDAC.

Fast Growing Furious Races for Targeting Fibroblast Growth Factor Receptors.

Targeting fibroblast growth factor receptors (FGFRs) has been slow compared to other targeted cancer therapies for receptor tyrosine kinases, such as epidermal growth factor receptors. The low efficacy and variable response have limited the growth of FGFR inhibitors in clinical use. Nevertheless, recent systematic and genomic approaches have identified the biological conditions for effectively targeting FGFRs and can accelerate the development of targeted drugs. Under clinical and preclinical trials, the inhibitors started fast growing furious races to target FGFRs. Finally, FGFRs will be more actionable and targetable with more precise and effective drugs at the end of the race, passing the finish line.

Challenges in genomic analysis of model systems and primary tumors of pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is characterized by aggressive tumor behavior and poor prognosis. Recent next-generation sequencing (NGS)-based genomic studies have provided novel treatment modes for pancreatic cancer via the identification of cancer driver variants and molecular subtypes in PDAC. Genome-wide approaches have been extended to model systems such as patient-derived xenografts (PDXs), organoids, and cell lines for pre-clinical purposes. However, the genomic characteristics vary in the model systems, which is mainly attributed to the clonal evolution of cancer cells during their construction and culture. Moreover, fundamental limitations such as low tumor cellularity and the complex tumor microenvironment of PDAC hinder the confirmation of genomic features in the primary tumor and model systems. The occurrence of these phenomena and their associated complexities may lead to false insights into the understanding of mechanisms and dynamics in tumor tissues of patients. In this review, we describe various model systems and discuss differences in the results based on genomics and transcriptomics between primary tumors and model systems. Finally, we introduce practical strategies to improve the accuracy of genomic analysis of primary tissues and model systems.

PHF20 is crucial for epigenetic control of starvation-induced autophagy through enhancer activation.

Autophagy is a catabolic pathway that maintains cellular homeostasis under various stress conditions, including conditions of nutrient deprivation. To elevate autophagic flux to a sufficient level under stress conditions, transcriptional activation of autophagy genes occurs to replenish autophagy components. Thus, the transcriptional and epigenetic control of the genes regulating autophagy is essential for cellular homeostasis. Here, we applied integrated transcriptomic and epigenomic profiling to reveal the roles of plant homeodomain finger protein 20 (PHF20), which is an epigenetic reader possessing methyl binding activity, in controlling the expression of autophagy genes. Phf20 deficiency led to impaired autophagic flux and autophagy gene expression under glucose starvation. Interestingly, the genome-wide characterization of chromatin states by Assay for Transposase-Accessible Chromatin (ATAC)-sequencing revealed that the PHF20-dependent chromatin remodelling occurs in enhancers that are co-occupied by dimethylated lysine 36 on histone H3 (H3K36me2). Importantly, the recognition of H3K36me2 by PHF20 was found to be highly correlated with increased levels of H3K4me1/2 at the enhancer regions. Collectively, these results indicate that PHF20 regulates autophagy genes through enhancer activation via H3K36me2 recognition as an epigenetic reader. Our findings emphasize the importance of nuclear events in the regulation of autophagy.

Substrate-binding loop interactions with pseudouridine trigger conformational changes that promote catalytic efficiency of pseudouridine kinase PUKI.

Pseudouridine, one major RNA modification, is catabolized into uracil and ribose-5'-phosphate by two sequential enzymatic reactions. In the first step, pseudouridine kinase (PUKI) phosphorylates pseudouridine to pseudouridine 5'-monophosphate. High-fidelity catalysis of pseudouridine by PUKI prevents possible disturbance of in vivo pyrimidine homeostasis. However, the molecular basis of how PUKI selectively phosphorylates pseudouridine over uridine with >100-fold greater efficiency despite minor differences in their Km values has not been elucidated. To investigate this selectivity, in this study we determined the structures of PUKI from Escherichia coli strain B (EcPUKI) in various ligation states. The structure of EcPUKI was determined to be similar to PUKI from Arabidopsis thaliana, including an α/ß core domain and ß-stranded small domain, with dimerization occurring via the ß-stranded small domain. In a binary complex, we show that Ser30 in the substrate-binding loop of the small domain mediates interactions with the hallmark N1 atom of pseudouridine nucleobase, causing conformational changes in its quaternary structure. Kinetic and fluorescence spectroscopic analyses also showed that the Ser30-mediated interaction is a prerequisite for conformational changes and subsequent catalysis by EcPUKI. Furthermore, S30A mutation or EcPUKI complexed with other nucleosides homologous to pseudouridine but lacking the pseudouridine-specific N1 atom did not induce such conformational changes, demonstrating the catalytic significance of the proposed Ser30-mediated interaction. These analyses provide structural and functional evidence for a pseudouridine-dependent conformational change of EcPUKI and its functional linkage to catalysis.

Influenza vaccination in the elderly boosts antibodies against conserved viral proteins and egg-produced glycans.

Seasonal influenza vaccination elicits a diminished adaptive immune response in the elderly, and the mechanisms of immunosenescence are not fully understood. Using Ig-Seq, we found a marked increase with age in the prevalence of cross-reactive (CR) serum antibodies that recognize both the H1N1 (vaccine-H1) and H3N2 (vaccine-H3) components of an egg-produced split influenza vaccine. CR antibodies accounted for 73% ± 18% of the serum vaccine responses in a cohort of elderly donors, 65% ± 15% in late middle-aged donors, and only 13% ± 5% in persons under 35 years of age. The antibody response to non-HA antigens was boosted by vaccination. Recombinant expression of 19 vaccine-H1+H3 CR serum monoclonal antibodies (s-mAbs) revealed that they predominantly bound to non-HA influenza proteins. A sizable fraction of vaccine-H1+H3 CR s-mAbs recognized with high affinity the sulfated glycans, in particular sulfated type 2 N-acetyllactosamine (Galß1-4GalNAcß), which is found on egg-produced proteins and thus unlikely to contribute to protection against influenza infection in humans. Antibodies against sulfated glycans in egg-produced vaccine had been identified in animals but were not previously characterized in humans. Collectively, our results provide a quantitative basis for how repeated exposure to split influenza vaccine correlates with unintended focusing of serum antibody responses to non-HA antigens that may result in suboptimal immunity against influenza.

Diffusion-weighted MR imaging in pancreatic ductal adenocarcinoma: prediction of next-generation sequencing-based tumor cellularity and prognosis after surgical resection.

PURPOSE: To identify features on preoperative MR imaging with diffusion-weighted imaging (DWI) for predicting next-generation sequencing (NGS)-based tumor cellularity and patient outcome after surgical resection of pancreatic ductal adenocarcinoma (PDAC). METHODS: This retrospective study included 105 patients with surgically resected PDAC who underwent preoperative MR imaging with DWI. Tumor cellularity was measured using molecular techniques and bioinformatics methods. Clinico-pathologic findings including tumor T stage for predicting disease-free survival (DFS) and overall survival (OS) were identified using Cox proportional hazards model. Important MR imaging findings including apparent diffusion coefficient (ADC) value of PDAC and modified ADC value (the ratio of the ADC value of PDAC to the ADC value of the spleen) for predicting higher tumor cellularity (≥ 30%) and poor prognosis were also identified. RESULTS: The median DFS and OS were 12.0 months [95% confidence interval (CI), 8.0-17.0] and 22.0 months (95% CI, 18.0-29.0), respectively. Higher T stage (T3/4) [hazard ratio (HR), 7.720, (95% CI 1.072, 55.612); p = 0.048] and higher tumor cellularity [HR, 1.599 (95% CI, 1.003-2.548); p = 0.048] were significantly associated with worse DFS. Among MR imaging features, the modified ADC value was significantly associated with tumor cellularity [odds ratio, 0.068 (95% CI, 0.012-0.372); p = 0.002], and PDAC with lower modified ADC value [≤ 1.40 (cutoff value)] showed significantly shorter median DFS than PDAC with higher modified ADC value [8 months (95% CI, 4-12) vs. 16 months (95% CI, 10-29); HR, 1.713 (95% CI, 1.073-2.735), log-rank p = 0.024]. CONCLUSION: Higher NGS-based tumor cellularity may be a negative prognostic factor in pancreatic cancer after resection, and modified ADC value derived from DWI could be helpful in predicting tumor cellularity and patient surgical outcome with regard to recurrence.

Telomeres reforged with non-telomeric sequences in mouse embryonic stem cells.

Telomeres are part of a highly refined system for maintaining the stability of linear chromosomes. Most telomeres rely on simple repetitive sequences and telomerase enzymes to protect chromosomal ends; however, in some species or telomerase-defective situations, an alternative lengthening of telomeres (ALT) mechanism is used. ALT mainly utilises recombination-based replication mechanisms and the constituents of ALT-based telomeres vary depending on models. Here we show that mouse telomeres can exploit non-telomeric, unique sequences in addition to telomeric repeats. We establish that a specific subtelomeric element, the mouse template for ALT (mTALT), is used for repairing telomeric DNA damage as well as for composing portions of telomeres in ALT-dependent mouse embryonic stem cells. Epigenomic and proteomic analyses before and after ALT activation reveal a high level of non-coding mTALT transcripts despite the heterochromatic nature of mTALT-based telomeres. After ALT activation, the increased HMGN1, a non-histone chromosomal protein, contributes to the maintenance of telomere stability by regulating telomeric transcription. These findings provide a molecular basis to study the evolution of new structures in telomeres.

The chromatin-binding protein PHF6 functions as an E3 ubiquitin ligase of H2BK120 via H2BK12Ac recognition for activation of trophectodermal genes.

Epigenetic regulation is important for establishing lineage-specific gene expression during early development. Although signaling pathways have been well-studied for regulation of trophectoderm reprogramming, epigenetic regulation of trophectodermal genes with histone modification dynamics have been poorly understood. Here, we identify that plant homeodomain finger protein 6 (PHF6) is a key epigenetic regulator for activation of trophectodermal genes using RNA-sequencing and ChIP assays. PHF6 acts as an E3 ubiquitin ligase for ubiquitination of H2BK120 (H2BK120ub) via its extended plant homeodomain 1 (PHD1), while the extended PHD2 of PHF6 recognizes acetylation of H2BK12 (H2BK12Ac). Intriguingly, the recognition of H2BK12Ac by PHF6 is important for exerting its E3 ubiquitin ligase activity for H2BK120ub. Together, our data provide evidence that PHF6 is crucial for epigenetic regulation of trophectodermal gene expression by linking H2BK12Ac to H2BK120ub modification.

Combinatorial Inhibition of Cell Surface Receptors Using Dual Aptamer-Functionalized Nanoconstructs for Cancer Treatment.

Membrane receptors overexpressed in diseased states are considered novel therapeutic targets. However, the single targeting approach faces several fundamental issues, such as poor efficacy, resistance, and toxicity. Here, we report a dual-targeting strategy to enhance anti-cancer efficacy via synergistic proximity interactions between therapeutics and two receptor proteins. Importantly, we report the first finding of an interaction between c-Met and nucleolin and demonstrate the therapeutic value of targeting the interaction between them. Bispecific nanocarriers densely grafted with anti-c-Met and -nucleolin aptamer increased the local concentration of aptamers at the target sites, in addition to inducing target receptor clustering. It was also demonstrated that the simultaneous targeting of c-Met and nucleolin inhibited the cellular functions of the receptors and increased anti-cancer efficacy by altering the cell cycle. Our findings pave the way for the development of an effective combinatorial treatment based on nanoconstruct-mediated interaction between receptors.

iNrich, Rapid and Robust Method to Enrich N-Terminal Proteome in a Highly Multiplexed Platform.

The field of terminal proteomics is limited in that it is optimized for large-scale analysis via multistep processes involving liquid chromatography. Here, we present an integrated N-terminal peptide enrichment method (iNrich) that can handle as little as 25 µg of cell lysate via a single-stage encapsulated solid-phase extraction column. iNrich enables simple, rapid, and reproducible sample processing, treatment of a wide range of protein amounts (25 µg ∼ 1 mg), multiplexed parallel sample preparation, and in-stage sample prefractionation using a mixed-anion-exchange filter. We identified ∼5000 N-terminal peptides (Nt-peptides) from only 100 µg of human cell lysate including Nt-formyl peptides. Multiplexed sample preparation facilitated quantitative and robust enrichment of N-terminome with dozens of samples simultaneously. We further developed the method to incorporate isobaric tags such as a tandem mass tag (TMT) and used it to discover novel peptides during ER stress analysis. The iNrich facilitated high-throughput N-terminomics and degradomics at a low cost using commercially available reagents and apparatus, without requiring arduous procedures.

Streamlined selection of cancer antigens for vaccine development through integrative multi-omics and high-content cell imaging.

Identification of tumor antigens that induce cytotoxic T lymphocytes (CTLs) is crucial for cancer-vaccine development. Despite their predictive ability, current algorithmic approaches and human leukocyte antigen (HLA)-peptidomic analysis allow limited selectivity. Here, we optimized a method to rapidly screen and identify highly immunogenic epitopes that trigger CTL responses. We used a combined application of this method involving immune-specific signature analysis and HLA-associated peptidomics using samples from six patients with triple-negative breast cancer (TNBC) in order to select immunogenic HLA epitopes for in vitro testing. Additionally, we applied high-throughput imaging at the single-cell level in order to confirm the immunoreactivity of the selected peptides. The results indicated that this method enabled identification of promising CTL peptides capable of inducing antitumor immunity. This platform combining high-resolution computational analysis, HLA-peptidomics, and high-throughput immunogenicity testing allowed rapid and robust identification of highly immunogenic epitopes and represents a powerful technique for cancer-vaccine development.

Determinants governing T cell receptor α/ß-chain pairing in repertoire formation of identical twins.

The T cell repertoire in each individual includes T cell receptors (TCRs) of enormous sequence diversity through the pairing of diverse TCR α- and ß-chains, each generated by somatic recombination of paralogous gene segments. Whether the TCR repertoire contributes to susceptibility to infectious or autoimmune diseases in concert with disease-associated major histocompatibility complex (MHC) polymorphisms is unknown. Due to a lack in high-throughput technologies to sequence TCR α-ß pairs, current studies on whether the TCR repertoire is shaped by host genetics have so far relied only on single-chain analysis. Using a high-throughput single T cell sequencing technology, we obtained the largest paired TCRαß dataset so far, comprising 965,523 clonotypes from 15 healthy individuals including 6 monozygotic twin pairs. Public TCR α- and, to a lesser extent, TCR ß-chain sequences were common in all individuals. In contrast, sharing of entirely identical TCRαß amino acid sequences was very infrequent in unrelated individuals, but highly increased in twins, in particular in CD4 memory T cells. Based on nucleotide sequence identity, a subset of these shared clonotypes appeared to be the progeny of T cells that had been generated during fetal development and had persisted for more than 50 y. Additional shared TCRαß in twins were encoded by different nucleotide sequences, implying that genetic determinants impose structural constraints on thymic selection that favor the selection of TCR α-ß pairs with entire sequence identities.

The histone variant H2A.Z in yeast is almost exclusively incorporated into the +1 nucleosome in the direction of transcription.

Transcription start sites (TSS) in eukaryotes are characterized by a nucleosome-depleted region (NDR), which appears to be flanked upstream and downstream by strongly positioned nucleosomes incorporating the histone variant H2A.Z. H2A.Z associates with both active and repressed TSS and is important for priming genes for rapid transcriptional activation. However, the determinants of H2A.Z occupancy at specific nucleosomes and its relationship to transcription initiation remain unclear. To further elucidate the specificity of H2A.Z, we determined its genomic localization at single nucleosome resolution, as well as the localization of its chromatin remodelers Swr1 and Ino80. By analyzing H2A.Z occupancy in conjunction with RNA expression data that captures promoter-derived antisense initiation, we find that H2A.Z's bimodal incorporation on either side of the NDR is not a general feature of TSS, but is specifically a marker for bidirectional transcription, such that the upstream flanking -1 H2A.Z-containing nucleosome is more appropriately considered as a +1 H2A.Z nucleosome for antisense transcription. The localization of H2A.Z almost exclusively at the +1 nucleosome suggests that a transcription-initiation dependent process could contribute to its specific incorporation.