Metabolic regulation of quiescence in plants.

Quiescence is a crucial survival attribute in which cell division is repressed in a reversible manner. Although quiescence has long been viewed as an inactive state, recent studies have shown that it is an actively monitored process that is influenced by environmental stimuli. Here, we provide a perspective of the quiescent state and discuss how this process is tuned by energy, nutrient and oxygen status, and the pathways that sense and transmit these signals. We not only highlight the governance of canonical regulators and signalling mechanisms that respond to changes in nutrient and energy status, but also consider the central significance of mitochondrial functions and cues as key regulators of nuclear gene expression. Furthermore, we discuss how reactive oxygen species and the associated redox processes, which are intrinsically linked to energy carbohydrate metabolism, also play a key role in the orchestration of quiescence.

Redox regulation of meristem quiescence: Outside/ in.

Quiescence is an essential property of meristematic cells, which restrains the cell cycle while retaining the capacity to divide. This crucial process not only facilitates life-long tissue homeostasis and regenerative capacity but also provides protection against adverse environmental conditions, enabling cells to conserve the proliferative potency while minimising DNA damage. As a survival attribute, quiescence is inherently regulated by the products of aerobic life, in particular reactive oxygen species (ROS) and the redox (reduction/ oxidation) mechanisms that plant have evolved to channel these into pervasive signals. Adaptive responses allow quiescent cells to compensate for reduced oxygen tension (hypoxia) in a reversible manner, while the regulated production of the superoxide anion (.O2-) facilitates cell division and the maintenance of stem cells. Here we discuss the role of ROS and redox reactions in the control of the quiescent state in plant meristems, and how this process is integrated with cellular energy and hormone biochemistry. We consider the pathways that sense and transmit redox signals with a focus on the central significance of redox regulation in the mitochondria and nucleus, which is a major regulator of quiescence in meristems. We discuss recent studies that suggest ROS are a critical component of the feedback loops that control stem cell identity and fate and suggest that the ROS/hypoxia interface is an important "outside/ in" positional cue for plant cells, particularly in meristems.

Integration of reactive oxygen species and nutrient signalling to shape root system architecture.

Yield losses due to nutrient deficiency are estimated as the primary cause of the yield gap worldwide. Understanding how plant roots perceive external nutrient status and elaborate morphological adaptations in response to it is necessary to develop reliable strategies to increase crop yield. In the last decade, reactive oxygen species (ROS) were shown to be key players of the mechanisms underlying root responses to nutrient limitation. ROS contribute in multiple ways to shape the root system in response to nutritional cues, both as direct effectors acting on cell wall architecture and as second messengers in signalling pathways. Here, we review the mutual interconnections existing between perception and signalling of the most common forms of the major macronutrients (nitrogen, phosphorus and potassium), and ROS in shaping plant root system architecture. We discuss recent advances in dissecting the integration of these elements and their impact on morphological traits of the root system, highlighting the functional ductility of ROS and enzymes implied in ROS metabolism, such as class III peroxidases.

Oxygen and reactive oxygen species-dependent regulation of plant growth and development.

Oxygen and reactive oxygen species (ROS) have been co-opted during evolution into the regulation of plant growth, development, and differentiation. ROS and oxidative signals arising from metabolism or phytohormone-mediated processes control almost every aspect of plant development from seed and bud dormancy, liberation of meristematic cells from the quiescent state, root and shoot growth, and architecture, to flowering and seed production. Moreover, the phytochrome and phytohormone-dependent transmissions of ROS waves are central to the systemic whole plant signaling pathways that integrate root and shoot growth. The sensing of oxygen availability through the PROTEOLYSIS 6 (PRT6) N-degron pathway functions alongside ROS production and signaling but how these pathways interact in developing organs remains poorly understood. Considerable progress has been made in our understanding of the nature of hydrogen peroxide sensors and the role of thiol-dependent signaling networks in the transmission of ROS signals. Reduction/oxidation (redox) changes in the glutathione (GSH) pool, glutaredoxins (GRXs), and thioredoxins (TRXs) are important in the control of growth mediated by phytohormone pathways. Although, it is clear that the redox states of proteins involved in plant growth and development are controlled by the NAD(P)H thioredoxin reductase (NTR)/TRX and reduced GSH/GRX systems of the cytosol, chloroplasts, mitochondria, and nucleus, we have only scratched the surface of this multilayered control and how redox-regulated processes interact with other cell signaling systems.

