ETV7 limits antiviral gene expression and control of influenza viruses

Smoking-associated gene expression alterations in nasal epithelium reveal immune impairment linked to lung cancer risk

BACKGROUND

Lung cancer is the leading cause of cancer-related death in the world. In contrast to many other cancers, a direct connection to modifiable lifestyle risk in the form of tobacco smoke has long been established. More than 50% of all smoking-related lung cancers occur in former smokers, 40% of which occur more than 15 years after smoking cessation. Despite extensive research, the molecular processes for persistent lung cancer risk remain unclear. We thus set out to examine whether risk stratification in the clinic and in the general population can be improved upon by the addition of genetic data and to explore the mechanisms of the persisting risk in former smokers.

METHODS

We analysed transcriptomic data from accessible airway tissues of 487 subjects, including healthy volunteers and clinic patients of different smoking statuses. We developed a computational model to assess smoking-associated gene expression changes and their reversibility after smoking is stopped, comparing healthy subjects to clinic patients with and without lung cancer.

RESULTS

We find persistent smoking-associated immune alterations to be a hallmark of the clinic patients. Integrating previous GWAS data using a transcriptional network approach, we demonstrate that the same immune- and interferon-related pathways are strongly enriched for genes linked to known genetic risk factors, demonstrating a causal relationship between immune alteration and lung cancer risk. Finally, we used accessible airway transcriptomic data to derive a non-invasive lung cancer risk classifier.

CONCLUSIONS

Our results provide initial evidence for germline-mediated personalized smoke injury response and risk in the general population, with potential implications for managing long-term lung cancer incidence and mortality.

Epigenetic variation impacts individual differences in the transcriptional response to influenza infection

Humans display remarkable interindividual variation in their immune response to identical challenges. Yet, our understanding of the genetic and epigenetic factors contributing to such variation remains limited. Here we performed in-depth genetic, epigenetic and transcriptional profiling on primary macrophages derived from individuals of European and African ancestry before and after infection with influenza A virus. We show that baseline epigenetic profiles are strongly predictive of the transcriptional response to influenza A virus across individuals. Quantitative trait locus (QTL) mapping revealed highly coordinated genetic effects on gene regulation, with many cis-acting genetic variants impacting concomitantly gene expression and multiple epigenetic marks. These data reveal that ancestry-associated differences in the epigenetic landscape can be genetically controlled, even more than gene expression. Lastly, among QTL variants that colocalized with immune-disease loci, only 7% were gene expression QTL, while the remaining genetic variants impact epigenetic marks, stressing the importance of considering molecular phenotypes beyond gene expression in disease-focused studies.

Identification of early gene expression profiles associated with long-lasting antibody responses to the Ebola vaccine Ad26.ZEBOV/MVA-BN-Filo

Ebola virus disease is a severe hemorrhagic fever with a high fatality rate. We investigate transcriptome profiles at 3 h, 1 day, and 7 days after vaccination with Ad26.ZEBOV and MVA-BN-Filo. 3 h after Ad26.ZEBOV injection, we observe an increase in genes related to antigen presentation, sensing, and T and B cell receptors. The highest response occurs 1 day after Ad26.ZEBOV injection, with an increase of the gene expression of interferon-induced antiviral molecules, monocyte activation, and sensing receptors. This response is regulated by the HESX1, ATF3, ANKRD22, and ETV7 transcription factors. A plasma cell signature is observed on day 7 post-Ad26.ZEBOV vaccination, with an increase of CD138, MZB1, CD38, CD79A, and immunoglobulin genes. We have identified early expressed genes correlated with the magnitude of the antibody response 21 days after the MVA-BN-Filo and 364 days after Ad26.ZEBOV vaccinations. Our results provide early gene signatures that correlate with vaccine-induced Ebola virus glycoprotein-specific antibodies.

CELLULAR RNA INTERACTS WITH MAVS TO PROMOTE ANTIVIRAL SIGNALING

