A methylation clock model of mild SARS-CoV-2 infection provides insight into immune dysregulation.

DNA methylation comprises a cumulative record of lifetime exposures superimposed on genetically determined markers. Little is known about methylation dynamics in humans following an acute perturbation, such as infection. We characterized the temporal trajectory of blood epigenetic remodeling in 133 participants in a prospective study of young adults before, during, and after asymptomatic and mildly symptomatic SARS-CoV-2 infection. The differential methylation caused by asymptomatic or mildly symptomatic infections was indistinguishable. While differential gene expression largely returned to baseline levels after the virus became undetectable, some differentially methylated sites persisted for months of follow-up, with a pattern resembling autoimmune or inflammatory disease. We leveraged these responses to construct methylation-based machine learning models that distinguished samples from pre-, during-, and postinfection time periods, and quantitatively predicted the time since infection. The clinical trajectory in the young adults and in a diverse cohort with more severe outcomes was predicted by the similarity of methylation before or early after SARS-CoV-2 infection to the model-defined postinfection state. Unlike the phenomenon of trained immunity, the postacute SARS-CoV-2 epigenetic landscape we identify is antiprotective.

SARS-CoV-2 Transmission among Marine Recruits during Quarantine.

BACKGROUND: The efficacy of public health measures to control the transmission of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has not been well studied in young adults. METHODS: We investigated SARS-CoV-2 infections among U.S. Marine Corps recruits who underwent a 2-week quarantine at home followed by a second supervised 2-week quarantine at a closed college campus that involved mask wearing, social distancing, and daily temperature and symptom monitoring. Study volunteers were tested for SARS-CoV-2 by means of quantitative polymerase-chain-reaction (qPCR) assay of nares swab specimens obtained between the time of arrival and the second day of supervised quarantine and on days 7 and 14. Recruits who did not volunteer for the study underwent qPCR testing only on day 14, at the end of the quarantine period. We performed phylogenetic analysis of viral genomes obtained from infected study volunteers to identify clusters and to assess the epidemiologic features of infections. RESULTS: A total of 1848 recruits volunteered to participate in the study; within 2 days after arrival on campus, 16 (0.9%) tested positive for SARS-CoV-2, 15 of whom were asymptomatic. An additional 35 participants (1.9%) tested positive on day 7 or on day 14. Five of the 51 participants (9.8%) who tested positive at any time had symptoms in the week before a positive qPCR test. Of the recruits who declined to participate in the study, 26 (1.7%) of the 1554 recruits with available qPCR results tested positive on day 14. No SARS-CoV-2 infections were identified through clinical qPCR testing performed as a result of daily symptom monitoring. Analysis of 36 SARS-CoV-2 genomes obtained from 32 participants revealed six transmission clusters among 18 participants. Epidemiologic analysis supported multiple local transmission events, including transmission between roommates and among recruits within the same platoon. CONCLUSIONS: Among Marine Corps recruits, approximately 2% who had previously had negative results for SARS-CoV-2 at the beginning of supervised quarantine, and less than 2% of recruits with unknown previous status, tested positive by day 14. Most recruits who tested positive were asymptomatic, and no infections were detected through daily symptom monitoring. Transmission clusters occurred within platoons. (Funded by the Defense Health Agency and others.).

Skeletal muscle transcriptome response to a bout of endurance exercise in physically active and sedentary older adults.

