Hepatitis C virus RNA is 5'-capped with flavin adenine dinucleotide.

RNA viruses have evolved elaborate strategies to protect their genomes, including 5' capping. However, until now no RNA 5' cap has been identified for hepatitis C virus1,2 (HCV), which causes chronic infection, liver cirrhosis and cancer3. Here we demonstrate that the cellular metabolite flavin adenine dinucleotide (FAD) is used as a non-canonical initiating nucleotide by the viral RNA-dependent RNA polymerase, resulting in a 5'-FAD cap on the HCV RNA. The HCV FAD-capping frequency is around 75%, which is the highest observed for any RNA metabolite cap across all kingdoms of life4-8. FAD capping is conserved among HCV isolates for the replication-intermediate negative strand and partially for the positive strand. It is also observed in vivo on HCV RNA isolated from patient samples and from the liver and serum of a human liver chimeric mouse model. Furthermore, we show that 5'-FAD capping protects RNA from RIG-I mediated innate immune recognition but does not stabilize the HCV RNA. These results establish capping with cellular metabolites as a novel viral RNA-capping strategy, which could be used by other viruses and affect anti-viral treatment outcomes and persistence of infection.

Riboswitch-Mediated Gene Regulation: Novel RNA Architectures Dictate Gene Expression Responses.

Riboswitches are RNA elements that act on the mRNA with which they are cotranscribed to modulate expression of that mRNA. These elements are widely found in bacteria, where they have a broad impact on gene expression. The defining feature of riboswitches is that they directly recognize a physiological signal, and the resulting shift in RNA structure affects gene regulation. The majority of riboswitches respond to cellular metabolites, often in a feedback loop to repress synthesis of the enzymes used to produce the metabolite. Related elements respond to the aminoacylation status of a specific tRNA or to a physical parameter, such as temperature or pH. Recent studies have identified new classes of riboswitches and have revealed new insights into the molecular mechanisms of signal recognition and gene regulation. Application of structural and biophysical approaches has complemented previous genetic and biochemical studies, yielding new information about how different riboswitches operate.

New tRNA contacts facilitate ligand binding in a Mycobacterium smegmatis T box riboswitch.

T box riboswitches are RNA regulatory elements widely used by organisms in the phyla Firmicutes and Actinobacteria to regulate expression of amino acid-related genes. Expression of T box family genes is down-regulated by transcription attenuation or inhibition of translation initiation in response to increased charging of the cognate tRNA. Three direct contacts with tRNA have been described; however, one of these contacts is absent in a subclass of T box RNAs and the roles of several structural domains conserved in most T box RNAs are unknown. In this study, structural elements of a Mycobacterium smegmatis ileS T box riboswitch variant with an Ultrashort (US) Stem I were sequentially deleted, which resulted in a progressive decrease in binding affinity for the tRNAIle ligand. Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) revealed structural changes in conserved riboswitch domains upon interaction with the tRNA ligand. Cross-linking and mutational analyses identified two interaction sites, one between the S-turn element in Stem II and the T arm of tRNAIle and the other between the Stem IIA/B pseudoknot and the D loop of tRNAIle These newly identified RNA contacts add information about tRNA recognition by the T box riboswitch and demonstrate a role for the S-turn and pseudoknot elements, which resemble structural elements that are common in many cellular RNAs.

Analytical evaluation of the clonoSEQ Assay for establishing measurable (minimal) residual disease in acute lymphoblastic leukemia, chronic lymphocytic leukemia, and multiple myeloma.

