Characterization of Hsp17, a Novel Small Heat Shock Protein, in Sphingomonas melonis TY under Heat Stress

Bacteria are constantly exposed to a variety of environmental stresses. Temperature is considered one of the most important environmental factors affecting microbial growth and survival. As ubiquitous environmental microorganisms, Sphingomonas species play essential roles in the biodegradation of organic contaminants, plant protection, and environmental remediation. Understanding the mechanism by which they respond to heat shock will help further improve cell resistance by applying synthetic biological strategies. Here, we assessed the transcriptomic and proteomic responses of Sphingomonas melonis TY to heat shock and found that stressful conditions caused significant changes in functional genes related to protein synthesis at the transcriptional level. The most notable changes observed were increases in the transcription (1,857-fold) and protein expression (11-fold) of Hsp17, which belongs to the small heat shock protein family, and the function of Hsp17 in heat stress was further investigated in this study. We found that the deletion of hsp17 reduced the capacity of the cells to tolerate high temperatures, whereas the overexpression of hsp17 significantly enhanced the ability of the cells to withstand high temperatures. Moreover, the heterologous expression of hsp17 in Escherichia coli DH5α conferred to the bacterium the ability to resist heat stress. Interestingly, its cells were elongated and formed connected cells following the increase in temperature, while hsp17 overexpression restored their normal morphology under high temperature. In general, these results indicate that the novel small heat shock protein Hsp17 greatly contributes to maintaining cell viability and morphology under stress conditions. IMPORTANCE Temperature is generally considered the most important factor affecting metabolic functions and the survival of microbes. As molecular chaperones, small heat shock proteins can prevent damaged protein aggregation during abiotic stress, especially heat stress. Sphingomonas species are widely distributed in nature, and they can frequently be found in various extreme environments. However, the role of small heat shock proteins in Sphingomonas under high-temperature stress has not been elucidated. This study greatly enhances our understanding of a novel identified protein, Hsp17, in S. melonis TY in terms of its ability to resist heat stress and maintain cell morphology under high temperature, leading to a broader understanding of how microbes adapt to environmental extremes. Furthermore, our study will provide potential heat resistance elements for further enhancing cellular resistance as well as the synthetic biological applications of Sphingomonas.

Long-term cold, freezing and drought: overlapping and specific regulatory mechanisms and signal transduction in tea plant (Camellia sinensis (L.) Kuntze)

Introduction

Low temperatures and drought are two main environmental constraints reducing the yield and geographical distribution of horticultural crops worldwide. Understanding the genetic crosstalk between stress responses has potential importance for crop improvement.

Methods

In this study, Illumina RNA-seq and Pac-Bio genome resequencing were used to annotate genes and analyze transcriptome dynamics in tea plants under long-term cold, freezing, and drought.

Results

The highest number of differentially expressed genes (DEGs) was identified under long-term cold (7,896) and freezing (7,915), with 3,532 and 3,780 upregulated genes, respectively. The lowest number of DEGs was observed under 3-day drought (47) and 9-day drought (220), with five and 112 genes upregulated, respectively. The recovery after the cold had 6.5 times greater DEG numbers as compared to the drought recovery. Only 17.9% of cold-induced genes were upregulated by drought. In total, 1,492 transcription factor genes related to 57 families were identified. However, only 20 transcription factor genes were commonly upregulated by cold, freezing, and drought. Among the 232 common upregulated DEGs, most were related to signal transduction, cell wall remodeling, and lipid metabolism. Co-expression analysis and network reconstruction showed 19 genes with the highest co-expression connectivity: seven genes are related to cell wall remodeling (GATL7, UXS4, PRP-F1, 4CL, UEL-1, UDP-Arap, and TBL32), four genes are related to calcium-signaling (PXL1, Strap, CRT, and CIPK6), three genes are related to photo-perception (GIL1, CHUP1, and DnaJ11), two genes are related to hormone signaling (TTL3 and GID1C-like), two genes are involved in ROS signaling (ERO1 and CXE11), and one gene is related to the phenylpropanoid pathway (GALT6).

