Peripheral Blood Gene Expression Profiles Research Articles

Blood transcriptional signatures of gallstone pathogenesis in patients with cholelithiasis and choledocholithiasis Running Title: Blood gene expression gallstone

Objective Disease-specific molecular signalling in peripheral blood has the potential to inform on pathophysiological mechanisms of diseases. Here, we aimed to investigate the blood-based gene expression profiles that reflect disease-specific pathogenic mechanisms of gallstones, such as antioxidant defence, heat-shock responses, DNA damage and repair, ABC pump mechanisms, mucin signals, mitochondrial dysfunction and apoptosis in patients with Cholelithiasis (CHL) and Choledocholithiasis (CHDL) Materials and Methods The relative fold-change in the mRNA expression levels of 73 genes were analysed using a gallstone-related qRT-PCR array in 10 control individuals, 24 CHL patients and 23 CHDL patients. Serum malondialdehyde levels was determined by thiobarbituric-acid reactive substance assays, and 8-hydroxy-2’-deoxyguanosine levels was determined by ELISA assays. Results Our results showed that peripheral whole blood gene expression profile strongly reflects tissue specific molecular signalling in gallstone pathogenesis. The present findings of altered gene expressions including antioxidant defence (CuZn-SOD, CAT), heat shock protein (HSP70), DNA repair (MLH1 and RAD18), pro-apoptotic (P53, BAX, Cyt-c and Caspase 3), ABC transporter (ABCB1, ABCC2, and LRP1) and Mucin signals (MUCIN4, MUCIN5AC and MUCIN5B) point out to DNA damages via oxidative stress as well as deteriorating ABC types pump mechanisms and mucin signals in CHL and CHDL patients. Our findings may also suggest that activation of mitophagy activator, DRP1/hFIS1/PINK1 axis induced by oxidative stress and DNA damage may have a role in the pathogenesis of CHL and CHDL. Conclusions Our results indicate that a blood-based gene expression signature has promising accuracy for monitorize pathogenesis of disease in CHL patients, CHDL patients and unaffected controls.

Transcriptomic score: A biomarker index of Alzheimer’s risk and resilience

AbstractBackgroundOne definition of resilience is the absence of disease despite high levels of risk. And in most cases, risk itself is best quantified by a multivariate profile rather than a single risk factor. Thus, multivariate risk profiles may be most appropriate for defining risk and resilience. RNA expression may serve as a risk biomarker, but no method has been proposed to build a whole‐transcriptome risk index of Alzheimer’s disease (AD) that could identify resilient individuals at the highest transcriptomic risk. Here, we derived a novel method for transcriptome‐wide gene expression data called “transcriptomic risk score” (TRS). We also employed the algorithm to build a “transcriptomic resilience score” (TReS) that offsets that risk.MethodAnalogous to a polygenic risk score, TRS is a weighted rank‐sum statistic that represents a subject’s relative expression level across the transcriptome, where the standardized beta coefficients are used as weights. We leveraged the TRS approach on peripheral blood gene expression profiles of 364 normal controls and 188 AD cases from ADNI and AddNeuroMed batch 1. Using TRSas a means of risk quantification, we then identified resilient normal controls as those with the highest percentiles of risk. High‐risk AD cases with matched TRSs range were retained for comparison. Expression levels of individual transcripts not included in the best‐fit TRS were contrasted between the high‐risk normal controls and high‐risk AD cases using mixed‐effect linear models. A TReS was derived using the residual transcriptomic variation that buffers against the transcriptomic risk. The association between TReS and case status was evaluated using a logistic regression model.ResultThe derived TReS differentiated resilient subjects (high‐risk normal controls) and high‐risk AD cases in independent replication samples (216 normal controls and 223 AD cases), with significant differentiation seen at all evaluated p‐value bins and a maximum of 3% variance accounted for at a p‐value threshold of 0.1.ConclusionThis analysis indicated that TRS and TReS analysis of blood has the potential of facilitating the understanding of phenotypic variability and improving clinically useful biomarker profiles for AD.

