ABSTRACT
Acute Respiratory Distress Syndrome (ARDS) is an illness that typically develops in people who are significantly ill or have serious injuries. ARDS is characterized by fluid build-up that occurs in the alveoli. T-cells are implicated as playing a role in the modulation of the aberrant response leading to excessive tissue damage and, eventually, ARDS. Complementarity Determining Region 3 (CDR3) sequences derived from T-cells are key players in the adaptive immune response. This response is governed by an elaborate specificity for distinct molecules and the ability to recognize and vigorously respond to repeated exposures to the same molecules. Most of the diversity in T-cell receptors (TCRs) is contained in the CDR3 regions of the heterodimeric cell-surface receptors. For this study, we employed the novel technology of immune sequencing to assess lung edema fluid. Our goal was to explore the landscape of CDR3 clonal sequences found within these samples. We obtained more than 3615 CDR3 sequences across samples in the study. Our data demonstrate that: (1) CDR3 sequences from lung edema fluid exhibit distinct clonal populations, and (2) CDR3 sequences can be further characterized based on biochemical features. Analysis of these CDR3 sequences offers insight into the CDR3-driven T-cell repertoire of ARDS. These findings represent the first step towards applications of this technology with these types of biological samples in the context of ARDS.
Subject(s)
Complementarity Determining Regions , Respiratory Distress Syndrome , Humans , Complementarity Determining Regions/genetics , Receptors, Antigen, T-Cell, alpha-beta/genetics , Receptors, Antigen, T-Cell , Respiratory Distress Syndrome/genetics , EdemaABSTRACT
Immune repertoire deep sequencing allows comprehensive characterization of antigen receptor-encoding genes in a lymphocyte population. We hypothesized that this method could enable a novel approach to diagnose disease by identifying antigen receptor sequence patterns associated with clinical phenotypes. In this study, we developed statistical classifiers of T-cell receptor (TCR) repertoires that distinguish tumor tissue from patient-matched healthy tissue of the same organ. The basis of both classifiers was a biophysicochemical motif in the complementarity determining region 3 (CDR3) of TCRß chains. To develop each classifier, we extracted 4-mers from every TCRß CDR3 and represented each 4-mer using biophysicochemical features of its amino acid sequence combined with quantification of 4-mer (or receptor) abundance. This representation was scored using a logistic regression model. Unlike typical logistic regression, the classifier is fitted and validated under the requirement that at least 1 positively labeled 4-mer appears in every tumor repertoire and no positively labeled 4-mers appear in healthy tissue repertoires. We applied our method to publicly available data in which tumor and adjacent healthy tissue were collected from each patient. Using a patient-holdout cross-validation, our method achieved classification accuracy of 93% and 94% for colorectal and breast cancer, respectively. The parameter values for each classifier revealed distinct biophysicochemical properties for tumor-associated 4-mers within each cancer type. We propose that such motifs might be used to develop novel immune-based cancer screening assays. SIGNIFICANCE: This study presents a novel computational approach to identify T-cell repertoire differences between normal and tumor tissue.See related commentary by Zoete and Coukos, p. 1299.
Subject(s)
Lymphocytes, Tumor-Infiltrating , Receptors, Antigen, T-Cell, alpha-beta/chemistry , Complementarity Determining Regions/chemistry , Humans , Receptors, Antigen, T-Cell/genetics , T-Lymphocytes/immunologyABSTRACT
Background: Recent technological advances in immune repertoire sequencing have created tremendous potential for advancing our understanding of adaptive immune response dynamics in various states of health and disease. Immune repertoire sequencing produces large, highly complex data sets, however, which require specialized methods and software tools for their effective analysis and interpretation. Results: VDJServer is a cloud-based analysis portal for immune repertoire sequence data that provide access to a suite of tools for a complete analysis workflow, including modules for preprocessing and quality control of sequence reads, V(D)J gene segment assignment, repertoire characterization, and repertoire comparison. VDJServer also provides sophisticated visualizations for exploratory analysis. It is accessible through a standard web browser via a graphical user interface designed for use by immunologists, clinicians, and bioinformatics researchers. VDJServer provides a data commons for public sharing of repertoire sequencing data, as well as private sharing of data between users. We describe the main functionality and architecture of VDJServer and demonstrate its capabilities with use cases from cancer immunology and autoimmunity. Conclusion: VDJServer provides a complete analysis suite for human and mouse T-cell and B-cell receptor repertoire sequencing data. The combination of its user-friendly interface and high-performance computing allows large immune repertoire sequencing projects to be analyzed with no programming or software installation required. VDJServer is a web-accessible cloud platform that provides access through a graphical user interface to a data management infrastructure, a collection of analysis tools covering all steps in an analysis, and an infrastructure for sharing data along with workflows, results, and computational provenance. VDJServer is a free, publicly available, and open-source licensed resource.
Subject(s)
Cloud Computing , Computational Biology/methods , Genomics/methods , VDJ Exons/immunology , Animals , Computing Methodologies , Humans , Information Dissemination , Mice , Software , User-Computer Interface , Web Browser , WorkflowABSTRACT
BACKGROUND: Pre-processing of high-throughput sequencing data for immune repertoire profiling is essential to insure high quality input for downstream analysis. VDJPipe is a flexible, high-performance tool that can perform multiple pre-processing tasks with just a single pass over the data files. RESULTS: Processing tasks provided by VDJPipe include base composition statistics calculation, read quality statistics calculation, quality filtering, homopolymer filtering, length and nucleotide filtering, paired-read merging, barcode demultiplexing, 5' and 3' PCR primer matching, and duplicate reads collapsing. VDJPipe utilizes a pipeline approach whereby multiple processing steps are performed in a sequential workflow, with the output of each step passed as input to the next step automatically. The workflow is flexible enough to handle the complex barcoding schemes used in many immunosequencing experiments. Because VDJPipe is designed for computational efficiency, we evaluated this by comparing execution times with those of pRESTO, a widely-used pre-processing tool for immune repertoire sequencing data. We found that VDJPipe requires <10% of the run time required by pRESTO. CONCLUSIONS: VDJPipe is a high-performance tool that is optimized for pre-processing large immune repertoire sequencing data sets.
Subject(s)
B-Lymphocytes/metabolism , High-Throughput Nucleotide Sequencing/methods , Immunoglobulin G/genetics , Software , Animals , DNA Primers , Humans , Mice , Time FactorsABSTRACT
BACKGROUND: The genes that produce antibodies and the immune receptors expressed on lymphocytes are not germline encoded; rather, they are somatically generated in each developing lymphocyte by a process called V(D)J recombination, which assembles specific, independent gene segments into mature composite genes. The full set of composite genes in an individual at a single point in time is referred to as the immune repertoire. V(D)J recombination is the distinguishing feature of adaptive immunity and enables effective immune responses against an essentially infinite array of antigens. Characterization of immune repertoires is critical in both basic research and clinical contexts. Recent technological advances in repertoire profiling via high-throughput sequencing have resulted in an explosion of research activity in the field. This has been accompanied by a proliferation of software tools for analysis of repertoire sequencing data. Despite the widespread use of immune repertoire profiling and analysis software, there is currently no standardized format for output files from V(D)J analysis. Researchers utilize software such as IgBLAST and IMGT/High V-QUEST to perform V(D)J analysis and infer the structure of germline rearrangements. However, each of these software tools produces results in a different file format, and can annotate the same result using different labels. These differences make it challenging for users to perform additional downstream analyses. RESULTS: To help address this problem, we propose a standardized file format for representing V(D)J analysis results. The proposed format, VDJML, provides a common standardized format for different V(D)J analysis applications to facilitate downstream processing of the results in an application-agnostic manner. The VDJML file format specification is accompanied by a support library, written in C++ and Python, for reading and writing the VDJML file format. CONCLUSIONS: The VDJML suite will allow users to streamline their V(D)J analysis and facilitate the sharing of scientific knowledge within the community. The VDJML suite and documentation are available from https://vdjserver.org/vdjml/ . We welcome participation from the community in developing the file format standard, as well as code contributions.
Subject(s)
Genomics/methods , Receptors, Immunologic/genetics , Software , V(D)J Recombination , Humans , Information DisseminationABSTRACT
BACKGROUND: The Bacillus cereus sensu lato group contains ubiquitous facultative anaerobic soil-borne Gram-positive spore-forming bacilli. Molecular phylogeny and comparative genome sequencing have suggested that these organisms should be classified as a single species. While clonal in nature, there do not appear to be species-specific clonal lineages, excepting B. anthracis, in spite of the wide array of phenotypes displayed by these organisms. RESULTS: We compared the protein-coding content of 201 B. cereus sensu lato genomes to characterize differences and understand the consequences of these differences on biological function. From this larger group we selected a subset consisting of 25 whole genomes for deeper analysis. Cluster analysis of orthologous proteins grouped these genomes into five distinct clades. Each clade could be characterized by unique genes shared among the group, with consequences for the phenotype of each clade. Surprisingly, this population structure recapitulates our recent observations on the divergence of the generalized stress response (SigB) regulons in these organisms. Divergence of the SigB regulon among these organisms is primarily due to the placement of SigB-dependent promoters that bring genes from a common gene pool into/out of the SigB regulon. CONCLUSIONS: Collectively, our observations suggest the hypothesis that the evolution of these closely related bacteria is a consequence of two distinct processes. Horizontal gene transfer, gene duplication/divergence and deletion dictate the underlying coding capacity in these genomes. Regulatory divergence overlays this protein coding reservoir and shapes the expression of both the unique and shared coding capacity of these organisms, resulting in phenotypic divergence. Data from other organisms suggests that this is likely a common pattern in prokaryotic evolution.
Subject(s)
Bacillus cereus/genetics , Bacterial Proteins/genetics , Bacillus cereus/classification , Bacillus cereus/metabolism , Cluster Analysis , Evolution, Molecular , Genome, Bacterial , Phenotype , Phylogeny , RegulonABSTRACT
SCB34 is a sequence type 131, highly invasive, multidrug-resistant Escherichia coli isolate that produced neonatal bacteremia. Whole-genome sequencing was performed using a 250-bp library on the Illumina MiSeq platform; 5,910,264 reads were assembled de novo using the A5 assembly pipeline. The total contig length was 5,227,742 bp; the RAST server was used for annotation.
ABSTRACT
We hypothesized that hypoxia would activate epidermal growth factor receptor (EGFR) tyrosine kinase, leading to increased arginase expression and resulting in proliferation of human pulmonary microvascular endothelial cell (hPMVEC). To test this hypothesis, hPMVEC were incubated in normoxia (20% O(2), 5% CO(2)) or hypoxia (1% O(2), 5% CO(2)). Immunoblotting for EGFR and proliferating cell nuclear antigen was done, and protein levels of both total EGFR and proliferating cell nuclear antigen were greater in hypoxic hPMVEC than in normoxic hPMVEC. Furthermore, hypoxic hPMVEC had greater levels of EGFR activity than did normoxic hPMVEC. Hypoxic hPMVEC had a twofold greater level of proliferation compared with normoxic controls, and this increase in proliferation was prevented by the addition of AG-1478 (a pharmacological inhibitor of EGFR). Immunoblotting for arginase I and arginase II demonstrated a threefold induction in arginase II protein levels in hypoxia, with little change in arginase I protein levels. The hypoxic induction of arginase II protein was prevented by treatment with AG-1478. Proliferation assays were performed in the presence of arginase inhibitors, and hypoxia-induced proliferation was also prevented by arginase inhibition. Finally, treatment with an EGFR small interfering RNA prevented hypoxia-induced proliferation and urea production. These findings demonstrate that hypoxia activates EGFR tyrosine kinase, leading to arginase expression and thereby promoting proliferation in hPMVEC.