Rare Genome-Wide Copy Number Variation and Expression of Schizophrenia in 22q11.2 Deletion Syndrome

OBJECTIVE

Chromosome 22q11.2 deletion syndrome (22q11.2DS) is associated with a more than 20-fold increased risk for developing schizophrenia. The aim of this study was to identify additional genetic factors (i.e., "second hits") that may contribute to schizophrenia expression.

METHOD

Through an international consortium, the authors obtained DNA samples from 329 psychiatrically phenotyped subjects with 22q11.2DS. Using a high-resolution microarray platform and established methods to assess copy number variation (CNV), the authors compared the genome-wide burden of rare autosomal CNV, outside of the 22q11.2 deletion region, between two groups: a schizophrenia group and those with no psychotic disorder at age ≥25 years. The authors assessed whether genes overlapped by rare CNVs were overrepresented in functional pathways relevant to schizophrenia.

RESULTS

Rare CNVs overlapping one or more protein-coding genes revealed significant between-group differences. For rare exonic duplications, six of 19 gene sets tested were enriched in the schizophrenia group; genes associated with abnormal nervous system phenotypes remained significant in a stepwise logistic regression model and showed significant interactions with 22q11.2 deletion region genes in a connectivity analysis. For rare exonic deletions, the schizophrenia group had, on average, more genes overlapped. The additional rare CNVs implicated known (e.g., GRM7, 15q13.3, 16p12.2) and novel schizophrenia risk genes and loci.

CONCLUSIONS

The results suggest that additional rare CNVs overlapping genes outside of the 22q11.2 deletion region contribute to schizophrenia risk in 22q11.2DS, supporting a multigenic hypothesis for schizophrenia. The findings have implications for understanding expression of psychotic illness and herald the importance of whole-genome sequencing to appreciate the overall genomic architecture of schizophrenia.

Nested Inversion Polymorphisms Predispose Chromosome 22q11.2 to Meiotic Rearrangements

Inversion polymorphisms between low-copy repeats (LCRs) might predispose chromosomes to meiotic non-allelic homologous recombination (NAHR) events and thus lead to genomic disorders. However, for the 22q11.2 deletion syndrome (22q11.2DS), the most common genomic disorder, no such inversions have been uncovered as of yet. Using fiber-FISH, we demonstrate that parents transmitting the de novo 3 Mb LCR22A-D 22q11.2 deletion, the reciprocal duplication, and the smaller 1.5 Mb LCR22A-B 22q11.2 deletion carry inversions of LCR22B-D or LCR22C-D. Hence, the inversions predispose chromosome 22q11.2 to meiotic rearrangements and increase the individual risk for transmitting rearrangements. Interestingly, the inversions are nested or flanking rather than coinciding with the deletion or duplication sizes. This finding raises the possibility that inversions are a prerequisite not only for 22q11.2 rearrangements but also for all NAHR-mediated genomic disorders.

Novel STIL Compound Heterozygous Mutations Cause Severe Fetal Microcephaly and Centriolar Lengthening

STIL (SCL/TAL1 interrupting locus) is a core component of the centriole duplication process. STIL mutations have been associated with both autosomal recessive primary microcephaly (MCPH) and holoprosencephaly. In this report, we describe a family with multiple miscarriages and 2 terminations of pregnancy due to marked fetal microcephaly, delayed cortical gyrification, and dysgenesis of the corpus callosum. Whole exome sequencing allowed us to identify novel compound heterozygous mutations in STIL. The mutations lie, respectively, in the CPAP/CENPJ and the hsSAS6 interacting domains of STIL. M-phase synchronized amniocytes from both affected fetuses did not display an aberrant number of centrioles, as shown previously for either STIL-depleted or overexpressing cells. However, we observed an elongation of at least 1 centriole for each duplicated centrosome. These preliminary results may point to a novel mechanism causing MCPH and embryonic lethality in humans.

Modulation of Aneuploidy in Leishmania donovani during Adaptation to Different In Vitro and In Vivo Environments and Its Impact on Gene Expression

