Network and biosignature analysis for the integration of transcriptomic and metabolomic data to characterize leaf senescence process in sunflower

Aging and obesity prime the methylome and transcriptome of adipose stem cells for disease and dysfunction

The epigenome of stem cells occupies a critical interface between genes and environment, serving to regulate expression through modification by intrinsic and extrinsic factors. We hypothesized that aging and obesity, which represent major risk factors for a variety of diseases, synergistically modify the epigenome of adult adipose stem cells (ASCs). Using integrated RNA- and targeted bisulfite-sequencing in murine ASCs from lean and obese mice at 5- and 12-months of age, we identified global DNA hypomethylation with either aging or obesity, and a synergistic effect of aging combined with obesity. The transcriptome of ASCs in lean mice was relatively stable to the effects of age, but this was not true in obese mice. Functional pathway analyses identified a subset of genes with critical roles in progenitors and in diseases of obesity and aging. Specifically, Mapt, Nr3c2, App, and Ctnnb1 emerged as potential hypomethylated upstream regulators in both aging and obesity (AL vs. YL and AO vs. YO), and App, Ctnnb1, Hipk2, Id2, and Tp53 exhibited additional effects of aging in obese animals. Furthermore, Foxo3 and Ccnd1 were potential hypermethylated upstream regulators of healthy aging (AL vs. YL), and of the effects of obesity in young animals (YO vs. YL), suggesting that these factors could play a role in accelerated aging with obesity. Finally, we identified candidate driver genes that appeared recurrently in all analyses and comparisons undertaken. Further mechanistic studies are needed to validate the roles of these genes capable of priming ASCs for dysfunction in aging- and obesity-associated pathologies.

The Eucalyptus grandis chloroplast proteome: Seasonal variations in leaf development

Chloroplast metabolism is very sensitive to environmental fluctuations and is intimately related to plant leaf development. Characterization of the chloroplast proteome dynamics can contribute to a better understanding on plant adaptation to different climate scenarios and leaf development processes. Herein, we carried out a discovery-driven analysis of the Eucalyptus grandis chloroplast proteome during leaf maturation and throughout different seasons of the year. The chloroplast proteome from young leaves differed the most from all assessed samples. Most upregulated proteins identified in mature and young leaves were those related to catabolic-redox signaling and biogenesis processes, respectively. Seasonal dynamics revealed unique proteome features in the fall and spring periods. The most abundant chloroplast protein in humid (wet) seasons (spring and summer) was a small subunit of RuBisCO, while in the dry periods (fall and winter) the proteins that showed the most pronounced accumulation were associated with photo-oxidative damage, Calvin cycle, shikimate pathway, and detoxification. Our investigation of the chloroplast proteome dynamics during leaf development revealed significant alterations in relation to the maturation event. Our findings also suggest that transition seasons induced the most pronounced chloroplast proteome changes over the year. This study contributes to a more comprehensive understanding on the subcellular mechanisms that lead to plant leaf adaptation and ultimately gives more insights into Eucalyptus grandis phenology.

Genome-wide and comparative phylogenetic analysis of senescence-associated NAC transcription factors in sunflower (Helianthus annuus)

BACKGROUND

Leaf senescence delay impacts positively in grain yield by maintaining the photosynthetic area during the reproductive stage and during grain filling. Therefore a comprehensive understanding of the gene families associated with leaf senescence is essential. NAC transcription factors (TF) form a large plant-specific gene family involved in regulating development, senescence, and responses to biotic and abiotic stresses. The main goal of this work was to identify sunflower NAC TF (HaNAC) and their association with senescence, studying their orthologous to understand possible functional relationships between genes of different species.

RESULTS

To clarify the orthologous relationships, we used an in-depth comparative study of four divergent taxa, in dicots and monocots, with completely sequenced genomes (Arabidopsis thaliana, Vitis vinifera, Musa acuminata and Oryza sativa). These orthologous groups provide a curated resource for large scale protein sequence annotation of NAC TF. From the 151 HaNAC genes detected in the latest version of the sunflower genome, 50 genes were associated with senescence traits. These genes showed significant differential expression in two contrasting lines according to an RNAseq assay. An assessment of overexpressing the Arabidopsis line for HaNAC001 (a gene of the same orthologous group of Arabidopsis thaliana ORE1) revealed that this line displayed a significantly higher number of senescent leaves and a pronounced change in development rate.

