Gut flora metagenomic analysis coupled with metabolic and deep immune profiling in chronic kidney disease.

BACKGROUND: Perturbation of gut microbiota has been linked to chronic kidney disease (CKD), which was correlated with a sophisticated milieu of metabolic and immune dysregulation. METHODS: To clarify the underlying host-microbe interaction in CKD, we performed multi-omics measurements, including systems-level gut microbiome, targeted serum metabolome and deep immunotyping, in a cohort of patients and non-CKD controls. RESULTS: Our analyses on functional profiles of the gut microbiome showed a decrease in the diversity and abundance of carbohydrate-active enzyme (CAZyme) genes but an increase in the abundance of antibiotic resistance, nitrogen cycling enzyme and virulence factor genes in CKD. Moreover, models generated using measurements of serum metabolites (amino acids, bile acids and short-chain fatty acids) or immunotypes were predictive of renal impairment but less so than many of the functional profiles derived from gut microbiota, with the CAZyme genes being the top-performing model to accurately predict the early stage of diseases. In addition, co-occurrence analyses revealed coordinated host-microbe relationships in CKD. Specifically, the highest fractions of significant correlations were identified with circulating metabolites by several taxonomic and functional profiles of gut microbiome, while immunotype features were moderately associated with the abundance of microbiome-encoded metabolic pathways and serum levels of amino acids (e.g. B cell cluster tryptophan and B cell cluster tryptophan metabolism). CONCLUSION: Overall, our multi-omics integration revealed several signatures of systems-level gut microbiome in robust associations with host-microbe co-metabolites and renal function, which may have aetiological and diagnostic implications in CKD.

O-GlcNAcylation of master growth repressor DELLA by SECRET AGENT modulates multiple signaling pathways in Arabidopsis.

The DELLA family of transcription regulators functions as master growth repressors in plants by inhibiting phytohormone gibberellin (GA) signaling in response to developmental and environmental cues. DELLAs also play a central role in mediating cross-talk between GA and other signaling pathways via antagonistic direct interactions with key transcription factors. However, how these crucial protein-protein interactions can be dynamically regulated during plant development remains unclear. Here, we show that DELLAs are modified by the O-linked N-acetylglucosamine (O-GlcNAc) transferase (OGT) SECRET AGENT (SEC) in Arabidopsis. O-GlcNAcylation of the DELLA protein REPRESSOR OF ga1-3 (RGA) inhibits RGA binding to four of its interactors-PHYTOCHROME-INTERACTING FACTOR3 (PIF3), PIF4, JASMONATE-ZIM DOMAIN1, and BRASSINAZOLE-RESISTANT1 (BZR1)-that are key regulators in light, jasmonate, and brassinosteroid signaling pathways, respectively. Consistent with this, the sec-null mutant displayed reduced responses to GA and brassinosteroid and showed decreased expression of several common target genes of DELLAs, BZR1, and PIFs. Our results reveal a direct role of OGT in repressing DELLA activity and indicate that O-GlcNAcylation of DELLAs provides a fine-tuning mechanism in coordinating multiple signaling activities during plant development.

Phenotype-specific signatures of systems-level gut microbiome associated with childhood airway allergies.

BACKGROUND: Perturbation of gut symbiosis has been linked to childhood allergic diseases. However, the underlying host-microbe interaction connected with specific phenotypes is poorly understood. METHODS: To address this, integrative analyses of stool metagenomic and metabolomic profiles associated with IgE reactions in 56 children with mite-sensitized airway allergies (25 with rhinitis and 31 with asthma) and 28 nonallergic healthy controls were conducted. RESULTS: We noted a decrease in the number and abundance of gut microbiome-encoded carbohydrate-active enzyme (CAZyme) genes, accompanied with a reduction in species richness, in the asthmatic gut microflora but not in that from allergic rhinitis. Such loss of CAZymes was consistent with the observation that a CAZyme-linked decrease in fecal butyrate was found in asthmatics and negatively correlated with mite-specific IgE responses. Different from the CAZymes, we demonstrated an increase in α diversity at the virulome levels in asthmatic gut microbiota and identified phenotype-specific variations of gut virulome. Moreover, use of fecal metagenomic and metabolomic signatures resulted in distinct effects on differentiating rhinitis and asthma from nonallergic healthy controls. CONCLUSION: Overall, our integrative analyses reveal several signatures of systems-level gut microbiome in robust associations with fecal metabolites and disease phenotypes, which may be of etiological and diagnostic implications in childhood airway allergies.

