Influenza A virus (IAV) infection remodels cellular processes to support viral replication. The modulation of host factors by the virus drives pathogenesis during infection, and these factors may serve as therapeutic targets. Here, we infect mice with two IAV strains, H1N1 and H5N1, and analyze lung tissue with multi-omics. Using network propagation analysis, we identify twenty-four distinct host modules altered by infection, encompassing 2920 genes/proteins. Independently, we develop a computational pipeline, MidTOD, which integrates metabolomic data with other OMICs data-types, linking metabolites to gene/protein alterations. Combining datasets from both approaches reveals alterations in mitochondrial and peroxisomal metabolism in IAV-infected cells and identifies arginine:glycine amidinotransferase (GATM) as a host dependency factor in both human cells and mice. Knockdown of this enzyme reduces IAV-mediated pathology and host inflammatory responses after infection. Collectively, this work provides an integrated systems-level view of host changes during infection and identifies an abundance of IAV-host factors.

Hybrid immunity (vaccination + natural infection) to SARS-CoV-2 provides superior protection to re-infection. We performed immune profiling studies during breakthrough infections in mRNA-vaccinated hamsters to evaluate hybrid immunity induction. The mRNA vaccine, BNT162b2, was dosed to induce binding antibody titers against ancestral spike, but inefficient serum virus neutralization of ancestral SARS-CoV-2 or variants of concern (VoCs). Vaccination reduced morbidity and controlled lung virus titers for ancestral virus and Alpha but allowed breakthrough infections in Beta, Delta and Mu-challenged hamsters. Vaccination primed for T cell responses that were boosted by infection. Infection back-boosted neutralizing antibody responses against ancestral virus and VoCs. Hybrid immunity resulted in more cross-reactive sera, reflected by smaller antigenic cartography distances. Transcriptomics post-infection reflects both vaccination status and disease course and suggests a role for interstitial macrophages in vaccine-mediated protection. Therefore, protection by vaccination, even in the absence of high titers of neutralizing antibodies in the serum, correlates with recall of broadly reactive B- and T-cell responses.

Many machine learning applications in bioinformatics currently rely on matching gene identities when analyzing input gene signatures and fail to take advantage of preexisting knowledge about gene functions. To further enable comparative analysis of OMICS datasets, including target deconvolution and mechanism of action studies, we develop an approach that represents gene signatures projected onto their biological functions, instead of their identities, similar to how the word2vec technique works in natural language processing. We develop the Functional Representation of Gene Signatures (FRoGS) approach by training a deep learning model and demonstrate that its application to the Broad Institute’s L1000 datasets results in more effective compound-target predictions than models based on gene identities alone. By integrating additional pharmacological activity data sources, FRoGS significantly increases the number of high-quality compound-target predictions relative to existing approaches, many of which are supported by in silico and/or experimental evidence. These results underscore the general utility of FRoGS in machine learning-based bioinformatics applications. Prediction networks pre-equipped with the knowledge of gene functions may help uncover new relationships among gene signatures acquired by large-scale OMICs studies on compounds, cell types, disease models, and patient cohorts.

