Single-cell landscape of idiopathic multicentric Castleman disease in identical twins.

ABSTRACT: Idiopathic multicentric Castleman disease (iMCD) is a rare cytokine-driven disorder characterized by systemic inflammation, generalized lymphadenopathy, and organ dysfunction. Here, we present an unusual occurrence of iMCD in identical twins and examined the immune milieu within the affected lymphoid organs and the host circulation using multiomic high-dimensional profiling. Using spatial enhanced resolution omics sequencing (Stereo-seq) transcriptomic profiling, we performed unsupervised spatially constrained clustering to identify different anatomic structures, mapping the follicles and interfollicular regions. After a cell segmentation approach, interleukin 6 (IL-6) pathway genes significantly colocalized with endothelial cells and fibroblastic reticular cells, confirming observations using a single-cell sequencing approach (10× Chromium). Furthermore, single-cell sequencing of peripheral blood mononuclear cells revealed an "inflammatory" peripheral monocytosis enriched for the expression of S100A family genes in both twins. In summary, we provided evidence of the putative cell-of-origin of IL-6 signals in iMCD and described a distinct monocytic host immune response phenotype through a unique identical twin model.

Next-generation integrated sequencing identifies poor prognostic factors in patients with MYD88-mutated chronic lymphocytic leukemia in Taiwan.

INTRODUCTION: Chronic lymphocytic leukemia (CLL) is the most common type of leukemia in the western countries and is very rare in Asia. METHODS: Peripheral blood or bone marrow mononuclear cells obtained at initial diagnosis from 215 patients with CLL were analyzed by using next-generation sequencing to investigate the ethnic differences in genetic abnormalities. RESULTS: Whole-genome sequencing and whole-exome sequencing analyses on 30 cases showed that 9 genes, including IGLL5, MYD88, TCHH, DSCAM, AXDND1, BICRA, KMT2D, MYT1L, RBM43 were more frequently mutated in our Taiwanese cohort compared with those of the western cohorts. IGLL5, MYD88, and KMT2D genes were further analyzed by targeted sequencing in another 185 CLL patients, unravelling frequencies of 29.3%, 20.9%, and 15.0%, respectively. The most frequent positional mutation of MYD88 was V217F (26/45, 57.8%), followed by L265P (9/45, 20.0%). MYD88 mutations were significantly associated with IGLL5 mutations (P = 0.0004), mutated IGHV (P < 0.0001) and 13q deletion (P = 0.0164). CLL patients with co-occurrence of MYD88 mutations with KMT2D or/and IGLL5 mutations were associated with a significantly inferior survival compared to those with MYD88 mutation alone (not reached vs. 131.8 months, P = 0.007). In multivariate analysis, MYD88 mutation without KMT2D or IGLL5 mutations was an independently favorable predictor. CONCLUSIONS: IGLL5, MYD88, and KMT2D mutations were enriched in Taiwanese CLL, and co-occurrence of MYD88 mutations with KMT2D or/and IGLL5 mutations was associated with a poorer prognosis.

Super-enhancer-driven TOX2 mediates oncogenesis in Natural Killer/T Cell Lymphoma.

