Robust detection of translocations in lymphoma FFPE samples using targeted locus capture-based sequencing.

In routine diagnostic pathology, cancer biopsies are preserved by formalin-fixed, paraffin-embedding (FFPE) procedures for examination of (intra-) cellular morphology. Such procedures inadvertently induce DNA fragmentation, which compromises sequencing-based analyses of chromosomal rearrangements. Yet, rearrangements drive many types of hematolymphoid malignancies and solid tumors, and their manifestation is instructive for diagnosis, prognosis, and treatment. Here, we present FFPE-targeted locus capture (FFPE-TLC) for targeted sequencing of proximity-ligation products formed in FFPE tissue blocks, and PLIER, a computational framework that allows automated identification and characterization of rearrangements involving selected, clinically relevant, loci. FFPE-TLC, blindly applied to 149 lymphoma and control FFPE samples, identifies the known and previously uncharacterized rearrangement partners. It outperforms fluorescence in situ hybridization (FISH) in sensitivity and specificity, and shows clear advantages over standard capture-NGS methods, finding rearrangements involving repetitive sequences which they typically miss. FFPE-TLC is therefore a powerful clinical diagnostics tool for accurate targeted rearrangement detection in FFPE specimens.

TEG001 Insert Integrity from Vector Producer Cells until Medicinal Product.

Despite extensive usage of gene therapy medicinal products (GTMPs) in clinical studies and recent approval of chimeric antigen receptor (CAR) T cell therapy, little information has been made available on the precise molecular characterization and possible variations in terms of insert integrity and vector copy numbers of different GTMPs during the complete production chain. Within this context, we characterize αßT cells engineered to express a defined γδT cell engineered to express a defined γδT receptor (TEG) currently used in a first-in-human clinical study (NTR6541). Utilizing targeted locus amplification in combination with next generation sequencing for the vector producer clone and TEG001 products, we report on five single-nucleotide variants and nine intact vector copies integrated in the producer clone. The vector copy number in TEG001 cells was on average a factor 0.72 (SD 0.11) below that of the producer cell clone. All nucleotide variants were transferred to TEG001 without having an effect on cellular proliferation during extensive in vitro culture. Based on an environmental risk assessment of the five nucleotide variants present in the non-coding viral region of the TEG001 insert, there was no altered environmental impact of TEG001 cells. We conclude that TEG001 cells do not have an increased risk for malignant transformation in vivo.

Enhanced CHO Clone Screening: Application of Targeted Locus Amplification and Next-Generation Sequencing Technologies for Cell Line Development.

Early analytical clone screening is important during Chinese hamster ovary (CHO) cell line development of biotherapeutic proteins to select a clonally derived cell line with most favorable stability and product quality. Sensitive sequence confirmation methods using mass spectrometry have limitations in throughput and turnaround time. Next-generation sequencing (NGS) technologies emerged as alternatives for CHO clone analytics. We report an efficient NGS workflow applying the targeted locus amplification (TLA) strategy for genomic screening of antibody expressing CHO clones. In contrast to previously reported RNA sequencing approaches, TLA allows for targeted sequencing of genomic integrated transgenic DNA without prior locus information, robust detection of single-nucleotide variants (SNVs) and transgenic rearrangements. During clone selection, TLA/NGS revealed CHO clones with high-level SNVs within the antibody gene and we report in another case the utility of TLA/NGS to identify rearrangements at transgenic DNA level. We also determined detection limits for SNVs calling and the potential to identify clone contaminations by TLA/NGS. TLA/NGS also allows to identify genetically identical clones. In summary, we demonstrate that TLA/NGS is a robust screening method useful for routine clone analytics during cell line development with the potential to process up to 24 CHO clones in less than 7 workdays.

Efficient mapping of transgene integration sites and local structural changes in Cre transgenic mice using targeted locus amplification.

Cre/LoxP technology is widely used in the field of mouse genetics for spatial and/or temporal regulation of gene function. For Cre lines generated via pronuclear microinjection of a Cre transgene construct, the integration site is random and in most cases not known. Integration of a transgene can disrupt an endogenous gene, potentially interfering with interpretation of the phenotype. In addition, knowledge of where the transgene is integrated is important for planning of crosses between animals carrying a conditional allele and a given Cre allele in case the alleles are on the same chromosome. We have used targeted locus amplification (TLA) to efficiently map the transgene location in seven previously published Cre and CreERT2 transgenic lines. In all lines, transgene insertion was associated with structural changes of variable complexity, illustrating the importance of testing for rearrangements around the integration site. In all seven lines the exact integration site and breakpoint sequences were identified. Our methods, data and genotyping assays can be used as a resource for the mouse community and our results illustrate the power of the TLA method to not only efficiently map the integration site of any transgene, but also provide additional information regarding the transgene integration events.

