Whole-Genome Sequencing of Cytogenetically Balanced Chromosome Translocations Identifies Potentially Pathological Gene Disruptions and Highlights the Importance of Microhomology in the Mechanism of Formation

Genome sequencing differentiates a paracentric inversion from a balanced insertion enabling more accurate preimplantation genetic testing

INTRODUCTION

Distinguishing paracentric inversions (PAIs) from chromosomal insertions has traditionally relied on fluorescent in situ hybridization (FISH) techniques, but recent advancements in high-throughput sequencing have enabled the use of genome sequencing for such differentiation. In this study, we present a 38-year-old male carrier of a paracentric inversion on chromosome 2q, inv (2)(q31.2q34), whose partner experienced recurrent miscarriages.

MATERIAL AND METHODS

FISH analysis confirmed the inversion, and genome sequencing was employed for detailed characterization.

RESULTS

Preimplantation genetic testing (PGT) revealed that all assessed embryos were balanced, consistent with the low risk of unbalanced offspring associated with PAIs. While PAI carriers traditionally exhibit low risk of producing unbalanced offspring, exceptions exist due to crossover events within the inversion loop. Although the sample size was limited, the findings align with existing sperm study data, supporting the rare occurrence of unbalanced progeny in PAI carriers.

CONCLUSIONS

This study highlights the possibility of characterizing PAIs using genome sequencing to enable correct reproductive counseling and PGT decisions. Detailed characterization of a PAI is crucial for understanding potential outcomes and guiding PGT strategies, as accurate knowledge of the inversion size is essential for appropriate method selection in PGT. Given the very low risk of unbalanced offspring in PAI carriers, routine PGT may not be warranted but should be considered in specific cases with a history of unbalanced progeny or recurrent miscarriages. This study contributes to our understanding of PAI segregation and its implications for reproductive outcomes.

An Integrated Transcriptomics and Genomics Approach Detects an X/Autosome Translocation in a Female with Duchenne Muscular Dystrophy

Duchenne and Becker muscular dystrophies, caused by pathogenic variants in DMD, are the most common inherited neuromuscular conditions in childhood. These diseases follow an X-linked recessive inheritance pattern, and mainly males are affected. The most prevalent pathogenic variants in the DMD gene are copy number variants (CNVs), and most patients achieve their genetic diagnosis through Multiplex Ligation-dependent Probe Amplification (MLPA) or exome sequencing. Here, we investigated a female patient presenting with muscular dystrophy who remained genetically undiagnosed after MLPA and exome sequencing. RNA sequencing (RNAseq) from the patient's muscle biopsy identified an 85% reduction in DMD expression compared to 116 muscle samples included in the cohort. A de novo balanced translocation between chromosome 17 and the X chromosome (t(X;17)(p21.1;q23.2)) disrupting the DMD and BCAS3 genes was identified through trio whole genome sequencing (WGS). The combined analysis of RNAseq and WGS played a crucial role in the detection and characterisation of the disease-causing variant in this patient, who had been undiagnosed for over two decades. This case illustrates the diagnostic odyssey of female DMD patients with complex structural variants that are not detected by current panel or exome sequencing analysis.

Inverted triplications formed by iterative template switches generate structural variant diversity at genomic disorder loci

The duplication-triplication/inverted-duplication (DUP-TRP/INV-DUP) structure is a complex genomic rearrangement (CGR). Although it has been identified as an important pathogenic DNA mutation signature in genomic disorders and cancer genomes, its architecture remains unresolved. Here, we studied the genomic architecture of DUP-TRP/INV-DUP by investigating the DNA of 24 patients identified by array comparative genomic hybridization (aCGH) on whom we found evidence for the existence of 4 out of 4 predicted structural variant (SV) haplotypes. Using a combination of short-read genome sequencing (GS), long-read GS, optical genome mapping, and single-cell DNA template strand sequencing (strand-seq), the haplotype structure was resolved in 18 samples. The point of template switching in 4 samples was shown to be a segment of ∼2.2-5.5 kb of 100% nucleotide similarity within inverted repeat pairs. These data provide experimental evidence that inverted low-copy repeats act as recombinant substrates. This type of CGR can result in multiple conformers generating diverse SV haplotypes in susceptible dosage-sensitive loci.

Glass syndrome derived from chromosomal breakage downstream region of SATB2

BACKGROUND

Glass syndrome, derived from chromosomal 2q33.1 microdeletions, manifests with intellectual disability, microcephaly, epilepsy, and distinctive features, including micrognathia, down-slanting palpebral fissures, cleft palate, and crowded teeth. Recently, SATB2 located within the deletion region, was identified as the causative gene responsible for Glass syndrome. Numerous disease-causing variants within the SATB2 coding region have been reported.

OBJECTIVE

Given the presentation of intellectual disability and multiple congenital anomalies in a patient with a de novo reciprocal translocation between chromosomes 1 and 2, disruption of the causative gene(s) was suspected. This study sought to identify the causative gene in the patient.

METHODS

Long-read whole-genome sequencing was performed, and the expression level of the candidate gene was analyzed.

RESULTS

The detection of breakpoints was successful. While the breakpoint on chromosome 1 disrupted RNF220, it was not deemed to be a genetic cause. Conversely, SATB2 is located in the approximately 100-kb telomeric region of the breakpoint on chromosome 2. The patient's clinical features resembled those of previously reported cases of Glass syndrome, despite the lack of confirmed reduced SATB2 expression.

CONCLUSION

The patient was diagnosed with Glass syndrome due to the similarity in clinical features. This led us to hypothesize that disruption in the downstream region of SATB2 could result in Glass syndrome. The microhomologies identified in the breakpoint junctions indicate a potential molecular mechanism involving microhomology-mediated break-induced repair mechanism or template switching.

Mind the gap: the relevance of the genome reference to resolve rare and pathogenic inversions

