Gene Network Modules Research Articles

Integrative Multiomics in the Lung Reveals a Protective Role of Asporin in Pulmonary Arterial Hypertension.

Integrative multiomics can elucidate pulmonary arterial hypertension (PAH) pathobiology, but procuring human PAH lung samples is rare. We leveraged transcriptomic profiling and deep phenotyping of the largest multicenter PAH lung biobank to date (96 disease and 52 control) by integration with clinicopathologic data, genome-wide association studies, Bayesian regulatory networks, single-cell transcriptomics, and pharmacotranscriptomics. We identified 2 potentially protective gene network modules associated with vascular cells, and we validated ASPN, coding for asporin, as a key hub gene that is upregulated as a compensatory response to counteract PAH. We found that asporin is upregulated in lungs and plasma of multiple independent PAH cohorts and correlates with reduced PAH severity. We show that asporin inhibits proliferation and transforming growth factor-β/phosphorylated SMAD2/3 signaling in pulmonary artery smooth muscle cells from PAH lungs. We demonstrate in Sugen-hypoxia rats that ASPN knockdown exacerbated PAH and recombinant asporin attenuated PAH. Our integrative systems biology approach to dissect the PAH lung transcriptome uncovered asporin as a novel protective target with therapeutic potential in PAH.

Unveiling Gene Interactions in Alzheimer's Disease by Integrating Genetic and Epigenetic Data with a Network-Based Approach.

Alzheimer's Disease (AD) is a complex disease and the leading cause of dementia in older people. We aimed to uncover aspects of AD's pathogenesis that may contribute to drug repurposing efforts by integrating DNA methylation and genetic data. Implementing the network-based tool, a dense module search of genome-wide association studies (dmGWAS), we integrated a large-scale GWAS dataset with DNA methylation data to identify gene network modules associated with AD. Our analysis yielded 286 significant gene network modules. Notably, the foremost module included the BIN1 gene, showing the largest GWAS signal, and the GNAS gene, the most significantly hypermethylated. We conducted Web-based Cell-type-Specific Enrichment Analysis (WebCSEA) on genes within the top 10% of dmGWAS modules, highlighting monocyte as the most significant cell type (p < 5 × 10-12). Functional enrichment analysis revealed Gene Ontology Biological Process terms relevant to AD pathology (adjusted p < 0.05). Additionally, drug target enrichment identified five FDA-approved targets (p-value = 0.03) for further research. In summary, dmGWAS integration of genetic and epigenetic signals unveiled new gene interactions related to AD, offering promising avenues for future studies.

Dissecting gene expression networks in the developing hippocampus through the lens of NEIL3 depletion

Gene regulation in the hippocampus is fundamental for its development, synaptic plasticity, memory formation, and adaptability. Comparisons of gene expression among different developmental stages, distinct cell types, and specific experimental conditions have identified differentially expressed genes contributing to the organization and functionality of hippocampal circuits. The NEIL3 DNA glycosylase, one of the DNA repair enzymes, plays an important role in hippocampal maturation and neuron functionality by shaping transcription. While differential gene expression (DGE) analysis has identified key genes involved, broader gene expression patterns crucial for high-order hippocampal functions remain uncharted. By utilizing the weighted gene co-expression network analysis (WGCNA), we mapped gene expression networks in immature (p8-neonatal) and mature (3 m-adult) hippocampal circuits in wild-type and NEIL3-deficient mice. Our study unveiled intricate gene network structures underlying hippocampal maturation, delineated modules of co-expressed genes, and pinpointed highly interconnected hub genes specific to the maturity of hippocampal subregions. We investigated variations within distinct gene network modules following NEIL3 depletion, uncovering NEIL3-targeted hub genes that influence module connectivity and specificity. By integrating WGCNA with DGE, we delve deeper into the NEIL3-dependent molecular intricacies of hippocampal maturation. This study provides a comprehensive systems-level analysis for assessing the potential correlation between gene connectivity and functional connectivity within the hippocampal network, thus shaping hippocampal function throughout development.

Modulation of neural gene networks by estradiol in old rhesus macaque females.