The bud dormancy disconnect: latent buds of grapevine are dormant during summer despite a high metabolic rate.

Grapevine (Vitis vinifera L.) displays wide plasticity to climate; however, the physiology of dormancy along a seasonal continuum is poorly understood. Here we investigated the apparent disconnect between dormancy and the underlying respiratory physiology and transcriptome of grapevine buds, from bud set in summer to bud burst in spring. The establishment of dormancy in summer was pronounced and reproducible; however, this was coupled with little or no change in physiology, indicated by respiration, hydration, and tissue oxygen tension. The release of dormancy was biphasic; the depth of dormancy declined substantially by mid-autumn, while the subsequent decline towards spring was moderate. Observed changes in physiology failed to explain the first phase of dormancy decline, in particular. Transcriptome data contrasting development from summer through to spring also indicated that dormancy was poorly reflected by metabolic quiescence during summer and autumn. Gene Ontology and enrichment data revealed the prevailing influence of abscisic acid (ABA)-related gene expression during the transition from summer to autumn, and promoter motif analysis suggested that photoperiod may play an important role in regulating ABA functions during the establishment of dormancy. Transcriptomic data from later transitions reinforced the importance of oxidation and hypoxia as physiological cues to regulate the maintenance of quiescence and resumption of growth. Collectively these data reveal a novel disconnect between growth and metabolic quiescence in grapevine following bud set, which requires further experimentation to explain the phenology and dormancy relationships.

Stress effects on the reactive oxygen species-dependent regulation of plant growth and development.

Plant growth is mediated by cell proliferation and expansion. Both processes are controlled by a network of endogenous factors such as phytohormones, reactive oxygen species (ROS), sugars, and other signals, which influence gene expression and post-translational regulation of proteins. Stress resilience requires rapid and appropriate responses in plant growth and development as well as defence. Regulation of ROS accumulation in different cellular compartments influences growth responses to abiotic and biotic stresses. While ROS are essential for growth, they are also implicated in the stress-induced cessation of growth and, in some cases, programmed cell death. It is widely accepted that redox post-translational modifications of key proteins determine the growth changes and cell fate responses to stress, but the molecular pathways and factors involved remain poorly characterized. Here we discuss ROS as a signalling molecule, the mechanisms of ROS-dependent regulation that influence protein-protein interactions, protein function, and turnover, together with the relocation of key proteins to different intracellular compartments in a manner that can alter cell fate. Understanding how the redox interactome responds to stress-induced increases in ROS may provide a road map to tailoring the dynamic ROS interactions that determine growth and cell fate in order to enhance stress resilience.

Root system architecture, physiological and transcriptional traits of soybean (Glycine max L.) in response to water deficit: A review.

Drought stress is the main limiting factor for global soybean growth and production. Genetic improvement for water and nutrient uptake efficiency is critical to advance tolerance and enable more sustainable and resilient production, underpinning yield growth. The identification of quantitative traits and genes related to water and nutrient uptake will enhance our understanding of the mechanisms of drought tolerance in soybean. This review summarizes drought stress in the context of the physiological traits that enable effective acclimation, with a particular focus on roots. Genes controlling root system architecture play an important role in water and nutrient availability, and therefore important targets for breeding strategies to improve drought tolerance. This review highlights the candidate genes that have been identified as regulators of important root traits and responses to water stress. Progress in our understanding of the function of particular genes, including GmACX1, GmMS and GmPEPCK are discussed in the context of developing a system-based platform for genetic improvement of drought tolerance in soybean.

Integrated single-cell and bulk gene expression and ATAC-seq reveals heterogeneity and early changes in pathways associated with resistance to cetuximab in HNSCC-sensitive cell lines.