BACKGROUND: Extranodal natural killer/T-cell lymphoma (NKTL) is an aggressive type of non-Hodgkin lymphoma with dismal outcome. A better understanding of disease biology and key oncogenic process is necessary for the development of targeted therapy. Super-enhancers (SEs) have been shown to drive pivotal oncogenes in various malignancies. However, the landscape of SEs and SE-associated oncogenes remain elusive in NKTL. METHODS: We used Nano-ChIP-seq of the active enhancer marker histone H3 lysine 27 acetylation (H3K27ac) to profile unique SEs NKTL primary tumor samples. Integrative analysis of RNA-seq and survival data further pinned down high value, novel SE oncogenes. We utilized shRNA knockdown, CRISPR-dCas9, luciferase reporter assay, ChIP-PCR to investigate the regulation of transcription factor (TF) on SE oncogenes. Multi-color immunofluorescence (mIF) staining was performed on an independent cohort of clinical samples. Various function experiments were performed to evaluate the effects of TOX2 on the malignancy of NKTL in vitro and in vivo. RESULTS: SE landscape was substantially different in NKTL samples in comparison with normal tonsils. Several SEs at key transcriptional factor (TF) genes, including TOX2, TBX21(T-bet), EOMES, RUNX2, and ID2, were identified. We confirmed that TOX2 was aberrantly overexpressed in NKTL relative to normal NK cells and high expression of TOX2 was associated with worse survival. Modulation of TOX2 expression by shRNA, CRISPR-dCas9 interference of SE function impacted on cell proliferation, survival and colony formation ability of NKTL cells. Mechanistically, we found that RUNX3 regulates TOX2 transcription by binding to the active elements of its SE. Silencing TOX2 also impaired tumor formation of NKTL cells in vivo. Metastasis-associated phosphatase PRL-3 has been identified and validated as a key downstream effector of TOX2-mediated oncogenesis. CONCLUSIONS: Our integrative SE profiling strategy revealed the landscape of SEs, novel targets and insights into molecular pathogenesis of NKTL. The RUNX3-TOX2-SE-TOX2-PRL-3 regulatory pathway may represent a hallmark of NKTL biology. Targeting TOX2 could be a valuable therapeutic intervene for NKTL patients and warrants further study in clinic.

A composite single-nucleotide polymorphism prediction signature for extranodal natural killer/T-cell lymphoma.

Current prognostic scoring systems based on clinicopathologic variables are inadequate in predicting the survival and treatment response of extranodal natural killer/T-cell lymphoma (ENKTL) patients undergoing nonanthracyline-based treatment. We aimed to construct a classifier based on single-nucleotide polymorphisms (SNPs) for improving predictive accuracy and guiding clinical decision making. Data from 722 patients with ENKTL from international centers were analyzed. A 7-SNP-based classifier was constructed using LASSO Cox regression in the training cohort (n = 336) and further validated in the internal testing cohort (n = 144) and in 2 external validation cohorts (n = 142 and n = 100). The 7-SNP-based classifier showed good prognostic predictive efficacy in the training cohort and the 3 validation cohorts. Patients with high- and low-risk scores calculated by the classifier exhibited significantly different progression-free survival (PFS) and overall survival (OS) (all P < .001). The 7-SNP-based classifier was further proved to be an independent prognostic factor by multivariate analysis, and its predictive accuracy was significantly better than clinicopathological risk variables. Application of the 7-SNP-based classifier was not affected by sample types. Notably, chemotherapy combined with radiotherapy significantly improved PFS and OS vs radiotherapy alone in high-risk Ann Arbor stage I patients, whereas there was no statistical difference between the 2 therapeutic modalities among low-risk patients. A nomogram was constructed comprising the classifier and clinicopathological variables; it showed remarkably better predictive accuracy than either variable alone. The 7-SNP-based classifier is a complement to existing risk-stratification systems in ENKTL, which could have significant implications for clinical decision making for patients with ENKTL.

Analysis of N-nitrosodimethylamine in metformin hydrochloride products by high-resolution accurate mass gas chromatography mass spectrometry.

