Interpreting omics data with pathway enrichment analysis.

Pathway enrichment analysis is indispensable for interpreting omics datasets and generating hypotheses. However, the foundations of enrichment analysis remain elusive to many biologists. Here, we discuss best practices in interpreting different types of omics data using pathway enrichment analysis and highlight the importance of considering intrinsic features of various types of omics data. We further explain major components that influence the outcomes of a pathway enrichment analysis, including defining background sets and choosing reference annotation databases. To improve reproducibility, we describe how to standardize reporting methodological details in publications. This article aims to serve as a primer for biologists to leverage the wealth of omics resources and motivate bioinformatics tool developers to enhance the power of pathway enrichment analysis.

Plant Metabolic Network 15: A resource of genome-wide metabolism databases for 126 plants and algae.

To understand and engineer plant metabolism, we need a comprehensive and accurate annotation of all metabolic information across plant species. As a step towards this goal, we generated genome-scale metabolic pathway databases of 126 algal and plant genomes, ranging from model organisms to crops to medicinal plants (https://plantcyc.org). Of these, 104 have not been reported before. We systematically evaluated the quality of the databases, which revealed that our semi-automated validation pipeline dramatically improves the quality. We then compared the metabolic content across the 126 organisms using multiple correspondence analysis and found that Brassicaceae, Poaceae, and Chlorophyta appeared as metabolically distinct groups. To demonstrate the utility of this resource, we used recently published sorghum transcriptomics data to discover previously unreported trends of metabolism underlying drought tolerance. We also used single-cell transcriptomics data from the Arabidopsis root to infer cell type-specific metabolic pathways. This work shows the quality and quantity of our resource and demonstrates its wide-ranging utility in integrating metabolism with other areas of plant biology.

METACLUSTER-an R package for context-specific expression analysis of metabolic gene clusters.

SUMMARY: Plants and microbes produce numerous compounds to cope with their environments but the biosynthetic pathways for most of these compounds have yet to be elucidated. Some biosynthetic pathways are encoded by enzymes collocated in the chromosome. To facilitate a more comprehensive condition and tissue-specific expression analysis of metabolic gene clusters, we developed METACLUSTER, a probabilistic framework for characterizing metabolic gene clusters using context-specific gene expression information. AVAILABILITY AND IMPLEMENTATION: METACLUSTER is freely available at https://github.com/mbanf/METACLUSTER. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

SECONDARY WALL ASSOCIATED MYB1 is a positive regulator of secondary cell wall thickening in Brachypodium distachyon and is not found in the Brassicaceae.

Grass biomass is comprised chiefly of secondary walls that surround fiber and xylem cells. A regulatory network of interacting transcription factors in part regulates cell wall thickening. We identified Brachypodium distachyon SECONDARY WALL ASSOCIATED MYB1 (SWAM1) as a potential regulator of secondary cell wall biosynthesis based on gene expression, phylogeny, and transgenic plant phenotypes. SWAM1 interacts with cellulose and lignin gene promoters with preferential binding to AC-rich sequence motifs commonly found in the promoters of cell wall-related genes. SWAM1 overexpression (SWAM-OE) lines had greater above-ground biomass with only a slight change in flowering time while SWAM1 dominant repressor (SWAM1-DR) plants were severely dwarfed with a striking reduction in lignin of sclerenchyma fibers and stem epidermal cell length. Cellulose, hemicellulose, and lignin genes were significantly down-regulated in SWAM1-DR plants and up-regulated in SWAM1-OE plants. There was no reduction in bioconversion yield in SWAM1-OE lines; however, it was significantly increased for SWAM1-DR samples. Phylogenetic and syntenic analyses strongly suggest that the SWAM1 clade was present in the last common ancestor between eudicots and grasses, but is not in the Brassicaceae. Collectively, these data suggest that SWAM1 is a transcriptional activator of secondary cell wall thickening and biomass accumulation in B. distachyon.

Comparative genomic analysis of the R2R3 MYB secondary cell wall regulators of Arabidopsis, poplar, rice, maize, and switchgrass.

BACKGROUND: R2R3 MYB proteins constitute one of the largest plant transcription factor clades and regulate diverse plant-specific processes. Several R2R3 MYB proteins act as regulators of secondary cell wall (SCW) biosynthesis in Arabidopsis thaliana (At), a dicotyledenous plant. Relatively few studies have examined SCW R2R3 MYB function in grasses, which may have diverged from dicots in terms of SCW regulatory mechanisms, as they have in cell wall composition and patterning. Understanding cell wall regulation is especially important for improving lignocellulosic bioenergy crops, such as switchgrass. RESULTS: Here, we describe the results of applying phylogenic, OrthoMCL, and sequence identity analyses to classify the R2R3 MYB family proteins from the annotated proteomes of Arabidposis, poplar, rice, maize and the initial genome (v0.0) and translated transcriptome of switchgrass (Panicum virgatum). We find that the R2R3 MYB proteins of the five species fall into 48 subgroups, including three dicot-specific, six grass-specific, and two panicoid grass-expanded subgroups. We observe four classes of phylogenetic relationships within the subgroups of known SCW-regulating MYB proteins between Arabidopsis and rice, ranging from likely one-to-one orthology (for AtMYB26, AtMYB103, AtMYB69) to no homologs identifiable (for AtMYB75). Microarray data for putative switchgrass SCW MYBs indicate that many maintain similar expression patterns with the Arabidopsis SCW regulators. However, some of the switchgrass-expanded candidate SCW MYBs exhibit differences in gene expression patterns among paralogs, consistent with subfunctionalization. Furthermore, some switchgrass representatives of grass-expanded clades have gene expression patterns consistent with regulating SCW development. CONCLUSIONS: Our analysis suggests that no single comparative genomics tool is able to provide a complete picture of the R2R3 MYB protein family without leaving ambiguities, and establishing likely false-negative and -positive relationships, but that used together a relatively clear view emerges. Generally, we find that most R2R3 MYBs that regulate SCW in Arabidopsis are likely conserved in the grasses. This comparative analysis of the R2R3 MYB family will facilitate transfer of understanding of regulatory mechanisms among species and enable control of SCW biosynthesis in switchgrass toward improving its biomass quality.

Persistent postoperative hypotension caused by subclinical empty sella syndrome after a simple surgery: A case report.

BACKGROUND: Empty sella is an anatomical and radiological finding of the herniation of the subarachnoid space into the pituitary fossa leading to a flattened pituitary gland. Patients with empty sella may present with various symptoms, including headache due to intracranial hypertension and endocrine symptoms related to the specific pituitary hormones affected. Here, we report a female patient who developed persistent postoperative hypotension caused by subclinical empty sella syndrome after a simple surgery. CASE SUMMARY: A 47-year-old woman underwent vocal cord polypectomy under general anesthesia with endotracheal intubation. She denied any medical history, and her vital signs were normal before the surgery. Anesthesia and surgery were uneventful. However, she developed dizziness, headache and persistent hypotension in the ward. Thus, intravenous dopamine was started to maintain normal blood pressure, which improved her symptoms. However, she remained dependent on dopamine for over 24 h without any obvious anesthesia- and surgery-related complications. An endocrine etiology was then suspected, and further examination showed a high prolactin level, a low normal adrenocorticotropic hormone level and a low cortisol level. Magnetic resonance imaging of the brain revealed an empty sella. Therefore, she was diagnosed with empty sella syndrome and secondary adrenal insufficiency. Her symptoms disappeared one week later after daily glucocorticoid supplement. CONCLUSION: Endocrine etiologies such as pituitary and adrenal-related dysfunction should be considered in patients showing persistent postoperative hypotension when anesthesia- and surgery-related factors are excluded.

Omics-guided metabolic pathway discovery in plants: Resources, approaches, and opportunities.

Plants produce a vast array of metabolites, the biosynthetic routes of which remain largely undetermined. Genome-scale enzyme and pathway annotations and omics technologies have revolutionized research to decrypt plant metabolism and produced a growing list of functionally characterized metabolic genes and pathways. However, what is known is still a tiny fraction of the metabolic capacity harbored by plants. Here, we review plant enzyme and pathway annotation resources and cutting-edge omics approaches to guide discovery and characterization of plant metabolic pathways. We also discuss strategies for improving enzyme function prediction by integrating protein 3D structure information and single cell omics. This review aims to serve as a primer for plant biologists to leverage omics datasets to facilitate understanding and engineering plant metabolism.

A novel bivalent chromatin associates with rapid induction of camalexin biosynthesis genes in response to a pathogen signal in Arabidopsis.

Temporal dynamics of gene expression underpin responses to internal and environmental stimuli. In eukaryotes, regulation of gene induction includes changing chromatin states at target genes and recruiting the transcriptional machinery that includes transcription factors. As one of the most potent defense compounds in Arabidopsis thaliana, camalexin can be rapidly induced by bacterial and fungal infections. Though several transcription factors controlling camalexin biosynthesis genes have been characterized, how the rapid activation of genes in this pathway upon a pathogen signal is enabled remains unknown. By combining publicly available epigenomic data with in vivo chromatin modification mapping, we found that camalexin biosynthesis genes are marked with two epigenetic modifications with opposite effects on gene expression, trimethylation of lysine 27 of histone 3 (H3K27me3) (repression) and acetylation of lysine 18 of histone 3 (H3K18ac) (activation), to form a previously uncharacterized type of bivalent chromatin. Mutants with reduced H3K27me3 or H3K18ac suggested that both modifications were required to determine the timing of gene expression and metabolite accumulation at an early stage of the stress response. Our study indicates that the H3K27me3-H3K18ac bivalent chromatin, which we name as kairostat, plays an important role in controlling the timely induction of gene expression upon stress stimuli in plants.

In the fight against harmful fungi and bacteria, plants have an arsenal of chemicals at their disposal. For instance, species in the crucifer family ­ which includes mustard, cabbages and the model plant Arabidopsis thaliana ­ can defend themselves with camalexin, a compound produced soon after the plant receives signals from its attacker. What controls this precise timing, however, is still unclear. For the genes that rule the production of camalexin to be 'read', interpreted, and ultimately converted into proteins, their DNA sequences first need to be physically accessible to the cell. This availability is controlled, in part, by adding or removing chemical groups onto histones, the spool-like structures which DNA wraps around. These precisely controlled modifications ultimately help to activate or repress a gene. Sometimes, activating and inhibiting chemical groups can be present in the same location, creating what is known as a bivalent chromatin domain. Zhao et al. investigated whether histone modifications regulate when A. thaliana produces camalexin in response to an attack. A combination of bioinformatics and experimental approaches highlighted two chemical modifications (one repressive, the other activating) which were present on the same histone, forming a previously unknown bivalent chromatin domain. Mutant plants which did not carry these modifications could not produce camalexin at the right time. Further experiments showed that under normal conditions, both histone modifications were present. However, when the plant was under attack, the level of repressive and activating modifications respectively decreased and increased, leading to gene activation. Together, the results by Zhao et al. suggest that both histone modifications are required for camalexin genes to respond appropriately to signals from a harmful agent. A deeper understanding of this new mechanism could, in turn, allow scientists to engineer crops that are better at resisting disease.

Interdependent iron and phosphorus availability controls photosynthesis through retrograde signaling.

Iron deficiency hampers photosynthesis and is associated with chlorosis. We recently showed that iron deficiency-induced chlorosis depends on phosphorus availability. How plants integrate these cues to control chlorophyll accumulation is unknown. Here, we show that iron limitation downregulates photosynthesis genes in a phosphorus-dependent manner. Using transcriptomics and genome-wide association analysis, we identify two genes, PHT4;4 encoding a chloroplastic ascorbate transporter and bZIP58, encoding a nuclear transcription factor, which prevent the downregulation of photosynthesis genes leading to the stay-green phenotype under iron-phosphorus deficiency. Joint limitation of these nutrients induces ascorbate accumulation by activating expression of an ascorbate biosynthesis gene, VTC4, which requires bZIP58. Furthermore, we demonstrate that chloroplastic ascorbate transport prevents the downregulation of photosynthesis genes under iron-phosphorus combined deficiency through modulation of ROS homeostasis. Our study uncovers a ROS-mediated chloroplastic retrograde signaling pathway to adapt photosynthesis to nutrient availability.

First plant cell atlas workshop report.

The societal challenges posed by a growing human population and climate change necessitate technical advances in plant science. Plant research makes vital contributions to society by advancing technologies that improve agricultural food production, biological energy capture and conversion, and human health. However, the plant biology community lacks a comprehensive understanding of molecular machinery, including their locations within cells, distributions and variations among different cell types, and real-time dynamics. Fortunately, rapid advances in molecular methods, imaging, proteomics, and metabolomics made in the last decade afford unprecedented opportunities to develop a molecular-level map of plant cells with high temporal and spatial resolution. The Plant Cell Atlas (PCA) initiative aims to generate a resource that will provide fresh insight into poorly understood aspects of plant cell structure and organization and enable the discovery of new cellular compartments and features. The PCA will be a community resource (www.plantcellatlas.org/) that describes the state of various plant cell types and integrates high-resolution spatio-temporal information of nucleic acids, proteins, and metabolites within plant cells. This first PCA initiative workshop convened scientists passionate about developing a comprehensive PCA to brainstorm about the state of the field, recent advances, the development of tools, and the future directions of this initiative. The workshop featured invited talks to share initial data, along with broader ideas for the PCA. Additionally, breakout sessions were organized around topics including the conceptual goals of the PCA, technical challenges, and community wants and needs. These activities connected scientists with diverse expertise and sparked important discussions about how to leverage and extend leading-edge technologies and develop new techniques. A major outcome of the workshop was that the community wishes to redefine concepts of plant cell types and tissues quantitatively. A long-term goal is to delineate all molecules within the cell at high spatio-temporal resolution, obtain information about interacting molecular networks, and identify the contribution of these networks to development of the organism as a whole. As a first step, we wish to create comprehensive cellular and subcellular biomolecular maps of transcripts, proteins, and metabolites, track the dynamic interactions of these molecules intra- and intercellularly, discern complete states and transitions of specialized cell types, and integrate these disparate data points to generate testable models of cellular function. Ultimately, the PCA initiative will have a substantial positive impact by empowering a broad, diverse group of scientists to forge exciting paths in the field of plant science, facilitating connections with interested stakeholders beyond the scientific community, and enabling new agricultural technologies for a sustainable future.

Genome-enabled discovery of anthraquinone biosynthesis in Senna tora.

Senna tora is a widely used medicinal plant. Its health benefits have been attributed to the large quantity of anthraquinones, but how they are made in plants remains a mystery. To identify the genes responsible for plant anthraquinone biosynthesis, we reveal the genome sequence of S. tora at the chromosome level with 526 Mb (96%) assembled into 13 chromosomes. Comparison among related plant species shows that a chalcone synthase-like (CHS-L) gene family has lineage-specifically and rapidly expanded in S. tora. Combining genomics, transcriptomics, metabolomics, and biochemistry, we identify a CHS-L gene contributing to the biosynthesis of anthraquinones. The S. tora reference genome will accelerate the discovery of biologically active anthraquinone biosynthesis pathways in medicinal plants.

Rice Genome-Scale Network Integration Reveals Transcriptional Regulators of Grass Cell Wall Synthesis.

Grasses have evolved distinct cell wall composition and patterning relative to dicotyledonous plants. However, despite the importance of this plant family, transcriptional regulation of its cell wall biosynthesis is poorly understood. To identify grass cell wall-associated transcription factors, we constructed the Rice Combined mutual Ranked Network (RCRN). The RCRN covers >90% of annotated rice (Oryza sativa) genes, is high quality, and includes most grass-specific cell wall genes, such as mixed-linkage glucan synthases and hydroxycinnamoyl acyltransferases. Comparing the RCRN and an equivalent Arabidopsis network suggests that grass orthologs of most genetically verified eudicot cell wall regulators also control this process in grasses, but some transcription factors vary significantly in network connectivity between these divergent species. Reverse genetics, yeast-one-hybrid, and protoplast-based assays reveal that OsMYB61a activates a grass-specific acyltransferase promoter, which confirms network predictions and supports grass-specific cell wall synthesis genes being incorporated into conserved regulatory circuits. In addition, 10 of 15 tested transcription factors, including six novel Wall-Associated regulators (WAP1, WACH1, WAHL1, WADH1, OsMYB13a, and OsMYB13b), alter abundance of cell wall-related transcripts when transiently expressed. The results highlight the quality of the RCRN for examining rice biology, provide insight into the evolution of cell wall regulation, and identify network nodes and edges that are possible leads for improving cell wall composition.

Perioperative Single-Donor Platelet Apheresis and Red Blood Cell Transfusion Impact on 90-Day and Overall Survival in Living Donor Liver Transplantation.

BACKGROUND:: Although many previous studies have confirmed that perioperative blood transfusion is associated with poor outcomes after liver transplantation (LT), few studies described the influence of single-donor platelet apheresis transfusion in living donor LT (LDLT). This study aimed to assess the effect of blood products on outcomes for LDLT recipients, focusing on apheresis platelets. METHODS:: This retrospective study included 126 recipients who underwent their first adult-to-adult LDLT. Twenty-four variables including consumption of blood products of 126 LDLT recipients were assessed for their link to short-term outcomes and overall survival. Kaplan-Meier survival curve and the log-rank test were used for recipient survival analysis. A multivariate Cox proportional-hazard model and a propensity score analysis were applied to adjust confounders after potential risk factors were identified by a univariate Cox analysis. RESULTS: Patients who received apheresis platelet transfusion had a lower 90-day cumulative survival (78.9% vs. 94.2%, P = 0.009), but had no significant difference in overall survival in the Cox model, compared with those without apheresis platelet transfusion. Units of apheresis platelet transfusion (hazard ratio [HR] = 3.103, 95% confidence interval [CI]: 1.720-5.600, P < 0.001) and preoperative platelet count (HR = 0.170, 95% CI: 0.040-0.730, P = 0.017) impacted 90-day survival independently. Multivariate Cox regression analysis also found that units of red blood cell (RBC) transfusion (HR = 1.036, 95% CI: 1.006-1.067, P = 0.018), recipient's age (HR = 1.045, 95% CI: 1.005-1.086, P = 0.025), and ABO blood group comparison (HR = 2.990, 95% CI: 1.341-6.669, P = 0.007) were independent risk factors for overall survival after LDLT. CONCLUSIONS:: This study suggested that apheresis platelets were only associated with early mortality but had no impact on overall survival in LDLT. Units of RBC, recipient's age, and ABO group comparison were independent predictors of long-term outcomes.

Proteomics Coupled with Metabolite and Cell Wall Profiling Reveal Metabolic Processes of a Developing Rice Stem Internode.

Internodes of grass stems function in mechanical support, transport, and, in some species, are a major sink organ for carbon in the form of cell wall polymers. This study reports cell wall composition, proteomic, and metabolite analyses of the rice elongating internode. Cellulose, lignin, and xylose increase as a percentage of cell wall material along eight segments of the second rice internode (internode II) at booting stage, from the younger to the older internode segments, indicating active cell wall synthesis. Liquid-chromatography tandem mass spectrometry (LC-MS/MS) of trypsin-digested proteins from this internode at booting reveals 2,547 proteins with at least two unique peptides in two biological replicates. The dataset includes many glycosyltransferases, acyltransferases, glycosyl hydrolases, cell wall-localized proteins, and protein kinases that have or may have functions in cell wall biosynthesis or remodeling. Phospho-enrichment of internode II peptides identified 21 unique phosphopeptides belonging to 20 phosphoproteins including a leucine rich repeat-III family receptor like kinase. GO over-representation and KEGG pathway analyses highlight the abundances of proteins involved in biosynthetic processes, especially the synthesis of secondary metabolites such as phenylpropanoids and flavonoids. LC-MS/MS of hot methanol-extracted secondary metabolites from internode II at four stages (booting/elongation, early mature, mature, and post mature) indicates that internode secondary metabolites are distinct from those of roots and leaves, and differ across stem maturation. This work fills a void of in-depth proteomics and metabolomics data for grass stems, specifically for rice, and provides baseline knowledge for more detailed studies of cell wall synthesis and other biological processes characteristic of internode development, toward improving grass agronomic properties.

[Effects of enflurane on myocardial reperfusion injury during cardiac surgery with cardiopulmonary bypass].

OBJECTIVE: To study the effects of inhalation enflurane (Enf) before aortic clamping on myocardial reperfusion injury in cardiac surgery with cardiopulmonary bypass (CPB). METHODS: hirty patents undergoing selective cardiac valve replacement were randomly allocated to three groups. Group I and group II inhaled 1.0 MAC and 0.5 MAC Enf before clamping aorta, respectively. Group III was the control group interval administration with Fentanyl. RESULTS: Immediately upon aortic clamp release (T2), the value of CK-MB, MDA and SOD of all the groups was significantly increased, however,their concentration did not peak significantly until T3 and T4(10 and 30 min after clamp aorta release). The levels at 60 min (T5) and 24 hours (T6) aorta were lower than T4 but still higher than T(0). At T3 and T4, CK-MB levels in group I were significantly higher than those in II and III groups (P=0.0220, 0.0108 and 0.0202, 0.0295). At T6, the CK-MB level of group II was significantly higher than that of group III (P<0.0001). At T4 and T5, the MDA value of group I was higher than that of group II (P=0.0060 and 0.0364). Meanwhile, the SOD level in group I was also higher than that of group II and group III at the T4 point (P<0.0001 and 0.0084). There was a correlation between the CK-MB value and the aorta clamping time,correlation coefficient range being 0.55 - 0.81,(P<0.05). However, there was no correlation between the CK-MB and MDA, SOD. CONCLUSION: There is ischemia reperfusion injury during cardiac surgery CPB with the increase of OFR production and elevation of the antioxidant reserve. Inhalation of large dose of enflurane may result in increased myocardial ischemia reperfusion injury manifested by elevated levels of myocardial enzymes and OFR production.