Age-related declines in cardiorespiratory fitness and physical function are mitigated by regular endurance exercise in older adults. This may be due, in part, to changes in the transcriptional program of skeletal muscle following repeated bouts of exercise. However, the impact of chronic exercise training on the transcriptional response to an acute bout of endurance exercise has not been clearly determined. Here, we characterized baseline differences in muscle transcriptome and exercise-induced response in older adults who were active/endurance trained or sedentary. RNA-sequencing was performed on vastus lateralis biopsy specimens obtained before, immediately after, and 3 h following a bout of endurance exercise (40 min of cycling at 60%-70% of heart rate reserve). Using a recently developed bioinformatics approach, we found that transcript signatures related to type I myofibers, mitochondria, and endothelial cells were higher in active/endurance-trained adults and were associated with key phenotypic features including V̇o2peak, ATPmax, and muscle fiber proportion. Immune cell signatures were elevated in the sedentary group and linked to visceral and intermuscular adipose tissue mass. Following acute exercise, we observed distinct temporal transcriptional signatures that were largely similar among groups. Enrichment analysis revealed catabolic processes were uniquely enriched in the sedentary group at the 3-h postexercise timepoint. In summary, this study revealed key transcriptional signatures that distinguished active and sedentary adults, which were associated with difference in oxidative capacity and depot-specific adiposity. The acute response signatures were consistent with beneficial effects of endurance exercise to improve muscle health in older adults irrespective of exercise history and adiposity.NEW & NOTEWORTHY Muscle transcript signatures associated with oxidative capacity and immune cells underlie important phenotypic and clinical characteristics of older adults who are endurance trained or sedentary. Despite divergent phenotypes, the temporal transcriptional signatures in response to an acute bout of endurance exercise were largely similar among groups. These data provide new insight into the transcriptional programs of aging muscle and the beneficial effects of endurance exercise to promote healthy aging in older adults.

SARS-CoV-2 Outbreak Dynamics in an Isolated US Military Recruit Training Center With Rigorous Prevention Measures.

BACKGROUND: Marine recruits training at Parris Island experienced an unexpectedly high rate of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, despite preventive measures including a supervised, 2-week, pre-entry quarantine. We characterize SARS-CoV-2 transmission in this cohort. METHODS: Between May and November 2020, we monitored 2,469 unvaccinated, mostly male, Marine recruits prospectively during basic training. If participants tested negative for SARS-CoV-2 by quantitative polymerase chain reaction (qPCR) at the end of quarantine, they were transferred to the training site in segregated companies and underwent biweekly testing for 6 weeks. We assessed the effects of coronavirus disease 2019 (COVID-19) prevention measures on other respiratory infections with passive surveillance data, performed phylogenetic analysis, and modeled transmission dynamics and testing regimens. RESULTS: Preventive measures were associated with drastically lower rates of other respiratory illnesses. However, among the trainees, 1,107 (44.8%) tested SARS-CoV-2-positive, with either mild or no symptoms. Phylogenetic analysis of viral genomes from 580 participants revealed that all cases but one were linked to five independent introductions, each characterized by accumulation of mutations across and within companies, and similar viral isolates in individuals from the same company. Variation in company transmission rates (mean reproduction number R 0 ; 5.5 [95% confidence interval [CI], 5.0, 6.1]) could be accounted for by multiple initial cases within a company and superspreader events. Simulations indicate that frequent rapid-report testing with case isolation may minimize outbreaks. CONCLUSIONS: Transmission of wild-type SARS-CoV-2 among Marine recruits was approximately twice that seen in the community. Insights from SARS-CoV-2 outbreak dynamics and mutations spread in a remote, congregate setting may inform effective mitigation strategies.

Differential analysis of chromatin accessibility and gene expression profiles identifies cis-regulatory elements in rat adipose and muscle.

Chromatin accessibility is a key factor influencing gene expression. We optimized the Omni-ATAC-seq protocol and used it together with RNA-seq to investigate cis-regulatory elements in rat white adipose and skeletal muscle, two tissues with contrasting metabolic functions. While promoter accessibility correlated with RNA expression, integration of the two datasets identified tissue-specific differentially accessible regions (DARs) that predominantly localized in intergenic and intron regions. DARs were mapped to differentially expressed (DE) genes enriched in distinct biological processes in each tissue. Randomly selected DE genes were validated by qPCR. Top enriched motifs in DARs predicted binding sites for transcription factors (TFs) showing tissue-specific up-regulation. The correlation between differential chromatin accessibility at a given TF binding motif and differential expression of target genes further supported the functional relevance of that motif. Our study identified cis-regulatory regions that likely play a major role in the regulation of tissue-specific gene expression in adipose and muscle.

Innate Immune Response to Influenza Virus at Single-Cell Resolution in Human Epithelial Cells Revealed Paracrine Induction of Interferon Lambda 1.