RATIONALE: The high resolving power of the Orbitrap mass spectrometer in a high-resolution accurate mass gas chromatography (HRAM-GC-MS) system provides greater selectivity and sensitivity for the identification and quantification of volatile analytes at low parts per billion (ppb) levels. Hence, it can be applied for the analysis of pharmaceutical impurities like N-nitrosodimethylamine (NDMA) in metformin hydrochloride products (METs). METHODS: Different METs extracted by a dichloromethane/aqueous system were analyzed by HRAM-GC-MS under softer electron ionization (EI) at 30 eV. The accurate masses of NDMA and its internal standard NDMA-d6 were analyzed by full scan and targeted selected ion monitoring modes under 60 000 and 30 000 full width at half maximum at m/z 200, respectively. Data acquisition and processing were managed by Xcalibur and Trace Finder software, respectively. RESULTS: Limits of detection (LOD) and quantification (LOQ) at 10 and 20 ng/g were achieved, which is below the allowed daily intake of 32 ng/g. The mass errors measured from experimental data were within ±2 ppm of the theoretical values over a period of a week. Sample analysis showed that 180 out of 212 samples (85%) were below LOD and 15 out of 212 samples (7 %) were within LOD and LOQ. Only 17 samples (8%) were found to be above LOQ, comprising one active pharmaceutical ingredient (API), five immediate-release METs and 11 extended-released METs. Amongst these, seven extended-release METs and one API exceeded the daily allowed intake, 32 ng/g. CONCLUSIONS: The validated method has been successfully applied for NDMA analysis in various forms of METs. The method is rather straightforward without an additional clean-up step. The scope can also be extended to other volatile impurities in finished pharmaceutical products.

Avelumab for the treatment of relapsed or refractory extranodal NK/T-cell lymphoma: an open-label phase 2 study.

This study aimed to assess the efficacy and safety of treatment with avelumab, an anti-programmed death ligand 1 (PD-L1) antibody, in patients with relapsed or refractory extranodal natural killer/T-cell lymphoma (ENKTL). In this phase 2 trial, 21 patients with relapsed or refractory ENKTL were treated with 10 mg/kg of avelumab on days 1 and 15 of a 28-day cycle. The primary end point was the complete response (CR) rate based on the best response. Targeted sequencing and immunohistochemistry were performed using pretreatment tumor tissue, and blood samples were drawn before and after treatment for measurement of cytokines and soluble programmed cell death protein 1 (PD1), PD-L1, and PD-L2. The CR rate was 24% (5 of 21), and the overall response rate was 38% (8 of 21). Although nonresponders showed early progression, 5 responders currently continue to receive treatment and have maintained their response. Most treatment-related adverse events were grade 1 or 2; no grade 4 adverse events were observed. Treatment responses did not correlate with mutation profiles, tumor mutation burden, serum levels of cytokines, or soluble PD1/PD-L1 and PD-L2. However, the response to avelumab was significantly associated with the expression of PD-L1 by tumor tissue (P = .001). Therefore, all patients achieving CR showed high PD-L1 expression, and their tumor subtyping based on PD-L1 expression correlated with treatment response. In summary, avelumab showed single-agent activity in a subset of patients with relapsed or refractory ENKTL. The assessment of PD-L1 expression on tumor cells might be helpful for identifying responders to avelumab. This trial was registered at www.clinicaltrials.gov as #NCT03439501.

A genomic-augmented multivariate prognostic model for the survival of natural-killer/T-cell lymphoma patients from an international cohort.

With lowering costs of sequencing and genetic profiling techniques, genetic drivers can now be detected readily in tumors but current prognostic models for Natural-killer/T cell lymphoma (NKTCL) have yet to fully leverage on them for prognosticating patients. Here, we used next-generation sequencing to sequence 260 NKTCL tumors, and trained a genomic prognostic model (GPM) with the genomic mutations and survival data from this retrospective cohort of patients using LASSO Cox regression. The GPM is defined by the mutational status of 13 prognostic genes and is weakly correlated with the risk-features in International Prognostic Index (IPI), Prognostic Index for Natural-Killer cell lymphoma (PINK), and PINK-Epstein-Barr virus (PINK-E). Cox-proportional hazard multivariate regression also showed that the new GPM is independent and significant for both progression-free survival (PFS, HR: 3.73, 95% CI 2.07-6.73; p < .001) and overall survival (OS, HR: 5.23, 95% CI 2.57-10.65; p = .001) with known risk-features of these indices. When we assign an additional risk-score to samples, which are mutant for the GPM, the Harrell's C-indices of GPM-augmented IPI, PINK, and PINK-E improved significantly (p < .001, χ2 test) for both PFS and OS. Thus, we report on how genomic mutational information could steer toward better prognostication of NKTCL patients.

In-depth characterization of the cisplatin mutational signature in human cell lines and in esophageal and liver tumors.

