Farm animal exposure, respiratory illnesses, and nasal cell gene expression.

BACKGROUND: Farm exposures in early life reduce the risks for childhood allergic diseases and asthma. There is less information about how farm exposures relate to respiratory illnesses and mucosal immune development. OBJECTIVE: We hypothesized that children raised in farm environments have a lower incidence of respiratory illnesses over the first 2 years of life than nonfarm children. We also analyzed whether farm exposures or respiratory illnesses were related to patterns of nasal cell gene expression. METHODS: The Wisconsin Infant Study Cohort included farm (n = 156) and nonfarm (n = 155) families with children followed to age 2 years. Parents reported prenatal farm and other environmental exposures. Illness frequency and severity were assessed using illness diaries and periodic surveys. Nasopharyngeal cell gene expression in a subset of 64 children at age 2 years was compared to farm exposure and respiratory illness history. RESULTS: Farm versus nonfarm children had nominally lower rates of respiratory illnesses (rate ratio 0.82 [95% CI, 0.69, 0.97]) with a stepwise reduction in illness rates in children exposed to 0, 1, or ≥2 animal species, but these trends were nonsignificant in a multivariable model. Farm exposures and preceding respiratory illnesses were positively related to nasal cell gene signatures for mononuclear cells and innate and antimicrobial responses. CONCLUSIONS: Maternal and infant exposure to farms and farm animals was associated with nonsignificant trends for reduced respiratory illnesses. Nasal cell gene expression in a subset of children suggests that farm exposures and respiratory illnesses in early life are associated with distinct patterns of mucosal immune expression.

Proteogenomic and V(D)J Analysis of Human Decidual T Cells Highlights Unique Transcriptional Programming and Clonal Distribution.

Immunological tolerance toward the semiallogeneic fetus is one of many maternal adaptations required for a successful pregnancy. T cells are major players of the adaptive immune system and balance tolerance and protection at the maternal-fetal interface; however, their repertoire and subset programming are still poorly understood. Using emerging single-cell RNA sequencing technologies, we simultaneously obtained transcript, limited protein, and receptor repertoire at the single-cell level, from decidual and matched maternal peripheral human T cells. The decidua maintains a tissue-specific distribution of T cell subsets compared with the periphery. We find that decidual T cells maintain a unique transcriptome programming, characterized by restraint of inflammatory pathways by overexpression of negative regulators (DUSP, TNFAIP3, ZFP36) and expression of PD-1, CTLA-4, TIGIT, and LAG3 in some CD8 clusters. Finally, analyzing TCR clonotypes demonstrated decreased diversity in specific decidual T cell populations. Overall, our data demonstrate the power of multiomics analysis in revealing regulation of fetal-maternal immune coexistence.

Dynamic regulatory module networks for inference of cell type-specific transcriptional networks.

Changes in transcriptional regulatory networks can significantly alter cell fate. To gain insight into transcriptional dynamics, several studies have profiled bulk multi-omic data sets with parallel transcriptomic and epigenomic measurements at different stages of a developmental process. However, integrating these data to infer cell type-specific regulatory networks is a major challenge. We present dynamic regulatory module networks (DRMNs), a novel approach to infer cell type-specific cis-regulatory networks and their dynamics. DRMN integrates expression, chromatin state, and accessibility to predict cis-regulators of context-specific expression, where context can be cell type, developmental stage, or time point, and uses multitask learning to capture network dynamics across linearly and hierarchically related contexts. We applied DRMNs to study regulatory network dynamics in three developmental processes, each showing different temporal relationships and measuring a different combination of regulatory genomic data sets: cellular reprogramming, liver dedifferentiation, and forward differentiation. DRMN identified known and novel regulators driving cell type-specific expression patterns, showing its broad applicability to examine dynamics of gene regulatory networks from linearly and hierarchically related multi-omic data sets.

Multiomic analysis reveals decidual-specific transcriptional programing of MAIT cells.