MetaSMC: a coalescent-based shotgun sequence simulator for evolving microbial populations.

MOTIVATION: High-throughput sequencing technology has revolutionized the study of metagenomics and cancer evolution. In a relatively simple environment, a metagenomics sequencing data is dominated by a few species. By analyzing the alignment of reads from microbial species, single nucleotide polymorphisms can be discovered and the evolutionary history of the populations can be reconstructed. The ever-increasing read length will allow more detailed analysis about the evolutionary history of microbial or tumor cell population. A simulator of shotgun sequences from such populations will be helpful in the development or evaluation of analysis algorithms. RESULTS: Here, we described an efficient algorithm, MetaSMC, which simulates reads from evolving microbial populations. Based on the coalescent theory, our simulator supports all evolutionary scenarios supported by other coalescent simulators. In addition, the simulator supports various substitution models, including Jukes-Cantor, HKY85 and generalized time-reversible models. The simulator also supports mutator phenotypes by allowing different mutation rates and substitution models in different subpopulations. Our algorithm ignores unnecessary chromosomal segments and thus is more efficient than standard coalescent when recombination is frequent. We showed that the process behind our algorithm is equivalent to Sequentially Markov Coalescent with an incomplete sample. The accuracy of our algorithm was evaluated by summary statistics and likelihood curves derived from Monte Carlo integration over large number of random genealogies. AVAILABILITY AND IMPLEMENTATION: MetaSMC is written in C. The source code is available at https://github.com/tarjxvf/metasmc. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Integration of metagenomics-metabolomics reveals specific signatures and functions of airway microbiota in mite-sensitized childhood asthma.

BACKGROUND: Childhood asthma is a multifactorial inflammatory condition of the airways, associated with specific changes in respiratory microbiome and circulating metabolome. METHODS: To explore the functional capacity of asthmatic microbiome and its intricate connection with the host, we performed shotgun sequencing of airway microbiome and untargeted metabolomics profiling of serum samples in a cohort of children with mite-sensitized asthma and non-asthmatic controls. RESULTS: We observed higher gene counts and sample-to-sample dissimilarities in asthmatic microbiomes, indicating a more heterogeneous community structure and functionality among the cases than in controls. Moreover, we identified airway microbial species linked to changes in circulating metabolites and IgE responses of the host, including a positive correlation between Prevotella sp oral taxon 306 and dimethylglycine that were both decreased in patients. Several control-enriched species (Eubacterium sulci, Prevotella pallens, and Prevotella sp oral taxon 306) were inversely correlated with total and allergen-specific IgE levels. Genes related to microbial carbohydrate, amino acid, and lipid metabolism were differentially enriched, suggesting that changes in microbial metabolism may contribute to respiratory health in asthmatics. Pathway modules relevant to allergic responses were differentially abundant in asthmatic microbiome, such as enrichments for biofilm formation by Pseudomonas aeruginosa, membrane trafficking, histidine metabolism, and glycosaminoglycan degradation, and depletions for polycyclic aromatic hydrocarbon degradation. Further, we identified metagenomic and metabolomic markers (eg, Eubacterium sulci) to discriminate cases from the non-asthmatic controls. CONCLUSIONS: Our dual-omics data reveal the connections between respiratory microbes and circulating metabolites perturbed in mite-sensitized pediatric asthma, which may be of etiological and diagnostic implications.

The Arabidopsis O-fucosyltransferase SPINDLY activates nuclear growth repressor DELLA.