Aneuploidy is usually deleterious in multicellular organisms but appears to be tolerated and potentially beneficial in unicellular organisms, including pathogens. Leishmania, a major protozoan parasite, is emerging as a new model for aneuploidy, since in vitro-cultivated strains are highly aneuploid, with interstrain diversity and intrastrain mosaicism. The alternation of two life stages in different environments (extracellular promastigotes and intracellular amastigotes) offers a unique opportunity to study the impact of environment on aneuploidy and gene expression. We sequenced the whole genomes and transcriptomes of Leishmania donovani strains throughout their adaptation to in vivo conditions mimicking natural vertebrate and invertebrate host environments. The nucleotide sequences were almost unchanged within a strain, in contrast to highly variable aneuploidy. Although high in promastigotes in vitro, aneuploidy dropped significantly in hamster amastigotes, in a progressive and strain-specific manner, accompanied by the emergence of new polysomies. After a passage through a sand fly, smaller yet consistent karyotype changes were detected. Changes in chromosome copy numbers were correlated with the corresponding transcript levels, but additional aneuploidy-independent regulation of gene expression was observed. This affected stage-specific gene expression, downregulation of the entire chromosome 31, and upregulation of gene arrays on chromosomes 5 and 8. Aneuploidy changes in Leishmania are probably adaptive and exploited to modulate the dosage and expression of specific genes; they are well tolerated, but additional mechanisms may exist to regulate the transcript levels of other genes located on aneuploid chromosomes. Our model should allow studies of the impact of aneuploidy on molecular adaptations and cellular fitness.IMPORTANCE Aneuploidy is usually detrimental in multicellular organisms, but in several microorganisms, it can be tolerated and even beneficial. Leishmania-a protozoan parasite that kills more than 30,000 people each year-is emerging as a new model for aneuploidy studies, as unexpectedly high levels of aneuploidy are found in clinical isolates. Leishmania lacks classical regulation of transcription at initiation through promoters, so aneuploidy could represent a major adaptive strategy of this parasite to modulate gene dosage in response to stressful environments. For the first time, we document the dynamics of aneuploidy throughout the life cycle of the parasite, in vitro and in vivo We show its adaptive impact on transcription and its interaction with regulation. Besides offering a new model for aneuploidy studies, we show that further genomic studies should be done directly in clinical samples without parasite isolation and that adequate methods should be developed for this.

Copy number variation analysis in adults with catatonia confirms haploinsufficiency of SHANK3 as a predisposing factor

BACKGROUND

Catatonia is a motor dysregulation syndrome co-occurring with a variety of psychiatric and medical disorders. Response to treatment with benzodiazepines and electroconvulsive therapy suggests a neurobiological background. The genetic etiology however remains largely unexplored. Copy Number Variants (CNV), known to predispose to neurodevelopmental disorders, may play a role in the etiology of catatonia.

METHODS

This study is exploring the genetic field of catatonia through CNV analysis in a cohort of psychiatric patients featuring intellectual disability and catatonia. Fifteen adults admitted to a psychiatric inpatient unit and diagnosed with catatonia were selected for array Comparative Genomic Hybridization analysis at 200 kb resolution. We introduced a CNV interpretation algorithm to define detected CNVs as benign, unclassified, likely pathogenic or causal with regard to catatonia.

RESULTS

Co-morbid psychiatric diagnoses in these patients were autism, psychotic or mood disorders. Eight patients were found to carry rare CNVs, which could not be classified as benign, comprising 6 duplications and 2 deletions. Microdeletions on 22q13.3, considered causal for catatonia, were detected in 2 patients. Duplications on 16p11.2 and 22q11.2 were previously implicated in psychiatric disorders, but not in catatonia, and were therefore considered likely pathogenic. Driven by the identification of a rare 14q11.2 duplication in one catatonic patient, additional patients with overlapping duplications were gathered to delineate a novel susceptibility locus for intellectual disability and psychiatric disorders on 14q11.2, harboring the gene SUPT16H. Three remaining variants respectively on 2q36.1, 16p13.13 and 17p13.3 were considered variants of unknown significance.