CONCLUSIONS

This finding suggests HaNAC001 as an interesting candidate to explore the molecular regulation of senescence in sunflower.

Stem girdling affects the onset of autumn senescence in aspen in interaction with metabolic signals

Autumn senescence in aspen (Populus tremula) is precisely timed every year to relocate nutrients from leaves to storage organs before winter. Here we demonstrate how stem girdling, which leads to the accumulation of photosynthates in the crown, influences senescence. Girdling resulted in an early onset of senescence, but the chlorophyll degradation was slower and nitrogen more efficiently resorbed than during normal autumn senescence. Girdled stems accumulated or retained anthocyanins potentially providing photoprotection in senescing leaves. Girdling of one stem in a clonal stand sharing the same root stock did not affect senescence in the others, showing that the stems were autonomous in this respect. One girdled stem with unusually high chlorophyll and nitrogen contents maintained low carbon-to-nitrogen (C/N) ratio and did not show early senescence or depleted chlorophyll level unlike the other girdled stems suggesting that the responses depended on the genotype or its carbon and nitrogen status. Metabolite analysis highlighted that the tricarboxylic acid (TCA) cycle, salicylic acid pathway, and redox homeostasis are involved in the regulation of girdling-induced senescence. We propose that disrupted sink-source relation and C/N status can provide cues through the TCA cycle and phytohormone signaling to override the phenological control of autumn senescence in the girdled stems.

Dietary switch to Western diet induces hypothalamic adaptation associated with gut microbiota dysbiosis in rats

BACKGROUND

Early hyperphagia and hypothalamic inflammation encountered after Western diet (WD) are linked to rodent propensity to obesity. Inflammation in several brain structures has been associated with gut dysbiosis. Since gut microbiota is highly sensitive to dietary changes, we hypothesised that immediate gut microbiota adaptation to WD in rats is involved in inflammation-related hypothalamic modifications.

METHODS

We evaluated short-term impact of WD consumption (2 h, 1, 2 and 4 days) on hypothalamic metabolome and caecal microbiota composition and metabolome. Data integration analyses were performed to uncover potential relationships among these three datasets. Finally, changes in hypothalamic gene expression in absence of gut microbiota were evaluated in germ-free rats fed WD for 2 days.

RESULTS

WD quickly and profoundly affected the levels of several hypothalamic metabolites, especially oxidative stress markers. In parallel, WD consumption reduced caecal microbiota diversity, modified its composition towards pro-inflammatory profile and changed caecal metabolome. Data integration identified strong correlations between gut microbiota sub-networks, unidentified caecal metabolites and hypothalamic oxidative stress metabolites. Germ-free rats displayed reduced energy intake and no changes in redox homoeostasis machinery expression or pro-inflammatory cytokines after 2 days of WD, in contrast to conventional rats, which exhibited increased SOD2, GLRX and IL-6 mRNA levels.

CONCLUSION

A potentially pro-inflammatory gut microbiota and an early hypothalamic oxidative stress appear shortly after WD introduction. Tripartite data integration highlighted putative links between gut microbiota sub-networks and hypothalamic oxidative stress. Together with the absence of hypothalamic modifications in germ-free rats, this strongly suggests the involvement of the microbiota-hypothalamus axis in rat adaptation to WD introduction and in energy homoeostasis regulation.

Gene Coexpression Network Reveals Insights into the Origin and Evolution of a Theanine-Associated Regulatory Module in Non-Camellia and Camellia Species

