Metabolic Dynamics and Prediction of Gestational Age and Time to Delivery in Pregnant Women.

Metabolism during pregnancy is a dynamic and precisely programmed process, the failure of which can bring devastating consequences to the mother and fetus. To define a high-resolution temporal profile of metabolites during healthy pregnancy, we analyzed the untargeted metabolome of 784 weekly blood samples from 30 pregnant women. Broad changes and a highly choreographed profile were revealed: 4,995 metabolic features (of 9,651 total), 460 annotated compounds (of 687 total), and 34 human metabolic pathways (of 48 total) were significantly changed during pregnancy. Using linear models, we built a metabolic clock with five metabolites that time gestational age in high accordance with ultrasound (R = 0.92). Furthermore, two to three metabolites can identify when labor occurs (time to delivery within two, four, and eight weeks, AUROC ≥ 0.85). Our study represents a weekly characterization of the human pregnancy metabolome, providing a high-resolution landscape for understanding pregnancy with potential clinical utilities.

Tryptophan Metabolites Target Transmembrane and Immunoglobulin Domain-Containing 1 Signaling to Augment Renal Tubular Injury.

Chronic kidney disease (CKD) is characterized by the accumulation of uremic toxins and renal tubular damage. Tryptophan-derived uremic toxins [indoxyl sulfate (IS) and kynurenine (Kyn)] are well-characterized tubulotoxins. Emerging evidence suggests that transmembrane and immunoglobulin domain-containing 1 (TMIGD1) protects tubular cells and promotes survival. However, the direct molecular mechanism(s) underlying how these two opposing pathways crosstalk remains unknown. We posited that IS and Kyn mediate tubular toxicity through TMIGD1 and the loss of TMIGD1 augments tubular injury. Results from the current study showed that IS and Kyn suppressed TMIGD1 transcription in tubular cells in a dose-dependent manner. The wild-type CCAAT enhancer-binding protein ß (C/EBPß) enhanced, whereas a dominant-negative C/EBPß suppressed, TMIGD1 promoter activity. IS down-regulated C/EBPß in primary human renal tubular cells. The adenine-induced CKD, unilateral ureteric obstruction, and deoxycorticosterone acetate salt unilateral nephrectomy models showed reduced TMIGD1 expression in the renal tubules, which correlated with C/EBPß expression. C/EBPß levels negatively correlated with the IS and Kyn levels. Inactivation of TMIGD1 in mice significantly lowered acetylated tubulin, decreased tubular cell proliferation, caused severe tubular damage, and worsened renal function. Thus, the current results demonstrate that TMIGD1 protects renal tubular cells from renal injury in different models of CKD and uncovers a novel mechanism of tubulotoxicity of tryptophan-based uremic toxins.

Parasitoid-host eavesdropping reveals temperature coupling of preferences to communication signals without genetic coupling.

Receivers of acoustic communication signals evaluate signal features to identify conspecifics. Changes in the ambient temperature can alter these features, rendering species recognition a challenge. To maintain effective communication, temperature coupling-changes in receiver signal preferences that parallel temperature-induced changes in signal parameters-occurs among genetically coupled signallers and receivers. Whether eavesdroppers of communication signals exhibit temperature coupling is unknown. Here, we investigate if the parasitoid fly Ormia ochracea, an eavesdropper of cricket calling songs, exhibits song pulse rate preferences that are temperature coupled. We use a high-speed treadmill system to record walking phonotaxis at three ambient temperatures (21, 25, and 30°C) in response to songs that varied in pulse rates (20 to 90 pulses per second). Total walking distance, peak steering velocity, angular heading, and the phonotaxis performance index varied with song pulse rates and ambient temperature. The peak of phonotaxis performance index preference functions became broader and shifted to higher pulse rate values at higher temperatures. Temperature-related changes in cricket songs between 21 and 30°C did not drastically affect the ability of flies to recognize cricket calling songs. These results confirm that temperature coupling can occur in eavesdroppers that are not genetically coupled with signallers.

Behind the mask(ing): how frogs cope with noise.

Albert Feng was a pioneer in the field of auditory neuroethology who used frogs to investigate the neural basis of spectral and temporal processing and directional hearing. Among his many contributions was connecting neural mechanisms for sound pattern recognition and localization to the problems of auditory masking that frogs encounter when communicating in noisy, real-world environments. Feng's neurophysiological studies of auditory processing foreshadowed and inspired subsequent behavioral investigations of auditory masking in frogs. For frogs, vocal communication frequently occurs in breeding choruses, where males form dense aggregations and produce loud species-specific advertisement calls to attract potential mates and repel competitive rivals. In this review, we aim to highlight how Feng's research advanced our understanding of how frogs cope with noise. We structure our narrative around three themes woven throughout Feng's research-spectral, temporal, and directional processing-to illustrate how frogs can mitigate problems of auditory masking by exploiting frequency separation between signals and noise, temporal fluctuations in noise amplitude, and spatial separation between signals and noise. We conclude by proposing future research that would build on Feng's considerable legacy to advance our understanding of hearing and sound communication in frogs and other vertebrates.