Triclosan (TCS) is an antimicrobial toxicant found in a myriad of consumer products and has been detected in human tissues, including breastmilk. We have evaluated the impact of lactational TCS on UDP-glucuronosyltransferase 1A1 (UGT1A1) expression and bilirubin metabolism in humanized UGT1 ( hUGT1 ) neonatal mice. In hUGT1 mice, expression of the hepatic UGT1A1 gene is developmentally delayed resulting in elevated total serum bilirubin (TSB) levels. We found that newborn hUGT1 mice breastfed or orally treated with TCS presented lower TSB levels along with induction of hepatic UGT1A1. Lactational and oral treatment by gavage with TCS leads to the activation of hepatic nuclear receptors constitutive androstane receptor (CAR), peroxisome proliferator-activated receptor alpha (PPARa), and stress sensor, activating transcription factor 4 (ATF4). When CAR-deficient hUGT1 mice ( hUGT1/Car - / - ) were treated with TCS, TSB levels were reduced with a robust induction of hepatic UGT1A1, leaving us to conclude that CAR is not tied to UGT1A1 induction. Alternatively, when PPARa-deficient hUGT1 mice ( hUGT1/ Ppara-/-) were treated with TCS, hepatic UGT1A1 was not induced. Additionally, we had previously demonstrated that TCS is a potent inducer of ATF4, a transcriptional factor linked to the integrated stress response. When ATF4 was deleted in liver of hUGT1 mice ( hUGT1/Atf4 D Hep ) and these mice treated with TCS, we observed superinduction of hepatic UGT1A1. Oxidative stress genes in livers of hUGT1/Atf4DHeptreated with TCS were increased, suggesting that ATF4 protects liver from excessive oxidative stress. The increase oxidative stress may be associated with superinduction of UGT1A1. The expression of ATF4 in neonatal hUGT1 hepatic tissue may play a role in the developmental repression of UGT1A1.

Patterns of transcriptional activity are encoded in our genome through regulatory elements such as promoters or enhancers that, paradoxically, contain similar assortments of sequence-specific transcription factor (TF) binding sites1-3. Knowledge of how these sequence motifs encode multiple, often overlapping, gene expression programs is central to understanding gene regulation and how mutations in non-coding DNA manifest in disease4,5. Here, by studying gene regulation from the perspective of individual transcription start sites (TSSs), using natural genetic variation, perturbation of endogenous TF protein levels and massively parallel analysis of natural and synthetic regulatory elements, we show that the effect of TF binding on transcription initiation is position dependent. Analysing TF-binding-site occurrences relative to the TSS, we identified several motifs with highly preferential positioning. We show that these patterns are a combination of a TF's distinct functional profiles-many TFs, including canonical activators such as NRF1, NFY and Sp1, activate or repress transcription initiation depending on their precise position relative to the TSS. As such, TFs and their spacing collectively guide the site and frequency of transcription initiation. More broadly, these findings reveal how similar assortments of TF binding sites can generate distinct gene regulatory outcomes depending on their spatial configuration and how DNA sequence polymorphisms may contribute to transcription variation and disease and underscore a critical role for TSS data in decoding the regulatory information of our genome.

Background: Newborns can be exposed to inorganic arsenic (iAs) through contaminated drinking water, formula, and other infant foods. Epidemiological studies have demonstrated a positive association between urinary iAs levels and the risk of developing nonalcoholic fatty liver disease (NAFLD) among U.S. adolescents and adults. Objectives: The present study examined how oral iAs administration to neonatal mice impacts the intestinal tract, which acts as an early mediator for NAFLD. Methods: Neonatal mice were treated with a single dose of iAs via oral gavage. Effects on the small intestine were determined by histological examination, RNA sequencing, and biochemical analysis. Serum lipid profiling was analyzed by fast protein liquid chromatography (FPLC), and hepatosteatosis was characterized histologically and biochemically. Liver X receptor-alpha (LXRα) knockout (Lxrα−/−) mice and liver-specific activating transcription factor 4 (ATF4)-deficient (Atf4ΔHep) mice were used to define their roles in iAs-induced effects during the neonatal stage. Results: Neonatal mice exposed to iAs via oral gavage exhibited accumulation of dietary fat in enterocytes, with higher levels of enterocyte triglycerides and free fatty acids. These mice also showed accelerated enterocyte maturation and a longer small intestine. This was accompanied by higher levels of liver-derived very low-density lipoprotein and low-density lipoprotein triglycerides, and a lower level of high-density lipoprotein cholesterol in the serum. Mice exposed during the neonatal period to oral iAs also developed hepatosteatosis. Compared with the control group, iAs-induced fat accumulation in enterocytes became more significant in neonatal Lxrα−/− mice, accompanied by accelerated intestinal growth, hypertriglyceridemia, and hepatosteatosis. In contrast, regardless of enterocyte fat accumulation, hepatosteatosis was largely reduced in iAs-treated neonatal Atf4ΔHep mice. Conclusion: Exposure to iAs in neonatal mice resulted in excessive accumulation of fat in enterocytes, disrupting lipid homeostasis in the serum and liver, revealing the importance of the gut–liver axis and endoplasmic reticulum stress in mediating iAs-induced NAFLD at an early age. https://doi.org/10.1289/EHP12381