BACKGROUND: The clonoSEQ® Assay (Adaptive Biotechnologies Corporation, Seattle, USA) identifies and tracks unique disease-associated immunoglobulin (Ig) sequences by next-generation sequencing of IgH, IgK, and IgL rearrangements and IgH-BCL1/2 translocations in malignant B cells. Here, we describe studies to validate the analytical performance of the assay using patient samples and cell lines. METHODS: Sensitivity and specificity were established by defining the limit of detection (LoD), limit of quantitation (LoQ) and limit of blank (LoB) in genomic DNA (gDNA) from 66 patients with multiple myeloma (MM), acute lymphoblastic leukemia (ALL), or chronic lymphocytic leukemia (CLL), and three cell lines. Healthy donor gDNA was used as a diluent to contrive samples with specific DNA masses and malignant-cell frequencies. Precision was validated using a range of samples contrived from patient gDNA, healthy donor gDNA, and 9 cell lines to generate measurable residual disease (MRD) frequencies spanning clinically relevant thresholds. Linearity was determined using samples contrived from cell line gDNA spiked into healthy gDNA to generate 11 MRD frequencies for each DNA input, then confirmed using clinical samples. Quantitation accuracy was assessed by (1) comparing clonoSEQ and multiparametric flow cytometry (mpFC) measurements of ALL and MM cell lines diluted in healthy mononuclear cells, and (2) analyzing precision study data for bias between clonoSEQ MRD results in diluted gDNA and those expected from mpFC based on original, undiluted samples. Repeatability of nucleotide base calls was assessed via the assay's ability to recover malignant clonotype sequences across several replicates, process features, and MRD levels. RESULTS: LoD and LoQ were estimated at 1.903 cells and 2.390 malignant cells, respectively. LoB was zero in healthy donor gDNA. Precision ranged from 18% CV (coefficient of variation) at higher DNA inputs to 68% CV near the LoD. Variance component analysis showed MRD results were robust, with expected laboratory process variations contributing ≤3% CV. Linearity and accuracy were demonstrated for each disease across orders of magnitude of clonal frequencies. Nucleotide sequence error rates were extremely low. CONCLUSIONS: These studies validate the analytical performance of the clonoSEQ Assay and demonstrate its potential as a highly sensitive diagnostic tool for selected lymphoid malignancies.

A Diverse Lipid Antigen-Specific TCR Repertoire Is Clonally Expanded during Active Tuberculosis.

Human T cells that recognize lipid Ags presented by highly conserved CD1 proteins often express semi-invariant TCRs, but the true diversity of lipid Ag-specific TCRs remains unknown. We use CD1b tetramers and high-throughput immunosequencing to analyze thousands of TCRs from ex vivo-sorted or in vitro-expanded T cells specific for the mycobacterial lipid Ag, glucose monomycolate. Our results reveal a surprisingly diverse repertoire resulting from editing of germline-encoded gene rearrangements analogous to MHC-restricted TCRs. We used a distance-based metric (TCRDist) to show how this diverse TCR repertoire builds upon previously reported conserved motifs by including subject-specific TCRs. In a South African cohort, we show that TCRDist can identify clonal expansion of diverse glucose monomycolate-specific TCRs and accurately distinguish patients with active tuberculosis from control subjects. These data suggest that similar mechanisms govern the selection and expansion of peptide and lipid Ag-specific T cells despite the nonpolymorphic nature of CD1.

Longitudinal immunosequencing in healthy people reveals persistent T cell receptors rich in highly public receptors.

BACKGROUND: The adaptive immune system maintains a diversity of T cells capable of recognizing a broad array of antigens. Each T cell's specificity for antigens is determined by its T cell receptors (TCRs), which together across all T cells form a repertoire of millions of unique receptors in each individual. Although many studies have examined how TCR repertoires change in response to disease or drugs, few have explored the temporal dynamics of the TCR repertoire in healthy individuals. RESULTS: Here we report immunosequencing of TCR ß chains (TCRß) from the blood of three healthy individuals at eight time points over one year. TCRß repertoires of all peripheral-blood T cells and sorted memory T cells clustered clearly by individual, systematically demonstrating that TCRß repertoires are specific to individuals across time. This individuality was absent from TCRßs from naive T cells, suggesting that the differences resulted from an individual's antigen exposure history, not genetic background. Many characteristics of the TCRß repertoire (e.g., diversity, clonality) were stable across time, although we found evidence of T cell expansion dynamics even within healthy individuals. We further identified a subset of "persistent" TCRßs present across all time points. These receptors were rich in clonal and highly public receptors and may play a key role in immune system maintenance. CONCLUSIONS: Our results highlight the importance of longitudinal sampling of the immune system, providing a much-needed baseline for TCRß dynamics in healthy individuals. Such a baseline will improve interpretation of changes in the TCRß repertoire during disease or treatment.

