A regulatory archipelago controls Hox genes transcription in digits.

The evolution of digits was an essential step in the success of tetrapods. Among the key players, Hoxd genes are coordinately regulated in developing digits, where they help organize growth and patterns. We identified the distal regulatory sites associated with these genes by probing the three-dimensional architecture of this regulatory unit in developing limbs. This approach, combined with in vivo deletions of distinct regulatory regions, revealed that the active part of the gene cluster contacts several enhancer-like sequences. These elements are dispersed throughout the nearby gene desert, and each contributes either quantitatively or qualitatively to Hox gene transcription in presumptive digits. We propose that this genetic system, which we call a "regulatory archipelago," provides an inherent flexibility that may partly underlie the diversity in number and morphology of digits across tetrapods, as well as their resilience to drastic variations.

CTCF Binding Polarity Determines Chromatin Looping.

CCCTC-binding factor (CTCF) is an architectural protein involved in the three-dimensional (3D) organization of chromatin. In this study, we assayed the 3D genomic contact profiles of a large number of CTCF binding sites with high-resolution 4C-seq. As recently reported, our data also suggest that chromatin loops preferentially form between CTCF binding sites oriented in a convergent manner. To directly test this, we used CRISPR/Cas9 genome editing to delete core CTCF binding sites in three loci, including the CTCF site in the Sox2 super-enhancer. In all instances, CTCF and cohesin recruitment were lost, and chromatin loops with distal, convergent CTCF sites were disrupted or destabilized. Re-insertion of oppositely oriented CTCF recognition sequences restored CTCF and cohesin recruitment, but did not re-establish chromatin loops. We conclude that CTCF binding polarity plays a functional role in the formation of higher-order chromatin structure.

Large-scale discovery of mouse transgenic integration sites reveals frequent structural variation and insertional mutagenesis.

Transgenesis has been a mainstay of mouse genetics for over 30 yr, providing numerous models of human disease and critical genetic tools in widespread use today. Generated through the random integration of DNA fragments into the host genome, transgenesis can lead to insertional mutagenesis if a coding gene or an essential element is disrupted, and there is evidence that larger scale structural variation can accompany the integration. The insertion sites of only a tiny fraction of the thousands of transgenic lines in existence have been discovered and reported, due in part to limitations in the discovery tools. Targeted locus amplification (TLA) provides a robust and efficient means to identify both the insertion site and content of transgenes through deep sequencing of genomic loci linked to specific known transgene cassettes. Here, we report the first large-scale analysis of transgene insertion sites from 40 highly used transgenic mouse lines. We show that the transgenes disrupt the coding sequence of endogenous genes in half of the lines, frequently involving large deletions and/or structural variations at the insertion site. Furthermore, we identify a number of unexpected sequences in some of the transgenes, including undocumented cassettes and contaminating DNA fragments. We demonstrate that these transgene insertions can have phenotypic consequences, which could confound certain experiments, emphasizing the need for careful attention to control strategies. Together, these data show that transgenic alleles display a high rate of potentially confounding genetic events and highlight the need for careful characterization of each line to assure interpretable and reproducible experiments.

Targeted locus amplification to detect molecular markers in mantle cell and follicular lymphoma.

Minimal residual disease (MRD) monitoring by PCR methods is a strong and standardized predictor of clinical outcome in mantle cell lymphoma (MCL) and follicular lymphoma (FL). However, about 20% of MCL and 40% of FL patients lack a reliable molecular marker, being thus not eligible for MRD studies. Recently, targeted locus amplification (TLA), a next-generation sequencing (NGS) method based on the physical proximity of DNA sequences for target selection, identified novel gene rearrangements in leukemia. The aim of this study was to test TLA in MCL and FL diagnostic samples lacking a classical, PCR-detectable, t(11; 14) MTC (BCL1/IGH), or t(14; 18) major breakpoint region and minor cluster region (BCL2/IGH) rearrangements. Overall, TLA was performed on 20 MCL bone marrow (BM) or peripheral blood (PB) primary samples and on 20 FL BM, identifying a novel BCL1 or BCL2/IGH breakpoint in 16 MCL and 8 FL patients (80% and 40%, respectively). These new breakpoints (named BCL1-TLA and BCL2-TLA) were validated by ASO primers design and compared as MRD markers to classical IGH rearrangements in eight MCL: overall, MRD results by BCL1-TLA were superimposable (R Pearson = 0.76) to the standardized IGH-based approach. Moreover, MRD by BCL2-TLA reached good sensitivity levels also in FL and was predictive of a primary refractory case. In conclusion, this study offers the proof of principle that TLA is a promising and reliable NGS-based technology for the identification of novel molecular markers, suitable for further MRD analysis in previously not traceable MCL and FL patients.