Effects of Deoxynivalenol and Fumonisins Fed in Combination to Beef Cattle: Immunotoxicity and Gene Expression.

We evaluated the effects of a treatment diet contaminated with 1.7 mg deoxynivalenol and 3.5 mg fumonisins (B1, B2 and B3) per kg ration on immune status and peripheral blood gene expression profiles in finishing-stage Angus steers. The mycotoxin treatment diet was fed for a period of 21 days followed by a two-week washout period during which time all animals consumed the control diet. Whole-blood leukocyte differentials were performed weekly throughout the experimental and washout period. Comparative profiles of CD4+ and CD8+ T cells, along with bactericidal capacity of circulating neutrophils and monocytes were evaluated at 0, 7, 14, 21 and 35 days. Peripheral blood gene expression was measured at 0, 7, 21 and 35 days via RNA sequencing. Significant increases in the percentage of CD4−CD8+ T cells were observed in treatment-fed steers after two weeks of treatment and were associated with decreased CD4:CD8 T-cell ratios at this same timepoint (p ≤ 0.10). No significant differences were observed as an effect of treatment in terms of bactericidal capacity at any timepoint. Dietary treatments induced major changes in transcripts associated with endocrine, metabolic and infectious diseases; protein digestion and absorption; and environmental information processing (inhibition of signaling and processing), as evaluated by dynamic impact analysis. DAVID analysis also suggested treatment effects on oxygen transport, extra-cellular signaling, cell membrane structure and immune system function. These results indicate that finishing-stage beef cattle are susceptible to the immunotoxic and transcript-inhibitory effects of deoxynivalenol and fumonisins at levels which may be realistically encountered in feedlot situations.

Transcriptomic profiling of peripheral blood mononuclear cells in patients with multiple primary melanoma.

e21558 Background: Patients (PTs) with cutaneous melanoma (MEL) are at increased risk of developing subsequent MEL. The primary objective of the present study was to determine whether the peripheral blood gene expression profile (GEP) of PTs with multiple primary melanoma (MPM) differs from that of patients with single primary melanoma (SPM) and healthy controls (HCs). A secondary objective is to assess if the peripheral blood immune transcriptional profile varies in relation with clinical and immunopathological characteristics of MEL. Methods: Cases with MPM, SPM and HCs were retrieved from the University of Pittsburgh Human Biological Sample and Nevus Image Banking and Analysis Repository. NanoString PanCancer Immune Profiling of 770-immune related genes was used to define GEP in peripheral blood mononuclear cells (PBMCs). Quality control (QC), data normalization, and differential expression analyses were performed using the ROSALIND Analysis System. Genes with fold change ≥1.5 and p-value < 0.05 were considered as differentially-expressed between cohorts. Results: In total, 47 PTs including 10 with SPM, 40 with MPM, and 10 HCs, were included. GEP data from 3 cases did not meet QC requirements and were discarded. The majority of PTs had stage I/II disease (87.2%). The median Breslow thickness was 1 mm, ulceration was present in 6 (12.2%) MELs. Tumor-infiltrating lymphocytes (TILs) were absent, non-brisk, and brisk in 20%, 65%, 15% of cases, respectively. The expression of IFIT1, MX1, ISG15 was decreased in PTs with MPM in comparison to HCs. In comparison to SPMs, PTs with MPM demonstrated downregulation of B-cell receptor signaling pathway genes (CD19, CD22, CD79A, BLK). 14 genes involved in TNF-alpha effects on cytokine activity, cell motility, and apoptosis, and 12 genes in the IL18 signaling pathway were differentially expressed in PTs with Breslow thickness of >1mm vs ≤1mm. The top downregulated genes in >1mm thick MELs were CD83, VEGFA and CXCL3. Ulcerated MEL had decreased expression of CCL2, CCL3, CCL3L1, CCL4. Stage III MEL had decreased expression of GZMH, MHC Class II genes HLA-DMB, HLA-DRB3, lectin receptor family genes KLRC2 and KLRD1, as well as Type-1 immune-associated biomarkers IRF1, SYK, TBX21 in comparison to Stage I/II cases. Tumors deficient in TILs exhibited decreased expression of Granzyme A and B, T-box 21, and IL18RAP (enhances IL18R-mediated signaling) in peripheral blood. Tumors with brisk TILs had increased expression of the B-cell associated transcripts CD19, CD22, CD79A in peripheral blood. Conclusions: The GEP of PBMC's from PTs with MPM was similar overall to that of PTs with SPM/HC, with only a subset of genes exhibiting differentially expression in PTs with MPM vs. PTs with SPM and HCs. GEP varied in relation to Breslow thickness, stage, ulceration, and TIL content in MELs.

