Using VisANT to Analyze Networks

Characterization of Accessible Chromatin Regions in Cattle Rumen Epithelial Tissue during Weaning

Weaning in ruminants is characterized by the transition from a milk-based diet to a solid diet, which drives a critical gastrointestinal tract transformation. Understanding the regulatory control of this transformation during weaning can help to identify strategies to improve rumen health. This study aimed to identify regions of accessible chromatin in rumen epithelial tissue in pre- and post-weaning calves and investigate differentially accessible regions (DARs) to uncover regulatory elements in cattle rumen development using the ATAC-seq approach. A total of 126,071 peaks were identified, covering 1.15% of the cattle genome. From these accessible regions, 2766 DARs were discovered. Gene ontology enrichment resulted in GO terms related to the cell adhesion, anchoring junction, growth, cell migration, motility, and morphogenesis. In addition, putative regulatory canonical pathways were identified (TGFβ, integrin-linked kinase, integrin signaling, and regulation of the epithelial–mesenchymal transition). Canonical pathways integrated with co-expression results showed that TGFβ and ILK signaling pathways play essential roles in rumen development through the regulation of cellular adhesions. In this study, DARs during weaning were identified, revealing enhancers, transcription factors, and candidate target genes that represent potential biomarkers for the bovine rumen development, which will serve as a molecular tool for rumen development studies.

GWENA: gene co-expression networks analysis and extended modules characterization in a single Bioconductor package

BACKGROUND

Network-based analysis of gene expression through co-expression networks can be used to investigate modular relationships occurring between genes performing different biological functions. An extended description of each of the network modules is therefore a critical step to understand the underlying processes contributing to a disease or a phenotype. Biological integration, topology study and conditions comparison (e.g. wild vs mutant) are the main methods to do so, but to date no tool combines them all into a single pipeline.

RESULTS

Here we present GWENA, a new R package that integrates gene co-expression network construction and whole characterization of the detected modules through gene set enrichment, phenotypic association, hub genes detection, topological metric computation, and differential co-expression. To demonstrate its performance, we applied GWENA on two skeletal muscle datasets from young and old patients of GTEx study. Remarkably, we prioritized a gene whose involvement was unknown in the muscle development and growth. Moreover, new insights on the variations in patterns of co-expression were identified. The known phenomena of connectivity loss associated with aging was found coupled to a global reorganization of the relationships leading to expression of known aging related functions.

CONCLUSION

