IMGT/HighV QUEST paradigm for T cell receptor IMGT clonotype diversity and next generation repertoire immunoprofiling.

T cell repertoire diversity and clonotype follow-up in vaccination, cancer, infectious and immune diseases represent a major challenge owing to the enormous complexity of the data generated. Here we describe a next generation methodology, which combines 5'RACE PCR, 454 sequencing and, for analysis, IMGT, the international ImMunoGeneTics information system (IMGT), IMGT/HighV-QUEST web portal and IMGT-ONTOLOGY concepts. The approach is validated in a human case study of T cell receptor beta (TRB) repertoire, by chronologically tracking the effects of influenza vaccination on conventional and regulatory T cell subpopulations. The IMGT/HighV-QUEST paradigm defines standards for genotype/haplotype analysis and characterization of IMGT clonotypes for clonal diversity and expression and achieves a degree of resolution for next generation sequencing verifiable by the user at the sequence level, while providing a normalized reference immunoprofile for human TRB.

Fusobacterium nucleatum infection is prevalent in human colorectal carcinoma.

An estimated 15% or more of the cancer burden worldwide is attributable to known infectious agents. We screened colorectal carcinoma and matched normal tissue specimens using RNA-seq followed by host sequence subtraction and found marked over-representation of Fusobacterium nucleatum sequences in tumors relative to control specimens. F. nucleatum is an invasive anaerobe that has been linked previously to periodontitis and appendicitis, but not to cancer. Fusobacteria are rare constituents of the fecal microbiota, but have been cultured previously from biopsies of inflamed gut mucosa. We obtained a Fusobacterium isolate from a frozen tumor specimen; this showed highest sequence similarity to a known gut mucosa isolate and was confirmed to be invasive. We verified overabundance of Fusobacterium sequences in tumor versus matched normal control tissue by quantitative PCR analysis from a total of 99 subjects (p = 2.5 × 10(-6)), and we observed a positive association with lymph node metastasis.

The sensitivity of massively parallel sequencing for detecting candidate infectious agents associated with human tissue.

Massively parallel sequencing technology now provides the opportunity to sample the transcriptome of a given tissue comprehensively. Transcripts at only a few copies per cell are readily detectable, allowing the discovery of low abundance viral and bacterial transcripts in human tissue samples. Here we describe an approach for mining large sequence data sets for the presence of microbial sequences. Further, we demonstrate the sensitivity of this approach by sequencing human RNA-seq libraries spiked with decreasing amounts of an RNA-virus. At a modest depth of sequencing, viral transcripts can be detected at frequencies less than 1 in 1,000,000. With current sequencing platforms approaching outputs of one billion reads per run, this is a highly sensitive method for detecting putative infectious agents associated with human tissues.

Exhaustive T-cell repertoire sequencing of human peripheral blood samples reveals signatures of antigen selection and a directly measured repertoire size of at least 1 million clonotypes.

Massively parallel sequencing is a useful approach for characterizing T-cell receptor diversity. However, immune receptors are extraordinarily difficult sequencing targets because any given receptor variant may be present in very low abundance and may differ legitimately by only a single nucleotide. We show that the sensitivity of sequence-based repertoire profiling is limited by both sequencing depth and sequencing accuracy. At two timepoints, 1 wk apart, we isolated bulk PBMC plus naïve (CD45RA+/CD45RO-) and memory (CD45RA-/CD45RO+) T-cell subsets from a healthy donor. From T-cell receptor beta chain (TCRB) mRNA we constructed and sequenced multiple libraries to obtain a total of 1.7 billion paired sequence reads. The sequencing error rate was determined empirically and used to inform a high stringency data filtering procedure. The error filtered data yielded 1,061,522 distinct TCRB nucleotide sequences from this subject which establishes a new, directly measured, lower limit on individual T-cell repertoire size and provides a useful reference set of sequences for repertoire analysis. TCRB nucleotide sequences obtained from two additional donors were compared to those from the first donor and revealed limited sharing (up to 1.1%) of nucleotide sequences among donors, but substantially higher sharing (up to 14.2%) of inferred amino acid sequences. For each donor, shared amino acid sequences were encoded by a much larger diversity of nucleotide sequences than were unshared amino acid sequences. We also observed a highly statistically significant association between numbers of shared sequences and shared HLA class I alleles.

Profiling the T-cell receptor beta-chain repertoire by massively parallel sequencing.