Immune profiles of male giant panda (Ailuropoda melanoleuca) during the breeding season

BackgroundThe giant panda (Ailuropoda melanoleuca) is a threatened endemic Chinese species and a flagship species of national and global conservation concern. Life history theory proposes that reproduction and immunity can be mutually constraining and interrelated. Knowledge of immunity changes of male giant pandas during the breeding season is limited.ResultsHere, we researched peripheral blood gene expression profiles associated with immunity. Thirteen captive giant pandas, ranging from 9 to 11 years old, were divided into two groups based on their reproductive status. We identified 318 up-regulated DEGs and 43 down-regulated DEGs, which were enriched in 87 GO terms and 6 KEGG pathways. Additionally, we obtained 45 immune-related genes with altered expression, mostly up-regulated, and identified four hub genes HSPA4, SUGT1, SOD1, and IL1B in PPI analysis. These 45 genes were related to pattern recognition receptors, autophagy, peroxisome, proteasome, natural killer cell, antigen processing and presentation. SUGT1 and IL1B were related to pattern recognition receptors. HSP90AA1 was the most up-regulated gene and is a member of heat shock protein 90 family. HSP90 contributes to the translocation of extracellular antigen. KLRD1 encodes CD94, whose complex is an inhibitor of the cytotoxic activity of NK cells, was down-regulated. IGIP, which has the capability of inducing IgA production by B cells, was down-regulated, suggesting low concentration of IgA in male giant pandas. Our results suggest that most immune-related genes were up-regulated and more related to innate immune than adaptive immune.ConclusionsOur results indicated that breeding male giant pandas presented an immunoenhancement in innate immunity, enhanced antigen presentation and processing in cellular immunity compared to non-breeding males. The humoral immunity of male giant pandas may show a tendency to decrease during the breeding season. This study will provide a foundation for further studies of immunity and reproduction in male giant pandas.

Correlation and integration of circulating miRNA and peripheral whole blood gene expression profiles in patients with venous thromboembolism

ABSTRACT The main aim of this work was to evaluate differential expression and biological functions of circulating miRNA and whole peripheral blood (PB) genes in patients affected by venous thromboembolism (VTE) and in healthy subjects. Circulating miRNA sequences and PB expression profiles were obtained from GEO datasets. Ten miRNAs with the most significant differential expression rate (dif-miRNA) were subjected to miRbase to confirm their identity. Dif-miRNA targets were predicted by TargetScan and aligned with differentially expressed genes to obtain overlapping co-genes. Biological functions of co-genes were analyzed by Gene Ontology and KEGG analysis. Interaction network of dif-miRNAs, co-genes, and their downstream pathways were studied by analyzing protein–protein interaction (PPI) clusters (STRING) and determining the crucial hubs (Cytoscape). MiR-522-3p and miR-134 dif-miRNAs are involved in protein translation and apoptosis by regulating their respective co-genes in PB. Co-genes are present in nucleolus and extracellular exosomes and are involved in oxidative phosphorylation and ribosome/poly(A)-RNA organization. The predicted PPI network covered 107 clustered genes and 220 marginal joints, where ten hub genes participating in PPIs were found. All these hub genes were down-regulated in VTE patients. Our study identifies new miRNAs as potential biological markers and therapeutic targets for VTE.

Revising the Diagnosis of Idiopathic Uveitis by Peripheral Blood Transcriptomics