Theanine (thea) is one of the most important plant-derived characteristic secondary metabolites and a major healthcare product because of its beneficial biological activities, such as being an antianxiety agent, promoting memory, and lowering blood pressure. Thea mostly accumulates in Camellia plants and is especially rich in Camellia sinensis (tea plant). Although some functional genes (e.g., TS, GOGAT, and GS) attributed to thea accumulation have been separately well explored in tea plants, the evolution of a regulatory module (highly interacting gene group) related to thea metabolism remains to be elaborated. Herein, a thea-associated regulatory module (TARM) was mined by using a comprehensive analysis of a weighted gene coexpression network in Camellia and non-Camellia species. Comparative genomic analysis of 84 green plant species revealed that TARM originated from the ancestor of green plants (algae) and that TARM genes were recruited from different evolutionary nodes with the most gene duplication events at the early stage. Among the TARM genes, two core transcription factors of NAC080 and LBD38 were deduced, which may play a crucial role in regulating the biosynthesis of thea. Our findings provide the first insights into the origin and evolution of TARM and indicate a promising paradigm for identifying vital regulatory genes involved in thea metabolism.

Transcriptomic profiling of Solanum peruvianum LA3858 revealed a Mi-3-mediated hypersensitive response to Meloidogyne incognita

Background

The Mi-1 gene was the first identified and cloned gene that provides resistance to root-knot nematodes (RKNs) in cultivated tomato. However, owing to its temperature sensitivity, this gene does not meet the need for breeding disease-resistant plants that grow under high temperature. In this study, Mi-3 was isolated from the wild species PI 126443 (LA3858) and was shown to display heat-stable resistance to RKNs. However, the mechanism that regulates this resistance remains unknown.

Results

In this study, 4760, 1024 and 137 differentially expressed genes (DEGs) were enriched on the basis of pairwise comparisons (34 °C vs. 25 °C) at 0 (before inoculation), 3 and 6 days post-inoculation (dpi), respectively. A total of 7035 DEGs were identified from line LA3858 in the respective groups under the different soil temperature treatments. At 3 dpi, most DEGs were enriched in Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways related to plant biotic responses, such as “plant-pathogen interaction” and “plant hormone signal transduction”. Significantly enriched DEGs were found to encode key proteins such as R proteins and heat-shock proteins (HSPs). Moreover, other DEGs were found to participate in Ca²⁺ signal transduction; the production of ROS; DEGs encoding transcription factors (TFs) from the bHLH, TGA, ERF, heat-shock transcription factor (HSF) and WRKY families were highly expressed, which contribute to be involved into the formation of phytohormones, such as salicylic acid (SA), jasmonic acid (JA) and ethylene (ET), the expression of most was upregulated at 3 dpi at the 25 °C soil temperature compared with the 34 °C soil temperature.

Conclusion

Taken together, the results of our study revealed reliable candidate genes from wild materials LA3858, that are related to Mi-3-mediate resistance to Meloidogyne incognita. A large number of vital pathways and DEGs were expressed specifically in accession LA3858 grown at 34 °C and 25 °C soil temperatures at 3 dpi. Upon infection by RKNs, pattern-recognition receptors (PRRs) specifically recognized conserved pathogen-associated molecular patterns (PAMPs) as a result of pathogen-triggered immunity (PTI), and the downstream defensive signal transduction pathway was likely activated through Ca²⁺ signal channels. The expression of various TFs was induced to synthesize phytohormones and activate R proteins related to resistance, resulting in the development of effector-triggered immunity (ETI). Last, a hypersensitive response in the roots occurred, which was probably induced by the accumulation of ROS.

Leaf Senescence: The Chloroplast Connection Comes of Age

Leaf senescence is a developmental process critical for plant fitness, which involves genetically controlled cell death and ordered disassembly of macromolecules for reallocating nutrients to juvenile and reproductive organs. While natural leaf senescence is primarily associated with aging, it can also be induced by environmental and nutritional inputs including biotic and abiotic stresses, darkness, phytohormones and oxidants. Reactive oxygen species (ROS) are a common thread in stress-dependent cell death and also increase during leaf senescence. Involvement of chloroplast redox chemistry (including ROS propagation) in modulating cell death is well supported, with photosynthesis playing a crucial role in providing redox-based signals to this process. While chloroplast contribution to senescence received less attention, recent findings indicate that changes in the redox poise of these organelles strongly affect senescence timing and progress. In this review, the involvement of chloroplasts in leaf senescence execution is critically assessed in relation to available evidence and the role played by environmental and developmental cues such as stress and phytohormones. The collected results indicate that chloroplasts could cooperate with other redox sources (e.g., mitochondria) and signaling molecules to initiate the committed steps of leaf senescence for a best use of the recycled nutrients in plant reproduction.