T box riboswitches in Actinobacteria: translational regulation via novel tRNA interactions.

The T box riboswitch regulates many amino acid-related genes in Gram-positive bacteria. T box riboswitch-mediated gene regulation was shown previously to occur at the level of transcription attenuation via structural rearrangements in the 5' untranslated (leader) region of the mRNA in response to binding of a specific uncharged tRNA. In this study, a novel group of isoleucyl-tRNA synthetase gene (ileS) T box leader sequences found in organisms of the phylum Actinobacteria was investigated. The Stem I domains of these RNAs lack several highly conserved elements that are essential for interaction with the tRNA ligand in other T box RNAs. Many of these RNAs were predicted to regulate gene expression at the level of translation initiation through tRNA-dependent stabilization of a helix that sequesters a sequence complementary to the Shine-Dalgarno (SD) sequence, thus freeing the SD sequence for ribosome binding and translation initiation. We demonstrated specific binding to the cognate tRNA(Ile) and tRNA(Ile)-dependent structural rearrangements consistent with regulation at the level of translation initiation, providing the first biochemical demonstration, to our knowledge, of translational regulation in a T box riboswitch.

High-throughput sequencing of T-cell receptors reveals a homogeneous repertoire of tumour-infiltrating lymphocytes in ovarian cancer.

The cellular adaptive immune system mounts a response to many solid tumours mediated by tumour-infiltrating T lymphocytes (TILs). Basic measurements of these TILs, including total count, show promise as prognostic markers for a variety of cancers, including ovarian and colorectal. In addition, recent therapeutic advances are thought to exploit this immune response to effectively fight melanoma, with promising studies showing efficacy in additional cancers. However, many of the basic properties of TILs are poorly understood, including specificity, clonality, and spatial heterogeneity of the T-cell response. We utilize deep sequencing of rearranged T-cell receptor beta (TCRB) genes to characterize the basic properties of TILs in ovarian carcinoma. Due to somatic rearrangement during T-cell development, the TCR beta chain sequence serves as a molecular tag for each T-cell clone. Using these sequence tags, we assess similarities and differences between infiltrating T cells in discretely sampled sections of large tumours and compare to T cells from peripheral blood. Within the limits of sensitivity of our assay, the TIL repertoires show strong similarity throughout each tumour and are distinct from the circulating T-cell repertoire. We conclude that the cellular adaptive immune response within ovarian carcinomas is spatially homogeneous and distinct from the T-cell compartment of peripheral blood.

Tumor-infiltrating lymphocytes in colorectal tumors display a diversity of T cell receptor sequences that differ from the T cells in adjacent mucosal tissue.

Tumors from colorectal cancer (CRC) are generally immunogenic and commonly infiltrated with T lymphocytes. However, the details of the adaptive immune reaction to these tumors are poorly understood. We have accrued both colon tumor samples and adjacent healthy mucosal samples from 15 CRC patients to study lymphocytes infiltrating these tissues. We apply a method for detailed sequencing of T-cell receptor (TCR) sequences from tumor-infiltrating lymphocytes (TILs) in CRC tumors at high throughput to probe T-cell clones in comparison with the TCRs from adjacent healthy mucosal tissue. In parallel, we captured TIL counts using standard immunohistochemistry. The variation in diversity of the TIL repertoire was far wider than the variation of T-cell clones in the healthy mucosa, and the oligoclonality was higher on average in the tumors. However, the diversity of the T-cell repertoire in both CRC tumors and healthy mucosa was on average 100-fold lower than in peripheral blood. Using the TCR sequences to identify and track clones between mucosal and tumor samples, we determined that the immune response in the tumor is different than in the adjacent mucosal tissue, and the number of shared clones is not dependent on distance between the samples. Together, these data imply that CRC tumors induce a specific adaptive immune response, but that this response differs widely in strength and breadth between patients.