Transcriptional dysregulation in developing trigeminal sensory neurons in the LgDel mouse model of DiGeorge 22q11.2 deletion syndrome.

LgDel mice, which model the heterozygous deletion of genes at human chromosome 22q11.2 associated with DiGeorge/22q11.2 deletion syndrome (22q11DS), have cranial nerve and craniofacial dysfunction as well as disrupted suckling, feeding and swallowing, similar to key 22q11DS phenotypes. Divergent trigeminal nerve (CN V) differentiation and altered trigeminal ganglion (CNgV) cellular composition prefigure these disruptions in LgDel embryos. We therefore asked whether a distinct transcriptional state in a specific population of early differentiating LgDel cranial sensory neurons, those in CNgV, a major source of innervation for appropriate oropharyngeal function, underlies this departure from typical development. LgDel versus wild-type (WT) CNgV transcriptomes differ significantly at E10.5 just after the ganglion has coalesced. Some changes parallel altered proportions of cranial placode versus cranial neural crest-derived CNgV cells. Others are consistent with a shift in anterior-posterior patterning associated with divergent LgDel cranial nerve differentiation. The most robust quantitative distinction, however, is statistically verifiable increased variability of expression levels for most of the over 17 000 genes expressed in common in LgDel versus WT CNgV. Thus, quantitative expression changes of functionally relevant genes and increased stochastic variation across the entire CNgV transcriptome at the onset of CN V differentiation prefigure subsequent disruption of cranial nerve differentiation and oropharyngeal function in LgDel mice.

Immunoarray Measurements of Parathyroid Hormone-Related Peptides Combined with Other Biomarkers to Diagnose Aggressive Prostate Cancer.

Parathyroid hormone-related peptide (PTHrP) is related to bone metastasis and hypercalcemia in prostate and breast cancers and should be an excellent biomarker for aggressive forms of these cancers. Current clinical detection protocols for PTHrP are immunoradiometric assay and radioimmunoassay but are not sensitive enough to detect PTHrPs at early stages. We recently evaluated a prostate cancer biomarker panel, including serum monocyte differentiation antigen (CD-14), ETS-related gene protein, pigment epithelial-derived factor, and insulin-like growth factor-1, with promise for identifying aggressive prostate cancers. This panel predicted the need for patient biopsy better than PSA alone. In the present paper, we report an ultrasensitive microfluidic assay for PTHrPs and evaluate their diagnostic value and the value of including them with our prior biomarker panel to diagnose aggressive forms of prostate cancer. The immunoarray features screen-printed carbon sensor electrodes coated with 5 nm glutathione gold nanoparticles with capture antibodies attached. PTHrPs are bound to a secondary antibody attached to a polyhorseradish peroxidase label and delivered to the sensors to provide high sensitivity when activated by H2O2 and a mediator. We obtained an unprecedented 0.3 fg mL-1 limit of detection for PTHrP bioactive moieties PTHrP 1-173 and PTHrP 1-86. We also report the first study of PTHrPs in a large serum pool to identify aggressive malignancies. In assays of 130 human patient serum samples, PTHrP levels distinguished between aggressive and indolent prostate cancers with 83-91% clinical sensitivity and 78-96% specificity. Logistic regression identified the best predictive model as a combination of PTHrP 1-86 and vascular endothelial growth factor-D. PTHrP 1-173 alone also showed a high ability to differentiate aggressive and indolent cancers.

Indoleamine 2,3-dioxygenase-1, a Novel Therapeutic Target for Post-Vascular Injury Thrombosis in CKD.

