RNAseq of INOCA patients identifies innate, invariant, and acquired immune changes: potential autoimmune microvascular dysfunction.

Background: Ischemia with non-obstructive coronary arteries (INOCA) is a major clinical entity that involves potentially 20%-30% of patients with chest pain. INOCA is typically attributed either to coronary microvascular disease and/or vasospasm, but is likely distinct from classical coronary artery disease (CAD). Objectives: To gain insights into the etiology of INOCA and CAD, RNA sequencing of whole blood from patients undergoing both stress testing and elective invasive coronary angiography (ICA) was conducted. Methods: Stress testing and ICA of 177 patients identified 40 patients (23%) with INOCA compared to 39 controls (stress-, ICA-). ICA+ patients divided into 38 stress- and 60 stress+. RNAseq was performed by Illumina with ribosomal RNA depletion. Transcriptome changes were analyzed by DeSeq2 and curated by manual and automated methods. Results: Differentially expressed genes for INOCA were associated with elevated levels of transcripts related to mucosal-associated invariant T (MAIT) cells, plasmacytoid dendritic cells (pcDC), and memory B cells, and were associated with autoimmune diseases such as rheumatoid arthritis. Decreased transcripts were associated with neutrophils, but neutrophil transcripts, per se, were not less abundant in INOCA. CAD transcripts were more related to T cell functions. Conclusions: Elevated transcripts related to pcDC, MAIT, and memory B cells suggest an autoimmune component to INOCA. Reduced neutrophil transcripts are likely attributed to chronic activation leading to increased translation and degradation. Thus, INOCA could result from stimulation of B cell, pcDC, invariant T cell, and neutrophil activation that compromises cardiac microvascular function.

COVID-19 point-of-care tests can identify low-antibody individuals: In-depth immunoanalysis of boosting benefits in a healthy cohort.

The recommended COVID-19 booster vaccine uptake is low. At-home lateral flow assay (LFA) antigen tests are widely accepted for detecting infection during the pandemic. Here, we present the feasibility and potential benefits of using LFA-based antibody tests as a means for individuals to detect inadequate immunity and make informed decisions about COVID-19 booster immunization. In a health care provider cohort, we investigated the changes in the breadth and depth of humoral and T cell immune responses following mRNA vaccination and boosting in LFA-positive and LFA-negative antibody groups. We show that negative LFA antibody tests closely reflect the lack of functional humoral immunity observed in a battery of sophisticated immune assays, while positive results do not necessarily reflect adequate immunity. After booster vaccination, both groups gain depth and breadth of systemic antibodies against evolving SARS-CoV-2 and related viruses. Our findings show that LFA-based antibody tests can alert individuals about inadequate immunity against COVID-19, thereby increasing booster shots and promoting herd immunity.

Blood RNA Biomarkers Identify Bacterial and Biofilm Coinfections in COVID-19 Intensive Care Patients.

Purpose: Secondary opportunistic coinfections are a significant contributor to morbidity and mortality in intensive care unit (ICU) patients, but can be difficult to identify. Presently, new blood RNA biomarkers were tested in ICU patients to diagnose viral, bacterial, and biofilm coinfections. Methods: COVID-19 ICU patients had whole blood drawn in RNA preservative and stored at -80°C. Controls and subclinical infections were also studied. Droplet digital polymerase chain reaction (ddPCR) quantified 6 RNA biomarkers of host neutrophil activation to bacterial (DEFA1), biofilm (alkaline phosphatase [ALPL], IL8RB/CXCR2), and viral infections (IFI27, RSAD2). Viral titer in blood was measured by ddPCR for SARS-CoV2 (SCV2). Results: RNA biomarkers were elevated in ICU patients relative to controls. DEFA1 and ALPL RNA were significantly higher in severe versus incidental/moderate cases. SOFA score was correlated with white blood cell count (0.42), platelet count (-0.41), creatinine (0.38), and lactate dehydrogenase (0.31). ALPL RNA (0.59) showed the best correlation with SOFA score. IFI27 (0.52) and RSAD2 (0.38) were positively correlated with SCV2 viral titer. Overall, 57.8% of COVID-19 patients had a positive RNA biomarker for bacterial or biofilm infection. Conclusions: RNA biomarkers of host neutrophil activation indicate the presence of bacterial and biofilm coinfections in most COVID-19 patients. Recognizing coinfections may help to guide the treatment of ICU patients.