Discussion

Based on our results, several important overlapping mechanisms of long-term stress responses include cell wall remodeling through lignin biosynthesis, o-acetylation of polysaccharides, pectin biosynthesis and branching, and xyloglucan and arabinogalactan biosynthesis. This study provides new insight into long-term stress responses in woody crops, and a set of new target candidate genes were identified for molecular breeding aimed at tolerance to abiotic stresses.

Metabolomics and Transcriptomics Analyses Reveal Regulatory Networks Associated with Fatty Acid Accumulation in Pecan Kernels

Pecans are a globally important tree nut crop. Pecan nuts are rich in fatty acids (FAs), proteins, and flavonoids in addition to thiamine and numerous micronutrients. Although several of these nutriments have been studied in this plant, the comprehensive metabolite variations and molecular mechanisms associated with them have not been fully elucidated. In this study, untargeted metabolomics and transcriptomics were integrated to reveal the metabolite accumulation patterns and their associated molecular mechanisms during pecan kernel development. In total, 4260 (under positive mode) and 2726 (under negative mode) high quality features were retained. Overall, 163 differentially accumulated metabolites were identified. Most components were classified into the categories "organic acids and derivatives" and "lipids and lipid-like molecules." The accumulation patterns of amino acids, FAs, carbohydrates, organic acids, vitamins, flavonoids, and phenylpropanoids alongside embryo development were determined. Furthermore, transcriptomes from four pecan kernel developmental stages were used to assess transcript expression levels. Coexpression analyses were performed between FAs and their related genes. This study provides a comprehensive overview of the metabolic changes and regulations during pecan kernel development. We believe that the identification of nutriment accumulation trends and hub genes associated with the biosynthesis of the components will be valuable for genetically improving this plant.

Optimization of Heterologous Glucoraphanin Production In Planta

Glucoraphanin is a plant specialized metabolite found in cruciferous vegetables that has long been a target for production in a heterologous host because it can subsequently be hydrolyzed to form the chemopreventive compound sulforaphane before and during consumption. However, previous studies have only been able to produce small amounts of glucoraphanin in heterologous plant and microbial systems compared to the levels found in glucoraphanin-producing plants, suggesting that there may be missing auxiliary genes that play a role in improving production in planta. In an effort to identify auxiliary genes required for high glucoraphanin production, we leveraged transient expression in Nicotiana benthamiana to screen a combination of previously uncharacterized coexpressed genes and rationally selected genes alongside the glucoraphanin biosynthetic pathway. This strategy alleviated metabolic bottlenecks, which improved glucoraphanin production by 4.74-fold. Our optimized glucoraphanin biosynthetic pathway provides a pathway amenable for high glucoraphanin production.

Epilepsy Combined With Multiple Gene Heterozygous Mutation

The fast pace of gene discovery has resulted in groundbreaking advances in the field of epilepsy genetics. Clinical testing using comprehensive gene panels, exomes, or genomes is now increasingly available and has significantly increased the diagnostic yield for early-onset epilepsies and enabled precision medicine approaches. In this paper, we report a case of epilepsy in a pedigree. The proband had heterozygous mutations in KCNC1 (NM_001112741.1:c.959G>A, p. Arg320His), CAPN3 (NM_000070.2:c.526G>A, p. Val176Met), and NEFH (NM_021076.3:c. 2595 delC, p. Lys866Argfs^*51). Sanger sequencing verification was consistent with the results of whole-exome sequencing. The KCNC1 mutation was a de novo mutation, and the CAPN3 and NEFH mutations were inherited from their father and mother, respectively. Based on the American College of Medical Genetics and Genomics (ACMG) guidelines, a heterozygous mutation was found for APOB (NM_000384.2: c.10579C > T, p. Arg3527Trp). The heterozygous mutation at this site was inherent in the pedigree. Coexpression analysis indicated that heterozygous mutations of KCNC1, CAPN3, NEFH, and APOB were closely related to the clinical phenotypes of the patient, and the clinical phenotypic heterogeneity of the disease may be the result of the interaction of multiple genes.