Immune signaling needs to be well-regulated to promote clearance of pathogens, while preventing aberrant inflammation. Interferons (IFNs) and antiviral genes are activated by the detection of viral RNA by RIG-I-like receptors (RLRs). Signal transduction downstream of RLRs proceeds through a multi-protein complex organized around the central adaptor protein MAVS. Recent work has shown that protein complex function can be modulated by RNA molecules providing allosteric regulation or acting as molecular guides or scaffolds. Thus, we hypothesized that RNA plays a role in organizing MAVS signaling platforms. Here, we show that MAVS, through its central intrinsically disordered domain, directly interacts with the 3' untranslated regions of cellular mRNAs. Importantly, elimination of RNA by RNase treatment disrupts the MAVS signalosome, including newly identified regulators of RLR signaling, and inhibits phosphorylation of the transcription factor IRF3. This supports the hypothesis that RNA molecules scaffold proteins in the MAVS signalosome to induce IFNs. Together, this work uncovers a function for cellular RNA in promoting signaling through MAVS and highlights a generalizable principle of RNA regulatory control of cytoplasmic immune signaling complexes.

MEF2A suppresses stress responses that trigger DDX41-dependent IFN production

Cellular stress in the form of disrupted transcription, loss of organelle integrity, or damage to nucleic acids can elicit inflammatory responses by activating signaling cascades canonically tasked with controlling pathogen infections. These stressors must be kept in check to prevent unscheduled activation of interferon, which contributes to autoinflammation. This study examines the role of the transcription factor myocyte enhancing factor 2A (MEF2A) in setting the threshold of transcriptional stress responses to prevent R-loop accumulation. Increases in R-loops lead to the induction of interferon and inflammatory responses in a DEAD-box helicase 41 (DDX41)-, cyclic GMP-AMP synthase (cGAS)-, and stimulator of interferon genes (STING)-dependent manner. The loss of MEF2A results in the activation of ATM and RAD3-related (ATR) kinase, which is also necessary for the activation of STING. This study identifies the role of MEF2A in sustaining transcriptional homeostasis and highlights the role of ATR in positively regulating R-loop-associated inflammatory responses.

Gene expression of peripheral blood mononuclear cells and CD8+ T cells from gilts after PRRSV infection

Porcine reproductive and respiratory syndrome virus (PRRSV) is a positive-stranded RNA virus, which emerged in Europe and U.S.A. in the late 1980s and has since caused huge economic losses. Infection with PRRSV causes mild to severe respiratory and reproductive clinical symptoms in pigs. Alteration of the host immune response by PRRSV is associated with the increased susceptibility to secondary viral and bacterial infections resulting in more serious and chronic disease. However, the expression profiles underlying innate and adaptive immune responses to PRRSV infection are yet to be further elucidated. In this study, we investigated gene expression profiles of PBMCs and CD8⁺ T cells after PRRSV AUT15-33 infection. We identified the highest number of differentially expressed genes in PBMCs and CD8⁺ T cells at 7 dpi and 21 dpi, respectively. The gene expression profile of PBMCs from infected animals was dominated by a strong innate immune response at 7 dpi which persisted through 14 dpi and 21 dpi and was accompanied by involvement of adaptive immunity. The gene expression pattern of CD8⁺ T cells showed a strong adaptive immune response to PRRSV, leading to the formation of highly differentiated CD8⁺ T cells starting from 14 dpi. The hallmark of the CD8⁺ T-cell response was the increased expression of effector and cytolytic genes (PRF1, GZMA, GZMB, GZMK, KLRK1, KLRD1, FASL, NKG7), with the highest levels observed at 21 dpi. Temporal clustering analysis of DEGs of PBMCs and CD8⁺ T cells from PRRSV-infected animals revealed three and four clusters, respectively, suggesting tight transcriptional regulation of both the innate and the adaptive immune response to PRRSV. The main cluster of PBMCs was related to the innate immune response to PRRSV, while the main clusters of CD8⁺ T cells represented the initial transformation and differentiation of these cells in response to the PRRSV infection. Together, we provided extensive transcriptomics data explaining gene signatures of the immune response of PBMCs and CD8⁺ T cells after PRRSV infection. Additionally, our study provides potential biomarker targets useful for vaccine and therapeutics development.

RNAseq profiling of blood from patients with coronary artery disease: Signature of a T cell imbalance

Background

Cardiovascular disease had a global prevalence of 523 million cases and 18.6 million deaths in 2019. The current standard for diagnosing coronary artery disease (CAD) is coronary angiography either by invasive catheterization (ICA) or computed tomography (CTA). Prior studies employed single-molecule, amplification-independent RNA sequencing of whole blood to identify an RNA signature in patients with angiographically confirmed CAD. The present studies employed Illumina RNAseq and network co-expression analysis to identify systematic changes underlying CAD.

Methods

