Computational Prediction of Position Effects of Apparently Balanced Human Chromosomal Rearrangements

Complex balanced intrachromosomal rearrangement involving PITX2 identified as a cause of Axenfeld-Rieger Syndrome

Axenfeld-Rieger Syndrome (ARS) type 1 is a rare autosomal dominant condition characterized by anterior chamber anomalies, umbilical defects, dental hypoplasia, and craniofacial anomalies, with Meckel's diverticulum in some individuals. Here, we describe a clinically ascertained female of childbearing age with ARS for whom clinical targeted sequencing and deletion/duplication analysis followed by clinical exome and genome sequencing resulted in no pathogenic variants or variants of unknown significance in PITX2 or FOXC1. Advanced bioinformatic analysis of the genome data identified a complex, balanced rearrangement disrupting PITX2. This case is the first reported intrachromosomal rearrangement leading to ARS, illustrating that for patients with compelling clinical phenotypes but negative genomic testing, additional bioinformatic analysis are essential to identify subtle genomic abnormalities in target genes.

Balanced chromosomal rearrangements implicate YIPF5 and SPATC1L in non-obstructive oligoasthenozoospermia and oligozoospermia and of a derivative chromosome 22 in recurrent miscarriage

Naturally occurring balanced, unbalanced, and complex chromosomal rearrangements have been reported to cause pathogenic genomic or genetic variants leading to infertility and recurrent miscarriage. Therefore, balanced chromosomal rearrangements were used as genomic signposts for identification of candidate genes or genomic loci associated with male infertility due to defects of spermatogenesis, or with recurrent miscarriage. In three male probands, structural chromosomal variants and copy number variants were identified at nucleotide resolution by long-insert genome sequencing approaches and Sanger sequencing. The pathogenic potential of these and affected candidate genes was assessed based on convergent genomic and genotype-phenotype correlation data. Identification of balanced chromosomal rearrangement breakpoints and interpretation in the context of their genomic background of structural and copy number variants led us to conclude that the infertility due to oligoasthenozoospermia and oligozoospermia is most likely associated with a position effect on YIPF5 and SPATC1L, respectively. In a third proband with intellectual disability and recurrent miscarriage, disruption of CAMK2B causing autosomal dominant, intellectual developmental disorder 54 and increased meiotic segregation during gametogenesis of a der(22) are responsible for the reported phenotype. Our data further support the existence of loci at 5q23 and 21q22.3 for these spermatogenesis defects and highlight the importance of the naturally occurring balanced chromosomal rearrangements for assessment of the pathogenic mechanisms. Furthermore, we show comorbidities due to the same balanced chromosomal rearrangement caused by different pathogenic mechanisms.

Disruption of regulatory domains and novel transcripts as disease-causing mechanisms

Deletions, duplications, insertions, inversions, and translocations, collectively called structural variations (SVs), affect more base pairs of the genome than any other sequence variant. The recent technological advancements in genome sequencing have enabled the discovery of tens of thousands of SVs per human genome. These SVs primarily affect non-coding DNA sequences, but the difficulties in interpreting their impact limit our understanding of human disease etiology. The functional annotation of non-coding DNA sequences and methodologies to characterize their three-dimensional (3D) organization in the nucleus have greatly expanded our understanding of the basic mechanisms underlying gene regulation, thereby improving the interpretation of SVs for their pathogenic impact. Here, we discuss the various mechanisms by which SVs can result in altered gene regulation and how these mechanisms can result in rare genetic disorders. Beyond changing gene expression, SVs can produce novel gene-intergenic fusion transcripts at the SV breakpoints.

Premature ovarian insufficiency is associated with global alterations in the regulatory landscape and gene expression in balanced X-autosome translocations

BACKGROUND

Patients with balanced X-autosome translocations and premature ovarian insufficiency (POI) constitute an interesting paradigm to study the effect of chromosome repositioning. Their breakpoints are clustered within cytobands Xq13-Xq21, 80% of them in Xq21, and usually, no gene disruption can be associated with POI phenotype. As deletions within Xq21 do not cause POI, and since different breakpoints and translocations with different autosomes lead to this same gonadal phenotype, a "position effect" is hypothesized as a possible mechanism underlying POI pathogenesis.

OBJECTIVE AND METHODS

To study the effect of the balanced X-autosome translocations that result in POI, we fine-mapped the breakpoints in six patients with POI and balanced X-autosome translocations and addressed gene expression and chromatin accessibility changes in four of them.

RESULTS

We observed differential expression in 85 coding genes, associated with protein regulation, multicellular regulation, integrin signaling, and immune response pathways, and 120 differential peaks for the three interrogated histone marks, most of which were mapped in high-activity chromatin state regions. The integrative analysis between transcriptome and chromatin data pointed to 12 peaks mapped less than 2 Mb from 11 differentially expressed genes in genomic regions not related to the patients' chromosomal rearrangement, suggesting that translocations have broad effects on the chromatin structure.

CONCLUSION

Since a wide impact on gene regulation was observed in patients, our results observed in this study support the hypothesis of position effect as a pathogenic mechanism for premature ovarian insufficiency associated with X-autosome translocations. This work emphasizes the relevance of chromatin changes in structural variation, since it advances our knowledge of the impact of perturbations in the regulatory landscape within interphase nuclei, resulting in the position effect pathogenicity.