Cisplatin reacts with DNA and thereby likely generates a characteristic pattern of somatic mutations, called a mutational signature. Despite widespread use of cisplatin in cancer treatment and its role in contributing to secondary malignancies, its mutational signature has not been delineated. We hypothesize that cisplatin's mutational signature can serve as a biomarker to identify cisplatin mutagenesis in suspected secondary malignancies. Knowledge of which tissues are at risk of developing cisplatin-induced secondary malignancies could lead to guidelines for noninvasive monitoring for secondary malignancies after cisplatin chemotherapy. We performed whole genome sequencing of 10 independent clones of cisplatin-exposed MCF-10A and HepG2 cells and delineated the patterns of single and dinucleotide mutations in terms of flanking sequence, transcription strand bias, and other characteristics. We used the mSigAct signature presence test and nonnegative matrix factorization to search for cisplatin mutagenesis in hepatocellular carcinomas and esophageal adenocarcinomas. All clones showed highly consistent patterns of single and dinucleotide substitutions. The proportion of dinucleotide substitutions was high: 8.1% of single nucleotide substitutions were part of dinucleotide substitutions, presumably due to cisplatin's propensity to form intra- and interstrand crosslinks between purine bases in DNA. We identified likely cisplatin exposure in nine hepatocellular carcinomas and three esophageal adenocarcinomas. All hepatocellular carcinomas for which clinical data were available and all esophageal cancers indeed had histories of cisplatin treatment. We experimentally delineated the single and dinucleotide mutational signature of cisplatin. This signature enabled us to detect previous cisplatin exposure in human hepatocellular carcinomas and esophageal adenocarcinomas with high confidence.

Beyond fitness tracking: The use of consumer-grade wearable data from normal volunteers in cardiovascular and lipidomics research.

The use of consumer-grade wearables for purposes beyond fitness tracking has not been comprehensively explored. We generated and analyzed multidimensional data from 233 normal volunteers, integrating wearable data, lifestyle questionnaires, cardiac imaging, sphingolipid profiling, and multiple clinical-grade cardiovascular and metabolic disease markers. We show that subjects can be stratified into distinct clusters based on daily activity patterns and that these clusters are marked by distinct demographic and behavioral patterns. While resting heart rates (RHRs) performed better than step counts in being associated with cardiovascular and metabolic disease markers, step counts identified relationships between physical activity and cardiac remodeling, suggesting that wearable data may play a role in reducing overdiagnosis of cardiac hypertrophy or dilatation in active individuals. Wearable-derived activity levels can be used to identify known and novel activity-modulated sphingolipids that are in turn associated with insulin sensitivity. Our findings demonstrate the potential for wearables in biomedical research and personalized health.

Integrated paired-end enhancer profiling and whole-genome sequencing reveals recurrent CCNE1 and IGF2 enhancer hijacking in primary gastric adenocarcinoma.

OBJECTIVE: Genomic structural variations (SVs) causing rewiring of cis-regulatory elements remain largely unexplored in gastric cancer (GC). To identify SVs affecting enhancer elements in GC (enhancer-based SVs), we integrated epigenomic enhancer profiles revealed by paired-end H3K27ac ChIP-sequencing from primary GCs with tumour whole-genome sequencing (WGS) data (PeNChIP-seq/WGS). DESIGN: We applied PeNChIP-seq to 11 primary GCs and matched normal tissues combined with WGS profiles of >200 GCs. Epigenome profiles were analysed alongside matched RNA-seq data to identify tumour-associated enhancer-based SVs with altered cancer transcription. Functional validation of candidate enhancer-based SVs was performed using CRISPR/Cas9 genome editing, chromosome conformation capture assays (4C-seq, Capture-C) and Hi-C analysis of primary GCs. RESULTS: PeNChIP-seq/WGS revealed ~150 enhancer-based SVs in GC. The majority (63%) of SVs linked to target gene deregulation were associated with increased tumour expression. Enhancer-based SVs targeting CCNE1, a key driver of therapy resistance, occurred in 8% of patients frequently juxtaposing diverse distal enhancers to CCNE1 proximal regions. CCNE1-rearranged GCs were associated with high CCNE1 expression, disrupted CCNE1 topologically associating domain (TAD) boundaries, and novel TAD interactions in CCNE1-rearranged primary tumours. We also observed IGF2 enhancer-based SVs, previously noted in colorectal cancer, highlighting a common non-coding genetic driver alteration in gastric and colorectal malignancies. CONCLUSION: Integrated paired-end NanoChIP-seq and WGS of gastric tumours reveals tumour-associated regulatory SV in regions associated with both simple and complex genomic rearrangements. Genomic rearrangements may thus exploit enhancer-hijacking as a common mechanism to drive oncogene expression in GC.