Inorganic arsenic (iAs) is an environmental toxicant that can lead to severe health consequences, which can be exacerbated if exposure occurs early in development. Here, we evaluated the impact of oral iAs treatment on UDP-glucuronosyltransferase 1A1 (UGT1A1) expression and bilirubin metabolism in hu-manized UGT1 (hUGT1) mice. We found that oral adminis-tration of iAs to neonatal hUGT1 mice that display severe neonatal hyperbilirubinemia leads to induction of intestinal UGT1A1 and a reduction in total serum bilirubin values. Oral iAs administration accelerates neonatal intestinal maturation, an event that is directly associated with UGT1A1 induction. As a reactive oxygen species producer, oral iAs treatment activated the Keap-Nrf2 pathway in the intestinal tract and liver. When Nrf2-deficient hUGT1 mice (hUGT1/Nrf2-/-) were treated with iAs, it was shown that activated Nrf2 contributed significantly toward intestinal maturation and UGT1A1 induction. Howev-er, hepatic UGT1A1 was not induced upon iAs exposure. We previously demonstrated that the nuclear receptor PXR re-presses liver UGT1A1 in neonatal hUGT1 mice. When PXR was deleted in hUGT1 mice (hUGT1/Pxr-/-), derepression of UGT1A1 was evident in both liver and intestinal tissue in ne-onates. Furthermore, when neonatal hUGT1/Pxr-/- mice were treated with iAs, UGT1A1 was superinduced in both tissues, confirming PXR release derepressed key regulatory elements on the gene that could be activated by iAs exposure. With iAs capable of generating reactive oxygen species in both liver and intestinal tissue, we conclude that PXR deficiency in neonatal hUGT1/Pxr-/- mice allows greater access of activated tran-scriptional modifiers such as Nrf2 leading to superinduction of UGT1A1.

When mRNA vaccines were rolled out under EUA, it was observed that a single vaccination, already protected vaccinees from severe disease, even before neutralizing antibodies appearded. To better understand the underlying mechanism of protection, we vaccinated hamsters once with a suboptimal dose of Pfizer/BioNTech COVID-19 mRNA vaccine (5mg/animal). Four weeks post vaccination, animals were challenged with different variants of SARS-CoV-2. Vaccination resulted in antibody ELISA titers against SARS-CoV-2 spike protein, but antibody titers were too low for efficient neutralization of antigenically matching USA-WA1/2020 (WA1) or drifted variants of concern (Alpha, Beta, Delta and Mu). Nevertheless, suboptimal vaccination resulted in reduced morbidity for WA1 and Alpha-challenged animals and complete control of lung virus titers for WA1, Alpha and Delta but with breakthrough infection for the antigenically more distant Beta and Mu. Moreover, T cell responses were higher in vaccinated animals compared to unvaccinated ones at 5 days post infection, suggesting a recall of vaccine-primed T cells. Infection with different SARS-CoV-2 variants also back-boosted neutralizing antibody responses against the challenge virus and in the case of the variants of concern, against the antigenically distant but vaccine-matched WA1. Transcriptomic analysis of host immune responses to infection reflects both vaccination status and disease course, is further compared with lung pathology and suggests a role for interstitial lung macrophages in vaccine-mediated protection. In summary, protection after suboptimal vaccination against SARS-CoV-2 challenge in hamsters correlates with recall of both vaccine-induced B- and T-cell responses. This study was also partly funded by CRIPT (Center for Research on Influenza Pathogenesis and Transmission), a NIH NIAID funded Center of Excellence for Influenza Research and Response (CEIRR, contract number 75N93021C00014) to A.G.S and M.S. and by R01 AI160706 to M.S.