TADeus2: a web server facilitating the clinical diagnosis by pathogenicity assessment of structural variations disarranging 3D chromatin structure

In recent years great progress has been made in identification of structural variants (SV) in the human genome. However, the interpretation of SVs, especially located in non-coding DNA, remains challenging. One of the reasons stems in the lack of tools exclusively designed for clinical SVs evaluation acknowledging the 3D chromatin architecture. Therefore, we present TADeus2 a web server dedicated for a quick investigation of chromatin conformation changes, providing a visual framework for the interpretation of SVs affecting topologically associating domains (TADs). This tool provides a convenient visual inspection of SVs, both in a continuous genome view as well as from a rearrangement’s breakpoint perspective. Additionally, TADeus2 allows the user to assess the influence of analyzed SVs within flaking coding/non-coding regions based on the Hi-C matrix. Importantly, the SVs pathogenicity is quantified and ranked using TADA, ClassifyCNV tools and sampling-based P-value. TADeus2 is publicly available at https://tadeus2.mimuw.edu.pl.

Multidisciplinary approaches for elucidating genetics and molecular pathogenesis of urinary tract malformations

Advances in clinical diagnostics and molecular tools have improved our understanding of the genetically heterogeneous causes underlying congenital anomalies of kidney and urinary tract (CAKUT). However, despite a sharp incline of CAKUT reports in the literature within the past 2 decades, there remains a plateau in the genetic diagnostic yield that is disproportionate to the accelerated ability to generate robust genome-wide data. Explanations for this observation include (i) diverse inheritance patterns with incomplete penetrance and variable expressivity, (ii) rarity of single-gene drivers such that large sample sizes are required to meet the burden of proof, and (iii) multigene interactions that might produce either intra- (e.g., copy number variants) or inter- (e.g., effects in trans) locus effects. These challenges present an opportunity for the community to implement innovative genetic and molecular avenues to explain the missing heritability and to better elucidate the mechanisms that underscore CAKUT. Here, we review recent multidisciplinary approaches at the intersection of genetics, genomics, in vivo modeling, and in vitro systems toward refining a blueprint for overcoming the diagnostic hurdles that are pervasive in urinary tract malformation cohorts. These approaches will not only benefit clinical management by reducing age at molecular diagnosis and prompting early evaluation for comorbid features but will also serve as a springboard for therapeutic development.

Factors affecting clinical manifestation of chromosomal imbalance in carriers of segmental autosomal mosaicism: differential impact of gender

Mosaicism for unbalanced chromosomal rearrangements segmental mosaicism (SM) is rare, both in patients referred for cytogenetic testing and in prenatal diagnoses. In contrast, in preimplantation embryos SM is a frequent finding and, therefore, is even more challenging. However, there is no consistency among results of published studies on the clinical outcomes of embryos with SM, primarily due to the small number of reported cases. Moreover, there is the problem of predicting the potential for the optimal development of a mosaic embryo to a healthy individual. Therefore, we suggested comparing factors predisposing to favorable and poor prognoses, identified in postnatal and prenatal cohorts of SM carriers, with those obtained from studies on preimplantation embryos. We analyzed 580 published cases of SM including (i) postnatally diagnosed affected carriers, (ii) clinically asymptomatic carriers, (iii) prenatally diagnosed carriers, and (iv) miscarriages. We observed a concordance with preimplantation diagnoses regarding the clinical significance of the extent of mosaicism as well as a predominance of deletions over other types of rearrangements. However, there is no concordance regarding excessive involvement of chromosomes 1, 5, and 9 in unbalanced rearrangements and a preferential involvement of larger chromosomes compared to short ones. Paternal age was not found to be associated with SM in postnatally disease-defined individuals. We have identified maternal age and preferential involvement of chromosome 18 in rearrangements associated with clinical manifestations. Male predominance was found among normal pregnancy outcomes and among disease-defined carriers of rearrangements resulting in a gain of genomic material. Female predominance was found among abnormal pregnancy outcomes, among disease-defined carriers of loss and gain/loss rearrangements, and among transmitting carriers of gonadal SM, both affected and asymptomatic. According to data obtained from “post-embryo” studies, clinical manifestations of chromosomal imbalance are associated with a high proportion of abnormal cells, female gender, the type of rearrangement and involved chromosome(s), and maternal age. We believe these data are instructive in the challenging medical genetic counseling of parents faced with no option other than transfer of an embryo with segmental mosaicism.

SVInterpreter: A Comprehensive Topologically Associated Domain-Based Clinical Outcome Prediction Tool for Balanced and Unbalanced Structural Variants