Multimodal, broadly neutralizing antibodies against SARS-CoV-2 identified by high-throughput native pairing of BCRs from bulk B cells.

TruAB Discovery is an approach that integrates cellular immunology, high-throughput immunosequencing, bioinformatics, and computational biology in order to discover naturally occurring human antibodies for prophylactic or therapeutic use. We adapted our previously described pairSEQ technology to pair B cell receptor heavy and light chains of SARS-CoV-2 spike protein-binding antibodies derived from enriched antigen-specific memory B cells and bulk antibody-secreting cells. We identified approximately 60,000 productive, in-frame, paired antibody sequences, from which 2,093 antibodies were selected for functional evaluation based on abundance, isotype and patterns of somatic hypermutation. The exceptionally diverse antibodies included RBD-binders with broad neutralizing activity against SARS-CoV-2 variants, and S2-binders with broad specificity against betacoronaviruses and the ability to block membrane fusion. A subset of these RBD- and S2-binding antibodies demonstrated robust protection against challenge in hamster and mouse models. This high-throughput approach can accelerate discovery of diverse, multifunctional antibodies against any target of interest.

m6A modification of U6 snRNA modulates usage of two major classes of pre-mRNA 5' splice site.

Alternative splicing of messenger RNAs is associated with the evolution of developmentally complex eukaryotes. Splicing is mediated by the spliceosome, and docking of the pre-mRNA 5' splice site into the spliceosome active site depends upon pairing with the conserved ACAGA sequence of U6 snRNA. In some species, including humans, the central adenosine of the ACAGA box is modified by N6 methylation, but the role of this m6A modification is poorly understood. Here, we show that m6A modified U6 snRNA determines the accuracy and efficiency of splicing. We reveal that the conserved methyltransferase, FIONA1, is required for Arabidopsis U6 snRNA m6A modification. Arabidopsis fio1 mutants show disrupted patterns of splicing that can be explained by the sequence composition of 5' splice sites and cooperative roles for U5 and U6 snRNA in splice site selection. U6 snRNA m6A influences 3' splice site usage. We generalise these findings to reveal two major classes of 5' splice site in diverse eukaryotes, which display anti-correlated interaction potential with U5 snRNA loop 1 and the U6 snRNA ACAGA box. We conclude that U6 snRNA m6A modification contributes to the selection of degenerate 5' splice sites crucial to alternative splicing.

All the information necessary to build the proteins that perform the biological processes required for life is encoded in the DNA of an organism. Making these proteins requires the DNA sequence of a gene to be transcribed into a 'messenger RNA' (mRNA), which is then processed into a final, mature form. This blueprint is then translated to assemble the corresponding protein. When an mRNA is processed, segments of the sequence that do not code for protein are removed and the remaining coding sequences are joined together in the right order. An intricate molecular machine known as the spliceosome controls this mechanism by recognising the 'splice sites' where coding and non-coding sequences meet. Depending on external conditions, the spliceosome can 'pick-and-mix' the coding sequences to create different processed mRNAs (and therefore proteins) from a single gene. This alternative splicing mechanism is often used to regulate when certain biological processes take place based on environmental cues; for example, the splicing of genes which control the timing of plant flowering is sensitive to ambient temperatures. To investigate this mechanism, Parker et al. focused on Arabidopsis thaliana, a plant that blooms later when temperatures are low. This precise timing partly relies on a gene whose mRNA is efficiently spliced in the cold, resulting in an active form of its protein that blocks blooming. Parker et al. grew and screened many A. thaliana plants to find individuals that could flower early in the cold, in which splicing of this gene was disrupted. A mutant fitting these criteria was identified and subjected to further investigation, which revealed that it could not produce FIONA1. In non-mutant plants, this enzyme chemically modifies one of the components of the spliceosome, a small nuclear RNA known as U6. Parker et al found that there are two types of splice site ­ one more likely to interact with U6 and another that preferentially interacts with another small nuclear RNA, U5. When FIONA1 is inactive (such as in the mutant identified by Parker et al.), splice sites that tend to strongly interact with U5 are selected. However, when the enzyme is active, splice sites that tend to bind with the chemically modified U6 are used instead. Further work by Parker et al. showed that these two types of splice sites ('preferring' either U5 or U6) are found in equal proportions in the genomes of many species, including humans. This suggests that Parker et al. have uncovered an essential feature of how genomes are organised and splicing is controlled.