Transcription factor programs mediating the immune response to coronavirus disease 2019 (COVID-19) are not fully understood. Capturing active transcription initiation from cis-regulatory elements such as enhancers and promoters by capped small RNA sequencing (csRNA-seq), in contrast to capturing steady-state tran-scripts by conventional RNA-seq, allows unbiased identification of the underlying transcription factor activity and regulatory pathways. Here, we profile transcription initiation in critically ill COVID-19 patients, identifying transcription factor motifs that correlate with clinical lung injury and disease severity. Unbiased clustering re-veals distinct subsets of cis-regulatory elements that delineate the cell type, pathway-specific, and combi-natorial transcription factor activity. We find evidence of critical roles of regulatory networks, showing that STAT/BCL6 and E2F/MYB regulatory programs from myeloid cell populations are activated in patients with poor disease outcomes and associated with COVID-19 susceptibility genetic variants. More broadly, we demonstrate how capturing acute, disease-mediated changes in transcription initiation can provide insight into the underlying molecular mechanisms and stratify patient disease severity.

Substance abuse and addiction represent a significant public health problem that impacts multiple dimensions of society, including healthcare, the economy, and the workforce. In 2021, over 100,000 drug overdose deaths were reported in the US, with an alarming increase in fatalities related to opioids and psychostimulants. Understanding the fundamental gene regulatory mechanisms underlying addiction and related behaviors could facilitate more effective treatments. To explore how repeated drug exposure alters gene regulatory networks in the brain, we combined capped small (cs)RNA-seq, which accurately captures nascent-like initiating transcripts from total RNA, with Hi-C and single nuclei (sn)ATAC-seq. We profiled initiating transcripts in two addiction-related brain regions, the prefrontal cortex (PFC) and the nucleus accumbens (NAc), from rats that were never exposed to drugs or were subjected to prolonged abstinence after oxycodone or cocaine intravenous self-administration (IVSA). Interrogating over 100,000 active transcription start regions (TSRs) revealed that most TSRs had hallmarks of bonafide enhancers and highlighted the KLF/SP1, RFX, and AP1 transcription factors families as central to establishing brain-specific gene regulatory programs. Analysis of rats with addiction-like behaviors versus controls identified addiction-associated repression of transcription at regulatory enhancers recognized by nuclear receptor subfamily 3 group C (NR3C) factors, including glucocorticoid receptors. Cell-type deconvolution analysis using snATAC-seq uncovered a potential role of glial cells in driving the gene regulatory programs associated with addiction-related phenotypes. These findings highlight the power of advanced transcriptomics methods to provide insight into how addiction perturbs gene regulatory programs in the brain.

Compositional oscillations of the gut microbiome are essential for normal peripheral circadian rhythms, both of which are disrupted in diet-induced obesity (DIO). Although time-restricted feeding (TRF) maintains circadian synchrony and protects against DIO, its impact on the dynamics of the cecal gut microbiome is modest. Thus, other regions of the gut, particularly the ileum, the nexus for incretin and bile acid signaling, may play an important role in entraining peripheral circadian rhythms. We demonstrate the effect of diet and feeding rhythms on the ileal microbiome composition and transcriptome in mice. The dynamic rhythms of ileal microbiome composition and transcriptome are dampened in DIO. TRF partially restores diurnal rhythms of the ileal microbiome and transcriptome, increases GLP-1 release, and alters the ileal bile acid pool and farnesoid X receptor (FXR) signaling, which could explain how TRF exerts its metabolic benefits. Finally, we provide a web resource for exploration of ileal microbiome and transcriptome circadian data.

