Key Pathways and Genes of Arabidopsis thaliana and Arabidopsis halleri Roots under Cadmium Stress Responses: Differences and Similarities.

Cadmium (Cd) is among the world's major health concerns, as it renders soils unsuitable and unsafe for food and feed production. Phytoremediation has the potential to remediate Cd-polluted soils, but efforts are still needed to develop a deep understanding of the processes underlying it. In this study, we performed a comprehensive analysis of the root response to Cd stress in A. thaliana, which can phytostabilize Cd, and in A. halleri, which is a Cd hyperaccumulator. Suitable RNA-seq data were analyzed by WGCNA to identify modules of co-expressed genes specifically associated with Cd presence. The results evidenced that the genes of the hyperaccumulator A. halleri mostly associated with the Cd presence are finely regulated (up- and downregulated) and related to a general response to chemical and other stimuli. Additionally, in the case of A. thaliana, which can phytostabilize metals, the genes upregulated during Cd stress are related to a general response to chemical and other stimuli, while downregulated genes are associated with functions which, affecting root growth and development, determine a deep modification of the organ both at the cellular and physiological levels. Furthermore, key genes of the Cd-associated modules were identified and confirmed by differentially expressed gene (DEG) detection and external knowledge. Together, key functions and genes shed light on differences and similarities among the strategies that the plants use to cope with Cd and may be considered as possible targets for future research.

Network-Based Analysis to Identify Hub Genes Involved in Spatial Root Response to Mechanical Constrains.

Previous studies report that the asymmetric response, observed along the main poplar woody bent root axis, was strongly related to both the type of mechanical forces (compression or tension) and the intensity of force displacement. Despite a large number of targets that have been proposed to trigger this asymmetry, an understanding of the comprehensive and synergistic effect of the antistress spatially related pathways is still lacking. Recent progress in the bioinformatics area has the potential to fill these gaps through the use of in silico studies, able to investigate biological functions and pathway overlaps, and to identify promising targets in plant responses. Presently, for the first time, a comprehensive network-based analysis of proteomic signatures was used to identify functions and pivotal genes involved in the coordinated signalling pathways and molecular activities that asymmetrically modulate the response of different bent poplar root sectors and sides. To accomplish this aim, 66 candidate proteins, differentially represented across the poplar bent root sides and sectors, were grouped according to their abundance profile patterns and mapped, together with their first neighbours, on a high-confidence set of interactions from STRING to compose specific cluster-related subnetworks (I-VI). Successively, all subnetworks were explored by a functional gene set enrichment analysis to identify enriched gene ontology terms. Subnetworks were then analysed to identify the genes that are strongly interconnected with other genes (hub gene) and, thus, those that have a pivotal role in the bent root asymmetric response. The analysis revealed novel information regarding the response coordination, communication, and potential signalling pathways asymmetrically activated along the main root axis, delegated mainly to Ca2+ (for new lateral root formation) and ROS (for gravitropic response and lignin accumulation) signatures. Furthermore, some of the data indicate that the concave side of the bent sector, where the mechanical forces are most intense, communicates to the other (neighbour and distant) sectors, inducing spatially related strategies to ensure water uptake and accompanying cell modification. This information could be critical for understanding how plants maintain and improve their structural integrity-whenever and wherever it is necessary-in natural mechanical stress conditions.

Root and shoot biology of Arabidopsis halleri dissected by WGCNA: an insight into the organ pivotal pathways and genes of an hyperaccumulator.

Arabidopsis halleri is a hyperaccumulating pseudo-metallophyte and an emerging model to explore molecular basis of metal tolerance and hyperaccumulation. In this regard, understanding of interacting genes can be a crucial aspect as these interactions regulate several biological functions at molecular level in response to multiple signals. In this current study, we applied a weighted gene co-expression network analysis (WGCNA) on root and shoot RNA-seq data of A. halleri to predict the related scale-free organ specific co-expression networks, for the first time. A total of 19,653 genes of root and 18,081 genes of shoot were grouped into 14 modules and subjected to GO and KEGG enrichment analysis. "Photosynthesis" and "photosynthesis-antenna proteins" were identified as the most enriched and common pathway to both root and shoot. Whereas "glucosinolate biosynthesis," "autophagy," and "SNARE interactions in vesicular transport" were specific to root, and "circadian rhythm" was found to be enriched only in shoot. Later, hub and bottleneck genes were identified in each module by using cytoHubba plugin based on Cytoscape and scoring the relevance of each gene to the topology of the network. The modules with the most significant differential expression pattern across control and treatment (Cd-Zn treatment) were selected and their hub and bottleneck genes were screened to validate their possible involvement in heavy metal stress. Moreover, we combined the analysis of co-expression modules together with protein-protein interactions (PPIs), confirming some genes as potential candidates in plant heavy metal stress and as biomarkers. The results from this analysis shed the light on the pivotal functions to the hyperaccumulative trait of A. halleri, giving perspective to new paths for future research on this species.