Whole blood RNA was depleted of ribosomal RNA (rRNA) and analyzed by Illumina total RNA sequencing (RNAseq) to identify transcripts associated with CAD in 177 patients presenting for elective invasive coronary catheterization. The resulting transcript counts were compared between groups to identify differentially expressed genes (DEGs) and to identify patterns of changes through whole genome co-expression network analysis (WGCNA).

Results

The correlation between Illumina amplified RNAseq and the prior SeqLL unamplified RNAseq was quite strong (r = 0.87), but there was only 9 % overlap in the DEGs identified. Consistent with the prior RNAseq, the majority (93 %) of DEGs were down-regulated ~1.7-fold in patients with moderate to severe CAD (>20 % stenosis). DEGs were predominantly related to T cells, consistent with known reductions in Tregs in CAD. Network analysis did not identify pre-existing modules with a strong association with CAD, but patterns of T cell dysregulation were evident. DEGs were enriched for transcripts associated with ciliary and synaptic transcripts, consistent with changes in the immune synapse of developing T cells.

Conclusions

These studies confirm and extend a novel mRNA signature of a Treg-like defect in CAD. The pattern of changes is consistent with stress-related changes in the maturation of T and Treg cells, possibly due to changes in the immune synapse.

ETV7 reduces inflammatory responses in breast cancer cells by repressing the TNFR1/NF-κB axis

The transcription factor ETV7 is an oncoprotein that is up-regulated in all breast cancer (BC) types. We have recently demonstrated that ETV7 promoted breast cancer progression by increasing cancer cell proliferation and stemness and was also involved in the development of chemo- and radio-resistance. However, the roles of ETV7 in breast cancer inflammation have yet to be studied. Gene ontology analysis previously performed on BC cells stably over-expressing ETV7 demonstrated that ETV7 was involved in the suppression of innate immune and inflammatory responses. To better decipher the involvement of ETV7 in these signaling pathways, in this study, we identified TNFRSF1A, encoding for the main receptor of TNF-α, TNFR1, as one of the genes down-regulated by ETV7. We demonstrated that ETV7 directly binds to the intron I of this gene, and we showed that the ETV7-mediated down-regulation of TNFRSF1A reduced the activation of NF-κB signaling. Furthermore, in this study, we unveiled a potential crosstalk between ETV7 and STAT3, another master regulator of inflammation. While it is known that STAT3 directly up-regulates the expression of TNFRSF1A, here we demonstrated that ETV7 reduces the ability of STAT3 to bind to the TNFRSF1A gene via a competitive mechanism, recruiting repressive chromatin remodelers, which results in the repression of its transcription. The inverse correlation between ETV7 and TNFRSF1A was confirmed also in different cohorts of BC patients. These results suggest that ETV7 can reduce the inflammatory responses in breast cancer through the down-regulation of TNFRSF1A.

ETV3 and ETV6 enable monocyte differentiation into dendritic cells by repressing macrophage fate commitment

In inflamed tissues, monocytes differentiate into macrophages (mo-Macs) or dendritic cells (mo-DCs). In chronic nonresolving inflammation, mo-DCs are major drivers of pathogenic events. Manipulating monocyte differentiation would therefore be an attractive therapeutic strategy. However, how the balance of mo-DC versus mo-Mac fate commitment is regulated is not clear. In the present study, we show that the transcriptional repressors ETV3 and ETV6 control human monocyte differentiation into mo-DCs. ETV3 and ETV6 inhibit interferon (IFN)-stimulated genes; however, their action on monocyte differentiation is independent of IFN signaling. Instead, we find that ETV3 and ETV6 directly repress mo-Mac development by controlling MAFB expression. Mice deficient for Etv6 in monocytes have spontaneous expression of IFN-stimulated genes, confirming that Etv6 regulates IFN responses in vivo. Furthermore, these mice have impaired mo-DC differentiation during inflammation and reduced pathology in an experimental autoimmune encephalomyelitis model. These findings provide information about the molecular control of monocyte fate decision and identify ETV6 as a therapeutic target to redirect monocyte differentiation in inflammatory disorders.

Genome-wide transcriptional profiling of pulmonary functional sequelae in ARDS- secondary to SARS-CoV-2 infection

BACKGROUND

Up to 80% of patients surviving acute respiratory distress syndrome (ARDS) secondary to SARS-CoV-2 infection present persistent anomalies in pulmonary function after hospital discharge. There is a limited understanding of the mechanistic pathways linked to post-acute pulmonary sequelae.

AIM