The postmenopausal decrease in circulating estradiol (E2) levels has been shown to contribute to several adverse physiological and psychiatric effects. To elucidate the molecular effects of E2 on the brain, we examined differential gene expression and DNA methylation (DNAm) patterns in the nonhuman primate brain following ovariectomy (Ov) and subsequent subcutaneous bioidentical E2 chronic treatment. We identified several dysregulated molecular networks, including MAPK signaling and dopaminergic synapse response, that are associated with ovariectomy and shared across two different brain areas, the occipital cortex (OC) and prefrontal cortex (PFC). The finding that hypomethylation (p = 1.6 × 10-51) and upregulation (p = 3.8 × 10-3) of UBE2M across both brain regions provide strong evidence for molecular differences in the brain induced by E2 depletion. Additionally, differential expression (p = 1.9 × 10-4; interaction p = 3.5 × 10-2) of LTBR in the PFC provides further support for the role E2 plays in the brain, by demonstrating that the regulation of some genes that are altered by ovariectomy may also be modulated by Ov followed by hormone replacement therapy (HRT). These results present real opportunities to understand the specific biological mechanisms that are altered with depleted E2. Given E2's potential role in cognitive decline and neuroinflammation, our findings could lead to the discovery of novel therapeutics to slow cognitive decline. Together, this work represents a major step toward understanding molecular changes in the brain that are caused by ovariectomy and how E2 treatment may revert or protect against the negative neuro-related consequences caused by a depletion in estrogen as women approach menopause.

Complex regulatory networks influence pluripotent cell state transitions in human iPSCs

Stem cells exist in vitro in a spectrum of interconvertible pluripotent states. Analyzing hundreds of hiPSCs derived from different individuals, we show the proportions of these pluripotent states vary considerably across lines. We discover 13 gene network modules (GNMs) and 13 regulatory network modules (RNMs), which are highly correlated with each other suggesting that the coordinated co-accessibility of regulatory elements in the RNMs likely underlie the coordinated expression of genes in the GNMs. Epigenetic analyses reveal that regulatory networks underlying self-renewal and pluripotency are more complex than previously realized. Genetic analyses identify thousands of regulatory variants that overlapped predicted transcription factor binding sites and are associated with chromatin accessibility in the hiPSCs. We show that the master regulator of pluripotency, the NANOG-OCT4 Complex, and its associated network are significantly enriched for regulatory variants with large effects, suggesting that they play a role in the varying cellular proportions of pluripotency states between hiPSCs. Our work bins tens of thousands of regulatory elements in hiPSCs into discrete regulatory networks, shows that pluripotency and self-renewal processes have a surprising level of regulatory complexity, and suggests that genetic factors may contribute to cell state transitions in human iPSC lines.

Intrinsically dysregulated cellular stress signaling genes and gene networks in postpartum depression.

Postpartum depression (PPD) is a leading cause of morbidity and mortality among women. Clinically, the administration and withdrawal of supraphysiologic estradiol and progesterone (E2 + P) can cause affective symptom reoccurrence in women with a history of PPD, but not matched controls. To investigate the cellular basis underlying this differential affective response, lymphoblastoid cell lines (LCLs) were derived from women with and without past PPD and compared transcriptomically in hormone conditions mimicking pregnancy and parturition: supraphysiologic E2 + P-addback; supraphysiologic E2 + P-withdrawal; and no added E2 + P (Baseline). RNA-sequencing identified unique differentially expressed genes (DEGs) in all hormone conditions, but the majority tended to be downregulated in PPD and observed in E2 + P-addback. Two of these DEGs were evolutionarily conserved cellular stress regulators: IMPACT, an integrative response protein maintaining translational homeostasis, and WWTR1, a transcriptional coactivator in the 'Hippo' pathway mediating cell proliferation and survival. Correspondingly, significant gene network modules were linked to cell cycle progression, estrogen response, and immune dysregulation, suggesting innate differences in intracellular signaling in PPD. In certain hormone conditions, PPD LCLs displayed increased GATA3 expression (an upstream regulator of IMPACT and WWTR1) and differentially phosphorylated eiF2α (the ultimate downstream target of IMPACT). Taken together, these transcriptomic data primarily implicate innately dysregulated cellular responses as potentially influencing mood and/or escalating PPD risk. Furthermore, the intrinsic downregulation of IMPACT's translation and WWTR1's transcription networks may suggest a novel link between PPD and a compromised ability to maintain homeostasis in the context of cellular stress occurring during pregnancy and parturition.