Early interactions of influenza A virus (IAV) with respiratory epithelium might determine the outcome of infection. The study of global cellular innate immune responses often masks multiple aspects of the mechanisms by which populations of cells work as organized and heterogeneous systems to defeat virus infection, and how the virus counteracts these systems. In this study, we experimentally dissected the dynamics of IAV and human epithelial respiratory cell interaction during early infection at the single-cell level. We found that the number of viruses infecting a cell (multiplicity of infection [MOI]) influences the magnitude of virus antagonism of the host innate antiviral response. Infections performed at high MOIs resulted in increased viral gene expression per cell and stronger antagonist effect than infections at low MOIs. In addition, single-cell patterns of expression of interferons (IFN) and IFN-stimulated genes (ISGs) provided important insights into the contributions of the infected and bystander cells to the innate immune responses during infection. Specifically, the expression of multiple ISGs was lower in infected than in bystander cells. In contrast with other IFNs, IFN lambda 1 (IFNL1) showed a widespread pattern of expression, suggesting a different cell-to-cell propagation mechanism more reliant on paracrine signaling. Finally, we measured the dynamics of the antiviral response in primary human epithelial cells, which highlighted the importance of early innate immune responses at inhibiting virus spread.IMPORTANCE Influenza A virus (IAV) is a respiratory pathogen of high importance to public health. Annual epidemics of seasonal IAV infections in humans are a significant public health and economic burden. IAV also causes sporadic pandemics, which can have devastating effects. The main target cells for IAV replication are epithelial cells in the respiratory epithelium. The cellular innate immune responses induced in these cells upon infection are critical for defense against the virus, and therefore, it is important to understand the complex interactions between the virus and the host cells. In this study, we investigated the innate immune response to IAV in the respiratory epithelium at the single-cell level, providing a better understanding on how a population of epithelial cells functions as a complex system to orchestrate the response to virus infection and how the virus counteracts this system.

Combination of Gene Expression Signature and Model for End-Stage Liver Disease Score Predicts Survival of Patients With Severe Alcoholic Hepatitis.

BACKGROUND & AIMS: Patients with severe alcoholic hepatitis (AH) have a high risk of death within 90 days. Corticosteroids, which can cause severe adverse events, are the only treatment that increases short-term survival. It is a challenge to predict outcomes of patients with severe AH. Therefore, we developed a scoring system to predict patient survival, integrating baseline molecular and clinical variables. METHODS: We obtained fixed liver biopsy samples from 71 consecutive patients diagnosed with severe AH and treated with corticosteroids from July 2006 through December 2013 in Brussels, Belgium (derivation cohort). Gene expression patterns were analyzed by microarrays and clinical data were collected for 180 days. We identified gene expression signatures and clinical data that are associated with survival without liver transplantation at 90 and 180 days after initiation of corticosteroid therapy. Findings were validated using liver biopsies from 48 consecutive patients with severe AH treated with corticosteroids, collected from March 2010 through February 2015 at hospitals in Belgium and Switzerland (validation cohort 1) and in liver biopsies from 20 patients (9 received corticosteroid treatment), collected from January 2012 through May 2015 in the United States (validation cohort 2). RESULTS: We integrated data on expression patterns of 123 genes and the model for end-stage liver disease (MELD) scores to assign patients to groups with poor survival (29% survived 90 days and 26% survived 180 days) and good survival (76% survived 90 days and 65% survived 180 days) (P < .001) in the derivation cohort. We named this assignment system the gene signature-MELD (gs-MELD) score. In validation cohort 1, the gs-MELD score discriminated patients with poor survival (43% survived 90 days) from those with good survival (96% survived 90 days) (P < .001). The gs-MELD score also discriminated between patients with a poor survival at 180 days (34% survived) and a good survival at 180 days (84% survived) (P < .001). The time-dependent area under the receiver operator characteristic curve for the score was 0.86 (95% confidence interval 0.73-0.99) for survival at 90 days, and 0.83 (95% confidence interval 0.71-0.96) for survival at 180 days. This score outperformed other clinical models to predict survival of patients with severe AH in validation cohort 1. In validation cohort 2, the gs-MELD discriminated patients with a poor survival at 90 days (12% survived) from those with a good survival at 90 days (100%) (P < .001). CONCLUSIONS: We integrated data on baseline liver gene expression pattern and the MELD score to create the gs-MELD scoring system, which identifies patients with severe AH, treated or not with corticosteroids, most and least likely to survive for 90 and 180 days.