Exploring gene networks in two sunflower lines with contrasting leaf senescence phenotype using a system biology approach

BACKGROUND

Leaf senescence is a complex process, controlled by multiple genetic and environmental variables. In sunflower, leaf senescence is triggered abruptly following anthesis thereby limiting the capacity of plants to keep their green leaf area during grain filling, which subsequently has a strong impact on crop yield. Recently, we performed a selection of contrasting sunflower inbred lines for the progress of leaf senescence through a physiological, cytological and molecular approach. Here we present a large scale transcriptomic analysis using RNA-seq and its integration with metabolic profiles for two contrasting sunflower inbred lines, R453 and B481-6 (early and delayed senescence respectively), with the aim of identifying metabolic pathways associated to leaf senescence.

RESULTS

Gene expression profiles revealed a higher number of differentially expressed genes, as well as, higher expression levels in R453, providing evidence for early activation of the senescence program in this line. Metabolic pathways associated with sugars and nutrient recycling were differentially regulated between the lines. Additionally, we identified transcription factors acting as hubs in the co-expression networks; some previously reported as senescence-associated genes in model species but many are novel candidate genes.

CONCLUSIONS

Understanding the onset and the progress of the senescence process in crops and the identification of these new candidate genes will likely prove highly useful for different management strategies to mitigate the impact of senescence on crop yield. Functional characterization of candidate genes will help to develop molecular tools for biotechnological applications in breeding crop yield.

Transcription Factors Associated with Leaf Senescence in Crops

Leaf senescence is a complex mechanism controlled by multiple genetic and environmental variables. Different crops present a delay in leaf senescence with an important impact on grain yield trough the maintenance of the photosynthetic leaf area during the reproductive stage. Additionally, because of the temporal gap between the onset and phenotypic detection of the senescence process, candidate genes are key tools to enable the early detection of this process. In this sense and given the importance of some transcription factors as hub genes in senescence pathways, we present a comprehensive review on senescence-associated transcription factors, in model plant species and in agronomic relevant crops. This review will contribute to the knowledge of leaf senescence process in crops, thus providing a valuable tool to assist molecular crop breeding.

Molecular profiling of mucinous epithelial ovarian cancer by weighted gene co-expression network analysis

PURPOSE

In order to identify the molecular characteristics and improve the efficacy of early diagnosis of mucinous epithelial ovarian cancer (mEOC), here, the transcriptome profiling by weighted gene co-expression network analysis (WGCNA) has been proposed as an effective method.

METHODS

The gene expression dataset GSE26193 was reanalyzed with a systematical approach, WGCNA. mEOC-related gene co-expression modules were detected and the functional enrichments of these modules were performed at GO and KEGG terms. Ten hub genes in the mEOC-related modules were validated using two independent datasets GSE44104 and GSE30274.

RESULTS

11 co-expressed gene modules were identified by WGCNA based on 4917 genes and 99 epithelial ovarian cancer samples. The turquoise module was found to be significantly associated with the subtype of mEOC. KEGG pathway enrichment analysis showed genes in the turquoise module significantly enriched in metabolism of xenobiotics by cytochrome P450 and steroid hormone biosynthesis. Ten hub genes (LIPH, BCAS1, FUT3, ZG16B, PTPRH, SLC4A4, MUC13, TFF1, HNF4G and TFF2) in the turquoise module were validated to be highly expressed in mEOC using two independent gene expression datasets GSE44104 and GSE30274.

CONCLUSION