BACKGROUND: CKD, characterized by retained uremic solutes, is a strong and independent risk factor for thrombosis after vascular procedures . Urem ic solutes such as indoxyl sulfate (IS) and kynurenine (Kyn) mediate prothrombotic effect through tissue factor (TF). IS and Kyn biogenesis depends on multiple enzymes, with therapeutic implications unexplored. We examined the role of indoleamine 2,3-dioxygenase-1 (IDO-1), a rate-limiting enzyme of kynurenine biogenesis, in CKD-associated thrombosis after vascular injury. METHODS: IDO-1 expression in mice and human vessels was examined. IDO-1-/- mice, IDO-1 inhibitors, an adenine-induced CKD, and carotid artery injury models were used. RESULTS: Both global IDO-1-/- CKD mice and IDO-1 inhibitor in wild-type CKD mice showed reduced blood Kyn levels, TF expression in their arteries, and thrombogenicity compared with respective controls. Several advanced IDO-1 inhibitors downregulated TF expression in primary human aortic vascular smooth muscle cells specifically in response to uremic serum. Further mechanistic probing of arteries from an IS-specific mouse model, and CKD mice, showed upregulation of IDO-1 protein, which was due to inhibition of its polyubiquitination and degradation by IS in vascular smooth muscle cells. In two cohorts of patients with advanced CKD, blood IDO-1 activity was significantly higher in sera of study participants who subsequently developed thrombosis after endovascular interventions or vascular surgery. CONCLUSION: Leveraging genetic and pharmacologic manipulation in experimental models and data from human studies implicate IS as an inducer of IDO-1 and a perpetuator of the thrombotic milieu and supports IDO-1 as an antithrombotic target in CKD.

Prostate Cancer Diagnosis in the Clinic Using an 8-Protein Biomarker Panel.

The inability to distinguish aggressive from indolent prostate cancer is a longstanding clinical problem. Prostate specific antigen (PSA) tests and digital rectal exams cannot differentiate these forms. Because only ∼10% of diagnosed prostate cancer cases are aggressive, existing practice often results in overtreatment including unnecessary surgeries that degrade patients' quality of life. Here, we describe a fast microfluidic immunoarray optimized to determine 8-proteins simultaneously in 5 µL of blood serum for prostate cancer diagnostics. Using polymeric horseradish peroxidase (poly-HRP, 400 HRPs) labels to provide large signal amplification and limits of detection in the sub-fg mL-1 range, a protocol was devised for the optimization of the fast, accurate assays of 100-fold diluted serum samples. Analysis of 130 prostate cancer patient serum samples revealed that some members of the protein panel can distinguish aggressive from indolent cancers. Logistic regression was used to identify a subset of the panel, combining biomarker proteins ETS-related gene protein (ERG), insulin-like growth factor-1 (IGF-1), pigment epithelial-derived factor (PEDF), and serum monocyte differentiation antigen (CD-14) to predict whether a given patient should be referred for biopsy, which gave a much better predictive accuracy than PSA alone. This represents the first prostate cancer blood test that can predict which patients will have a high biopsy Gleason score, a standard pathology score used to grade tumors.

Metabolites in a mouse cancer model enhance venous thrombogenicity through the aryl hydrocarbon receptor-tissue factor axis.

Patients with malignancy are at 4- to 7-fold higher risk of venous thromboembolism (VTE), a potentially fatal, yet preventable complication. Although general mechanisms of thrombosis are enhanced in these patients, malignancy-specific triggers and their therapeutic implication remain poorly understood. Here we examined a colon cancer-specific VTE model and probed a set of metabolites with prothrombotic propensity in the inferior vena cava (IVC) ligation model. Athymic mice injected with human colon adenocarcinoma cells exhibited significantly higher IVC clot weights, a biological readout of venous thrombogenicity, compared with the control mice. Targeted metabolomics analysis of plasma of mice revealed an increase in the blood levels of kynurenine and indoxyl sulfate (tryptophan metabolites) in xenograft-bearing mice, which correlated positively with the increase in the IVC clot size. These metabolites are ligands of aryl hydrocarbon receptor (AHR) signaling. Accordingly, plasma from the xenograft-bearing mice activated the AHR pathway and augmented tissue factor (TF) and plasminogen activator inhibitor 1 (PAI-1) levels in venous endothelial cells in an AHR-dependent manner. Consistent with these findings, the endothelium from the IVC of xenograft-bearing animals revealed nuclear AHR and upregulated TF and PAI-1 expression, telltale signs of an activated AHR-TF/PAI-1 axis. Importantly, pharmacological inhibition of AHR activity suppressed TF and PAI-1 expression in endothelial cells of the IVC and reduced clot weights in both kynurenine-injected and xenograft-bearing mice. Together, these data show dysregulated tryptophan metabolites in a mouse cancer model, and they reveal a novel link between these metabolites and the control of the AHR-TF/PAI-1 axis and VTE in cancer.

Temporal and tissue-specific activation of aryl hydrocarbon receptor in discrete mouse models of kidney disease.