GWENA is an R package available through Bioconductor ( https://bioconductor.org/packages/release/bioc/html/GWENA.html ) that has been developed to perform extended analysis of gene co-expression networks. Thanks to biological and topological information as well as differential co-expression, the package helps to dissect the role of genes relationships in diseases conditions or targeted phenotypes. GWENA goes beyond existing packages that perform co-expression analysis by including new tools to fully characterize modules, such as differential co-expression, additional enrichment databases, and network visualization.

A novel LncRNA-based prognostic score reveals TP53-dependent subtype of lung adenocarcinoma with poor survival

The prognostic signatures play an essential role in the era of personalised therapy for cancer patients including lung adenocarcinoma (LUAD). Long noncoding RNA (LncRNA), a relatively novel class of RNA, has shown to play a crucial role in all the areas of cancer biology. Here, we developed and validated a robust LncRNA-based prognostic signature for LUAD patients using three different cohorts. In the discovery cohort, four LncRNAs were identified with 10% false discovery rate and a hazard ratio of >10 using univariate Cox regression analysis. A risk score, generated from the four LncRNAs' expression, was found to be a significant predictor of survival in the discovery and validation cohort (p = 9.97 × 10 ^-8 and 1.41 × 10 ^-3 , respectively). Further optimisation of four LncRNAs signature in the validation cohort, generated a three LncRNAs prognostic score (LPS), which was found to be an independent predictor of survival in both the cohorts ( p = 1.00 × 10 ^-6 and 7.27 × 10 ^-4 , respectively). The LPS also significantly divided survival in clinically important subsets, including Stage I ( p = 9.00 × 10 ^-4 and 4.40 × 10 ^-2 , respectively), KRAS wild-type (WT), KRAS mutant ( p = 4.00 × 10 ^-3 and 4.30 × 10 ^-2 , respectively) and EGFR WT ( p = 2.00 × 10 ^-4 ). In multivariate analysis LPS outperformed, eight previous prognosticators. Further, individual members of LPS showed a significant correlation with survival in microarray data sets. Mutation analysis showed that high-LPS patients have a higher mutation rate and inactivation of the TP53 pathway. In summary, we identified and validated a novel LncRNA signature LPS for LUAD.

Whole-transcriptome analysis delineates the human placenta gene network and its associations with fetal growth

BACKGROUND

The placenta is the principal organ regulating intrauterine growth and development, performing critical functions on behalf of the developing fetus. The delineation of functional networks and pathways driving placental processes has the potential to provide key insight into intrauterine perturbations that result in adverse birth as well as later life health outcomes.

RESULTS

We generated the transcriptome-wide profile of 200 term human placenta using the Illumina HiSeq 2500 platform and characterized the functional placental gene network using weighted gene coexpression network analysis (WGCNA). We identified 17 placental coexpression network modules that were dominated by functional processes including growth, organ development, gas exchange and immune response. Five network modules, enriched for processes including cellular respiration, amino acid transport, hormone signaling, histone modifications and gene expression, were associated with birth weight; hub genes of all five modules (CREB3, DDX3X, DNAJC14, GRHL1 and C21orf91) were significantly associated with fetal growth restriction, and one hub gene (CREB3) was additionally associated with fetal overgrowth.

CONCLUSIONS

In this largest RNA-Seq based transcriptome-wide profiling study of human term placenta conducted to date, we delineated a placental gene network with functional relevance to fetal growth using a network-based approach with superior scale reduction capacity. Our study findings not only implicate potential molecular mechanisms underlying fetal growth but also provide a reference placenta gene network to inform future studies investigating placental dysfunction as a route to future disease endpoints.

Intra- and interregional coregulation of opioid genes: broken symmetry in spinal circuits

Regulation of the formation and rewiring of neural circuits by neuropeptides may require coordinated production of these signaling molecules and their receptors that may be established at the transcriptional level. Here, we address this hypothesis by comparing absolute expression levels of opioid peptides with their receptors, the largest neuropeptide family, and by characterizing coexpression (transcriptionally coordinated) patterns of these genes. We demonstrated that expression patterns of opioid genes highly correlate within and across functionally and anatomically different areas. Opioid peptide genes, compared with their receptor genes, are transcribed at much greater absolute levels, which suggests formation of a neuropeptide cloud that covers the receptor-expressed circuits. Surprisingly, we found that both expression levels and the proportion of opioid receptors are strongly lateralized in the spinal cord, interregional coexpression patterns are side specific, and intraregional coexpression profiles are affected differently by left- and right-side unilateral body injury. We propose that opioid genes are regulated as interconnected components of the same molecular system distributed between distinct anatomic regions. The striking feature of this system is its asymmetric coexpression patterns, which suggest side-specific regulation of selective neural circuits by opioid neurohormones.-Kononenko, O., Galatenko, V., Andersson, M., Bazov, I., Watanabe, H., Zhou, X. W., Iatsyshyna, A., Mityakina, I., Yakovleva, T., Sarkisyan, D., Ponomarev, I., Krishtal, O., Marklund, N., Tonevitsky, A., Adkins, D. L., Bakalkin, G. Intra- and interregional coregulation of opioid genes: broken symmetry in spinal circuits.

Integrative Analysis of Genomics and Transcriptome Data to Identify Potential Functional Genes of BMDs in Females

Osteoporosis is known to be highly heritable. However, to date, the findings from more than 20 genome-wide association studies (GWASs) have explained less than 6% of genetic risks. Studies suggest that the missing heritability data may be because of joint effects among genes. To identify novel heritability for osteoporosis, we performed a system-level study on bone mineral density (BMD) by weighted gene coexpression network analysis (WGCNA), using the largest GWAS data set for BMD in the field, Genetic Factors for Osteoporosis Consortium (GEFOS-2), and a transcriptomic gene expression data set generated from transiliac bone biopsies in women. A weighted gene coexpression network was generated for 1574 genes with GWAS nominal evidence of association (p ≤ 0.05) based on dissimilarity measurement on the expression data. Twelve distinct gene modules were identified, and four modules showed nominally significant associations with BMD (p ≤ 0.05), but only one module, the yellow module, demonstrated a good correlation between module membership (MM) and gene significance (GS), suggesting that the yellow module serves an important biological role in bone regulation. Interestingly, through characterization of module content and topology, the yellow module was found to be significantly enriched with contractile fiber part (GO:044449), which is widely recognized as having a close relationship between muscle and bone. Furthermore, detailed submodule analyses of important candidate genes (HOMER1, SPTBN1) by all edges within the yellow module implied significant enrichment of functional connections between bone and cytoskeletal protein binding. Our study yielded novel information from system genetics analyses of GWAS data jointly with transcriptomic data. The findings highlighted a module and several genes in the model as playing important roles in the regulation of bone mass in females, which may yield novel insights into the genetic basis of osteoporosis. © 2016 American Society for Bone and Mineral Research.

Visualization of Metabolic Interaction Networks in Microbial Communities Using VisANT 5.0

The complexity of metabolic networks in microbial communities poses an unresolved visualization and interpretation challenge. We address this challenge in the newly expanded version of a software tool for the analysis of biological networks, VisANT 5.0. We focus in particular on facilitating the visual exploration of metabolic interaction between microbes in a community, e.g. as predicted by COMETS (Computation of Microbial Ecosystems in Time and Space), a dynamic stoichiometric modeling framework. Using VisANT's unique metagraph implementation, we show how one can use VisANT 5.0 to explore different time-dependent ecosystem-level metabolic networks. In particular, we analyze the metabolic interaction network between two bacteria previously shown to display an obligate cross-feeding interdependency. In addition, we illustrate how a putative minimal gut microbiome community could be represented in our framework, making it possible to highlight interactions across multiple coexisting species. We envisage that the "symbiotic layout" of VisANT can be employed as a general tool for the analysis of metabolism in complex microbial communities as well as heterogeneous human tissues. VisANT is freely available at: http://visant.bu.edu and COMETS at http://comets.bu.edu.

Frontotemporal dementia: insights into the biological underpinnings of disease through gene co-expression network analysis

Background

In frontotemporal dementia (FTD) there is a critical lack in the understanding of biological and molecular mechanisms involved in disease pathogenesis. The heterogeneous genetic features associated with FTD suggest that multiple disease-mechanisms are likely to contribute to the development of this neurodegenerative condition.

We here present a systems biology approach with the scope of i) shedding light on the biological processes potentially implicated in the pathogenesis of FTD and ii) identifying novel potential risk factors for FTD. We performed a gene co-expression network analysis of microarray expression data from 101 individuals without neurodegenerative diseases to explore regional-specific co-expression patterns in the frontal and temporal cortices for 12 genes (MAPT, GRN, CHMP2B, CTSC, HLA-DRA, TMEM106B, C9orf72, VCP, UBQLN2, OPTN, TARDBP and FUS) associated with FTD and we then carried out gene set enrichment and pathway analyses, and investigated known protein-protein interactors (PPIs) of FTD-genes products.

Results

Gene co-expression networks revealed that several FTD-genes (such as MAPT and GRN, CTSC and HLA-DRA, TMEM106B, and C9orf72, VCP, UBQLN2 and OPTN) were clustering in modules of relevance in the frontal and temporal cortices. Functional annotation and pathway analyses of such modules indicated enrichment for: i) DNA metabolism, i.e. transcription regulation, DNA protection and chromatin remodelling (MAPT and GRN modules); ii) immune and lysosomal processes (CTSC and HLA-DRA modules), and; iii) protein meta/catabolism (C9orf72, VCP, UBQLN2 and OPTN, and TMEM106B modules). PPI analysis supported the results of the functional annotation and pathway analyses.

Conclusions

This work further characterizes known FTD-genes and elaborates on their biological relevance to disease: not only do we indicate likely impacted regional-specific biological processes driven by FTD-genes containing modules, but also do we suggest novel potential risk factors among the FTD-genes interactors as targets for further mechanistic characterization in hypothesis driven cell biology work.

Electronic supplementary material

The online version of this article (doi:10.1186/s13024-016-0085-4) contains supplementary material, which is available to authorized users.