Chromosome evolution and the genetic basis of agronomically important traits in greater yam.

The nutrient-rich tubers of the greater yam, Dioscorea alata L., provide food and income security for millions of people around the world. Despite its global importance, however, greater yam remains an orphan crop. Here, we address this resource gap by presenting a highly contiguous chromosome-scale genome assembly of D. alata combined with a dense genetic map derived from African breeding populations. The genome sequence reveals an ancient allotetraploidization in the Dioscorea lineage, followed by extensive genome-wide reorganization. Using the genomic tools, we find quantitative trait loci for resistance to anthracnose, a damaging fungal pathogen of yam, and several tuber quality traits. Genomic analysis of breeding lines reveals both extensive inbreeding as well as regions of extensive heterozygosity that may represent interspecific introgression during domestication. These tools and insights will enable yam breeders to unlock the potential of this staple crop and take full advantage of its adaptability to varied environments.

Widespread premature transcription termination of Arabidopsis thaliana NLR genes by the spen protein FPA.

Genes involved in disease resistance are some of the fastest evolving and most diverse components of genomes. Large numbers of nucleotide-binding, leucine-rich repeat (NLR) genes are found in plant genomes and are required for disease resistance. However, NLRs can trigger autoimmunity, disrupt beneficial microbiota or reduce fitness. It is therefore crucial to understand how NLRs are controlled. Here, we show that the RNA-binding protein FPA mediates widespread premature cleavage and polyadenylation of NLR transcripts, thereby controlling their functional expression and impacting immunity. Using long-read Nanopore direct RNA sequencing, we resolved the complexity of NLR transcript processing and gene annotation. Our results uncover a co-transcriptional layer of NLR control with implications for understanding the regulatory and evolutionary dynamics of NLRs in the immune responses of plants.

Nanopore direct RNA sequencing maps the complexity of Arabidopsis mRNA processing and m6A modification.

Understanding genome organization and gene regulation requires insight into RNA transcription, processing and modification. We adapted nanopore direct RNA sequencing to examine RNA from a wild-type accession of the model plant Arabidopsis thaliana and a mutant defective in mRNA methylation (m6A). Here we show that m6A can be mapped in full-length mRNAs transcriptome-wide and reveal the combinatorial diversity of cap-associated transcription start sites, splicing events, poly(A) site choice and poly(A) tail length. Loss of m6A from 3' untranslated regions is associated with decreased relative transcript abundance and defective RNA 3' end formation. A functional consequence of disrupted m6A is a lengthening of the circadian period. We conclude that nanopore direct RNA sequencing can reveal the complexity of mRNA processing and modification in full-length single molecule reads. These findings can refine Arabidopsis genome annotation. Further, applying this approach to less well-studied species could transform our understanding of what their genomes encode.

A Public Database of Memory and Naive B-Cell Receptor Sequences.