Plant development requires coordination among complex signaling networks to enhance the plant's adaptation to changing environments. DELLAs, transcription regulators originally identified as repressors of phytohormone gibberellin signaling, play a central role in integrating multiple signaling activities via direct protein interactions with key transcription factors. Here, we found that DELLA is mono-O-fucosylated by the novel O-fucosyltransferase SPINDLY (SPY) in Arabidopsis thaliana. O-fucosylation activates DELLA by promoting its interaction with key regulators in brassinosteroid- and light-signaling pathways, including BRASSINAZOLE-RESISTANT1 (BZR1), PHYTOCHROME-INTERACTING-FACTOR3 (PIF3) and PIF4. Moreover, spy mutants displayed elevated responses to gibberellin and brassinosteroid, and increased expression of common target genes of DELLAs, BZR1 and PIFs. Our study revealed that SPY-dependent protein O-fucosylation plays a key role in regulating plant development. This finding may have broader importance because SPY orthologs are conserved in prokaryotes and eukaryotes, thus suggesting that intracellular O-fucosylation may regulate a wide range of biological processes in diverse organisms.

HUA ENHANCER1 is involved in posttranscriptional regulation of positive and negative regulators in Arabidopsis photomorphogenesis.

Light regulates growth and developmental processes in plants via global transcriptome adjustment, translational control, and multilayered posttranslational modification of proteins. The transcriptional activation and repression of light-responsive genes has been well documented; however, the impact of posttranscriptional regulation on conveying light signals has been less addressed. Here, we examined whether optimal photomorphogenesis in Arabidopsis thaliana requires the proper biogenesis of small regulatory RNAs that play pivotal roles in the posttranscriptional regulation of gene expression. Arabidopsis carrying a mutation in HUA ENHANCER1 (HEN1), required for stabilization of small regulatory RNAs, showed defects in multiple aspects of photomorphogenic and skotomorphogenic development. HEN1 negatively regulated Arabidopsis photomorphogenesis. Light-activated HEN1 expression depended on the photoreceptors phytochrome A (phyA), phyB, cryptochrome 1 (cry1), and cry2 and key transcriptional regulators ELONGATED HYPOCOTYL5 (HY5) and HY5-HOMOLOG. We also demonstrate the involvement of the small regulatory RNAs miR157d and miR319 in modulating the expression of a positive regulator, HY5, and negative regulators TEOSINTE BRANCHED1, CYCLOIDEA AND PCF family proteins, respectively, for optimal photomorphogenic development in Arabidopsis.

The CrdRS two-component system in Helicobacter pylori responds to nitrosative stress.

Helicobacter pylori inhabits the gastric mucosa where it senses and responds to various stresses via a two-component systems (TCSs) that enable its persistent colonization. The aim of this study was to investigate whether any of the three paired TCSs (ArsRS, FleRS and CrdRS) in H. pylori respond to nitrosative stress. The results showed that the expression of crdS was significantly increased upon exposure to nitric oxide (NO). crdS-knockout (ΔcrdS) and crdR/crdS-knockout (ΔcrdRS) H. pylori, but not arsS-knockout (ΔarsS) or fleS-knockout (ΔfleS) H. pylori, showed a significant loss of viability upon exposure to NO compared with wild-type strain. Knockin crdS (ΔcrdS-in) significantly restored viability in the presence of NO. Global transcriptional profiling analysis of wild-type and ΔcrdS H. pylori in the presence or absence of NO showed that 101 genes were differentially expressed, including copper resistance determinant A (crdA), transport, binding and envelope proteins. The CrdR binding motifs were investigated by competitive electrophoretic mobility shift assay, which revealed that the two AC-rich regions in the crdA promoter region are required for binding. These results demonstrate that CrdR-crdA interaction enables H. pylori to survive under nitrosative stress.

Causal inference of gene regulation with subnetwork assembly from genetical genomics data.