Our work proposed an applicable framework of molecular characteristics for patients with mEOC, which may help us to obtain a precise and comprehensive understanding on the molecular complexities of mEOC. The hub genes identified in our study, as potential specific biomarkers of mEOC, may be applied in the early diagnosis of mEOC in the future.

Metabolomics of mammalian brain reveals regional differences

Background

The mammalian brain is organized into regions with specific biological functions and properties. These regions have distinct transcriptomes, but little is known whether they may also differ in their metabolome. The metabolome, a collection of small molecules or metabolites, is at the intersection of the genetic background of a given cell or tissue and the environmental influences that affect it. Thus, the metabolome directly reflects information about the physiologic state of a biological system under a particular condition. The objective of this study was to investigate whether various brain regions have diverse metabolome profiles, similarly to their genetic diversity. The answer to this question would suggest that not only the genome but also the metabolome may contribute to the functional diversity of brain regions.

Methods

We investigated the metabolome of four regions of the mouse brain that have very distinct functions: frontal cortex, hippocampus, cerebellum, and olfactory bulb. We utilized gas- and liquid- chromatography mass spectrometry platforms and identified 215 metabolites.

Results

Principal component analysis, an unsupervised multivariate analysis, clustered each brain region based on its metabolome content, thus providing the unique metabolic profile of each region. A pathway-centric analysis indicated that olfactory bulb and cerebellum had most distinct metabolic profiles, while the cortical parenchyma and hippocampus were more similar in their metabolome content. Among the notable differences were distinct oxidative-anti-oxidative status and region-specific lipid profiles. Finally, a global metabolic connectivity analysis using the weighted correlation network analysis identified five hub metabolites that organized a unique metabolic network architecture within each examined brain region. These data indicate the diversity of global metabolome corresponding to specialized regional brain function and provide a new perspective on the underlying properties of brain regions.

Conclusion

In summary, we observed many differences in the metabolome among the various brain regions investigated. All four brain regions in our study had a unique metabolic signature, but the metabolites came from all categories and were not pathway-centric.

Electronic supplementary material

The online version of this article (10.1186/s12918-018-0644-0) contains supplementary material, which is available to authorized users.

Full-Length Transcriptome Sequencing and Modular Organization Analysis of the Naringin/Neoeriocitrin-Related Gene Expression Pattern in Drynaria roosii

Drynaria roosii (Nakaike) is a traditional Chinese medicinal fern, known as 'GuSuiBu'. The effective components, naringin and neoeriocitrin, share a highly similar chemical structure and medicinal function. Our HPLC-tandem mass spectrometry (MS/MS) results showed that the accumulation of naringin/neoeriocitrin depended on specific tissues or ages. However, little was known about the expression patterns of naringin/neoeriocitrin-related genes involved in their regulatory pathways. Due to a lack of basic genetic information, we applied a combination of single molecule real-time (SMRT) sequencing and second-generation sequencing (SGS) to generate the complete and full-length transcriptome of D. roosii. According to the SGS data, the differentially expressed gene (DEG)-based heat map analysis revealed that naringin/neoeriocitrin-related gene expression exhibited obvious tissue- and time-specific transcriptomic differences. Using the systems biology method of modular organization analysis, we clustered 16,472 DEGs into 17 gene modules and studied the relationships between modules and tissue/time point samples, as well as modules and naringin/neoeriocitrin contents. We found that naringin/neoeriocitrin-related DEGs distributed in nine distinct modules, and DEGs in these modules showed significantly different patterns of transcript abundance to be linked to specific tissues or ages. Moreover, weighted gene co-expression network analysis (WGCNA) results further identified that PAL, 4CL and C4H, and C3H and HCT acted as the major hub genes involved in naringin and neoeriocitrin synthesis, respectively, and exhibited high co-expression with MYB- and basic helix-leucine-helix (bHLH)-regulated genes. In this work, modular organization and co-expression networks elucidated the tissue and time specificity of the gene expression pattern, as well as hub genes associated with naringin/neoeriocitrin synthesis in D. roosii. Simultaneously, the comprehensive transcriptome data set provided important genetic information for further research on D. roosii.