Transcriptome and morphological analyses of double flower formation in Dianthus chinensis

The double flower developmental process is regulated via a complex transcriptional regulatory network. To understand this highly dynamic and complex developmental process of Dianthus spp., we performed a comparative analysis of floral morphology and transcriptome dynamics in simple flowers and double flowers. We found that the primordium of double flowers of ‘X’ was larger in size compared to that of simple flowers of ‘L’ in Dianthus chinensis. RNA-seq and Weighted Gene Co-expression Network Analysis (WGCNA) during flower development, identified stage-specific gene network modules. Expression analysis by RNA-seq indicated that a group of genes related to floral meristem identity, primordia position and polarity were highly expressed in double flowers genotypes compared to simple flowers genotypes, suggesting their roles in double-petal formation. A total of 21 DEGs related to petal number were identified between simple and double flowers. The experiments of in situ hybridization revealed that DcaAP2L, DcaLFY and DcaUFO genes were expressed in the intra-sepal boundary and petal boundary. We proposed a potential transcriptional regulatory network for simple and double flower development. This study provides novel insights into the molecular mechanism underlying double flower formation in Dianthus spp.

Identifying novel regulators of placental development using time-series transcriptome data.

The placenta serves as a connection between the mother and the fetus during pregnancy, providing the fetus with oxygen, nutrients, and growth hormones. However, the regulatory mechanisms and dynamic gene interaction networks underlying early placental development are understudied. Here, we generated RNA-sequencing data from mouse fetal placenta at embryonic days 7.5, 8.5, and 9.5 to identify genes with timepoint-specific expression, then inferred gene interaction networks to analyze highly connected network modules. We determined that timepoint-specific gene network modules were associated with distinct developmental processes, and with similar expression profiles to specific human placental cell populations. From each module, we identified hub genes and their direct neighboring genes, which were predicted to govern placental functions. We confirmed that four novel candidate regulators identified through our analyses regulate cell migration in the HTR-8/SVneo cell line. Overall, we predicted several novel regulators of placental development expressed in specific placental cell types using network analysis of bulk RNA-sequencing data. Our findings and analysis approaches will be valuable for future studies investigating the transcriptional landscape of early development.

WGCNA Analysis Revealed the Hub Genes Related to Soil Cadmium Stress in Maize Kernel (Zea mays L.).

Soil contamination by heavy metals has become a prevalent topic due to their widespread release from industry, agriculture, and other human activities. Great progress has been made in elucidating the uptake and translocation of cadmium (Cd) accumulation in rice. However, there is still little known about corresponding progress in maize. In the current study, we performed a comparative RNA-Seq-based approach to identify differentially expressed genes (DEGs) of maize immature kernel related to Cd stress. In total, 55, 92, 22, and 542 DEGs responsive to high cadmium concentration soil were identified between XNY22-CHS-8 vs. XNY22-YA-8, XNY22-CHS-24 vs. XNY22-YA-24, XNY27-CHS-8 vs. XNY27-YA-8, and XNY27-CHS-24 vs. XNY27-YA-24, respectively. The weighted gene co-expression network analysis (WGCNA) categorized the 9599 Cd stress-responsive hub genes into 37 different gene network modules. Combining the hub genes and DEGs, we obtained 71 candidate genes. Gene Ontology (GO) enrichment analysis of genes in the greenyellow module in XNY27-YA-24 and connectivity genes of these 71 candidate hub genes showed that the responses to metal ion, inorganic substance, abiotic stimulus, hydrogen peroxide, oxidative stress, stimulus, and other processes were enrichment. Moreover, five candidate genes that were responsive to Cd stress in maize kernel were detected. These results provided the putative key genes and pathways to response to Cd stress in maize kernel, and a useful dataset for unraveling the underlying mechanism of Cd accumulation in maize kernel.

Transcriptomic Analysis Reveals the Correlation between End-of-Day Far Red Light and Chilling Stress in Setaria viridis.