To identify the molecular underpinnings associated with severe lung diffusion involvement in survivors of SARS-CoV-2-induced ARDS.

METHODS

Survivors attended to a complete pulmonary evaluation 3 months after hospital discharge. RNA sequencing (RNA-seq) was performed using Illumina technology in whole-blood samples from 50 patients with moderate to severe diffusion impairment (D_LCO<60%) and age- and sex-matched individuals with mild-normal lung function (D_LCO≥60%). A transcriptomic signature for optimal classification was constructed using random forest. Transcriptomic data were analyzed for biological pathway enrichment, cellular deconvolution, cell/tissue-specific gene expression and candidate drugs.

RESULTS

RNA-seq identified 1357 differentially expressed transcripts. A model composed of 14 mRNAs allowed the optimal discrimination of survivors with severe diffusion impairment (AUC=0.979). Hallmarks of lung sequelae involved cell death signaling, cytoskeleton reorganization, cell growth and differentiation and the immune response. Resting natural killer (NK) cells were the most important immune cell subtype for the prediction of severe diffusion impairment. Components of the signature correlated with neutrophil, lymphocyte and monocyte counts. A variable expression profile of the transcripts was observed in lung cell subtypes and bodily tissues. One upregulated gene, TUBB4A, constitutes a target for FDA-approved drugs.

CONCLUSIONS

This work defines the transcriptional programme associated with post-acute pulmonary sequelae and provides novel insights for targeted interventions and biomarker development.

In Vivo Profiling of Individual Multiciliated Cells during Acute Influenza A Virus Infection

Influenza virus infections are thought to be initiated in a small number of cells; however, the heterogeneity across the cellular responses of the epithelial cells during establishment of disease is incompletely understood. Here, we used an H1N1 influenza virus encoding a fluorescent reporter gene, a cell lineage-labeling transgenic mouse line, and single-cell RNA sequencing to explore the range of responses in a susceptible epithelial cell population during an acute influenza A virus (IAV) infection. Focusing on multiciliated cells, we identified a subpopulation that basally expresses interferon-stimulated genes (ISGs), which we hypothesize may be important for the early response to infection. We subsequently found that a population of infected ciliated cells produce most of the ciliated cell-derived inflammatory cytokines, and nearly all bystander ciliated cells induce a broadly antiviral state. From these data together, we propose that variable preexisting gene expression patterns in the initial cells targeted by the virus may ultimately affect the establishment of viral disease. IMPORTANCE Influenza A virus poses a significant threat to public health, and each year, millions of people in the United States alone are exposed to the virus. We do not currently, however, fully understand why some individuals clear the infection asymptomatically and others become severely ill. Understanding how these divergent phenotypes arise could eventually be leveraged to design therapeutics that prevent severe disease. As a first step toward understanding these different infection states, we used a technology that allowed us to determine how thousands of individual murine lung epithelial cells behaved before and during IAV infection. We found that small subsets of epithelial cells exhibited an antiviral state prior to infection, and similarly, some cells made high levels of inflammatory cytokines during infection. We propose that different ratios of these individual cellular responses may contribute to the broader antiviral state of the lung and may ultimately affect disease severity.

Genetic Associations and Differential mRNA Expression Levels of Host Genes Suggest a Viral Trigger for Endemic Pemphigus Foliaceus

The long search for the environmental trigger of the endemic pemphigus foliaceus (EPF, fogo selvagem) has not yet resulted in any tangible findings. Here, we searched for genetic associations and the differential expression of host genes involved in early viral infections and innate antiviral defense. Genetic variants could alter the structure, expression sites, or levels of the gene products, impacting their functions. By analyzing 3063 variants of 166 candidate genes in 227 EPF patients and 194 controls, we found 12 variants within 11 genes associated with differential susceptibility (p < 0.005) to EPF. The products of genes TRIM5, TPCN2, EIF4E, EIF4E3, NUP37, NUP50, NUP88, TPR, USP15, IRF8, and JAK1 are involved in different mechanisms of viral control, for example, the regulation of viral entry into the host cell or recognition of viral nucleic acids and proteins. Only two of nine variants were also associated in an independent German cohort of sporadic PF (75 patients, 150 controls), aligning with our hypothesis that antiviral host genes play a major role in EPF due to a specific virus−human interaction in the endemic region. Moreover, CCL5, P4HB, and APOBEC3G mRNA levels were increased (p < 0.001) in CD4+ T lymphocytes of EPF patients. Because there is limited or no evidence that these genes are involved in autoimmunity, their crucial role in antiviral responses and the associations that we observed support the hypothesis of a viral trigger for EPF, presumably a still unnoticed flavivirus. This work opens new frontiers in searching for the trigger of EPF, with the potential to advance translational research that aims for disease prevention and treatment.