BACKGROUND: Identifying potential resistance mechanisms while tumour cells still respond to therapy is critical to delay acquired resistance. METHODS: We generated the first comprehensive multi-omics, bulk and single-cell data in sensitive head and neck squamous cell carcinoma (HNSCC) cells to identify immediate responses to cetuximab. Two pathways potentially associated with resistance were focus of the study: regulation of receptor tyrosine kinases by TFAP2A transcription factor, and epithelial-to-mesenchymal transition (EMT). RESULTS: Single-cell RNA-seq demonstrates heterogeneity, with cell-specific TFAP2A and VIM expression profiles in response to treatment and also with global changes to various signalling pathways. RNA-seq and ATAC-seq reveal global changes within 5 days of therapy, suggesting early onset of mechanisms of resistance; and corroborates cell line heterogeneity, with different TFAP2A targets or EMT markers affected by therapy. Lack of TFAP2A expression is associated with HNSCC decreased growth, with cetuximab and JQ1 increasing the inhibitory effect. Regarding the EMT process, short-term cetuximab therapy has the strongest effect on inhibiting migration. TFAP2A silencing does not affect cell migration, supporting an independent role for both mechanisms in resistance. CONCLUSION: Overall, we show that immediate adaptive transcriptional and epigenetic changes induced by cetuximab are heterogeneous and cell type dependent; and independent mechanisms of resistance arise while tumour cells are still sensitive to therapy.

Correction: Integrated single-cell and bulk gene expression and ATAC-seq reveals heterogeneity and early changes in pathways associated with resistance to cetuximab in HNSCC-sensitive cell lines.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

The initiation of bud burst in grapevine features dynamic regulation of the apoplastic pore size.

The physiological constraints on bud burst in woody perennials, including vascular development and oxygenation, remain unresolved. Both light and tissue oxygen status have emerged as important cues for vascular development in other systems; however, grapevine buds have only a facultative light requirement, and data on the tissue oxygen status have been confounded by the spatial variability within the bud. Here, we analysed apoplastic development at early stages of grapevine bud burst and combined molecular modelling with histochemical techniques to determine the pore size of cell walls in grapevine buds. The data demonstrate that quiescent grapevine buds were impermeable to apoplastic dyes (acid fuchsin and eosin Y) until after bud burst was established. The molecular exclusion size was calculated to be 2.1 nm, which would exclude most macromolecules except simple sugars and phytohormones until after bud burst. We used micro-computed tomography to demonstrate that tissue oxygen partial pressure data correlated well with structural heterogeneity of the bud and differences in tissue density, confirming that the primary bud complex becomes rapidly and preferentially oxygenated during bud burst. Taken together, our results reveal that the apoplastic porosity is highly regulated during the early stages of bud burst, suggesting a role for vascular development in the initial, rapid oxygenation of the primary bud complex.

Tankyrase inhibition promotes a stable human naïve pluripotent state with improved functionality.

The derivation and maintenance of human pluripotent stem cells (hPSCs) in stable naïve pluripotent states has a wide impact in human developmental biology. However, hPSCs are unstable in classical naïve mouse embryonic stem cell (ESC) WNT and MEK/ERK signal inhibition (2i) culture. We show that a broad repertoire of conventional hESC and transgene-independent human induced pluripotent stem cell (hiPSC) lines could be reverted to stable human preimplantation inner cell mass (ICM)-like naïve states with only WNT, MEK/ERK, and tankyrase inhibition (LIF-3i). LIF-3i-reverted hPSCs retained normal karyotypes and genomic imprints, and attained defining mouse ESC-like functional features, including high clonal self-renewal, independence from MEK/ERK signaling, dependence on JAK/STAT3 and BMP4 signaling, and naïve-specific transcriptional and epigenetic configurations. Tankyrase inhibition promoted a stable acquisition of a human preimplantation ICM-like ground state via modulation of WNT signaling, and was most efficacious in efficiently reprogrammed conventional hiPSCs. Importantly, naïve reversion of a broad repertoire of conventional hiPSCs reduced lineage-primed gene expression and significantly improved their multilineage differentiation capacities. Stable naïve hPSCs with reduced genetic variability and improved functional pluripotency will have great utility in regenerative medicine and human disease modeling.

Splice Expression Variation Analysis (SEVA) for inter-tumor heterogeneity of gene isoform usage in cancer.