Low temperature and end-of-day far-red (EOD-FR) light signaling are two key factors limiting plant production and geographical location worldwide. However, the transcriptional dynamics of EOD-FR light conditions during chilling stress remain poorly understood. Here, we performed a comparative RNA-Seq-based approach to identify differentially expressed genes (DEGs) related to EOD-FR and chilling stress in Setaria viridis. A total of 7911, 324, and 13431 DEGs that responded to low temperature, EOD-FR and these two stresses were detected, respectively. Further DEGs analysis revealed that EOD-FR may enhance cold tolerance in plants by regulating the expression of genes related to cold tolerance. The result of weighted gene coexpression network analysis (WGCNA) using 13431 nonredundant DEGs exhibited 15 different gene network modules. Interestingly, a CO-like transcription factor named BBX2 was highly expressed under EOD-FR or chilling conditions. Furthermore, we could detect more expression levels when EOD-FR and chilling stress co-existed. Our dataset provides a valuable resource for the regulatory network involved in EOD-FR signaling and chilling tolerance in C4 plants.

Identification of Key Gene Network Modules and Hub Genes Associated with Wheat Response to Biotic Stress Using Combined Microarray Meta-analysis and WGCN Analysis.

Wheat (Triticum aestivum) is one of the major crops worldwide and a primary source of calories for human food. Biotic stresses such as fungi, bacteria, and diseases limit wheat production. Althoughplant breeding and genetic engineeringforbiotic stress resistance have been suggested as promising solutions to handle losses caused by biotic stress factors, a comprehensive understanding of molecular mechanisms and identifying key genes is a critical step to obtaining success. Here, a network-based meta-analysis approach based on two main statistical methods was used to identify key genes and molecular mechanisms of the wheat response to biotic stress. A total of 163 samples (21,792 genes) from 10 datasets were analyzed. Fisher Z test based on the p-value and REM method based on effect size resulted in 533 differentially expressed genes (p < 0.001 and FDR < 0.001). WGCNA analysis using a dynamic tree-cutting algorithm was used to construct a co-expression network and three significant modules were detected. The modules were significantly enriched by 16 BP terms and 4 KEGG pathways (Benjamini-Hochberg FDR < 0.001). A total of nine hub genes (a top 1.5% of genes with the highest degree) were identified from the constructed network. The identification of DE genes, gene-gene co-expressing network, and hub genes may contribute to uncovering the molecular mechanisms of the wheat response to biotic stress.

Within- and cross-tissue gene regulations were disrupted by PM2.5 nitrate exposure and associated with respiratory functions

BackgroundPathogenesis of complex diseases often involves multiple organs/tissue-types. To date, the PM2.5 exposure's toxic effects and induced disease risks were not studied at multi-tissue level. MethodsC57BL/6 mice (n = 40) were exposed to PM2.5 NO3− and clean air, respectively, and afterwards assessed respiratory functions and transcriptome in relevant tissues: blood and lung. We constructed within- and cross-tissue gene regulation networks and identified network modules associated with exposure and respiratory functions. ResultsPM2.5 NO3− exposure elevated naïve B cells proportion in blood (p = 0.0028). Among the 6000 highest expressed genes in blood, 18.8 % (1133 genes) were altered by exposure at p ≤ 0.05 level, among which 763 genes were also associated with respiratory function (enrichment folds = 7.63, p = 2.7E-189). The exposure disrupted blood genes were primarily in the immunoregulation pathways. Both within- and cross-tissue gene network modules were perturbed by exposure and associated with respiratory function. An immunodeficiency related cross-tissue module of 555 genes was affected by exposure (p = 0.0023) and strongly correlated with FEV0.05/FVC (r = 0.61 and p = 3E-5). ConclusionsThis study aims to fill in a major knowledge gap and investigated the effect of PM2.5 exposure simultaneously in multiple tissues. We provided novel evidence that PM2.5 NO3− exposure profoundly perturbed within- and cross-tissue gene regulations, and highlighted their roles in the etiology of respiratory decline.

FUNAGE-Pro: comprehensive web server for gene set enrichment analysis of prokaryotes.