Genome-wide identification of seed storage protein gene regulators in wheat through coexpression analysis

Only a few transcriptional regulators of seed storage protein (SSP) genes have been identified in common wheat (Triticum aestivum L.). Coexpression analysis could be an efficient approach to characterize novel transcriptional regulators at the genome-scale considering the correlated expression between transcriptional regulators and target genes. As the A genome donor of common wheat, Triticum urartu is more suitable for coexpression analysis than common wheat considering the diploid genome and single gene copy. In this work, the transcriptome dynamics in endosperm of T. urartu throughout grain filling were revealed by RNA-Seq analysis. In the coexpression analysis, a total of 71 transcription factors (TFs) from 23 families were found to be coexpressed with SSP genes. Among these TFs, TuNAC77 enhanced the transcription of SSP genes by binding to cis-elements distributed in promoters. The homolog of TuNAC77 in common wheat, TaNAC77, shared an identical function, and the total SSPs were reduced by about 24% in common wheat when TaNAC77 was knocked down. This is the first genome-wide identification of transcriptional regulators of SSP genes in wheat, and the newly characterized transcriptional regulators will undoubtedly expand our knowledge of the transcriptional regulation of SSP synthesis.

De novo Transcriptome Sequencing Coupled With Co-expression Analysis Reveal the Transcriptional Regulation of Key Genes Involved in the Formation of Active Ingredients in Peucedanum praeruptorum Dunn Under Bolting Period

Peucedanum praeruptorum Dunn is a perennial and one-off flowering plant of the Peucedanum genus in Umbelliferae. The cultivated P. praeruptorum Dunn usually grows nutritionally in the first year and then moves into the reproductive growth in the second year. The lignification of the roots caused by bolting leads to the quality decline of crude materials. Since most of the previous studies have dealt with coumarin biosynthesis and identification of functional genes in P. praeruptorum, the scientific connotation of the inability that the bolted P. praeruptorum cannot be used medically is still unclear. Here, we employed a transcriptome sequencing combined with coexpression analysis to unearth the regulation mechanism of key genes related to coumarin synthesis in pre- and postbolting period, and to explore the mechanisms underlying the effects of bolting on the formation and transport of coumarins between the annual and biennial plants. Six cDNA libraries were constructed, and the transcripts were sequenced and assembled by Illumina Hiseq platform. A total of 336,505 unigenes were obtained from 824,129 non-redundant spliced transcripts. Unigenes (114,488) were annotated to the NCBI nr database, 119,017 and 10,475 unigenes were aligned to Gene Ontology (GO) functional groups and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, respectively. Differential expression analysis screened out a series of upregulated and downregulated genes related to the phenylpropanoid pathway. The heatmap clustering showed that the similar expression patterns were both observed in groups C vs. D and groups C vs. F. The WGCNA-based coexpression was performed to elucidate the module and trait relationship to unearth important genes related to the bolting process. Seven pivotal modules on the KEGG functional annotations suggested these genes were mainly enriched in the process of plant–pathogen interaction, plant hormone signal transduction, MAPK signaling pathway, α-linolenic acid metabolism, circadian rhythm, and phenylpropanoid pathway. Further analysis provided clues that the key genes of the phenylpropanoid pathway, the ABC transporters, the apoptosis-related and circadian rhythm regulatory genes may play pivotal roles in regulating bolting signaling, biosynthesis, and transportation of coumarins.

Comprehensive Analysis of Cell Population Dynamics and Related Core Genes During Vitiligo Development

