Pan- and core- gene association networks: Integrative approaches to understanding biological regulation

Co-expression pan-network reveals genes involved in complex traits within maize pan-genome

BACKGROUND

With the advances in the high throughput next generation sequencing technologies, genome-wide association studies (GWAS) have identified a large set of variants associated with complex phenotypic traits at a very fine scale. Despite the progress in GWAS, identification of genotype-phenotype relationship remains challenging in maize due to its nature with dozens of variants controlling the same trait. As the causal variations results in the change in expression, gene expression analyses carry a pivotal role in unraveling the transcriptional regulatory mechanisms behind the phenotypes.

RESULTS

To address these challenges, we incorporated the gene expression and GWAS-driven traits to extend the knowledge of genotype-phenotype relationships and transcriptional regulatory mechanisms behind the phenotypes. We constructed a large collection of gene co-expression networks and identified more than 2 million co-expressing gene pairs in the GWAS-driven pan-network which contains all the gene-pairs in individual genomes of the nested association mapping (NAM) population. We defined four sub-categories for the pan-network: (1) core-network contains the highest represented ~ 1% of the gene-pairs, (2) near-core network contains the next highest represented 1-5% of the gene-pairs, (3) private-network contains ~ 50% of the gene pairs that are unique to individual genomes, and (4) the dispensable-network contains the remaining 50-95% of the gene-pairs in the maize pan-genome. Strikingly, the private-network contained almost all the genes in the pan-network but lacked half of the interactions. We performed gene ontology (GO) enrichment analysis for the pan-, core-, and private- networks and compared the contributions of variants overlapping with genes and promoters to the GWAS-driven pan-network.

CONCLUSIONS

Gene co-expression networks revealed meaningful information about groups of co-regulated genes that play a central role in regulatory processes. Pan-network approach enabled us to visualize the global view of the gene regulatory network for the studied system that could not be well inferred by the core-network alone.

The spatial transcriptomic landscape of non-small cell lung cancer brain metastasis

Brain metastases (BrMs) are a common occurrence in lung cancer with a dismal outcome. To understand the mechanism of metastasis to inform prognosis and treatment, here we analyze primary and metastasized tumor specimens from 44 non-small cell lung cancer patients by spatial RNA sequencing, affording a whole transcriptome map of metastasis resolved with morphological markers for the tumor core, tumor immune microenvironment (TIME), and tumor brain microenvironment (TBME). Our data indicate that the tumor microenvironment (TME) in the brain, including the TIME and TBME, undergoes extensive remodeling to create an immunosuppressive and fibrogenic niche for the BrMs. Specifically, the brain TME is characterized with reduced antigen presentation and B/T cell function, increased neutrophils and M2-type macrophages, immature microglia, and reactive astrocytes. Differential gene expression and network analysis identify fibrosis and immune regulation as the major functional modules disrupted in both the lung and brain TME. Besides providing systems-level insights into the mechanism of lung cancer brain metastasis, our study uncovers potential prognostic biomarkers and suggests that therapeutic strategies should be tailored to the immune and fibrosis status of the BrMs.

Brain metastases (BrMs) in non-small cell lung cancer (NSCLC) are associated with dismal outcomes, and are possibly sustained by the brain microenvironment. Here, the authors analyse NSCLC BrMs using Digital Spatial Profiling and reveal fibrosis, immune suppression, and cell reprogramming in the BrM microenvironment.

Comparative transcriptome analyses between cultivated and wild grapes reveal conservation of expressed genes but extensive rewiring of co-expression networks