Sensitive Monogenic Noninvasive Prenatal Diagnosis by Targeted Haplotyping.

During pregnancy, cell-free DNA (cfDNA) in maternal blood encompasses a small percentage of cell-free fetal DNA (cffDNA), an easily accessible source for determination of fetal disease status in risk families through non-invasive procedures. In case of monogenic heritable disease, background maternal cfDNA prohibits direct observation of the maternally inherited allele. Non-invasive prenatal diagnostics (NIPD) of monogenic diseases therefore relies on parental haplotyping and statistical assessment of inherited alleles from cffDNA, techniques currently unavailable for routine clinical practice. Here, we present monogenic NIPD (MG-NIPD), which requires a blood sample from both parents, for targeted locus amplification (TLA)-based phasing of heterozygous variants selectively at a gene of interest. Capture probes-based targeted sequencing of cfDNA from the pregnant mother and a tailored statistical analysis enables predicting fetal gene inheritance. MG-NIPD was validated for 18 pregnancies, focusing on CFTR, CYP21A2, and HBB. In all cases we could predict the inherited alleles with >98% confidence, even at relatively early stages (8 weeks) of pregnancy. This prediction and the accuracy of parental haplotyping was confirmed by sequencing of fetal material obtained by parallel invasive procedures. MG-NIPD is a robust method that requires standard instrumentation and can be implemented in any clinic to provide families carrying a severe monogenic disease with a prenatal diagnostic test based on a simple blood draw.

The pluripotent genome in three dimensions is shaped around pluripotency factors.

It is becoming increasingly clear that the shape of the genome importantly influences transcription regulation. Pluripotent stem cells such as embryonic stem cells were recently shown to organize their chromosomes into topological domains that are largely invariant between cell types. Here we combine chromatin conformation capture technologies with chromatin factor binding data to demonstrate that inactive chromatin is unusually disorganized in pluripotent stem-cell nuclei. We show that gene promoters engage in contacts between topological domains in a largely tissue-independent manner, whereas enhancers have a more tissue-restricted interaction profile. Notably, genomic clusters of pluripotency factor binding sites find each other very efficiently, in a manner that is strictly pluripotent-stem-cell-specific, dependent on the presence of Oct4 and Nanog protein and inducible after artificial recruitment of Nanog to a selected chromosomal site. We conclude that pluripotent stem cells have a unique higher-order genome structure shaped by pluripotency factors. We speculate that this interactome enhances the robustness of the pluripotent state.

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.

The inactive X chromosome adopts a unique three-dimensional conformation that is dependent on Xist RNA.

Three-dimensional topology of DNA in the cell nucleus provides a level of transcription regulation beyond the sequence of the linear DNA. To study the relationship between the transcriptional activity and the spatial environment of a gene, we used allele-specific chromosome conformation capture-on-chip (4C) technology to produce high-resolution topology maps of the active and inactive X chromosomes in female cells. We found that loci on the active X form multiple long-range interactions, with spatial segregation of active and inactive chromatin. On the inactive X, silenced loci lack preferred interactions, suggesting a unique random organization inside the inactive territory. However, escapees, among which is Xist, are engaged in long-range contacts with each other, enabling identification of novel escapees. Deletion of Xist results in partial refolding of the inactive X into a conformation resembling the active X without affecting gene silencing or DNA methylation. Our data point to a role for Xist RNA in shaping the conformation of the inactive X chromosome at least partially independent of transcription.