To test the hypothesis that idiopathic uveitis can be categorized into subtypes based on gene expression from blood. Case control study. We applied RNA-Seq to peripheral blood from patients with uveitis associated with 1 of 4 systemic diseases, including axial spondyloarthritis (n= 17), sarcoidosis (n= 13), inflammatory bowel disease (n= 12), tubulo-interstitial nephritis with uveitis (n= 10), or idiopathic uveitis (n= 38) as well as 18 healthy control subjects evaluated predominantly at Oregon Health and Science University. A high-dimensional negative binomial regression model implemented in the edgeR R package compared each disease group with the control subjects. The 20 most distinctive genes for each diagnosis were extracted. Of 80 genes, there were 75 unique genes. A classification algorithm was developed by fitting a gradient boosting tree with 5-fold cross-validation. Messenger RNA from subjects with idiopathic uveitis were analyzed to see if any fit clinically and by gene expression pattern with one of the diagnosable entities. For uveitis associated with a diagnosable systemic disease, gene expression profiling achieved an overall accuracy of 85% (balanced average of sensitivity plus specificity, P < .001). Although most patients with idiopathic uveitis presumably have none of these 4 associated systemic diseases, gene expression profiles helped to reclassify 11 of 38 subjects. Peripheral blood gene expression profiling is a potential adjunct in accurate differential diagnosis of the cause of uveitis. Validation of these results and characterization of the gene expression profile from additional discrete diagnoses could enhance the value of these observations.

Use of peripheral blood transcriptomic biomarkers to distinguish high-grade cervical squamous intraepithelial lesions from low-grade lesions.

It is crucial to classify cervical lesions into high-grade squamous intraepithelial lesions (HSILs) and low-grade SILs (LSILs), as LSILs are conservatively treated by observation, based on an expectation of natural regression, whereas HSILs usually require electrosurgical excision. In the present study, peripheral blood gene expression profiles were analyzed to identify transcriptomic biomarkers distinguishing HSILs from LSILs. A total of 102 blood samples were collected from women with cervical SILs (66 HSIL and 36 LSIL) for microarray hybridization. Candidate gene signatures were identified using AdaBoost algorithms, and a predictive model was constructed using logistic regression to differentiate HSILs from LSILs. To correct for possible bias as a result of the limited sample size and to verify the stability of the predictive model, a two-fold cross validation and null set analysis was conducted over 1,000 iterations. The functions of the transcriptomic biomarkers were then analyzed to elucidate the pathogenesis of cervical SIL. A total of 10 transcriptomic genes (STMN3, TRPC4AP, DYRK2, AGK, KIAA0319L, GRPEL1, ZFC3H1, LYL1, ITGB1 and ARHGAP18) were identified. The predictive model based on the 10-gene panel exhibited well-discriminated power. A cross validation process using known disease status exhibited almost the same performance as that of the predictive model, whereas null-set analysis with randomly reassigned disease status exhibited much lower predictive performance for distinguishing HSILs from LSILs. These biomarkers were involved in the ‘Rho GTPase cycle’, ‘mitochondrial protein import’, ‘oncogenic MAPK signaling’, ‘integrin cell surface interaction’ and ‘signaling by BRAF and RAF fusions’. In conclusion, peripheral blood gene expression analysis is a promising method for distinguishing HSILs from LSILs. The present study proposes 10 candidate genes that could be used in the future as diagnostic biomarkers and potential therapeutic targets for cervical SILs. A simple, non-invasive blood test would be clinically useful in the diagnosis and classification of patients with cervical SILs.

Whole blood transcriptome profile at hospital admission discriminates between patients with ST-segment elevation and non-ST-segment elevation acute myocardial infarction