CONCLUSION

The identification of catatonia-related copy number changes in this study, underscores the importance of genetic research in patients with catatonia. We confirmed that 22q13.3 deletions, affecting the gene SHANK3, predispose to catatonia, and we uncover 14q11.2 duplications as a novel susceptibility factor for intellectual and psychiatric disorders.

Variant discovery and breakpoint region prediction for studying the human 22q11.2 deletion using BAC clone and whole genome sequencing analysis

Velo-cardio-facial syndrome/DiGeorge syndrome/22q11.2 deletion syndrome (22q11.2DS) is caused by meiotic non-allelic homologous recombination events between flanking low copy repeats termed LCR22A and LCR22D, resulting in a 3 million base pair (Mb) deletion. Due to their complex structure, large size and high sequence identity, genetic variation within LCR22s among different individuals has not been well characterized. In this study, we sequenced 13 BAC clones derived from LCR22A/D and aligned them with 15 previously available BAC sequences to create a new genetic variation map. The thousands of variants identified by this analysis were not uniformly distributed in the two LCR22s. Moreover, shared single nucleotide variants between LCR22A and LCR22D were enriched in the Breakpoint Cluster Region pseudogene (BCRP) block, suggesting the existence of a possible recombination hotspot there. Interestingly, breakpoints for atypical 22q11.2 rearrangements have previously been located to BCRPs To further explore this finding, we carried out in-depth analyses of whole genome sequence (WGS) data from two unrelated probands harbouring a de novo 3Mb 22q11.2 deletion and their normal parents. By focusing primarily on WGS reads uniquely mapped to LCR22A, using the variation map from our BAC analysis to help resolve allele ambiguity, and by performing PCR analysis, we infer that the deletion breakpoints were most likely located near or within the BCRP module. In summary, we found a high degree of sequence variation in LCR22A and LCR22D and a potential recombination breakpoint near or within the BCRP block, providing a starting point for future breakpoint mapping using additional trios.

Molecular Karyotyping: Array CGH Quality Criteria for Constitutional Genetic Diagnosis

Array CGH (comparative genomic hybridization) enables the identification of chromosomal copy number changes. The availability of clone sets covering the human genome opens the possibility for the widespread use of array CGH for both research and diagnostic purposes. In this manuscript we report on the parameters that were critical for successful implementation of the technology, assess quality criteria, and discuss the potential benefits and pitfalls of the technology for improved pre- and postnatal constitutional genetic diagnosis. We propose to name the genome-wide array CGH “molecular karyotyping,” in analogy with conventional karyotyping that uses staining methods to visualize chromosomes.

Identification of Intellectual Disability Genes in Female Patients with a Skewed X-Inactivation Pattern

Intellectual disability (ID) is a heterogeneous disorder with an unknown molecular etiology in many cases. Previously, X-linked ID (XLID) studies focused on males because of the hemizygous state of their X chromosome. Carrier females are generally unaffected because of the presence of a second normal allele, or inactivation of the mutant X chromosome in most of their cells (skewing). However, in female ID patients, we hypothesized that the presence of skewing of X-inactivation would be an indicator for an X chromosomal ID cause. We analyzed the X-inactivation patterns of 288 females with ID, and found that 22 (7.6%) had extreme skewing (>90%), which is significantly higher than observed in the general population (3.6%; P = 0.029). Whole-exome sequencing of 19 females with extreme skewing revealed causal variants in six females in the XLID genes DDX3X, NHS, WDR45, MECP2, and SMC1A. Interestingly, variants in genes escaping X-inactivation presumably cause both XLID and skewing of X-inactivation in three of these patients. Moreover, variants likely accounting for skewing only were detected in MED12, HDAC8, and TAF9B. All tested candidate causative variants were de novo events. Hence, extreme skewing is a good indicator for the presence of X-linked variants in female patients.

Zygotes segregate entire parental genomes in distinct blastomere lineages causing cleavage-stage chimerism and mixoploidy

