Integration of Hi-C with short and long-read genome sequencing reveals the structure of germline rearranged genomes

Here the authors characterize structural variations (SVs) in a cohort of individuals with complex genomic rearrangements, identifying breakpoints by employing short- and long-read genome sequencing and investigate their impact on gene expression and the three-dimensional chromatin architecture. They find breakpoints are enriched in inactive regions and can result in chromatin domain fusions.Structural variants are a common cause of disease and contribute to a large extent to inter-individual variability, but their detection and interpretation remain a challenge. Here, we investigate 11 individuals with complex genomic rearrangements including germline chromothripsis by combining short- and long-read genome sequencing (GS) with Hi-C. Large-scale genomic rearrangements are identified in Hi-C interaction maps, allowing for an independent assessment of breakpoint calls derived from the GS methods, resulting in >300 genomic junctions. Based on a comprehensive breakpoint detection and Hi-C, we achieve a reconstruction of whole rearranged chromosomes. Integrating information on the three-dimensional organization of chromatin, we observe that breakpoints occur more frequently than expected in lamina-associated domains (LADs) and that a majority reshuffle topologically associating domains (TADs). By applying phased RNA-seq, we observe an enrichment of genes showing allelic imbalanced expression (AIG) within 100 kb around the breakpoints. Interestingly, the AIGs hit by a breakpoint (19/22) display both up- and downregulation, thereby suggesting different mechanisms at play, such as gene disruption and rearrangements of regulatory information. However, the majority of interpretable genes located 200 kb around a breakpoint do not show significant expression changes. Thus, there is an overall robustness in the genome towards large-scale chromosome rearrangements

De novo unbalanced translocations have a complex history/aetiology

We investigated 52 cases of de novo unbalanced translocations, consisting in a terminally deleted or inverted-duplicated deleted (inv-dup del) 46th chromosome to which the distal portion of another chromosome or its opposite end was transposed. Array CGH, whole-genome sequencing, qPCR, FISH, and trio genotyping were applied. A biparental origin of the deletion and duplication was detected in 6 cases, whereas in 46, both imbalances have the same parental origin. Moreover, the duplicated region was of maternal origin in more than half of the cases, with 25% of them showing two maternal and one paternal haplotype. In all these cases, maternal age was increased. These findings indicate that the primary driver for the occurrence of the de novo unbalanced translocations is a maternal meiotic non-disjunction, followed by partial trisomy rescue of the supernumerary chromosome present in the trisomic zygote. In contrast, asymmetric breakage of a dicentric chromosome, originated either at the meiosis or postzygotically, in which the two resulting chromosomes, one being deleted and the other one inv-dup del, are repaired by telomere capture, appears at the basis of all inv-dup del translocations. Notably, this mechanism also fits with the origin of some simple translocations in which the duplicated region was of paternal origin. In all cases, the signature at the translocation junctions was that of non-homologous end joining (NHEJ) rather than non-allelic homologous recombination (NAHR). Our data imply that there is no risk of recurrence in the following pregnancies for any of the de novo unbalanced translocations we discuss here

Integration of Hi-C with short and long-read genome sequencing reveals the structure of germline rearranged genomes

Structural variants are a common cause of disease and contribute to a large extent to inter-individual variability, but their detection and interpretation remain a challenge. Here, we investigate 11 individuals with complex genomic rearrangements including germline chromothripsis by combining short- and long-read genome sequencing (GS) with Hi-C. Large-scale genomic rearrangements are identified in Hi-C interaction maps, allowing for an independent assessment of breakpoint calls derived from the GS methods, resulting in >300 genomic junctions. Based on a comprehensive breakpoint detection and Hi-C, we achieve a reconstruction of whole rearranged chromosomes. Integrating information on the three-dimensional organization of chromatin, we observe that breakpoints occur more frequently than expected in lamina-associated domains (LADs) and that a majority reshuffle topologically associating domains (TADs). By applying phased RNA-seq, we observe an enrichment of genes showing allelic imbalanced expression (AIG) within 100 kb around the breakpoints. Interestingly, the AIGs hit by a breakpoint (19/22) display both up- and downregulation, thereby suggesting different mechanisms at play, such as gene disruption and rearrangements of regulatory information. However, the majority of interpretable genes located 200 kb around a breakpoint do not show significant expression changes. Thus, there is an overall robustness in the genome towards large-scale chromosome rearrangements