Whether ST-segment (STEMI) and non-ST-segment elevation myocardial infarction (NSTEMI) should be regarded as distinct pathophysiological entities is a matter of debate. We tested the hypothesis that peripheral blood gene-expression profiles at presentation distinguish STEMI from NSTEMI. We performed a case-control study collecting whole-blood from 60 STEMI and 58 NSTEMI (defined according to the third universal definition of MI) consecutive patients on hospital admission. We used RNA-sequencing for the discovery phase, comparing 15 STEMI vs. 15 NSTEMI patients, matched for age, sex, and cardiovascular risk factors, and quantitative PCR in the remaining unmatched patients for validating top-significant genes. Gene-level differential expression analysis identified significant differences in the expression of 323 genes: 153 genes withstood correction for admission cardiac troponin I (cTnI), differentiating the two conditions independently of myocardial necrosis extent. Functional annotation analysis uncovered divergent modulation in leukocyte and platelet activation, cell migration, and mitochondrial respiratory processes. Linear regression analysis revealed gene expression patterns on admission predicting infarct size, as indexed by cTnI peak (R2 = 0.58–0.75). Our results unveil distinctive pathological traits for these two MI subtypes and provide insights into the early assessment of injury extent. This could translate into RNA-based disease-specific biomarkers for precision diagnosis and risk stratification.

Serial Peripheral Blood Gene Expression Profiling to Assess Immune Quiescence in Kidney Transplant Recipients with Stable Renal Function.

BackgroundTruGraf is a blood-based biomarker test that measures differential expression of a collection of genes that have been shown to correlate with surveillance biopsy results. However, in the majority of U.S. transplant centers, surveillance biopsies are not performed. The objectives of this study were to evaluate the clinical validity of TruGraf in stable kidney transplant recipients and to demonstrate the potential clinical utility of serial TruGraf testing in a center not utilizing surveillance biopsies.Material/MethodsSerum creatinine levels, TruGraf testing at multiple time points, and subsequent clinical follow-up were obtained for 28 patients.ResultsOverall concordance of TruGraf results, when compared with independent clinical assessment of testing, was 77% (54/70) for all tests; 79% (22/28) for test 1, 75% (21/28) for test 2, and 79% (11/14) for test 3. The negative predictive value (NPV) was 98.0%. Analysis of clinical utility indicated that 77% of TruGraf results would have been useful in patient management.ConclusionsOur results indicate the value of serial TruGraf testing in those transplant centers that do not perform surveillance biopsies as part of their standard of care. The high negative predictive value indicates the ability of TruGraf to confirm immune quiescence with a high degree of probability in patients with a Transplant eXcellence (TX) result, without the need to perform a surveillance biopsy.

Downregulation of peripheral PTGS2/COX-2 in response to valproate treatment in patients with epilepsy

Antiepileptic drug therapy has significant inter-patient variability in response towards it. The current study aims to understand this variability at the molecular level using microarray-based analysis of peripheral blood gene expression profiles of patients receiving valproate (VA) monotherapy. Only 10 unique genes were found to be differentially expressed in VA responders (n = 15) and 6 genes in the non-responders (n = 8) (fold-change >2, p < 0.05). PTGS2 which encodes cyclooxygenase-2, COX-2, showed downregulation in the responders compared to the non-responders. PTGS2/COX-2 mRNA profiles in the two groups corresponded to their plasma profiles of the COX-2 product, prostaglandin E2 (PGE2). Since COX-2 is believed to regulate P-glycoprotein (P-gp), a multidrug efflux transporter over-expressed at the blood-brain barrier (BBB) in drug-resistant epilepsy, the pathway connecting COX-2 and P-gp was further explored in vitro. Investigation of the effect of VA upon the brain endothelial cells (hCMEC/D3) in hyperexcitatory conditions confirmed suppression of COX-2-dependent P-gp upregulation by VA. Our findings suggest that COX-2 downregulation by VA may suppress seizure-mediated P-gp upregulation at the BBB leading to enhanced drug delivery to the brain in the responders. Our work provides insight into the association of peripheral PTGS2/COX-2 expression with VA efficacy and the role of COX-2 as a potential therapeutic target for developing efficacious antiepileptic treatment.

Gene Expression Profiles in Children With Suspected Sepsis