Genetic Screening Test to Detect Translocations in Acute Leukemias by Use of Targeted Locus Amplification.

BACKGROUND: Over 500 translocations have been identified in acute leukemia. To detect them, most diagnostic laboratories use karyotyping, fluorescent in situ hybridization, and reverse transcription PCR. Targeted locus amplification (TLA), a technique using next-generation sequencing, now allows detection of the translocation partner of a specific gene, regardless of its chromosomal origin. We present a TLA multiplex assay as a potential first-tier screening test for detecting translocations in leukemia diagnostics. METHODS: The panel includes 17 genes involved in many translocations present in acute leukemias. Procedures were optimized by using a training set of cell line dilutions and 17 leukemia patient bone marrow samples and validated by using a test set of cell line dilutions and a further 19 patient bone marrow samples. Per gene, we determined if its region was involved in a translocation and, if so, the translocation partner. To balance sensitivity and specificity, we introduced a gray zone showing indeterminate translocation calls needing confirmation. We benchmarked our method against results from the 3 standard diagnostic tests. RESULTS: In patient samples passing QC, we achieved a concordance with benchmarking tests of 81% in the training set and 100% in the test set, after confirmation of 4 and nullification of 3 gray zone calls (in total). In cell line dilutions, we detected translocations in 10% aberrant cells at several genetic loci. CONCLUSIONS: Multiplex TLA shows promising results as an acute leukemia screening test. It can detect cryptic and other translocations in selected genes. Further optimization may make this assay suitable for diagnostic use.

Quantitative analysis of chromatin interaction changes upon a 4.3 Mb deletion at mouse 4E2.

BACKGROUND: Circular chromosome conformation capture (4C) has provided important insights into three dimensional (3D) genome organization and its critical impact on the regulation of gene expression. We developed a new quantitative framework based on polymer physics for the analysis of paired-end sequencing 4C (PE-4Cseq) data. We applied this strategy to the study of chromatin interaction changes upon a 4.3 Mb DNA deletion in mouse region 4E2. RESULTS: A significant number of differentially interacting regions (DIRs) and chromatin compaction changes were detected in the deletion chromosome compared to a wild-type (WT) control. Selected DIRs were validated by 3D DNA FISH experiments, demonstrating the robustness of our pipeline. Interestingly, significant overlaps of DIRs with CTCF/Smc1 binding sites and differentially expressed genes were observed. CONCLUSIONS: Altogether, our PE-4Cseq analysis pipeline provides a comprehensive characterization of DNA deletion effects on chromatin structure and function.

The complex transcription regulatory landscape of our genome: control in three dimensions.

The non-coding part of our genome contains sequence motifs that can control gene transcription over distance. Here, we discuss functional genomics studies that uncover and characterize these sequences across the mammalian genome. The picture emerging is of a genome being a complex regulatory landscape. We explore the principles that underlie the wiring of regulatory DNA sequences and genes. We argue transcriptional control over distance can be understood when considering action in the context of the folded genome. Genome topology is expected to differ between individual cells, and this may cause variegated expression. High-resolution three-dimensional genome topology maps, ultimately of single cells, are required to understand the cis-regulatory networks that underlie cellular transcriptomes.

Robust 4C-seq data analysis to screen for regulatory DNA interactions.

Regulatory DNA elements can control the expression of distant genes via physical interactions. Here we present a cost-effective methodology and computational analysis pipeline for robust characterization of the physical organization around selected promoters and other functional elements using chromosome conformation capture combined with high-throughput sequencing (4C-seq). Our approach can be multiplexed and routinely integrated with other functional genomics assays to facilitate physical characterization of gene regulation.

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.

Nuclear organization of active and inactive chromatin domains uncovered by chromosome conformation capture-on-chip (4C).