Recent advances in the field of high throughput (meta-)transcriptomics and proteomics call for easy and rapid methods enabling to explore not only single genes or proteins but also extended biological systems. Gene set enrichment analysis is commonly used to find relations in a set of genes and helps to uncover the biological meaning in results derived from high-throughput data. The basis for gene set enrichment analysis is a solid functional classification of genes. Here, we describe a comprehensive database containing multiple functional classifications of genes of all (>55 000) publicly available complete bacterial genomes. In addition to the most common functional classes such as COG and GO, also KEGG, InterPro, PFAM, eggnog and operon classes are supported. As classification data for features is often not available, we offer fast annotation and classification of proteins in any newly sequenced bacterial genome. The web server FUNAGE-Pro enables fast functional analysis on single gene sets, multiple experiments, time series data, clusters, and gene network modules for any prokaryote species or strain. FUNAGE-Pro is freely available at http://funagepro.molgenrug.nl.

Transcriptional Activation, Deactivation and Rebound Patterns in Cortex, Hippocampus and Amygdala in Response to Ketamine Infusion in Rats.

Ketamine, an N-methyl-D-aspartate (NMDA)-receptor antagonist, is a recently revitalized treatment for pain and depression, yet its actions at the molecular level remain incompletely defined. In this molecular-pharmacological investigation in the rat, we used short- and longer-term infusions of high dose ketamine to stimulate neuronal transcription processes. We hypothesized that a progressively stronger modulation of neuronal gene networks would occur over time in cortical and limbic pathways. A continuous intravenous administration paradigm for ketamine was developed in rat consisting of short (1 h) and long duration (10 h, and 10 h + 24 h recovery) infusions of anesthetic concentrations to activate or inhibit gene transcription in a pharmacokinetically controlled fashion. Transcription was measured by RNA-Seq in three brain regions: frontal cortex, hippocampus, and amygdala. Cellular level gene localization was performed with multiplex fluorescent in situ hybridization. Induction of a shared transcriptional regulatory network occurred within 1 h in all three brain regions consisting of (a) genes involved in stimulus-transcription factor coupling that are induced during altered synaptic activity (immediate early genes, IEGs, such as c-Fos, 9–12 significant genes per brain region, p < 0.01 per gene) and (b) the Nrf2 oxidative stress-antioxidant response pathway downstream from glutamate signaling (Nuclear Factor Erythroid-Derived 2-Like 2) containing 12–25 increasing genes (p < 0.01) per brain region. By 10 h of infusion, the acute results were further reinforced and consisted of more and stronger gene alterations reflecting a sustained and accentuated ketamine modulation of regional excitation and plasticity. At the cellular level, in situ hybridization localized up-regulation of the plasticity-associated gene Bdnf, and the transcription factors Nr4a1 and Fos, in cortical layers III and V. After 24 h recovery, we observed overshoot of transcriptional processes rather than a smooth return to homeostasis suggesting an oscillation of plasticity occurs during the transition to a new phase of neuronal regulation. These data elucidate critical molecular regulatory actions during and downstream of ketamine administration that may contribute to the unique drug actions of this anesthetic agent. These molecular investigations point to pathways linked to therapeutically useful attributes of ketamine.

Key gene network related to primary ciliary dyskinesia in hippocampus of patients with Alzheimer’s disease revealed by weighted gene co-expression network analysis

BackgroundAlzheimer’s disease (AD) is closely related to aging, showing an increasing incidence rate for years. As one of the main brain regions involved in AD, hippocampus has been extensively studied due to its association with many human diseases. However, little is known about its association with primary ciliary dyskinesia (PCD).Material and MethodsThe microarray data of hippocampus on AD were retrieved from the Gene Expression Omnibus (GEO) database to construct the co-expression network by weighted gene co-expression network analysis (WGCNA). The gene network modules associated with AD screened with the common genes were further annotated based on Gene Ontology (GO) database and enriched based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. The protein-protein interaction (PPI) network was constructed based on STRING database to identify the hub genes in the network.ResultsGenes involved in PCD were identified in the hippocampus of AD patients. Functional analysis revealed that these genes were mainly enriched in ciliary tissue, ciliary assembly, axoneme assembly, ciliary movement, microtubule based process, microtubule based movement, organelle assembly, axoneme dynamin complex, cell projection tissue, and microtubule cytoskeleton tissue. A total of 20 central genes, e.g., DYNLRB2, ZMYND10, DRC1, DNAH5, WDR16, TTC25, and ARMC4 were identified as hub genes related to PCD in hippocampus of AD patients.ConclusionOur study demonstrated that AD and PCD have common metabolic pathways. These common pathways provide novel evidence for further investigation of the pathophysiological mechanism and the hub genes suggest new therapeutic targets for the diagnosis and treatment of AD and PCD.SubjectsBioinformatics, Cell Biology, Molecular Biology, Neurology.