Low-variance RNAs identify Parkinson's disease molecular signature in blood.

The diagnosis of Parkinson's disease (PD) is usually not established until advanced neurodegeneration leads to clinically detectable symptoms. Previous blood PD transcriptome studies show low concordance, possibly resulting from the use of microarray technology, which has high measurement variation. The Leucine-rich repeat kinase 2 (LRRK2) G2019S mutation predisposes to PD. Using preclinical and clinical studies, we sought to develop a novel statistically motivated transcriptomic-based approach to identify a molecular signature in the blood of Ashkenazi Jewish PD patients, including LRRK2 mutation carriers. Using a digital gene expression platform to quantify 175 messenger RNA (mRNA) markers with low coefficients of variation (CV), we first compared whole-blood transcript levels in mouse models (1) overexpressing wild-type (WT) LRRK2, (2) overexpressing G2019S LRRK2, (3) lacking LRRK2 (knockout), and (4) and in WT controls. We then studied an Ashkenazi Jewish cohort of 34 symptomatic PD patients (both WT LRRK2 and G2019S LRRK2) and 32 asymptomatic controls. The expression profiles distinguished the four mouse groups with different genetic background. In patients, we detected significant differences in blood transcript levels both between individuals differing in LRRK2 genotype and between PD patients and controls. Discriminatory PD markers included genes associated with innate and adaptive immunity and inflammatory disease. Notably, gene expression patterns in levodopa-treated PD patients were significantly closer to those of healthy controls in a dose-dependent manner. We identify whole-blood mRNA signatures correlating with LRRK2 genotype and with PD disease state. This approach may provide insight into pathogenesis and a route to early disease detection.

Protocol for high-throughput DNA methylation profiling in rat tissues using automated reduced representation bisulfite sequencing.

Although reduced representation bisulfite sequencing (RRBS) measures DNA methylation (DNAme) across CpG-rich genomic regions with high sensitivity, the assay can be time-consuming and prone to batch effects. Here, we present a high-throughput, automated RRBS protocol starting with DNA extraction from frozen rat tissues. We describe steps for RRBS library preparation, library quality control, and sequencing. We also detail an optimized pipeline for sequencing data processing. This protocol has been applied successfully to DNAme profiling across multiple rat tissues. For complete details on the use and execution of this protocol, please refer to Nair et al.1.

Molecular adaptations in response to exercise training are associated with tissue-specific transcriptomic and epigenomic signatures.

Regular exercise has many physical and brain health benefits, yet the molecular mechanisms mediating exercise effects across tissues remain poorly understood. Here we analyzed 400 high-quality DNA methylation, ATAC-seq, and RNA-seq datasets from eight tissues from control and endurance exercise-trained (EET) rats. Integration of baseline datasets mapped the gene location dependence of epigenetic control features and identified differing regulatory landscapes in each tissue. The transcriptional responses to 8 weeks of EET showed little overlap across tissues and predominantly comprised tissue-type enriched genes. We identified sex differences in the transcriptomic and epigenomic changes induced by EET. However, the sex-biased gene responses were linked to shared signaling pathways. We found that many G protein-coupled receptor-encoding genes are regulated by EET, suggesting a role for these receptors in mediating the molecular adaptations to training across tissues. Our findings provide new insights into the mechanisms underlying EET-induced health benefits across organs.

Localization of BRCA1 protein in breast cancer tissue and cell lines with mutations.