Deciphering the causal networks of gene interactions is critical for identifying disease pathways and disease-causing genes. We introduce a method to reconstruct causal networks based on exploring phenotype-specific modules in the human interactome and including the expression quantitative trait loci (eQTLs) that underlie the joint expression variation of each module. Closely associated eQTLs help anchor the orientation of the network. To overcome the inherent computational complexity of causal network reconstruction, we first deduce the local causality of individual subnetworks using the selected eQTLs and module transcripts. These subnetworks are then integrated to infer a global causal network using a random-field ranking method, which was motivated by animal sociology. We demonstrate how effectively the inferred causality restores the regulatory structure of the networks that mediate lymph node metastasis in oral cancer. Network rewiring clearly characterizes the dynamic regulatory systems of distinct disease states. This study is the first to associate an RXRB-causal network with increased risks of nodal metastasis, tumor relapse, distant metastases and poor survival for oral cancer. Thus, identifying crucial upstream drivers of a signal cascade can facilitate the discovery of potential biomarkers and effective therapeutic targets.

Quantitative trait transcripts for nicotine resistance in Drosophila melanogaster.

Although most genetic association studies are performed with the intention of detecting nucleotide polymorphisms that are correlated with a complex trait, transcript abundance should also be expected to associate with diseases or phenotypes. We performed a scan for such quantitative trait transcripts in adult female heads of the fruit fly (Drosophila melanogaster) that might explain variation for nicotine resistance. The strongest association was seen for abundance of ornithine aminotransferase transcripts, implicating detoxification and neurotransmitter biosynthesis as mediators of the quantitative response to the drug. Subsequently, genetic analysis and metabolite profiling confirmed a complex role for ornithine and GABA levels in modification of survival time upon chronic nicotine exposure. Differences between populations from North Carolina and California suggest that the resistance mechanism may be an evolved response to environmental exposure.

Arabidopsis bZIP16 transcription factor integrates light and hormone signaling pathways to regulate early seedling development.

Transcriptomic adjustment plays an important role in Arabidopsis thaliana seed germination and deetiolation in response to environmental light signals. The G-box cis-element is commonly present in promoters of genes that respond positively or negatively to the light signal. In pursuing additional transcriptional regulators that modulate light-mediated transcriptome changes, we identified bZIP16, a basic region/Leu zipper motif transcription factor, by G-box DNA affinity chromatography. We confirmed that bZIP16 has G-box-specific binding activity. Analysis of bzip16 mutants revealed that bZIP16 is a negative regulator in light-mediated inhibition of cell elongation but a positive regulator in light-regulated seed germination. Transcriptome analysis supported that bZIP16 is primarily a transcriptional repressor regulating light-, gibberellic acid (GA)-, and abscisic acid (ABA)-responsive genes. Chromatin immunoprecipitation analysis revealed that bZIP16 could directly target ABA-responsive genes and RGA-like2, a DELLA gene in the GA signaling pathway. bZIP16 could also indirectly repress the expression of phytochrome interacting factoR3-like5, which encodes a basic helix-loop-helix protein coordinating hormone responses during seed germination. By repressing the expression of these genes, bZIP16 functions to promote seed germination and hypocotyl elongation during the early stages of Arabidopsis seedling development.

The role of Mss11 in Candida albicans biofilm formation.

Candida albicans is an opportunistic human pathogen that can form a biofilm on biotic or inert surfaces such as epithelia and clinical devices. In this study, we examine the formation of C. albicans biofilm by establishing a key gene-centered network based on protein-protein interaction (PPI) and gene expression datasets. Starting from C. albicans Cph1 and Efg1, transcription factors associated with morphogenesis of biofilm formation, a network elucidates the complex cellular process and predicts potential unknown components related to biofilm formation. Subsequently, we analyzed the functions of Mss11 among these identified proteins to test the efficiency of the proposed computational approach. MSS11-deleted mutants were compared with a wild-type strain, indicating that the mutant is defective in forming a mature biofilm and partially attenuates the virulence of C. albicans in an infected mouse model. Finally, a DNA microarray analysis was conducted to identify the potential target genes of C. albicans Mss11. The findings of this study clarify complex gene or protein interaction during the biofilm formation process of C. albicans, supporting the application of a systems biology approach to study fungal pathogenesis.

The chromatin modification by SUMO-2/3 but not SUMO-1 prevents the epigenetic activation of key immune-related genes during Kaposi's sarcoma associated herpesvirus reactivation.

