List of fold change proteomic data sets of dFCM- 39 vs. 12Day and105 vs. 12Day, cFCM- 40 vs. 12Day and115 vs. 12Day , mouse heart 90 vs. 1 day after selecting mitochondrial proteins from the scaffold software for Ingenuity Pathway Analysis (IPA Qiagen)

Renal cell carcinoma (RCC), particularly the clear cell subtype (ccRCC), poses a significant global health concern due to its increasing prevalence and resistance to conventional therapies. Early detection of ccRCC remains challenging, resulting in poor patient survival rates. In this study, we employed a bioinformatic approach to identify potential prognostic biomarkers for kidney renal clear cell carcinoma (KIRC). By analyzing RNA sequencing data from the TCGA-KIRC project, differentially expressed genes (DEGs) associated with ccRCC were identified. Pathway analysis utilizing the Qiagen Ingenuity Pathway Analysis (IPA) tool elucidated key pathways and genes involved in ccRCC dysregulation. Prognostic value assessment was conducted through survival analysis, including Cox univariate proportional hazards (PH) modeling and Kaplan–Meier plotting. This analysis unveiled several promising biomarkers, such as MMP9, PIK3R6, IFNG, and PGF, exhibiting significant associations with overall survival and relapse-free survival in ccRCC patients. Cox multivariate PH analysis, considering gene expression and age at diagnosis, further confirmed the prognostic potential of MMP9, IFNG, and PGF genes. These findings enhance our understanding of ccRCC and provide valuable insights into potential prognostic biomarkers that can aid healthcare professionals in risk stratification and treatment decision-making. The study also establishes a foundation for future research, validation, and clinical translation of the identified prognostic biomarkers, paving the way for personalized approaches in the management of KIRC.

TSPO-associated protein 1 (TSPOAP1) is a cytoplasmic protein and is closely associated with its mitochondrial transmembrane protein partner translocator protein (TSPO). To decipher the canonical signalling pathways of TSPOAP1, its role in human diseases and disorders, and relationship with TSPO; expression analyses of TSPOAP1- and TSPO-associated human genes were performed by Qiagen Ingenuity Pathway Analysis (IPA). In the expression analysis, necroptosis and sirtuin signalling pathways, mitochondrial dysfunction, and inflammasome were the top canonical pathways for both TSPOAP1 and TSPO, confirming the close relationship between these two proteins. A distribution analysis of common proteins in all the canonical pathways predicted for TSPOAP1 revealed that tumor necrosis factor receptor 1 (TNFR1), vascular cell adhesion molecule 1 (VCAM1), cyclic AMP response element-binding protein 1 (CREB1), T-cell receptor (TCR), nucleotide-binding oligomerization domain, leucine-rich repeat and pyrin domain containing 3 (NLRP3), DNA-dependent protein kinase (DNA-PK or PRKDC), and mitochondrial permeability transition pore (mPTP) were the major interaction partners of TSPOAP1, highlighting the role of TSPOAP1 in inflammation, particularly neuroinflammation. An analysis of the overlap between TSPO and TSPOAP1 Homo sapiens genes and top-ranked canonical pathways indicated that TSPO and TSPOAP1 interact via voltage-dependent anion-selective channels (VDAC1/2/3). A heat map analysis indicated that TSPOAP1 has critical roles in inflammatory, neuroinflammatory, psychiatric, and metabolic diseases and disorders, and cancer. Taken together, this information improves our understanding of the mechanism of action and biological functions of TSPOAP1 as well as its relationship with TSPO; furthermore, these results could provide new directions for in-depth functional studies of TSPOAP1 aimed at unmasking its detailed functions.