With the advent of genomic sequencing, a number of balanced and unbalanced structural variants (SVs) can be detected per individual. Mainly due to incompleteness and the scattered nature of the available annotation data of the human genome, manual interpretation of the SV's clinical significance is laborious and cumbersome. Since bioinformatic tools developed for this task are limited, a comprehensive tool to assist clinical outcome prediction of SVs is warranted. Herein, we present SVInterpreter, a free Web application, which analyzes both balanced and unbalanced SVs using topologically associated domains (TADs) as genome units. Among others, gene-associated data (as function and dosage sensitivity), phenotype similarity scores, and copy number variants (CNVs) scoring metrics are retrieved for an informed SV interpretation. For evaluation, we retrospectively applied SVInterpreter to 97 balanced (translocations and inversions) and 125 unbalanced (deletions, duplications, and insertions) previously published SVs, and 145 SVs identified from 20 clinical samples. Our results showed the ability of SVInterpreter to support the evaluation of SVs by (1) confirming more than half of the predictions of the original studies, (2) decreasing 40% of the variants of uncertain significance, and (3) indicating several potential position effect events. To our knowledge, SVInterpreter is the most comprehensive TAD-based tool to identify the possible disease-causing candidate genes and to assist prediction of the clinical outcome of SVs. SVInterpreter is available at http://dgrctools-insa.min-saude.pt/cgi-bin/SVInterpreter.py.

Position effects of 22q13 rearrangements on candidate genes in Phelan-McDermid syndrome

Phelan-McDermid syndrome (PMS) is a multi-system disorder characterized by significant variability in clinical presentation. The genetic etiology is also variable with differing sizes of deletions in the chromosome 22q13 region and types of genetic abnormalities (e.g., terminal or interstitial deletions, translocations, ring chromosomes, or SHANK3 variants). Position effects have been shown to affect gene expression and function and play a role in the clinical presentation of various genetic conditions. This study employed a topologically associating domain (TAD) analysis approach to investigate position effects of chromosomal rearrangements on selected candidate genes mapped to 22q13 in 81 individuals with PMS. Data collected were correlated with clinical information from these individuals and with expression and metabolic profiles of lymphoblastoid cells from selected cases. The data confirmed TAD predictions for genes encompassed in the deletions and the clinical and molecular data indicated clear differences among individuals with different 22q13 deletion sizes. The results of the study indicate a positive correlation between deletion size and phenotype severity in PMS and provide evidence of the contribution of other genes to the clinical variability in this developmental disorder by reduced gene expression and altered metabolomics.

Congenital chromoanagenesis in the routine postnatal chromosomal microarray analyses

Chromosomal microarray analyses (CMA) have greatly increased both the yield and diagnostic accuracy of postnatal analysis; it has been used as a first-tier cytogenetic test in patients with intellectual disability, autism spectrum disorder, and multiple congenital abnormalities. During the last 15 years, we performed CMA in approximately 8,000 patients with neurodevelopmental and/or congenital disorders, of which 13 (0.16%) genetically catastrophic complex chromosomal rearrangements were identified. These ultrarare rearrangements showed clustering of breakpoints, characteristic of chromoanagenesis events. Al1 13 complex events display underlying formation mechanisms, originating either by a synchronization of the shattering of clustered chromosome regions in which regional asynchrony of DNA replication may be one of the main causes of disruption. We provide an overview of the copy number profiling in these patients. Although several previous studies have suggested that chromoanagenesis is often a genetic disease source in postnatal diagnostic screening, due to either the challenge of clinical interpretation of these complex rearrangements or the limitation of microarray resolution relative to the small size and complexity of chromogenic induced chromosome abnormalities, bringing further attention and to study its occurrence in the clinical setting is extremely important.

What a finding of gene copy number variation can add to the diagnosis of developmental neuropsychiatric disorders

Among medical disciplines, diagnosis in psychiatry depends highly upon descriptive signs and symptoms, rather than biomarkers. Clear descriptions of specific genetic etiologies have been lacking; genomic technologies, however, are rapidly changing that landscape. Notably, chromosomal microarrays-which detect gene copy number variants (CNVs)-are a recommended standard of care for neurodevelopmental disorders. As a result, an increasing number of patients now receive a clinical diagnosis based on the Diagnostic and Statistical Manual of Mental Disorders (DSM-5) and an identified genetic etiological variant. However, psychiatric and genetic diagnoses are frequently communicated and managed as two disconnected diagnostic parameters. Here, we advocate for a transition model, allowing the integration of genetic etiological information-starting with diagnostically proven CNVs-within the DSM-5 classification framework.

Limited diagnostic impact of duplications <1 Mb of uncertain clinical significance: a 10-year retrospective analysis of reporting practices at the Mayo Clinic

PURPOSE

Copy-number variants (CNVs) of uncertain clinical significance are routinely reported in a clinical setting only when exceeding predetermined reporting thresholds, typically based on CNV size. Given that very few genes are associated with triplosensitive phenotypes, it is not surprising that many interstitial duplications <1 Mb are found to be inherited and anticipated to be of limited or no clinical significance.

METHODS

In an effort to further refine our reporting criteria to maximize diagnostic yield while minimizing the return of uncertain variants, we performed a retrospective analysis of all clinical microarray cases reported in a 10-year window. A total of 1112 reported duplications had parental follow-up, and these were compared by size, RefSeq gene content, and inheritance pattern. De novo origin was used as a rough proxy for pathogenicity.

RESULTS