BACKGROUND: SUMOylation, as part of the epigenetic regulation of transcription, has been intensively studied in lower eukaryotes that contain only a single SUMO protein; however, the functions of SUMOylation during mammalian epigenetic transcriptional regulation are largely uncharacterized. Mammals express three major SUMO paralogues: SUMO-1, SUMO-2, and SUMO-3 (normally referred to as SUMO-1 and SUMO-2/3). Herpesviruses, including Kaposi's sarcoma associated herpesvirus (KSHV), seem to have evolved mechanisms that directly or indirectly modulate the SUMO machinery in order to evade host immune surveillance, thus advancing their survival. Interestingly, KSHV encodes a SUMO E3 ligase, K-bZIP, with specificity toward SUMO-2/3 and is an excellent model for investigating the global functional differences between SUMO paralogues. RESULTS: We investigated the effect of experimental herpesvirus reactivation in a KSHV infected B lymphoma cell line on genomic SUMO-1 and SUMO-2/3 binding profiles together with the potential role of chromatin SUMOylation in transcription regulation. This was carried out via high-throughput sequencing analysis. Interestingly, chromatin immunoprecipitation sequencing (ChIP-seq) experiments showed that KSHV reactivation is accompanied by a significant increase in SUMO-2/3 modification around promoter regions, but SUMO-1 enrichment was absent. Expression profiling revealed that the SUMO-2/3 targeted genes are primarily highly transcribed genes that show no expression changes during viral reactivation. Gene ontology analysis further showed that these genes are involved in cellular immune responses and cytokine signaling. High-throughput annotation of SUMO occupancy of transcription factor binding sites (TFBS) pinpointed the presence of three master regulators of immune responses, IRF-1, IRF-2, and IRF-7, as potential SUMO-2/3 targeted transcriptional factors after KSHV reactivation. CONCLUSION: Our study is the first to identify differential genome-wide SUMO modifications between SUMO paralogues during herpesvirus reactivation. Our findings indicate that SUMO-2/3 modification near protein-coding gene promoters occurs in order to maintain host immune-related gene unaltered during viral reactivation.

Rhb1 regulates the expression of secreted aspartic protease 2 through the TOR signaling pathway in Candida albicans.

Candida albicans is a major fungal pathogen in humans. In C. albicans, secreted aspartyl protease 2 (Sap2) is the most highly expressed secreted aspartic protease in vitro and is a virulence factor. Recent research links the small GTPase Rhb1 to C. albicans target of rapamycin (TOR) signaling in response to nitrogen availability. The results of this study show that Rhb1 is related to cell growth through the control of SAP2 expression when protein is the major nitrogen source. This process involves various components of the TOR signaling pathway, including Tor1 kinase and its downstream effectors. TOR signaling not only controls SAP2 transcription but also affects Sap2 protein levels, possibly through general amino acid control. DNA microarray analysis identifies other target genes downstream of Rhb1 in addition to SAP2. These findings provide new insight into nutrients, Rhb1-TOR signaling, and expression of C. albicans virulence factor.

Deep immune profiling of patients with renal impairment unveils distinct immunotypes associated with disease severity.

Background: Chronic kidney disease (CKD) is pathologically correlated with a sophisticated milieu of innate and adaptive immune dysregulation, but the underlying immunological disturbances remain poorly understood. Methods: To address this, we comprehensively interrogated cellular and soluble elements of the immune system by using high-dimensional flow cytometry to analyze peripheral blood mononuclear cells and performing cytokine/chemokine profiling of serum samples, respectively, in a cohort of 69 patients and 19 non-CKD controls. Results: Altered serum levels of several cytokines/chemokines were identified, among which concentrations of stem cell factor (SCF) were found to be elevated with the progression of CKD and inversely correlated with estimated glomerular filtration rate (eGFR). Deep immunophenotyping analyses reveal a global change in immune modulation associated with CKD severity. Specifically, a decrease in the subsets of CD56dim natural killer (NK) cells (KLRG-1+CD38+CD64+CD15+CD197+) and monocytes (KLRG-1+CD38+PD-1+) was detected in severe CKD compared with controls and mild CKD. In addition, comparisons between mild and severe CKD demonstrated a loss of a mature B cell population (PD-1+CD197+IgD+HLA-DR+) in the advanced stages of disease. Further, we identified immunophenotypic markers to discriminate mild CKD from the controls, among which the portion of CD38+ monocytes was of particular value in early diagnosis. Conclusions: Our data unveil severity-specific immunological signatures perturbed in CKD patients.