Transcriptome analysis of Aedes albopictus midguts infected by dengue virus identifies a gene network module highly associated with temperature

BackgroundDengue is prevalent worldwide and is transmitted by Aedes mosquitoes. Temperature is a strong driver of dengue transmission. However, little is known about the underlying mechanisms.MethodsAedes albopictus mosquitoes exposed or not exposed to dengue virus serotype 2 (DENV-2) were reared at 23 °C, 28 °C and 32 °C, and midguts and residual tissues were evaluated at 7 days after infection. RNA sequencing of midgut pools from the control group, midgut breakthrough group and midgut nonbreakthrough group at different temperatures was performed. The transcriptomic profiles were analyzed using the R package, followed by weighted gene correlation network analysis (WGCNA) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis to identify the important molecular mechanisms regulated by temperature.ResultsThe midgut infection rate and midgut breakthrough rate at 28 °C and 32 °C were significantly higher than those at 23 °C, which indicates that high temperature facilitates DENV-2 breakthrough in the Ae. albopictus midgut. Transcriptome sequencing was performed to investigate the antiviral mechanism in the midgut. The midgut gene expression datasets clustered with respect to temperature, blood-feeding and midgut breakthrough. Over 1500 differentially expressed genes were identified by pairwise comparisons of midguts at different temperatures. To assess key molecules regulated by temperature, we used WGCNA, which identified 28 modules of coexpressed genes; the ME3 module correlated with temperature. KEGG analysis indicated that RNA degradation, Toll and immunodeficiency factor signaling and other pathways are regulated by temperature.ConclusionsTemperature affects the infection and breakthrough of Ae. albopictus midguts invaded by DENV-2, and Ae. albopictus midgut transcriptomes change with temperature. The candidate genes and key pathways regulated by temperature provide targets for the prevention and control of dengue.Graphical abstract

Abstract WP241: Transcriptome Analysis Of The Liver Reveals Novel Mechanisms By Which Prenatal Alcohol Exposure Contributes To Worse Stoke Outcome