RNAseq profiling of blood from patients with coronary artery disease: Signature of a T cell imbalance.

Background: Cardiovascular disease had a global prevalence of 523 million cases and 18.6 million deaths in 2019. The current standard for diagnosing coronary artery disease (CAD) is coronary angiography either by invasive catheterization (ICA) or computed tomography (CTA). Prior studies employed single-molecule, amplification-independent RNA sequencing of whole blood to identify an RNA signature in patients with angiographically confirmed CAD. The present studies employed Illumina RNAseq and network co-expression analysis to identify systematic changes underlying CAD. Methods: Whole blood RNA was depleted of ribosomal RNA (rRNA) and analyzed by Illumina total RNA sequencing (RNAseq) to identify transcripts associated with CAD in 177 patients presenting for elective invasive coronary catheterization. The resulting transcript counts were compared between groups to identify differentially expressed genes (DEGs) and to identify patterns of changes through whole genome co-expression network analysis (WGCNA). Results: The correlation between Illumina amplified RNAseq and the prior SeqLL unamplified RNAseq was quite strong (r = 0.87), but there was only 9 % overlap in the DEGs identified. Consistent with the prior RNAseq, the majority (93 %) of DEGs were down-regulated ~1.7-fold in patients with moderate to severe CAD (>20 % stenosis). DEGs were predominantly related to T cells, consistent with known reductions in Tregs in CAD. Network analysis did not identify pre-existing modules with a strong association with CAD, but patterns of T cell dysregulation were evident. DEGs were enriched for transcripts associated with ciliary and synaptic transcripts, consistent with changes in the immune synapse of developing T cells. Conclusions: These studies confirm and extend a novel mRNA signature of a Treg-like defect in CAD. The pattern of changes is consistent with stress-related changes in the maturation of T and Treg cells, possibly due to changes in the immune synapse.

Longitudinal Analysis of Humoral and Cellular Immune Response Following SARS-CoV-2 Vaccination Supports Utilizing Point-Of-Care Tests to Enhance COVID-19 Booster Uptake.

Individuals with weaker neutralizing responses show reduced protection with SARS-CoV-2 variants. Booster vaccines are recommended for vaccinated individuals, but the uptake is low. We present the feasibility of utilizing point-of-care tests (POCT) to support evidence-based decision-making around COVID-19 booster vaccinations. Using infectious virus neutralization, ACE2 blocking, spike binding, and TCR sequencing assays, we investigated the dynamics of changes in the breadth and depth of blood and salivary antibodies as well as T-cell clonal response following mRNA vaccination in a cohort of healthcare providers. We evaluated the accuracy of two POCTs utilizing either blood or saliva to identify those in whom humoral immunity was inadequate. >4 months after two doses of mRNA vaccine, SARS-CoV-2 binding and neutralizing Abs (nAbs) and T-cell clones declined 40-80%, and 2/3rd lacked Omicron nAbs. After the third mRNA booster, binding and neutralizing Abs increased overall in the systemic compartment; notably, individuals with previously weak nAbs gained sharply. The third dose failed to stimulate secretory IgA, but salivary IgG closely tracked systemic IgG levels. Vaccine boosting increased Ab breadth against a divergent bat sarbecovirus, SHC014, although the TCR-beta sequence breadth was unchanged. Post 3rd booster dose, Ab avidity increased for the Wuhan and Delta strains, while avidity against Omicron and SHC014 increased to levels seen for Wuhan after the second dose. Negative results on POCTs strongly correlated with a lack of functional humoral immunity. The third booster dose helps vaccinees gain depth and breadth of systemic Abs against evolving SARS-CoV-2 and related viruses. Our findings show that POCTs are useful and easy-to-access tools to inform inadequate humoral immunity accurately. POCTs designed to match the circulating variants can help individuals with booster vaccine decisions and could serve as a population-level screening platform to preserve herd immunity. One Sentence Summary: SARS-CoV-2 point-of-care antibody tests are valuable and easy-to-access tools to inform inadequate humoral immunity and to support informed decision-making regarding the current and future booster vaccination.