BACKGROUND: The breast and ovarian cancer susceptibility gene (BRCA1) encodes a tumor suppressor. The BRCA1 protein is found primarily in cell nuclei and plays an important role in the DNA damage response and transcriptional regulation. Deficiencies in DNA repair capabilities have been associated with higher histopathological grade and worse prognosis in breast cancer. METHODS: In order to investigate the subcellular distribution of BRCA1 in tumor tissue we randomly selected 22 breast carcinomas and tested BRCA1 protein localization in frozen and contiguous formalin-fixed, paraffin embedded (FFPE) tissue, using pressure cooker antigen-retrieval and the MS110 antibody staining. To assess the impact of BRCA1 germline mutations on protein localization, we retrospectively tested 16 of the tumor specimens to determine whether they contained the common Ashkenazi Jewish founder mutations in BRCA1 (185delAG, 5382insC), and BRCA2 (6174delT). We also compared co-localization of BRCA1 and nucleolin in MCF7 cells (wild type) and a mutant BRCA1 cell line, HCC1937 (5382insC). RESULTS: In FFPE tissue, with MS110 antibody staining, we frequently found reduced BRCA1 nuclear staining in breast tumor tissue compared to normal tissue, and less BRCA1 staining with higher histological grade in the tumors. However, in the frozen sections, BRCA1 antibody staining showed punctate, intra-nuclear granules in varying numbers of tumor, lactating, and normal cells. Two mutation carriers were identified and were confirmed by gene sequencing. We have also compared co-localization of BRCA1 and nucleolin in MCF7 cells (wild type) and a mutant BRCA1 cell line, HCC1937 (5382insC) and found altered sub-nuclear and nucleolar localization patterns consistent with a functional impact of the mutation on protein localization. CONCLUSIONS: The data presented here support a role for BRCA1 in the pathogenesis of sporadic and inherited breast cancers. The use of well-characterized reagents may lead to further insights into the function of BRCA1 and possibly the further development of targeted therapeutics.

Activating Transcription Factor 3 Stimulates Follicle-Stimulating Hormone-ß Expression In Vitro But Is Dispensable for Follicle-Stimulating Hormone Production in Murine Gonadotropes In Vivo.

Follicle-stimulating hormone (FSH), a dimeric glycoprotein produced by pituitary gonadotrope cells, regulates spermatogenesis in males and ovarian follicle growth in females. Hypothalamic gonadotropin-releasing hormone (GnRH) stimulates FSHß subunit gene (Fshb) transcription, though the underlying mechanisms are poorly understood. To address this gap in knowledge, we examined changes in pituitary gene expression in GnRH-deficient mice (hpg) treated with a regimen of exogenous GnRH that increases pituitary Fshb but not luteinizing hormone ß (Lhb) messenger RNA levels. Activating transcription factor 3 (Atf3) was among the most upregulated genes. Activating transcription factor 3 (ATF3) can heterodimerize with members of the activator protein 1 family to regulate gene transcription. Co-expression of ATF3 with JunB stimulated murine Fshb, but not Lhb, promoter-reporter activity in homologous LßT2b cells. ATF3 also synergized with a constitutively active activin type I receptor to increase endogenous Fshb expression in these cells. Nevertheless, FSH production was intact in gonadotrope-specific Atf3 knockout [conditional knockout (cKO)] mice. Ovarian follicle development, ovulation, and litter sizes were equivalent between cKOs and controls. Testis weights and sperm counts did not differ between genotypes. Following gonadectomy, increases in LH secretion were enhanced in cKO animals. Though FSH levels did not differ between genotypes, post-gonadectomy increases in pituitary Fshb and gonadotropin α subunit expression were more pronounced in cKO than control mice. These data indicate that ATF3 can selectively stimulate Fshb expression in vitro but is not required for FSH production in vivo.

Integrated single-cell multiome analysis reveals muscle fiber-type gene regulatory circuitry modulated by endurance exercise.

Endurance exercise is an important health modifier. We studied cell-type specific adaptations of human skeletal muscle to acute endurance exercise using single-nucleus (sn) multiome sequencing in human vastus lateralis samples collected before and 3.5 hours after 40 min exercise at 70% VO2max in four subjects, as well as in matched time of day samples from two supine resting circadian controls. High quality same-cell RNA-seq and ATAC-seq data were obtained from 37,154 nuclei comprising 14 cell types. Among muscle fiber types, both shared and fiber-type specific regulatory programs were identified. Single-cell circuit analysis identified distinct adaptations in fast, slow and intermediate fibers as well as LUM-expressing FAP cells, involving a total of 328 transcription factors (TFs) acting at altered accessibility sites regulating 2,025 genes. These data and circuit mapping provide single-cell insight into the processes underlying tissue and metabolic remodeling responses to exercise.

Multi-omic identification of key transcriptional regulatory programs during endurance exercise training.