Insertional translocation involving an additional nonchromothriptic chromosome in constitutional chromothripsis: Rule or exception?

Background: Chromothripsis, which is the local massive shattering of one or more chromosomes and their reassembly in a disordered array with frequent loss of some fragments, has been mainly reported in association with abnormal phenotypes. We report three unrelated healthy persons, two of which parenting a child with some degree of intellectual disability, carrying a chromothripsis involving respectively one, two, and three chromosomes, which was detected only after whole-genome sequencing. Unexpectedly, in all three cases a fragment from one of the chromothripsed chromosomes resulted to be inserted within a nonchromothripsed one. Methods: Conventional cytogenetic techniques, paired-end whole-genome sequencing, polymerase chain reaction, and Sanger sequencing were used to characterize complex rearrangements, copy-number variations, and breakpoint sequences in all three families. Results: In two families, one parent was carrier of a balanced chromothripsis causing in the index case a deletion and a noncontiguous duplication at 3q in case 1, and a t(6;14) translocation associated with interstitial 14q deletion in case 2. In the third family, an unbalanced chromothripsis involving chromosomes 6, 7, and 15 was inherited to the proband by the mosaic parent. In all three parents, the chromothripsis was concurrent with an insertional translocation of a portion of one of the chromothriptic chromosomes within a further chromosome that was not involved in the chromothripsis event. Conclusion: Our findings show that (a) both simple and complex unbalanced rearrangements may result by the recombination of a cryptic parental balanced chromothripsis and that (b) insertional translocations are the spy of more complex rearrangements and not simply a three-breakpoint event

Chromothripsis and ring chromosome 22: a paradigm of genomic complexity in the Phelan-McDermid syndrome (22q13 deletion syndrome)

INTRODUCTION: Phelan-McDermid syndrome (PMS) is caused by SHANK3 haploinsufficiency. Its wide phenotypic variation is attributed partly to the type and size of 22q13 genomic lesion (deletion, unbalanced translocation, ring chromosome), partly to additional undefined factors. We investigated a child with severe global neurodevelopmental delay (NDD) compatible with her distal 22q13 deletion, complicated by bilateral perisylvian polymicrogyria (BPP) and urticarial rashes, unreported in PMS. METHODS: Following the cytogenetic and array-comparative genomic hybridization (CGH) detection of a r(22) with SHANK3 deletion and two upstream duplications, whole-genome sequencing (WGS) in blood and whole-exome sequencing (WES) in blood and saliva were performed to highlight potential chromothripsis/chromoanagenesis events and any possible BPP-associated variants, even in low-level mosaicism. RESULTS: WGS confirmed the deletion and highlighted inversion and displaced order of eight fragments, three of them duplicated. The microhomology-mediated insertion of partial Alu-elements at one breakpoint junction disrupted the topological associating domain joining NFAM1 to the transcriptional coregulator TCF20. WES failed to detect BPP-associated variants. CONCLUSIONS: Although we were unable to highlight the molecular basis of BPP, our data suggest that SHANK3 haploinsufficiency and TCF20 misregulation, both associated with intellectual disability, contributed to the patient's NDD, while NFAM1 interruption likely caused her skin rashes, as previously reported. We provide the first example of chromoanasynthesis in a constitutional ring chromosome and reinforce the growing evidence that chromosomal rearrangements may be more complex than estimated by conventional diagnostic approaches and affect the phenotype by global alteration of the topological chromatin organisation rather than simply by deletion or duplication of dosage-sensitive genes

MBD5 molecular screening on Smith Magenis-like Syndrome patients without the typical 17p11.2 deletion

