Cryptic breakpoint identified by whole-genome mate-pair sequencing in a rare paternally inherited complex chromosomal rearrangement

Abstract Background Precise characterization of apparently balanced complex chromosomal rearrangements in non-affected individuals is crucial as they may result in reproductive failure, recurrent miscarriages or affected offspring. Case presentation We present a family, where the non-affected father and daughter were found, using FISH and karyotyping, to be carriers of a three-way complex chromosomal rearrangement [t(6;7;10)(q16.2;q34;q26.1), de novo in the father]. The family suffered from two stillbirths, one miscarriage, and has a son with severe intellectual disability. In the present study, the family was revisited using whole-genome mate-pair sequencing. Interestingly, whole-genome mate-pair sequencing revealed a cryptic breakpoint on derivative (der) chromosome 6 rendering the rearrangement even more complex. FISH using a chromosome (chr) 6 custom-designed probe and a chr10 control probe confirmed that the interstitial chr6 segment, created by the two chr6 breakpoints, was translocated onto der(10). Breakpoints were successfully validated with Sanger sequencing, and small imbalances as well as microhomology were identified. Finally, the complex chromosomal rearrangement breakpoints disrupted the SIM1, GRIK2, CNTNAP2, and PTPRE genes without causing any phenotype development. Conclusions In contrast to the majority of maternally transmitted complex chromosomal rearrangement cases, our study investigated a rare case where a complex chromosomal rearrangement, which most probably resulted from a Type IV hexavalent during the pachytene stage of meiosis I, was stably transmitted from a fertile father to his non-affected daughter. Whole-genome mate-pair sequencing proved highly successful in identifying cryptic complexity, which consequently provided further insight into the meiotic segregation of chromosomes and the increased reproductive risk in individuals carrying the specific complex chromosomal rearrangement. We propose that such complex rearrangements should be characterized in detail using a combination of conventional cytogenetic and NGS-based approaches to aid in better prenatal preimplantation genetic diagnosis and counseling in couples with reproductive problems

Enrichment of megabase-sized DNA molecules for single-molecule optical mapping and next-generation sequencing

Abstract Next-generation sequencing (NGS) has caused a revolution, yet left a gap: long-range genetic information from native, non-amplified DNA fragments is unavailable. It might be obtained by optical mapping of megabase-sized DNA molecules. Frequently only a specific genomic region is of interest, so here we introduce a method for selection and enrichment of megabase-sized DNA molecules intended for single-molecule optical mapping: DNA from a human cell line is digested by the NotI rare-cutting enzyme and size-selected by pulsed-field gel electrophoresis. For demonstration, more than 600 sub-megabase- to megabase-sized DNA molecules were recovered from the gel and analysed by denaturation-renaturation optical mapping. Size-selected molecules from the same gel were sequenced by NGS. The optically mapped molecules and the NGS reads showed enrichment from regions defined by NotI restriction sites. We demonstrate that the unannotated genome can be characterized in a locus-specific manner via molecules partially overlapping with the annotated genome. The method is a promising tool for investigation of structural variants in enriched human genomic regions for both research and diagnostic purposes. Our enrichment method could potentially work with other genomes or target specified regions by applying other genomic editing tools, such as the CRISPR/Cas9 system

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

Balanced chromosomal rearrangements offer insights into coding and noncoding genomic features associated with developmental disorders

Balanced chromosomal rearrangements (BCRs), including inversions, translocations, and insertions, reorganize large sections of the genome and contribute substantial risk for developmental disorders (DDs). However, the rarity and lack of systematic screening for BCRs in the population has precluded unbiased analyses of the genomic features and mechanisms associated with risk for DDs versus normal developmental outcomes. Here, we sequenced and analyzed 1,420 BCR breakpoints across 710 individuals, including 406 DD cases and the first large-scale collection of 304 control BCR carriers. We found that BCRs were not more likely to disrupt genes in DD cases than controls, but were seven-fold more likely to disrupt genes associated with dominant DDs (21.3% of cases vs. 3.4% of controls; P = 1.60×10$^{−12}$). Moreover, BCRs that did not disrupt a known DD gene were significantly enriched for breakpoints that altered topologically associated domains (TADs) containing dominant DD genes in cases compared to controls (odds ratio [OR] = 1.43, P = 0.036). We discovered six TADs enriched for noncoding BCRs (false discovery rate < 0.1) that contained known DD genes (MEF2C, FOXG1, SOX9, BCL11A, BCL11B, and SATB2) and represent candidate pathogenic long-range positional effect (LRPE) loci. These six TADs were collectively disrupted in 7.4% of the DD cohort. Phased Hi-C analyses of five cases with noncoding BCR breakpoints localized to one of these putative LRPEs, the 5q14.3 TAD encompassing MEF2C, confirmed extensive disruption to local 3D chromatin structures and reduced frequency of contact between the MEF2C promoter and annotated enhancers. We further identified six genomic features enriched in TADs preferentially disrupted by noncoding BCRs in DD cases versus controls and used these features to build a model to predict TADs at risk for LRPEs across the genome. These results emphasize the potential impact of noncoding structural variants to cause LRPEs in unsolved DD cases, as well as the complex interaction of features associated with predicting three-dimensional chromatin structures intolerant to disruption