Molecular dissection of germline chromothripsis in a developmental context using patient-derived iPS cells.

BACKGROUND: Germline chromothripsis causes complex genomic rearrangements that are likely to affect multiple genes and their regulatory contexts. The contribution of individual rearrangements and affected genes to the phenotypes of patients with complex germline genomic rearrangements is generally unknown. METHODS: To dissect the impact of germline chromothripsis in a relevant developmental context, we performed trio-based RNA expression analysis on blood cells, induced pluripotent stem cells (iPSCs), and iPSC-derived neuronal cells from a patient with de novo germline chromothripsis and both healthy parents. In addition, Hi-C and 4C-seq experiments were performed to determine the effects of the genomic rearrangements on transcription regulation of genes in the proximity of the breakpoint junctions. RESULTS: Sixty-seven genes are located within 1 Mb of the complex chromothripsis rearrangements involving 17 breakpoints on four chromosomes. We find that three of these genes (FOXP1, DPYD, and TWIST1) are both associated with developmental disorders and differentially expressed in the patient. Interestingly, the effect on TWIST1 expression was exclusively detectable in the patient's iPSC-derived neuronal cells, stressing the need for studying developmental disorders in the biologically relevant context. Chromosome conformation capture analyses show that TWIST1 lost genomic interactions with several enhancers due to the chromothripsis event, which likely led to deregulation of TWIST1 expression and contributed to the patient's craniosynostosis phenotype. CONCLUSIONS: We demonstrate that a combination of patient-derived iPSC differentiation and trio-based molecular profiling is a powerful approach to improve the interpretation of pathogenic complex genomic rearrangements. Here we have applied this approach to identify misexpression of TWIST1, FOXP1, and DPYD as key contributors to the complex congenital phenotype resulting from germline chromothripsis rearrangements.

The Genomic Scrapheap Challenge; Extracting Relevant Data from Unmapped Whole Genome Sequencing Reads, Including Strain Specific Genomic Segments, in Rats.

Unmapped next-generation sequencing reads are typically ignored while they contain biologically relevant information. We systematically analyzed unmapped reads from whole genome sequencing of 33 inbred rat strains. High quality reads were selected and enriched for biologically relevant sequences; similarity-based analysis revealed clustering similar to previously reported phylogenetic trees. Our results demonstrate that on average 20% of all unmapped reads harbor sequences that can be used to improve reference genomes and generate hypotheses on potential genotype-phenotype relationships. Analysis pipelines would benefit from incorporating the described methods and reference genomes would benefit from inclusion of the genomic segments obtained through these efforts.

Genomic landscape of rat strain and substrain variation.

BACKGROUND: Since the completion of the rat reference genome in 2003, whole-genome sequencing data from more than 40 rat strains have become available. These data represent the broad range of strains that are used in rat research including commonly used substrains. Currently, this wealth of information cannot be used to its full extent, because the variety of different variant calling algorithms employed by different groups impairs comparison between strains. In addition, all rat whole genome sequencing studies to date used an outdated reference genome for analysis (RGSC3.4 released in 2004). RESULTS: Here we present a comprehensive, multi-sample and uniformly called set of genetic variants in 40 rat strains, including 19 substrains. We reanalyzed all primary data using a recent version of the rat reference assembly (RGSC5.0 released in 2012) and identified over 12 million genomic variants (SNVs, indels and structural variants) among the 40 strains. 28,318 SNVs are specific to individual substrains, which may be explained by introgression from other unsequenced strains and ongoing evolution by genetic drift. Substrain SNVs may have a larger predicted functional impact compared to older shared SNVs. CONCLUSIONS: In summary we present a comprehensive catalog of uniformly analyzed genetic variants among 40 widely used rat inbred strains based on the RGSC5.0 assembly. This represents a valuable resource, which will facilitate rat functional genomic research. In line with previous observations, our genome-wide analyses do not show evidence for contribution of multiple ancestral founder rat subspecies to the currently used rat inbred strains, as is the case for mouse. In addition, we find that the degree of substrain variation is highly variable between strains, which is of importance for the correct interpretation of experimental data from different labs.

