High-throughput Sequencing Assay Research Articles

A Complementary Role of High Throughput Sequencing and Multiparameter Cytometry for Minimal Residual Disease (MRD) Detection in Chronic Lymphocytic Leukemia (CLL):an European Research Initiative (ERIC) Study

BACKGROUND. The detection of minimal residual disease (MRD) at the level of 0.01%/10-4 or above is a strong independent predictor of reduced progression-free (PFS) and overall survival (OS) in patients with CLL treated with chemoimmunotherapy. Although newer agents such as B-cell receptor pathway inhibitors can result in prolonged survival without achieving complete response, there remains a important role for MRD analysis in assessing therapeutic strategies aimed at disease eradication and cure. This is particularly important in front-line trials for fit patients which now require at least five years of follow-up if PFS is used as an endpoint. The feasibility of using MRD as a surrogate or intermediate endpoint for accelerated approval of new treatments is under review by regulatory agencies but further prospective validation is required. At the same time technology is rapidly evolving and high-throughput sequencing (HTS) technologies now detect MRD at the 0.0001%/10-6 level. It is therefore important to determine the most effective approaches for quantifying MRD that are compatible with previous studies but sufficiently sensitive for current treatments.AIMS. This collaborative project had two objectives. First, to identify the simplest and most flexible flow cytometry panel capable of detecting MRD at the 0.01%/10-4 or lower, that is compatible with published data and independent of instrument/reagent manufacturer. Second, to compare the flow cytometry approach with HTS analysis using the ClonoSEQ assay (Adaptive Biotechnologies, Seattle, WA).METHODS AND RESULTS. A core panel of antibodies for MRD detection was identified by testing an 8-marker combination in 52 samples (27 post-treatment and 25 dilution study) and re-analysing data with serial exclusion of single markers to determine redundancy. A 1-tube core panel of CD19, CD20, CD5, CD43, CD79b, and CD81 was identified and validated against the previously published 2-tube 6-marker and 4-tube 4-marker ERIC panels in 76 samples (19 post-treatment and 57 dilution study). The results showed good concordance (for log-transformed data above the LoQ, linearity=0.977, Pearson correlation co-efficient=0.983, average difference=0.026 log, 95% limit of agreement 0.20log) and a limit of detection of 0.001%/10-5 for the 1-tube core panel. Inter-operator variation was similar to CML MRD monitoring with both experienced operators, or inexperienced cytometrists after ~1 hour of specific education, achieving a 95% limit of agreement less than 0.3log in samples with MRD levels ranging from 0.0001 – 100%.The flow cytometry approach was compared with the ClonoSEQ HTS assay in 109 samples (21 dilution study and 88 post-treatment samples, complete data currently available on 13/88). The assay was applicable to the vast majority CLL patients, often with two clonal markers. There was 94% concordance at the 0.01% (10-4) threshold (15 samples with ≥0.01% CLL by both methods, 14 samples with <0.01% by both methods, 1 sample with 0.03% CLL by HTS and <0.003% CLL by flow cytometry, and 1 sample with 0.005% CLL by HTS and 0.012% by flow cytometry. HTS detected CLL IGH sequences in 22% (7/31) samples with no detectable CLL cells by flow cytometry (i.e. CLL level 0.0001-0.001%, 3/13 patient samples and 4/18 dilution samples). HTS demonstrated a relatively high variability in quantification, as seen in previous studies, but with a clear superiority in the limit of detection and good linearity (linearity=0.905, Pearson correlation co-efficient=0.870, average difference=0.078 log, 95% limit of agreement 1.5 log).CONCLUSIONS. The 1-tube 6-marker flow cytometry core panel is compatible with published studies, manufacturer-independent and flexible, providing directly quantitative results to 0.001%/10-5 without requiring a pre-treatment sample. HTS requires further work to standardise the quantitative analysis and prospective validation but the ClonoSEQ assay is applicable to >95% of CLL patients, does not require viable cells and is extremely sensitive, detecting residual disease in a significant proportion of cases with <0.01% CLL. The results indicate that flow cytometry and HTS are complementary technologies with a combined approach offering the most reliable way of quantifying CLL at the 0.01%/10-4 threshold while allowing higher sensitivity in clinical trials aimed at disease eradication. DisclosuresRawstron:Roche: Honoraria; Biogen Idec: Consultancy; Gilead: Consultancy, Honoraria; Abbvie: Honoraria; BD Biosciences: Intrasure reagent Patents & Royalties; Celgene: Honoraria; GSK: Honoraria. Williamson:Adaptive Biotechnologies: Employment, Equity Ownership. Sanders:Adaptive Biotechnologies: Employment, Equity Ownership. Robins:Adaptive Biotechnologies: Consultancy, Equity Ownership, Patents & Royalties. Hallek:Celgene: Honoraria, Research Funding; Mundipharma: Honoraria, Research Funding; Roche: Honoraria, Research Funding; Janssen: Honoraria, Research Funding; GSK: Honoraria; Gilead: Honoraria. Hillmen:Roche: Honoraria, Research Funding; GSK: Honoraria, Research Funding; Janssen: Honoraria, Research Funding; Pharmacyclics: Honoraria, Research Funding; Gilead: Honoraria, Research Funding.

Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2.

In comparative high-throughput sequencing assays, a fundamental task is the analysis of count data, such as read counts per gene in RNA-seq, for evidence of systematic changes across experimental conditions. Small replicate numbers, discreteness, large dynamic range and the presence of outliers require a suitable statistical approach. We present DESeq2, a method for differential analysis of count data, using shrinkage estimation for dispersions and fold changes to improve stability and interpretability of estimates. This enables a more quantitative analysis focused on the strength rather than the mere presence of differential expression. The DESeq2 package is available at http://www.bioconductor.org/packages/release/bioc/html/DESeq2.html.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-014-0550-8) contains supplementary material, which is available to authorized users.

Abstract 4260: Standardized decision support in next-generation sequencing (NGS) reports of somatic cancer variants

Abstract Introduction: Of hundreds to thousands of somatic mutations that exist in each cancer genome, a large number are unique and non-recurrent variants. Many such variants occur in tumor genes that have well-established biological and clinical relevance and are putative targets of molecular therapy, however, most variants are still of unknown significance. Prioritizing and reporting genetic variants identified via NGS technologies remains a major challenge. Methods: Based on a systematic framework for cancer variant annotation/prioritization, we propose a structured molecular pathology report using standardized terminology in order to best inform oncology clinical practice. In parallel with the variant interpretation pipeline, we developed a comprehensive knowledge database (Kdb) that integrates tumor types, genes, variants, response/resistance patterns to approved and experimental agents and PubMed identifiers. Results: As physicians rely heavily on limited subsets of data and frequently apply rules based on strength of evidence in the decision-making process, NGS results are presented in a very discrete manner, avoiding detailed descriptions of each genomic aberration. Reportable variants are grouped in three categories: (i) actionable, which support treatment recommendation, enrollment in clinical trials and/or have prognostic or diagnostic implications; (ii) biologically relevant but not clearly actionable, including novel variants lacking functional preclinical validation in known cancer genes; and (iii) variants of unknown significance. Predictive associations are reported according to tumor type and classified in a hierarchical way based on the strength of evidence: (i) late trials; (ii) early trials; (iii) case reports; and (iv) preclinical data being explored in clinical trials. In our internal Kdb, tumor types with the largest number of emerging predictive associations are acute leukemias, non-small cell lung cancer, brain tumors, melanoma, colorectal, ovarian and breast cancers, comprising more than 100 unique gene - drug interactions. Further interpretation of gene variants requires a team with strong background in cancer biology, careful evaluation of the published literature and ability to match a patient's tumor genotype to clinical trials in the context of multiple aberrations. Conclusions: With large amounts of data being generated as high throughput sequencing assays enter the clinical realm, there is a growing need to better communicate relevant findings in a timely manner while remaining cognizant of the potential consequences of misuse or overinterpretation of genomic information. The scientific knowledge on actionable events should be presented the report, so that physicians can make evidence-based decisions and patients may benefit from treatment with matched targeted agents. We hope that our experience in this process will help other institutions implement similar programs. Citation Format: Rodrigo Dienstmann, Fei Dong, Darrell Borger, Dora Dias Santagata, Leif W. Ellisen, Long P. Le, A. John Iafrate. Standardized decision support in next-generation sequencing (NGS) reports of somatic cancer variants. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 4260. doi:10.1158/1538-7445.AM2014-4260

N-terminal domains of human DNA polymerase lambda promote primer realignment during translesion DNA synthesis.

The X-family DNA polymerases λ (Polλ) and β (Polβ) possess similar 5'-2-deoxyribose-5-phosphate lyase (dRPase) and polymerase domains. Besides these domains, Polλ also possesses a BRCA1 C-terminal (BRCT) domain and a proline-rich domain at its N terminus. However, it is unclear how these non-enzymatic domains contribute to the unique biological functions of Polλ. Here, we used primer extension assays and a newly developed high-throughput short oligonucleotide sequencing assay (HT-SOSA) to compare the efficiency of lesion bypass and fidelity of human Polβ, Polλ and two N-terminal deletion constructs of Polλ during the bypass of either an abasic site or an 8-oxo-7,8-dihydro-2'-deoxyguanosine (8-oxodG) lesion. We demonstrate that the BRCT domain of Polλ enhances the efficiency of abasic site bypass by approximately 1.6-fold. In contrast, deletion of the N-terminal domains of Polλ did not affect the efficiency of 8-oxodG bypass relative to nucleotide incorporations opposite undamaged dG. HT-SOSA analysis demonstrated that Polλ and Polβ preferentially generated -1 or -2 frameshift mutations when bypassing an abasic site and the single or double base deletion frequency was highly sequence dependent. Interestingly, the BRCT and proline-rich domains of Polλ cooperatively promoted the generation of -2 frameshift mutations when the abasic site was situated within a sequence context that was susceptible to homology-driven primer realignment. Furthermore, both N-terminal domains of Polλ increased the generation of -1 frameshift mutations during 8-oxodG bypass and influenced the frequency of substitution mutations produced by Polλ opposite the 8-oxodG lesion. Overall, our data support a model wherein the BRCT and proline-rich domains of Polλ act cooperatively to promote primer/template realignment between DNA strands of limited sequence homology. This function of the N-terminal domains may facilitate the role of Polλ as a gap-filling polymerase within the non-homologous end joining pathway.

Mutagenic potential of 8-oxo-7,8-dihydro-2'-deoxyguanosine bypass catalyzed by human Y-family DNA polymerases.

One of the most common lesions induced by oxidative DNA damage is 8-oxo-7,8-dihydro-2′-deoxyguanosine (8-oxodG). Replicative DNA polymerases poorly traverse this highly mutagenic lesion, suggesting that the replication fork may switch to a polymerase specialized for translesion DNA synthesis (TLS) to catalyze 8-oxodG bypass in vivo. Here, we systematically compared the 8-oxodG bypass efficiencies and fidelities of the TLS-specialized, human Y-family DNA polymerases eta (hPolη), iota (hPolι), kappa (hPolκ), and Rev1 (hRev1) either alone or in combination. Primer extension assays revealed that the times required for hPolη, hRev1, hPolκ, and hPolι to bypass 50% of the 8-oxodG lesions encountered (t50bypass) were 0.58, 0.86, 108, and 670 s, respectively. Although hRev1 bypassed 8-oxodG efficiently, hRev1 failed to catalyze the extension step of TLS, and a second polymerase was required to extend the lesion bypass products. A high-throughput short oligonucleotide sequencing assay (HT-SOSA) was used to quantify the types and frequencies of incorporation errors produced by the human Y-family DNA polymerases at and near the 8-oxodG site. Although hPolη bypassed 8-oxodG most efficiently, hPolη correctly incorporated dCTP opposite 8-oxodG within only 54.5% of the sequences analyzed. In contrast, hPolι bypassed the lesion least efficiently but correctly incorporated dCTP at a frequency of 65.8% opposite the lesion. The combination of hRev1 and hPolκ was most accurate opposite 8-oxodG (92.3%), whereas hPolκ alone was the least accurate (18.5%). The t50bypass value and correct dCTP incorporation frequency in the presence of an equal molar concentration of all four Y-family enzymes were 0.60 s and 43.5%, respectively. These values are most similar to those of hPolη alone, suggesting that hPolη outcompetes the other three Y-family polymerases to catalyze 8-oxodG bypass in vitro and possibly in vivo.

STATegra EMS: an Experiment Management System for complex next-generation omics experiments.

High-throughput sequencing assays are now routinely used to study different aspects of genome organization. As decreasing costs and widespread availability of sequencing enable more laboratories to use sequencing assays in their research projects, the number of samples and replicates in these experiments can quickly grow to several dozens of samples and thus require standardized annotation, storage and management of preprocessing steps. As a part of the STATegra project, we have developed an Experiment Management System (EMS) for high throughput omics data that supports different types of sequencing-based assays such as RNA-seq, ChIP-seq, Methyl-seq, etc, as well as proteomics and metabolomics data. The STATegra EMS provides metadata annotation of experimental design, samples and processing pipelines, as well as storage of different types of data files, from raw data to ready-to-use measurements. The system has been developed to provide research laboratories with a freely-available, integrated system that offers a simple and effective way for experiment annotation and tracking of analysis procedures.

Compression of Structured High-Throughput Sequencing Data

Large biological datasets are being produced at a rapid pace and create substantial storage challenges, particularly in the domain of high-throughput sequencing (HTS). Most approaches currently used to store HTS data are either unable to quickly adapt to the requirements of new sequencing or analysis methods (because they do not support schema evolution), or fail to provide state of the art compression of the datasets. We have devised new approaches to store HTS data that support seamless data schema evolution and compress datasets substantially better than existing approaches. Building on these new approaches, we discuss and demonstrate how a multi-tier data organization can dramatically reduce the storage, computational and network burden of collecting, analyzing, and archiving large sequencing datasets. For instance, we show that spliced RNA-Seq alignments can be stored in less than 4% the size of a BAM file with perfect data fidelity. Compared to the previous compression state of the art, these methods reduce dataset size more than 40% when storing exome, gene expression or DNA methylation datasets. The approaches have been integrated in a comprehensive suite of software tools (http://goby.campagnelab.org) that support common analyses for a range of high-throughput sequencing assays.

Targeted Next-Generation Sequencing Of 2,761 Genes Detects Copy Number States and Molecular Mutations In a Single Approach

IntroductionIn acute myeloid leukemia (AML), the karyotype and molecular mutation profile are the strongest determinants for prognosis and biological subclassification. Yet, diagnostic analyses rely on chromosome banding technique and sequencing of a constantly growing number of genes. AimsIn an era of novel high-throughput sequencing assays becoming viable options for diagnostic implementation we aimed to evaluate whether the application of targeted exome sequencing can reliably identify copy number states and molecular mutations in a single-step procedure. Patients and MethodsThe pilot cohort included four AML cases with a complex karyotype with known chromosomal alterations as detected by chromosome banding analysis, 24-color FISH and array CGH (12x270K microarrays, NimbleGen, Madison, WI). The size of the aberrant clone was determined by suitable probes using interphase-FISH on bone marrow smears. For sequencing analysis genomic DNA was extracted from mononuclear cells and 50 ng were processed using the TruSight Rapid Capture kit (Illumina, San Diego, CA). Sequencing was performed on a MiSeq instrument using the 2x150 bp paired-end read chemistry targeting a subset of the human exome (2,761 genes; 37,366 exons). This exome enrichment library contained >50,000 probes (7.75 Mb) focusing on disease-causing variants in specific inherited conditions (Illumina). Data analysis was performed applying default settings of the on-board MiSeq Reporter Software version 2.2.29 using the Burrows-Wheeler Aligner to align the reads against the hg19 reference genome. Further processing to delineate copy number states was performed using the ExomeCNV package. ResultsEach patient was analyzed in a single MiSeq run and in median 22,022,240 (range 19,233,134 - 23,507,016) reads were generated. The median coverage per target region was in the range of 74-186 reads. Coverage uniformity was assessed according to the manufacturer's recommendations. Over 98% of bases were covered at 0.12X mean coverage for each sample. Next, two data analysis pipelines were triggered, i.e. copy number states and mutation analysis. With respect to copy number alterations (CNA), in total 65 CNA were detected by chromosome banding analysis/array CGH. Of these, 21 were gains, 44 were losses. The size of the deletions ranged between 378,377 and 141,048,720 bp (median 10,731,680 bp), the size of the gains ranged between 281,608 and 46,404,876 bp (median 4,947,125 bp), respectively.In total, 63/65 (96.9%) copy number alterations were correctly identified by targeted exome sequencing. The NGS assay was able to detect copy number alterations that were present in only 23% of cells as determined by interphase-FISH. In detail, one of the deletions was homozygous with a larger deletion on the long arm of chromosome 17 (size: 1,070,162 bp) and a small intragenic deletion within the NF1 gene. This homozygous deletion was detected by array-CGH and by exome sequencing. Interestingly, the higher resolution of the exome sequencing assay in this area enabled the exact localization (exons 37 to 58) and size determination (78,415 bp) of the deletion. Overall, only 2 gains escaped detection. These were two small gained regions on a highly rearranged chr. 19.Secondly, with respect to mutation analysis, the same assay detected 19, 20, 21 and 28 mutations in the four analyzed patients. This pipeline took only putative variants into account that were not present in the control sample, were having a coverage ≥30 reads with a mutation load ≥10%, and had a confirmed COSMIC mutation entry (v66). 12/2,761 (0.4%) genes harbored mutations in at least 2/4 patients. This included genes known to be involved in leukemogenesis. TP53 mutations were detected in all four cases and all were confirmed by Sanger sequencing. ConclusionsA targeted exome sequencing assay allowed to robustly assess copy number states in AML at diagnosis at a resolution greater than current conventional array CGH analyses. Moreover, exome sequencing read data also can be used to delineate mutation profiles. Thus, this workflow enabled to call gene mutations and copy number states in a single assay and is a promising option for a routine diagnostics assay in the future. The gene panel has to be further optimized by adding genes known to be mutated in hematological malignancies. More data is necessary to precisely determine the detection limit and to optimize software tools for a routine use. Disclosures:Kohlmann:MLL Munich Leukemia Laboratory: Employment. Roller:MLL Munich Leukemia Laboratory: Employment. Weissmann:MLL Munich Leukemia Laboratory: Employment. Kuznia:MLL Munich Leukemia Laboratory: Employment. Zenger:MLL Munich Leukemia Laboratory: Employment. Kern:MLL Munich Leukemia Laboratory: Employment, Equity Ownership. Schnittger:MLL Munich Leukemia Laboratory: Employment, Equity Ownership. Haferlach:MLL Munich Leukemia Laboratory: Employment, Equity Ownership. Haferlach:MLL Munich Leukemia Laboratory: Employment, Equity Ownership.

Zinc finger protein ZBTB20 promotes toll-like receptor-triggered innate immune responses by repressing IκBα gene transcription

Toll-like receptor (TLR) signaling is critical in innate response against invading pathogens. However, the molecular mechanisms for full activation of TLR-triggered innate immunity need to be fully elucidated. The broad complex tramtrack bric-a-brac/poxvirus and zinc finger (BTB/POZ) family is a class of transcription factors involved in many biological processes. However, few BTB/POZ proteins were reported to function in innate immune response. Zinc finger and BTB domain-containing 20 (ZBTB20), a member of BTB/POZ family, functions in neurogenesis and represses α-fetoprotein gene transcription in liver. However, the immunological functions of ZBTB20 remain unknown. Here, we found that myeloid cell-specific ZBTB20 KO mice were resistant to endotoxin shock and Escherichia coli-caused sepsis. ZBTB20 deficiency attenuated TLR-triggered production of proinflammatory cytokines and type I IFN in macrophages, which attributed to higher abundance of IκBα protein and impaired activity of NF-κB. Furthermore, ChIP and next generation high-throughput DNA sequencing assay showed that ZBTB20 specifically bound to IκBα gene promoter (+1 to +60 region) after TLR activation. ZBTB20 could inhibit IκBα gene transcription, govern IκBα protein expression, and then promote NF-κB activation. Therefore, transcriptional repressor ZBTB20 is needed to promote full activation of TLR signaling and TLR-triggered innate immune response by selectively suppressing the suppressor IκBα gene transcription.

Desktop Transcriptome Sequencing From Archival Tissue to Identify Clinically Relevant Translocations

Somatic mutations, often translocations or single nucleotide variations, are pathognomonic for certain types of cancers and are increasingly of clinical importance for diagnosis and prediction of response to therapy. Conventional clinical assays only evaluate 1 mutation at a time, and targeted tests are often constrained to identify only the most common mutations. Genome-wide or transcriptome-wide high-throughput sequencing (HTS) of clinical samples offers an opportunity to evaluate for all clinically significant mutations with a single test. Recently a "desktop version" of HTS has become available, but most of the experience to date is based on data obtained from high-quality DNA from frozen specimens. In this study, we demonstrate, as a proof of principle, that translocations in sarcomas can be diagnosed from formalin-fixed paraffin-embedded (FFPE) tissue with desktop HTS. Using the first generation MiSeq platform, full transcriptome sequencing was performed on FFPE material from archival blocks of 3 synovial sarcomas, 3 myxoid liposarcomas, 2 Ewing sarcomas, and 1 clear cell sarcoma. Mapping the reads to the "sarcomatome" (all known 83 genes involved in translocations and mutations in sarcoma) and using a novel algorithm for ranking fusion candidates, the pathognomonic fusions and the exact breakpoints were identified in all cases of synovial sarcoma, myxoid liposarcoma, and clear cell sarcoma. The Ewing sarcoma fusion gene was detectable in FFPE material only with a sequencing platform that generates greater sequencing depth. The results show that a single transcriptome HTS assay, from FFPE, has the potential to replace conventional molecular diagnostic techniques for the evaluation of clinically relevant mutations in cancer.

Minimizing PCR bias in high-throughput adaptive immune repertoire sequencing assays (P3389)

Abstract The human adaptive immune system is a primary defense against the vast world of potential pathogens. To counter this, B and T cells rearrange a diverse repertoire of surface receptors. This system functions in part by clonal expansion, so measuring changes in the abundance of both unique and total clones is necessary to understand the immune response. We developed immune profiling assays to amplify and sequence rearranged receptors. To efficiently capture receptor diversity we use a multiplex PCR approach. However, small variations in annealing kinetics can have large impacts on amplification efficiency causing PCR amplification bias (PAB). To address this issue, we developed a method to quantify and correct PAB introduced by multiplex PCR. We demonstrate this method using our T cell receptor γ (TCRG) assay. Our method uses a quantitated, synthetic immune repertoire standard consisting of every possible Variable(V) and Joining(J) gene combination. Following multiplex TCRG PCR amplification, we assessed the difference between the synthetic VJ composition before and after amplification and quantified PAB. We adjusted primer concentrations and redesigned primers to generate an unbiased multiplex PCR assay. Using this method we reduced PAB for the TCRG assay from a 340-fold to a 4.5-fold dynamic range. Now the abundance of receptor sequences is representative of the frequency of each receptor chain in the input sample, allowing for accurate quantification of the TCRG repertoire.

A high-throughput and quantitative method to assess the mutagenic potential of translesion DNA synthesis

Cellular genomes are constantly damaged by endogenous and exogenous agents that covalently and structurally modify DNA to produce DNA lesions. Although most lesions are mended by various DNA repair pathways in vivo, a significant number of damage sites persist during genomic replication. Our understanding of the mutagenic outcomes derived from these unrepaired DNA lesions has been hindered by the low throughput of existing sequencing methods. Therefore, we have developed a cost-effective high-throughput short oligonucleotide sequencing assay that uses next-generation DNA sequencing technology for the assessment of the mutagenic profiles of translesion DNA synthesis catalyzed by any error-prone DNA polymerase. The vast amount of sequencing data produced were aligned and quantified by using our novel software. As an example, the high-throughput short oligonucleotide sequencing assay was used to analyze the types and frequencies of mutations upstream, downstream and at a site-specifically placed cis–syn thymidine–thymidine dimer generated individually by three lesion-bypass human Y-family DNA polymerases.

Modeling Long Chromatin Fibers based on In-Vivo Nucleosome Positioning Maps

In the nucleus of eukaryotes DNA is wrapped around histone proteins forming so-called nucleosomes, the basic unit of chromatin. This packaging controls DNA accessibility, thus influences directly gene expression, DNA repair and recombination. The precise structure of chromatin in-vivo is still under debate, because experimental probing of chromatin in the nucleus remains difficult and evidence for proposed models is rare. The majority of theoretical models imply a uniform nucleosome repeat length. However, in-vivo nucleosomes exhibit more irregular patterns. It has been shown that the internucleosomal distance has a great impact on the local and global structure of the chromatin arrangement.Here we investigate the influence of nucleosome positioning on the three-dimensional chromatin structure by modeling chromatin fibers comprising more than 1000 nucleosomes. The spacing of nucleosomes was based on in-vivo nucleosome positioning maps from mouse embryonic stem cells and differentiated cells derived from these. We determined non-overlapping nucleosome configurations from MNase assisted high throughput sequencing assays by a simulation procedure. In the next step, we derived a three-dimensional coarse-grained chromatin model from nucleosome positioning data, including fundamental physical properties like elasticity, electrostatics and nucleosome interactions. To improve the understanding of how structural differences and function are linked, we simulate models of gene-clusters at different stages of cell differentiation.

Detection of Multiple Myeloma Cells in Peripheral Blood Using High-Throughput Sequencing Assay

AbstractAbstract 321 Background:Multiple myeloma (MM) is characterized by the presence of monoclonal protein (M-protein) in serum and/or urine, clonal plasma cell accumulation in bone marrow (BM), and related organ or tissue impairment. MM patients are monitored during therapy and posttherapy using immunoglobulin, M-protein and free light chain assays. Assessing pathological myeloma cells using flow cytometry and RT-PCR has been shown to have superior prognostic value. However, the sensitivity of these techniques has generally limited their use to assessment of BM. In order to determine whether myeloma minimal residual disease could be detected in peripheral blood (PB), we assessed a cohort of MM patients using a sequencing based platform, LymphoSIGHT, with a sensitivity of 1 cancer cell per million leukocytes. Methods:We obtained from 4 sources (UCSF, NYU, Washington Univ, commercial source) pairs of BM and PB samples from 60 MM patients at different points of disease. BM samples were used to identify the clonal MM sequence and detection of that sequence in the PB was assessed. For some patients BM/PB sample pairs were obtained from >1 time point resulting in a total of 78 pairs. In addition blood and bone marrow plasma samples were available for 44 and 6 patients, respectively, to assess presence of the myeloma clonotype in cell free DNA. Altogether there were 206 samples. BM samples were either available as BM mononuclear cells (BMMC) or bead enriched CD138+ cells. Using universal primer sets, we amplified IgH@ variable (V), diversity (D), and joining (J) gene segments from genomic DNA and/or RNA samples, the incomplete IgH rearrangement (DJ), and IgK from genomic DNA. Amplified products were sequenced to obtain >1 million reads and analyzed using standardized algorithms for clonotype determination. Myeloma-specific IgH, IgK, and DJ clonotypes were identified for each patient based on their high frequency in BM samples. The presence of the myeloma clonotype was then assessed in all PBMC (DNA and RNA), BM plasma (DNA), and PB plasma (DNA) samples using the same IgH and in some samples using the IgK sequencing assays. A quantitative and standardized measure of clone level among all leukocytes in each PB or BM sample was determined using internal reference DNA. Here we describe data on 46/60 patients; data from all 60 patients will be presented. Results:In BM samples, we detected the myeloma clonal rearrangement of at least one receptor (“calibrating receptor”) in 34/46 (74%) of MM patients (Table 1). The calibration rate varied by receptor, with 30/46 (65%) patients calibrating with IgH, 14/43 (33%) with IgH DJ, and 22/43 (51%) with IgK (Table 1). Identification of myeloma-specific clonal rearrangement is based on presence at high frequency and may not occur in samples from patients with low disease load (e.g., post-treatment). Of the 12 non-calibrating patients, only 3 had high disease load. The myeloma clonotype that was identified in the BM was also detected in PBMC in 22/30 (73%) and 28/30 (93%) patients with the DNA and RNA IgH analysis, respectively (Table 2). IgK DNA analysis showed the presence of the myeloma clonotype in 9/10 PBMC samples, all of which were concordant with IgH results. The myeloma clonotype that was identified in the BM was also detected in the cell-free BM and PB samples in 5/5 and 7/11 patients, respectively, using the IgH DNA assay. The evaluation of blood plasma and PBMC were at times complementary in detecting the myeloma. Conclusions:Results from the application of a high-throughput sequencing method for detection of myeloma-specific clonotypes in 46 MM patients are shown. A clonal rearrangement was detected in 74% of MM BM samples. Importantly, 93% of peripheral blood samples from 30 patients showed evidence of circulating myeloma in PBMC. Analysis of BM and PB samples from 14 additional MM patients as well as association of the level of myeloma in PBMC and BM with clinical measures is ongoing.Table 1Calibration rates for IgH and IgK receptors in MM patients.Sample sourceIgH (complete)IgKIgH (incomplete)Any receptorUCSF8/126/123/129/12NYU15/2113/2110/2117/21Commercial source7/133/101/108/13Total30/4622/4314/4334/46Table 2Detection of circulating myeloma material in MM patients.Sample sourceIgH clone present in blood DNAIgH clone present in blood RNAIgK clone present in blood DNAUCSF4/87/82/3NYU13/1515/157/7Commercial source5/76/7N/ATotal22/3028/309/10 Disclosures:Faham:Sequenta, Inc.: Employment, Equity Ownership, Research Funding. Klinger:Sequenta, Inc.: Employment, Equity Ownership, Research Funding.

Forkhead box protein p1 is a transcriptional repressor of immune signaling in the CNS: implications for transcriptional dysregulation in Huntington disease

Forkhead box protein p1 (Foxp1), a transcription factor showing highly enriched expression in the striatum, has been implicated in central nervous system (CNS) development, but its role in the mature brain is unknown. In order to ascertain functional roles for Foxp1 in the CNS, we have identified gene targets for Foxp1 both in vitro and in vivo using genome-wide expression microarrays and chromatin-immunoprecipitation followed by high-throughput sequencing (ChIP-seq) assays. We found that mouse Foxp1 overexpression in striatal cells elicited expression changes of genes related to immune signaling, transcriptional regulation and a manually curated Huntington's disease (HD)-signaling pathway. Similar results were found when the gene expression data set was integrated with Foxp1-binding data determined from ChIP-seq analysis. In vivo lentiviral-mediated overexpression of human FOXP1 in the context of mutant huntingtin (Htt) protein resulted in a robust downregulation of glial cell-associated, immune genes, including those encoding a variety of cytokines and chemokines. Furthermore, Foxp1-induced expression changes were significantly negatively correlated with those changes elicited by mutant Htt protein in several different HD mouse models, and most significantly in post-mortem caudate from human HD subjects. We finally show that Foxp1 interacts with mutant Htt protein in mouse brain and is present in nuclear Htt aggregates in the striatum of R6/1 transgenic mice. These findings implicate Foxp1 as a key repressor of immune signaling in the CNS and suggest that the loss of Foxp1-mediated gene regulation in HD contributes to the immune dysfunction in this disease. We further suggest that Foxp1-regulated pathways might be important mediators of neuronal-glial cell communication.

Differential expression analysis for sequence count data.

High-throughput sequencing assays such as RNA-Seq, ChIP-Seq or barcode counting provide quantitative readouts in the form of count data. To infer differential signal in such data correctly and with good statistical power, estimation of data variability throughout the dynamic range and a suitable error model are required. We propose a method based on the negative binomial distribution, with variance and mean linked by local regression and present an implementation, DESeq, as an R/Bioconductor package.

High-Throughput Gene Sequencing Assay Development for Hereditary Nonpolyposis Colon Cancer

Hereditary nonpolyposis colorectal cancer (HNPCC) is the most common hereditary colon cancer syndrome and is responsible for as many as 10% of all colorectal cancers. Hereditary nonpolyposis colorectal cancer is autosomally dominant with a prevalence of 1 in 200-2000 and exhibits incomplete penetrance. Affected individuals have an approximately 70% lifetime risk of colon cancer with a mean age of onset of 44 years and an approximately 40% lifetime risk of endometrial cancer. At least 5 mismatch repair genes (MLH1, MSH2, MSH6, PMS1, PMS2) have been implicated in HNPCC; however, no predominant mutations were found in these genes. Mutation detection by direct sequencing has proven to be the most sensitive method. We have developed high-throughput full-length sequencing assays of the MLH1, MSH2, and MSH6 genes. These 3 genes account for approximately 90% of all germline mutations found in HNPCC. In our assays, 19 exons of MLH1, 16 exons of MSH2, 10 exons of MSH6, and the adjacent splice sites were amplified using polymerase chain reaction and loaded onto a capillary sequencing machine. Results were analyzed using sequence analysis software and stored in a relational database. Our assay method was validated using 15 affected patients and normal controls. It is anticipated that our high-throughput assay technique will provide accurate diagnoses for patients at risk for HNPCC and thereby facilitate early curative intervention.

A semi-automated system for analysis and storage of SNPs.

The discovery of single nucleotide polymorphisms ( SNPs) is currently pursued with a tremendous effort. SNPs represent a rich source for molecular markers, since estimations predict six to seven million of these DNA variations in the human genome. A subset of these genetic variants is thought to have a pervasive impact on modern medicine, be it for the elucidation of differential pharmacological response or for the facilitated identification of genes involved in monogenetic and complex human diseases. Here we describe the overall process that leads to the set up of a SNP database. We describe a high-throughput sequencing assay for SNP discovery, automation of the dataflow from the DNA sequencer to the SNP analysis, and the tools to facilitate it. At the end of the process, a web-accessible interface collects the SNP information, which is processed in order to be written into the SNP database and to be available for end users who would like to select appropriate SNPs for their special screening needs.