Approximately 6% of duplications 500 kb-1 Mb were de novo observations, compared with approximately 14% for 1-2 Mb duplications (p = 0.0005). On average, de novo duplications had higher gene counts than inherited duplications.

CONCLUSION

Our data reveal limited diagnostic utility for duplications of uncertain significance <1 Mb. Considerations for revised reporting criteria are discussed and are applicable to CNVs detected by any genome-wide exploratory methodology, including exome/genome sequencing.

Rare or Overlooked? Structural Disruption of Regulatory Domains in Human Neurocristopathies

In the last few years, the role of non-coding regulatory elements and their involvement in human disease have received great attention. Among the non-coding regulatory sequences, enhancers are particularly important for the proper establishment of cell type-specific gene-expression programs. Furthermore, the disruption of enhancers can lead to human disease through two main mechanisms: (i) Mutations or copy number variants can directly alter the enhancer sequences and thereby affect expression of their target genes; (ii) structural variants can provoke changes in 3-D chromatin organization that alter neither the enhancers nor their target genes, but rather the physical communication between them. In this review, these pathomechanisms are mostly discussed in the context of neurocristopathies, congenital disorders caused by defects that occur during neural crest development. We highlight why, due to its contribution to multiple tissues and organs, the neural crest represents an important, yet understudied, cell type involved in multiple congenital disorders. Moreover, we discuss currently available resources and experimental models for the study of human neurocristopathies. Last, we provide some practical guidelines that can be followed when investigating human neurocristopathies caused by structural variants. Importantly, these guidelines can be useful not only to uncover the etiology of human neurocristopathies, but also of other human congenital disorders in which enhancer disruption is involved.

Genotoxic potential of a novel PDE-4B inhibitor Apremilast by chromosomal aberration and micronucleus assay in mice

Objective

Researchers have confirmed that chronic administration of drugs at high doses causes genotoxicity which serve as first step in development of cancers. Apremilast, a phosphodiesterase-4 inhibitor is Food and Drug Administration (FDA) approved drug for Psoriatic Arthritis. The present study designed to conduct genotoxicity testing using the genotoxic study which give simple, sensitive, economical and fast tools for the assessment of damage of genetic material.

Methods

To conduct genotoxicity study of Apremilast, 60 Swiss albino male mice divided into 6 groups (n = 10). Group1 served as a normal control group without any treatment, Group 2 treated as a disease control and administered with cyclophosphamide 40 mg/kg, IP. Group 3, 4, 5 and 6 treated as test groups and received 10, 20, 40 and 80 mg/kg/day Apremilast respectively. The total duration of study was 13 weeks. At termination day animals were sacrificed and chromosomal aberration assay (BMCAA) and micronucleus assay (BMMNA) were performed to know the genotoxicity potential of Apremilast.

Results

The results indicates significant rise in chromosomal aberrations (CA) frequency in bone marrow cells and decrease in the MI of the disease control animals as well as Apremilast treated groups. Further significant (p < 0.001; p < 0.0001) increase in score of micronucleated polychromatic erythrocytes (MNPCEs) and percentage of micronucleated PCEs per 1000 PCEs and decrease in the ratio of polychromatic/normochromatic erythrocytes (PCE/NCE) was observed in micronucleus assay. Genotoxic effect increases with the increase of Apremilast dose. Conclusion: Finding of present indicates that Apremilast shows genotoxic potential on high administration although further detailed toxicity studies required for confirmations.

Low-pass genome sequencing versus chromosomal microarray analysis: implementation in prenatal diagnosis

PURPOSE

Emerging studies suggest that low-pass genome sequencing (GS) provides additional diagnostic yield of clinically significant copy-number variants (CNVs) compared with chromosomal microarray analysis (CMA). However, a prospective back-to-back comparison evaluating accuracy, efficacy, and incremental yield of low-pass GS compared with CMA is warranted.

METHODS

A total of 1023 women undergoing prenatal diagnosis were enrolled. Each sample was subjected to low-pass GS and CMA for CNV analysis in parallel. CNVs were classified according to guidelines of the American College of Medical Genetics and Genomics.

RESULTS

Low-pass GS not only identified all 124 numerical disorders or pathogenic or likely pathogenic (P/LP) CNVs detected by CMA in 121 cases (11.8%, 121/1023), but also defined 17 additional and clinically relevant P/LP CNVs in 17 cases (1.7%, 17/1023). In addition, low-pass GS significantly reduced the technical repeat rate from 4.6% (47/1023) for CMA to 0.5% (5/1023) and required less DNA (50 ng) as input.

CONCLUSION

In the context of prenatal diagnosis, low-pass GS identified additional and clinically significant information with enhanced resolution and increased sensitivity of detecting mosaicism as compared with the CMA platform used. This study provides strong evidence for applying low-pass GS as an alternative prenatal diagnostic test.

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.

Quantitative prediction of enhancer-promoter interactions

Recent experimental and computational efforts have provided large data sets describing three-dimensional organization of mouse and human genomes and showed the interconnection between the expression profile, epigenetic state, and spatial interactions of loci. These interconnections were utilized to infer the spatial organization of chromatin, including enhancer–promoter contacts, from one-dimensional epigenetic marks. Here, we show that the predictive power of some of these algorithms is overestimated due to peculiar properties of the biological data. We propose an alternative approach, which provides high-quality predictions of chromatin interactions using information on gene expression and CTCF-binding alone. Using multiple metrics, we confirmed that our algorithm could efficiently predict the three-dimensional architecture of both normal and rearranged genomes.