Diagnostic accuracy of novel mRNA blood biomarkers of infection to predict outcomes in emergency department patients with undifferentiated abdominal pain.

Abdominal pain represents greater than 20% of US Emergency Department (ED) visits due to a wide range of illnesses. There are currently no reliable blood biomarkers to predict serious outcomes in patients with abdominal pain. Our previous studies have identified three mRNA transcripts related to innate immune activation: alkaline phosphatase (ALPL), interleukin-8 receptor-ß (IL8RB), and defensin-1 (DEFA1) as promising candidates to detect an intra-abdominal infection. The objective of this study was to evaluate the accuracy of these mRNA biomarkers to predict likely infection, hospitalization and surgery in Emergency Department patients with undifferentiated abdominal pain. We prospectively enrolled Emergency Department patients with undifferentiated abdominal pain who received an abdominal CT scan as part of their evaluation. Clinical outcomes were abstracted from the CT scan and medical records. mRNA biomarker levels were calculated independent of the clinical outcomes and their accuracy was assessed to predict infectious diagnoses, surgery and hospital admission. 89 patients were enrolled; 21 underwent surgery; 47 underwent hospital admission; and, no deaths were observed within 30 days. In identifying which cases were likely infectious, mRNA biomarkers' AUC values were: ALPL, 0.83; DEFA1 0.51; IL8RB, 0.74; and ALPL + IL8RB, 0.79. In predicting which Emergency Department patients would receive surgery, the AUC values were: ALPL, 0.75; DEFA1, 0.58; IL8RB, 0.75; and ALPL + IL8RB, 0.76. In predicting hospital admission, the AUC values were: ALPL, 0.78; DEFA1, 0.52; IL8RB, 0.74; and, ALPL + IL8RB, 0.77. For predicting surgery, ALPL + IL8RB's positive likelihood ratio (LR) was 3.97; negative LR (NLR) was 0.70. For predicting hospital admission, the same marker's positive LR was 2.80 with an NLR of 0.45. Where the primary cause for admission was a potentially infectious disorder, 33 of 34 cases (97%) had positive RNA scores. In a pragmatic, prospective diagnostic accuracy trial in Emergency Department patients with undifferentiated abdominal pain, mRNA biomarkers showed good accuracy to identify patients with potential infection, as well as those needing surgery or hospital admission.

RNA Sequencing in COVID-19 patients identifies neutrophil activation biomarkers as a promising diagnostic platform for infections.

Infection with the SARS-CoV2 virus can vary from asymptomatic, or flu-like with moderate disease, up to critically severe. Severe disease, termed COVID-19, involves acute respiratory deterioration that is frequently fatal. To understand the highly variable presentation, and identify biomarkers for disease severity, blood RNA from COVID-19 patient in an intensive care unit was analyzed by whole transcriptome RNA sequencing. Both SARS-CoV2 infection and the severity of COVID-19 syndrome were associated with up to 25-fold increased expression of neutrophil-related transcripts, such as neutrophil defensin 1 (DEFA1), and 3-5-fold reductions in T cell related transcripts such as the T cell receptor (TCR). The DEFA1 RNA level detected SARS-CoV2 viremia with 95.5% sensitivity, when viremia was measured by ddPCR of whole blood RNA. Purified CD15+ neutrophils from COVID-19 patients were increased in abundance and showed striking increases in nuclear DNA staining by DAPI. Concurrently, they showed >10-fold higher elastase activity than normal controls, and correcting for their increased abundance, still showed 5-fold higher elastase activity per cell. Despite higher CD15+ neutrophil elastase activity, elastase activity was extremely low in plasma from the same patients. Collectively, the data supports the model that increased neutrophil and decreased T cell activity is associated with increased COVID-19 severity, and suggests that blood DEFA1 RNA levels and neutrophil elastase activity, both involved in neutrophil extracellular traps (NETs), may be informative biomarkers of host immune activity after viral infection.

