nf-core/airrflow: an adaptive immune receptor repertoire analysis workflow employing the Immcantation framework.

Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) is a valuable experimental tool to study the immune state in health and following immune challenges such as infectious diseases, (auto)immune diseases, and cancer. Several tools have been developed to reconstruct B cell and T cell receptor sequences from AIRR-seq data and infer B and T cell clonal relationships. However, currently available tools offer limited parallelization across samples, scalability or portability to high-performance computing infrastructures. To address this need, we developed nf-core/airrflow, an end-to-end bulk and single-cell AIRR-seq processing workflow which integrates the Immcantation Framework following BCR and TCR sequencing data analysis best practices. The Immcantation Framework is a comprehensive toolset, which allows the processing of bulk and single-cell AIRR-seq data from raw read processing to clonal inference. nf-core/airrflow is written in Nextflow and is part of the nf-core project, which collects community contributed and curated Nextflow workflows for a wide variety of analysis tasks. We assessed the performance of nf-core/airrflow on simulated sequencing data with sequencing errors and show example results with real datasets. To demonstrate the applicability of nf-core/airrflow to the high-throughput processing of large AIRR-seq datasets, we validated and extended previously reported findings of convergent antibody responses to SARS-CoV-2 by analyzing 97 COVID-19 infected individuals and 99 healthy controls, including a mixture of bulk and single-cell sequencing datasets. Using this dataset, we extended the convergence findings to 20 additional subjects, highlighting the applicability of nf-core/airrflow to validate findings in small in-house cohorts with reanalysis of large publicly available AIRR datasets. nf-core/airrflow is available free of charge, under the MIT license on GitHub (https://github.com/nf-core/airrflow). Detailed documentation and example results are available on the nf-core website at (https://nf-co.re/airrflow).

Antibody gene features associated with binding and functional activity in vaccine-derived human mAbs targeting malaria parasites.

Adjuvants have been essential to malaria vaccine development, but their impact on the vaccine-induced antibody repertoire is poorly understood. Here, we used cDNA sequences from antigen-specific single memory B cells to express 132 recombinant human anti-Pfs230 monoclonal antibodies (mAbs). Alhydrogel®-induced mAbs demonstrated higher binding to Pfs230D1, although functional activity was similar between adjuvants. All Alhydrogel® mAbs using IGHV1-69 gene bound to recombinant Pfs230D1, but none blocked parasite transmission to mosquitoes; similarly, no AS01 mAb using IGHV1-69 blocked transmission. Functional mAbs from both Alhydrogel® and AS01 vaccines used IGHV3-21 and IGHV3-30 genes. Antibodies with the longest CDR3 sequences were associated with binding but not functional activity. This study assesses adjuvant effects on antibody clonotype diversity during malaria vaccination.

Mpox vaccine and infection-driven human immune signatures: an immunological analysis of an observational study.

BACKGROUND: Monkeypox virus has recently infected more than 88 000 people, raising concerns about our preparedness against this emerging viral pathogen. Licensed and approved for mpox, the JYNNEOS vaccine has fewer side-effects than previous smallpox vaccines and has shown immunogenicity against monkeypox in animal models. This study aims to elucidate human immune responses to JYNNEOS vaccination compared with mpox-induced immunity. METHODS: Peripheral blood mononuclear cells and sera were obtained from ten individuals vaccinated with one or two doses of JYNNEOS and six individuals diagnosed with monkeypox virus infection. Samples were obtained from seven individuals before vaccination to serve as a baseline. We examined the polyclonal serum (ELISA) and single B-cell (heavy chain gene and transcriptome data) antibody repertoires and T-cell responses (activation-induced marker and intracellular cytokine staining assays) induced by the JYNNEOS vaccine versus monkeypox virus infection. FINDINGS: All participants were men between the ages of 21 and 60 years, except for one woman in the group of mpox-convalescent individuals, and none had previous orthopoxvirus exposure. All mpox cases were mild. Vaccinee samples were collected 6-33 days after the first dose and 5-40 days after the second dose. Mpox-convalescent samples were collected 20-102 days after infection. In vaccine recipients, gene-level plasmablast and antibody responses were negligible and sera displayed moderate binding to recombinant orthopoxviral proteins (A29L, A35R, E8L, A30L, A27L, A33R, B18R, and L1R) and native proteins from the 2022 monkeypox outbreak strain. By contrast, recent monkeypox virus infection (within 20-102 days) induced robust serum antibody responses to monkeypox virus proteins and to native monkeypox virus proteins from a viral isolate obtained during the 2022 outbreak. JYNNEOS vaccine recipients presented robust orthopoxviral CD4+ and CD8+ T-cell responses. INTERPRETATION: Infection with monkeypox virus resulted in robust B-cell and T-cell responses, whereas immunisation with JYNNEOS elicited more robust T-cell responses. These data can help to inform vaccine design and policies for preventing mpox in humans. FUNDING: National Cancer Institute (National Institutes of Health), National Institute of Allergy and Infectious Diseases (National Institutes of Health), and Icahn School of Medicine.

Mpox vaccine and infection-driven human immune signatures.

Background: Mpox (formerly known as monkeypox) outbreaks outside endemic areas peaked in July 2022, infecting > 85,000 people and raising concerns about our preparedness against this emerging viral pathogen. Licensed and approved for mpox, the JYNNEOS vaccine has fewer side effects than previous smallpox vaccines and demonstrated efficacy against mpox infection in humans. Comparing JYNNEOS vaccine- and mpox-induced immunity is imperative to evaluate JYNNEOS' immunogenicity and inform vaccine administration and design. Methods: We examined the polyclonal serum (ELISA) and single B cell (heavy chain gene and transcriptome data) antibody repertoires and T cells (AIM and ICS assays) induced by the JYNNEOS vaccine as well as mpox infection. Findings: Gene-level plasmablast and antibody responses were negligible and JYNNEOS vaccinee sera displayed minimal binding to recombinant mpox proteins and native proteins from the 2022 outbreak strain. In contrast, recent mpox infection (within 20-102 days) induced robust serum antibody responses to A29L, A35R, A33R, B18R, and A30L, and to native mpox proteins, compared to vaccinees. JYNNEOS vaccine recipients presented comparable CD4 and CD8 T cell responses against orthopox peptides to those observed after mpox infection. Interpretation: JYNNEOS immunization does not elicit a robust B cell response, and its immunogenicity may be mediated by T cells. Funding: Research reported in this publication was supported, in part, by the National Cancer Institute of the National Institutes of Health under Award Number U54CA267776, U19AI168631(VS), as well as institutional funds from the Icahn School of Medicine.

Adaptive Immune Receptor Repertoire (AIRR) Community Guide to TR and IG Gene Annotation.

High-throughput sequencing of adaptive immune receptor repertoires (AIRR, i.e., IG and TR) has revolutionized the ability to carry out large-scale experiments to study the adaptive immune response. Since the method was first introduced in 2009, AIRR sequencing (AIRR-Seq) has been applied to survey the immune state of individuals, identify antigen-specific or immune-state-associated signatures of immune responses, study the development of the antibody immune response, and guide the development of vaccines and antibody therapies. Recent advancements in the technology include sequencing at the single-cell level and in parallel with gene expression, which allows the introduction of multi-omics approaches to understand in detail the adaptive immune response. Analyzing AIRR-seq data can prove challenging even with high-quality sequencing, in part due to the many steps involved and the need to parameterize each step. In this chapter, we outline key factors to consider when preprocessing raw AIRR-Seq data and annotating the genetic origins of the rearranged receptors. We also highlight a number of common difficulties with common AIRR-seq data processing and provide strategies to address them.

Adaptive Immune Receptor Repertoire (AIRR) Community Guide to Repertoire Analysis.

Adaptive immune receptor repertoires (AIRRs) are rich with information that can be mined for insights into the workings of the immune system. Gene usage, CDR3 properties, clonal lineage structure, and sequence diversity are all capable of revealing the dynamic immune response to perturbation by disease, vaccination, or other interventions. Here we focus on a conceptual introduction to the many aspects of repertoire analysis and orient the reader toward the uses and advantages of each. Along the way, we note some of the many software tools that have been developed for these investigations and link the ideas discussed to chapters on methods provided elsewhere in this volume.

Bulk Sequencing from mRNA with UMI for Evaluation of B-Cell Isotype and Clonal Evolution: A Method by the AIRR Community.

During the course of an immune response to a virus such as influenza, B cells undergo activation, clonal expansion, isotype switching, and somatic hypermutation (SHM). Members of an antigen-experienced B-cell clone can have different sequence features including SHM in the immunoglobulin heavy-chain V (IGHV) gene and can use the same IGVH gene in combination with different constant regions or isotypes (e.g., IgM, IgG, IgA). To study these features of expanded clones in an immune response by AIRR-seq, we provide a bulk RNA-based sequencing experimental procedure with unique molecular identifiers (UMIs) and the accompanying bioinformatics analytical workflow.

Muscular carnosine is a marker for cardiorespiratory fitness and cardiometabolic risk factors in men with type 1 diabetes.

PURPOSE: Muscle is an essential organ for glucose metabolism and can be influenced by metabolic disorders and physical activity. Elevated muscle carnosine levels have been associated with insulin resistance and cardiometabolic risk factors. Little is known about muscle carnosine in type 1 diabetes (T1D) and how it is influenced by physical activity. The aim of this study was to characterize muscle carnosine in vivo by proton magnetic resonance spectroscopy (1H MRS) and evaluate the relationship with physical activity, clinical characteristics and lipoprotein subfractions. METHODS: 16 men with T1D (10 athletes/6 sedentary) and 14 controls without diabetes (9/5) were included. Body composition by DXA, cardiorespiratory capacity (VO2peak) and serum lipoprotein profile by proton nuclear magnetic resonance (1H NMR) were obtained. Muscle carnosine scaled to water (carnosineW) and to creatine (carnosineCR), creatine and intramyocellular lipids (IMCL) were quantified in vivo using 1H MRS in a 3T MR scanner in soleus muscle. RESULTS: Subjects with T1D presented higher carnosine CR levels compared to controls. T1D patients with a lower VO2peak presented higher carnosineCR levels compared to sedentary controls, but both T1D and control groups presented similar levels of carnosineCR at high VO2peak levels. CarnosineW followed the same trend. Integrated correlation networks in T1D demonstrated that carnosineW and carnosineCR were associated with cardiometabolic risk factors including total and abdominal fat, pro-atherogenic lipoproteins (very low-density lipoprotein subfractions), low VO2peak, and IMCL. CONCLUSIONS: Elevated muscle carnosine levels in persons with T1D and their effect on atherogenic lipoproteins can be modulated by physical activity.

Benchmarking immunoinformatic tools for the analysis of antibody repertoire sequences.

SUMMARY: Antibody repertoires reveal insights into the biology of the adaptive immune system and empower diagnostics and therapeutics. There are currently multiple tools available for the annotation of antibody sequences. All downstream analyses such as choosing lead drug candidates depend on the correct annotation of these sequences; however, a thorough comparison of the performance of these tools has not been investigated. Here, we benchmark the performance of commonly used immunoinformatic tools, i.e. IMGT/HighV-QUEST, IgBLAST and MiXCR, in terms of reproducibility of annotation output, accuracy and speed using simulated and experimental high-throughput sequencing datasets.We analyzed changes in IMGT reference germline database in the last 10 years in order to assess the reproducibility of the annotation output. We found that only 73/183 (40%) V, D and J human genes were shared between the reference germline sets used by the tools. We found that the annotation results differed between tools. In terms of alignment accuracy, MiXCR had the highest average frequency of gene mishits, 0.02 mishit frequency and IgBLAST the lowest, 0.004 mishit frequency. Reproducibility in the output of complementarity determining three regions (CDR3 amino acids) ranged from 4.3% to 77.6% with preprocessed data. In addition, run time of the tools was assessed: MiXCR was the fastest tool for number of sequences processed per unit of time. These results indicate that immunoinformatic analyses greatly depend on the choice of bioinformatics tool. Our results support informed decision-making to immunoinformaticians based on repertoire composition and sequencing platforms. AVAILABILITY AND IMPLEMENTATION: All tools utilized in the paper are free for academic use. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Phenotypic and Ig Repertoire Analyses Indicate a Common Origin of IgD-CD27- Double Negative B Cells in Healthy Individuals and Multiple Sclerosis Patients.

IgD-CD27- double negative (DN) B cells with proinflammatory characteristics are abnormally elevated in a proportion of multiple sclerosis (MS) patients. In this study, the origin and selection characteristics of DN B cells were studied in MS patients and healthy controls (HC). Expression of developmental markers on peripheral blood DN, IgD-CD27+ class-switched memory (CSM) and IgD+CD27- naive B cells of HC (n = 48) and MS patients (n = 96) was determined by flow cytometry. High-throughput adaptive immune receptor repertoire sequencing was performed on peripheral blood DN and CSM B cells of HC and MS patients (n = 3 each). DN B cells from HC and MS patients showed similar phenotypic and Ig repertoire characteristics. Phenotypic analysis indicated a mature state of DN B cells by low CD5, CD10, and CD38 expression. However, the frequency of CD95+ and IgA+ cells was lower in DN versus CSM B cells. DN B cells are Ag experienced, as shown by somatic hypermutation of their Ig genes in adaptive immune receptor repertoire sequencing, although they showed a lower mutation load than CSM B cells. Shared clones were found between DN and CSM B cells, although >95% of the clones were unique to each population, and differences in V(D)J usage and CDR3 physicochemical properties were found. Thus, DN B cells arise in HC and MS patients via a common developmental pathway that is probably linked to immune aging. However, DN and CSM B cells develop through unique differentiation pathways, with most DN B cells representing an earlier maturation state.

The Med31 Conserved Component of the Divergent Mediator Complex in Tetrahymena thermophila Participates in Developmental Regulation.

Mediator is a large protein complex required for basal and regulated expression of most RNA polymerase II (RNAP II)-transcribed genes, in part due to its interaction with and phosphorylation of the conserved C-terminal domain (CTD) of Rpb1 [1, 2]. Mediator has been implicated in many aspects of gene expression including chromatin looping [3], higher-order chromatin folding [4], mRNA processing [5] and export [6], and transcriptional memory [7]. Mediator is thought to have played a major role during eukaryotic diversification [8, 9], although its function remains unknown in evolutionarily deep branching eukaryotes lacking canonical CTD heptad repeats. We used the ciliate protozoan Tetrahymena thermophila as a model organism whose genome encodes a highly divergent Rpb1 lacking canonical CTD heptad repeats. We endogenously tagged the Med31 subunit of the Mediator complex and performed affinity purification coupled with mass spectrometry (AP-MS) to identify Mediator subunits. We found that Med31 physically interacts with a large number of proteins (>20), several of which share similarities to canonical Mediator subunits in yeast and humans as well as Tetrahymena-specific proteins. Furthermore, Med31 ChIP-seq analysis suggested a global role for Mediator in transcription regulation. We demonstrated that MED31 knockdown in growing Tetrahymena results in the ectopic expression of developmental genes important for programmed DNA rearrangements. In addition, indirect immunofluorescence revealed Med31 localization in meiotic micronuclei, implicating Mediator in RNAPII-dependent ncRNA transcription. Our results reveal structural and functional insights and implicate Mediator as an ancient cellular machinery for transcription regulation with a possible involvement in global transcription of ncRNAs.

Identification of Subject-Specific Immunoglobulin Alleles From Expressed Repertoire Sequencing Data.

The adaptive immune receptor repertoire (AIRR) contains information on an individuals' immune past, present and potential in the form of the evolving sequences that encode the B cell receptor (BCR) repertoire. AIRR sequencing (AIRR-seq) studies rely on databases of known BCR germline variable (V), diversity (D), and joining (J) genes to detect somatic mutations in AIRR-seq data via comparison to the best-aligning database alleles. However, it has been shown that these databases are far from complete, leading to systematic misidentification of mutated positions in subsets of sample sequences. We previously presented TIgGER, a computational method to identify subject-specific V gene genotypes, including the presence of novel V gene alleles, directly from AIRR-seq data. However, the original algorithm was unable to detect alleles that differed by more than 5 single nucleotide polymorphisms (SNPs) from a database allele. Here we present and apply an improved version of the TIgGER algorithm which can detect alleles that differ by any number of SNPs from the nearest database allele, and can construct subject-specific genotypes with minimal prior information. TIgGER predictions are validated both computationally (using a leave-one-out strategy) and experimentally (using genomic sequencing), resulting in the addition of three new immunoglobulin heavy chain V (IGHV) gene alleles to the IMGT repertoire. Finally, we develop a Bayesian strategy to provide a confidence estimate associated with genotype calls. All together, these methods allow for much higher accuracy in germline allele assignment, an essential step in AIRR-seq studies.

AIRR Community Standardized Representations for Annotated Immune Repertoires.

Increased interest in the immune system's involvement in pathophysiological phenomena coupled with decreased DNA sequencing costs have led to an explosion of antibody and T cell receptor sequencing data collectively termed "adaptive immune receptor repertoire sequencing" (AIRR-seq or Rep-Seq). The AIRR Community has been actively working to standardize protocols, metadata, formats, APIs, and other guidelines to promote open and reproducible studies of the immune repertoire. In this paper, we describe the work of the AIRR Community's Data Representation Working Group to develop standardized data representations for storing and sharing annotated antibody and T cell receptor data. Our file format emphasizes ease-of-use, accessibility, scalability to large data sets, and a commitment to open and transparent science. It is composed of a tab-delimited format with a specific schema. Several popular repertoire analysis tools and data repositories already utilize this AIRR-seq data format. We hope that others will follow suit in the interest of promoting interoperable standards.

Generation of Long-Lived Bone Marrow Plasma Cells Secreting Antibodies Specific for the HIV-1 gp41 Membrane-Proximal External Region in the Absence of Polyreactivity.

UNLABELLED: An effective preventive vaccine is highly sought after in order to stem the current HIV-1 pandemic. Both conservation of contiguous gp41 membrane-proximal external region (MPER) amino acid sequences across HIV-1 clades and the ability of anti-MPER broadly neutralizing antibodies (BNAbs) to block viral hemifusion/fusion establish the MPER as a prime vaccination target. In earlier studies, we described the development of an MPER vaccine formulation that takes advantage of liposomes to array the MPER on a lipid bilayer surface, paralleling its native configuration on the virus membrane while also incorporating molecular adjuvant and CD4 T cell epitope cargo. Here we demonstrate that several immunizations with MPER/liposomes induce high levels of bone marrow long-lived plasma cell (LLPC) antibody production. Single-cell immunoglobulin gene retrieval analysis shows that these plasma cells are derived from a germ line repertoire of B cells with a diverse representation of immunoglobulin genes, exhibiting antigen-driven positive selection. Characterization of LLPC recombinant monoclonal antibodies (rMAbs) indicates that antigen recognition is achieved through convergence on a common epitopic focus by utilizing various complementarity-determining region H3 (CDRH3) lengths. Importantly, the vast majority of rMAbs produced from these cells lack polyreactivity yet manifest antigen specificity in the context of lipids, shaping MPER-specific paratopes through selective pressure. Taken together, these findings demonstrate that the MPER is a vaccine target with minimal risk of generating off-target autoimmunity. IMPORTANCE: A useful vaccine must generate desired long-term, antigen-specific antibody responses devoid of polyreactivity or autoreactivity. The common polyreactive features of some HIV-1 BNAbs have raised concern about elicitation of anti-MPER antibodies. Utilizing single-LLPC repertoire analysis and biophysical characterization of anti-MPER rMAbs, we show that their fine specificities require a structural fitness of the antibody combining site involving heavy and light chain variable domains shaped by somatic hypermutation and affinity maturation of B cells in the germinal center. Perhaps more importantly, our results demonstrate that the majority of MPER-specific antibodies are not inherently polyspecific and/or autoreactive, suggesting that polyreactivity of MPER-specific antibodies is separable from their antigen specificity.

Mapping of functionally characterized cell classes onto canonical circuit operations in primate prefrontal cortex.

Microcircuits are composed of multiple cell classes that likely serve unique circuit operations. But how cell classes map onto circuit functions is largely unknown, particularly for primate prefrontal cortex during actual goal-directed behavior. One difficulty in this quest is to reliably distinguish cell classes in extracellular recordings of action potentials. Here we surmount this issue and report that spike shape and neural firing variability provide reliable markers to segregate seven functional classes of prefrontal cells in macaques engaged in an attention task. We delineate an unbiased clustering protocol that identifies four broad spiking (BS) putative pyramidal cell classes and three narrow spiking (NS) putative inhibitory cell classes dissociated by how sparse, bursty, or regular they fire. We speculate that these functional classes map onto canonical circuit functions. First, two BS classes show sparse, bursty firing, and phase synchronize their spiking to 3-7 Hz (theta) and 12-20 Hz (beta) frequency bands of the local field potential (LFP). These properties make cells flexibly responsive to network activation at varying frequencies. Second, one NS and two BS cell classes show regular firing and higher rate with only marginal synchronization preference. These properties are akin to setting tonically the excitation and inhibition balance. Finally, two NS classes fired irregularly and synchronized to either theta or beta LFP fluctuations, tuning them potentially to frequency-specific subnetworks. These results suggest that a limited set of functional cell classes emerges in macaque prefrontal cortex (PFC) during attentional engagement to not only represent information, but to subserve basic circuit operations.

Conserved Asf1-importin ß physical interaction in growth and sexual development in the ciliate Tetrahymena thermophila.

How the eukaryotic cell specifies distinct chromatin domains is a central problem in molecular biology. The ciliate protozoan Tetrahymena thermophila features a separation of structurally and functionally distinct germ-line and somatic chromatin into two distinct nuclei, the micronucleus (MIC) and macronucleus (MAC) respectively. To address questions about how distinct chromatin states are assembled in the MAC and MIC, we have initiated studies to define protein-protein interactions for T. thermophila chromatin-related proteins. Affinity purification followed by mass spectrometry analysis of the conserved Asf1 histone chaperone in T. thermophila revealed that it forms a complex with an importin ß, ImpB6. Furthermore, these proteins co-localized to both the MAC and MIC in growth and development. We suggest that newly synthesized histones H3 and H4 in T. thermophila are transported via Asf1-ImpB6 in an evolutionarily conserved pathway to both nuclei where they then enter nucleus-specific chromatin assembly pathways. These studies set the stage for further use of functional proteomics to elucidate details of the characterization and functional analysis of the unique chromatin domains in T. thermophila. BIOLOGICAL SIGNIFICANCE: Asf1 is an evolutionarily conserved chaperone of H3 and H4 histones that functions in replication dependent and independent chromatin assembly. Although Asf1 has been well studied in humans and yeast (members of the Opisthokonta lineage of eukaryotes), questions remain concerning its mechanism of function. To obtain additional insight into the Asf1 function we have initiated a proteomic analysis in the ciliate protozoan T. thermophila, a member of the Alveolata lineage of eukaryotes. Our results suggest that an evolutionarily conserved function of Asf1 is mediating the nuclear transport of newly synthesized histones H3 and H4.

Antigen-specific, antibody-coated, exosome-like nanovesicles deliver suppressor T-cell microRNA-150 to effector T cells to inhibit contact sensitivity.

BACKGROUND: T-cell tolerance of allergic cutaneous contact sensitivity (CS) induced in mice by high doses of reactive hapten is mediated by suppressor cells that release antigen-specific suppressive nanovesicles. OBJECTIVE: We sought to determine the mechanism or mechanisms of immune suppression mediated by the nanovesicles. METHODS: T-cell tolerance was induced by means of intravenous injection of hapten conjugated to self-antigens of syngeneic erythrocytes and subsequent contact immunization with the same hapten. Lymph node and spleen cells from tolerized or control donors were harvested and cultured to produce a supernatant containing suppressive nanovesicles that were isolated from the tolerized mice for testing in active and adoptive cell-transfer models of CS. RESULTS: Tolerance was shown due to exosome-like nanovesicles in the supernatants of CD8(+) suppressor T cells that were not regulatory T cells. Antigen specificity of the suppressive nanovesicles was conferred by a surface coat of antibody light chains or possibly whole antibody, allowing targeted delivery of selected inhibitory microRNA (miRNA)-150 to CS effector T cells. Nanovesicles also inhibited CS in actively sensitized mice after systemic injection at the peak of the responses. The role of antibody and miRNA-150 was established by tolerizing either panimmunoglobulin-deficient JH(-/-) or miRNA-150(-/-) mice that produced nonsuppressive nanovesicles. These nanovesicles could be made suppressive by adding antigen-specific antibody light chains or miRNA-150, respectively. CONCLUSIONS: This is the first example of T-cell regulation through systemic transit of exosome-like nanovesicles delivering a chosen inhibitory miRNA to target effector T cells in an antigen-specific manner by a surface coating of antibody light chains.

Reconstruction of regulatory networks through temporal enrichment profiling and its application to H1N1 influenza viral infection.

BACKGROUND: H1N1 influenza viruses were responsible for the 1918 pandemic that caused millions of deaths worldwide and the 2009 pandemic that caused approximately twenty thousand deaths. The cellular response to such virus infections involves extensive genetic reprogramming resulting in an antiviral state that is critical to infection control. Identifying the underlying transcriptional network driving these changes, and how this program is altered by virally-encoded immune antagonists, is a fundamental challenge in systems immunology. RESULTS: Genome-wide gene expression patterns were measured in human monocyte-derived dendritic cells (DCs) infected in vitro with seasonal H1N1 influenza A/New Caledonia/20/1999. To provide a mechanistic explanation for the timing of gene expression changes over the first 12 hours post-infection, we developed a statistically rigorous enrichment approach integrating genome-wide expression kinetics and time-dependent promoter analysis. Our approach, TIme-Dependent Activity Linker (TIDAL), generates a regulatory network that connects transcription factors associated with each temporal phase of the response into a coherent linked cascade. TIDAL infers 12 transcription factors and 32 regulatory connections that drive the antiviral response to influenza. To demonstrate the generality of this approach, TIDAL was also used to generate a network for the DC response to measles infection. The software implementation of TIDAL is freely available at http://tsb.mssm.edu/primeportal/?q=tidal_prog. CONCLUSIONS: We apply TIDAL to reconstruct the transcriptional programs activated in monocyte-derived human dendritic cells in response to influenza and measles infections. The application of this time-centric network reconstruction method in each case produces a single transcriptional cascade that recapitulates the known biology of the response with high precision and recall, in addition to identifying potentially novel antiviral factors. The ability to reconstruct antiviral networks with TIDAL enables comparative analysis of antiviral responses, such as the differences between pandemic and seasonal influenza infections.

Antiviral response dictated by choreographed cascade of transcription factors.

The dendritic cell (DC) is a master regulator of immune responses. Pathogenic viruses subvert normal immune function in DCs through the expression of immune antagonists. Understanding how these antagonists interact with the host immune system requires knowledge of the underlying genetic regulatory network that operates during an uninhibited antiviral response. To isolate and identify this network, we studied DCs infected with Newcastle disease virus, which is able to stimulate innate immunity and DC maturation through activation of RIG-I signaling, but lacks the ability to evade the human IFN response. To analyze this experimental model, we developed a new approach integrating genome-wide expression kinetics and time-dependent promoter analysis. We found that the genetic program underlying the antiviral cell-state transition during the first 18 h postinfection could be explained by a single convergent regulatory network. Gene expression changes were driven by a stepwise multifactor cascading control mechanism, where the specific transcription factors controlling expression changed over time. Within this network, most individual genes were regulated by multiple factors, indicating robustness against virus-encoded immune evasion genes. In addition to effectively recapitulating current biological knowledge, we predicted, and validated experimentally, antiviral roles for several novel transcription factors. More generally, our results show how a genetic program can be temporally controlled through a single regulatory network to achieve the large-scale genetic reprogramming characteristic of cell-state transitions.