Transcription factors (TFs) play a key role in regulating gene expression and responses to stimuli. We conducted an integrated analysis of chromatin accessibility, DNA methylation, and RNA expression across eight rat tissues following endurance exercise training (EET) to map epigenomic changes to transcriptional changes and determine key TFs involved. We uncovered tissue-specific changes and TF motif enrichment across all omic layers, differentially accessible regions (DARs), differentially methylated regions (DMRs), and differentially expressed genes (DEGs). We discovered distinct routes of EET-induced regulation through either epigenomic alterations providing better access for TFs to affect target genes, or via changes in TF expression or activity enabling target gene response. We identified TF motifs enriched among correlated epigenomic and transcriptomic alterations, DEGs correlated with exercise-related phenotypic changes, and EET-induced activity changes of TFs enriched for DEGs among their gene targets. This analysis elucidates the unique transcriptional regulatory mechanisms mediating diverse organ effects of EET.

Blood RNA alternative splicing events as diagnostic biomarkers for infectious disease.

Assays detecting blood transcriptome changes are studied for infectious disease diagnosis. Blood-based RNA alternative splicing (AS) events, which have not been well characterized in pathogen infection, have potential normalization and assay platform stability advantages over gene expression for diagnosis. Here, we present a computational framework for developing AS diagnostic biomarkers. Leveraging a large prospective cohort of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection and whole-blood RNA sequencing (RNA-seq) data, we identify a major functional AS program switch upon viral infection. Using an independent cohort, we demonstrate the improved accuracy of AS biomarkers for SARS-CoV-2 diagnosis compared with six reported transcriptome signatures. We then optimize a subset of AS-based biomarkers to develop microfluidic PCR diagnostic assays. This assay achieves nearly perfect test accuracy (61/62 = 98.4%) using a naive principal component classifier, significantly more accurate than a gene expression PCR assay in the same cohort. Therefore, our RNA splicing computational framework enables a promising avenue for host-response diagnosis of infection.

Mapping disease regulatory circuits at cell-type resolution from single-cell multiomics data.

Resolving chromatin-remodeling-linked gene expression changes at cell-type resolution is important for understanding disease states. Here we describe MAGICAL (Multiome Accessibility Gene Integration Calling and Looping), a hierarchical Bayesian approach that leverages paired single-cell RNA sequencing and single-cell transposase-accessible chromatin sequencing from different conditions to map disease-associated transcription factors, chromatin sites, and genes as regulatory circuits. By simultaneously modeling signal variation across cells and conditions in both omics data types, MAGICAL achieved high accuracy on circuit inference. We applied MAGICAL to study Staphylococcus aureus sepsis from peripheral blood mononuclear single-cell data that we generated from subjects with bloodstream infection and uninfected controls. MAGICAL identified sepsis-associated regulatory circuits predominantly in CD14 monocytes, known to be activated by bacterial sepsis. We addressed the challenging problem of distinguishing host regulatory circuit responses to methicillin-resistant and methicillin-susceptible S. aureus infections. Although differential expression analysis failed to show predictive value, MAGICAL identified epigenetic circuit biomarkers that distinguished methicillin-resistant from methicillin-susceptible S. aureus infections.

Optimization of the Omni-ATAC protocol to chromatin accessibility profiling in snap-frozen rat adipose and muscle tissues.

ATAC-seq is a fast and sensitive method for the epigenomic profiling of open chromatin and for mapping of transcription factor binding sites [1]. Despite the development of the Omni-ATAC protocol for the profiling of chromatin accessibility in frozen tissues [2], studies in adipose tissue have been restricted due to technical challenges including the high lipid content of adipocytes and reproducibility issues between replicates. Here, we provide a modified Omni-ATAC protocol that achieves high data reproducibility in various tissue types from rat, including adipose and muscle tissues [3].•This protocol describes a methodology that enables chromatin accessibility profiling from snap-frozen rat adipose and muscle tissues.•The technique comprises an optimized bead-based tissue homogenization process that substitutes to Dounce homogenization, reduces variability in the experimental procedure, and is adaptable to various tissue types.•In comparison with the Omni-ATAC protocol, the method described here results in improved ATAC-seq data quality that complies with ENCODE quality standards.