Histone demethylase JMJD2A regulates Kaposi's sarcoma-associated herpesvirus replication and is targeted by a viral transcriptional factor.

The switch between the latency and lytic cycles of Kaposi's sarcoma-associated herpesvirus (KSHV) is accompanied by specific alterations of histone codes. Recently, comprehensive analysis of histone modifications of KSHV showed the deposition of H3K27me3 across the KSHV genome with two specific regions occupied by the heterochromatin marker H3K9me3. Here, we show that knockdown of JMJD2A, an H3K9me3 demethylase, attenuates viral titers, whereas its overexpression increases KSHV reactivation. JMJD2A is localized in regions of latent viral chromosomes that are deficient in the H3K9me3 mark, indicating that JMJD2A may be responsible for the low level of this mark on viral chromatin. The presence of JMJD2A on the latent genome maintains H3K9 in unmethylated form and signals the readiness of specific sets of viral genes to be reactivated. The demethylase activity of JMJD2A is important for KSHV reactivation, because a demethylase-deficient mutant cannot restore the JMJD2A knockdown phenotype. Interestingly, we found that the KSHV encoded K-bZIP associated with JMJD2A, resulting in the inhibition of demethylase activity of JMJD2A both in vivo and in vitro. Inhibition of JMJD2A by K-bZIP is likely due to a physical interaction which blocks substrate accessibility. A consequence of such an inhibition is increasing global levels of H3K9me3 and gene silencing. Consistently, K-bZIP overexpression resulted in a repression of ∼80% of the ≥2-fold differentially regulated genes compared to results for the uninduced control cells. The consequences of K-bZIP targeting JMJD2A during viral replication will be discussed. To our knowledge, this is the first description of a viral product shown to be a potent inhibitor of a host cellular histone demethylase.

Functional analysis reveals the possible role of the C-terminal sequences and PI motif in the function of lily (Lilium longiflorum) PISTILLATA (PI) orthologues.

Two lily (Lilium longiflorum) PISTILLATA (PI) genes, Lily MADS Box Gene 8 and 9 (LMADS8/9), were characterized. LMADS9 lacked 29 C-terminal amino acids including the PI motif that was present in LMADS8. Both LMADS8/9 mRNAs were prevalent in the first and second whorl tepals during all stages of development and were expressed in the stamen only in young flower buds. LMADS8/9 could both form homodimers, but the ability of LMADS8 homodimers to bind to CArG1 was relatively stronger than that of LMADS9 homodimers. 35S:LMADS8 completely, and 35S:LMADS9 only partially, rescued the second whorl petal formation and partially converted the first whorl sepal into a petal-like structure in Arabidopsis pi-1 mutants. Ectopic expression of LMADS8-C (with deletion of the 29 amino acids of the C-terminal sequence) or LMADS8-PI (with only the PI motif deleted) only partially rescued petal formation in pi mutants, which was similar to what was observed in 35S:LMADS9/pi plants. In contrast, 35:LMADS9+L8C (with the addition of the 29 amino acids of the LMADS8 C-terminal sequence) or 35S:LMADS9+L8PI (with the addition of the LMADS8 PI motif) demonstrated an increased ability to rescue petal formation in pi mutants, which was similar to what was observed in 35S:LMADS8/pi plants. Furthermore, ectopic expression of LMADS8-M (with the MADS domain truncated) generated more severe dominant negative phenotypes than those seen in 35S:LMADS9-M flowers. These results revealed that the 29 amino acids including the PI motif in the C-terminal region of the lily PI orthologue are valuable for its function in regulating perianth organ formation.

Kernel density weighted loess normalization improves the performance of detection within asymmetrical data.