Introduction: Susceptibility to adult-onset disease has been linked to early life adversity. In this study, we assessed the impact of a common early adverse experience, prenatal alcohol exposure (PAE), on neurobehavioral outcomes following cerebrovascular ischemic stroke. We also assessed impact of PAE and stroke on the transcriptome of intermediary organs of the gut-brain axis. Specifically, mesenteric adipose and liver, which are supplied by portal circulation from the gut are hypothesized to mediate inflammatory and neuroprotective adaptations to stroke. Methods: Pregnant Sprague Dawley rats were exposed to either vaporized ethanol or room air from gestation day 8 to 19. At 5 months, progeny were subjected to unilateral endothelin-1 induced occlusion of the middle cerebral artery and outcomes evaluated after 2 days. Stroke induced disabilities were assessed from behavioral assays (adhesive removal, Vibrissae evoked fore-limb placement and circling) and infarct size. The liver transcriptome was assessed by sequencing, from age-matched stroke-exposed and naïve offspring. Results: PAE rats had worse neurological scores 2 days post stroke compared to control (p&lt;0.05). We observed that the liver was substantially more transcriptionally responsive to PAE status than adipose tissue. Moreover, PAE offspring exhibited a 1.6-fold increase in upregulated genes and a 2.4-fold increase in down-regulated liver gene transcripts following stroke, compared to control offspring. Only 448 differentially expressed genes (DEGs) were common between the control and PAE groups. DEGs unique to PAE contribute to triglyceride metabolism, steroid metabolism, and cholesterol biosynthesis pathways. Weighted correlation network analysis identified 24 gene network modules, one of which, a liver proinflammatory network, was strongly correlated with post stroke neurological function (Pearson's r=-0.63, p&lt;0.01). Conclusion: The present study shows that PAE predisposes middle-aged rat offspring for worse outcomes following cerebrovascular ischemic stroke. This effect is associated with persistent reprogramming of the liver transcriptome and an exacerbated proinflammatory response.

Transcriptomic and physiological analysis identifies a gene network module highly associated with brassinosteroid regulation in hybrid sweetgum tissues differing in the capability of somatic embryogenesis.

Somatic embryogenesis is a preferred method for large-scale production of forest trees due to its high propagation efficiency. In this study, hybrid sweetgum leaves with phase changes from mature to embryogenic state were selected as experimental material to study somatic embryo initiation. Embryogenicity ranged from high to low, i.e. from 45%, 25%, and 12.5% to 0, with the samples of embryogenic callus (EC), whiten leaf edge (WLI), whiten leaf (WLII), and green leaf (GL) respectively. High correlations existed between embryogenicity and endogenous brassinosteroids (BRs) (r = 0.95, p < 0.05). Similarly, concentrations of endogenous BRs of the sample set correlated positively (r = 0.93, 0.99, 0.87, 0.99, 0.96 respectively, P < 0.05) to expression of somatic embryo (SE)-related genes, i.e. BBM, LEC2, ABI3, PLT2, and WOX2. Hierarchical cluster and weighted gene coexpression network analysis identified modules of coexpressed genes and network in 4820 differentially expressed genes (DEGs) from All-BR-Regulated Genes (ABRG). Moreover, exogenously-supplemented epiBR, together with 2,4-D and 6-BA, increased embryogenicity of GL-sourced callus, and expression of SE- and auxin-related genes, while brassinazole (BRZ), a BR biosynthesis inhibitor, reduced embryogenicity. Evidences obtained in this study revealed that BRs involved in phase change of leaf explants and may function in regulating gene expression and enhancing auxin effects. This study successfully established protocols for inducing somatic embryogenesis from leaf explants in hybrid sweetgum, which could facilitate the propagation process greatly, and provide theoretical basis for manipulating SE competence of explants in ornamental woody plants.

Toxicity Analysis of Pentachlorophenol Data with a Bioinformatics Tool Set.

Rapid progress in technologies opened the new era of computer-leaded analytics, leaving humans more space for experimental design and decision making. Here we demonstrate the machine learning analysis workflow represented by spectral clustering, elucidation of evolutionary conserved transcription regulation, and network analysis using reverse engineering. Analysis of genes induced by the Pentachlorophenol toxic chemical revealed two subnetworks, one orchestrated by Interferon and another by Nuclear receptor factor 2 (NRF2) gene. Furthermore, network-inference based analysis identified a gene network module composed of genes associated with interferon signaling and their regulatory interaction with downstream genes, especially TRIM family proteins involved in responses of innate immune systems.

Placental gene network modules are associated with maternal stress during pregnancy and infant temperament.

Maternal psychosocial stress during pregnancy (MPSP) is a known contributor to maladaptive neurobehavioral development of the offspring; however, the underlying molecular mechanisms linking MPSP with childhood outcome remain largely unknown. Transcriptome-wide gene expression data were generated using RNA-seq from placenta samples collected in a multi-ethnic urban birth cohort in New York City (n=129). Weighted gene co-expression network analysis (WGCNA) was used to characterize placental co-expression modules, which were then evaluated for their associations with MPSP and infant temperament. WGCNA revealed 16 gene coexpression modules. One module, enriched for regulation of chromosome organization/gene expression, was positively associated with MPSP and negatively associated with Regulatory Capacity (REG), a component of infant temperament. Two other modules, enriched for cotranslational protein targeting and cell cycle regulation, respectively, displayed negative associations with MPSP and positive associations with REG. A module enriched with oxidative phosphorylation/mitochondrial translation was positively associated with REG. These findings support the notion that the placenta provides a functional in utero link between MPSP and infant temperament, possibly through transcriptional regulation of placental gene expression.