Dramatic genome dynamics, such as chromosome instability, contribute to the remarkable genomic heterogeneity among the blastomeres comprising a single embryo during human preimplantation development. This heterogeneity, when compatible with life, manifests as constitutional mosaicism, chimerism, and mixoploidy in live-born individuals. Chimerism and mixoploidy are defined by the presence of cell lineages with different parental genomes or different ploidy states in a single individual, respectively. Our knowledge of their mechanistic origin results from indirect observations, often when the cell lineages have been subject to rigorous selective pressure during development. Here, we applied haplarithmisis to infer the haplotypes and the copy number of parental genomes in 116 single blastomeres comprising entire preimplantation bovine embryos (n = 23) following in vitro fertilization. We not only demonstrate that chromosome instability is conserved between bovine and human cleavage embryos, but we also discovered that zygotes can spontaneously segregate entire parental genomes into different cell lineages during the first post-zygotic cleavage division. Parental genome segregation was not exclusively triggered by abnormal fertilizations leading to triploid zygotes, but also normally fertilized zygotes can spontaneously segregate entire parental genomes into different cell lineages during cleavage of the zygote. We coin the term "heterogoneic division" to indicate the events leading to noncanonical zygotic cytokinesis, segregating the parental genomes into distinct cell lineages. Persistence of those cell lines during development is a likely cause of chimerism and mixoploidy in mammals.

Reciprocal 22q11.2 Deletion and Duplication in Siblings with Karyotypically Normal Parents

The 22q11.2 locus is known to harbor a high risk for structural variation caused by non-allelic homologous recombination, resulting in deletions and duplications. Here, we describe the first family with one sibling carrying the 22q11 deletion and the other carrying the reciprocal duplication. FISH and SNP array analysis of the parents show a maternal origin for both deletion and duplication, without indications of balanced deletions/duplications or mosaicism. We hypothesize that germline mosaicism in the mother underlies the deletion and duplication, which would implicate a high recurrence risk for her offspring.

Polymerase specific error rates and profiles identified by single molecule sequencing

DNA polymerases have an innate error rate which is polymerase and DNA context specific. Historically the mutational rate and profiles have been measured using a variety of methods, each with their own technical limitations. Here we used the unique properties of single molecule sequencing to evaluate the mutational rate and profiles of six DNA polymerases at the sequence level. In addition to accurately determining mutations in double strands, single molecule sequencing also captures direction specific transversions and transitions through the analysis of heteroduplexes. Not only did the error rates vary, but also the direction specific transitions differed among polymerases.

Pulmonary transcriptome analysis in the surgically induced rabbit model of diaphragmatic hernia treated with fetal tracheal occlusion

Congenital diaphragmatic hernia (CDH) is a malformation leading to pulmonary hypoplasia, which can be treated in utero by fetal tracheal occlusion (TO). However, the changes of gene expression induced by TO remain largely unknown but could be used to further improve the clinically used prenatal treatment of this devastating malformation. Therefore, we aimed to investigate the pulmonary transcriptome changes caused by surgical induction of diaphragmatic hernia (DH) and additional TO in the fetal rabbit model. Induction of DH was associated with 378 upregulated genes compared to controls when allowing a false-discovery rate (FDR) of 0.1 and a fold change (FC) of 2. Those genes were again downregulated by consecutive TO. But DH+TO was associated with an upregulation of 157 genes compared to DH and controls. When being compared to control lungs, 106 genes were downregulated in the DH group and were not changed by TO. Therefore, the overall pattern of gene expression in DH+TO is more similar to the control group than to the DH group. In this study, we further provide a database of gene expression changes induced by surgical creation of DH and consecutive TO in the rabbit model. Future treatment strategies could be developed using this dataset. We also discuss the most relevant genes that are involved in CDH.

Summary: Rabbit fetuses with induced diaphragmatic hernia and treated with prenatal tracheal occlusion have a similar pulmonary transcriptome as unaffected controls. This study describes a valuable database of gene expressions in this model.

Clinical implementation of NIPT - technical and biological challenges