Comparative Investigation of Gene Regulatory Processes Underlying Avian Influenza Viruses in Chicken and Duck

Simple Summary

Avian influenza poses a great risk to gallinaceous poultry, while mallard ducks can withstand most virus strains. To date, the mechanisms underlying the susceptibility of chicken and the effective immune response of duck have not been completely understood. In this study, our aim is to investigate the transcriptional gene regulation governing the expression of important avian-influenza-induced genes and to reveal the master regulators stimulating an effective immune response after virus infection in ducks while dysfunctioning in chicken.

Abstract

The avian influenza virus (AIV) mainly affects birds and not only causes animals’ deaths, but also poses a great risk of zoonotically infecting humans. While ducks and wild waterfowl are seen as a natural reservoir for AIVs and can withstand most virus strains, chicken mostly succumb to infection with high pathogenic avian influenza (HPAI). To date, the mechanisms underlying the susceptibility of chicken and the effective immune response of duck have not been completely unraveled. In this study, we investigate the transcriptional gene regulation underlying disease progression in chicken and duck after AIV infection. For this purpose, we use a publicly available RNA-sequencing dataset from chicken and ducks infected with low-pathogenic avian influenza (LPAI) H5N2 and HPAI H5N1 (lung and ileum tissues, 1 and 3 days post-infection). Unlike previous studies, we performed a promoter analysis based on orthologous genes to detect important transcription factors (TFs) and their cooperation, based on which we apply a systems biology approach to identify common and species-specific master regulators. We found master regulators such as EGR1, FOS, and SP1, specifically for chicken and ETS1 and SMAD3/4, specifically for duck, which could be responsible for the duck’s effective and the chicken’s ineffective immune response.

Comprehensive Analysis Identified ETV7 as a Potential Prognostic Biomarker in Bladder Cancer

Background

The tumor microenvironment (TME) plays a crucial role in the initiation and progression of cancer. Bladder cancer (BLCA) is a malignant tumor of the genitourinary system. Its heterogeneity results in significant differences in the prognosis of patients. To date, this is still a huge challenge for clinical treatment. In recent years, more and more evidence showed that dysregulation of transcription factors (TFs) plays an important role in tumor progression, invasion, and metastasis. Unfortunately, the role of TFs on the tumor microenvironment in bladder cancer is unclear.

Methods

The original data of BLCA and corresponding adjacent tissues were obtained from The Cancer Genome Atlas (TCGA) database. TFs were downloaded from the Animal Transcription Factor DataBase (Animal TFDB). Intersection analysis was used to obtain TFs that were differentially expressed between tumor and adjacent tissues. Gene Set Cancer Analysis (GSCALite) and CIBERSORT software were used to reveal the key differentially expressed TFs (DE-TFs). Subsequently, UALCAN and Human Protein Atlas (HPA) databases were used to disclose the expression of key DE-TFs in BLCA. The K-M curve divulged the relationship between the key DE-TFs and the patient's overall survival (OS), and the univariate and multivariate Cox regression analyses were conducted to explore independent prognostic factors. The cluster profiler package and Gene Set Enrichment Analysis (GSEA) were used for functional enrichment of genes related to the key DE-TFs. Finally, CIBERSORT software analyzed the immune landscape of BLCA.

Results

We obtained a total of 117 BLCA-related DE-TFs. Among them, ETV7 was identified as the key DE-TFs due to its association with the autophagy activation pathway and various immune cells in cancer. Online databases of UALCAN and HPA indicated that ETV7 was overexpressed in tumors and negatively correlated with tumor severity. The K-M curve showed that the OS of patients with high expression of ETV7 was poor, which indicated that it was an independent prognostic factor. Functional enrichment of 87 DEGs between ETV7-high and -low expression groups indicated that it was closely related to the immune response and the functions of a variety of immune cells. Finally, CIBERSORT results proved that the high and low expression of ETV7 also caused significant differences in the tumor immune microenvironment of patients.

Conclusion

Overall, we proved that the transcription factor ETV7 was a novel prognostic factor, which may improve the individualized outcome prediction in BLCA by regulating the tumor immune microenvironment.