A Retrospective Analysis of 10-Year Data Assessed the Diagnostic Accuracy and Efficacy of Cytogenomic Abnormalities in Current Prenatal and Pediatric Settings

Background: Array comparative genomic hybridization (aCGH), karyotyping and fluorescence in situ hybridization (FISH) analyses have been used in a clinical cytogenetic laboratory. A systematic analysis on diagnostic findings of cytogenomic abnormalities in current prenatal and pediatric settings provides approaches for future improvement.

Methods: A retrospective analysis was performed on abnormal findings by aCGH, karyotyping, and FISH from 3,608 prenatal cases and 4,509 pediatric cases during 2008–2017. The diagnostic accuracy was evaluated by comparing the abnormality detection rate (ADR) and the relative frequency (RF) of different types of cytogenomic abnormalities between prenatal and pediatric cases. A linear regression correlation between known prevalence and ADR of genomic disorders was used to extrapolate the prevalence of other genomic disorders. The diagnostic efficacy was estimated as percentage of detected abnormal cases by expected abnormal cases from served population.

Results: The composite ADR for numerical chromosome abnormalities, structural chromosome abnormalities, recurrent genomic disorders, and sporadic pathogenic copy number variants (pCNVs) in prenatal cases were 13.03%, 1.77%, 1.69%, and 0.9%, respectively, and were 5.13%, 2.84%, 7.08%, and 2.69% in pediatric cases, respectively. The chromosomal abnormalities detected in prenatal cases (14.80%) were significantly higher than that of pediatric cases (7.97%) (p < 0.05), while the pCNVs detected in prenatal cases (2.59%) were significantly lower than that of pediatric cases (9.77%) (p < 0.05). The prevalence of recurrent genomic disorders and total pCNVs was estimated to be 1/396 and 1/291, respectively. Approximately, 29% and 35% of cytogenomic abnormalities expected from the population served were detected in current prenatal and pediatric diagnostic practice, respectively.

Conclusion: For chromosomal abnormalities, effective detection of Down syndrome (DS) and Turner syndrome (TS) and under detection of sex chromosome numerical abnormalities in both prenatal and pediatric cases were noted. For pCNVs, under detection of pCNVs in prenatal cases and effective detection of DiGeorge syndrome (DGS) and variable efficacy in detecting other pCNVs in pediatric cases were noted. Extend aCGH analysis to more prenatal cases with fetal ultrasonographic anomalies, enhanced non-invasive prenatal (NIPT) testing screening for syndromic genomic disorders, and better clinical indications for pCNVs are approaches that could improve diagnostic yield of cytogenomic abnormalities.

Structural variation in the 3D genome

Structural and quantitative chromosomal rearrangements, collectively referred to as structural variation (SV), contribute to a large extent to the genetic diversity of the human genome and thus are of high relevance for cancer genetics, rare diseases and evolutionary genetics. Recent studies have shown that SVs can not only affect gene dosage but also modulate basic mechanisms of gene regulation. SVs can alter the copy number of regulatory elements or modify the 3D genome by disrupting higher-order chromatin organization such as topologically associating domains. As a result of these position effects, SVs can influence the expression of genes distant from the SV breakpoints, thereby causing disease. The impact of SVs on the 3D genome and on gene expression regulation has to be considered when interpreting the pathogenic potential of these variant types.

Case for genome sequencing in infants and children with rare, undiagnosed or genetic diseases

Up to 350 million people worldwide suffer from a rare disease, and while the individual diseases are rare, in aggregate they represent a substantial challenge to global health systems. The majority of rare disorders are genetic in origin, with children under the age of five disproportionately affected. As these conditions are difficult to identify clinically, genetic and genomic testing have become the backbone of diagnostic testing in this population. In the last 10 years, next-generation sequencing technologies have enabled testing of multiple disease genes simultaneously, ranging from targeted gene panels to exome sequencing (ES) and genome sequencing (GS). GS is quickly becoming a practical first-tier test, as cost decreases and performance improves. A growing number of studies demonstrate that GS can detect an unparalleled range of pathogenic abnormalities in a single laboratory workflow. GS has the potential to deliver unbiased, rapid and accurate molecular diagnoses to patients across diverse clinical indications and complex presentations. In this paper, we discuss clinical indications for testing and historical testing paradigms. Evidence supporting GS as a diagnostic tool is supported by superior genomic coverage, types of pathogenic variants detected, simpler laboratory workflow enabling shorter turnaround times, diagnostic and reanalysis yield, and impact on healthcare.

The Iceberg under Water: Unexplored Complexity of Chromoanagenesis in Congenital Disorders