The vast diversity of B-cell receptors (BCR) and secreted antibodies enables the recognition of, and response to, a wide range of epitopes, but this diversity has also limited our understanding of humoral immunity. We present a public database of more than 37 million unique BCR sequences from three healthy adult donors that is many fold deeper than any existing resource, together with a set of online tools designed to facilitate the visualization and analysis of the annotated data. We estimate the clonal diversity of the naive and memory B-cell repertoires of healthy individuals, and provide a set of examples that illustrate the utility of the database, including several views of the basic properties of immunoglobulin heavy chain sequences, such as rearrangement length, subunit usage, and somatic hypermutation positions and dynamics.

T-cell receptor sequencing reveals the clonal diversity and overlap of colonic effector and FOXP3+ T cells in ulcerative colitis.

BACKGROUND: FOXP3 regulatory T cell prevent inflammation but are paradoxically increased in ulcerative colitis (UC). Local T-cell activation has been hypothesized to account for increased FOXP3 expression in colon lamina propria (LP) T cells. METHODS: To see if human FOXP3 LP T cells are an activated fraction of otherwise FOXP3 effector T cells and explore their clonal diversity in health and disease, we deep sequenced clonally unique T-cell receptor hypervariable regions of FOXP3 and FOXP3CD4 T-cell subpopulations from inflamed versus noninflamed colon LP or mesenteric lymph nodes of patients with or without UC. RESULTS: The clonal diversity of each LP T-cell population was not different between patients with versus without UC. Repertoire overlap was only seen between a minority of FOXP3 and FOXP3 cells, including recently activated CD38 cells and Th17-like CD161 effector T cells, but this repertoire overlap was not different between patients with versus without UC and was no larger than the overlap between Helios and Helios FOXP3 cells. CONCLUSIONS: Thus, at steady state, only a minority of FOXP3, and particularly Helios, T cells share a T-cell receptor sequence with FOXP3 effector populations in the colon LP, even in UC, revealing distinct clonal origins for LP regulatory T cell and effector T cells in humans.

Annotation of pseudogenic gene segments by massively parallel sequencing of rearranged lymphocyte receptor loci.

BACKGROUND: The adaptive immune system generates a remarkable range of antigen-specific T-cell receptors (TCRs), allowing the recognition of a diverse set of antigens. Most of this diversity is encoded in the complementarity determining region 3 (CDR3) of the ß chain of the αß TCR, which is generated by somatic recombination of noncontiguous variable (V), diversity (D), and joining (J) gene segments. Deletion and non-templated insertion of nucleotides at the D-J and V-DJ junctions further increases diversity. Many of these gene segments are annotated as non-functional owing to defects in their primary sequence, the absence of motifs necessary for rearrangement, or chromosomal locations outside the TCR locus. METHODS: We sought to utilize a novel method, based on high-throughput sequencing of rearranged TCR genes in a large cohort of individuals, to evaluate the use of functional and non-functional alleles. We amplified and sequenced genomic DNA from the peripheral blood of 587 healthy volunteers using a multiplexed polymerase chain reaction assay that targets the variable region of the rearranged TCRß locus, and we determined the presence and the proportion of productive rearrangements for each TCRß V gene segment in each individual. We then used this information to annotate the functional status of TCRß V gene segments in this cohort. RESULTS: For most TCRß V gene segments, our method agrees with previously reported functional annotations. However, we identified novel non-functional alleles for several gene segments, some of which were used exclusively in our cohort to the detriment of reported functional alleles. We also saw that some gene segments reported to have both functional and non-functional alleles consistently behaved in our cohort as either functional or non-functional, suggesting that some reported alleles were not present in the population studied. CONCLUSIONS: In this proof-of-principle study, we used high-throughput sequencing of the TCRß locus of a large cohort of healthy volunteers to evaluate the use of functional and non-functional alleles of individual TCRß V gene segments. With some modifications, our method has the potential to be extended to gene segments in the α, γ, and δ TCR loci, as well as the genes encoding for B-cell receptor chains.

High-throughput pairing of T cell receptor α and ß sequences.