Multi-omics profiling of single nuclei from frozen archived postmortem human pituitary tissue.

Concomitant profiling of transcriptome and chromatin accessibility in isolated nuclei can reveal gene regulatory control mechanisms in health and disease. We report a single nucleus multi-omics analysis protocol optimized for frozen archived postmortem human pituitaries that is also effective for frozen ovine and murine pituitaries and human skeletal muscle biopsies. Its main advantages are that (1) it is not limited to fresh tissue, (2) it avoids tissue dissociation-induced transcriptional changes, and (3) it includes a novel, automated quality control pipeline. For complete details on the use and execution of this protocol, please refer to Ruf-Zamojski et al. (2021) and Zhang et al. (2022).

Transcriptome signatures preceding the induction of anti-stalk antibodies elicited after universal influenza vaccination.

A phase 1 clinical trial to test the immunogenicity of a chimeric group 1 HA (cHA) universal influenza virus vaccine targeting the conserved stalk domain of the hemagglutinin of influenza viruses was carried out. Vaccination with adjuvanted-inactivated vaccines induced high anti-stalk antibody titers. We sought to identify gene expression signatures that correlate with such induction. Messenger-RNA sequencing in whole blood was performed on the peripheral blood of 53 vaccinees. We generated longitudinal data on the peripheral blood of 53 volunteers, at early (days 3 and 7) and late (28 days) time points after priming and boosting with cHAs. Differentially expressed gene analysis showed no differences between placebo and live-attenuated vaccine groups. However, an upregulation of genes involved in innate immune responses and type I interferon signaling was found at day 3 after vaccination with inactivated adjuvanted formulations. Cell type deconvolution analysis revealed a significant enrichment for monocyte markers and different subsets of dendritic cells as mediators for optimal B cell responses and significant increase of anti-stalk antibodies in sera. A significant upregulation of immunoglobulin-related genes was only observed after administration of adjuvanted vaccines (either as primer or booster) with specific induction of anti-stalk IGVH1-69. This approach informed of specific immune signatures that correlate with robust anti-stalk antibody responses, while also helping to understand the regulation of gene expression induced by cHA proteins under different vaccine regimens.

Asymptomatic SARS-CoV-2 Infection Is Associated With Higher Levels of Serum IL-17C, Matrix Metalloproteinase 10 and Fibroblast Growth Factors Than Mild Symptomatic COVID-19.

Young adults infected with SARS-CoV-2 are frequently asymptomatic or develop only mild disease. Because capturing representative mild and asymptomatic cases require active surveillance, they are less characterized than moderate or severe cases of COVID-19. However, a better understanding of SARS-CoV-2 asymptomatic infections might shed light into the immune mechanisms associated with the control of symptoms and protection. To this aim, we have determined the temporal dynamics of the humoral immune response, as well as the serum inflammatory profile, of mild and asymptomatic SARS-CoV-2 infections in a cohort of 172 initially seronegative prospectively studied United States Marine recruits, 149 of whom were subsequently found to be SARS-CoV-2 infected. The participants had blood samples taken, symptoms surveyed and PCR tests for SARS-CoV-2 performed periodically for up to 105 days. We found similar dynamics in the profiles of viral load and in the generation of specific antibody responses in asymptomatic and mild symptomatic participants. A proteomic analysis using an inflammatory panel including 92 analytes revealed a pattern of three temporal waves of inflammatory and immunoregulatory mediators, and a return to baseline for most of the inflammatory markers by 35 days post-infection. We found that 23 analytes were significantly higher in those participants that reported symptoms at the time of the first positive SARS-CoV-2 PCR compared with asymptomatic participants, including mostly chemokines and cytokines associated with inflammatory response or immune activation (i.e., TNF-α, TNF-ß, CXCL10, IL-8). Notably, we detected 7 analytes (IL-17C, MMP-10, FGF-19, FGF-21, FGF-23, CXCL5 and CCL23) that were higher in asymptomatic participants than in participants with symptoms; these are known to be involved in tissue repair and may be related to the control of symptoms. Overall, we found a serum proteomic signature that differentiates asymptomatic and mild symptomatic infections in young adults, including potential targets for developing new therapies and prognostic tests.