PROBLEM: Mucosal-Associated Invariant T (MAIT) cells have been recently identified at the maternal-fetal interface. However, transcriptional programming of decidual MAIT cells in pregnancy remains poorly understood. METHOD OF STUDY: We employed a multiomic approach to address this question. Mononuclear cells from the decidua basalis and parietalis, and control PBMCs, were analyzed via flow cytometry to investigate MAIT cells in the decidua and assess their transcription factor expression. In a separate study, both decidual and matched peripheral MAIT cells were analyzed using Cellular Indexing of Transcriptomes and Epitopes by Sequencing (CITE-seq) coupled with gene expression analysis. Lastly, decidual MAIT cells were stimulated with E.coli and expression of MR1 by antigen presenting cells was measured to evaluate decidual MAIT cell function. RESULTS: First, we identified MAIT cells in both the decidua basalis and parietalis. CITE-seq, coupled with scRNA-seq gene expression analysis, highlighted transcriptional programming differences between decidual and matched peripheral MAIT cells at a single cell resolution. Transcription factor expression analysis further highlighted transcriptional differences between decidual MAIT cells and non-matched peripheral MAIT cells. Functionally, MAIT cells are skewed towards IFNγ and TNFα production upon stimulation, with E.coli leading to IFNγ production. Lastly, we demonstrate that MR1, the antigen presenting molecule restricting MAIT cells, is expressed by decidual APCs. CONCLUSION: MAIT cells are present in the decidua basalis and obtain a unique gene expression profile. The presence of MR1 on APCs coupled with in vitro activation by E.coli suggests that MAIT cells might be involved in tissue-repair mechanisms at the maternal-fetal interface.

Deciphering the Role of 3D Genome Organization in Breast Cancer Susceptibility.

Cancer risk by environmental exposure is modulated by an individual's genetics and age at exposure. This age-specific period of susceptibility is referred to as the "Window of Susceptibility" (WOS). Rats have a similar WOS for developing breast cancer. A previous study in rat identified an age-specific long-range regulatory interaction for the cancer gene, Pappa, that is associated with breast cancer susceptibility. However, the global role of three-dimensional genome organization and downstream gene expression programs in the WOS is not known. Therefore, we generated Hi-C and RNA-seq data in rat mammary epithelial cells within and outside the WOS. To systematically identify higher-order changes in 3D genome organization, we developed NE-MVNMF that combines network enhancement followed by multitask non-negative matrix factorization. We examined three-dimensional genome organization dynamics at the level of individual loops as well as higher-order domains. Differential chromatin interactions tend to be associated with differentially up-regulated genes with the WOS and recapitulate several human SNP-gene interactions associated with breast cancer susceptibility. Our approach identified genomic blocks of regions with greater overall differences in contact count between the two time points when the cluster assignments change and identified genes and pathways implicated in early carcinogenesis and cancer treatment. Our results suggest that WOS-specific changes in 3D genome organization are linked to transcriptional changes that may influence susceptibility to breast cancer.

Identification of FMR1-regulated molecular networks in human neurodevelopment.

RNA-binding proteins (RNA-BPs) play critical roles in development and disease to regulate gene expression. However, genome-wide identification of their targets in primary human cells has been challenging. Here, we applied a modified CLIP-seq strategy to identify genome-wide targets of the FMRP translational regulator 1 (FMR1), a brain-enriched RNA-BP, whose deficiency leads to Fragile X Syndrome (FXS), the most prevalent inherited intellectual disability. We identified FMR1 targets in human dorsal and ventral forebrain neural progenitors and excitatory and inhibitory neurons differentiated from human pluripotent stem cells. In parallel, we measured the transcriptomes of the same four cell types upon FMR1 gene deletion. We discovered that FMR1 preferentially binds long transcripts in human neural cells. FMR1 targets include genes unique to human neural cells and associated with clinical phenotypes of FXS and autism. Integrative network analysis using graph diffusion and multitask clustering of FMR1 CLIP-seq and transcriptional targets reveals critical pathways regulated by FMR1 in human neural development. Our results demonstrate that FMR1 regulates a common set of targets among different neural cell types but also operates in a cell type-specific manner targeting distinct sets of genes in human excitatory and inhibitory neural progenitors and neurons. By defining molecular subnetworks and validating specific high-priority genes, we identify novel components of the FMR1 regulation program. Our results provide new insights into gene regulation by a critical neuronal RNA-BP in human neurodevelopment.