Vitiligo is a common immune-related depigmentation condition, and its pathogenesis remains unclear. This study used a combination of bioinformatics methods and expression analysis techniques to explore the relationship between immune cell infiltration and gene expression in vitiligo. Previously reported gene expression microarray data from the skin (GSE53146 and GSE75819) and peripheral blood (GSE80009 and GSE90880) of vitiligo patients and healthy controls was used in the analysis. R software was used to filter the differentially expressed genes (DEGs) in each dataset, and the KOBAS 2.0 server was used to perform functional enrichment analysis. Compared with healthy controls, the upregulated genes in skin lesions and peripheral blood leukocytes of vitiligo patents were highly enriched in immune response pathways and inflammatory response signaling pathways. Immunedeconv software and the EPIC method were used to analyze the expression levels of marker genes to obtain the immune cell population in the samples. In the lesional skin of vitiligo patients, the proportions of macrophages, B cells and NK cells were increased compared with healthy controls. In the peripheral blood of vitiligo patients, CD8+ T cells and macrophages were significantly increased. A coexpression analysis of the cell populations and DEGs showed that differentially expressed immune and inflammation response genes had a strong positive correlation with macrophages. The TLR4 receptor pathway, interferon gamma-mediated signaling pathway and lipopolysaccharide-related pathway were positively correlated with CD4+ T cells. Regarding immune response-related genes, the overexpression of IFITM2, TNFSF10, GZMA, ADAMDEC1, NCF2, ADAR, SIGLEC16, and WIPF2 were related to macrophage abundance, while the overexpression of ICOS, GPR183, RGS1, ILF2 and CD28 were related to CD4+ T cell abundance. GZMA and CXCL10 expression were associated with CD8+ T cell abundance. Regarding inflammatory response-related genes, the overexpression of CEBPB, ADAM8, CXCR3, and TNIP3 promoted macrophage infiltration. Only ADORA1 expression was associated with CD4+ T cell infiltration. ADAM8 and CXCL10 expression were associated with CD8+ T cell abundance. The overexpression of CCL18, CXCL10, FOS, NLRC4, LY96, HCK, MYD88, and KLRG1, which are related to inflammation and immune responses, were associated with macrophage abundance. We also found that immune cells infiltration in vitiligo was associated with antigen presentation-related genes expression. The genes and pathways identified in this study may point to new directions for vitiligo treatment.

Coexpression analysis of a large-scale transcriptome identified a calmodulin-like protein regulating the development of adventitious roots in poplar

Poplars are important woody plants, and the ability to form adventitious roots (ARs) is the key factor for their cultivation because most poplars are propagated by cloning. In previous studies, Ca2+ was confirmed to regulate AR formation in poplar. In this study, wild-type poplar cuttings grown in 1.0 mM Ca2+ solution showed the best visible performance of AR development. Coexpression analysis of a large-scale RNA-Seq transcriptome was conducted to identify Ca2+-related genes that regulate AR development in poplar. A total of 15 coexpression modules (CMs) were identified, and two CMs showed high association with AR development. Functional analysis identified a number of biological pathways, including 'oxidation-reduction process', 'response to biotic stimulus' and 'metabolic process', in tissues of AR development. The Ca2+-related pathway was specifically selected, and its regulation in poplar AR development was predicted. A Ca2+ sensor, PdeCML23-1, which is a member of the calmodulin-like protein (CML) family, was found to promote AR development by phenotypic assay of overexpressed PdeCML23-1 transgenic lines at various growing conditions. By measuring cytosolic Ca2+ in AR tips, PdeCML23-1 seemed to play a role in decreasing cytosolic Ca2+ concentration. Additionally, the expression profiles of some genes and phytohormone indole acetic acid (IAA) were also changed in the overexpressed PdeCML23-1 transgenic lines. According to this study, we were able to provide a global view of gene regulation for poplar AR development. Moreover, we also observed the regulation of cytosolic Ca2+ concentration by PdeCML23-1, and this regulation was involved in AR development in poplar. We also predicted that PdeCML23-1 possibly regulates AR development by modulating IAA content in poplar.