Evaluation of the PIK3 pathway in peripheral T-cell lymphoma and NK/T-cell lymphoma.

Peripheral T-cell lymphomas (PTCL) and natural killer (NK)/T-cell lymphomas (NKTCL) are a heterogeneous group of aggressive malignancies with dismal outcomes and limited treatment options. While the phosphatidylinositol 3-kinase (PIK3) pathway has been shown to be highly activated in many B-cell lymphomas, its therapeutic relevance in PTCL and NKTCL remains unclear. The aim of this study is to investigate the expression of PIK3 and phosphatase and tensin homolog (PTEN) in these subtypes of lymphoma and to identify potential therapeutic targets for clinical testing. Therefore, the expression of PIK3α, PIK3ß, PIK3γ, PIK3δ and PTEN was analyzed in 88 cases of PTCL and NKTCL samples by immunohistochemistry. All PTCL and NKTCL samples demonstrated high expression of PIK3 isoforms. In particular, high PIK3α expression was significantly associated with poor survival, even after adjustment for age, International Prognostic Index (IPI) score and anthracycline-based chemotherapy in first line. Notably, copanlisib, a pan-class I inhibitor with predominant activities towards PIK3α and PIK3δ isoforms, effectively inhibited phosphorylation of AKT, 4E-BP-1 and STAT3, causing G0 /G1 cell cycle arrest and resulting in suppression of tumour cell growth in vitro and in vivo. This study provides evidence that targeting the PIK3 pathway, particularly simultaneous inhibition of PIK3α and δ, could be a promising approach for the treatment of PTCL and NKTCL.

Morphologic and genetic heterogeneity in breast fibroepithelial lesions-a comprehensive mapping study.

Breast fibroepithelial lesions (FELs) encompass the common fibroadenoma (FA) and relatively rare phyllodes tumour (PT); the latter entity is usually classified as benign, borderline or malignant. Intratumoural heterogeneity is frequently present in these tumours, making accurate histologic evaluation challenging. Despite their rarity, PTs are an important clinical problem due to their propensity for recurrence and, in the case of malignant PT, metastasis. Surgical excision is the mainstay of management. Recent work has uncovered myriad genetic alterations in breast FELs. In this study, exome sequencing was performed on seven cases of morphologically heterogeneous breast FELs, including FAs, PTs of all grades, and a case of metaplastic spindle cell carcinoma arising in PT, in order to elucidate their intratumoural genetic repertoire. Gene mutations identified encompassed cell signalling, tumour suppressor, DNA repair and cell cycle regulating pathways. Mutations common to multiple tumour regions generally showed higher variant allele frequency. Frequent mutations included MED12, TP53, RARA and PIK3CA. Histological observations of increased cellular density and pleomorphism correlated with mutational burden. Phylogenetic analyses revealed disparate pathways of possible tumour progression. In summary, histological heterogeneity correlated with genetic changes in breast FELs.

Oncogenic activation of the STAT3 pathway drives PD-L1 expression in natural killer/T-cell lymphoma.

