Modeling Genome-Wide Dynamic Regulatory Network in Mouse Lungs with Influenza Infection Using High-Dimensional Ordinary Differential Equations

The immune response to viral infection is regulated by an intricate network of many genes and their products. The reverse engineering of gene regulatory networks (GRNs) using mathematical models from time course gene expression data collected after influenza infection is key to our understanding of the mechanisms involved in controlling influenza infection within a host. A five-step pipeline: detection of temporally differentially expressed genes, clustering genes into co-expressed modules, identification of network structure, parameter estimate refinement, and functional enrichment analysis, is developed for reconstructing high-dimensional dynamic GRNs from genome-wide time course gene expression data. Applying the pipeline to the time course gene expression data from influenza-infected mouse lungs, we have identified 20 distinct temporal expression patterns in the differentially expressed genes and constructed a module-based dynamic network using a linear ODE model. Both intra-module and inter-module annotations and regulatory relationships of our inferred network show some interesting findings and are highly consistent with existing knowledge about the immune response in mice after influenza infection. The proposed method is a computationally efficient, data-driven pipeline bridging experimental data, mathematical modeling, and statistical analysis. The application to the influenza infection data elucidates the potentials of our pipeline in providing valuable insights into systematic modeling of complicated biological processes.     

Modulation of the ComA-dependent quorum response in Bacillus subtilis by multiple Rap proteins and Phr peptides

In Bacillus subtilis, extracellular peptide signaling regulates several biological processes. Secreted Phr signaling peptides are imported into the cell and act intracellularly to antagonize the activity of regulators known as Rap proteins. B. subtilis encodes several Rap proteins and Phr peptides, and the processes regulated by many of these Rap proteins and Phr peptides are unknown. We used DNA microarrays to characterize the roles that several rap-phr signaling modules play in regulating gene expression. We found that rapK-phrK regulates the expression of a number of genes activated by the response regulator ComA. ComA activates expression of genes involved in competence development and the production of several secreted products. Two Phr peptides, PhrC and PhrF, were previously known to stimulate the activity of ComA. We assayed the roles that PhrC, PhrF, and PhrK play in regulating gene expression and found that these three peptides stimulate ComA-dependent gene expression to different levels and are all required for full expression of genes activated by ComA. The involvement of multiple Rap proteins and Phr peptides allows multiple physiological cues to be integrated into a regulatory network that modulates the timing and magnitude of the ComA response.     

Viruses associated with human cancer

It is estimated that viral infections contribute to 15-20% of all human cancers. As obligatory intracellular parasites, viruses encode proteins that reprogram host cellular signaling pathways that control proliferation, differentiation, cell death, genomic integrity, and recognition by the immune system. These cellular processes are governed by complex and redundant regulatory networks and are surveyed by sentinel mechanisms that ensure that aberrant cells are removed from the proliferative pool. Given that the genome size of a virus is highly restricted to ensure packaging within an infectious structure, viruses must target cellular regulatory nodes with limited redundancy and need to inactivate surveillance mechanisms that would normally recognize and extinguish such abnormal cells. In many cases, key proteins in these same regulatory networks are subject to mutation in non-virally associated diseases and cancers. Oncogenic viruses have thus served as important experimental models to identify and molecularly investigate such cellular networks. These include the discovery of oncogenes and tumor suppressors, identification of regulatory networks that are critical for maintenance of genomic integrity, and processes that govern immune surveillance.     

Signature patterns revealed by microarray analyses of mice infected with influenza virus A and Streptococcus pneumoniae

We used cDNA microarrays to identify differentially expressed genes in mice in response to infections with influenza virus A/PR/8/34 (H1N1) and Streptococcus pneumoniae. Expression microarray analysis showed up-regulation and down-regulation of many genes involved in the defense, inflammatory response and intracellular signaling pathways including chemokine, apoptosis, MAPK, Notch, Jak-STAT, T-cell receptor and complement and coagulation cascades. We have revealed signature patterns of gene expression in mice infected with two different classes of pathogens: influenza virus A and S. pneumoniae. Quantitative real-time RT-PCR results confirmed microarray results for most of the genes tested. These studies document clear differences in gene expression profiles between mice infected with influenza virus A and S. pneumoniae. Identification of genes that are differentially expressed after respiratory infections can provide insights into the mechanisms by which the host interacts with different pathogens, useful information about stage of diseases and selection of suitable targets for early diagnosis and treatments. The advantage of this novel approach is that the detection of pathogens is based on the differences in host gene expression profiles in response to different pathogens instead of detecting pathogens directly.     