RNAseq analysis of treatment-dependent signaling changes during inflammation in a mouse cutaneous wound healing model.

BACKGROUND: Despite proven therapeutic effects in inflammatory conditions, the specific mechanisms of phytochemical therapies are not well understood. The transcriptome effects of Traumeel (Tr14), a multicomponent natural product, and diclofenac, a non-selective cyclooxygenase (COX) inhibitor, were compared in a mouse cutaneous wound healing model to identify both known and novel pathways for the anti-inflammatory effect of plant-derived natural products. METHODS: Skin samples from abraded mice were analyzed by single-molecule, amplification-free RNAseq transcript profiling at 7 points between 12 and 192 h after injury. Immediately after injury, the wounds were treated with either diclofenac, Tr14, or placebo control (n = 7 per group/time). RNAseq levels were compared between treatment and control at each time point using a systems biology approach. RESULTS: At early time points (12-36 h), both control and Tr14-treated wounds showed marked increase in the inducible COX2 enzyme mRNA, while diclofenac-treated wounds did not. Tr14, in contrast, modulated lipoxygenase transcripts, especially ALOX12/15, and phospholipases involved in arachidonate metabolism. Notably, Tr14 modulated a group of cell-type specific markers, including the T cell receptor, that could be explained by an overarching effect on the type of cells that were recruited into the wound tissue. CONCLUSIONS: Tr14 and diclofenac had very different effects on the COX/LOX synthetic pathway after cutaneous wounding. Tr14 allowed normal autoinduction of COX2 mRNA, but suppressed mRNA levels for key enzymes in the leukotriene synthetic pathway. Tr14 appeared to have a broad 'phytocellular' effect on the wound transcriptome by altering the balance of cell types present in the wound.

RNA sequencing of blood in coronary artery disease: involvement of regulatory T cell imbalance.

BACKGROUND: Cardiovascular disease had a global prevalence of 523 million cases and 18.6 million deaths in 2019. The current standard for diagnosing coronary artery disease (CAD) is coronary angiography. Surprisingly, despite well-established clinical indications, up to 40% of the one million invasive cardiac catheterizations return a result of 'no blockage'. The present studies employed RNA sequencing of whole blood to identify an RNA signature in patients with angiographically confirmed CAD. METHODS: Whole blood RNA was depleted of ribosomal RNA (rRNA) and analyzed by single-molecule sequencing of RNA (RNAseq) to identify transcripts associated with CAD (TRACs) in a discovery group of 96 patients presenting for elective coronary catheterization. The resulting transcript counts were compared between groups to identify differentially expressed genes (DEGs). RESULTS: Surprisingly, 98% of DEGs/TRACs were down-regulated ~ 1.7-fold in patients with mild to severe CAD (> 20% stenosis). The TRACs were independent of comorbid risk factors for CAD, such as sex, hypertension, and smoking. Bioinformatic analysis identified an enrichment in transcripts such as FoxP1, ICOSLG, IKZF4/Eos, SMYD3, TRIM28, and TCF3/E2A that are likely markers of regulatory T cells (Treg), consistent with known reductions in Tregs in CAD. A validation cohort of 80 patients confirmed the overall pattern (92% down-regulation) and supported many of the Treg-related changes. TRACs were enriched for transcripts associated with stress granules, which sequester RNAs, and ciliary and synaptic transcripts, possibly consistent with changes in the immune synapse of developing T cells. CONCLUSIONS: These studies identify a novel mRNA signature of a Treg-like defect in CAD patients and provides a blueprint for a diagnostic test for CAD. The pattern of changes is consistent with stress-related changes in the maturation of T and Treg cells, possibly due to changes in the immune synapse.