Sepsis recognition is a clinical challenge in children. We aim to determine whether peripheral blood gene expression profiles are associated with pathogen type and sepsis severity in children with suspected sepsis. This was a prospective pilot observational study in a tertiary pediatric emergency department with a convenience sample of children enrolled. Participants were older than 56 days and younger than 18 years, had suspected sepsis, and had not received broad-spectrum antibiotics in the previous 4 hours. Primary outcome was source pathogen, defined as confirmed bacterial source from sterile body fluid or confirmed viral source. Secondary outcome was sepsis severity, defined as maximum therapy required for shock reversal in the first 3 hospital days. We drew peripheral blood for ribonucleic acid isolation at the sepsis protocol activation, obtained gene expression measures with the GeneChip Human Gene 2.0 ST Array, and conducted differential expression analysis. We collected ribonucleic acid samples from a convenience sample of 122 children with suspected sepsis and 12 healthy controls. We compared the 66 children (54%) with confirmed bacterial or viral infection and found 558 differentially expressed genes, many related to interferon signaling or viral immunity. We did not find statistically significant gene expression differences in patients according to sepsis severity. The study demonstrates feasibility of evaluating gene expression profiling data in children evaluated for sepsis in the pediatric emergency department setting. Our results suggest that gene expression profiling may facilitate identification of source pathogen in children with suspected sepsis, which could ultimately lead to improved tailoring of sepsis treatment and antimicrobial stewardship.

An age-independent gene signature for monitoring acute rejection in kidney transplantation.

Acute rejection (AR) remains a significant problem that negatively impacts long-term renal allograft survival. Numerous therapies are used to prevent AR that differ by center and recipient age. This variability confounds diagnostic methods.Methods: To develop an age-independent gene signature for AR effective across a broad array of immunosuppressive regimens, we compiled kidney transplant biopsy (n=1091) and peripheral blood (n=392) gene expression profiles from 12 independent public datasets. After removing genes differentially expressed in pediatric and adult patients, we compared gene expression profiles from biopsy and peripheral blood samples of patients with AR to those who were stable (STA), using Mann-Whitney U Tests with validation in independent testing datasets. We confirmed this signature in pediatric and adult patients (42 AR and 47 STA) from our institutional biorepository.Results: We identified a novel age-independent gene network that identified AR from both kidney and blood samples. We developed a 90-probe set signature targeting 76 genes that differentiated AR from STA and found an 8 gene subset (DIP2C, ENOSF1, FBXO21, KCTD6, PDXDC1, REXO2, HLA-E, and RAB31) that was associated with AR.Conclusion: We used publicly available datasets to create a gene signature of AR that identified AR irrespective of immunosuppression regimen or recipient age. This study highlights a novel model to screen and validate biomarkers across multiple treatment regimens.

RNA-Seq analysis of ileocecal valve and peripheral blood from Holstein cattle infected with Mycobacterium avium subsp. paratuberculosis revealed dysregulation of the CXCL8/IL8 signaling pathway

Paratuberculosis is chronic granulomatous enteritis of ruminants caused by Mycobacterium avium subsp. paratuberculosis (MAP). Whole RNA-sequencing (RNA-Seq) is a promising source of novel biomarkers for early MAP infection and disease progression in cattle. Since the blood transcriptome is widely used as a source of biomarkers, we analyzed whether it recapitulates, at least in part, the transcriptome of the ileocecal valve (ICV), the primary site of MAP colonization. Total RNA was prepared from peripheral blood (PB) and ICV samples, and RNA-Seq was used to compare gene expression between animals with focal or diffuse histopathological lesions in gut tissues versus control animals with no detectable signs of infection. Our results demonstrated both shared, and PB and ICV-specific gene expression in response to a natural MAP infection. As expected, the number of differentially expressed (DE) genes was larger in the ICV than in the PB samples. Among the DE genes in the PB and ICV samples, there were some common genes irrespective of the type of lesion including the C-X-C motif chemokine ligand 8 (CXCL8/IL8), apolipoprotein L (APOLD1), and the interferon inducible protein 27 (IFI27). The biological processes (BP) enriched in the PB gene expression profiles from the cows with diffuse lesions included the killing of cells of other organism, defense response, immune response and the regulation of neutrophil chemotaxis. Two of these BP, the defense and immune response, were also enriched in the ICV from the cows with diffuse lesions. Metabolic analysis of the DE genes revealed that the N-glycan biosynthesis, bile secretion, one-carbon pool by folate and purine metabolism were significantly enriched in the ICV from the cows with focal lesions. In the ICV from cows with diffuse lesions; the valine, leucine and isoleucine degradation route, purine metabolism, vitamin digestion and absorption and the cholesterol routes were enriched. Some of the identified DE genes, BP and metabolic pathways will be studied further to develop novel diagnostic tools, vaccines and immunotherapeutics.