Both long-read genome sequencing (lrGS) and the recently published Telomere to Telomere (T2T) reference genome provide increased coverage and resolution across repetitive regions promising heightened structural variant detection and improved mapping. Inversions (INV), intrachromosomal segments which are rotated 180° and inserted back into the same chromosome, are a class of structural variants particularly challenging to detect due to their copy-number neutral state and association with repetitive regions. Inversions represent about 1/20 of all balanced structural chromosome aberrations and can lead to disease by gene disruption or altering regulatory regions of dosage sensitive genes in cis . Here we remapped the genome data from six individuals carrying unsolved cytogenetically detected inversions. An INV6 and INV10 were resolved using GRCh38 and T2T-CHM13. Finally, an INV9 required optical genome mapping, de novo assembly of lrGS data and T2T-CHM13. This inversion disrupted intron 25 of EHMT1, confirming a diagnosis of Kleefstra syndrome 1 (MIM#610253). These three inversions, only mappable in specific references, prompted us to investigate the presence and population frequencies of differential reference regions (DRRs) between T2T-CHM13, GRCh37, GRCh38, the chimpanzee and bonobo, and hundreds of megabases of DRRs were identified. Our results emphasize the significance of the chosen reference genome and the added benefits of lrGS and optical genome mapping in solving rearrangements in challenging regions of the genome. This is particularly important for inversions and may impact clinical diagnostics.

Evaluation of genetic risk of apparently balanced chromosomal rearrangement carriers by breakpoint characterization

PURPOSE

To report genetic characteristics and associated risk of chromosomal breaks due to chromosomal rearrangements in large samples.

METHODS

MicroSeq, a technique that combines chromosome microdissection and next-generation sequencing, was used to identify chromosomal breakpoints. Long-range PCR and Sanger sequencing were used to precisely characterize 100 breakpoints in 50 ABCR carriers.

RESULTS

In addition to the recurrent regions of balanced rearrangement breaks in 8q24.13, 11q11.23, and 22q11.21 that had been documented, we have discovered a 10-Mb region of 12q24.13-q24.3 that could potentially be a sparse region of balanced rearrangement breaks. We found that 898 breakpoints caused gene disruption and a total of 188 breakpoints interrupted genes recorded in OMIM. The percentage of breakpoints that disrupted autosomal dominant genes recorded in OMIM was 25.53% (48/188). Fifty-four of the precisely characterized breakpoints had 1-8-bp microhomologous sequences.

CONCLUSION

Our findings provide a reference for the evaluation of the pathogenicity of mutations in related genes that cause protein truncation in clinical practice. According to the characteristics of breakpoints, non-homologous end joining and microhomology-mediated break-induced replication may be the main mechanism for ABCRs formation.

Break-induced replication underlies formation of inverted triplications and generates unexpected diversity in haplotype structures

Background

The duplication-triplication/inverted-duplication (DUP-TRP/INV-DUP) structure is a type of complex genomic rearrangement (CGR) hypothesized to result from replicative repair of DNA due to replication fork collapse. It is often mediated by a pair of inverted low-copy repeats (LCR) followed by iterative template switches resulting in at least two breakpoint junctions in cis . Although it has been identified as an important mutation signature of pathogenicity for genomic disorders and cancer genomes, its architecture remains unresolved and is predicted to display at least four structural variation (SV) haplotypes.

Results

Here we studied the genomic architecture of DUP-TRP/INV-DUP by investigating the genomic DNA of 24 patients with neurodevelopmental disorders identified by array comparative genomic hybridization (aCGH) on whom we found evidence for the existence of 4 out of 4 predicted SV haplotypes. Using a combination of short-read genome sequencing (GS), long- read GS, optical genome mapping and StrandSeq the haplotype structure was resolved in 18 samples. This approach refined the point of template switching between inverted LCRs in 4 samples revealing a DNA segment of ∼2.2-5.5 kb of 100% nucleotide similarity. A prediction model was developed to infer the LCR used to mediate the non-allelic homology repair.

Conclusions

These data provide experimental evidence supporting the hypothesis that inverted LCRs act as a recombinant substrate in replication-based repair mechanisms. Such inverted repeats are particularly relevant for formation of copy-number associated inversions, including the DUP-TRP/INV-DUP structures. Moreover, this type of CGR can result in multiple conformers which contributes to generate diverse SV haplotypes in susceptible loci .

Combination of trio-based whole exome sequencing and optical genome mapping reveals a cryptic balanced translocation that causes unbalanced chromosomal rearrangements in a family with multiple anomalies

Background: Balanced translocation (BT) carriers can produce imbalanced gametes and experience recurrent spontaneous abortions (RSAs) and even give birth to a child with complex chromosomal disorders. Here, we report a cryptic BT, t(5; 6) (p15.31; p25.1), in the proband's grandmother, which caused unbalanced chromosomal rearrangements and various anomalies in the two subsequent generations. We also provide a thorough overview of the application of optical genome mapping (OGM) to identify chromosomal structural variants (SVs). Methods: Trio-based whole exome sequencing (Trio-WES) was conducted to explore the genetic basis of the phenotype of the proband and her mother. High-resolution karyotype analysis and OGM detection were performed on the proband's grandparents to trace the origin of the unbalanced rearrangements between chromosomes 5 and 6. A PubMed search was conducted with the following keywords: "OGM" and "SVs." Then, relevant studies were collected and systematically reviewed. Results: The proband and her mother presented with various anomalies, whereas the grandmother was healthy but had a history of four abnormal pregnancies. Trio-WES revealed a heterozygous duplication on the terminal region of chromosome 5p and a heterozygous deletion on the proximal end of chromosome 6p in the proband and her mother. High-resolution karyotype analysis revealed no aberrant karyotypes in either grandparent, whereas OGM detection revealed a cryptic BT, t(5; 6)(p15.31; p25.1), in the proband's grandmother. An overwhelming majority of research publications have verified the clinical utility of OGM in detecting SVs. Conclusion: The results of this study revealed that the unbalanced chromosomal rearrangements and many anomalies observed in multiple members of the family were attributable to the cryptic BT carried by the proband's grandmother. This study supports that OGM has a unique advantage for detecting cryptic BTs, and can be used as a first-tier genetic test for the etiological diagnosis of infertility, RSAs, and other complex genetic disorders.

Analysis of chromosomal structural variations in patients with recurrent spontaneous abortion using optical genome mapping

Background and aims: Certain chromosomal structural variations (SVs) in biological parents can lead to recurrent spontaneous abortions (RSAs). Unequal crossing over during meiosis can result in the unbalanced rearrangement of gamete chromosomes such as duplication or deletion. Unfortunately, routine techniques such as karyotyping, fluorescence in situ hybridization (FISH), chromosomal microarray analysis (CMA), and copy number variation sequencing (CNV-seq) cannot detect all types of SVs. In this study, we show that optical genome mapping (OGM) quickly and accurately detects SVs for RSA patients with a high resolution and provides more information about the breakpoint regions at gene level. Methods: Seven couples who had suffered RSA with unbalanced chromosomal rearrangements of aborted embryos were recruited, and ultra-high molecular weight (UHMW) DNA was isolated from their peripheral blood. The consensus genome map was created by de novo assembly on the Bionano Solve data analysis software. SVs and breakpoints were identified via alignments of the reference genome GRCh38/hg38. The exact breakpoint sequences were verified using either Oxford Nanopore sequencing or Sanger sequencing. Results: Various SVs in the recruited couples were successfully detected by OGM. Also, additional complex chromosomal rearrangement (CCRs) and four cryptic balanced reciprocal translocations (BRTs) were revealed, further refining the underlying genetic causes of RSA. Two of the disrupted genes identified in this study, FOXK2 [46,XY,t(7; 17)(q31.3; q25)] and PLXDC2 [46,XX,t(10; 16)(p12.31; q23.1)], had been previously shown to be associated with male fertility and embryo transit. Conclusion: OGM accurately detects chromosomal SVs, especially cryptic BRTs and CCRs. It is a useful complement to routine human genetic diagnostics, such as karyotyping, and detects cryptic BRTs and CCRs more accurately than routine genetic diagnostics.

Resolving complex structural variants via nanopore sequencing

The recent development of high-throughput sequencing platforms provided impressive insights into the field of human genetics and contributed to considering structural variants (SVs) as the hallmark of genome instability, leading to the establishment of several pathologic conditions, including neoplasia and neurodegenerative and cognitive disorders. While SV detection is addressed by next-generation sequencing (NGS) technologies, the introduction of more recent long-read sequencing technologies have already been proven to be invaluable in overcoming the inaccuracy and limitations of NGS technologies when applied to resolve wide and structurally complex SVs due to the short length (100-500 bp) of the sequencing read utilized. Among the long-read sequencing technologies, Oxford Nanopore Technologies developed a sequencing platform based on a protein nanopore that allows the sequencing of "native" long DNA molecules of virtually unlimited length (typical range 1-100 Kb). In this review, we focus on the bioinformatics methods that improve the identification and genotyping of known and novel SVs to investigate human pathological conditions, discussing the possibility of introducing nanopore sequencing technology into routine diagnostics.

Breakpoint analysis for cytogenetically balanced translocation revealed unexpected complex structural abnormalities and suggested the position effect for MEF2C

Many disease-causing genes have been identified by determining the breakpoints of balanced chromosomal translocations. Recent progress in genomic analysis has accelerated the analysis of chromosomal translocation-breakpoints at the nucleotide level. Using a long-read whole-genome sequence, we analyzed the breakpoints of the cytogenetically balanced chromosomal translocation t(5;15)(q21;26.3), which was confirmed to be of de novo origin, in a patient with a neurodevelopmental disorder. The results showed complex rearrangements with seven fragments consisting of five breakpoint-junctions (BJs). Four of the five BJs showed microhomologies of 1-3-bp, and only one BJ displayed a signature of blunt-end ligation, indicating chromothripsis as the underlying mechanism. Although the BJs did not disrupt any disease-causing gene, the clinical features of the patient were compatible with MEF2C haploinsufficiency syndrome. Complex rearrangements were located approximately 2.5-Mb downstream of MEF2C. Therefore, position effects were considered the mechanism of the occurrence of MEF2C haploinsufficiency syndrome.

Nanopore sequencing for detecting reciprocal translocation carrier status in preimplantation genetic testing

BACKGROUND

Balanced reciprocal translocation (BRT) is one of the most common chromosomal abnormalities that causes infertility, recurrent miscarriage, and birth defects. Preimplantation genetic testing (PGT) is widely used to select euploid embryos for BRT carriers to increase the chance of a healthy live birth. Several strategies can be used to distinguish reciprocal translocation carrier embryos from those with a normal karyotype; however, these techniques are time-consuming and difficult to implement in clinical laboratories. In this study, nanopore sequencing was performed in two reciprocal translocation carriers, and the results were validated using the next-generation sequencing-based method named, "Mapping Allele with Resolved Carrier Status" (MaReCs).

RESULTS

The translocation breakpoints in both reciprocal translocation carriers were accurately identified by nanopore sequencing and were in accordance with the results obtained using MaReCs. More than one euploid non-balanced translocation carrier embryo was identified in both patients. Amniocentesis results revealed normal karyotypes, consistent with the findings by MaReCs and nanopore sequencing.

CONCLUSION

Our results suggest that nanopore sequencing is a powerful strategy for accurately distinguishing non-translocation embryos from translocation carrier embryos and precisely localizing translocation breakpoints, which is essential for PGT and aids in reducing the propagation of reciprocal translocation in the population.

Exploring the Genetic Causality of Discordant Phenotypes in Familial Apparently Balanced Translocation Cases Using Whole Exome Sequencing

Familial apparently balanced translocations (ABTs) are usually not associated with a phenotype; however, rarely, ABTs segregate with discordant phenotypes in family members carrying identical rearrangements. The current study was a follow-up investigation of four familial ABTs, where whole exome sequencing (WES) was implemented as a diagnostic tool to identify the underlying genetic aetiology of the patients' phenotypes. Data were analysed using an in-house bioinformatics pipeline alongside VarSome Clinical. WES findings were validated with Sanger sequencing, while the impact of splicing and missense variants was assessed by reverse-transcription PCR and in silico tools, respectively. Novel candidate variants were identified in three families. In family 1, it was shown that the de novo pathogenic STXBP1 variant (NM_003165.6:c.1110+2T>G) affected splicing and segregated with the patient's phenotype. In family 2, a likely pathogenic TUBA1A variant (NM_006009.4:c.875C>T, NP_006000.2:p.(Thr292Ile)) could explain the patient's symptoms. In family 3, an SCN1A variant of uncertain significance (NM_006920.6:c.5060A>G, NP_008851.3:p.(Glu1687Gly)) required additional evidence to sufficiently support causality. This first report of WES application in familial ABT carriers with discordant phenotypes supported our previous findings describing such rearrangements as coincidental. Thus, WES can be recommended as a complementary test to find the monogenic cause of aberrant phenotypes in familial ABT carriers.

Expanding the genetic spectrum of the pyruvate carboxylase deficiency with novel missense, deep intronic and structural variants

Introduction

Pathogenic variants in the pyruvate carboxylase (PC) gene cause a wide spectrum of recessive phenotypes, ranging from the early-onset fatal encephalopathy to the adult-onset benign form.

Results

Patient 1 is a 6 y.o. boy with ataxia, hypoglycemia and episodes of lactic acidosis. WGS revealed the novel heterozygous missense variant c.1372A > G (p.Asn458Asp) in the PC gene. Additional analysis revealed discordant reads mapped to chromosomes 11 and 1, so a reciprocal translocation disrupted the PC gene was suspected. The translocation was validated via FISH-analysis and Sanger sequencing of its boundaries.

Patient 2 is a 13 y.o. girl with psychomotor delay, episodes of lactic acidosis and ketonuria. WES revealed the novel homozygous intronic variant c.1983-116C > T. The PC's mRNA analysis demonstrated the exonization of several intron 16 sequences and some residual amount of WT mRNA isoform.

Two other patients had more severe course of the disease. Their genotype represents missense variants in compound heterozygous and homozygous state (c.1876C > T (p.Arg626Trp), c.2606G > C (p.Gly869Ala), c.2435C > A (p.Ala812Asp).

Conclusion

In patients with metabolic crises, lactic acidosis and hypoglycemia analysis of PC gene is recommended. WGS with deep bioinformatic analysis should be taken into consideration when none or the only one pathogenic variant in the PC gene is found.

•

In patients with metabolic crises with lactic acidosis and hypoglycemia WGS with deep bioinformatic analysis should be taken when none or the only one heterozygous pathogenic variant in the PC gene is found.

•

PC gene may have more deep intronic pathogenic variants than included in known mutation databases.

Multi-Omic Investigations of a 17-19 Translocation Links MINK1 Disruption to Autism, Epilepsy and Osteoporosis

Balanced structural variants, such as reciprocal translocations, are sometimes hard to detect with sequencing, especially when the breakpoints are located in repetitive or insufficiently mapped regions of the genome. In such cases, long-range information is required to resolve the rearrangement, identify disrupted genes and, in symptomatic carriers, pinpoint the disease-causing mechanisms. Here, we report an individual with autism, epilepsy and osteoporosis and a de novo balanced reciprocal translocation: t(17;19) (p13;p11). The genomic DNA was analyzed by short-, linked- and long-read genome sequencing, as well as optical mapping. Transcriptional consequences were assessed by transcriptome sequencing of patient-specific neuroepithelial stem cells derived from induced pluripotent stem cells (iPSC). The translocation breakpoints were only detected by long-read sequencing, the first on 17p13, located between exon 1 and exon 2 of MINK1 (Misshapen-like kinase 1), and the second in the chromosome 19 centromere. Functional validation in induced neural cells showed that MINK1 expression was reduced by >50% in the patient’s cells compared to healthy control cells. Furthermore, pathway analysis revealed an enrichment of changed neural pathways in the patient’s cells. Altogether, our multi-omics experiments highlight MINK1 as a candidate monogenic disease gene and show the advantages of long-read genome sequencing in capturing centromeric translocations.

Destabilizing Effects of Ionizing Radiation on Chromosomes: Sizing up the Damage

For long-term survival and evolution, all organisms have depended on a delicate balance between processes involved in maintaining stability of their genomes and opposing processes that lead toward destabilization. At the level of mammalian somatic cells in renewal tissues, events or conditions that can tip this balance toward instability have attracted special interest in connection with carcinogenesis. Mutations affecting DNA (and its subsequent repair) would, of course, be a major consideration here. These may occur spontaneously through endogenous cellular processes or as a result of exposure to mutagenic environmental agents. It is in this context that we discuss the rather unique destabilizing effects of ionizing radiation (IR) in terms of its ability to cause large-scale structural rearrangements to the genome. We present arguments supporting the conclusion that these and other important effects of IR originate largely from microscopically visible chromosome aberrations.

Deciphering balanced translocations in infertile males by next-generation sequencing to identify candidate genes for spermatogenesis disorders

Male infertility affects about 7% of the general male population. Balanced structural chromosomal rearrangements are observed in 0.4-1.4% of infertile males and are considered as a well-established cause of infertility. However, underlying pathophysiological mechanisms still need to be clarified. A strategy combining standard and high throughput cytogenetic and molecular technologies was applied in order to identify the candidate genes that might be implicated in the spermatogenesis defect in three male carriers of different balanced translocations. Fluorescence in situ hybridization (FISH) and whole-genome paired-end sequencing were used to characterize translocation breakpoints at the molecular level while exome sequencing was performed in order to exclude the presence of any molecular event independent from the chromosomal rearrangement in the patients. All translocation breakpoints were characterized in the three patients. We identified four variants: a position effect on LACTB2 gene in Patient 1, a heterozygous CTDP1 gene disruption in Patient 2, two single-nucleotide variations (SNVs) in DNAH5 gene and a heterozygous 17q12 deletion in Patient 3. The variants identified in this study need further validation to assess their roles in male infertility. This study shows that beside the mechanical effect of structural rearrangement on meiosis, breakpoints could result in additional alterations such as gene disruption or position effect. Moreover, additional SNVs or copy number variations may be fortuitously present and could explain the variable impact of chromosomal rearrangements on spermatogenesis. In conclusion, this study confirms the relevance of combining different cytogenetic and molecular techniques to investigate patients with spermatogenesis disorders and structural rearrangements on genomic scale.

Classical, Molecular, and Genomic Cytogenetics of the Pig, a Clinical Perspective

The chromosomes of the domestic pig (Sus scrofa domesticus) are known to be prone to reciprocal chromosome translocations and other balanced chromosome rearrangements with concomitant fertility impairment of carriers. In response to the remarkable prevalence of chromosome rearrangements in swine herds, clinical cytogenetics laboratories have been established in several countries in order to screen young boars for chromosome rearrangements prior to service. At present, clinical cytogenetics laboratories typically apply classical cytogenetics techniques such as giemsa-trypsin (GTG)-banding to produce high-quality karyotypes and reveal large-scale chromosome ectopic exchanges. Further refinements to clinical cytogenetics practices have led to the implementation of molecular cytogenetics techniques such as fluorescent in-situ hybridization (FISH), allowing for rearrangements to be visualized and breakpoints refined using fluorescently labelled painting probes. The next-generation of clinical cytogenetics include the implementation of DNA microarrays, and next-generation sequencing (NGS) technologies such as DNA sequencing to better explore tentative genome architecture changes. The implementation of these cytogenomics techniques allow the genomes of rearrangement carriers to be deciphered at the highest resolution, allowing rearrangements to be detected; breakpoints to be delineated; and, most importantly, potential gene implications of those chromosome rearrangements to be interrogated. Clinical cytogenetics has become an integral tool in the livestock industry, identifying rearrangements and allowing breeders to make informed breeding decisions.

Robust Performance of Potentially Functional SNPs in Machine Learning Models for the Prediction of Atorvastatin-Induced Myalgia

Statins can cause muscle symptoms resulting in poor adherence to therapy and increased cardiovascular risk. We hypothesize that combinations of potentially functional SNPs (pfSNPs), rather than individual SNPs, better predict myalgia in patients on atorvastatin. This study assesses the value of potentially functional single nucleotide polymorphisms (pfSNPs) and employs six machine learning algorithms to identify the combination of SNPs that best predict myalgia.

Methods: Whole genome sequencing of 183 Chinese, Malay and Indian patients from Singapore was conducted to identify genetic variants associated with atorvastatin induced myalgia. To adjust for confounding factors, demographic and clinical characteristics were also examined for their association with myalgia. The top factor, sex, was then used as a covariate in the whole genome association analyses. Variants that were highly associated with myalgia from this and previous studies were extracted, assessed for potential functionality (pfSNPs) and incorporated into six machine learning models. Predictive performance of a combination of different models and inputs were compared using the average cross validation area under ROC curve (AUC). The minimum combination of SNPs to achieve maximum sensitivity and specificity as determined by AUC, that predict atorvastatin-induced myalgia in most, if not all the six machine learning models was determined.

Results: Through whole genome association analyses using sex as a covariate, a larger proportion of pfSNPs compared to non-pf SNPs were found to be highly associated with myalgia. Although none of the individual SNPs achieved genome wide significance in univariate analyses, machine learning models identified a combination of 15 SNPs that predict myalgia with good predictive performance (AUC >0.9). SNPs within genes identified in this study significantly outperformed SNPs within genes previously reported to be associated with myalgia. pfSNPs were found to be more robust in predicting myalgia, outperforming non-pf SNPs in the majority of machine learning models tested.

Conclusion: Combinations of pfSNPs that were consistently identified by different machine learning models to have high predictive performance have good potential to be clinically useful for predicting atorvastatin-induced myalgia once validated against an independent cohort of patients.

Resolving Breakpoints of Chromosomal Rearrangements at the Nucleotide Level Using Sanger Sequencing

Novel cytogenetic tools are increasingly based on genome sequencing for detecting chromosomal abnormalities. Different sequence-based techniques optimized for diagnosis of structural variants can be useful for narrowing down the localization of breakpoints of chromosomal abnormalities, but do not offer nucleotide resolution of breakpoints for proper interpretation of gene disruption. This protocol presents the characterization of structural variants at nucleotide resolution using Sanger sequencing after low-pass large-insert genome sequencing or other long-molecule methods. © 2020 Wiley Periodicals LLC. Basic Protocol 1: Primer design for junction amplification at translocations and inversions Basic Protocol 2: Amplification of derivative chromosomes using a long-range polymerase Alternate Protocol: Amplification of derivative chromosomes using a hot-start polymerase Basic Protocol 3: Preparation of DNA for Sanger sequencing Basic Protocol 4: Interpretation and reporting of breakpoints based on Sanger sequencing.

Cytogenetically visible inversions are formed by multiple molecular mechanisms

Cytogenetically detected inversions are generally assumed to be copy number and phenotypically neutral events. While nonallelic homologous recombination is thought to play a major role, recent data suggest the involvement of other molecular mechanisms in inversion formation. Using a combination of short-read whole-genome sequencing (WGS), 10X Genomics Chromium WGS, droplet digital polymerase chain reaction and array comparative genomic hybridization we investigated the genomic structure of 18 large unique cytogenetically detected chromosomal inversions and achieved nucleotide resolution of at least one chromosomal inversion junction for 13/18 (72%). Surprisingly, we observed that seemingly copy number neutral inversions can be accompanied by a copy-number gain of up to 350 kb and local genomic complexities (3/18, 17%). In the resolved inversions, the mutational signatures are consistent with nonhomologous end-joining (8/13, 62%) or microhomology-mediated break-induced replication (5/13, 38%). Our study indicates that short-read 30x coverage WGS can detect a substantial fraction of chromosomal inversions. Moreover, replication-based mechanisms are responsible for approximately 38% of those events leading to a significant proportion of inversions that are actually accompanied by additional copy-number variation potentially contributing to the overall phenotypic presentation of those patients.

TSC1 intragenic deletion transmitted from a mosaic father to two siblings with cardiac rhabdomyomas: Identification of two aberrant transcripts

Tuberous sclerosis complex (TSC) is a rare autosomal dominant disorder characterized by non-cancerous tumors in multiple organs including the brain, kidney, lung, heart, and skin. We encountered a Japanese family consisting of two siblings (a four-year-old boy and a one-year-old girl) with multiple cardiac rhabdomyomas conveying a high risk of TSC and apparently unaffected sibling (a two-year-old girl) and parents. Whole exome sequencing and application of Integrative Genomic Viewer revealed an identical intragenic TSC1 deletion with the breakpoints on intron 15 and exon 19 in the affected siblings, but not in the apparently unaffected sibling and parents. Subsequently, PCR-based analyses were performed using primers flanking the deletion, showing that the deletion was also present in the father and that the deletion occurred between chr9:135,777,038 (bp) and chr9:135,780,540 (bp) in association with a one bp overlap. Furthermore, RT-PCR analyses were carried out using lymphoblastoid cell lines, revealing a major in-frame insertion/deletion transcript produced by aberrant splicing using a cryptic ″ag″ splice acceptor motif at intron 15 (r.1998_2438delinsTTCATTAGGTGG) and a minor frameshift transcript generated by aberrant splicing between exon 15 and exon 20 (r.1998_2502del, p.Lys666Asnfs*15) in the affected siblings. These findings imply that the intragenic deletion producing two aberrant transcripts was generated as a somatic pathogenic variant involving the germline in the father and was transmitted to the affected siblings.

Whole Genome Low-Coverage Sequencing Concurrently Detecting Copy Number Variations and Their Underlying Complex Chromosomal Rearrangements by Systematic Breakpoint Mapping in Intellectual Deficiency/Developmental Delay Patients

Simple copy number variations (CNVs) detected by chromosomal microarray (CMA) can result from complex structural changes. Therefore, it is necessary to characterize potential structural changes that cause pathogenic CNVs. We applied whole-genome low-coverage sequencing (WGLCS) to concurrently detect pathogenic CNVs and their associated chromosomal rearrangements in 15 patients. All the patients had an average of 2–3 pathogenic CNVs involving 1–2 chromosomes. WGLCS identified all the 34 pathogenic CNVs found by microarray. By identifying chimeric read pairs, WGLCS mapped 70 breakpoints in these patients, of which 47 were finely mapped at the nucleotide level and confirmed by subsequent PCR amplification and Sanger sequencing of the junction fragments. In 15 patients, structural rearrangements were defined at molecular level in 13 patients. In 13 patients, WGLCS reveal no additional results in two patients. In another 11 patients, WGLCS revealed new breakpoints or finely mapped the genes disrupted by breakpoints or 1–6 bp microhomology and/or short insertion (4–70 bp) in the breakpoints junctions. However, structural changes in the other two patients still remained unclear after WGLCS was performed. The structural alteration identified in the 13 patients could be divided into the following categories: (1) interstitial inverted duplication with concomitant terminal deletion (inv dup del) (P1,P4,P9,P11); (2) the product of pericentric inversion (P5); (3) ring chromosome (P8); (4) interstitial duplication and/or triplication (P6, P7); and (5) +der(22)t(11;22) (P2,P15); (6) complex structural rearrangements (P3,P12,P14). WGLCS displayed the ability to discover CNVs and define breakpoints and its disrupted genes and its surrounding sequences in one experiment at base-pair-resolution, which help us to learn more about the mechanisms of formation of observed genomic rearrangements, and in which DNA replicative/repair mechanism might contribute to the formation of complex rearrangements in 11 patients. Clear karyotype at molecular level could help provide an accurate evaluation of recurrent risk and guide prenatal diagnosis or reproductive planning.

Targeted next-generation sequencing identifies the disruption of the SHANK3 and RYR2 genes in a patient carrying a de novo t(1;22)(q43;q13.3) associated with signs of Phelan-McDermid syndrome

Background

It has been known for more than 30 years that balanced translocations, especially if de novo, can associate with congenital malformations and / or neurodevelopmental disorders, following the disruption of a disease gene or its cis-regulatory elements at one or both breakpoints.

Case presentation

We describe a 10-year-old girl with a non-specific neurodevelopmental disorder characterized by moderate intellectual disability (ID), gross motor clumsiness, social and communication deficits. She carries a de novo reciprocal translocation between chromosomes 1q43 and 22q13.3, the latter suggesting the involvement of SHANK3. Indeed, its haploinsufficiency associates with Phelan-McDermid Syndrome, whose main symptoms are characterized by global developmental delay and absent or severely delayed expressive speech. A deep molecular approach, including next-generation sequencing of SHANK3 locus, allowed demonstrating the breakage of RYR2 and SHANK3 on the derivative chromosomes 1 and 22 respectively, and the formation of two fusion genes SHANK3-RYR2 and RYR2-SHANK3 with concomitant cryptic deletion of 3.6 and 4.1 kilobases at translocation junction of both derivatives chromosomes 22 and 1, respectively.

Conclusions

Although the interruption of SHANK3 accounts for the patient’s psychomotor retardation and autism-like behavior, we do not exclude that the interruption of RYR2 may also have a role on her disorder, or result in further pathogenicity in the future. Indeed, RYR2 that has a well-established role in the etiology of two autosomal dominant adulthood cardiac disorders (#600996 and #604772) is also expressed in the brain (cerebellum, hippocampus, and cerebral cortex) and about half of RYR2 mutation carriers present late onset primary generalized epilepsy without cardiac arrhythmogenic disorders. Moreover, RYR2 variants have also been sporadically reported in individuals with early onset schizophrenia or ID, and its constraint values suggest intolerance to loss-of-function. This study not only confirms the usefulness of the molecular mapping of de novo balanced rearrangements in symptomatic individuals, but also underscores the need for long-term clinical evaluation of the patients, for better evaluating the pathogenicity of the chromosomal breakpoints.

Whole-genome mate-pair sequencing of apparently balanced chromosome rearrangements reveals complex structural variations: two case studies

Background

Apparently balanced chromosome rearrangements (ABCRs) in non-affected individuals are well-known to possess high reproductive risks such as infertility, abnormal offspring, and pregnancy loss. However, caution should be exercised in genetic counseling and reproductive intervention because cryptic unbalanced defects and genome structural variations beyond the resolution of routine cytogenetics may not be detected.

Case presentation

Here, we studied two familial cases of ABCRs were recruited in this study. In family 1, the couple suffered two abortions pregnancies and underwent labor induction. Single nucleotide polymorphism (SNP) array analysis of the aborted sample from the second pregnancy revealed a 10.8 Mb heterozygous deletion at 10q26.13q26.3 and a 5.5 Mb duplication at 19q13.41-q13.43. The non-affected father was identified as a carrier of three-way complex chromosomal rearrangement [t (6;10;19)(p22;q26;q13)] by karyotyping. Whole-genome mate-pair sequencing revealed a cryptic breakpoint on the derivative chromosome 19 (der19), indicating that the karyotype was a more complex structural rearrangement comprising four breakpoints. Three genes, FAM24B, CACNG8, and KIAA0556, were disrupted without causing any abnormal phenotype in the carrier. In family 2, the couple suffered from a spontaneous miscarriage. This family had an affected child with multiple congenital deformities and an unbalanced karyotype, 46,XY,der (11) t (6;11)(q13;p11.2). The female partner was identified as a balanced translocation carrier with the karyotype 46,XX,t (6;11)(q13;p11.2) dn. Further SNP array and fluorescent in situ hybridization (FISH) indicated a cryptic insertion between chromosome 6 and chromosome 11. Finally, whole-genome mate-pair sequencing revealed an extremely complex genomic structural variation, including a cryptic deletion and 12 breakpoints on chromosome 11, and 1 breakpoint on chromosome 6 .

Conclusions

Our study investigated two rare cases of ABCRs and demonstrated the efficacy of whole-genome mate-pair sequencing in analyzing the genome complex structural variation. In case of ABCRs detected by conventional cytogenetic techniques, whole genome sequencing (WGS) based approaches should be considered for accurate diagnosis, effective genetic counseling, and correct reproductive intervention to avoid recurrence risks.

Whole-genome sequencing reveals complex chromosome rearrangement disrupting NIPBL in infant with Cornelia de Lange syndrome

Clinical laboratory diagnostic evaluation of the genomes of children with suspected genetic disorders, including chromosomal microarray and exome sequencing, cannot detect copy number neutral genomic rearrangements such as inversions, balanced translocations, and complex chromosomal rearrangements (CCRs). We describe an infant with a clinical diagnosis of Cornelia de Lange syndrome (CdLS) in whom chromosome analysis revealed a de novo complex balanced translocation, 46,XY,t(5;7;6)(q11.2;q32;q13)dn. Subsequent molecular characterization by whole-genome sequencing (WGS) identified 23 breakpoints, delineating segments derived from four chromosomes (5;6;7;21) in ancestral or inverted orientation. One of the breakpoints disrupted a known CdLS gene, NIPBL. Further investigation revealed paternal origin of the CCR allele, clustering of the breakpoint junctions, and molecular repair signatures suggestive of a single catastrophic event. Notably, very short DNA segments (25 and 41 bp) were included in the reassembled chromosomes, lending additional support that the DNA repair machinery can detect and repair such segments. Interestingly, there was an independent paternally derived miniscule complex rearrangement, possibly predisposing to subsequent genomic instability. In conclusion, we report a CCR causing a monogenic Mendelian disorder, urging WGS analysis of similar unsolved cases with suspected Mendelian disorders. Breakpoint analysis allowed for identification of the underlying molecular diagnosis and implicated chromoanagenesis in CCR formation.

NCALD affects drug resistance and prognosis by acting as a ceRNA of CX3CL1 in ovarian cancer

Drug resistance, an impenetrable barrier in the treatment of ovarian cancer (OC), is often associated with poor outcomes. Hence, it is urgent to discover new factors controlling drug resistance and survival. The association between neurocalcin delta (NCALD) and cancer drug resistance is poorly understood. Here, we reveal that NCALD messenger RNA expression, probably regulated by DNA methylation and microRNAs, was significantly downregulated in at least three independent microarrays covering 633 ovarian carcinomas and 16 normal controls, which includes the Cancer Genome Atlas (TCGA) ovarian cohort. In the sub-groups of the TCGA cohort, NCALD was suppressed in 90 platinum-resistant tissues vs in 197 sensitive tissues. It is consistent with the quantitative reverse transcription polymerase chain reaction results revealing gene downregulation in carboplatin-resistant SKOV3 and HeyA8 OC cells as compared with that in controls. Low expression of NCALD predicted poor overall survival (OS) in sub-groups of 1656 patients, progression-free survival (PFS) in 1435 patients, and post-progression survival (PPS) in 782 patients according to Kaplan-Meier plotter covering 1815 OC patients. Comprehensive bioinformatic analyses strongly implicated NCALD in the regulation of drug resistance, probably via competing for endogenous RNA (ceRNA) interactions with CX3CL1 and tumor immune-microenvironment. NCALD acted as a ceRNA for CX3CL1 in 21 different cancers includes OC according to Starbase. These two genes negatively correlated with tumor purity and positively correlated with infiltration levels of neutrophils and dendritic cells in OC. The combined low expression of NCALD and CX3CL1 showed better prognosis potential for OS, PFS, and PPS in the 1815 OC patients than any of the individually tested genes. In summary, NCALD acts as a ceRNA for CX3CL1, and its downregulation may affect drug resistance and prognosis in OC. Thus, NCALD could be a new therapeutic target for anticancer therapy and a new biomarker for survival prediction in OC.

Analysis of balanced reciprocal translocations in patients with subfertility using single-molecule optical mapping

PURPOSE

Approximately 1% of individuals who carry a balanced reciprocal translocation (BRT) are subfertile. Current karyotyping does not have the resolution to determine whether the breakpoints of the involved chromosomes perturb genes important for fertility. The aim of this study was to apply single-molecule optical mapping (SMOM) to patients presenting for IVF (in vitro fertilization) to ascertain whether the BRT disrupted any genes associated with normal fertility.

METHODS

Nine subfertile patients with different BRTs were recruited for the study. Methyltransferase enzyme DLE1 was used to fluorescently label their genomic DNA samples at the recognition motif CTTAAG. The SMOM was performed on the Bionano platform, and long molecules aligned against the reference genome hg19 to identify the breakpoint regions. Mate-pair and PCR-Sanger sequencing were used to confirm the precise breakpoint sequences.

RESULTS

Both breakpoint regions in each of the nine BRTs were finely mapped to small regions of approximately 10 Kb, and their positions were consistent with original cytogenetic banding patterns determined by karyotyping. In three BRTs, breakpoints disrupted genes known to be associated with male infertility, namely NUP155 and FNDC3A [46,XY,t(5;13)(p15;q22)], DPY19L1 [46,XY,t(1;7)(p36.3;p15), and BAI3 [46,XY,t(3;6)(p21;q16)].

CONCLUSIONS

The SMOM has potential clinical application as a rapid tool to screen patients with BRTs for underlying genetic causes of infertility and other diseases.

Genome sequencing in cytogenetics: Comparison of short-read and linked-read approaches for germline structural variant detection and characterization

BACKGROUND

Structural variants (SVs) include copy number variants (CNVs) and apparently balanced chromosomal rearrangements (ABCRs). Genome sequencing (GS) enables SV detection at base-pair resolution, but the use of short-read sequencing is limited by repetitive sequences, and long-read approaches are not yet validated for diagnosis. Recently, 10X Genomics proposed Chromium, a technology providing linked-reads to reconstruct long DNA fragments and which could represent a good alternative. No study has compared short-read to linked-read technologies to detect SVs in a constitutional diagnostic setting yet. The aim of this work was to determine whether the 10X Genomics technology enables better detection and comprehension of SVs than short-read WGS.

METHODS

We included 13 patients carrying various SVs. Whole genome analyses were performed using paired-end HiSeq X sequencing with (linked-read strategy) or without (short-read strategy) Chromium library preparation. Two different bioinformatic pipelines were used: Variants are called using BreakDancer for short-read strategy and LongRanger for long-read strategy. Variant interpretations were first blinded.

RESULTS

The short-read strategy allowed diagnosis of known SV in 10/13 patients. After unblinding, the linked-read strategy identified 10/13 SVs, including one (patient 7) missed by the short-read strategy.

CONCLUSION

In conclusion, regarding the results of this study, 10X Genomics solution did not improve the detection and characterization of SV.

Development of coupling controlled polymerizations by adapter-ligation in mate-pair sequencing for detection of various genomic variants in one single assay

The diversity of disease presentations warrants one single assay for detection and delineation of various genomic disorders. Herein, we describe a gel-free and biotin-capture-free mate-pair method through coupling Controlled Polymerizations by Adapter-Ligation (CP-AL). We first demonstrated the feasibility and ease-of-use in monitoring DNA nick translation and primer extension by limiting the nucleotide input. By coupling these two controlled polymerizations by a reported non-conventional adapter-ligation reaction 3′ branch ligation, we evidenced that CP-AL significantly increased DNA circularization efficiency (by 4-fold) and was applicable for different sequencing methods but at a faction of current cost. Its advantages were further demonstrated by fully elimination of small-insert-contaminated (by 39.3-fold) with a ∼50% increment of physical coverage, and producing uniform genome/exome coverage and the lowest chimeric rate. It achieved single-nucleotide variants detection with sensitivity and specificity up to 97.3 and 99.7%, respectively, compared with data from small-insert libraries. In addition, this method can provide a comprehensive delineation of structural rearrangements, evidenced by a potential diagnosis in a patient with oligo-atheno-terato-spermia. Moreover, it enables accurate mutation identification by integration of genomic variants from different aberration types. Overall, it provides a potential single-integrated solution for detecting various genomic variants, facilitating a genetic diagnosis in human diseases.

Genome Sequencing Explores Complexity of Chromosomal Abnormalities in Recurrent Miscarriage

Recurrent miscarriage (RM) affects millions of couples globally, and half of them have no demonstrated etiology. Genome sequencing (GS) is an enhanced and novel cytogenetic tool to define the contribution of chromosomal abnormalities in human diseases. In this study we evaluated its utility in RM-affected couples. We performed low-pass GS retrospectively for 1,090 RM-affected couples, all of whom had routine chromosome analysis. A customized sequencing and interpretation pipeline was developed to identify chromosomal rearrangements and deletions/duplications with confirmation by fluorescence in situ hybridization, chromosomal microarray analysis, and PCR studies. Low-pass GS yielded results in 1,077 of 1,090 couples (98.8%) and detected 127 chromosomal abnormalities in 11.7% (126/1,077) of couples; both members of one couple were identified with inversions. Of the 126 couples, 39.7% (50/126) had received former diagnostic results by karyotyping characteristic of normal human male or female karyotypes. Low-pass GS revealed additional chromosomal abnormalities in 50 (4.0%) couples, including eight with balanced translocations and 42 inversions. Follow-up studies of these couples showed a higher miscarriage/fetal-anomaly rate of 5/10 (50%) compared to 21/93 (22.6%) in couples with normal GS, resulting in a relative risk of 2.2 (95% confidence interval, 1.1 to 4.6). In these couples, this protocol significantly increased the diagnostic yield of chromosomal abnormalities per couple (11.7%) in comparison to chromosome analysis (8.0%, chi-square test p = 0.000751). In summary, low-pass GS identified underlying chromosomal aberrations in 1 in 9 RM-affected couples, enabling identification of a subgroup of couples with increased risk of subsequent miscarriage who would benefit from a personalized intervention.

From cytogenetics to cytogenomics: whole-genome sequencing as a first-line test comprehensively captures the diverse spectrum of disease-causing genetic variation underlying intellectual disability

Background

Since different types of genetic variants, from single nucleotide variants (SNVs) to large chromosomal rearrangements, underlie intellectual disability, we evaluated the use of whole-genome sequencing (WGS) rather than chromosomal microarray analysis (CMA) as a first-line genetic diagnostic test.

Methods

We analyzed three cohorts with short-read WGS: (i) a retrospective cohort with validated copy number variants (CNVs) (cohort 1, n = 68), (ii) individuals referred for monogenic multi-gene panels (cohort 2, n = 156), and (iii) 100 prospective, consecutive cases referred to our center for CMA (cohort 3). Bioinformatic tools developed include FindSV, SVDB, Rhocall, Rhoviz, and vcf2cytosure.

Results

First, we validated our structural variant (SV)-calling pipeline on cohort 1, consisting of three trisomies and 79 deletions and duplications with a median size of 850 kb (min 500 bp, max 155 Mb). All variants were detected. Second, we utilized the same pipeline in cohort 2 and analyzed with monogenic WGS panels, increasing the diagnostic yield to 8%. Next, cohort 3 was analyzed by both CMA and WGS. The WGS data was processed for large (> 10 kb) SVs genome-wide and for exonic SVs and SNVs in a panel of 887 genes linked to intellectual disability as well as genes matched to patient-specific Human Phenotype Ontology (HPO) phenotypes. This yielded a total of 25 pathogenic variants (SNVs or SVs), of which 12 were detected by CMA as well. We also applied short tandem repeat (STR) expansion detection and discovered one pathologic expansion in ATXN7. Finally, a case of Prader-Willi syndrome with uniparental disomy (UPD) was validated in the WGS data.

Important positional information was obtained in all cohorts. Remarkably, 7% of the analyzed cases harbored complex structural variants, as exemplified by a ring chromosome and two duplications found to be an insertional translocation and part of a cryptic unbalanced translocation, respectively.

Conclusion

The overall diagnostic rate of 27% was more than doubled compared to clinical microarray (12%). Using WGS, we detected a wide range of SVs with high accuracy. Since the WGS data also allowed for analysis of SNVs, UPD, and STRs, it represents a powerful comprehensive genetic test in a clinical diagnostic laboratory setting.

Expression, Localization, and Effect of High Salt Intake on Electroneutral Na+/HCO3 - Cotransporter NBCn2 in Rat Small Intestine: Implication in Intestinal NaCl Absorption

The electroneutral Na+/HCO3 - cotransporter NBCn2 (SLC4A10) of solute carrier family 4 (SLC4) plays important physiological and pathological roles in the body. Our previous study showed that NBCn2 is expressed on the protein level in the small intestine of rat. Here, by reverse-transcription polymerase chain reaction (PCR), we identified a novel full-length NBCn2 variant, i.e., NBCn2-K, from rat small intestine. By pHi measurement with Xenopus oocytes, the activity of NBCn2-K is not significantly different from NBCn2-G. By western blotting, NBCn2 and the Na+/H+ exchanger NHE3 (SLC9A3) are predominantly expressed in the jejunum of rat small intestine. By immunofluorescence, NBCn2 and NHE3 are localized at the apical domain of the jejunum. NaCl overload decreases the expression of NBCn2 by 56% and that of NHE3 by 40% in the small intestine. We propose that NBCn2 is involved in the transepithelial NaCl absorption in the small intestine, and that the down-regulation of NBCn2 by NaCl represents an adaptive response to high salt intake in rat.

Templated Insertions: A Smoking Gun for Polymerase Theta-Mediated End Joining

A recognized source of disease-causing genome alterations is erroneous repair of broken chromosomes, which can be executed by two distinct mechanisms: non-homologous end joining (NHEJ) and the recently discovered polymerase theta-mediated end joining (TMEJ) pathway. While TMEJ has previously been considered to act as an alternative mechanism backing up NHEJ, recent work points to a role for TMEJ in the repair of replication-associated DNA breaks that are excluded from repair through homologous recombination. Because of its mode of action, TMEJ is intrinsically mutagenic and sometimes leaves behind a recognizable genomic scar when joining chromosome break ends (i.e., 'templated insertions'). This review article focuses on the intriguing observation that this polymerase theta signature is frequently observed in disease alleles, arguing for a prominent role of this double-strand break repair pathway in genome diversification and disease-causing spontaneous mutagenesis in humans.

Copy number variant and runs of homozygosity detection by microarrays enabled more precise molecular diagnoses in 11,020 clinical exome cases

Background

Exome sequencing (ES) has been successfully applied in clinical detection of single nucleotide variants (SNVs) and small indels. However, identification of copy number variants (CNVs) using ES data remains challenging. The purpose of this study is to understand the contribution of CNVs and copy neutral runs of homozygosity (ROH) in molecular diagnosis of patients referred for ES.

Methods

In a cohort of 11,020 consecutive ES patients, an Illumina SNP array analysis interrogating mostly coding SNPs was performed as a quality control (QC) measurement and for CNV/ROH detection. Among these patients, clinical chromosomal microarray analysis (CMA) was performed at Baylor Genetics (BG) on 3229 patients, either before, concurrently, or after ES. We retrospectively analyzed the findings from CMA and the QC array.

Results

The QC array can detect ~ 70% of pathogenic/likely pathogenic CNVs (PCNVs) detectable by CMA. Out of the 11,020 ES cases, the QC array identified PCNVs in 327 patients and uniparental disomy (UPD) disorder-related ROH in 10 patients. The overall PCNV/UPD detection rate was 5.9% in the 3229 ES patients who also had CMA at BG; PCNV/UPD detection rate was higher in concurrent ES and CMA than in ES with prior CMA (7.2% vs 4.6%). The PCNVs/UPD contributed to the molecular diagnoses in 17.4% (189/1089) of molecularly diagnosed ES cases with CMA and were estimated to contribute in 10.6% of all molecularly diagnosed ES cases. Dual diagnoses with both PCNVs and SNVs were detected in 38 patients. PCNVs affecting single recessive disorder genes in a compound heterozygous state with SNVs were detected in 4 patients, and homozygous deletions (mostly exonic deletions) were detected in 17 patients. A higher PCNV detection rate was observed for patients with syndromic phenotypes and/or cardiovascular abnormalities.

Conclusions

Our clinical genomics study demonstrates that detection of PCNV/UPD through the QC array or CMA increases ES diagnostic rate, provides more precise molecular diagnosis for dominant as well as recessive traits, and enables more complete genetic diagnoses in patients with dual or multiple molecular diagnoses. Concurrent ES and CMA using an array with exonic coverage for disease genes enables most effective detection of both CNVs and SNVs and therefore is recommended especially in time-sensitive clinical situations.

Electronic supplementary material

The online version of this article (10.1186/s13073-019-0639-5) contains supplementary material, which is available to authorized users.

Whole genome paired-end sequencing elucidates functional and phenotypic consequences of balanced chromosomal rearrangement in patients with developmental disorders

BACKGROUND

Balanced chromosomal rearrangements associated with abnormal phenotype are rare events, but may be challenging for genetic counselling, since molecular characterisation of breakpoints is not performed routinely. We used next-generation sequencing to characterise breakpoints of balanced chromosomal rearrangements at the molecular level in patients with intellectual disability and/or congenital anomalies.

METHODS

Breakpoints were characterised by a paired-end low depth whole genome sequencing (WGS) strategy and validated by Sanger sequencing. Expression study of disrupted and neighbouring genes was performed by RT-qPCR from blood or lymphoblastoid cell line RNA.

RESULTS

Among the 55 patients included (41 reciprocal translocations, 4 inversions, 2 insertions and 8 complex chromosomal rearrangements), we were able to detect 89% of chromosomal rearrangements (49/55). Molecular signatures at the breakpoints suggested that DNA breaks arose randomly and that there was no major influence of repeated elements. Non-homologous end-joining appeared as the main mechanism of repair (55% of rearrangements). A diagnosis could be established in 22/49 patients (44.8%), 15 by gene disruption (KANSL1, FOXP1, SPRED1, TLK2, MBD5, DMD, AUTS2, MEIS2, MEF2C, NRXN1, NFIX, SYNGAP1, GHR, ZMIZ1) and 7 by position effect (DLX5, MEF2C, BCL11B, SATB2, ZMIZ1). In addition, 16 new candidate genes were identified. Systematic gene expression studies further supported these results. We also showed the contribution of topologically associated domain maps to WGS data interpretation.

CONCLUSION

Paired-end WGS is a valid strategy and may be used for structural variation characterisation in a clinical setting.

Comprehensive structural variation genome map of individuals carrying complex chromosomal rearrangements

Complex chromosomal rearrangements (CCRs) are rearrangements involving more than two chromosomes or more than two breakpoints. Whole genome sequencing (WGS) allows for outstanding high resolution characterization on the nucleotide level in unique sequences of such rearrangements, but problems remain for mapping breakpoints in repetitive regions of the genome, which are known to be prone to rearrangements. Hence, multiple complementary WGS experiments are sometimes needed to solve the structures of CCRs. We have studied three individuals with CCRs: Case 1 and Case 2 presented with de novo karyotypically balanced, complex interchromosomal rearrangements (46,XX,t(2;8;15)(q35;q24.1;q22) and 46,XY,t(1;10;5)(q32;p12;q31)), and Case 3 presented with a de novo, extremely complex intrachromosomal rearrangement on chromosome 1. Molecular cytogenetic investigation revealed cryptic deletions in the breakpoints of chromosome 2 and 8 in Case 1, and on chromosome 10 in Case 2, explaining their clinical symptoms. In Case 3, 26 breakpoints were identified using WGS, disrupting five known disease genes. All rearrangements were subsequently analyzed using optical maps, linked-read WGS, and short-read WGS. In conclusion, we present a case series of three unique de novo CCRs where we by combining the results from the different technologies fully solved the structure of each rearrangement. The power in combining short-read WGS with long-molecule sequencing or optical mapping in these unique de novo CCRs in a clinical setting is demonstrated.

Breakpoint mapping at nucleotide resolution in X-autosome balanced translocations associated with clinical phenotypes

Precise breakpoint mapping of balanced chromosomal rearrangements is crucial to identify disease etiology. Ten female patients with X-autosome balanced translocations associated with phenotypic alterations were evaluated, by mapping and sequencing their breakpoints. The rearrangements' impact on the expression of disrupted genes, and inferred mechanisms of formation in each case were assessed. For four patients that presented one of the chromosomal breaks in heterochromatic and highly repetitive segments, we combined cytogenomic methods and short-read sequencing to characterize, at nucleotide resolution, breakpoints that occurred in reference genome gaps. Most of rearrangements were possibly formed by non-homologous end joining and have breakpoints at repeat elements. Seven genes were found to be disrupted in six patients. Six of the affected genes showed altered expression, and the functional impairment of three of them were considered pathogenic. One gene disruption was considered potentially pathogenic, and three had uncertain clinical significance. Four patients presented no gene disruptions, suggesting other pathogenic mechanisms. Four genes were considered potentially affected by position effect and the expression abrogation of one of them was confirmed. This study emphasizes the importance of breakpoint-junction characterization at nucleotide resolution in balanced rearrangements to reveal genetic mechanisms associated with the patients' phenotypes, mechanisms of formation that originated the rearrangements, and genomic nature of disrupted DNA sequences.

Mapping of breakpoints in balanced chromosomal translocations by shallow whole-genome sequencing points to EFNA5, BAHD1 and PPP2R5E as novel candidates for genes causing human Mendelian disorders

BACKGROUND

Mapping the breakpoints in de novo balanced chromosomal translocations (BCT) in symptomatic individuals provides a unique opportunity to identify in an unbiased way the likely causative genetic defect and thus find novel human disease candidate genes. Our aim was to fine-map breakpoints of de novo BCTs in a case series of nine patients.

METHODS

Shallow whole-genome mate pair sequencing (SGMPS) together with long-range PCR and Sanger sequencing. In one case (BCT disrupting BAHD1 and RET) cDNA analysis was used to verify expression of a fusion transcript in cultured fibroblasts.

RESULTS

In all nine probands 11 disrupted genes were found, that is, EFNA5, EBF3, LARGE, PPP2R5E, TXNDC5, ZNF423, NIPBL, BAHD1, RET, TRPS1 and SLC4A10. Five subjects had translocations that disrupted genes with so far unknown (EFNA5, BAHD1, PPP2R5E, TXNDC5) or poorly delineated impact on the phenotype (SLC4A10, two previous reports of BCT disrupting the gene). The four genes with no previous disease associations (EFNA5, BAHD1, PPP2R5E, TXNDC5), when compared with all human genes by a bootstrap test, had significantly higher pLI (p<0.017) and DOMINO (p<0.02) scores indicating enrichment in genes likely to be intolerant to single copy damage. Inspection of individual pLI and DOMINO scores, and local topologically associating domain structure suggested that EFNA5, BAHD1 and PPP2R5E were particularly good candidates for novel disease loci. The pathomechanism for BAHD1 may involve deregulation of expression due to fusion with RET promoter.

CONCLUSION

SGMPS in symptomatic carriers of BCTs is a powerful approach to delineate novel human gene-disease associations.

Epigenetic maintenance of topological domains in the highly rearranged gibbon genome

The relationship between evolutionary genome remodeling and the three-dimensional structure of the genome remain largely unexplored. Here, we use the heavily rearranged gibbon genome to examine how evolutionary chromosomal rearrangements impact genome-wide chromatin interactions, topologically associating domains (TADs), and their epigenetic landscape. We use high-resolution maps of gibbon–human breaks of synteny (BOS), apply Hi-C in gibbon, measure an array of epigenetic features, and perform cross-species comparisons. We find that gibbon rearrangements occur at TAD boundaries, independent of the parameters used to identify TADs. This overlap is supported by a remarkable genetic and epigenetic similarity between BOS and TAD boundaries, namely presence of CpG islands and SINE elements, and enrichment in CTCF and H3K4me3 binding. Cross-species comparisons reveal that regions orthologous to BOS also correspond with boundaries of large (400–600 kb) TADs in human and other mammalian species. The colocalization of rearrangement breakpoints and TAD boundaries may be due to higher chromatin fragility at these locations and/or increased selective pressure against rearrangements that disrupt TAD integrity. We also examine the small portion of BOS that did not overlap with TAD boundaries and gave rise to novel TADs in the gibbon genome. We postulate that these new TADs generally lack deleterious consequences. Last, we show that limited epigenetic homogenization occurs across breakpoints, irrespective of their time of occurrence in the gibbon lineage. Overall, our findings demonstrate remarkable conservation of chromatin interactions and epigenetic landscape in gibbons, in spite of extensive genomic shuffling.

High-resolution detection of chromosomal rearrangements in leukemias through mate pair whole genome sequencing

The detection of recurrent somatic chromosomal rearrangements is standard of care for most leukemia types. Even though karyotype analysis-a low-resolution genome-wide chromosome analysis-is still the gold standard, it often needs to be complemented with other methods to increase resolution. To evaluate the feasibility and applicability of mate pair whole genome sequencing (MP-WGS) to detect structural chromosomal rearrangements in the diagnostic setting, we sequenced ten bone marrow samples from leukemia patients with recurrent rearrangements. Samples were selected based on cytogenetic and FISH results at leukemia diagnosis to include common rearrangements of prognostic relevance. Using MP-WGS and in-house bioinformatic analysis all sought rearrangements were successfully detected. In addition, unexpected complexity or additional, previously undetected rearrangements was unraveled in three samples. Finally, the MP-WGS analysis pinpointed the location of chromosome junctions at high resolution and we were able to identify the exact exons involved in the resulting fusion genes in all samples and the specific junction at the nucleotide level in half of the samples. The results show that our approach combines the screening character from karyotype analysis with the specificity and resolution of cytogenetic and molecular methods. As a result of the straightforward analysis and high-resolution detection of clinically relevant rearrangements, we conclude that MP-WGS is a feasible method for routine leukemia diagnostics of structural chromosomal rearrangements.

Balanced Chromosomal Rearrangement Detection by Low-Pass Whole-Genome Sequencing

Balanced chromosomal rearrangements (or balanced chromosome abnormalities, BCAs) are common chromosomal structural variants. Emerging studies have demonstrated the feasibility of using whole-genome sequencing (WGS) for detection of BCA-associated breakpoints, but the requirement for a priori knowledge of the rearranged regions from G-banded chromosome analysis limits its application. The protocols described here are based on low-pass WGS for detecting BCA events independent from chromosome analysis, and has been validated using genomic data from the 1000 Genomes Project. This approach adopts non-size-selected mate-pair library (3∼8 kb) with 2∼3 μg DNA as input, and requires only 30 million read-pairs (50 bp, equivalent to 1-fold base-coverage) for each sample. The complete procedure takes 13 days and the total cost is estimated to be less than $600 (USD) per sample. © 2018 by John Wiley & Sons, Inc.

Phenotypic interpretation of complex chromosomal rearrangements informed by nucleotide-level resolution and structural organization of chromatin

Molecular characterization of balanced chromosomal abnormalities constitutes a powerful tool in understanding the pathogenic mechanisms of complex genetic disorders. Here we report a male with severe global developmental delay in the presence of a complex karyotype and normal microarray and exome studies. The subject, referred to as DGAP294, has two de novo apparently balanced translocations involving chromosomes 1 and 14, and chromosomes 4 and 10, disrupting several different transcripts of adhesion G protein-coupled receptor L2 (ADGRL2) and protocadherin 15 (PCDH15). In addition, a maternally inherited inversion disrupts peptidyl arginine deiminase types 3 and 4 (PADI3 and PADI4) on chromosome 1. None of these gene disruptions explain the patient's phenotype. Using genome regulatory annotations and chromosome conformation data, we predict a position effect ~370 kb upstream of a translocation breakpoint located at 14q12. The position effect involves forkhead box G1 (FOXG1), mutations in which are associated with the congenital form of Rett syndrome and FOXG1 syndrome. We believe the FOXG1 position effect largely accounts for the clinical phenotype in DGAP294, which can be classified as FOXG1 syndrome like. Our findings emphasize the significance of not only analyzing disrupted genes by chromosomal rearrangements, but also evaluating potential long-range position effects in clinical diagnoses.

Historical and Clinical Perspectives on Chromosomal Translocations

Chromosomal translocations, rearrangements involving the exchange of segments between chromosomes, were documented in humans in 1959. The first accurately reported clinical phenotype resulting from a translocation was that of Down syndrome. In a small percentage of Down syndrome cases, an extra 21q is provided by a Robertsonian translocation chromosome, either occurring de novo or inherited from a phenotypically normal parent with the translocation chromosome and a balanced genome of 45 chromosomes. Balanced translocations, including both Robertsonian and reciprocal translocations, are typically benign, but meiosis in germ cells with balanced translocations may result in meiotic arrest and subsequent infertility, or in unbalanced gametes, with attendant risks of miscarriage and unbalanced progeny. Most reciprocal translocations are unique. A few to several percent of translocations disrupt haploinsufficient genes or their regulatory regions and result in clinical phenotypes. Balanced translocations from patients with clinical phenotypes have been valuable in mapping disease genes and in illuminating cis-regulatory regions. Mapping of discordant mate pairs from long-insert, low-pass genome sequencing now permits efficient and cost-effective discovery and nucleotide-level resolution of rearrangement breakpoints, information that is absolutely necessary for interpreting the etiology of clinical phenotypes in patients with rearrangements. Pathogenic translocations and other balanced chromosomal rearrangements constitute a class of typically highly penetrant mutation that is cryptic to both clinical microarray and exome sequencing. A significant proportion of rearrangements include additional complexity that is not visible by conventional karyotype analysis. Some proportion of patients with negative findings on exome/genome sequencing and clinical microarray will be found to have etiologic balanced rearrangements only discoverable by genome sequencing with analysis pipelines optimized to recover rearrangement breakpoints.

Induction of Chromosomal Translocations with CRISPR-Cas9 and Other Nucleases: Understanding the Repair Mechanisms That Give Rise to Translocations

Chromosomal translocations are associated with several tumor types, including hematopoietic malignancies, sarcomas, and solid tumors of epithelial origin, due to their activation of a proto-oncogene or generation of a novel fusion protein with oncogenic potential. In many cases, the availability of suitable human models has been lacking because of the difficulty in recapitulating precise expression of the fusion protein or other reasons. Further, understanding how translocations form mechanistically has been a goal, as it may suggest ways to prevent their occurrence. Chromosomal translocations arise when DNA ends from double-strand breaks (DSBs) on two heterologous chromosomes are improperly joined. This review provides a summary of DSB repair mechanisms and their contribution to translocation formation, the various programmable nuclease platforms that have been used to generate translocations, and the successes that have been achieved in this area.

Identification of balanced chromosomal rearrangements previously unknown among participants in the 1000 Genomes Project: implications for interpretation of structural variation in genomes and the future of clinical cytogenetics

Purpose

Recent studies demonstrate that whole-genome sequencing (WGS) enables detection of cryptic rearrangements in apparently balanced chromosomal rearrangements (also known as balanced chromosomal abnormalities, BCAs) previously identified by conventional cytogenetic methods. We aimed to assess our analytical tool for detecting BCAs in The 1000 Genomes Project without knowing affected bands.

Methods

The 1000 Genomes Project provides an unprecedented integrated map of structural variants in phenotypically normal subjects, but there is no information on potential inclusion of subjects with apparently BCAs akin to those traditionally detected in diagnostic cytogenetics laboratories. We applied our analytical tool to 1,166 genomes from the 1000 Genomes Project with sufficient physical coverage (8.25-fold).

Results

Our approach detected four reciprocal balanced translocations and four inversions ranging in size from 57.9 kb to 13.3 Mb, all of which were confirmed by cytogenetic methods and PCR studies. One of DNAs has a subtle translocation that is not readily identified by chromosome analysis due to similar banding patterns and size of exchanged segments, and another results in disruption of all transcripts of an OMIM gene.

Conclusions

Our study demonstrates the extension of utilizing low-coverage WGS for unbiased detection of BCAs including translocations and inversions previously unknown in the 1000 Genomes Project.