Ly6C hi inflammatory monocytes show high IFN responses, and contribute to both protective and pathogenic functions following influenza virus infection. In order to understand the significance of IFN responses in this subset, we examined monocytes during infection with a lethal H5N1 virus that induces high levels of IFN and a low-pathogenicity H1N1 virus that induces low levels of IFN. We show that H5N1 infection results in early recruitment of high numbers of Ly6C hi monocytes and induction of chemokines and Ifnb1 . Using unbiased transcriptomic and proteomic approaches, we also find that monocytes are significantly enriched during H5N1 infection and are associated with chemokine and IFN signatures in mice, and with severity of symptoms after influenza virus infection in humans. Recruited Ly6C hi monocytes subsequently become infected in the lung, produce IFN-β, and mature into FasL + monocyte-derived cells (FasL + MCs) expressing dendritic cell markers. Both Ccr2 -/- and Fasl gld mice are protected from lethal infection, indicating that monocytes contribute to pathogenesis. Global loss of type I and type III IFN signaling in Stat2 -/- mice results in loss of monocyte recruitment, likely reflecting a requirement for IFN-dependent chemokine induction. Here we show that IFN is not directly required for monocyte recruitment on an IFN-sufficient background, but is required for maturation to FasL + MCs. Loss of IFN signaling skews to a Ly6C lo phenotype associated with tissue repair, suggesting that IFN signaling in monocytes is a critical determinant of influenza virus pathogenesis.

A deficient interferon (IFN) response to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection has been implicated as a determinant of severe coronavirus disease 2019 (COVID-19). To identify the molecular effectors that govern IFN control of SARS-CoV-2 infection, we conducted a large-scale gain-of-function analysis that evaluated the impact of human IFN-stimulated genes (ISGs) on viral replication. A limited subset of ISGs were found to control viral infection, including endosomal factors inhibiting viral entry, RNA binding proteins suppressing viral RNA synthesis, and a highly enriched cluster of endoplasmic reticulum (ER)/Golgi-resident ISGs inhibiting viral assembly/egress. These included broad-acting antiviral ISGs and eight ISGs that specifically inhibited SARS-CoV-2 and SARS-CoV-1 replication. Among the broad-acting ISGs was BST2/tetherin, which impeded viral release and is antagonized by SARS-CoV-2 Orf7a protein. Overall, these data illuminate a set of ISGs that underlie innate immune control of SARS-CoV-2/SARS-CoV-1 infection, which will facilitate the understanding of host determinants that impact disease severity and offer potential therapeutic strategies for COVID-19.

The fate of influenza A virus (IAV) infection in the host cell depends on the balance between cellular defence mechanisms and viral evasion strategies. To illuminate the landscape of IAV cellular restriction, we generated and integrated global genetic loss-of-function screens with transcriptomics and proteomics data. Our multi-omics analysis revealed a subset of both IFN-dependent and independent cellular defence mechanisms that inhibit IAV replication. Amongst these, the autophagy regulator TBC1 domain family member 5 (TBC1D5), which binds Rab7 to enable fusion of autophagosomes and lysosomes, was found to control IAV replication in vitro and in vivo and to promote lysosomal targeting of IAV M2 protein. Notably, IAV M2 was observed to abrogate TBC1D5-Rab7 binding through a physical interaction with TBC1D5 via its cytoplasmic tail. Our results provide evidence for the molecular mechanism utilised by IAV M2 protein to escape lysosomal degradation and traffic to the cell membrane, where it supports IAV budding and growth.