Meristematic Connectome: A Cellular Coordinator of Plant Responses to Environmental Signals?

Mechanical stress in tree roots induces the production of reaction wood (RW) and the formation of new branch roots, both functioning to avoid anchorage failure and limb damage. The vascular cambium (VC) is the factor responsible for the onset of these responses as shown by their occurrence when all primary tissues and the root tips are removed. The data presented confirm that the VC is able to evaluate both the direction and magnitude of the mechanical forces experienced before coordinating the most fitting responses along the root axis whenever and wherever these are necessary. The coordination of these responses requires intense crosstalk between meristematic cells of the VC which may be very distant from the place where the mechanical stress is first detected. Signaling could be facilitated through plasmodesmata between meristematic cells. The mechanism of RW production also seems to be well conserved in the stem and this fact suggests that the VC could behave as a single structure spread along the plant body axis as a means to control the relationship between the plant and its environment. The observation that there are numerous morphological and functional similarities between different meristems and that some important regulatory mechanisms of meristem activity, such as homeostasis, are common to several meristems, supports the hypothesis that not only the VC but all apical, primary and secondary meristems present in the plant body behave as a single interconnected structure. We propose to name this structure "meristematic connectome" given the possibility that the sequence of meristems from root apex to shoot apex could represent a pluricellular network that facilitates long-distance signaling in the plant body. The possibility that the "meristematic connectome" could act as a single structure active in adjusting the plant body to its surrounding environment throughout the life of a plant is now proposed.

Reaction Wood Anatomical Traits and Hormonal Profiles in Poplar Bent Stem and Root.

Reaction wood (RW) formation is an innate physiological response of woody plants to counteract mechanical constraints in nature, reinforce structure and redirect growth toward the vertical direction. Differences and/or similarities between stem and root response to mechanical constraints remain almost unknown especially in relation to phytohormones distribution and RW characteristics. Thus, Populus nigra stem and root subjected to static non-destructive mid-term bending treatment were analyzed. The distribution of tension and compression forces was firstly modeled along the main bent stem and root axis; then, anatomical features, chemical composition, and a complete auxin and cytokinin metabolite profiles of the stretched convex and compressed concave side of three different bent stem and root sectors were analyzed. The results showed that in bent stems RW was produced on the upper stretched convex side whereas in bent roots it was produced on the lower compressed concave side. Anatomical features and chemical analysis showed that bent stem RW was characterized by a low number of vessel, poor lignification, and high carbohydrate, and thus gelatinous layer in fiber cell wall. Conversely, in bent root, RW was characterized by high vessel number and area, without any significant variation in carbohydrate and lignin content. An antagonistic interaction of auxins and different cytokinin forms/conjugates seems to regulate critical aspects of RW formation/development in stem and root to facilitate upward/downward organ bending. The observed differences between the response stem and root to bending highlight how hormonal signaling is highly organ-dependent.

PANEV: an R package for a pathway-based network visualization.

BACKGROUND: During the last decade, with the aim to solve the challenge of post-genomic and transcriptomic data mining, a plethora of tools have been developed to create, edit and analyze metabolic pathways. In particular, when a complex phenomenon is considered, the creation of a network of multiple interconnected pathways of interest could be useful to investigate the underlying biology and ultimately identify functional candidate genes affecting the trait under investigation. RESULTS: PANEV (PAthway NEtwork Visualizer) is an R package set for gene/pathway-based network visualization. Based on information available on KEGG, it visualizes genes within a network of multiple levels (from 1 to n) of interconnected upstream and downstream pathways. The network graph visualization helps to interpret functional profiles of a cluster of genes. CONCLUSIONS: The suite has no species constraints and it is ready to analyze genomic or transcriptomic outcomes. Users need to supply the list of candidate genes, specify the target pathway(s) and the number of interconnected downstream and upstream pathways (levels) required for the investigation. The package is available at https://github.com/vpalombo/PANEV.

Formation of Annual Ring Eccentricity in Coarse Roots within the Root Cage of Pinus ponderosa Growing on Slopes.