Motivation: Current bioinformatics methods to detect changes in gene isoform usage in distinct phenotypes compare the relative expected isoform usage in phenotypes. These statistics model differences in isoform usage in normal tissues, which have stable regulation of gene splicing. Pathological conditions, such as cancer, can have broken regulation of splicing that increases the heterogeneity of the expression of splice variants. Inferring events with such differential heterogeneity in gene isoform usage requires new statistical approaches. Results: We introduce Splice Expression Variability Analysis (SEVA) to model increased heterogeneity of splice variant usage between conditions (e.g. tumor and normal samples). SEVA uses a rank-based multivariate statistic that compares the variability of junction expression profiles within one condition to the variability within another. Simulated data show that SEVA is unique in modeling heterogeneity of gene isoform usage, and benchmark SEVA's performance against EBSeq, DiffSplice and rMATS that model differential isoform usage instead of heterogeneity. We confirm the accuracy of SEVA in identifying known splice variants in head and neck cancer and perform cross-study validation of novel splice variants. A novel comparison of splice variant heterogeneity between subtypes of head and neck cancer demonstrated unanticipated similarity between the heterogeneity of gene isoform usage in HPV-positive and HPV-negative subtypes and anticipated increased heterogeneity among HPV-negative samples with mutations in genes that regulate the splice variant machinery. These results show that SEVA accurately models differential heterogeneity of gene isoform usage from RNA-seq data. Availability and implementation: SEVA is implemented in the R/Bioconductor package GSReg. Contact: bahman@jhu.edu or favorov@sensi.org or ejfertig@jhmi.edu. Supplementary information: Supplementary data are available at Bioinformatics online.

Modelling predicts that soybean is poised to dominate crop production across Africa.

The superior agronomic and human nutritional properties of grain legumes (pulses) make them an ideal foundation for future sustainable agriculture. Legume-based farming is particularly important in Africa, where small-scale agricultural systems dominate the food production landscape. Legumes provide an inexpensive source of protein and nutrients to African households as well as natural fertilization for the soil. Although the consumption of traditionally grown legumes has started to decline, the production of soybeans (Glycine max Merr.) is spreading fast, especially across southern Africa. Predictions of future land-use allocation and production show that the soybean is poised to dominate future production across Africa. Land use models project an expansion of harvest area, whereas crop models project possible yield increases. Moreover, a seed change in farming strategy is underway. This is being driven largely by the combined cash crop value of products such as oils and the high nutritional benefits of soybean as an animal feed. Intensification of soybean production has the potential to reduce the dependence of Africa on soybean imports. However, a successful "soybean bonanza" across Africa necessitates an intensive research, development, extension, and policy agenda to ensure that soybean genetic improvements and production technology meet future demands for sustainable production.

CancerInSilico: An R/Bioconductor package for combining mathematical and statistical modeling to simulate time course bulk and single cell gene expression data in cancer.

Bioinformatics techniques to analyze time course bulk and single cell omics data are advancing. The absence of a known ground truth of the dynamics of molecular changes challenges benchmarking their performance on real data. Realistic simulated time-course datasets are essential to assess the performance of time course bioinformatics algorithms. We develop an R/Bioconductor package, CancerInSilico, to simulate bulk and single cell transcriptional data from a known ground truth obtained from mathematical models of cellular systems. This package contains a general R infrastructure for running cell-based models and simulating gene expression data based on the model states. We show how to use this package to simulate a gene expression data set and consequently benchmark analysis methods on this data set with a known ground truth. The package is freely available via Bioconductor: http://bioconductor.org/packages/CancerInSilico/.

DNA methylation genome-wide analysis in remnant and primary gastric cancers.

BACKGROUND: Although primary (PGC) and remnant gastric cancers (RGC) both originate from the same gastrointestinal organ, they have very distinct clinicopathological behaviors. We hypothesized that there would be distinct differences in DNA methylation patterns that would occur during carcinogenesis of RGC and PGC, and that the differences in methylation patterns may help identify the primary factor contributing to chronic inflammation in patients with RGC. METHODS: We investigated the genome-wide DNA methylation patterns of PGC and RGC tissues from 48 patients using the Infinium HumanMethylation450 Beadchip assay. The results were validated by quantitative methylation-specific PCR (qMSP) in separate, independent cohorts. RESULTS: We found that in our training cohort of 48 patients, the most variable genes from the gastric cancer tissues identified by the Infinium HumanMethylation450 Beadchip clustered the resultant heatmap into high and low methylation groups. On multivariate analysis, PGCs contributed significantly to the high methylation group (p = 0.004, OR 12.33), which suggested that the promoter methylation status in PGC is higher than that in RGC. Supporting this conclusion was the finding that in a separate qMSP analysis in a test cohort, the EPB41L3 gene, chosen because of its high ß value on microarray analysis in the gastric cancer tissues, had significantly higher DNA promoter methylation in cancer tissues in the validation PGC tissues than in RGC. CONCLUSIONS: This study demonstrated that promoter methylation status in PGC is higher than in RGC. This result may reflect the effects of the absence of Helicobacter pylori on the reduced DNA methylation in the remnant stomach.