Emerging evidence in animal models of chronic kidney disease (CKD) implicates Aryl Hydrocarbon Receptor (AHR) signaling as a mediator of uremic toxicity. However, details about its tissue-specific and time-dependent activation in response to various renal pathologies remain poorly defined. Here, a comprehensive analysis of AHR induction was conducted in response to discrete models of kidney diseases using a transgenic mouse line expressing the AHR responsive-promoter tethered to a ß-galactosidase reporter gene. Following validation using a canonical AHR ligand (a dioxin derivative), the transgenic mice were subjected to adenine-induced and ischemia/reperfusion-induced injury models representing CKD and acute kidney injury (AKI), respectively, in humans. Indoxyl sulfate was artificially increased in mice through the drinking water and by inhibiting its excretion into the urine. Adenine-fed mice showed a distinct and significant increase in ß-galactosidase in the proximal and distal renal tubules, cardiac myocytes, hepatocytes, and microvasculature in the cerebral cortex. The pattern of ß-galactosidase increase coincided with the changes in serum indoxyl sulfate levels. Machine-learning-based image quantification revealed positive correlations between indoxyl sulfate levels and ß-galactosidase expression in various tissues. This pattern of ß-galactosidase expression was recapitulated in the indoxyl sulfate-specific model. The ischemia/reperfusion injury model showed increase in ß-galactosidase in renal tubules that persisted despite reduction in serum indoxyl sulfate and blood urea nitrogen levels. Thus, our results demonstrate a relationship between AHR activation in various tissues of mice with CKD or AKI and the levels of indoxyl sulfate. This study demonstrates the use of a reporter gene mouse to probe tissue-specific manifestations of uremia in translationally relevant animal models and provide hypothesis-generating insights into the mechanism of uremic toxicity that warrant further investigation.

Lung-to-ear sound transmission does not improve directional hearing in green treefrogs (Hyla cinerea).

Amphibians are unique among extant vertebrates in having middle ear cavities that are internally coupled to each other and to the lungs. In frogs, the lung-to-ear sound transmission pathway can influence the tympanum's inherent directionality, but what role such effects might play in directional hearing remains unclear. In this study of the American green treefrog (Hyla cinerea), we tested the hypothesis that the lung-to-ear sound transmission pathway functions to improve directional hearing, particularly in the context of intraspecific sexual communication. Using laser vibrometry, we measured the tympanum's vibration amplitude in females in response to a frequency modulated sweep presented from 12 sound incidence angles in azimuth. Tympanum directionality was determined across three states of lung inflation (inflated, deflated, reinflated) both for a single tympanum in the form of the vibration amplitude difference (VAD) and for binaural comparisons in the form of the interaural vibration amplitude difference (IVAD). The state of lung inflation had negligible effects (typically less than 0.5 dB) on both VADs and IVADs at frequencies emphasized in the advertisement calls produced by conspecific males (834 and 2730 Hz). Directionality at the peak resonance frequency of the lungs (1558 Hz) was improved by ∼3 dB for a single tympanum when the lungs were inflated versus deflated, but IVADs were not impacted by the state of lung inflation. Based on these results, we reject the hypothesis that the lung-to-ear sound transmission pathway functions to improve directional hearing in frogs.

Accessible Telemedicine Diagnostics with ELISA in a 3D Printed Pipette Tip.

We report herein a novel pipet-based "ELISA in a tip" as a new versatile diagnostic tool featuring better sensitivity, shorter incubation time, accessibility, and low sample and reagent volumes compared to traditional ELISA. Capture and analysis of data by a cell phone facilitates electronic delivery of results to health care providers. Pipette tips were designed and 3D printed as adapters to fit most commercial 50-200 µL pipettes. Capture antibodies (Ab1) are immobilized on the inner walls of the pipet tip, which serves as the assay compartment where samples and reagents are moved in and out by pipetting. Signals are generated using colorimetric or chemiluminescent (CL) reagents and can be quantified using a cell phone, CCD camera, or plate reader. We utilized pipet-tip ELISA to detect four cancer biomarker proteins with detection limits similar to or lower than microplate ELISAs at 25% assay cost and time. Recoveries of these proteins from spiked human serum were 85-115% or better, depending slightly on detection mode. Using CCD camera quantification of CL with femto-luminol reagent gave limits of detection (LOD) as low as 0.5 pg/mL. Patient samples (13) were assayed for 3 biomarker proteins with results well correlated to conventional ELISA and an established microfluidic electrochemical immunoassay.

Role of Alternative Splicing in Prostate Cancer Aggressiveness and Drug Resistance in African Americans.