Integration of multi-omics techniques and physiological phenotyping within a holistic phenomics approach to study senescence in model and crop plants

The study of senescence in plants is complicated by diverse levels of temporal and spatial dynamics as well as the impact of external biotic and abiotic factors and crop plant management. Whereas the molecular mechanisms involved in developmentally regulated leaf senescence are very well understood, in particular in the annual model plant species Arabidopsis, senescence of other organs such as the flower, fruit, and root is much less studied as well as senescence in perennials such as trees. This review addresses the need for the integration of multi-omics techniques and physiological phenotyping into holistic phenomics approaches to dissect the complex phenomenon of senescence. That became feasible through major advances in the establishment of various, complementary 'omics' technologies. Such an interdisciplinary approach will also need to consider knowledge from the animal field, in particular in relation to novel regulators such as small, non-coding RNAs, epigenetic control and telomere length. Such a characterization of phenotypes via the acquisition of high-dimensional datasets within a systems biology approach will allow us to systematically characterize the various programmes governing senescence beyond leaf senescence in Arabidopsis and to elucidate the underlying molecular processes. Such a multi-omics approach is expected to also spur the application of results from model plants to agriculture and their verification for sustainable and environmentally friendly improvement of crop plant stress resilience and productivity and contribute to improvements based on postharvest physiology for the food industry and the benefit of its customers.

Sunflower Hybrid Breeding: From Markers to Genomic Selection

In sunflower, molecular markers for simple traits as, e.g., fertility restoration, high oleic acid content, herbicide tolerance or resistances to Plasmopara halstedii, Puccinia helianthi, or Orobanche cumana have been successfully used in marker-assisted breeding programs for years. However, agronomically important complex quantitative traits like yield, heterosis, drought tolerance, oil content or selection for disease resistance, e.g., against Sclerotinia sclerotiorum have been challenging and will require genome-wide approaches. Plant genetic resources for sunflower are being collected and conserved worldwide that represent valuable resources to study complex traits. Sunflower association panels provide the basis for genome-wide association studies, overcoming disadvantages of biparental populations. Advances in technologies and the availability of the sunflower genome sequence made novel approaches on the whole genome level possible. Genotype-by-sequencing, and whole genome sequencing based on next generation sequencing technologies facilitated the production of large amounts of SNP markers for high density maps as well as SNP arrays and allowed genome-wide association studies and genomic selection in sunflower. Genome wide or candidate gene based association studies have been performed for traits like branching, flowering time, resistance to Sclerotinia head and stalk rot. First steps in genomic selection with regard to hybrid performance and hybrid oil content have shown that genomic selection can successfully address complex quantitative traits in sunflower and will help to speed up sunflower breeding programs in the future. To make sunflower more competitive toward other oil crops higher levels of resistance against pathogens and better yield performance are required. In addition, optimizing plant architecture toward a more complex growth type for higher plant densities has the potential to considerably increase yields per hectare. Integrative approaches combining omic technologies (genomics, transcriptomics, proteomics, metabolomics and phenomics) using bioinformatic tools will facilitate the identification of target genes and markers for complex traits and will give a better insight into the mechanisms behind the traits.

Expression of a Plastid-Targeted Flavodoxin Decreases Chloroplast Reactive Oxygen Species Accumulation and Delays Senescence in Aging Tobacco Leaves

Leaf senescence is a concerted physiological process involving controlled degradation of cellular structures and reallocation of breakdown products to other plant organs. It is accompanied by increased production of reactive oxygen species (ROS) that are proposed to signal cell death, although both the origin and the precise role of ROS in the execution of this developmental program are still poorly understood. To investigate the contribution of chloroplast-associated ROS to natural leaf senescence, we used tobacco plants expressing a plastid-targeted flavodoxin, an electron shuttle flavoprotein present in prokaryotes and algae. When expressed in plants, flavodoxin specifically prevents ROS formation in chloroplasts during stress situations. Senescence symptoms were significantly mitigated in these transformants, with decreased accumulation of chloroplastic ROS and differential preservation of chlorophylls, carotenoids, protein contents, cell and chloroplast structures, membrane integrity and cell viability. Flavodoxin also improved maintenance of chlorophyll-protein complexes, photosynthetic electron flow, CO₂ assimilation, central metabolic routes and levels of bioactive cytokinins and auxins in aging leaves. Delayed induction of senescence-associated genes indicates that the entire genetic program of senescence was affected by flavodoxin. The results suggest that ROS generated in chloroplasts are involved in the regulation of natural leaf senescence.