Prioritization of genes driving congenital phenotypes of patients with de novo genomic structural variants

Background:Genomic structural variants (SVs) can affect many genes and regulatory elements. Therefore, the molecular mechanisms driving the phenotypes of patients carrying de novo SVs are frequently unknown. Methods:We applied a combination of systematic experimental and bioinformatic methods to improve the molecular diagnosis of 39 patients with multiple congenital abnormalities and/or intellectual disability harboring apparent de novo SVs, most with an inconclusive diagnosis after regular genetic testing. Results: In 7 of these cases (18%), whole-genome sequencing analysis revealed disease-relevant complexities of the SVs missed in routine microarray-based analyses. We developed a computational tool to predict the effects on genes directly affected by SVs and on genes indirectly affected likely due to the changes in chromatin organization and impact on regulatory mechanisms. By combining these functional predictions with extensive phenotype information, candidate driver genes were identified in 16/39 (41%) patients. In 8 cases, evidence was found for the involvement of multiple candidate drivers contributing to different parts of the phenotypes. Subsequently, we applied this computational method to two cohorts containing a total of 379 patients with previously detected and classified de novo SVs and identified candidate driver genes in 189 cases (50%), including 40 cases whose SVs were previously not classified as pathogenic. Pathogenic position effects were predicted in 28% of all studied cases with balanced SVs and in 11% of the cases with copy number variants. Conclusions:These results demonstrate an integrated computational and experimental approach to predict driver genes based on analyses of WGS data with phenotype association and chromatin organization datasets. These analyses nominate new pathogenic loci and have strong potential to improve the molecular diagnosis of patients with de novo SVs

Rare Pathogenic Variants in Genes Implicated in Glutamatergic Neurotransmission Pathway Segregate with Schizophrenia in Pakistani Families

Schizophrenia is a disabling neuropsychiatric disorder of adulthood onset with high heritability. Worldwide collaborations have identified an association of ~270 common loci, with small individual effects and hence weak clinical implications. The recent technological feasibility of exome sequencing enables the identification of rare variants of high penetrance that refine previous findings and improve risk assessment and prognosis. We recruited two multiplex Pakistani families, having 11 patients and 19 unaffected individuals in three generations. We performed genome-wide SNP genotyping, next-generation mate pairing and whole-exome sequencing of selected members to unveil genetic components. Candidate variants were screened in unrelated cohorts of 508 cases, 300 controls and fifteen families (with 51 affected and 47 unaffected individuals) of Pakistani origin. The structural impact of substituted residues was assessed through in silico modeling using iTASSER. In one family, we identified a rare novel microduplication (5q14.1_q14.2) encompassing critical genes involved in glutamate signaling, such as CMYA5, HOMER and RasGRF2. The second family segregates two ultra-rare, predicted pathogenic variants in the GRIN2A (NM_001134407.3: c.3505C>T, (p.R1169W) and in the NRG3 NM_001010848.4: c.1951G>A, (p.E651K). These genes encode for parts of AMPA and NMDA receptors of glutamatergic neurotransmission, respectively, and the variants are predicted to compromise protein function by destabilizing their structures. The variants were absent in the aforementioned cohorts. Our findings suggest that rare, highly penetrant variants of genes involved in glutamatergic neurotransmission are contributing to the etiology of schizophrenia in these families. It also highlights that genetic investigations of multiplex, multigenerational families could be a powerful approach to identify rare genetic variants involved in complex disorders