Structural variation, composed of balanced and unbalanced genomic rearrangements, is an important contributor to human genetic diversity with prominent roles in somatic and congenital disease. At the nucleotide level, structural variants (SVs) have been shown to frequently harbor additional breakpoints and copy-number imbalances, a complexity predicted to emerge wholly as a single-cell division event. Chromothripsis, chromoplexy, and chromoanasynthesis, collectively referred to as chromoanagenesis, are three major mechanisms that explain the occurrence of complex germline and somatic SVs. While chromothripsis and chromoplexy have been shown to be key signatures of cancer, chromoanagenesis has been detected in numerous cases of developmental disease and phenotypically normal individuals. Such observations advocate for a deeper study of the polymorphic and pathogenic properties of complex germline SVs, many of which go undetected by traditional clinical molecular and cytogenetic methods. This review focuses on congenital chromoanagenesis, mechanisms leading to occurrence of these complex rearrangements, and their impact on chromosome organization and genome function. We highlight future applications of routine screening of complex and balanced SVs in the clinic, as these represent a potential and often neglected genetic disease source, a true "iceberg under water."

Chromosomal Microarray Analysis versus Karyotyping in Fetuses with Increased Nuchal Translucency

We have carried out a retrospective study of chromosome anomalies associated with increased nuchal translucency (NT) in order to compare yield rates of karyotype, chromosome microarray analysis (CMA), and non-invasive prenatal testing (NIPT) in this condition. Presenting with increased NT or cystic hygroma ≥3.5 mm as an isolated sign, 249 fetuses underwent karyotype and/or CMA from 11 to 18 gestational weeks. Karyotype and fluorescence in situ hybridization (FISH) analyses detected 103 chromosomal anomalies including 95 aneuploidies and eight chromosomal rearrangements or derivatives. Further, seven pathogenic copy number variants (CNV), five likely pathogenic CNVs, and 15 variants of unknown significance (VOUS) were detected by CMA in fetuses with normal karyotype. Genetic testing is now facing new challenges due to results with uncertain clinical impacts. Additional investigations will be necessary to interpret these findings. More than 15% of the anomalies that we have diagnosed with invasive techniques could not be detected by NIPT. It is therefore definitely not recommended in the case of ultrasound anomalies. These results, while corroborating the use of CMA in fetuses with increased NT as a second tier after rapid aneuploidy testing, do not suggest a dismissal of karyotype analysis.

Chromatin features constrain structural variation across evolutionary timescales

Noncoding DNA sequences play crucial roles in gene regulation, including via three-dimensional genome organization where they define chromatin boundaries and segment the genome into a sequence of insulated neighborhoods. However, the relative importance of noncoding DNA elements, particularly in comparison with protein-coding DNA sequences, remains more poorly characterized. Here, we systematically test if chromatin boundary disruptions are under purifying selection. Our analyses uncover a genomewide depletion of structural variants that would have the potential to alter chromatin structure. This in turn has implications for predicting not only which variants are likely pathogenic in clinical genetics settings, but also which are likely key innovations in primate evolution, and argues for expanding the current gene-centric paradigm for interpreting structural variants.

The potential impact of structural variants includes not only the duplication or deletion of coding sequences, but also the perturbation of noncoding DNA regulatory elements and structural chromatin features, including topological domains (TADs). Structural variants disrupting TAD boundaries have been implicated both in cancer and developmental disease; this likely occurs via “enhancer hijacking,” whereby removal of the TAD boundary exposes enhancers to new target transcription start sites (TSSs). With this functional role, we hypothesized that boundaries would display evidence for negative selection. Here we demonstrate that the chromatin landscape constrains structural variation both within healthy humans and across primate evolution. In contrast, in patients with developmental delay, variants occur remarkably uniformly across genomic features, suggesting a potentially broad role for enhancer hijacking in human disease.

X-chromosomale Intelligenzminderung

X-chromosomale Intelligenzminderung („X-linked intellectual disability“, XLID) ist eine heterogene Krankheitsgruppe; inzwischen sind mehr als 100 XLID-Gene identifiziert worden. Das Fragile-X-Syndrom mit CGG-Repeatexpansion in der 5’-UTR des FMR1-Gens ist die häufigste monogene Ursache für Intelligenzminderung. Weitere X‑chromosomale Gene mit vergleichsweise hohen Mutationsprävalenzen sind ATRX, RPS6KA3, GPC3, SLC16A2, SLC6A8 und ARX. Die Ursachen für XLID verteilen sich zu ca. 90 % auf molekulargenetisch nachweisbare Mutationen und zu ca. 10 % auf chromosomale Kopienzahlvarianten („copy-number variants“, CNVs). Häufige CNVs sind Duplikationen in Xq28 unter Einschluss von MECP2 sowie das Xp11.22-Duplikations-Syndrom mit Überexpression von HUWE1. Mit den aktuellen Untersuchungsmethoden kann bei ca. 10 % der männlichen Patienten mit Intelligenzminderung eine X‑chromosomale Ursache nachgewiesen werden. Neue Erkenntnisse zu XLID sind für die nächsten Jahre am ehesten in den nicht kodierenden Regionen zu erwarten, wo wahrscheinlich ein weiterer Teil der Ursachen für das bislang nicht vollständig erklärte Überwiegen männlicher Patienten zu suchen ist.

Evolutionary stability of topologically associating domains is associated with conserved gene regulation

Background