Non-invasive prenatal testing (NIPT) for fetal aneuploidy detection is increasingly being offered in the clinical setting. Whereas the majority of tests only report fetal trisomies 21, 18 and 13, genome-wide analyses have the potential to detect other fetal, as well as maternal, aneuploidies. In this review, we discuss the technical and clinical advantages and challenges associated with genome-wide cell-free fetal DNA profiling.

Coffin-Siris and Nicolaides-Baraitser syndromes are a common well recognizable cause of intellectual disability

BACKGROUND

Nicolaides-Baraitser and Coffin-Siris syndromes are emerging conditions with overlapping clinical features including intellectual disability and typical somatic characteristics, especially sparse hair, low frontal hairline, large mouth with thick and everted lips, and hands and feet anomalies. Since 2012, mutations in genes encoding six proteins of the BAF complex were identified in both conditions.

METHODS AND RESULTS

We have clinically evaluated a cohort of 1161 patients with intellectual disability from three different Italian centers. A strong clinical suspicion of either Nicolaides-Baraitser syndrome or Coffin-Siris syndrome was proposed in 11 cases who were then molecularly confirmed: 8 having de novo missense mutations in SMARCA2, two frame-shift mutations in ARID1B and one missense mutation in SMARCB1. Given the high frequency of the condition we set up a one-step deep sequencing test for all 6 genes of the BAF complex.

CONCLUSIONS

These results prove that the frequency of these conditions may be as high as the most common syndromes with intellectual deficit (about 1%). Clinical geneticists should be well aware of this group of disorders in the clinical setting when ascertaining patients with intellectual deficit, the specific facial features being the major diagnostic handle. Finally, this work adds information on the clinical differences of the two conditions and presents a fast and sensitive test for the molecular diagnosis.

Clinical implementation of NIPT - technical and biological challenges

Non-invasive prenatal testing (NIPT) for fetal aneuploidy detection is increasingly being offered in the clinical setting. Whereas the majority of tests only report fetal trisomies 21, 18 and 13, genome-wide analyses have the potential to detect other fetal, as well as maternal, aneuploidies. In this review, we discuss the technical and clinical advantages and challenges associated with genome-wide cell-free fetal DNA profiling.

GBSX: a toolkit for experimental design and demultiplexing genotyping by sequencing experiments

Background

Massive parallel sequencing is a powerful tool for variant discovery and genotyping. To reduce costs, sequencing of restriction enzyme based reduced representation libraries can be utilized. This technology is generally referred to as Genotyping By Sequencing (GBS). To deal with GBS experimental design and initial processing specific bioinformatic tools are needed.

Results

GBSX is a package that assists in selecting the appropriate enzyme and the design of compatible in-line barcodes. Post sequencing, it performs optimized demultiplexing using these barcodes to create fastq files per barcode which can easily be plugged into existing variant analysis pipelines. Here we demonstrate the usability of the GBSX toolkit and demonstrate improved in-line barcode demultiplexing and trimming performance compared to existing tools.

Conclusions

GBSX provides an easy to use suite of tools for designing and demultiplexing of GBS experiments.

Non-invasive detection of genomic imbalances in Hodgkin/Reed-Sternberg cells in early and advanced stage Hodgkin's lymphoma by sequencing of circulating cell-free DNA: a technical proof-of-principle study

BACKGROUND

Hodgkin's lymphoma is one of the most common lymphoid neoplasms in young adults, but the low abundance of neoplastic Hodgkin/Reed-Sternberg cells in the tumour hampers the elucidation of its pathogenesis, biology, and diversity. After an incidental observation that genomic aberrations known to occur in Hodgkin's lymphoma were detectable in circulating cell-free DNA, this study was undertaken to investigate whether circulating cell-free DNA can be informative about genomic imbalances in Hodgkin's lymphoma.

METHODS

We applied massive parallel sequencing to circulating cell-free DNA in a prospective study of patients with biopsy proven nodular sclerosis Hodgkin's lymphoma. Genomic imbalances in Hodgkin/Reed-Sternberg cells were investigated by fluorescence in-situ hybridisation (FISH) on tumour specimens.