Network-based Transcriptome-wide Expression Study for Postmenopausal Osteoporosis

PURPOSE

Menopause is a crucial physiological transition during a woman's life, and it occurs with growing risks of health issues like osteoporosis. To identify postmenopausal osteoporosis-related genes, we performed transcriptome-wide expression analyses for human peripheral blood monocytes (PBMs) using Affymetrix 1.0 ST arrays in 40 Caucasian postmenopausal women with discordant bone mineral density (BMD) levels.

METHODS

We performed multiscale embedded gene coexpression network analysis (MEGENA) to study functionally orchestrating clusters of differentially expressed genes in the form of functional networks. Gene sets net correlations analysis (GSNCA) was applied to assess how the coexpression structure of a predefined gene set differs in high and low BMD groups. Bayesian network (BN) analysis was used to identify important regulation patterns between potential risk genes for osteoporosis. A small interfering ribonucleic acid (siRNA)-based gene silencing in vitro experiment was performed to validate the findings from BN analysis.

RESULT

MEGENA showed that the "T cell receptor signaling pathway" and the "osteoclast differentiation pathway" were significantly enriched in the identified compact network, which is significantly correlated with BMD variation. GSNCA revealed that the coexpression structure of the "Signaling by TGF-beta receptor complex pathway" is significantly different between the 2 BMD discordant groups; the hub genes in the postmenopausal low and high BMD group are FURIN and SMAD3 respectively. With siRNA in vitro experiments, we confirmed the regulation relationship of TGFBR2-SMAD7 and TGFBR1-SMURF2.

MAIN CONCLUSION

The present study suggests that biological signals involved in monocyte recruitment, monocyte/macrophage lineage development, osteoclast formation, and osteoclast differentiation might function together in PBMs that contribute to the pathogenesis of postmenopausal osteoporosis.

Meta-Analysis of the Alzheimer's Disease Human Brain Transcriptome and Functional Dissection in Mouse Models

We present a consensus atlas of the human brain transcriptome in Alzheimer's disease (AD), based on meta-analysis of differential gene expression in 2,114 postmortem samples. We discover 30 brain coexpression modules from seven regions as the major source of AD transcriptional perturbations. We next examine overlap with 251 brain differentially expressed gene sets from mouse models of AD and other neurodegenerative disorders. Human-mouse overlaps highlight responses to amyloid versus tau pathology and reveal age- and sex-dependent expression signatures for disease progression. Human coexpression modules enriched for neuronal and/or microglial genes broadly overlap with mouse models of AD, Huntington's disease, amyotrophic lateral sclerosis, and aging. Other human coexpression modules, including those implicated in proteostasis, are not activated in AD models but rather following other, unexpected genetic manipulations. Our results comprise a cross-species resource, highlighting transcriptional networks altered by human brain pathophysiology and identifying correspondences with mouse models for AD preclinical studies.

Identification of PIGU as the Hub Gene Associated with KRAS Mutation in Colorectal Cancer by Coexpression Analysis

Colorectal cancer (CRC) patients with KRAS mutation are refractory and usually have poor prognosis. We aimed to identify the hub gene associated with KRAS mutant CRCs. Weighted gene coexpression network analysis (WGCNA) was used to calculate the key module and the hub genes in GSE39582. Combined with the protein-protein interaction (PPI) network and survival analysis, the real hub gene was identified and further validated. With the highest module significance value and correlation coefficient, the blue module was selected as the key module, 19 genes were identified as the hub gene candidates. The above genes were significantly downregulated in KRAS mutant CRCs compared with the wild type. Four genes (AAR2, PSMA7, NELFCD, and PIGU) were further screened as the potential hub genes by the PPI network. Low expression of PIGU for KRAS mutant patients had a poor prognosis. Therefore, PIGU was identified as the hub gene. PIGU expression was also downregulated in other two CRC datasets. "MAPK SIGNALING PATHWAY" was enriched in PIGU lowly expressed samples. PIGU was identified and validated to be closely related to KRAS mutation. It could be a potential prognosis biomarker and a novel treatment target for KRAS mutant CRC patients.