Biomarker discovery in attention deficit hyperactivity disorder: RNA sequencing of whole blood in discordant twin and case-controlled cohorts.

BACKGROUND: A variety of DNA-based methods have been applied to identify genetic markers of attention deficit hyperactivity disorder (ADHD), but the connection to RNA-based gene expression has not been fully exploited. METHODS: Using well defined cohorts of discordant, monozygotic twins from the Michigan State University Twin Registry, and case-controlled ADHD cases in adolescents, the present studies utilized advanced single molecule RNA sequencing to identify expressed changes in whole blood RNA in ADHD. Multiple analytical strategies were employed to narrow differentially expressed RNA targets to a small set of potential biomarkers of ADHD. RESULTS: RNA markers common to both the discordant twin study and case-controlled subjects further narrowed the putative targets, some of which had been previously associated with ADHD at the DNA level. The potential role of several differentially expressed genes, including ABCB5, RGS2, GAK, GIT1 and 3 members of the galactose metabolism pathway (GALE, GALT, GALK1) are substantiated by prior associations to ADHD and by established mechanistic connections to molecular pathways relevant to ADHD and behavioral control. CONCLUSIONS: The convergence of DNA, RNA, and metabolic data suggests these may be promising targets for diagnostics and therapeutics in ADHD.

Benchmark Evaluation of True Single Molecular Sequencing to Determine Cystic Fibrosis Airway Microbiome Diversity.

Cystic fibrosis (CF) is an autosomal recessive disease associated with recurrent lung infections that can lead to morbidity and mortality. The impact of antibiotics for treatment of acute pulmonary exacerbations on the CF airway microbiome remains unclear with prior studies giving conflicting results and being limited by their use of 16S ribosomal RNA sequencing. Our primary objective was to validate the use of true single molecular sequencing (tSMS) and PathoScope in the analysis of the CF airway microbiome. Three control samples were created with differing amounts of Burkholderia cepacia, Pseudomonas aeruginosa, and Prevotella melaninogenica, three common bacteria found in cystic fibrosis lungs. Paired sputa were also obtained from three study participants with CF before and >6 days after initiation of antibiotics. Antibiotic resistant B. cepacia and P. aeruginosa were identified in concurrently obtained respiratory cultures. Direct sequencing was performed using tSMS, and filtered reads were aligned to reference genomes from NCBI using PathoScope and Kraken and unique clade-specific marker genes using MetaPhlAn. A total of 180-518 K of 6-12 million filtered reads were aligned for each sample. Detection of known pathogens in control samples was most successful using PathoScope. In the CF sputa, alpha diversity measures varied based on the alignment method used, but similar trends were found between pre- and post-antibiotic samples. PathoScope outperformed Kraken and MetaPhlAn in our validation study of artificial bacterial community controls and also has advantages over Kraken and MetaPhlAn of being able to determine bacterial strains and the presence of fungal organisms. PathoScope can be confidently used when evaluating metagenomic data to determine CF airway microbiome diversity.

Deep Sequencing Transcriptome Analysis of Murine Wound Healing: Effects of a Multicomponent, Multitarget Natural Product Therapy-Tr14.

Wound healing involves an orchestrated response that engages multiple processes, such as hemostasis, cellular migration, extracellular matrix synthesis, and in particular, inflammation. Using a murine model of cutaneous wound repair, the transcriptome was mapped from 12 h to 8 days post-injury, and in response to a multicomponent, multi-target natural product, Tr14. Using single-molecule RNA sequencing (RNA-seq), there were clear temporal changes in known transcripts related to wound healing pathways, and additional novel transcripts of both coding and non-coding genes. Tr14 treatment modulated >100 transcripts related to key wound repair pathways, such as response to wounding, wound contraction, and cytokine response. The results provide the most precise and comprehensive characterization to date of the transcriptome's response to skin damage, repair, and multicomponent natural product therapy. By understanding the wound repair process, and the effects of natural products, it should be possible to intervene more effectively in diseases involving aberrant repair.