The transcriptomes of wild and cultivated grapes consists of similar expressed genes but distinct wiring of co-expressed modules associated with environmental conditions. Grapevine is an important fruit crop worldwide, with high economic value and widespread distribution. Commercial production is based on Vitis vinifera, and, to a lesser extent, on hybrids with American grapes, such as V. labrusca. Wild grape relatives are important sources of resistance against biotic and abiotic factors; however, their global gene expression patterns remain poorly characterized. We associated genome-wide transcript profiling to phenotypic analyses to investigate the responses of cultivated and wild vines to vineyard conditions. The expressed genes in the Vitis reference transcriptome are largely shared by wild grapes, V. labrusca hybrids and vinifera cultivars. In contrast, significant differential regulation between wild and vinifera genotypes represents 80% of gene expression variation, regardless of the environment. In wild grapes, genes associated to regulatory processes are downregulated, whereas those involved in metabolic pathways are upregulated, in comparison to vinifera. Photosynthesis-related ontologies are overrepresented in the induced genes, in agreement with higher contents of chlorophyll in wild grapes. Co-regulated gene network analyses provide evidence of more complex transcriptome organization in vinifera. In wild grapes, genes involved in signaling pathways of stress-related hormones are overrepresented in modules associated with the environment. Consensus network analyses revealed high preservation within co-regulated gene modules between cultivated and wild grapes, but divergent relationships among the expression clusters. In conclusion, the distinct phenotypes of wild and cultivated grapes are underlain by differences in gene expression, but also by distinct higher-order organization of the transcriptome and contrasting association of co-expressed gene clusters with the environment.

Topological clustering of regulatory genes confers pathogenic tolerance to cassava brown streak virus (CBSV) in cassava

Robustness, a naïve property of biological systems, enables organisms to maintain functions during perturbation and is crucial for improving the resilience of crops to prevailing stress conditions and diseases, guaranteeing food security. Most studies of robustness in crops have focused on genetic superiority based upon individual genes, overlooking the collaborative actions of multiple responsive genes and the regulatory network topology. This research aims to uncover patterns of gene cooperation leading to organismal robustness by studying the topology of gene co-expression networks (GCNs) of both CBSV virus resistant and susceptible cassava cultivars. The resulting GCNs show higher topological clustering of cooperative genes in the resistant cultivar, suggesting that the network architecture is central to attaining robustness. Despite a reduction in the number of hub genes in the resistant cultivar following the perturbation, essential biological functions contained in the network were maintained through neighboring genes that withstood the shock. The susceptible cultivar seemingly coped by inducing more gene actions in the network but could not maintain the functions required for plant growth. These findings underscore the importance of regulatory network architecture in ensuring phenotypic robustness and deepen our understanding of transcriptional regulation.

Genome-wide analysis and prediction of genes involved in the biosynthesis of polysaccharides and bioactive secondary metabolites in high-temperature-tolerant wild Flammulina filiformis

Background

Flammulina filiformis (previously known as Asian F. velutipes) is a popular commercial edible mushroom. Many bioactive compounds with medicinal effects, such as polysaccharides and sesquiterpenoids, have been isolated and identified from F. filiformis, but their biosynthesis and regulation at the molecular level remains unclear. In this study, we sequenced the genome of the wild strain F. filiformis Liu355, predicted its biosynthetic gene clusters (BGCs) and profiled the expression of these genes in wild and cultivar strains and in different developmental stages of the wild F. filiformis strain by a comparative transcriptomic analysis.

Results

We found that the genome of the F. filiformis was 35.01 Mb in length and harbored 10,396 gene models. Thirteen putative terpenoid gene clusters were predicted and 12 sesquiterpene synthase genes belonging to four different groups and two type I polyketide synthase gene clusters were identified in the F. filiformis genome. The number of genes related to terpenoid biosynthesis was higher in the wild strain (119 genes) than in the cultivar strain (81 genes). Most terpenoid biosynthesis genes were upregulated in the primordium and fruiting body of the wild strain, while the polyketide synthase genes were generally upregulated in the mycelium of the wild strain. Moreover, genes encoding UDP-glucose pyrophosphorylase and UDP-glucose dehydrogenase, which are involved in polysaccharide biosynthesis, had relatively high transcript levels both in the mycelium and fruiting body of the wild F. filiformis strain.

Conclusions

F. filiformis is enriched in a number of gene clusters involved in the biosynthesis of polysaccharides and terpenoid bioactive compounds and these genes usually display differential expression between wild and cultivar strains, even in different developmental stages. This study expands our knowledge of the biology of F. filiformis and provides valuable data for elucidating the regulation of secondary metabolites in this unique F. filiformis strain.