Comparison of Metabolome and Transcriptome of Flavonoid Biosynthesis Pathway in a Purple-Leaf Tea Germplasm Jinmingzao and a Green-Leaf Tea Germplasm Huangdan reveals Their Relationship with Genetic Mechanisms of Color Formation

Purple-leaf tea is a phenotype with unique color because of its high anthocyanin content. The special flavor of purple-leaf tea is highly different from that of green-leaf tea, and its main ingredient is also of economic value. To probe the genetic mechanism of the phenotypic characteristics of tea leaf color, we conducted widely targeted metabolic and transcriptomic profiling. The metabolites in the flavonoid biosynthetic pathway of purple- and green-leaf tea were compared, and results showed that phenolic compounds, including phenolic acids, flavonoids, and tannins, accumulated in purple-leaf tea. The high expression of genes related to flavonoid biosynthesis (e.g., PAL and LAR) exhibits the specific expression of biosynthesis and the accumulation of these metabolites. Our result also shows that two CsUFGTs were positively related to the accumulation of anthocyanin. Moreover, genes encoding transcription factors that regulate flavonoids were identified by coexpression analysis. These results may help to identify the metabolic factors that influence leaf color differentiation and provide reference for future research on leaf color biology and the genetic improvement of tea.

Genome-Wide Identification and Functional Analysis of NADPH Oxidase Family Genes in Wheat During Development and Environmental Stress Responses

As the key producers of reactive oxygen species (ROS), NADPH oxidases (NOXs), also known as respiratory burst oxidase homologs (RBOHs), play crucial roles in various biological processes in plants with considerable evolutionary selection and functional diversity in the entire terrestrial plant kingdom. However, only limited resources are available on the phylogenesis and functions of this gene family in wheat. Here, a total of 46 NOX family genes were identified in the wheat genome, and these NOXs could be classified into three subgroups: typical TaNOXs, TaNOX-likes, and ferric reduction oxidases (TaFROs). Phylogenetic analysis indicated that the typical TaNOXs might originate from TaFROs during evolution, and the TaFROs located on Chr 2 might be the most ancient forms of TaNOXs. TaNOXs are highly expressed in wheat with distinct tissue or organ-specificity and stress-inducible diversity. A large-scale expression and/or coexpression analysis demonstrated that TaNOXs can be divided into four functional groups with different expression patterns under a broad range of environmental stresses. Different TaNOXs are coexpressed with different sets of other genes, which widely participate in several important intracellular processes such as cell wall biosynthesis, defence response, and signal transduction, suggesting their vital but diversity of roles in plant growth regulation and stress responses of wheat.

Overexpression of SbMyb60 in Sorghum bicolor impacts both primary and secondary metabolism

Few transcription factors have been identified in C₄ grasses that either positively or negatively regulate monolignol biosynthesis. Previously, the overexpression of SbMyb60 in sorghum (Sorghum bicolor) has been shown to induce monolignol biosynthesis, which leads to elevated lignin deposition and altered cell wall composition. To determine how SbMyb60 overexpression impacts other metabolic pathways, RNA-Seq and metabolite profiling were performed on stalks and leaves. 35S::SbMyb60 was associated with the transcriptional activation of genes involved in aromatic amino acid, S-adenosyl methionine (SAM) and folate biosynthetic pathways. The high coexpression values between SbMyb60 and genes assigned to these pathways indicate that SbMyb60 may directly induce their expression. In addition, 35S::SbMyb60 altered the expression of genes involved in nitrogen (N) assimilation and carbon (C) metabolism, which may redirect C and N towards monolignol biosynthesis. Genes linked to UDP-sugar biosynthesis and cellulose synthesis were also induced, which is consistent with the observed increase in cellulose deposition in the internodes of 35S::SbMyb60 plants. However, SbMyb60 showed low coexpression values with these genes and is not likely to be a direct regulator of cell wall polysaccharide biosynthesis. These findings indicate that SbMyb60 can activate pathways beyond monolignol biosynthesis, including those that synthesize the substrates and cofactors required for lignin biosynthesis.