In silico prediction of high-resolution Hi-C interaction matrices.

The three-dimensional (3D) organization of the genome plays an important role in gene regulation bringing distal sequence elements in 3D proximity to genes hundreds of kilobases away. Hi-C is a powerful genome-wide technique to study 3D genome organization. Owing to experimental costs, high resolution Hi-C datasets are limited to a few cell lines. Computational prediction of Hi-C counts can offer a scalable and inexpensive approach to examine 3D genome organization across multiple cellular contexts. Here we present HiC-Reg, an approach to predict contact counts from one-dimensional regulatory signals. HiC-Reg predictions identify topologically associating domains and significant interactions that are enriched for CCCTC-binding factor (CTCF) bidirectional motifs and interactions identified from complementary sources. CTCF and chromatin marks, especially repressive and elongation marks, are most important for HiC-Reg's predictive performance. Taken together, HiC-Reg provides a powerful framework to generate high-resolution profiles of contact counts that can be used to study individual locus level interactions and higher-order organizational units of the genome.

Inferring Regulatory Programs Governing Region Specificity of Neuroepithelial Stem Cells during Early Hindbrain and Spinal Cord Development.

Neuroepithelial stem cells (NSC) from different anatomical regions of the embryonic neural tube's rostrocaudal axis can differentiate into diverse central nervous system tissues, but the transcriptional regulatory networks governing these processes are incompletely understood. Here, we measure region-specific NSC gene expression along the rostrocaudal axis in a human pluripotent stem cell model of early central nervous system development over a 72-h time course, spanning the hindbrain to cervical spinal cord. We introduce Escarole, a probabilistic clustering algorithm for non-stationary time series, and combine it with prior-based regulatory network inference to identify genes that are regulated dynamically and predict their upstream regulators. We identify known regulators of patterning and neural development, including the HOX genes, and predict a direct regulatory connection between the transcription factor POU3F2 and target gene STMN2. We demonstrate that POU3F2 is required for expression of STMN2, suggesting that this regulatory connection is important for region specificity of NSCs.

Augmenting subnetwork inference with information extracted from the scientific literature.

Many biological studies involve either (i) manipulating some aspect of a cell or its environment and then simultaneously measuring the effect on thousands of genes, or (ii) systematically manipulating each gene and then measuring the effect on some response of interest. A common challenge that arises in these studies is to explain how genes identified as relevant in the given experiment are organized into a subnetwork that accounts for the response of interest. The task of inferring a subnetwork is typically dependent on the information available in publicly available, structured databases, which suffer from incompleteness. However, a wealth of potentially relevant information resides in the scientific literature, such as information about genes associated with certain concepts of interest, as well as interactions that occur among various biological entities. We contend that by exploiting this information, we can improve the explanatory power and accuracy of subnetwork inference in multiple applications. Here we propose and investigate several ways in which information extracted from the scientific literature can be used to augment subnetwork inference. We show that we can use literature-extracted information to (i) augment the set of entities identified as being relevant in a subnetwork inference task, (ii) augment the set of interactions used in the process, and (iii) support targeted browsing of a large inferred subnetwork by identifying entities and interactions that are closely related to concepts of interest. We use this approach to uncover the pathways involved in interactions between a virus and a host cell, and the pathways that are regulated by a transcription factor associated with breast cancer. Our experimental results demonstrate that these approaches can provide more accurate and more interpretable subnetworks. Integer program code, background network data, and pathfinding code are available at https://github.com/Craven-Biostat-Lab/subnetwork_inference.