The spatial organization of DNA in the cell nucleus is an emerging key contributor to genomic function. We developed 4C technology (chromosome conformation capture (3C)-on-chip), which allows for an unbiased genome-wide search for DNA loci that contact a given locus in the nuclear space. We demonstrate here that active and inactive genes are engaged in many long-range intrachromosomal interactions and can also form interchromosomal contacts. The active beta-globin locus in fetal liver preferentially contacts transcribed, but not necessarily tissue-specific, loci elsewhere on chromosome 7, whereas the inactive locus in fetal brain contacts different transcriptionally silent loci. A housekeeping gene in a gene-dense region on chromosome 8 forms long-range contacts predominantly with other active gene clusters, both in cis and in trans, and many of these intra- and interchromosomal interactions are conserved between the tissues analyzed. Our data demonstrate that chromosomes fold into areas of active chromatin and areas of inactive chromatin and establish 4C technology as a powerful tool to study nuclear architecture.

Diverse gene reprogramming events occur in the same spatial clusters of distal regulatory elements.

The spatial organization of genes in the interphase nucleus plays an important role in establishment and regulation of gene expression. Contradicting results have been reported to date, with little consensus about the dynamics of nuclear organization and the features of the contact loci. In this study, we investigated the properties and dynamics of genomic loci that are in contact with glucocorticoid receptor (GR)-responsive loci. We took a systematic approach, combining genome-wide interaction profiling by the chromosome conformation capture on chip (4C) technology with expression, protein occupancy, and chromatin accessibility profiles. This approach allowed a comprehensive analysis of how distinct features of the linear genome are organized in the three-dimensional nuclear space in the context of rapid gene regulation. We found that the transcriptional response to GR occurs without dramatic nuclear reorganization. Moreover, contrary to the view of transcription-driven organization, even genes with opposite transcriptional responses colocalize. Regions contacting GR-regulated genes are not particularly enriched for GR-regulated loci or for any functional group of genes, suggesting that these subnuclear environments are not organized to respond to a specific factor. The contact regions are, however, highly enriched for DNase I-hypersensitive sites that comprehensively mark cell-type-specific regulatory sites. These findings indicate that the nucleus is pre-organized in a conformation allowing rapid transcriptional reprogramming, and this organization is significantly correlated with cell-type-specific chromatin sites accessible to regulatory factors. Numerous open chromatin loci may be arranged in nuclear domains that are poised to respond to diverse signals in general and to permit efficient gene regulation.

Determining long-range chromatin interactions for selected genomic sites using 4C-seq technology: from fixation to computation.

Chromosome Conformation Capture (3C) and 3C-based technologies are constantly evolving in order to probe nuclear organization with higher depth and resolution. One such method is 4C-technology that allows the investigation of the nuclear environment of a locus of choice. The use of Illumina next generation sequencing as a detection platform for the analysis of 4C data has further improved the sensitivity and resolution of this method. Here we provide a step-by-step protocol for 4C-seq, describing the procedure from the initial template preparation until the final data analysis, interchanged with background information and considerations.

The beta-globin nuclear compartment in development and erythroid differentiation.

Efficient transcription of genes requires a high local concentration of the relevant trans-acting factors. Nuclear compartmentalization can provide an effective means to locally increase the concentration of rapidly moving trans-acting factors; this may be achieved by spatial clustering of chromatin-associated binding sites for such factors. Here we analyze the structure of an erythroid-specific spatial cluster of cis-regulatory elements and active beta-globin genes, the active chromatin hub (ACH; ref. 6), at different stages of development and in erythroid progenitors. We show, in mice and humans, that a core ACH is developmentally conserved and consists of the hypersensitive sites (HS1-HS6) of the locus control region (LCR), the upstream 5' HS-60/-62 and downstream 3' HS1. Globin genes switch their interaction with this cluster during development, correlating with the switch in their transcriptional activity. In mouse erythroid progenitors that are committed to but do not yet express beta-globin, only the interactions between 5' HS-60/-62, 3' HS1 and hypersensitive sites at the 5' side of the LCR are stably present. After induction of differentiation, these sites cluster with the rest of the LCR and the gene that is activated. We conclude that during erythroid differentiation, cis-regulatory DNA elements create a developmentally conserved nuclear compartment dedicated to RNA polymerase II transcription of beta-globin genes.

Targeted Locus Amplification and Haplotyping.