T-cell receptor (TCR) genomic loci undergo somatic V(D)J recombination, plus the addition/subtraction of nontemplated bases at recombination junctions, in order to generate the repertoire of structurally diverse T cells necessary for antigen recognition. TCR beta subunits can be unambiguously identified by their hypervariable CDR3 (Complement Determining Region 3) sequence. This is the site of V(D)J recombination encoding the principal site of antigen contact. The complexity and dynamics of the T-cell repertoire remain unknown because the potential repertoire size has made conventional sequence analysis intractable. Here, we use 5'-RACE, Illumina sequencing, and a novel short read assembly strategy to sample CDR3(beta) diversity in human T lymphocytes from peripheral blood. Assembly of 40.5 million short reads identified 33,664 distinct TCR(beta) clonotypes and provides precise measurements of CDR3(beta) length diversity, usage of nontemplated bases, sequence convergence, and preferences for TRBV (T-cell receptor beta variable gene) and TRBJ (T-cell receptor beta joining gene) gene usage and pairing. CDR3 length between conserved residues of TRBV and TRBJ ranged from 21 to 81 nucleotides (nt). TRBV gene usage ranged from 0.01% for TRBV17 to 24.6% for TRBV20-1. TRBJ gene usage ranged from 1.6% for TRBJ2-6 to 17.2% for TRBJ2-1. We identified 1573 examples of convergence where the same amino acid translation was specified by distinct CDR3(beta) nucleotide sequences. Direct sequence-based immunoprofiling will likely prove to be a useful tool for understanding repertoire dynamics in response to immune challenge, without a priori knowledge of antigen.

Steroid hormone control of cell death and cell survival: molecular insights using RNAi.

The insect steroid hormone ecdysone triggers programmed cell death of obsolete larval tissues during metamorphosis and provides a model system for understanding steroid hormone control of cell death and cell survival. Previous genome-wide expression studies of Drosophila larval salivary glands resulted in the identification of many genes associated with ecdysone-induced cell death and cell survival, but functional verification was lacking. In this study, we test functionally 460 of these genes using RNA interference in ecdysone-treated Drosophila l(2)mbn cells. Cell viability, cell morphology, cell proliferation, and apoptosis assays confirmed the effects of known genes and additionally resulted in the identification of six new pro-death related genes, including sorting nexin-like gene SH3PX1 and Sox box protein Sox14, and 18 new pro-survival genes. Identified genes were further characterized to determine their ecdysone dependency and potential function in cell death regulation. We found that the pro-survival function of five genes (Ras85D, Cp1, CG13784, CG32016, and CG33087), was dependent on ecdysone signaling. The TUNEL assay revealed an additional two genes (Kap-alpha3 and Smr) with an ecdysone-dependent cell survival function that was associated with reduced cell death. In vitro, Sox14 RNAi reduced the percentage of TUNEL-positive l(2)mbn cells (p<0.05) following ecdysone treatment, and Sox14 overexpression was sufficient to induce apoptosis. In vivo analyses of Sox14-RNAi animals revealed multiple phenotypes characteristic of aberrant or reduced ecdysone signaling, including defects in larval midgut and salivary gland destruction. These studies identify Sox14 as a positive regulator of ecdysone-mediated cell death and provide new insights into the molecular mechanisms underlying the ecdysone signaling network governing cell death and cell survival.

echinus, required for interommatidial cell sorting and cell death in the Drosophila pupal retina, encodes a protein with homology to ubiquitin-specific proteases.

BACKGROUND: Programmed cell death is used to remove excess cells between ommatidia in the Drosophila pupal retina. This death is required to establish the crystalline, hexagonal packing of ommatidia that characterizes the adult fly eye. In previously described echinus mutants, interommatidial cell sorting, which precedes cell death, occurred relatively normally. Interommatidial cell death was partially suppressed, resulting in adult eyes that contained excess pigment cells, and in which ommatidia were mildly disordered. These results have suggested that echinus functions in the pupal retina primarily to promote interommatidial cell death. RESULTS: We generated a number of new echinus alleles, some likely null mutants. Analysis of these alleles provides evidence that echinus has roles in cell sorting as well as cell death. echinus encodes a protein with homology to ubiquitin-specific proteases. These proteins cleave ubiquitin-conjugated proteins at the ubiquitin C-terminus. The echinus locus encodes multiple splice forms, including two proteins that lack residues thought to be critical for deubiquitination activity. Surprisingly, ubiquitous expression in the eye of versions of Echinus that lack residues critical for ubiquitin specific protease activity, as well as a version predicted to be functional, rescue the echinus loss-of-function phenotype. Finally, genetic interactions were not detected between echinus loss and gain-of-function and a number of known apoptotic regulators. These include Notch, EGFR, the caspases Dronc, Drice, Dcp-1, Dream, the caspase activators, Rpr, Hid, and Grim, the caspase inhibitor DIAP1, and Lozenge or Klumpfuss. CONCLUSION: The echinus locus encodes multiple splice forms of a protein with homology to ubiquitin-specific proteases, but protease activity is unlikely to be required for echinus function, at least when echinus is overexpressed. Characterization of likely echinus null alleles and genetic interactions suggests that echinus acts at a novel point(s) to regulate interommatidial cell sorting and/or cell death in the fly eye.