FINDINGS

By non-invasive prenatal testing, we observed several genomic imbalances in circulating cell-free DNA of a pregnant woman, who was subsequently diagnosed with early-stage nodular sclerosis Hodgkin's lymphoma stage IIA during gestation. FISH on tumour tissue confirmed corresponding genomic imbalances in Hodgkin/Reed-Sternberg cells. We prospectively studied circulating cell-free DNA of nine nodular sclerosis Hodgkin's lymphoma cases: eight at first diagnosis and one at first relapse. Seven patients had stage IIA disease and two had stage IVB disease. In eight, genomic imbalances were detected, including, among others, gain of chromosomes 2p and 9p, known to occur in Hodgkin's lymphoma. These gains and losses in circulating cell-free DNA were extensively validated by FISH on Hodgkin/Reed-Sternberg cells in biopsy samples. Initiation of chemotherapy induced normalisation of circulating cell-free DNA profiles within 2-6 weeks. The cell cycle indicator Ki67 and cleaved caspase-3 were detected in Hodgkin/Reed-Sternberg cells by immunohistochemistry, suggesting high turnover of Hodgkin/Reed-Sternberg cells.

INTERPRETATION

In early and advanced stage nodular sclerosis Hodgkin's lymphoma, genomic imbalances in Hodgkin/Reed-Sternberg cells can be identified by massive parallel sequencing of circulating cell-free DNA at diagnosis. The rapid normalisation of circulating cell-free DNA profiles on therapy initiation suggests a potential role for circulating cell-free DNA profiling in early response monitoring. This finding creates several new possibilities for exploring the diversity of Hodgkin's lymphoma, and has potential implications for the future clinical development of biomarkers and precision therapy for this malignancy.

FUNDING

KU Leuven-University of Leuven and University Hospitals Leuven.

Noninvasive prenatal testing using a novel analysis pipeline to screen for all autosomal fetal aneuploidies improves pregnancy management

Noninvasive prenatal testing by massive parallel sequencing of maternal plasma DNA has rapidly been adopted as a mainstream method for detection of fetal trisomy 21, 18 and 13. Despite the relative high accuracy of current NIPT testing, a substantial number of false-positive and false-negative test results remain. Here, we present an analysis pipeline, which addresses some of the technical as well as the biologically derived causes of error. Most importantly, it differentiates high z-scores due to fetal trisomies from those due to local maternal CNVs causing false positives. This pipeline was retrospectively validated for trisomy 18 and 21 detection on 296 samples demonstrating a sensitivity and specificity of 100%, and applied prospectively to 1350 pregnant women in the clinical diagnostic setting with a result reported in 99.9% of cases. In addition, values indicative for trisomy were observed two times for chromosome 7 and once each for chromosomes 15 and 16, and once for a segmental trisomy 18. Two of the trisomies were confirmed to be mosaic, one of which contained a uniparental disomy cell line. As placental trisomies pose a risk for low-grade fetal mosaicism as well as uniparental disomy, genome-wide noninvasive aneuploidy detection is improving prenatal management.

Expanding the phenotypic spectrum of PORCN variants in two males with syndromic microphthalmia

Variants in PORCN are a cause of Goltz-Gorlin syndrome or Focal Dermal Hypoplasia, an X-linked dominant disorder affecting heterozygous females and until now considered to be embryonic lethal in males. Exome sequencing was performed in a family in which two male siblings were characterized by microphthalmia and additional congenital anomalies including diaphragmatic hernia, spina bifida and cardiac defects. Surprisingly, we identified a maternally inherited variant in PORCN present in both males as well as in two female siblings. This represents the first finding of a PORCN variant in non-mosaic males affected with Goltz-Gorlin syndrome. The apparently asymptomatic mother showed extreme skewing of X-inactivation (90%), an asymptomatic female sibling showed skewing of 88%, and the second female sibling affected with cutis aplasia of the scalp showed X-inactivation considered within the normal range.