The T cell receptor (TCR) protein is a heterodimer composed of an α chain and a ß chain. TCR genes undergo somatic DNA rearrangements to generate the diversity of T cell binding specificities needed for effective immunity. Recently, high-throughput immunosequencing methods have been developed to profile the TCR α (TCRA) and TCR ß (TCRB) repertoires. However, these methods cannot determine which TCRA and TCRB chains combine to form a specific TCR, which is essential for many functional and therapeutic applications. We describe and validate a method called pairSEQ, which can leverage the diversity of TCR sequences to accurately pair hundreds of thousands of TCRA and TCRB sequences in a single experiment. Our TCR pairing method uses standard laboratory consumables and equipment without the need for single-cell technologies. We show that pairSEQ can be applied to T cells from both blood and solid tissues, such as tumors.

Detection of minimal residual disease in B lymphoblastic leukemia by high-throughput sequencing of IGH.

PURPOSE: High-throughput sequencing (HTS) of immunoglobulin heavy-chain genes (IGH) in unselected clinical samples for minimal residual disease (MRD) in B lymphoblastic leukemia (B-ALL) has not been tested. As current MRD-detecting methods such as flow cytometry or patient-specific qPCR are complex or difficult to standardize in the clinical laboratory, sequencing may enhance clinical prognostication. EXPERIMENTAL DESIGN: We sequenced IGH in paired pretreatment and day 29 post-treatment samples using residual material from consecutive, unselected samples from the Children's Oncology Group AALL0932 trial to measure MRD as compared with flow cytometry. We assessed the impact of ongoing recombination at IGH on MRD detection in post-treatment samples. Finally, we evaluated a subset of cases with discordant MRD results between flow cytometry and sequencing. RESULTS: We found clonal IGH rearrangements in 92 of 98 pretreatment patient samples. Furthermore, while ongoing recombination of IGH was evident, index clones typically prevailed in MRD-positive post-treatment samples, suggesting that clonal evolution at IGH does not contribute substantively to tumor fitness. MRD was detected by sequencing in all flow cytometry-positive cases with no false-negative results. In addition, in a subset of patients, MRD was detected by sequencing, but not by flow cytometry, including a fraction with MRD levels within the sensitivity of flow cytometry. We provide data that suggest that this discordance in some patients may be due to the phenotypic maturation of the transformed cell. CONCLUSION: Our results provide strong support for HTS of IGH to enhance clinical prognostication in B-ALL.

Estimating the ratio of CD4+ to CD8+ T cells using high-throughput sequence data.

Mature T cells express either CD8 or CD4, defining two physiologically distinct populations of T cells. CD8+ T cells, or killer T-cells, and CD4+ T cells, or helper T cells, effect different aspects of T cell mediated adaptive immunity. Currently, determining the ratio of CD4+ to CD8+ T cells requires flow cytometry or immunohistochemistry. The genomic T cell receptor locus is rearranged during T cell maturation, generating a highly variable T cell receptor locus in each mature T cell. As part of thymic maturation, T cells that will become CD4+ versus CD8+ are subjected to different selective pressures. In this study, we apply high-throughput next-generation sequencing to T cells from both a healthy cohort and a cohort with an autoimmune disease (multiple sclerosis) to identify sequence features in the variable CDR3 region of the rearranged T cell receptor gene that distinguish CD4+ from CD8+ T cells. We identify sequence features that differ between CD4+ and CD8+ T cells, including Variable gene usage and CDR3 region length. We implement a likelihood model to estimate relative proportions of CD4+ and CD8+ T cells using these features. Our model accurately estimates the proportion of CD4+ and CD8+ T cell sequences in samples from healthy and diseased immune systems, and simulations indicate that it can be applied to as few as 1000 T cell receptor sequences; we validate this model using in vitro mixtures of T cell sequences, and by comparing the results of our method to flow cytometry using peripheral blood samples. We believe our computational method for determining the CD4:CD8 ratio in T cell samples from sequence data will provide additional useful information for any samples on which high-throughput TCR sequencing is performed, potentially including some solid tumors.