Transcriptional and Functional Programming of Decidual Innate Lymphoid Cells.

A successful pregnancy requires many physiological adaptations from the mother, including the establishment of tolerance toward the semiallogeneic fetus. Innate lymphoid cells (ILCs) have arisen as important players in immune regulation and tissue homeostasis at mucosal and barrier surfaces. Dimensionality reduction and transcriptomic analysis revealed the presence of two novel CD56Bright decidual ILCs that express low T-bet and divergent Eomes levels. Transcriptional correlation with recently identified first trimester decidual dNKs suggests that these novel decidual ILCs might be present throughout pregnancy. Functional testing with permutation analysis revealed production of multiple factors by individual cells, with a preference for IFNγ and VEGF. Overall, our data suggests continuity of a unique decidual innate lymphocytes across pregnancy with a polyfunctional functional profile conducive for pregnancy.

Integrative Approaches for Inference of Genome-Scale Gene Regulatory Networks.

Transcriptional regulatory networks specify the regulatory proteins of target genes that control the context-specific expression levels of genes. With our ability to profile the different types of molecular components of cells under different conditions, we are now uniquely positioned to infer regulatory networks in diverse biological contexts such as different cell types, tissues, and time points. In this chapter, we cover two main classes of computational methods to integrate different types of information to infer genome-scale transcriptional regulatory networks. The first class of methods focuses on integrative methods for specifically inferring connections between transcription factors and target genes by combining gene expression data with regulatory edge-specific knowledge. The second class of methods integrates upstream signaling networks with transcriptional regulatory networks by combining gene expression data with protein-protein interaction networks and proteomic datasets. We conclude with a section on practical applications of a network inference algorithm to infer a genome-scale regulatory network.

Network inference reveals novel connections in pathways regulating growth and defense in the yeast salt response.

Cells respond to stressful conditions by coordinating a complex, multi-faceted response that spans many levels of physiology. Much of the response is coordinated by changes in protein phosphorylation. Although the regulators of transcriptome changes during stress are well characterized in Saccharomyces cerevisiae, the upstream regulatory network controlling protein phosphorylation is less well dissected. Here, we developed a computational approach to infer the signaling network that regulates phosphorylation changes in response to salt stress. We developed an approach to link predicted regulators to groups of likely co-regulated phospho-peptides responding to stress, thereby creating new edges in a background protein interaction network. We then use integer linear programming (ILP) to integrate wild type and mutant phospho-proteomic data and predict the network controlling stress-activated phospho-proteomic changes. The network we inferred predicted new regulatory connections between stress-activated and growth-regulating pathways and suggested mechanisms coordinating metabolism, cell-cycle progression, and growth during stress. We confirmed several network predictions with co-immunoprecipitations coupled with mass-spectrometry protein identification and mutant phospho-proteomic analysis. Results show that the cAMP-phosphodiesterase Pde2 physically interacts with many stress-regulated transcription factors targeted by PKA, and that reduced phosphorylation of those factors during stress requires the Rck2 kinase that we show physically interacts with Pde2. Together, our work shows how a high-quality computational network model can facilitate discovery of new pathway interactions during osmotic stress.

Integrative Genomic Analysis Predicts Causative Cis-Regulatory Mechanisms of the Breast Cancer-Associated Genetic Variant rs4415084.