Molecular signature of penumbra in acute ischemic stroke: a pilot transcriptomics study.

We aimed to characterize peripheral blood gene expression profile of penumbra defined as MRI perfusion–diffusion mismatch (PD MM) in peripheral blood of patients with acute ischemic stroke. We studied 23 patients. Perfusion–diffusion mismatch volume was observed to be associated and significantly correlated with the expression of 34 genes including those related to inflammation, SUMOylation, and coagulation; while lipopolysaccharide inhibition was identified to be a candidate upstream regulator of these processes (z‐score −2.38, P = 0.04). Penumbral volume is correlated with a specific gene expression profile in the peripheral blood characterized by overlap of inflammatory and neuroprotective pathways that are regulated by lipopolysaccharide inhibition.

Regulators of Salmonella-host interaction identified by peripheral blood transcriptome profiling: roles of TGFB1 and TRP53 in intracellular Salmonella replication in pigs

Peripheral blood transcriptome is an important intermediate data source for investigating the mechanism of Salmonella invasion, proliferation, and transmission. We challenged 4-week old piglets with Salmonella enterica serovar Typhimurium LT2 and investigated the peripheral blood gene expression profile before treatment (d0) and at 2 and 7 days post-inoculation (dpi) using deep sequencing. Regulator pathways were first predicted in silico and validated by wet-lab experiments. In total, 1255, 765, and 853 genes were differentially expressed between 2 dpi/d0, 7 dpi/d0, and 7 dpi/2 dpi, respectively. Additionally, 1333 genes showed a time effect during the investigated Salmonella infection period. Clustering analysis showed that the differentially expressed genes fell into six distinct expression clusters. Pathway annotation of these gene clusters showed that the innate immune system was first significantly upregulated at 2 dpi and then attenuated at 7 dpi. Toll-like receptor cascades, MyD88 cascade, phagosome pathway, cytokine signaling pathway, and lysosome pathway showed a similar expression pattern. Interestingly, we found that the ribosome pathway was significantly inhibited at 2 and 7 dpi. Gene expression regulation network enrichment analysis identified several candidate factors controlling the expression clusters. Further in vitro study showed that TGFB1 can inhibit Salmonella replication whereas TRP53 can promote Salmonella replication in porcine peripheral blood mononuclear cells and murine macrophages. These results provide new insights into the molecular mechanism of Salmonella-host interactions and clues for the genetic improvement of Salmonella infection resistance in pigs.

A crowdsourced analysis to identify ab initio molecular signatures predictive of susceptibility to viral infection

The response to respiratory viruses varies substantially between individuals, and there are currently no known molecular predictors from the early stages of infection. Here we conduct a community-based analysis to determine whether pre- or early post-exposure molecular factors could predict physiologic responses to viral exposure. Using peripheral blood gene expression profiles collected from healthy subjects prior to exposure to one of four respiratory viruses (H1N1, H3N2, Rhinovirus, and RSV), as well as up to 24 h following exposure, we find that it is possible to construct models predictive of symptomatic response using profiles even prior to viral exposure. Analysis of predictive gene features reveal little overlap among models; however, in aggregate, these genes are enriched for common pathways. Heme metabolism, the most significantly enriched pathway, is associated with a higher risk of developing symptoms following viral exposure. This study demonstrates that pre-exposure molecular predictors can be identified and improves our understanding of the mechanisms of response to respiratory viruses.

Longitudinal transcriptomic characterization of the immune response to acute hepatitis C virus infection in patients with spontaneous viral clearance.