: The coarse roots of Pinus ponderosa included in the cage are the ones most involved in tree stability. This study explored the variations in traits, such as volume, cross-sectional area, and radius length of cage roots, and used those data to develop a mathematical model to better understand the type of forces occurring for each shallow lateral root segment belonging to different quadrants of the three-dimensional (3D) root system architecture. The pattern and intensity of these forces were modelled along the root segment from the branching point to the cage edge. Data of root cage volume in the upper 30 cm of soil showed a higher value in the downslope and windward quadrant while, at a deeper soil depth (> 30 cm), we found higher values in both upslope and leeward quadrants. The analysis of radius length and the cross-sectional area of the shallow lateral roots revealed the presence of a considerable degree of eccentricity of the annual rings at the branching point and at the cage edge. This eccentricity is due to the formation of compression wood, and the eccentricity changes from the top portion at the branching point to the bottom portion at the cage edge, which we hypothesize may be a response to the variation in mechanical forces occurring in the various zones of the cage. This hypothesis is supported by a mathematical model that shows how the pattern and intensity of different types of mechanical forces are present within the various quadrants of the same root system from the taproot to the cage edge.

Phylo_dCor: distance correlation as a novel metric for phylogenetic profiling.

BACKGROUND: Elaboration of powerful methods to predict functional and/or physical protein-protein interactions from genome sequence is one of the main tasks in the post-genomic era. Phylogenetic profiling allows the prediction of protein-protein interactions at a whole genome level in both Prokaryotes and Eukaryotes. For this reason it is considered one of the most promising methods. RESULTS: Here, we propose an improvement of phylogenetic profiling that enables handling of large genomic datasets and infer global protein-protein interactions. This method uses the distance correlation as a new measure of phylogenetic profile similarity. We constructed robust reference sets and developed Phylo-dCor, a parallelized version of the algorithm for calculating the distance correlation that makes it applicable to large genomic data. Using Saccharomyces cerevisiae and Escherichia coli genome datasets, we showed that Phylo-dCor outperforms phylogenetic profiling methods previously described based on the mutual information and Pearson's correlation as measures of profile similarity. CONCLUSIONS: In this work, we constructed and assessed robust reference sets and propose the distance correlation as a measure for comparing phylogenetic profiles. To make it applicable to large genomic data, we developed Phylo-dCor, a parallelized version of the algorithm for calculating the distance correlation. Two R scripts that can be run on a wide range of machines are available upon request.

An Integrated Approach to Explore Composition and Dynamics of Cholesterol-rich Membrane Microdomains in Sexual Stages of Malaria Parasite.

Membrane microdomains that include lipid rafts, are involved in key physiological and pathological processes and participate in the entry of endocellular pathogens. These assemblies, enriched in cholesterol and sphingolipids, form highly dynamic, liquid-ordered phases that can be separated from the bulk membranes thanks to their resistance to solubilization by nonionic detergents. To characterize complexity and dynamics of detergent-resistant membranes of sexual stages of the rodent malaria parasite Plasmodium berghei, here we propose an integrated study of raft components based on proteomics, lipid analysis and bioinformatics. This analysis revealed unexpected heterogeneity and unexplored pathways associated with these specialized assemblies. Protein-protein relationships and protein-lipid co-occurrence were described through multi-component networks. The proposed approach can be widely applied to virtually every cell type in different contexts and perturbations, under physiological and/or pathological conditions.

The FAD-dependent glycerol-3-phosphate dehydrogenase of Giardia duodenalis: an unconventional enzyme that interacts with the g14-3-3 and it is a target of the antitumoral compound NBDHEX.

The flagellated protozoan Giardia duodenalis is a worldwide parasite causing giardiasis, an acute and chronic diarrheal disease. Metabolism in G. duodenalis has a limited complexity thus making metabolic enzymes ideal targets for drug development. However, only few metabolic pathways (i.e., carbohydrates) have been described so far. Recently, the parasite homolog of the mitochondrial-like glycerol-3-phosphate dehydrogenase (gG3PD) has been identified among the interactors of the g14-3-3 protein. G3PD is involved in glycolysis, electron transport, glycerophospholipids metabolism, and hyperosmotic stress response, and is emerging as promising target in tumor treatment. In this work, we demonstrate that gG3PD is a functional flavoenzyme able to convert glycerol-3-phosphate into dihydroxyacetone phosphate and that its activity and the intracellular glycerol level increase during encystation. Taking advantage of co-immunoprecipitation assays and deletion mutants, we provide evidence that gG3PD and g14-3-3 interact at the trophozoite stage, the intracellular localization of gG3PD is stage dependent and it partially co-localizes with mitosomes during cyst development. Finally, we demonstrate that the gG3PD activity is affected by the antitumoral compound 6-(7-nitro-2,1,3-benzoxadiazol-4-ylthio)hexanol, that results more effective in vitro at killing G. duodenalis trophozoites than the reference drug metronidazole. Overall, our results highlight the involvement of gG3PD in processes crucial for the parasite survival thus proposing this enzyme as target for novel antigiardial interventions.