Predictors of high on-aspirin platelet reactivity in elderly patients with coronary artery disease.

OBJECTIVES: Previous studies have illustrated the link between high on-aspirin platelet reactivity (HAPR) with increasing thrombotic risks. The aim of our study was to investigate relative risk factors of HAPR in elderly patients with coronary artery disease. METHODS: Elderly, hospitalized coronary artery disease patients on regular aspirin treatment were enrolled from January 2014 to September 2016. Medical records of each patient were collected, including demographic information, cardiovascular risk factors, concomitant drugs and routine biological parameters. Arachidonic acid (AA, 0.5 mg/mL) and adenosine diphosphate (ADP, 5 µmol/L) induced platelet aggregation were measured via light transmission assay (LTA) to evaluate antiplatelet responses, referred as LTA-AA and LTA-ADP. RESULTS: A total of 275 elderly patients were included, with mean age of 77.2±8.1 years, and males accounted for 81.8%. HAPR was defined as LTA-AA in the upper quartile of the enrolled population. HAPR patients tended to have lower renal function (P=0.052). Higher serum uric acid (SUA) level, as well as lower platelet count, hemoglobin and hematocrit were observed in HAPR patients, with a higher proportion of diuretics use (P<0.05). Multivariate analysis revealed that SUA (OR: 1.004, 95% CI: 1.000-1.007, P=0.048), platelet count (OR: 0.994, 95% CI: 0.989-1.000, P=0.045), hematocrit (OR: 0.921, 95% CI: 0.864-0.981, P=0.011) and concomitant P2Y12 receptor inhibitors use (OR: 1.965, 95% CI: 1.075-3.592, P=0.028) were correlated with HAPR. Spearman's correlation analysis demonstrated an inverse association of LTA-AA with hematocrit (r=-0.234, P<0.001), hemoglobin (r=-0.209, P<0.001) and estimated glomerular filtration rate (r=-0.132, P=0.031). CONCLUSION: SUA, platelet count, hematocrit and P2Y12 receptor inhibitors use were independently correlated with HAPR. These parameters might provide novel therapeutic targets for optimizing antiplatelet therapy.

Validation of aspirin response-related transcripts in patients with coronary artery disease and preliminary investigation on CMTM5 function.

Aspirin is widely used in the prevention of cardiovascular diseases, but the antiplatelet responses vary from one patient to another. To validate aspirin response related transcripts and illustrate their roles in predicting cardiovascular events, we have quantified the relative expression of 14 transcripts previously identified as related to high on-aspirin platelet reactivity (HAPR) in 223 patients with coronary artery disease (CAD) on regular aspirin treatment. All patients were followed up regularly for cardiovascular events (CVE). The mean age of our enrolled population was 75.80±8.57years. HAPR patients showed no significant differences in terms of co-morbidities and combined drugs. Besides, the relative expression of HLA-DQA1 was significantly lower in low on-aspirin platelet reactivity (LAPR) patients, when compared with HAPR and high normal (HN) group (p=0.028). What's more, the number of arteries involved, HAPR status and the relative expression of CLU, CMTM5 and SPARC were independent risk factors for CVE during follow up (p<0.05). In addition, overexpression of CMTM5 attenuated endothelial cells (ECs) migration and proliferation, with significantly decreased phosphorylated-Akt levels, while its inhibition promoted these processes in vitro (p<0.05).Our study provides evidence that circulating transcripts might be potential biomarkers in predicting cardiovascular events. CMTM5 might exert anti-atherosclerotic effects via suppressing migration and proliferation in the vessel wall. Nevertheless, larger-scale and long-term studies are still needed.