PatternMarkers & GWCoGAPS for novel data-driven biomarkers via whole transcriptome NMF.

SUMMARY: Non-negative Matrix Factorization (NMF) algorithms associate gene expression with biological processes (e.g. time-course dynamics or disease subtypes). Compared with univariate associations, the relative weights of NMF solutions can obscure biomarkers. Therefore, we developed a novel patternMarkers statistic to extract genes for biological validation and enhanced visualization of NMF results. Finding novel and unbiased gene markers with patternMarkers requires whole-genome data. Therefore, we also developed Genome-Wide CoGAPS Analysis in Parallel Sets (GWCoGAPS), the first robust whole genome Bayesian NMF using the sparse, MCMC algorithm, CoGAPS. Additionally, a manual version of the GWCoGAPS algorithm contains analytic and visualization tools including patternMatcher, a Shiny web application. The decomposition in the manual pipeline can be replaced with any NMF algorithm, for further generalization of the software. Using these tools, we find granular brain-region and cell-type specific signatures with corresponding biomarkers in GTEx data, illustrating GWCoGAPS and patternMarkers ascertainment of data-driven biomarkers from whole-genome data. AVAILABILITY AND IMPLEMENTATION: PatternMarkers & GWCoGAPS are in the CoGAPS Bioconductor package (3.5) under the GPL license. CONTACT: gsteinobrien@jhmi.edu or ccolantu@jhmi.edu or ejfertig@jhmi.edu. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Developmental control of hypoxia during bud burst in grapevine.

Dormant or quiescent buds of woody perennials are often dense and in the case of grapevine (Vitis vinifera L.) have a low tissue oxygen status. The precise timing of the decision to resume growth is difficult to predict, but once committed, the increase in tissue oxygen status is rapid and developmentally regulated. Here, we show that more than a third of the grapevine homologues of widely conserved hypoxia-responsive genes and nearly a fifth of all grapevine genes possessing a plant hypoxia-responsive promoter element were differentially regulated during bud burst, in apparent harmony with resumption of meristem identity and cell-cycle gene regulation. We then investigated the molecular and biochemical properties of the grapevine ERF-VII homologues, which in other species are oxygen labile and function in transcriptional regulation of hypoxia-responsive genes. Each of the 3 VvERF-VIIs were substrates for oxygen-dependent proteolysis in vitro, as a function of the N-terminal cysteine. Collectively, these data support an important developmental function of oxygen-dependent signalling in determining the timing and effective coordination bud burst in grapevine. In addition, novel regulators, including GASA-, TCP-, MYB3R-, PLT-, and WUS-like transcription factors, were identified as hallmarks of the orderly and functional resumption of growth following quiescence in buds.

Two clinical phenotypes in polycythemia vera.

BACKGROUND: Polycythemia vera is the ultimate phenotypic consequence of the V617F mutation in Janus kinase 2 (encoded by JAK2), but the extent to which this mutation influences the behavior of the involved CD34+ hematopoietic stem cells is unknown. METHODS: We analyzed gene expression in CD34+ peripheral-blood cells from 19 patients with polycythemia vera, using oligonucleotide microarray technology after correcting for potential confounding by sex, since the phenotypic features of the disease differ between men and women. RESULTS: Men with polycythemia vera had twice as many up-regulated or down-regulated genes as women with polycythemia vera, in a comparison of gene expression in the patients and in healthy persons of the same sex, but there were 102 genes with differential regulation that was concordant in men and women. When these genes were used for class discovery by means of unsupervised hierarchical clustering, the 19 patients could be divided into two groups that did not differ significantly with respect to age, neutrophil JAK2 V617F allele burden, white-cell count, platelet count, or clonal dominance. However, they did differ significantly with respect to disease duration; hemoglobin level; frequency of thromboembolic events, palpable splenomegaly, and splenectomy; chemotherapy exposure; leukemic transformation; and survival. The unsupervised clustering was confirmed by a supervised approach with the use of a top-scoring-pair classifier that segregated the 19 patients into the same two phenotypic groups with 100% accuracy. CONCLUSIONS: Removing sex as a potential confounder, we identified an accurate molecular method for classifying patients with polycythemia vera according to disease behavior, independently of their JAK2 V617F allele burden, and identified previously unrecognized molecular pathways in polycythemia vera outside the canonical JAK2 pathway that may be amenable to targeted therapy. (Funded by the Department of Defense and the National Institutes of Health.).