Alternative splicing, the process of removing introns and joining exons of pre-mRNA, is critical for growth, development, tissue homeostasis, and species diversity. Dysregulation of alternative splicing can initiate and drive disease. Aberrant alternative splicing has been shown to promote the "hallmarks of cancer" in both hematological and solid cancers. Of interest, recent work has focused on the role of alternative splicing in prostate cancer and prostate cancer health disparities. We will provide a review of prostate cancer health disparities involving the African American population, alternative RNA splicing, and alternative splicing in prostate cancer. Lastly, we will summarize our work on differential alternative splicing in prostate cancer disparities and its implications for disparate health outcomes and therapeutic targets.

Single-nucleotide polymorphisms of stemness genes predicted to regulate RNA splicing, microRNA and oncogenic signaling are associated with prostate cancer survival.

Prostate cancer (PCa) is a clinically and molecularly heterogeneous disease, with variation in outcomes only partially predicted by grade and stage. Additional tools to distinguish indolent from aggressive disease are needed. Phenotypic characteristics of stemness correlate with poor cancer prognosis. Given this correlation, we identified single-nucleotide polymorphisms (SNPs) of stemness-related genes and examined their associations with PCa survival. SNPs within stemness-related genes were analyzed for association with overall survival of PCa in the Prostate, Lung, Colorectal and Ovarian Cancer Screening Trial. Significant SNPs predicted to be functional were selected for linkage disequilibrium analysis and combined and stratified analyses. Identified SNPs were evaluated for association with gene expression. SNPs of CD44 (rs9666607), ABCC1 (rs35605 and rs212091) and GDF15 (rs1058587) were associated with PCa survival and predicted to be functional. A role for rs9666607 of CD44 and rs35605 of ABCC1 in RNA splicing regulation, rs212091 of ABCC1 in miRNA binding site activity and rs1058587 of GDF15 in causing an amino acid change was predicted. These SNPs represent potential novel prognostic markers for overall survival of PCa and support a contribution of the stemness pathway to PCa patient outcome.

Automated 3D-Printed Microfluidic Array for Rapid Nanomaterial-Enhanced Detection of Multiple Proteins.

We report here the fabrication and validation of a novel 3D-printed, automated immunoarray to detect multiple proteins with ultralow detection limits. This low cost, miniature immunoarray employs electrochemiluminescent (ECL) detection measured with a CCD camera and employs touch-screen control of a micropump to facilitate automated use. The miniaturized array features prefilled reservoirs to deliver sample and reagents to a paper-thin pyrolytic graphite microwell detection chip to complete sandwich immunoassays. The detection chip achieves high sensitivity by using single-wall carbon nanotube-antibody conjugates in the microwells and employing massively labeled antibody-decorated RuBPY-silica nanoparticles to generate ECL. The total cost of an array is $0.65, and an eight-protein assay can be done in duplicate for $0.14 per protein with limits of detection (LOD) as low as 78-110 fg mL-1 in diluted serum. The electronic control system costs $210 in components. Utility of the automated immunoarray was demonstrated by detecting an eight-protein prostate cancer biomarker panel in human serum samples in 25 min. The system is well suited to future clinical and point-of-care diagnostic testing and could be used in resource-limited environments.

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.

Hard to swallow: Developmental biological insights into pediatric dysphagia.

Pediatric dysphagia-feeding and swallowing difficulties that begin at birth, last throughout childhood, and continue into maturity--is one of the most common, least understood complications in children with developmental disorders. We argue that a major cause of pediatric dysphagia is altered hindbrain patterning during pre-natal development. Such changes can compromise craniofacial structures including oropharyngeal muscles and skeletal elements as well as motor and sensory circuits necessary for normal feeding and swallowing. Animal models of developmental disorders that include pediatric dysphagia in their phenotypic spectrum can provide mechanistic insight into pathogenesis of feeding and swallowing difficulties. A fairly common human genetic developmental disorder, DiGeorge/22q11.2 Deletion Syndrome (22q11DS) includes a substantial incidence of pediatric dysphagia in its phenotypic spectrum. Infant mice carrying a parallel deletion to 22q11DS patients have feeding and swallowing difficulties that approximate those seen in pediatric dysphagia. Altered hindbrain patterning, craniofacial malformations, and changes in cranial nerve growth prefigure these difficulties. Thus, in addition to craniofacial and pharyngeal anomalies that arise independently of altered neural development, pediatric dysphagia may result from disrupted hindbrain patterning and its impact on peripheral and central neural circuit development critical for feeding and swallowing. The mechanisms that disrupt hindbrain patterning and circuitry may provide a foundation to develop novel therapeutic approaches for improved clinical management of pediatric dysphagia.

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.