New BAR tools for mining expression data and exploring Cis-elements in Arabidopsis thaliana

Identifying sets of genes that are specifically expressed in certain tissues or in response to an environmental stimulus is useful for designing reporter constructs, generating gene expression markers, or for understanding gene regulatory networks. We have developed an easy-to-use online tool for defining a desired expression profile (a modification of our Expression Angler program), which can then be used to identify genes exhibiting patterns of expression that match this profile as closely as possible. Further, we have developed another online tool, Cistome, for predicting or exploring cis-elements in the promoters of sets of co-expressed genes identified by such a method, or by other methods. We present two use cases for these tools, which are freely available on the Bio-Analytic Resource at http://BAR.utoronto.ca.

Expression of p21-activated kinases 1 and 3 is altered in the brain of subjects with depression

The p21-activated kinases (PAKs) of group I are the main effectors for the small Rho GTPases, critically involved in neurodevelopment, plasticity and maturation of the nervous system. Moreover, the neuronal complexity controlled by PAK1/PAK3 signaling determines the postnatal brain size and synaptic properties. Stress induces alterations at the level of structural and functional synaptic plasticity accompanied by reductions in size and activity of the hippocampus and the prefrontal cortex (PFC). These abnormalities are likely to contribute to the pathology of depression and, in part, reflect impaired cytoskeleton remodeling pointing to the role of Rho GTPase signaling. Thus, the present study assessed the expression of the group I PAKs and their activators in the brain of depressed subjects. Using quantitative polymerase chain reaction (qPCR), mRNA levels and coexpression of the group I PAKs: PAK1, PAK2, and PAK3 as well as of their activators: RAC1, CDC42 and ARHGEF7 were examined in postmortem samples from the PFC (n=25) and the hippocampus (n=23) of subjects with depression and compared to control subjects (PFC n=24; hippocampus n=21). Results demonstrated that mRNA levels of PAK1 and PAK3, are significantly reduced in the brain of depressed subjects, with PAK1 being reduced in the PFC and PAK3 in the hippocampus. No differences were observed for the ubiquitously expressed PAK2. Following analysis of gene coexpression demonstrated disruption of coordinated gene expression in the brain of subjects with depression. Abnormalities in mRNA expression of PAK1 and PAK3 as well as their altered coexpression patterns were detected in the brain of subjects with depression.

Coexpression Network Analysis of Macronutrient Deficiency Response Genes in Rice

BACKGROUND

Macronutrients are pivotal elements for proper plant growth and development. Although extensive gene expression profiling revealed a large number of genes differentially expressed under various nutrient deprivation, characterization of these genes has never been fully explored especially in rice. Coexpression network analysis is a useful tool to elucidate the functional relationships of genes based on common expression. Therefore, we performed microarray analysis of rice shoot under nitrogen (N), phosphorus (P), and potassium (K) deficiency conditions. Moreover, we conducted a large scale coexpression analysis by integrating the data with previously generated gene expression profiles of organs and tissues at different developmental stages to obtain a global view of gene networks associated with plant response to nutrient deficiency.

RESULTS

We statistically identified 5400 differentially expressed genes under the nutrient deficiency treatments. Subsequent coexpression analysis resulted in the extraction of 6 modules (groups of highly interconnected genes) with distinct gene expression signatures. Three of these modules comprise mostly of downregulated genes under N deficiency associated with distinct functions such as development of immature organs, protein biosynthesis and photosynthesis in chloroplast of green tissues, and fundamental cellular processes in all organs and tissues. Furthermore, we identified one module containing upregulated genes under N and K deficiency conditions, and a number of genes encoding protein kinase, kinase-like domain containing protein and nutrient transporters. This module might be particularly involved in adaptation to nutrient deficiency via phosphorylation-mediated signal transduction and/or post-transcriptional regulation.