BACKGROUND: Normalization of gene expression data has been studied for many years and various strategies have been formulated to deal with various types of data. Most normalization algorithms rely on the assumption that the number of up-regulated genes and the number of down-regulated genes are roughly the same. However, the well-known Golden Spike experiment presents a unique situation in which differentially regulated genes are biased toward one direction, thereby challenging the conclusions of previous bench mark studies. RESULTS: This study proposes two novel approaches, KDL and KDQ, based on kernel density estimation to improve upon the basic idea of invariant set selection. The key concept is to provide various importance scores to data points on the MA plot according to their proximity to the cluster of the null genes under the assumption that null genes are more densely distributed than those that are differentially regulated. The comparison is demonstrated in the Golden Spike experiment as well as with simulation data using the ROC curves and compression rates. KDL and KDQ in combination with GCRMA provided the best performance among all approaches. CONCLUSIONS: This study determined that methods based on invariant sets are better able to resolve the problem of asymmetry. Normalization, either before or after expression summary for probesets, improves performance to a similar degree.

A novel method to identify cooperative functional modules: study of module coordination in the Saccharomyces cerevisiae cell cycle.

BACKGROUND: Identifying key components in biological processes and their associations is critical for deciphering cellular functions. Recently, numerous gene expression and molecular interaction experiments have been reported in Saccharomyces cerevisiae, and these have enabled systematic studies. Although a number of approaches have been used to predict gene functions and interactions, tools that analyze the essential coordination of functional components in cellular processes still need to be developed. RESULTS: In this work, we present a new approach to study the cooperation of functional modules (sets of functionally related genes) in a specific cellular process. A cooperative module pair is defined as two modules that significantly cooperate with certain functional genes in a cellular process. This method identifies cooperative module pairs that significantly influence a cellular process and the correlated genes and interactions that are essential to that process. Using the yeast cell cycle as an example, we identified 101 cooperative module associations among 82 modules, and importantly, we established a cell cycle-specific cooperative module network. Most of the identified module pairs cover cooperative pathways and components essential to the cell cycle. We found that 14, 36, 18, 15, and 20 cooperative module pairs significantly cooperate with genes regulated in early G1, late G1, S, G2, and M phase, respectively. Fifty-nine module pairs that correlate with Cdc28 and other essential regulators were also identified. These results are consistent with previous studies and demonstrate that our methodology is effective for studying cooperative mechanisms in the cell cycle. CONCLUSIONS: In this work, we propose a new approach to identifying condition-related cooperative interactions, and importantly, we establish a cell cycle-specific cooperation module network. These results provide a global view of the cell cycle and the method can be used to discover the dynamic coordination properties of functional components in other cellular processes.

Prognostic Significance of O-GlcNAc and PKM2 in Hormone Receptor-Positive and HER2-Nonenriched Breast Cancer.

Predictive metabolic biomarkers for the recurrent luminal breast cancer (BC) with hormone receptor (HR)-positive and human epidermal growth factor receptor type 2 (HER2)-negative are lacking. High levels of O-GlcNAcylation (O-GlcNAc) and pyruvate kinase isoenzyme M2 (PKM2) are associated with malignancy in BC; however, the association with the recurrence risk remains unclear. We first conduct survival analysis by using the METABRIC dataset to assess the correlation of PKM2 expression with BC clinical outcomes. Next, patients with HR+/HER2- luminal BC were recruited for PKM2/O-GlcNAc testing. Logistic regression and receiver operating characteristic curve analysis were performed to evaluate the 10-year DFS predicted outcome. Survival analysis of the METABRIC dataset revealed that high expression of PKM2 was significantly associated with worse overall survival in luminal BC. The high expression of O-GlcNAc or PKM2 was a significant independent marker for poor 10-year DFS using immunohistochemical analysis. The PKM2 or O-GlcNAc status was a significant predictor of DFS, with the combination of PKM2-O-GlcNAc status and T stage greatly enhancing the predictive outcome potential. In summary, O-GlcNAc, PKM2, and T stage serve as good prognostic discriminators in HR+/HER2- luminal BC.