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

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

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

A germline chromothripsis event stably segregating in 11 individuals through three generations

Purpose: Parentally transmitted germ-line chromothripsis (G-CTH) has been identified in only a few cases. Most of these rearrangements were stably transmitted, in an unbalanced form, from a healthy mother to her child with congenital abnormalities probably caused by de novo copy-number changes of dosage sensitive genes. We describe a G-CTH transmitted through three generations in 11 healthy carriers. Methods: Conventional cytogenetic analysis, mate-pair sequencing, and polymerase chain reaction (PCR) were used to identify the chromosome rearrangement and characterize the breakpoints in all three generations. Results: We identified an apparently balanced translocation t(3;5), later shown to be a G-CTH, in all individuals of a three-generation family. The G-CTH stably segregated without occurrence of additional rearrangements; however, several spontaneous abortions were reported, possibly due to unbalanced transmission. Although seven protein-coding genes are interrupted, no clinical features can be definitively attributed to the affected genes. However, it can be speculated that truncation of one of these genes, encoding ataxia-telangiectasia and Rad3-related protein kinase (ATR), a key component of the DNA damage response, may be related to G-CTH formation. Conclusion: G-CTH rearrangements are not always associated with abnormal phenotypes and may be misinterpreted as balanced two-way translocations, suggesting that G-CTH is an underdiagnosed phenomenon