Previous genome-wide association studies (GWAS) have identified several common genetic variants that may significantly modulate cancer susceptibility. However, the precise molecular mechanisms behind these associations remain largely unknown; it is often not clear whether discovered variants are themselves functional or merely genetically linked to other functional variants. Here, we provide an integrated method for identifying functional regulatory variants associated with cancer and their target genes by combining analyses of expression quantitative trait loci, a modified version of allele-specific expression that systematically utilizes haplotype information, transcription factor (TF)-binding preference, and epigenetic information. Application of our method to a breast cancer susceptibility region in 5p12 demonstrates that the risk allele rs4415084-T correlates with higher expression levels of the protein-coding gene mitochondrial ribosomal protein S30 (MRPS30) and lncRNA RP11-53O19.1 We propose an intergenic SNP rs4321755, in linkage disequilibrium (LD) with the GWAS SNP rs4415084 (r2 = 0.988), to be the predicted functional SNP. The risk allele rs4321755-T, in phase with the GWAS rs4415084-T, created a GATA3-binding motif within an enhancer, resulting in differential GATA3 binding and chromatin accessibility, thereby promoting transcription of MRPS30 and RP11-53O19.1. MRPS30 encodes a member of the mitochondrial ribosomal proteins, implicating the role of risk SNP in modulating mitochondrial activities in breast cancer. Our computational framework provides an effective means to integrate GWAS results with high-throughput genomic and epigenomic data and can be extended to facilitate rapid functional characterization of other genetic variants modulating cancer susceptibility.Significance: Unification of GWAS results with information from high-throughput genomic and epigenomic profiles provides a direct link between common genetic variants and measurable molecular perturbations. Cancer Res; 78(7); 1579-91. ©2018 AACR.

Inference of cell type specific regulatory networks on mammalian lineages.

Transcriptional regulatory networks are at the core of establishing cell type specific gene expression programs. In mammalian systems, such regulatory networks are determined by multiple levels of regulation, including by transcription factors, chromatin environment, and three-dimensional organization of the genome. Recent efforts to measure diverse regulatory genomic datasets across multiple cell types and tissues offer unprecedented opportunities to examine the context-specificity and dynamics of regulatory networks at a greater resolution and scale than before. In parallel, numerous computational approaches to analyze these data have emerged that serve as important tools for understanding mammalian cell type specific regulation. In this article, we review recent computational approaches to predict the expression and sequence-based regulators of a gene's expression level and examine long-range gene regulation. We highlight promising approaches, insights gained, and open challenges that need to be overcome to build a comprehensive picture of cell type specific transcriptional regulatory networks.

Physiological Responses and Gene Co-Expression Network of Mycorrhizal Roots under K+ Deprivation.

Arbuscular mycorrhizal (AM) associations enhance the phosphorous and nitrogen nutrition of host plants, but little is known about their role in potassium (K+) nutrition. Medicago truncatula plants were cocultured with the AM fungus Rhizophagus irregularis under high and low K+ regimes for 6 weeks. We determined how K+ deprivation affects plant development and mineral acquisition and how these negative effects are tempered by the AM colonization. The transcriptional response of AM roots under K+ deficiency was analyzed by whole-genome RNA sequencing. K+ deprivation decreased root biomass and external K+ uptake and modulated oxidative stress gene expression in M. truncatula roots. AM colonization induced specific transcriptional responses to K+ deprivation that seem to temper these negative effects. A gene network analysis revealed putative key regulators of these responses. This study confirmed that AM associations provide some tolerance to K+ deprivation to host plants, revealed that AM symbiosis modulates the expression of specific root genes to cope with this nutrient stress, and identified putative regulators participating in these tolerance mechanisms.

Integrating Transcriptomic and Proteomic Data Using Predictive Regulatory Network Models of Host Response to Pathogens.

Mammalian host response to pathogenic infections is controlled by a complex regulatory network connecting regulatory proteins such as transcription factors and signaling proteins to target genes. An important challenge in infectious disease research is to understand molecular similarities and differences in mammalian host response to diverse sets of pathogens. Recently, systems biology studies have produced rich collections of omic profiles measuring host response to infectious agents such as influenza viruses at multiple levels. To gain a comprehensive understanding of the regulatory network driving host response to multiple infectious agents, we integrated host transcriptomes and proteomes using a network-based approach. Our approach combines expression-based regulatory network inference, structured-sparsity based regression, and network information flow to infer putative physical regulatory programs for expression modules. We applied our approach to identify regulatory networks, modules and subnetworks that drive host response to multiple influenza infections. The inferred regulatory network and modules are significantly enriched for known pathways of immune response and implicate apoptosis, splicing, and interferon signaling processes in the differential response of viral infections of different pathogenicities. We used the learned network to prioritize regulators and study virus and time-point specific networks. RNAi-based knockdown of predicted regulators had significant impact on viral replication and include several previously unknown regulators. Taken together, our integrated analysis identified novel module level patterns that capture strain and pathogenicity-specific patterns of expression and helped identify important regulators of host response to influenza infection.