Lack of major genome instability in tumors of p53 null rats.

Tumorigenesis is often associated with loss of tumor suppressor genes (such as TP53), genomic instability and telomere lengthening. Previously, we generated and characterized a rat p53 knockout model in which the homozygous rats predominantly develop hemangiosarcomas whereas the heterozygous rats mainly develop osteosarcomas. Using genome-wide analyses, we find that the tumors that arise in the heterozygous and homozygous Tp53C273X mutant animals are also different in their genomic instability profiles. While p53 was fully inactivated in both heterozygous and homozygous knockout rats, tumors from homozygous animals show very limited aneuploidy and low degrees of somatic copy number variation as compared to the tumors from heterozygous animals. In addition, complex structural rearrangements such as chromothripsis and breakage-fusion-bridge cycles were never found in tumors from homozygous animals, while these were readily detectable in tumors from heterozygous animals. Finally, we measured telomere length and telomere lengthening pathway activity and found that tumors of homozygous animals have longer telomeres but do not show clear telomerase or alternative lengthening of telomeres (ALT) activity differences as compared to the tumors from heterozygous animals. Taken together, our results demonstrate that host p53 status in this rat p53 knockout model has a large effect on both tumor type and genomic instability characteristics, where full loss of functional p53 is not the main driver of large-scale structural variations. Our results also suggest that chromothripsis primarily occurs under p53 heterozygous rather than p53 null conditions.

Genomic and functional overlap between somatic and germline chromosomal rearrangements.

Genomic rearrangements are a common cause of human congenital abnormalities. However, their origin and consequences are poorly understood. We performed molecular analysis of two patients with congenital disease who carried de novo genomic rearrangements. We found that the rearrangements in both patients hit genes that are recurrently rearranged in cancer (ETV1, FOXP1, and microRNA cluster C19MC) and drive formation of fusion genes similar to those described in cancer. Subsequent analysis of a large set of 552 de novo germline genomic rearrangements underlying congenital disorders revealed enrichment for genes rearranged in cancer and overlap with somatic cancer breakpoints. Breakpoints of common (inherited) germline structural variations also overlap with cancer breakpoints but are depleted for cancer genes. We propose that the same genomic positions are prone to genomic rearrangements in germline and soma but that timing and context of breakage determines whether developmental defects or cancer are promoted.

Extensive localization of long noncoding RNAs to the cytosol and mono- and polyribosomal complexes.

BACKGROUND: Long noncoding RNAs (lncRNAs) form an abundant class of transcripts, but the function of the majority of them remains elusive. While it has been shown that some lncRNAs are bound by ribosomes, it has also been convincingly demonstrated that these transcripts do not code for proteins. To obtain a comprehensive understanding of the extent to which lncRNAs bind ribosomes, we performed systematic RNA sequencing on ribosome-associated RNA pools obtained through ribosomal fractionation and compared the RNA content with nuclear and (non-ribosome bound) cytosolic RNA pools. RESULTS: The RNA composition of the subcellular fractions differs significantly from each other, but lncRNAs are found in all locations. A subset of specific lncRNAs is enriched in the nucleus but surprisingly the majority is enriched in the cytosol and in ribosomal fractions. The ribosomal enriched lncRNAs include H19 and TUG1. CONCLUSIONS: Most studies on lncRNAs have focused on the regulatory function of these transcripts in the nucleus. We demonstrate that only a minority of all lncRNAs are nuclear enriched. Our findings suggest that many lncRNAs may have a function in cytoplasmic processes, and in particular in ribosome complexes.

Allelic exclusion of the immunoglobulin heavy chain locus is independent of its nuclear localization in mature B cells.

In developing B cells, the immunoglobulin heavy chain (IgH) locus is thought to move from repressive to permissive chromatin compartments to facilitate its scheduled rearrangement. In mature B cells, maintenance of allelic exclusion has been proposed to involve recruitment of the non-productive IgH allele to pericentromeric heterochromatin. Here, we used an allele-specific chromosome conformation capture combined with sequencing (4C-seq) approach to unambigously follow the individual IgH alleles in mature B lymphocytes. Despite their physical and functional difference, productive and non-productive IgH alleles in B cells and unrearranged IgH alleles in T cells share many chromosomal contacts and largely reside in active chromatin. In brain, however, the locus resides in a different repressive environment. We conclude that IgH adopts a lymphoid-specific nuclear location that is, however, unrelated to maintenance of allelic exclusion. We additionally find that in mature B cells-but not in T cells-the distal VH regions of both IgH alleles position themselves away from active chromatin. This, we speculate, may help to restrict enhancer activity to the productively rearranged VH promoter element.