The human genome is highly organized in the three-dimensional nucleus. Chromosomes fold locally into topologically associating domains (TADs) defined by increased intra-domain chromatin contacts. TADs contribute to gene regulation by restricting chromatin interactions of regulatory sequences, such as enhancers, with their target genes. Disruption of TADs can result in altered gene expression and is associated to genetic diseases and cancers. However, it is not clear to which extent TAD regions are conserved in evolution and whether disruption of TADs by evolutionary rearrangements can alter gene expression.

Results

Here, we hypothesize that TADs represent essential functional units of genomes, which are stable against rearrangements during evolution. We investigate this using whole-genome alignments to identify evolutionary rearrangement breakpoints of different vertebrate species. Rearrangement breakpoints are strongly enriched at TAD boundaries and depleted within TADs across species. Furthermore, using gene expression data across many tissues in mouse and human, we show that genes within TADs have more conserved expression patterns. Disruption of TADs by evolutionary rearrangements is associated with changes in gene expression profiles, consistent with a functional role of TADs in gene expression regulation.

Conclusions

Together, these results indicate that TADs are conserved building blocks of genomes with regulatory functions that are often reshuffled as a whole instead of being disrupted by rearrangements.

Electronic supplementary material

The online version of this article (10.1186/s12915-018-0556-x) contains supplementary material, which is available to authorized users.

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.

Risks and Recommendations in Prenatally Detected De Novo Balanced Chromosomal Rearrangements from Assessment of Long-Term Outcomes

The 6%-9% risk of an untoward outcome previously established by Warburton for prenatally detected de novo balanced chromosomal rearrangements (BCRs) does not account for long-term morbidity. We performed long-term follow-up (mean 17 years) of a registry-based nationwide cohort of 41 individuals carrying a prenatally detected de novo BCR with normal first trimester screening/ultrasound scan. We observed a significantly higher frequency of neurodevelopmental and/or neuropsychiatric disorders than in a matched control group (19.5% versus 8.3%, p = 0.04), which was increased to 26.8% upon clinical follow-up. Chromosomal microarray of 32 carriers revealed no pathogenic imbalances, illustrating a low prognostic value when fetal ultrasound scan is normal. In contrast, mate-pair sequencing revealed disrupted genes (ARID1B, NPAS3, CELF4), regulatory domains of known developmental genes (ZEB2, HOXC), and complex BCRs associated with adverse outcomes. Seven unmappable autosomal-autosomal BCRs with breakpoints involving pericentromeric/heterochromatic regions may represent a low-risk group. We performed independent phenotype-aware and blinded interpretation, which accurately predicted benign outcomes (specificity = 100%) but demonstrated relatively low sensitivity for prediction of the clinical outcome in affected carriers (sensitivity = 45%-55%). This sensitivity emphasizes the challenges associated with prenatal risk prediction for long-term morbidity in the absence of phenotypic data given the still immature annotation of the morbidity genome and poorly understood long-range regulatory mechanisms. In conclusion, we upwardly revise the previous estimates of Warburton to a morbidity risk of 27% and recommend sequencing of the chromosomal breakpoints as the first-tier diagnostic test in pregnancies with a de novo BCR.

Computational Prediction of Position Effects of Human Chromosome Rearrangements

Balanced and apparently balanced chromosome abnormalities (BCAs) have long been known to generate disease through position effects, either by altering local networks of gene regulation or positioning genes in architecturally different chromosome domains. Despite these observations, identification of distally affected genes by BCAs is oftentimes neglected, especially when predicted gene disruptions are found elsewhere in the genome. In this unit, we provide detailed instructions on how to run a computational pipeline that identifies relevant candidates of non-coding BCA position effects. This methodology facilitates quick identification of genes potentially involved in disease by non-coding BCAs and other types of rearrangements, and expands on the importance of considering the long-range consequences of genomic lesions.

Altered DNA methylation associated with a translocation linked to major mental illness

Recent work has highlighted a possible role for altered epigenetic modifications, including differential DNA methylation, in susceptibility to psychiatric illness. Here, we investigate blood-based DNA methylation in a large family where a balanced translocation between chromosomes 1 and 11 shows genome-wide significant linkage to psychiatric illness. Genome-wide DNA methylation was profiled in whole-blood-derived DNA from 41 individuals using the Infinium HumanMethylation450 BeadChip (Illumina Inc., San Diego, CA). We found significant differences in DNA methylation when translocation carriers (n = 17) were compared to related non-carriers (n = 24) at 13 loci. All but one of the 13 significant differentially methylated positions (DMPs) mapped to the regions surrounding the translocation breakpoints. Methylation levels of five DMPs were associated with genotype at SNPs in linkage disequilibrium with the translocation. Two of the five genes harbouring significant DMPs, DISC1 and DUSP10, have been previously shown to be differentially methylated in schizophrenia. Gene Ontology analysis revealed enrichment for terms relating to neuronal function and neurodevelopment among the genes harbouring the most significant DMPs. Differentially methylated region (DMR) analysis highlighted a number of genes from the MHC region, which has been implicated in psychiatric illness previously through genetic studies. We show that inheritance of a translocation linked to major mental illness is associated with differential DNA methylation at loci implicated in neuronal development/function and in psychiatric illness. As genomic rearrangements are over-represented in individuals with psychiatric illness, such analyses may be valuable more widely in the study of these conditions.

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.