Mature T-cell lymphomas, including peripheral T-cell lymphoma (PTCL) and extranodal NK/T-cell lymphoma (NKTL), represent a heterogeneous group of non-Hodgkin lymphomas with dismal outcomes and limited treatment options. To determine the extent of involvement of the JAK/STAT pathway in this malignancy, we performed targeted capture sequencing of 188 genes in this pathway in 171 PTCL and NKTL cases. A total of 272 nonsynonymous somatic mutations in 101 genes were identified in 73% of the samples, including 258 single-nucleotide variants and 14 insertions or deletions. Recurrent mutations were most frequently located in STAT3 and TP53 (15%), followed by JAK3 and JAK1 (6%) and SOCS1 (4%). A high prevalence of STAT3 mutation (21%) was observed specifically in NKTL. Novel STAT3 mutations (p.D427H, E616G, p.E616K, and p.E696K) were shown to increase STAT3 phosphorylation and transcriptional activity of STAT3 in the absence of cytokine, in which p.E616K induced programmed cell death-ligand 1 (PD-L1) expression by robust binding of activated STAT3 to the PD-L1 gene promoter. Consistent with these findings, PD-L1 was overexpressed in NKTL cell lines harboring hotspot STAT3 mutations, and similar findings were observed by the overexpression of p.E616K and p.E616G in the STAT3 wild-type NKTL cell line. Conversely, STAT3 silencing and inhibition decreased PD-L1 expression in STAT3 mutant NKTL cell lines. In NKTL tumors, STAT3 activation correlated significantly with PD-L1 expression. We demonstrated that STAT3 activation confers high PD-L1 expression, which may promote tumor immune evasion. The combination of PD-1/PD-L1 antibodies and STAT3 inhibitors might be a promising therapeutic approach for NKTL, and possibly PTCL.

Genomic characterisation of breast fibroepithelial lesions in an international cohort.

Fibroepithelial lesions (FELs) are a heterogeneous group of tumours comprising fibroadenomas (FAs) and phyllodes tumours (PTs). Here we used a 16-gene panel that was previously discovered to be implicated in pathogenesis and progression, to characterise a large international cohort of FELs via targeted sequencing. The study comprised 303 (38%) FAs and 493 (62%) PTs which were contributed by the International Fibroepithelial Consortium. There were 659 (83%) Asian and 109 (14%) non-Asian FELs, while the ethnicity of the rest was unknown. Genetic aberrations were significantly associated with increasing grade of PTs, and were detected more in PTs than FAs for MED12, TERT promoter, RARA, FLNA, SETD2, TP53, RB1, EGFR, and IGF1R. Most borderline and malignant PTs possessed ≥ 2 mutations, while there were more cases of FAs with ≤ 1 mutation compared to PTs. FELs with MED12 mutations had significantly higher rates of TERT promoter, RARA, SETD2, EGFR, ERBB4, MAP3K1, and IGF1R aberrations. However, FELs with wild-type MED12 were more likely to express TP53 and PIK3CA mutations. There were no significant differences observed between the mutational profiles of recurrent FAs, FAs with a history of subsequent ipsilateral recurrence or contralateral occurrence, and FAs without a history of subsequent events. We identified recurrent mutations which were more frequent in PTs than FAs, with borderline and malignant PTs harbouring cancer driver gene and multiple mutations. This study affirms the role of a set of genes in FELs, including its potential utility in classification based on mutational profiles. © 2019 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

An integrated package for bisulfite DNA methylation data analysis with Indel-sensitive mapping.