CONCLUSIONS

Our study demonstrated that large scale coexpression analysis is an efficient approach in characterizing the nutrient response genes based on biological functions and could provide new insights in understanding plant response to nutrient deficiency.

Altered expression of neuroplasticity-related genes in the brain of depressed suicides

BACKGROUND

Expression of the neuronal membrane glycoprotein M6a (GPM6A), the proteolipid protein (PLP/DM20) family member, is downregulated in the hippocampus of chronically stressed animals. Its neuroplastic function involves a role in neurite formation, filopodium outgrowth and synaptogenesis through an unknown mechanism. Disruptions in neuroplasticity mechanisms have been shown to play a significant part in the etiology of depression. Thus, the current investigation examined whether GPM6A expression is also altered in human depressed brain.

METHODS

Expression levels and coexpression patterns of GPM6A, GPM6B, and PLP1 (two other members of PLP/DM20 family) as well as of the neuroplasticity-related genes identified to associate with GPM6A were determined using quantitative polymerase chain reaction (qPCR) in postmortem samples from the hippocampus (n = 18) and the prefrontal cortex (PFC) (n = 25) of depressed suicide victims and compared with control subjects (hippocampus n = 18; PFC n = 25). Neuroplasticity-related proteins that form complexes with GPM6A were identified by coimmunoprecipitation technique followed by mass spectrometry.

RESULTS

Results indicated transcriptional downregulation of GPM6A and GPM6B in the hippocampus of depressed suicides. The expression level of calcium/calmodulin-dependent protein kinase II alpha (CAMK2A) and coronin1A (CORO1A) was also significantly decreased. Subsequent analysis of coexpression patterns demonstrated coordinated gene expression in the hippocampus and in the PFC indicating that the function of these genes might be coregulated in the human brain. However, in the brain of depressed suicides this coordinated response was disrupted.

CONCLUSIONS

Disruption of coordinated gene expression as well as abnormalities in GPM6A and GPM6B expression and expression of the components of GPM6A complexes were detected in the brain of depressed suicides.

Current challenges and future potential of tomato breeding using omics approaches

As tomatoes are one of the most important vegetables in the world, improvements in the quality and yield of tomato are strongly required. For this purpose, omics approaches such as metabolomics and transcriptomics are used not only for basic research to understand relationships between important traits and metabolism but also for the development of next generation breeding strategies of tomato plants, because an increase in the knowledge improves the taste and quality, stress resistance and/or potentially health-beneficial metabolites and is connected to improvements in the biochemical composition of tomatoes. Such omics data can be applied to network analyses to potentially reveal unknown cellular regulatory networks in tomato plants. The high-quality tomato genome that was sequenced in 2012 will likely accelerate the application of omics strategies, including next generation sequencing for tomato breeding. In this review, we highlight the current studies of omics network analyses of tomatoes and other plant species, in particular, a gene coexpression network. Key applications of omics approaches are also presented as case examples to improve economically important traits for tomato breeding.

Correlation networks visualization

New, in silico ways of generating hypotheses based on large data sets have emerged in the past decade. These data sets have been used to investigate different aspects of plant biology, especially at the level of transcriptome, from tissue-specific expression patterns to patterns in as little as a few cells. Such publicly available data are a boon to researchers for hypothesis generation by providing a guide for experimental work such as phenotyping or genetic analysis. More advanced computational methods can leverage these data via gene coexpression analysis, the results of which can be visualized and refined using network analysis. Other kinds of networks of, e.g., protein-protein interactions, can also be used to inform biology. These networks can be visualized and analyzed with additional information on gene expression levels, subcellular localization, etc., or with other emerging kinds information. Finally, cross-level correlation is an area that will become increasingly important. Visualizing these cross-level correlations will require new data visualization tools.