Contrasting Transcriptional Programs Control Postharvest Development of Apples (Malus x domestica Borkh.) Submitted to Cold Storage and Ethylene Blockage

Apple is commercially important worldwide. Favorable genomic contexts and postharvest technologies allow year-round availability. Although ripening is considered a unidirectional developmental process toward senescence, storage at low temperatures, alone or in combination with ethylene blockage, is effective in preserving apple properties. Quality traits and genome wide expression were integrated to investigate the mechanisms underlying postharvest changes. Development and conservation techniques were responsible for transcriptional reprogramming and distinct programs associated with quality traits. A large portion of the differentially regulated genes constitutes a program involved in ripening and senescence, whereas a smaller module consists of genes associated with reestablishment and maintenance of juvenile traits after harvest. Ethylene inhibition was associated with a reversal of ripening by transcriptional induction of anabolic pathways. Our results demonstrate that the blockage of ethylene perception and signaling leads to upregulation of genes in anabolic pathways. We also associated complex phenotypes to subsets of differentially regulated genes.

Integration of transcriptomic and metabolic data reveals hub transcription factors involved in drought stress response in sunflower (Helianthus annuus L.)

By integration of transcriptional and metabolic profiles we identified pathways and hubs transcription factors regulated during drought conditions in sunflower, useful for applications in molecular and/or biotechnological breeding. Drought is one of the most important environmental stresses that effects crop productivity in many agricultural regions. Sunflower is tolerant to drought conditions but the mechanisms involved in this tolerance remain unclear at the molecular level. The aim of this study was to characterize and integrate transcriptional and metabolic pathways related to drought stress in sunflower plants, by using a system biology approach. Our results showed a delay in plant senescence with an increase in the expression level of photosynthesis related genes as well as higher levels of sugars, osmoprotectant amino acids and ionic nutrients under drought conditions. In addition, we identified transcription factors that were upregulated during drought conditions and that may act as hubs in the transcriptional network. Many of these transcription factors belong to families implicated in the drought response in model species. The integration of transcriptomic and metabolomic data in this study, together with physiological measurements, has improved our understanding of the biological responses during droughts and contributes to elucidate the molecular mechanisms involved under this environmental condition. These findings will provide useful biotechnological tools to improve stress tolerance while maintaining crop yield under restricted water availability.

MatureP: prediction of secreted proteins with exclusive information from their mature regions

More than a third of the cellular proteome is non-cytoplasmic. Most secretory proteins use the Sec system for export and are targeted to membranes using signal peptides and mature domains. To specifically analyze bacterial mature domain features, we developed MatureP, a classifier that predicts secretory sequences through features exclusively computed from their mature domains. MatureP was trained using Just Add Data Bio, an automated machine learning tool. Mature domains are predicted efficiently with ~92% success, as measured by the Area Under the Receiver Operating Characteristic Curve (AUC). Predictions were validated using experimental datasets of mutated secretory proteins. The features selected by MatureP reveal prominent differences in amino acid content between secreted and cytoplasmic proteins. Amino-terminal mature domain sequences have enhanced disorder, more hydroxyl and polar residues and less hydrophobics. Cytoplasmic proteins have prominent amino-terminal hydrophobic stretches and charged regions downstream. Presumably, secretory mature domains comprise a distinct protein class. They balance properties that promote the necessary flexibility required for the maintenance of non-folded states during targeting and secretion with the ability of post-secretion folding. These findings provide novel insight in protein trafficking, sorting and folding mechanisms and may benefit protein secretion biotechnology.