A densely interconnected genome-wide network of microRNAs and oncogenic pathways revealed using gene expression signatures

MicroRNAs (miRNAs) are important components of cellular signaling pathways, acting either as pathway regulators or pathway targets. Currently, only a limited number of miRNAs have been functionally linked to specific signaling pathways. Here, we explored if gene expression signatures could be used to represent miRNA activities and integrated with genomic signatures of oncogenic pathway activity to identify connections between miRNAs and oncogenic pathways on a high-throughput, genome-wide scale. Mapping >300 gene expression signatures to >700 primary tumor profiles, we constructed a genome-wide miRNA-pathway network predicting the associations of 276 human miRNAs to 26 oncogenic pathways. The miRNA-pathway network confirmed a host of previously reported miRNA/pathway associations and uncovered several novel associations that were subsequently experimentally validated. Globally, the miRNA-pathway network demonstrates a small-world, but not scale-free, organization characterized by multiple distinct, tightly knit modules each exhibiting a high density of connections. However, unlike genetic or metabolic networks typified by only a few highly connected nodes ("hubs"), most nodes in the miRNA-pathway network are highly connected. Sequence-based computational analysis confirmed that highly-interconnected miRNAs are likely to be regulated by common pathways to target similar sets of downstream genes, suggesting a pervasive and high level of functional redundancy among coexpressed miRNAs. We conclude that gene expression signatures can be used as surrogates of miRNA activity. Our strategy facilitates the task of discovering novel miRNA-pathway connections, since gene expression data for multiple normal and disease conditions are abundantly available.     

RNAseq expression analysis of resistant and susceptible mice after influenza A virus infection identifies novel genes associated with virus replication and important for host resistance to infection

Background

The host response to influenza A infections is strongly influenced by host genetic factors. Animal models of genetically diverse mouse strains are well suited to identify host genes involved in severe pathology, viral replication and immune responses. Here, we have utilized a dual RNAseq approach that allowed us to investigate both viral and host gene expression in the same individual mouse after H1N1 infection.

Results

We performed a detailed expression analysis to identify (i) correlations between changes in expression of host and virus genes, (ii) host genes involved in viral replication, and (iii) genes showing differential expression between two mouse strains that strongly differ in resistance to influenza infections. These genes may be key players involved in regulating the differences in pathogenesis and host defense mechanisms after influenza A infections. Expression levels of influenza segments correlated well with the viral load and may thus be used as surrogates for conventional viral load measurements. Furthermore, we investigated the functional role of two genes, Reg3g and Irf7, in knock-out mice and found that deletion of the Irf7 gene renders the host highly susceptible to H1N1 infection.

Conclusions

Using RNAseq analysis we identified novel genes important for viral replication or the host defense. This study adds further important knowledge to host-pathogen-interactions and suggests additional candidates that are crucial for host susceptibility or survival during influenza A infections.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1867-8) contains supplementary material, which is available to authorized users.     

Multiple Signals Converge on a Differentiation MAPK Pathway

An important emerging question in the area of signal transduction is how information from different pathways becomes integrated into a highly coordinated response. In budding yeast, multiple pathways regulate filamentous growth, a complex differentiation response that occurs under specific environmental conditions. To identify new aspects of filamentous growth regulation, we used a novel screening approach (called secretion profiling) that measures release of the extracellular domain of Msb2p, the signaling mucin which functions at the head of the filamentous growth (FG) MAPK pathway. Secretion profiling of complementary genomic collections showed that many of the pathways that regulate filamentous growth (RAS, RIM101, OPI1, and RTG) were also required for FG pathway activation. This regulation sensitized the FG pathway to multiple stimuli and synchronized it to the global signaling network. Several of the regulators were required for MSB2 expression, which identifies the MSB2 promoter as a target “hub” where multiple signals converge. Accessibility to the MSB2 promoter was further regulated by the histone deacetylase (HDAC) Rpd3p(L), which positively regulated FG pathway activity and filamentous growth. Our findings provide the first glimpse of a global regulatory hierarchy among the pathways that control filamentous growth. Systems-level integration of signaling circuitry is likely to coordinate other regulatory networks that control complex behaviors.