Smith Magenis Syndrome (SMS) is a complex heterogeneous disorder, caused by RAI1 haploinsufficiency and triggered by 17p11.2 deletion or RAI1 gene mutation. Only 50% of patients with a clinical suspicion of SMS bear the known defects, making its diagnosis challenging and suggesting that other loci may underlie SMS-like phenotypes. Notably 2q23.1 deletion and brachydactyly mental retardation (BDMR) syndromes, caused by MBD5 and HDAC4 haploinsufficiency respectively, are two SMS phenotypical and molecular overlapping syndromes as the affected patients seem to show a reduced RAI1 transcript levels, implying that MBD5/HDAC4 regulate RAI1. Here we describe a cohort of 30 SMS-like patients without 17p11.2 deletion, and RAI1 microdeletion or mutation. Array-CGH analysis revealed a 29 kb deletion at 2q23.1 encompassing exons 1 and 2 of MBD5 in one patient, while MBD5 mutational screening identified in a second case an heterozygous missense mutation likely pathogenetic, p.A857T. Unexpectedly the deletion was found in a mosaic condition in the healthy father of the first patient, while the missense mutation, previously described in a child with mental retardation, was maternally inherited. The subsequent RTqPCR showed a MBD5 transcript downregulation in the patient with 2q23.1 deletion, normal RAI1 transcript levels in all the patients and unexpectedly in the SMS-like patient with MBD5 haploinsufficiency too. Our findings corroborate the hypothesis that other loci can be causative of SMS-like phenotypes and for this reason we are setting up a flow chart aiming at improve the molecular diagnosis in SMS-like patients and shed light on pathogenetic mechanisms

Small supernumerary marker chromosomes: A legacy of trisomy rescue?

We studied by a whole genomic approach and trios genotyping, 12 de novo, non-recurrent small supernumerary marker chromosomes (sSMC), detected as mosaics during pre- or postnatal diagnosis and associated with increased maternal age. Four sSMCs contained pericentromeric portions only, whereas eight had additional non-contiguous portions of the same chromosome, assembled together in a disordered fashion by repair-based mechanisms in a chromothriptic event. Maternal hetero/isodisomy was detected with a paternal origin of the sSMC in some cases, whereas in others two maternal alleles in the sSMC region and biparental haplotypes of the homologs were detected. In other cases the homologs were biparental while the sSMC had the same haplotype of the maternally inherited chromosome. These findings strongly suggest that most sSMCs are the result of a multiple-step mechanism, initiated by maternal meiotic non-disjunction followed by post-zygotic anaphase lagging of the supernumerary chromosome and its subsequent chromothripsis

A Genetics-First Approach Revealed Monogenic Disorders in Patients With ARM and VACTERL Anomalies

Background: The VATER/VACTERL association (VACTERL) is defined as the non-random occurrence of the following congenital anomalies: Vertebral, Anal, Cardiac, Tracheal-Esophageal, Renal, and Limb anomalies. As no unequivocal candidate gene has been identified yet, patients are diagnosed phenotypically. The aims of this study were to identify patients with monogenic disorders using a genetics-first approach, and to study whether variants in candidate genes are involved in the etiology of VACTERL or the individual features of VACTERL: Anorectal malformation (ARM) or esophageal atresia with or without trachea-esophageal fistula (EA/TEF). Methods: Using molecular inversion probes, a candidate gene panel of 56 genes was sequenced in three patient groups: VACTERL (n = 211), ARM (n = 204), and EA/TEF (n = 95). Loss-of-function (LoF) and additional likely pathogenic missense variants, were prioritized and validated using Sanger sequencing. Validated variants were tested for segregation and patients were clinically re-evaluated. Results: In 7 out of the 510 patients (1.4%), pathogenic or likely pathogenic variants were identified in SALL1, SALL4, and MID1, genes that are associated with Townes-Brocks, Duane-radial-ray, and Opitz-G/BBB syndrome. These syndromes always include ARM or EA/TEF, in combination with at least two other VACTERL features. We did not identify LoF variants in the remaining candidate genes. Conclusions: None of the other candidate genes were identified as novel unequivocal disease genes for VACTERL. However, a genetics-first approach allowed refinement of the clinical diagnosis in seven patients, in whom an alternative molecular-based diagnosis was found with important implications for the counseling of the families