The contribution of transcription factors (TFs) and gene regulatory programs in the immune response to COVID-19 and their relationship to disease outcome is not fully understood. Analysis of genome-wide changes in transcription at both promoter-proximal and distal cis-regulatory DNA elements, collectively termed the ’active cistrome,’ offers an unbiased assessment of TF activity identifying key pathways regulated in homeostasis or disease. Here, we profiled the active cistrome from peripheral leukocytes of critically ill COVID-19 patients to identify major regulatory programs and their dynamics during SARS-CoV-2 associated acute respiratory distress syndrome (ARDS). We identified TF motifs that track the severity of COVID- 19 lung injury, disease resolution, and outcome. We used unbiased clustering to reveal distinct cistrome subsets delineating the regulation of pathways, cell types, and the combinatorial activity of TFs. We found critical roles for regulatory networks driven by stimulus and lineage determining TFs, showing that STAT and E2F/MYB regulatory programs targeting myeloid cells are activated in patients with poor disease outcomes and associated with single nucleotide genetic variants implicated in COVID-19 susceptibility. Integration with single-cell RNA-seq found that STAT and E2F/MYB activation converged in specific neutrophils subset found in patients with severe disease. Collectively we demonstrate that cistrome analysis facilitates insight into disease mechanisms and provides an unbiased approach to evaluate global changes in transcription factor activity and stratify patient disease severity.

Influenza A virus (IAV) is a human respiratory pathogen that causes yearly global epidemics, as well as sporadic pandemics due to human adaptation of pathogenic strains. Efficient replication of IAV in different species is, in part, dictated by its ability to exploit the genetic environment of the host cell. To investigate IAV tropism in human cells, we evaluated the replication of IAV strains in a diverse subset of epithelial cell lines. HeLa cells were refractory to the growth of human H1N1 and H3N2 viruses and low-pathogenic avian influenza (LPAI) viruses. Interestingly, a human isolate of the highly pathogenic avian influenza (HPAI) H5N1 virus successfully propagated in HeLa cells to levels comparable to those in a human lung cell line. Heterokaryon cells generated by fusion of HeLa and permissive cells supported H1N1 virus growth, suggesting the absence of a host factor(s) required for the replication of H1N1, but not H5N1, viruses in HeLa cells. The absence of this factor(s) was mapped to reduced nuclear import, replication, and translation, as well as deficient viral budding. Using reassortant H1N1:H5N1 viruses, we found that the combined introduction of nucleoprotein (NP) and hemagglutinin (HA) from an H5N1 virus was necessary and sufficient to enable H1N1 virus growth. Overall, this study suggests that the absence of one or more cellular factors in HeLa cells results in abortive replication of H1N1, H3N2, and LPAI viruses, which can be circumvented upon the introduction of H5N1 virus NP and HA. Further understanding of the molecular basis of this restriction will provide important insights into the virus-host interactions that underlie IAV pathogenesis and tropism.IMPORTANCE Many zoonotic avian influenza A viruses have successfully crossed the species barrier and caused mild to life-threatening disease in humans. While human-to-human transmission is limited, there is a risk that these zoonotic viruses may acquire adaptive mutations enabling them to propagate efficiently and cause devastating human pandemics. Therefore, it is important to identify viral determinants that provide these viruses with a replicative advantage in human cells. Here, we tested the growth of influenza A virus in a subset of human cell lines and found that abortive replication of H1N1 viruses in HeLa cells can be circumvented upon the introduction of H5N1 virus HA and NP. Overall, this work leverages the genetic diversity of multiple human cell lines to highlight viral determinants that could contribute to H5N1 virus pathogenesis and tropism.

SBV is a newly emerging orthobunyavirus (family Peribunyaviridae ) that has spread rapidly across Europe since 2011, resulting in substantial economic losses in livestock farming. In this study, we performed unbiased genome-wide CRISPR-Cas9 screening and identified PAPST1, a sulfotransferase encoded by SLC35B2 , as a host entry factor for SBV. Consistent with its role in the synthesis of heparan sulfate, we show that this activity is required for efficient infection by SBV. A comparable dependency on heparan sulfate was also observed for La Crosse virus and Rift Valley fever virus, highlighting the importance of heparan sulfate for host cell infection by bunyaviruses. Thus, the present work provides crucial insights into virus-host interactions of important animal and human pathogens.