An efficient concordant integrative analysis of multiple large-scale two-sample expression data sets.

MOTIVATION: We have proposed a mixture model based approach to the concordant integrative analysis of multiple large-scale two-sample expression datasets. Since the mixture model is based on the transformed differential expression test P-values (z-scores), it is generally applicable to the expression data generated by either microarray or RNA-seq platforms. The mixture model is simple with three normal distribution components for each dataset to represent down-regulation, up-regulation and no differential expression. However, when the number of datasets increases, the model parameter space increases exponentially due to the component combination from different datasets. RESULTS: In this study, motivated by the well-known generalized estimating equations (GEEs) for longitudinal data analysis, we focus on the concordant components and assume that the proportions of non-concordant components follow a special structure. We discuss the exchangeable, multiset coefficient and autoregressive structures for model reduction, and their related expectation-maximization (EM) algorithms. Then, the parameter space is linear with the number of datasets. In our previous study, we have applied the general mixture model to three microarray datasets for lung cancer studies. We show that more gene sets (or pathways) can be detected by the reduced mixture model with the exchangeable structure. Furthermore, we show that more genes can also be detected by the reduced model. The Cancer Genome Atlas (TCGA) data have been increasingly collected. The advantage of incorporating the concordance feature has also been clearly demonstrated based on TCGA RNA sequencing data for studying two closely related types of cancer. AVAILABILITY AND IMPLEMENTATION: Additional results are included in a supplemental file. Computer program R-functions are freely available at http://home.gwu.edu/∼ylai/research/Concordance. CONTACT: ylai@gwu.edu. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Detecting discordance enrichment among a series of two-sample genome-wide expression data sets.

BACKGROUND: With the current microarray and RNA-seq technologies, two-sample genome-wide expression data have been widely collected in biological and medical studies. The related differential expression analysis and gene set enrichment analysis have been frequently conducted. Integrative analysis can be conducted when multiple data sets are available. In practice, discordant molecular behaviors among a series of data sets can be of biological and clinical interest. METHODS: In this study, a statistical method is proposed for detecting discordance gene set enrichment. Our method is based on a two-level multivariate normal mixture model. It is statistically efficient with linearly increased parameter space when the number of data sets is increased. The model-based probability of discordance enrichment can be calculated for gene set detection. RESULTS: We apply our method to a microarray expression data set collected from forty-five matched tumor/non-tumor pairs of tissues for studying pancreatic cancer. We divided the data set into a series of non-overlapping subsets according to the tumor/non-tumor paired expression ratio of gene PNLIP (pancreatic lipase, recently shown it association with pancreatic cancer). The log-ratio ranges from a negative value (e.g. more expressed in non-tumor tissue) to a positive value (e.g. more expressed in tumor tissue). Our purpose is to understand whether any gene sets are enriched in discordant behaviors among these subsets (when the log-ratio is increased from negative to positive). We focus on KEGG pathways. The detected pathways will be useful for our further understanding of the role of gene PNLIP in pancreatic cancer research. Among the top list of detected pathways, the neuroactive ligand receptor interaction and olfactory transduction pathways are the most significant two. Then, we consider gene TP53 that is well-known for its role as tumor suppressor in cancer research. The log-ratio also ranges from a negative value (e.g. more expressed in non-tumor tissue) to a positive value (e.g. more expressed in tumor tissue). We divided the microarray data set again according to the expression ratio of gene TP53. After the discordance enrichment analysis, we observed overall similar results and the above two pathways are still the most significant detections. More interestingly, only these two pathways have been identified for their association with pancreatic cancer in a pathway analysis of genome-wide association study (GWAS) data. CONCLUSIONS: This study illustrates that some disease-related pathways can be enriched in discordant molecular behaviors when an important disease-related gene changes its expression. Our proposed statistical method is useful in the detection of these pathways. Furthermore, our method can also be applied to genome-wide expression data collected by the recent RNA-seq technology.