BACKGROUND: DNA methylation plays crucial roles in most eukaryotic organisms. Bisulfite sequencing (BS-Seq) is a sequencing approach that provides quantitative cytosine methylation levels in genome-wide scope and single-base resolution. However, genomic variations such as insertions and deletions (indels) affect methylation calling, and the alignment of reads near/across indels becomes inaccurate in the presence of polymorphisms. Hence, the simultaneous detection of DNA methylation and indels is important for exploring the mechanisms of functional regulation in organisms. RESULTS: These problems motivated us to develop the algorithm BatMeth2, which can align BS reads with high accuracy while allowing for variable-length indels with respect to the reference genome. The results from simulated and real bisulfite DNA methylation data demonstrated that our proposed method increases alignment accuracy. Additionally, BatMeth2 can calculate the methylation levels of individual loci, genomic regions or functional regions such as genes/transposable elements. Additional programs were also developed to provide methylation data annotation, visualization, and differentially methylated cytosine/region (DMC/DMR) detection. The whole package provides new tools and will benefit bisulfite data analysis. CONCLUSION: BatMeth2 improves DNA methylation calling, particularly for regions close to indels. It is an autorun package and easy to use. In addition, a DNA methylation visualization program and a differential analysis program are provided in BatMeth2. We believe that BatMeth2 will facilitate the study of the mechanisms of DNA methylation in development and disease. BatMeth2 is an open source software program and is available on GitHub ( https://github.com/GuoliangLi-HZAU/BatMeth2 /).

BatAlign: an incremental method for accurate alignment of sequencing reads.

Structural variations (SVs) play a crucial role in genetic diversity. However, the alignments of reads near/across SVs are made inaccurate by the presence of polymorphisms. BatAlign is an algorithm that integrated two strategies called 'Reverse-Alignment' and 'Deep-Scan' to improve the accuracy of read-alignment. In our experiments, BatAlign was able to obtain the highest F-measures in read-alignments on mismatch-aberrant, indel-aberrant, concordantly/discordantly paired and SV-spanning data sets. On real data, the alignments of BatAlign were able to recover 4.3% more PCR-validated SVs with 73.3% less callings. These suggest BatAlign to be effective in detecting SVs and other polymorphic-variants accurately using high-throughput data. BatAlign is publicly available at https://goo.gl/a6phxB.

Determination of N-nitroso folic acid in folic acid and multivitamin supplements by LC-MS/MS.

Analysis of N-nitroso folic acid, a nitrosamine impurity found in folic acid, is challenging due to the complex sample matrices. Many of such supplements contain not only a variety of vitamins, including vitamin A, Bs, C, D and E, but also other ingredients such as minerals, docosahexaenoic acid (DHA), glucose syrup, sugar, and herbs. On the other hand, the strength of folic acid is typically low, ranging from 50 µg to 5 mg per unit. In this study, a highly selective and sensitive LC-MS/MS method was developed to accurately quantify N-nitroso folic acid in supplements containing folic acid. The sample was extracted by 0.1% ammonia solution: MeOH (9:1, v/v) containing 5 ng/mL of N-nitroso folic acid-d4 (Isotope internal standard). The quantification was performed by MRM in negative ionization mode. Mobile phases A and B were 0.1% formic acid in deionized water and methanol, respectively. The method was validated and found to have sufficient linearity (R2 > 0.995), accuracy (recovery 83-110%), precision (RSD 3%) and low LOD, LOQ (4 and 10 µg/g respectively, with respect to folic acid). The method was applied to the determination of N-nitroso folic acid in 40 supplements containing folic acid with different strengths and formulation. The content of N-nitroso folic acid was found to be up to 898 ng/unit (1794 µg/g with respect to folic acid). It enabled regulatory actions, such as product recall, to safeguard public health from unsafe products.

ASTER: A Method to Predict Clinically Relevant Synthetic Lethal Genetic Interactions.

A Synthetic Lethal (SL) interaction is a functional relationship between two genes or functional entities where the loss of either entity is viable but the loss of both is lethal. Such pairs can be used to develop targeted anticancer therapies with fewer side effects and reduced overtreatment. However, finding clinically relevant SL interactions remains challenging. Leveraging unified gene expression data of both disease-free and cancerous samples, we design a new technique based on statistical hypothesis testing, called ASTER, to identify SL pairs. We empirically find that the patterns of mutually exclusivity ASTER finds using genomic and transcriptomic data provides a strong signal of synthetic lethality. For large-scale multiple hypothesis testing, we develop an extension called ASTER++ that can utilize additional input gene features within the hypothesis testing framework. Our computational and functional experiments demonstrate the efficacy of ASTER in identifying SL pairs with potential therapeutic benefits.