Network-based approaches for analysis of complex biological systems.

Cells function and respond to changes in their environment by the coordinated activity of their molecular components, including mRNAs, proteins and metabolites. At the heart of proper cellular function are molecular networks connecting these components to process extra-cellular environmental signals and drive dynamic, context-specific cellular responses. Network-based computational approaches aim to systematically integrate measurements from high-throughput experiments to gain a global understanding of cellular function under changing environmental conditions. We provide an overview of recent methodological developments toward solving two major computational problems within this field in the past two years (2013-2015): network reconstruction and network-based interpretation. Looking forward, we envision development of methods that can predict phenotypes with high accuracy as well as provide biologically plausible mechanistic hypotheses.

Multi-task consensus clustering of genome-wide transcriptomes from related biological conditions.

MOTIVATION: Identifying the shared and pathogen-specific components of host transcriptional regulatory programs is important for understanding the principles of regulation of immune response. Recent efforts in systems biology studies of infectious diseases have resulted in a large collection of datasets measuring host transcriptional response to various pathogens. Computational methods to identify and compare gene expression modules across different infections offer a powerful way to identify strain-specific and shared components of the regulatory program. An important challenge is to identify statistically robust gene expression modules as well as to reliably detect genes that change their module memberships between infections. RESULTS: We present MULCCH (MULti-task spectral Consensus Clustering for Hierarchically related tasks), a consensus extension of a multi-task clustering algorithm to infer high-confidence strain-specific host response modules under infections from multiple virus strains. On simulated data, MULCCH more accurately identifies genes exhibiting pathogen-specific patterns compared to non-consensus and nonmulti-task clustering approaches. Application of MULCCH to mammalian transcriptional response to a panel of influenza viruses showed that our method identifies clusters with greater coherence compared to non-consensus methods. Further, MULCCH derived clusters are enriched for several immune system-related processes and regulators. In summary, MULCCH provides a reliable module-based approach to identify molecular pathways and gene sets characterizing commonality and specificity of host response to viruses of different pathogenicities. AVAILABILITY AND IMPLEMENTATION: The source code is available at https://bitbucket.org/roygroup/mulcch CONTACT: sroy@biostat.wisc.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A predictive modeling approach for cell line-specific long-range regulatory interactions.

Long range regulatory interactions among distal enhancers and target genes are important for tissue-specific gene expression. Genome-scale identification of these interactions in a cell line-specific manner, especially using the fewest possible datasets, is a significant challenge. We develop a novel computational approach, Regulatory Interaction Prediction for Promoters and Long-range Enhancers (RIPPLE), that integrates published Chromosome Conformation Capture (3C) data sets with a minimal set of regulatory genomic data sets to predict enhancer-promoter interactions in a cell line-specific manner. Our results suggest that CTCF, RAD21, a general transcription factor (TBP) and activating chromatin marks are important determinants of enhancer-promoter interactions. To predict interactions in a new cell line and to generate genome-wide interaction maps, we develop an ensemble version of RIPPLE and apply it to generate interactions in five human cell lines. Computational validation of these predictions using existing ChIA-PET and Hi-C data sets showed that RIPPLE accurately predicts interactions among enhancers and promoters. Enhancer-promoter interactions tend to be organized into subnetworks representing coordinately regulated sets of genes that are enriched for specific biological processes and cis-regulatory elements. Overall, our work provides a systematic approach to predict and interpret enhancer-promoter interactions in a genome-wide cell-type specific manner using a few experimentally tractable measurements.