Six female littermate piglets were used in an experiment to evaluate the mRNA expression in tissues from piglets given one or two 1 mL injections of iron dextran (200 mg Fe/mL). All piglets in the litter were administered the first 1 mL injection < 24 hours after birth. On d 7, piglets were paired by weight (mean BW = 1.72 ± 0.13 kg) and one piglet from each pair was randomly selected as control (CON) and the other received a second injection (+Fe). At weaning on d 22, each piglet was anesthetized, and samples of liver and duodenum were taken from the anesthetized piglets and preserved until mRNA extraction. Differential Gene Expression data were analyzed with a fold-change cutoff (FC) of |1.2| P < 0.05. Pathway analysis was conducted with Z-score cutoff of P < 0.05. In the duodenum 435 genes were significantly changed with a FC ≥ |1.2| P < 0.05. In the duodenum, Claudin 1 and Claudin 2 were inversely affected by +Fe. Claudin 1 (CLDN1) plays a key role in cell-to-cell adhesion in the epithelial cell sheets and was upregulated (FC = 4.48, P = 0.0423). Claudin 2 (CLDN2) is expressed in cation leaky epithelia, especially during disease or inflammation and was downregulated (FC = -1.41, P = 0.0097). In the liver, 362 genes were expressed with a FC ≥ |1.2| P < 0.05. The gene most affected by a second dose of 200 mg Fe was HAMP (Hepcidin Antimicrobial Peptide) with a FC of 40.8. HAMP is a liver-produced hormone that is the main circulating regulator of Fe absorption and distribution across tissues. It also controls the major flows of Fe into plasma by promoting endocytosis and degradation of ferroportin (SLC4A1). This leads to the retention of Fe in Fe-exporting cells and decreased flow of Fe into plasma. Metabolic pathway changes in the duodenum and liver provide evidence for the improved feed conversion and growth rates in piglets given two iron injections pre-weaning with contemporary pigs in a companion study. In the duodenum, there is a down regulation of gene clusters associated with gluconeogenesis (P < 0.05). Concurrently, there was a decrease in the mRNA expression of genes for enzymes required for urea production in the liver (P < 0.05). These observations suggest that there may be less need for gluconeogenesis, and possibly less urea production from deaminated amino acids. The genomic and pathway analyses provided empirical evidence linking gene expression with phenotypic observations of piglet health and growth improvements. Methods RNA Sequencing Samples were submitted to Zymo Research (Irvine, CA, USA) for total mRNA extraction, cDNA library preparation and RNA sequencing. Total RNA-Seq libraries were constructed from 250 ng of total RNA. To remove rRNA, a method described by Bogdanova et al., 2011 was followed. Libraries were prepared using the Zymo-Seq RiboFree Total RNA Library Prep Kit (Cat # R3000) according to the manufacturer’s instructions (Zymo-Research, 2022). The RNA- Seq libraries were sequenced on an Illumina NovaSeq to a sequencing depth of at least 30 million read pairs (150 bp paired-end sequencing) per sample. Detection of DGE and Bioinformatic Data Handling Differentially expressed genes were detected using GeneSpring software (Agilent, Santa Clara, CA) using selection criteria that accepted a DGE threshold of greater than a |1.2|-FC in expression level and statistical probability levels of P < 0.05. The filtered genes were then subjected to Ingenuity Pathway Analysis (IPA; QIAGEN Inc., Redwood City, CA; https://digitalinsights.qiagen.com/products/) to gain insights into canonical pathways, networks, and biological functions. Qiagen IPA uses algorithms, tools, and visualizations to combine the directional information from gene expression patterns (up- or down-regulation) with the expected causal effects of the genes, as reported in the published literature. A prediction for effects of a treatment on a particular biological pathway function or disease can then be made based upon the direction of change in gene expression and calculated Z-scores. Briefly, a Z-score is used to compare data that have different means and standard deviations. The Z-score is the distance of a point, such as a complete pathway, from the mean of the distribution in terms of the standard deviation (Corchete et al., 2020; Shao et al., 2020; Wieder, Lai and Ebbels, 2022).

Ordered cellular architecture and high concentrations of stable crystallins are required for the lens to maintain transparency. Here we investigate the molecular mechanism of cataractogenesis of the CRYGC c.119-123dupGCGGC (p.Cys42AlafsX63) (CRYGC5bpd) mutation. Lenses were extracted from wild-type and transgenic mice carrying the CRYGC5bpdup minigene and RNA was isolated and converted into cDNA. Expression of genes in the unfolded protein response (UPR) pathways was estimated by qRT-PCR and RNA seq and pathway analysis was carried out using the Qiagen IPA website. P3W Transgenic mice exhibited phenotypic diversity with a dimorphic population of severe and clear lenses. PCA of RNA seq data showed separate clustering of wild-type, clear CRYGC5bpd, and severe CRYGC5bpd lenses. Transgenic mice showed differential upregulation in Master regulator Grp78 (Hspa5) and downstream targets in the PERK-dependent UPR pathway including Atf4 and Chop (Ddit3), but not GADD34. Thus, high levels of CRYGC..., RNA isolation, cDNA synthesis, and qRT-PCR Total RNA was isolated using an RNA isolation kit (The RNeasy Plus Mini Kit; Qiagen, Valencia, CA) and quantified using a spectrophotometer (Nanodrop 2000C; ThermoFisher). A first-strand cDNA was synthesized from approximately 0.5mg of total RNA by cDNA synthesis kit (Super III first-strand synthesis for RT PCR kit; Invitrogen) according to the manufacturer's protocol. qRT-PCR was performed using Applied Biosystems ViiA7 Real-Time PCR system with the following amplification conditions: an initial incubation of the samples at 50°C for 2min and denaturation at 95°C 15min followed by 40 cycles of denaturation, annealing, and extension at 95°C 15sec, 60°C 30sec, and 72°C 30sec. Gapdh was used as an endogenous control for normalizing the target mRNA. The relative expression of each target gene was calculated using the 2^(∆∆Ct) method. The primers were standardized, and efficiencies were tested before performing qRT-PCR. RNA-Seq About 200ng of ..., , # The c.119-123dup5bp mutation in human gamma-C-crystallin destabilizes the protein and activates the unfolded protein response to cause highly variable cataracts

https://doi.org/10.5061/dryad.rn8pk0pmm

Description of the data and file structure

The RNASeq and qRT-PCR files included refer to mice transgenic for a CRYGC c.119-123dupGCGGC (p.Cys42AlafsX63) (CRYGC5bpd) mutation. Lenses were extracted from wild-type and transgenic mice carrying the CRYGC5bpdup minigene and RNA was isolated and converted into cDNA and submitted to Novogene for RNASeq analysis. The descriptions of the mice are given in Ma, Z. et al. Overexpression of human γC-crystallin 5bp duplication Disrupts Lens Morphology in Transgenic Mice. Invest Ophthalmol Vis Sci 52, 5269-5375 (2011).

Identification of CRYGC as the causative gene is described in Ren, Z. et al. A 5-base insertion in the γC-crystallin gene is associated with autosomal dominant variable zonular pulveru...

Up-regulated and down-regulated molecules are shown with up and down arrows, respectively.

A significant area of rare diseases research is the investigation of druggable protein encoding genes that contribute to pathogenesis. Niemann Pick Type C (NPC) disease can be considered as a challenging one, among all rare diseases, because it has been associated with a poor prognosis and an unclear molecular pathogenesis. In recent years, Proteomics analysis has become a functional and useful technology for profiling protein expression and find possible drugs targets. In the present study, hepatocytes derived from wild type and Npc1 deficient mice were analyzed by mass spectrometry-based proteomics, followed by pathway analysis and statistical interpretation, performed with the QIAGEN Ingenuity Pathway Analysis (IPA) software. Applications for protein function was built using IPA and gene ontology analysis. We identified and reliably quantified a total of 3833 proteins, among them 416 presented a significant p-value (<0.05) being classified 149 as upregulated (log2 fold change >1) and 6 as downregulated (Log2 fold change <-1). Our analysis revealed that the most significant changed proteins are related to liver damage, lipid metabolism and inflammatory response. In addition, in the group of up/downregulated proteins, 47 proteins were identified as lysosomal proteins and 22 as mitochondrial proteins. Importantly, we found that of those proteins CTSB, LIPA, DPP7, GLMP and DECR1 are related to liver fibrosis, liver damage and steatosis.

The project contains raw and result files from a proteomic profiling of a male breast cancer (MBC) case. Label-free quantification-mass spectrometry (LFQ-MS) and bioinformatics analysis were employed to investigate the differentially expressed proteins (DEPs) among distinct tissue samples: the primary breast tumor, axillary metastatic lymph nodes and the adjacent non-tumor breast tissue. An additional proteomic comparative analysis was performed with a primary breast tumor of a female patient. A number of Ingenuity® Pathway Analysis (IPA) (QIAGEN Inc.) and functional annotation tools were used to further analyze the DEPs. Altogether, our findings revealed deregulated proteins into signaling pathways involved in the cancer development and provided a landscape of proteomic data for the MBC research.

Macrophages play a key role in ozone-induced lung injury by regulating both the initiation and resolution of inflammation. These distinct activities are mediated by pro-inflammatory and anti-inflammatory/pro-resolution macrophages which sequentially accumulate in injured tissues. Macrophage activation is dependent, in part, on intracellular metabolism. Herein, we used RNA-sequencing (seq) to identify signaling pathways regulating macrophage immunometabolic activity following exposure of mice to ozone (0.8 ppm, 3 hr) or air control. Analysis of lung macrophages using an Agilent Seahorse showed that inhalation of ozone increased macrophage glycolytic activity and oxidative phosphorylation at 24 and 72 hr post exposure. An increase in the percentage of macrophages in the S phase of the cell cycle was observed 24 hr post ozone. RNA-seq revealed significant enrichment of pathways involved in innate immune signaling and cytokine production among differentially expressed genes at both 24 and 72 hr after ozone, while pathways involved in cell cycle regulation were upregulated at 24 hr and intracellular metabolism at 72 hr. An interaction network analysis identified tumor suppressor 53 (TP53), E2F family of transcription factors (E2Fs), Cyclin Dependent Kinase Inhibitor 1A (CDKN1a/p21), and Cyclin D1 (CCND1) as upstream regulators of cell cycle pathways at 24 hr and TP53, nuclear receptor subfamily 4 group a member 1 (NR4A1/Nur77), and estrogen receptor alpha (ESR1/ERα) as central upstream regulators of mitochondrial respiration pathways at 72 hr. These results highlight the complex interaction between cell cycle, intracellular metabolism, and macrophage activation which may be important in the initiation and resolution of inflammation following ozone exposure. Methods Total RNA was extracted as described above from 3 mice/treatment group. In a pilot study, we found that 3 mice were sufficient to identify a significant difference in Ptgs2 gene expression by qPCR at α = 0.05 and power = 80%. RNA integrity numbers (RINs) were confirmed to be ≥ 8.8 using a 2100 Bioanalyzer Instrument (Agilent, Santa Clara, CA). cDNA libraries were prepared using mouse TruSeq® Stranded Total RNA Library Prep kit (illumina, San Diego, CA) and quantified using a KAPA Library Quantification kit (Roche, Pleasanton, CA). cDNA libraries were sequenced (75 bp single-ended, ~35-44M reads per sample) on an Illumina NextSeq instrument. Raw reads in FastQ files were trimmed using Trimmomatic-0.39 (Bolger et al. 2014) and quality control of trimmed files performed using FastQC. Salmon was used to align reads in mapping-based mode with selective alignment against a decoy-aware transcriptome generated from mouse transcriptome GENCODE Release M23 (GRCm38.p6). Estimated counts per transcript were generated using the gcBias flag and normalized to transcript length to correct for potential changes in gene length across samples from differential isoform usage (Love et al. 2016; Patro et al. 2017). Transcript level quantitation data were aggregated to the gene-level using tximport (Soneson et al. 2015). Differential gene expression analysis was performed with air exposed mice as controls using DESeq2 with corrections for differences in library size (Love et al. 2014) in R version 4.0.3. Significantly enriched canonical pathways and upstream regulators were identified with Ingenuity IPA Version 65367011 (QIAGEN Inc, https://www.qiagenbioinformatics.com/products/ingenuity-pathway-analysis/) using a right-tailed Fisher’s Exact Test (Krämer et al. 2014). A less stringent criteria (fold change > 1.3 and experimental false discovery rate [padj] < 0.05) was used to augment the number of genes included in the pathway analysis (Bennett et al. 2024). Data were deposited NCBI’s Gene Expression Omnibus (Edgar et al. 2002) and are accessible through GEO Series accession number GSE237594 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE237594).

Summary:The datasets described here were gathered while investigating the molecular processes by which some human ductal carcinoma in situ (DCIS) lesions advance to the more aggressive form while others remain indolent.Data access:All RNA sequencing data have been deposited in the Gene Expression Omnibus with accession https://identifiers.org/geo:GSE143790 All Chip-Exo data have been deposited in the Gene Expression Omnibus with accession https://identifiers.org/geo:GSE143313 RPPA data is included together with this data record, in the file Supplementary Figure 3-RPPA.xlsx. The Raw RPPA and ANOVA tabs are the results, while the other tabs are IPA analysis performed by the authors. Permission to use figures and data generated using QIAGEN Ingenuity Pathway Analysis (IPA) is given in the file QIAGEN Ingenuity Product Support Permission letter for Dr. Behbod.pdf.The specific data underlying each figure and supplementary figure in the manuscript are provided as part of this data record, and are as follows:Figure 1-BCL9-STAT3 interaction.xlsxFigure 2-ChIP Exo.xlsxFigure 3-ChIP.xlsxFigure 4-MMP16 and avb3 MIND xenografts.xlsxFigure 5-MMP16 avb3 TMA analysis.xlsxFigure 6-Carnosic data.xlsxSupplementary Figure 3-RPPA.xlsxSupplementary Figure 4-STAT3 Reporter.xlsxSupplementary Figure 5-ChIP Exo Motifs.xlsxSupplementary Figure 6-integrin data.xlsxSupplementary Figure 7-MMP data.xlsxStudy approval and patient consent: Patients gave written informed consent for participation in the University of Kansas Medical Center Institutional Review Board–approved study allowing collection of additional biopsies and or surgical tissue for research. Animal experiments were conducted following protocols approved by the University of Kansas School of Medicine Animal Care and Use and Human Subjects Committee. Study aims and methodology: The aim of the related study was to determine the molecular processes underlying progression to invasion in DCIS using PDX DCIS MIND animal models. Using a novel intraductal model they identify downregulation of specific STAT3 targets to promote progression and use a purified component from rosemary extract to show that treatment in vivo decreases DCIS progression in patient derived DCIS and cell line models.

Pathways are arranged according to the rank for obese versus lean ZDF rats. The probability for a difference between ZDF obese versus lean and Wistar-STZ versus Wistar respectively, p(FDR), is corrected for multiple testing by a Benjamini-Hochberg false discovery rate (FDR) methodology and shown here is negative common logarithm, -lg p(FDR).

LC-MS/MS raw dataSpectrum matching and protein identification and validation were performed with MSFragger, and quantification of protein intensities with matching between runs was performed with IonQuant as components of the FragPipe analysis pipeline using the default settings of each module. The protein database used for the search was the Mus musculus reviewed sequence database downloaded from UniProt on June 1, 2023. The results were subsequently processed to filter out common contaminants, decoy hits from the reverse database, and protein groups identified by a single peptide. The data were filtered as follows: (a) binary expression of a protein (i.e., protein exclusively identified in either scLRP1+/+ or scLRP1-/-) was only considered relevant if all scLRP1+/+ samples or all scLRP1-/- samples expressed the protein. The missing values were imputed with the minimum intensity value for each specific data set; (b) for samples expressed in both scLRP1+/+ and scLRP1-/- tissue, the filtering process required 2 or more proteins to be detected in both scLRP1+/+ and scLRP1-/- samples. False discovery analysis was performed using the Benjamini, Krieger, and Yekutieli method using GraphPad Prism 10.0 software. Causal analysis of proteomic data was performed in IPA upstream analysis software (QIAGEN). For IPA, the binary values were imputed using local minimum intensities. Enrichment analyses for gene ontology (biological process) were performed using clusterProfiler 4.2.2 R package on R 4.1.0. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the data set identifier PXD038236. Global quantification of protein expression. Sciatic nerves from scLRP1+/+ and scLRP1-/- mice were rinsed in PBS to remove the blood and frozen with liquid nitrogen in cryogenic storage tubes (#5016-0001, Thermo Fisher Scientific). Fractured tissue was transferred to a 1 mL milliTUBE containing an AFA fiber (#520135, Covaris catalog) in 200 μL of 50 mM HEPES pH 8.5, 150 mM NaCl, and 2% Triton X-114 and sonicated with a M220 Focused-Ultrasonicator (#500295, Covaris catalog). Sonication parameters were temperature 15 °C, peak power 75 W, duty factor 26, cycles/burst = 1,000, and duration 600 seconds. Extracted proteins were clarified of insoluble material by centrifugation at 15,000 x g for 20 minutes at 4 °C. Protein concentrations were determined with the Micro BCA colorimetric assay (Pierce, Thermo Fisher Scientific) with the addition of SDS to a final concentration of 1% in the assay solvent to prevent detergent clouding. Aliquots containing approximately 5 μg of protein were processed using the SP3 protocol as described (89) with some modifications. Briefly, the sample aliquots were brought to 50 μL volume, and disulfide bonds were reduced and alkylated simultaneously with 10 mM TCEP, 40 mM 2-chloroacetamide in 50 mM HEPES pH 8.5, and 1% sodium deoxycholate at 70 °C for 10 minutes, then cooled on ice. Proteins were precipitated and captured following the addition of 10 μL of a washed 10 μg/μL suspension of SpeedBeads (Cytiva) and 400 μL of ethanol. After shaking for 10 minutes at room temperature, the beads were magnetically captured and washed 3 times with 200 μL of 80% ethanol in water. Proteins were digested on the beads in 50 μL of 50 mM HEPES pH 8.5, 1% sodium deoxycholate, and 10 ng/μL trypsin (Promega) overnight at room temperature with shaking sufficient to maintain the beads in suspension. The digest was diluted 10-fold with 80% acetonitrile and 1% formic acid, then separated from the beads magnetically, and the resulting peptides were captured on 2 mm discs of Empore Cation (CDS Analytical) fitted into 1,000 μL pipette tips (Sartorious catalog 791000). Detergents and other contaminants were removed by washing the tips serially with 1) ethyl acetate; 2) 80% acetonitrile, 1% formic acid; and 3) 10% acetonitrile, 0.2% formic acid. Peptides were eluted directly into injection vials with freshly prepared 80% acetonitrile and 5% ammonium hydroxide and immediately dried down in a centrifugal vacuum evaporator. One-fifth of the recovered peptides from each sample were subsequently analyzed by liquid chromatography-tandem mass spectrometry. In-house capillary columns were constructed from 360 μm OD and 100 μm internal diameter × 30 cm fused silica tubing (Molex) with laser-pulled tips (Sutter Instruments) and were packed with Reprosil-PUR 3 μm C18-AQ (Dr. Maisch GmBH). Solvents A and B consisted of 0.1% formic acid in water and 80% acetonitrile with 0.1% formic acid, respectively. A 180-minute linear gradient from 2% to 35% solvent B was used for chromatographic separation. Peptides were analyzed with an Orbitrap Elite (Thermo Fisher Scientific) mass spectrometer using nano-electrospray ionization with an applied voltage of 1,800 V. MS1 spectra were acquired at a resolution of 120,000, and the 15 most abundant precursor ions were selected for fragmentation by higher energy collision dissociation. MS2 spectra were acquired at a resolution of 15,000. Dynamic exclusion parameters were a list size of 500, a mass window of ±7 ppm, and a duration of 1 minute. Automatic gain control settings were MS1 target 1 × 106, maximum inject time 100 ms; MS2 target 4 × 104, maximum inject time 100 ms. Principal components of the 8 samples (2 groups: 4 scLRP1+/+, 4 scLRP1-/-) were analyzed. The centroid of each group, generated by the K-nearest neighbor (KNN) algorithm, was used to define each cluster. All samples from each group were restricted to the same cluster with no overlap.RNA-Seq dataL4 and L5 DRGs from the left and right sides of each mouse were acutely isolated, pooled, and snap-frozen in liquid nitrogen (n=3 per genotype). RNA was extracted with the AllPrep DNA/RNA Micro Kit (Qiagen, Inc.). RNA quality was assessed on an Agilent 2100 Bioanalyzer. Samples with RNA integrity numbers ≥8 were used for RNA sequencing (RNA-seq). RNA-seq was conducted at the NYU Genomic Core. All samples were processed in the same time period, following the same protocol to limit batch effects and other confounders using the Illumina RNA-seq platform. The cDNA library was prepared via the standard Illumina protocol. Three samples were pooled into each lane and sequenced by 75-bp paired-end sequencing on an Illumina HiSeq 2500 using standard protocols. A Phi-X positive control provided by Illumina was spiked into all lanes at a concentration of 0.3% to monitor sequencing quality. The sample error rate was < 2% and the distribution of reads per sample in a lane was within a reasonable tolerance. Data generated during sequencing runs were transferred to the high-performance computing (HPC) cluster, with individual base calls transferred for downstream analysis. The target for average reads per sample was approximately 25 million. The QC pipeline included: 1) quality check of the raw sequencing data using FastQC (v 0.11.9) and MultiQC (v 1.9); 2) mapping the sequencing reads to the human genome (build 102) using HISAT2 (v 2.2.1), followed by SAMtools (v 1.12) to convert BAM (Binary Alignment Map) into SAM (Sequence Alignment Map) files; 3) assembly of RNA-seq reads into transcripts using StringTie (v 2.1.4); and 4) calculation of expression levels from read counts, producing a gene count matrix. To improve the rigor of gene expression analysis, a minimum read depth of 10 gene counts in at least 90% of samples was required for retention in the matrix. The Ensembl identifiers (ID) of the gene counts were annotated to Entrez IDs using the EnrichmentBrowser (v.2.18.2) package in R. The Entrez IDs were annotated to gene symbols using Homo sapiens (v. 1.3.1). Gene expression was compared using DESeq2 (v. 1.32.0), which assumes a negative binomial distribution using the mean and variance estimated from gene counts and p-values calculated using the Wald test. GSEA was performed using GSEA 3.0 (http://www.broadinstitute.org/gsea/). Adjusted p < 0.05 and a false discovery rate (FDR) q < 0.25 were considered significant. DEGs were selected when presented with an adjusted p-value < 0.05 and a log2FC (fold change) > 1 for upregulated genes and < -1 for downregulated genes. The Database for Annotation, Visualization, and Integrated Discovery (DAVID, http://david.abcc.ncifcrf.gov/, version 6.8) was used to perform GO and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses. Set enrichment analysis was used for the pathway by selecting non-significant differentially expressed genes specified as the “background universe” and accounting for multiple testing using a false discovery rate of q < 0.1. Differentially altered pathways were evaluated by using the enrich plot package in R for visualization of functional enrichment (i.e., dot plot).

The project contains raw and result files from a comparative proteomic analysis of malignant [primary breast tumor (PT) and axillary metastatic lymph nodes (LN)] and non-tumor [contralateral (NCT) and adjacent breast (ANT)] tissues of patients diagnosed with invasive ductal carcinoma. A label-free mass spectrometry was conducted using nano-liquid chromatography coupled to electrospray ionization–mass spectrometry (LC-ESI-MS/MS) followed by functional annotation to reveal differentially expressed proteins and their predicted impacts on pathways and cellular functions in breast cancer. A total of 462 proteins was observed as differentially expressed (DEPs) among the groups of samples analyzed. Ingenuity Pathway Analysis software version 2.3 (QIAGEN Inc.) was employed to identify the most relevant signaling and metabolic pathways, diseases, biological functions and interaction networks affected by the deregulated proteins. Upstream regulator and biomarker analyses were also performed by IPA’s tools. Altogether, our findings revealed differential proteomic profiles that affected the associated and interconnected cancer signaling processes.

Mutually enriched IPA Canonical pathways that belong to sphingolipid signaling in GBMvsCTRL and MSvsCTRL contrasts.

Additional file 2. Table S1. List of somatic mutations identified in each patient. Table S2. List, information and literature supply for genes mutated or copy number altered described in the main text as belonging to DDR or MErT/EMT categories. Table S3. List of somatic copy number aberrations identified by EXCAVATOR2. Table S4. List of gains with at least two copies and losses with homozygous deletions. Table S5. List of recurrent amplified genes on chromosome 1. Table S6. List of recurrently amplified genes of chromosome 1, belonging to the two significant biological functions identified by Ingenuity Pathway Analysis (IPA®, Qiagen; Bioinformatics, Redwood City, CA, USA; http://www.qiagen.com/ingenuity ). For the entire name of the genes, reported as gene ID, see Table S5. Table S7. List of recurrent deleted genes on chromosome 6.

IPA results and Sphingolipid signaling pathway DEGs with potential targeting differentially expressed serum exosomal miRNAs.

Sheet A in S3 Table (S6C Fig): Gene list (with RNA-Seq values) used for pathway analysis in S6C Fig. Cellular genes upregulated in both glucose/IFN-γ/mock infected and glucose/IFN-γ/vaccinia virus-infected relative to matching conditions but in galactose media. Comparison made using DEGs from 1) glucose/galactose A549 cultured cells treated with IFN-γ and mock infected with 2) glucose/galactose A549 cultured cells treated with IFN-γ and infected with vaccinia virus (MOI = 0.01). RNA-Seq data were curated with the cut-off values of log2 fold-change ≥ 1 or log2 fold-change ≤ -1 with an adjusted p-value ≤ 0.01 under mock infection condition. Data were analyzed with QIAGEN Ingenuity Pathway Analysis (QIAGEN IPA). Sheet B in S3 Table (S6C Fig): Gene list (with RNA-Seq values) used for pathway analysis in S6C Fig. Cellular genes upregulated in both galactose/IFN-γ/mock infected and galactose/IFN-γ/vaccinia virus infected relative to matching conditions but in glucose media. Comparison made using DEGs from 1) glucose/galactose A549 cultured cells treated with IFN-γ and mock infected with 2) glucose/galactose A549 cultured cells treated with IFN-γ and infected with vaccinia virus (M.O.I. = 0.01). RNA-Seq data were curated with the cut-off values of log2 fold-change ≥ 1 or log2 fold-change ≤ -1 with an adjusted p-value ≤ 0.01 under mock infection condition. Data were analyzed with QIAGEN Ingenuity Pathway Analysis (QIAGEN IPA). (XLSX)

Five oppositely-regulated differentially expressed genes between GBMvsCTRL and MSvsCTRL contrasts.

LimmaVoom mRNA differential expression analysis results.

An excel file with the top six networks from frameshift and non-frameshift INDELs exclusively from the layer or broiler line. QIAGEN’s Ingenuity® Pathway Analysis (IPA®) software ( http://www.ingenuity.com/ ) with default parameters were used to find metabolic pathways of genes with relevant biological functions. (XLSX 10 kb)