Cell cycle arrest in plants: what distinguishes quiescence, dormancy and differentiated G1?

Background: Quiescence is a fundamental feature of plant life, which enables plasticity, renewal and fidelity of the somatic cell line. Cellular quiescence is defined by arrest in a particular phase of the cell cycle, typically G1 or G2; however, the regulation of quiescence and proliferation can also be considered across wider scales in space and time. As such, quiescence is a defining feature of plant development and phenology, from meristematic stem cell progenitors to terminally differentiated cells, as well as dormant or suppressed seeds and buds. While the physiology of each of these states differs considerably, each is referred to as 'cell cycle arrest' or 'G1 arrest'. Scope: Here the physiology and molecular regulation of (1) meristematic quiescence, (2) dormancy and (3) terminal differentiation (cell cycle exit) are considered in order to determine whether and how the molecular decisions guiding these nuclear states are distinct. A brief overview of the canonical cell cycle regulators is provided, and the genetic and genomic, as well as physiological, evidence is considered regarding two primary questions: (1) Are the canonical cell cycle regulators superior or subordinate in the regulation of quiescence? (2) Are these three modes of quiescence governed by distinct molecular controls? Conclusion: Meristematic quiescence, dormancy and terminal differentiation are each predominantly characterized by G1 arrest but regulated distinctly, at a level largely superior to the canonical cell cycle. Meristematic quiescence is intrinsically linked to non-cell-autonomous regulation of meristem cell identity, and particularly through the influence of ubiquitin-dependent proteolysis, in partnership with reactive oxygen species, abscisic acid and auxin. The regulation of terminal differentiation shares analogous features with meristematic quiescence, albeit with specific activators and a greater role for cytokinin signalling. Dormancy meanwhile appears to be regulated at the level of chromatin accessibility, by Polycomb group-type histone modifications of particular dormancy genes.

Characterization of functionally active gene fusions in human papillomavirus related oropharyngeal squamous cell carcinoma.

The Cancer Genome Atlas (TCGA) sequencing analysis of head and neck squamous cell carcinoma (HNSCC) recently reported on gene fusions, however, few human papillomavirus (HPV) positive samples were included, and the functional relevance of identified fusions was not explored. We therefore performed an independent analysis of gene fusions in HPV-positive oropharyngeal SCC (OPSCC). RNA sequencing was performed on 47 HPV-positive OPSCC primary tumors and 25 normal mucosal samples from cancer unaffected controls on an Illumina TruSeq platform. MapSplice2 was used for alignment and identification of fusion candidates. Putative fusions with less than five spanning reads, detected in normal tissues, or that mapped to the same gene were filtered out. Selected fusions were validated by RT-PCR and Sanger sequencing. Within 47 HPV-positive OPSCC tumors, 282 gene fusions were identified. Most fusions (85.1%) occurred in a single tumor, and the remaining fusions recurred in 2-16 tumors. Gene fusions were associated with significant up regulation of 16 genes (including EGFR and ERBB4) and down regulation of four genes (PTPRT, ZNF750, DLG2, SLCO5A1). Expression of these genes followed similar patterns of up regulation and down regulation in tumors without these fusions compared to normal tissue. Five of six gene fusions selected for validation were confirmed through RT-PCR and sequencing. This integrative analysis provides a method of prioritizing functionally relevant gene fusions that may be expanded to other tumor types. These results demonstrate that gene fusions may be one mechanism by which functionally relevant genes are altered in HPV-positive OPSCC.