Most individuals exposed to hepatitis C virus (HCV) become persistently infected while a minority spontaneously eliminate the virus. Although early immune events influence infection outcome, the cellular composition, molecular effectors, and timeframe of the host response active shortly after viral exposure remain incompletely understood. Employing specimens collected from people who inject drugs (PWID) with high risk of HCV exposure, we utilized RNA-Seq and blood transcriptome module (BTM) analysis to characterize immune function in peripheral blood mononuclear cells (PBMC) before, during, and after acute HCV infection resulting in spontaneous resolution. Our results provide a detailed description of innate immune programs active in peripheral blood during acute HCV infection, which include prominent type I interferon and inflammatory signatures. Innate immune gene expression rapidly returns to pre-infection levels upon viral clearance. Comparative analyses using peripheral blood gene expression profiles from other viral and vaccine studies demonstrate similarities in the immune responses to acute HCV and flaviviruses. Of note, both acute dengue virus (DENV) infection and acute HCV infection elicit similar innate antiviral signatures. However, while transient in DENV infection, this signature was sustained for many weeks in the response to HCV. These results represent the first longitudinal transcriptomic characterization of human immune function in PBMC during acute HCV infection and identify several dynamically regulated features of the complex response to natural HCV exposure.

Each pregnancy linearly changes immune gene expression in the blood of healthy women compared with breast cancer patients

BackgroundThere is a large body of evidence demonstrating long-lasting protective effect of each full-term pregnancy (FTP) on the development of breast cancer (BC) later in life, a phenomenon that could be related to both hormonal and immunological changes during pregnancies. In this work, we studied the pregnancy-associated differences in peripheral blood gene expression profiles between healthy women and women diagnosed with BC in a prospective design.MethodsUsing an integrated system epidemiology approach, we modeled BC incidence as a function of parity in the Norwegian Women and Cancer (NOWAC) cohort (165,000 women) and then tested the resulting mathematical model using gene expression profiles in blood in a nested case–control study (460 invasive case–control pairs) of women from the NOWAC postgenome cohort. Lastly, we undertook a gene set enrichment analysis for immunological gene sets.ResultsA linear trend fitted the dataset precisely showing an 8% decrease in risk of BC for each FTP, independent of stratification on other risk factors and lasting for decades after a woman’s last FTP. Women with six children demonstrated 48% reduction in the incidence of BC compared to nulliparous. When we looked at gene expression, we found that 756 genes showed linear trends in cancer-free controls (false discovery rate [FDR] 5%), but this was not the case for any of the genes in BC cases. Gene set enrichment analysis of immunologic gene sets (C7 collection in Molecular Signatures Database) revealed 215 significantly enriched human gene sets (FDR 5%).ConclusionWe found marked differences in gene expression and enrichment profiles of immunologic gene sets between BC cases and healthy controls, suggesting an important protective effect of the immune system on BC risk.

A peripheral blood transcriptomic signature predicts autoantibody development in infants at risk of type 1 diabetes.

Autoimmune-mediated destruction of pancreatic islet β cells results in type 1 diabetes (T1D). Serum islet autoantibodies usually develop in genetically susceptible individuals in early childhood before T1D onset, with multiple islet autoantibodies predicting diabetes development. However, most at-risk children remain islet-antibody negative, and no test currently identifies those likely to seroconvert. We sought a genomic signature predicting seroconversion risk by integrating longitudinal peripheral blood gene expression profiles collected in high-risk children included in the BABYDIET and DIPP cohorts, of whom 50 seroconverted. Subjects were followed for 10 years to determine time of seroconversion. Any cohort effect and the time of seroconversion were corrected to uncover genes differentially expressed (DE) in seroconverting children. Gene expression signatures associated with seroconversion were evident during the first year of life, with 67 DE genes identified in seroconverting children relative to those remaining antibody negative. These genes contribute to T cell-, DC-, and B cell-related immune responses. Near-birth expression of ADCY9, PTCH1, MEX3B, IL15RA, ZNF714, TENM1, and PLEKHA5, along with HLA risk score predicted seroconversion (AUC 0.85). The ubiquitin-proteasome pathway linked DE genes and T1D susceptibility genes. Therefore, a gene expression signature in infancy predicts risk of seroconversion. Ubiquitination may play a mechanistic role in diabetes progression.