Acute appendicitis: transcript profiling of blood identifies promising biomarkers and potential underlying processes.

BACKGROUND: The diagnosis of acute appendicitis can be surprisingly difficult without computed tomography, which carries significant radiation exposure. Circulating blood cells may carry informative changes in their RNA expression profile that would signal internal infection or inflammation of the appendix. METHODS: Genome-wide expression profiling was applied to whole blood RNA of acute appendicitis patients versus patients with other abdominal disorders, in order to identify biomarkers of appendicitis. From a large cohort of emergency patients, a discovery set of patients with surgically confirmed appendicitis, or abdominal pain from other causes, was identified. RNA from whole blood was profiled by microarrays, and RNA levels were filtered by a combined fold-change (>2) and p value (<0.05). A separate set of patients, including patients with respiratory infections, was used to validate a partial least squares discriminant (PLSD) prediction model. RESULTS: Transcript profiling identified 37 differentially expressed genes (DEG) in appendicitis versus abdominal pain patients. The DEG list contained 3 major ontologies: infection-related, inflammation-related, and ribosomal processing. Appendicitis patients had lower level of neutrophil defensin mRNA (DEFA1,3), but higher levels of alkaline phosphatase (ALPL) and interleukin-8 receptor-ß (CXCR2/IL8RB), which was confirmed in a larger cohort of 60 patients using droplet digital PCR (ddPCR). CONCLUSIONS: Patients with acute appendicitis have detectable changes in the mRNA expression levels of factors related to neutrophil innate defense systems. The low defensin mRNA levels suggest that appendicitis patient's immune cells are not directly activated by pathogens, but are primed by diffusible factors in the microenvironment of the infection. The detected biomarkers are consistent with prior evidence that biofilm-forming bacteria in the appendix may be an important factor in appendicitis.

Metataxonomic and Metagenomic Approaches vs. Culture-Based Techniques for Clinical Pathology.

Diagnoses that are both timely and accurate are critically important for patients with life-threatening or drug resistant infections. Technological improvements in High-Throughput Sequencing (HTS) have led to its use in pathogen detection and its application in clinical diagnoses of infectious diseases. The present study compares two HTS methods, 16S rRNA marker gene sequencing (metataxonomics) and whole metagenomic shotgun sequencing (metagenomics), in their respective abilities to match the same diagnosis as traditional culture methods (culture inference) for patients with ventilator associated pneumonia (VAP). The metagenomic analysis was able to produce the same diagnosis as culture methods at the species-level for five of the six samples, while the metataxonomic analysis was only able to produce results with the same species-level identification as culture for two of the six samples. These results indicate that metagenomic analyses have the accuracy needed for a clinical diagnostic tool, but full integration in diagnostic protocols is contingent on technological improvements to decrease turnaround time and lower costs.

Functional annotation of the vlinc class of non-coding RNAs using systems biology approach.

Functionality of the non-coding transcripts encoded by the human genome is the coveted goal of the modern genomics research. While commonly relied on the classical methods of forward genetics, integration of different genomics datasets in a global Systems Biology fashion presents a more productive avenue of achieving this very complex aim. Here we report application of a Systems Biology-based approach to dissect functionality of a newly identified vast class of very long intergenic non-coding (vlinc) RNAs. Using highly quantitative FANTOM5 CAGE dataset, we show that these RNAs could be grouped into 1542 novel human genes based on analysis of insulators that we show here indeed function as genomic barrier elements. We show that vlinc RNAs genes likely function in cisto activate nearby genes. This effect while most pronounced in closely spaced vlinc RNA-gene pairs can be detected over relatively large genomic distances. Furthermore, we identified 101 vlinc RNA genes likely involved in early embryogenesis based on patterns of their expression and regulation. We also found another 109 such genes potentially involved in cellular functions also happening at early stages of development such as proliferation, migration and apoptosis. Overall, we show that Systems Biology-based methods have great promise for functional annotation of non-coding RNAs.