Nucleosomal DNA binding drives the recognition of H3K36-methylated nucleosomes by the PSIP1-PWWP domain.

BACKGROUND: Recognition of histone modifications by specialized protein domains is a key step in the regulation of DNA-mediated processes like gene transcription. The structural basis of these interactions is usually studied using histone peptide models, neglecting the nucleosomal context. Here, we provide the structural and thermodynamic basis for the recognition of H3K36-methylated (H3K36me) nucleosomes by the PSIP1-PWWP domain, based on extensive mutational analysis, advanced nuclear magnetic resonance (NMR), and computational approaches. RESULTS: The PSIP1-PWWP domain binds H3K36me3 peptide and DNA with low affinity, through distinct, adjacent binding surfaces. PWWP binding to H3K36me nucleosomes is enhanced approximately 10,000-fold compared to a methylated peptide. Based on mutational analyses and NMR data, we derive a structure of the complex showing that the PWWP domain is bound to H3K36me nucleosomes through simultaneous interactions with both methylated histone tail and nucleosomal DNA. CONCLUSION: Concerted binding to the methylated histone tail and nucleosomal DNA underlies the high- affinity, specific recognition of H3K36me nucleosomes by the PSIP1-PWWP domain. We propose that this bipartite binding mechanism is a distinctive and general property in the recognition of histone modifications close to the nucleosome core.

Genetic basis of transcriptome differences between the founder strains of the rat HXB/BXH recombinant inbred panel.

BACKGROUND: With the advent of next generation sequencing it has become possible to detect genomic variation on a large scale. However, predicting which genomic variants are damaging to gene function remains a challenge, as knowledge of the effects of genomic variation on gene expression is still limited. Recombinant inbred panels are powerful tools to study the cis and trans effects of genetic variation on molecular phenotypes such as gene expression. RESULTS: We generated a comprehensive inventory of genomic differences between the two founder strains of the rat HXB/BXH recombinant inbred panel: SHR/OlaIpcv and BN-Lx/Cub. We identified 3.2 million single nucleotide variants, 425,924 small insertions and deletions, 907 copy number changes and 1,094 large structural genetic variants. RNA-sequencing analyses on liver tissue of the two strains identified 532 differentially expressed genes and 40 alterations in transcript structure. We identified both coding and non-coding variants that correlate with differential expression and alternative splicing. Furthermore, structural variants, in particular gene duplications, show a strong correlation with transcriptome alterations. CONCLUSIONS: We show that the panel is a good model for assessing the genetic basis of phenotypic heterogeneity and for providing insights into possible underlying molecular mechanisms. Our results reveal a high diversity and complexity underlying quantitative and qualitative transcriptional differences.

Variegated gene expression caused by cell-specific long-range DNA interactions.

Mammalian genomes contain numerous regulatory DNA sites with unknown target genes. We used mice with an extra ß-globin locus control region (LCR) to investigate how a regulator searches the genome for target genes. We find that the LCR samples a restricted nuclear subvolume, wherein it preferentially contacts genes controlled by shared transcription factors. No contacted gene is detectably upregulated except for endogenous ß-globin genes located on another chromosome. This demonstrates genetically that mammalian trans activation is possible, but suggests that it will be rare. Trans activation occurs not pan-cellularly, but in 'jackpot' cells enriched for the interchromosomal interaction. Therefore, cell-specific long-range DNA contacts can cause variegated expression.

Interactions among Polycomb domains are guided by chromosome architecture.

Polycomb group (PcG) proteins bind and regulate hundreds of genes. Previous evidence has suggested that long-range chromatin interactions may contribute to the regulation of PcG target genes. Here, we adapted the Chromosome Conformation Capture on Chip (4C) assay to systematically map chromosomal interactions in Drosophila melanogaster larval brain tissue. Our results demonstrate that PcG target genes interact extensively with each other in nuclear space. These interactions are highly specific for PcG target genes, because non-target genes with either low or high expression show distinct interactions. Notably, interactions are mostly limited to genes on the same chromosome arm, and we demonstrate that a topological rather than a sequence-based mechanism is responsible for this constraint. Our results demonstrate that many interactions among PcG target genes exist and that these interactions are guided by overall chromosome architecture.