Identification of OAF and PVRL1 as candidate genes for an ocular anomaly characterized by Peters anomaly type 2 and ectopia lentis

Keratolenticular dysgenesis (KLD) and ectopia lentis are congenital eye defects. The aim of this study is the identification of molecular genetic alterations responsible for those ocular anomalies with neurologic impairment in an individual with a de novo balanced chromosome translocation t(11;18)(q23.3;q11.2)dn. Disruption of OAF, the human orthologue of the Drosophila oaf, by the 11q23.3 breakpoint results in reduced expression of this transcriptional regulator. Furthermore, four most likely nonfunctional chimeric transcripts comprising up to OAF exon 3, derived from the der(11) allele, have also been identified. This locus has been implicated by publicly available genome-wide association data in corneal disease and corneal topography. The expression of the poliovirus receptor-related 1(PVRL1) or nectin cell adhesion molecule 1 (NECTIN1), a paralogue of nectin cell adhesion molecule 3 (PVRL3) associated with congenital ocular defects, situated 500 kb upstream from 11q23.3 breakpoint, is increased. The 18q11.2 breakpoint is localized between cutaneous T-cell lymphoma-associated antigen 1(CTAGE1) and retinoblastoma binding protein 8 (RBBP8) genes. Genomic imbalance that could contribute to the observed phenotype was excluded. Analysis of gene expression datasets throughout normal murine ocular lens embryogenesis suggests that OAF expression is significantly enriched in the lens from early stages of development through adulthood, whereas PVRL1 is lens-enriched until E12.5 and then down-regulated. This contrasts with the observation that the proposita's lymphoblastoid cell lines exhibit low OAF and high PVRL1 expression as compared to control, which offers further support that the alterations described above are most likely responsible for the clinical phenotype. Finally, gene interaction topology data for PVRL1 also agree with our proposal that disruption of OAF by the translocation breakpoint and misregulation of PVRL1 due to a position effect contribute to the observed ocular and neurological phenotype.

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.

Regulatory variants of FOXG1 in the context of its topological domain organisation

FOXG1 syndrome is caused by FOXG1 intragenic point mutations, or by long-range position effects (LRPE) of intergenic structural variants. However, the size of the FOXG1 regulatory landscape is uncertain, because the associated topologically associating domain (TAD) in fibroblasts is split into two domains in embryonic stem cells (hESC). Indeed, it has been suggested that the pathogenetic mechanism of deletions that remove the stem-cell-specific TAD boundary may be enhancer adoption due to ectopic activity of enhancer(s) located in the distal hESC-TAD. Herein we map three de novo translocation breakpoints to the proximal regulatory domain of FOXG1. The classical FOXG1 syndrome in these and in other translocation patients, and in a patient with an intergenic deletion that removes the hESC-specific TAD boundary, do not support the hypothesised enhancer adoption as a main contributor to the FOXG1 syndrome. Also, virtual 4 C and HiC-interaction data suggest that the hESC-specific TAD boundary may not be critical for FOXG1 regulation in a majority of human cells and tissues, including brain tissues and a neuronal progenitor cell line. Our data support the importance of a critical regulatory region (SRO) proximal to the hESC-specific TAD boundary. We further narrow this critical region by a deletion distal to the hESC-specific boundary, associated with a milder clinical phenotype. The distance from FOXG1 to the SRO ( > 500 kb) highlight a limitation of ENCODE DNase hypersensitivity data for functional prediction of LRPE. Moreover, the SRO has little overlap with a cluster of frequently associating regions (FIREs) located in the proximal hESC-TAD.

3C-PCR: a novel proximity ligation-based approach to phase chromosomal rearrangement breakpoints with distal allelic variants

Recent advances in molecular cytogenetics highlight the importance of noncoding structural variation in human disease. Genomic rearrangements can disrupt chromatin architecture, leading to long-range alterations in gene expression. With increasing ability to assess distal gene dysregulation comes new challenges in clinical interpretation of rearrangements. While haplotyping methods to determine compound heterozygosity in a single gene with two pathogenic variants are established, such methods are insufficient for phasing larger distances between a pathogenic variant and a genomic rearrangement breakpoint. Herein, we present an inexpensive and efficient proximity ligation-based method called 3C-PCR for phasing chromosomal rearrangement breakpoints with distal allelic variants. 3C-PCR uses canonical chromosome conformation capture (3C) libraries for targeted distal phasing by implementing a novel nested PCR strategy with primers anchored across the rearrangement breakpoints and subsequent Sanger sequencing. As a proof of concept, 3C-PCR was used to phase a highly variable region 1.3 Mb upstream of a chromosomal rearrangement breakpoint in a balanced translocation. We found that the nested PCR approach amplified the derivative chromosome substrate exclusively and identified the same haplotype by Sanger sequencing reliably. Given its efficacy and versatility, 3C-PCR is ideal for use in phasing chromosomal rearrangement breakpoints with allelic variants located at a genomic distance over a megabase.