Targeted locus amplification (TLA) allows for the detection of all genetic variation (including structural variation) in a genomic region of interest. As TLA is based on proximity ligation, variants can be linked to each other, thereby enabling allelic phasing and the generation of haplotypes. This allows for the study of genetic variants in an allele-specific manner. Here, we provide a step-by-step protocol for TLA sample preparation and a complete bioinformatics pipeline for the allelic phasing of TLA data. Additionally, to illustrate the protocol, we show the ability of TLA to re-sequence and haplotype the complete cystic fibrosis transmembrane (CFTR) gene (> 200 kb in size) from patient-derived intestinal organoids.

Formalin-Fixed, Paraffin-Embedded-Targeted Locus Capture: A Next-Generation Sequencing Technology for Accurate DNA-Based Gene Fusion Detection in Bone and Soft Tissue Tumors.

Chromosomal rearrangements are important drivers in cancer, and their robust detection is essential for diagnosis, prognosis, and treatment selection, particularly for bone and soft tissue tumors. Current diagnostic methods are hindered by limitations, including difficulties with multiplexing targets and poor quality of RNA. A novel targeted DNA-based next-generation sequencing method, formalin-fixed, paraffin-embedded-targeted locus capture (FFPE-TLC), has shown advantages over current diagnostic methods when applied on FFPE lymphomas, including the ability to detect novel rearrangements. We evaluated the utility of FFPE-TLC in bone and soft tissue tumor diagnostics. FFPE-TLC sequencing was successfully applied on noncalcified and decalcified FFPE samples (n = 44) and control samples (n = 19). In total, 58 rearrangements were identified in 40 FFPE tumor samples, including three previously negative samples, and none was identified in the FFPE control samples. In all five discordant cases, FFPE-TLC could identify gene fusions where other methods had failed due to either detection limits or poor sample quality. FFPE-TLC achieved a high specificity and sensitivity (no false positives and negatives). These results indicate that FFPE-TLC is applicable in cancer diagnostics to simultaneously analyze many genes for their involvement in gene fusions. Similar to the observation in lymphomas, FFPE-TLC is a good DNA-based alternative to the conventional methods for detection of rearrangements in bone and soft tissue tumors.

Targeted locus amplification to develop robust patient-specific assays for liquid biopsies in pediatric solid tumors.

Background: Liquid biopsies combine minimally invasive sample collection with sensitive detection of residual disease. Pediatric malignancies harbor tumor-driving copy number alterations or fusion genes, rather than recurrent point mutations. These regions contain tumor-specific DNA breakpoint sequences. We investigated the feasibility to use these breakpoints to design patient-specific markers to detect tumor-derived cell-free DNA (cfDNA) in plasma from patients with pediatric solid tumors. Materials and methods: Regions of interest (ROI) were identified through standard clinical diagnostic pipelines, using SNP array for CNAs, and FISH or RT-qPCR for fusion genes. Using targeted locus amplification (TLA) on tumor organoids grown from tumor material or targeted locus capture (TLC) on FFPE material, ROI-specific primers and probes were designed, which were used to design droplet digital PCR (ddPCR) assays. cfDNA from patient plasma at diagnosis and during therapy was analyzed. Results: TLA was performed on material from 2 rhabdomyosarcoma, 1 Ewing sarcoma and 3 neuroblastoma. FFPE-TLC was performed on 8 neuroblastoma tumors. For all patients, at least one patient-specific ddPCR was successfully designed and in all diagnostic plasma samples the patient-specific markers were detected. In the rhabdomyosarcoma and Ewing sarcoma patients, all samples after start of therapy were negative. In neuroblastoma patients, presence of patient-specific markers in cfDNA tracked tumor burden, decreasing during induction therapy, disappearing at complete remission and re-appearing at relapse. Conclusion: We demonstrate the feasibility to determine tumor-specific breakpoints using TLA/TLC in different pediatric solid tumors and use these for analysis of cfDNA from plasma. Considering the high prevalence of CNAs and fusion genes in pediatric solid tumors, this approach holds great promise and deserves further study in a larger cohort with standardized plasma sampling protocols.