High-resolution identification of balanced and complex chromosomal rearrangements by 4C technology.

Balanced chromosomal rearrangements can cause disease, but techniques for their rapid and accurate identification are missing. Here we demonstrate that chromatin conformation capture on chip (4C) technology can be used to screen large genomic regions for balanced and complex inversions and translocations at high resolution. The 4C technique can be used to detect breakpoints also in repetitive DNA sequences as it uniquely relies on capturing genomic fragments across the breakpoint. Using 4C, we uncovered LMO3 as a potentially leukemogenic translocation partner of TRB@. We developed multiplex 4C to simultaneously screen for translocation partners of multiple selected loci. We identified unsuspected translocations and complex rearrangements. Furthermore, using 4C we detected translocations even in small subpopulations of cells. This strategy opens avenues for the rapid fine-mapping of cytogenetically identified translocations and inversions, and the efficient screening for balanced rearrangements near candidate loci, even when rearrangements exist only in subpopulations of cells.

FISH-eyed and genome-wide views on the spatial organisation of gene expression.

Eukaryotic cells store their genome inside a nucleus, a dedicated organelle shielded by a double lipid membrane. Pores in these membranes allow the exchange of molecules between the nucleus and cytoplasm. Inside the mammalian cell nucleus, roughly 2 m of DNA, divided over several tens of chromosomes is packed. In addition, protein and RNA molecules functioning in DNA-metabolic processes such as transcription, replication, repair and the processing of RNA fill the nuclear space. While many of the nuclear proteins freely diffuse and display a more or less homogeneous distribution across the nuclear interior, some appear to preferentially cluster and form foci or bodies. A non-random structure is also observed for DNA: increasing evidence shows that selected parts of the genome preferentially contact each other, sometimes even at specific sites in the nucleus. Currently a lot of research is dedicated to understanding the functional significance of nuclear architecture, in particular with respect to the regulation of gene expression. Here we will evaluate evidence implying that the folding of DNA is important for transcriptional control in mammals and we will discuss novel high-throughput techniques expected to further boost our knowledge on nuclear organisation.

Maintenance of long-range DNA interactions after inhibition of ongoing RNA polymerase II transcription.

A relationship exists between nuclear architecture and gene activity and it has been proposed that the activity of ongoing RNA polymerase II transcription determines genome organization in the mammalian cell nucleus. Recently developed 3C and 4C technology allowed us to test the importance of transcription for nuclear architecture. We demonstrate that upon transcription inhibition binding of RNA polymerase II to gene regulatory elements is severely reduced. However, contacts between regulatory DNA elements and genes in the beta-globin locus are unaffected and the locus still interacts with the same genomic regions elsewhere on the chromosome. This is a general phenomenon since the great majority of intra- and interchromosomal interactions with the ubiquitously expressed Rad23a gene are also not affected. Our data demonstrate that without transcription the organization and modification of nucleosomes at active loci and the local binding of specific trans-acting factors is unaltered. We propose that these parameters, more than transcription or RNA polymerase II binding, determine the maintenance of long-range DNA interactions.

Three-dimensional organization of gene expression in erythroid cells.

The history of globin research is marked by a series of contributions seminal to our understanding of the genome, its function, and its relation to disease. For example, based on studies on hemoglobinopathies, it was understood that gene expression can be under the control of DNA elements that locate away from the genes on the linear chromosome template. Recent technological developments have allowed the demonstration that these regulatory DNA elements communicate with the genes through physical interaction, which loops out the intervening chromatin fiber. Subsequent studies showed that the spatial organization of the beta-globin locus dynamically changes in relation to differences in gene expression. Moreover, it was shown that the beta-globin locus adopts a different position in the nucleus during development and erythroid maturation. Here, we discuss the most recent insight into the three-dimensional organization of gene expression.

An evaluation of 3C-based methods to capture DNA interactions.

The shape of the genome is thought to play an important part in the coordination of transcription and other DNA-metabolic processes. Chromosome conformation capture (3C) technology allows us to analyze the folding of chromatin in the native cellular state at a resolution beyond that provided by current microscopy techniques. It has been used, for example, to demonstrate that regulatory DNA elements communicate with distant target genes through direct physical interactions that loop out the intervening chromatin fiber. Here we discuss the intricacies of 3C and new 3C-based methods including the 4C, 5C and ChIP-loop assay.