Complex rearrangements in patients with duplications of MECP2 can occur by fork stalling and template switching

Structural variant allelic heterogeneity in MECP2 duplication syndrome provides insight into clinical severity and variability of disease expression

BACKGROUND

MECP2 Duplication Syndrome, also known as X-linked intellectual developmental disorder Lubs type (MRXSL; MIM: 300260), is a neurodevelopmental disorder caused by copy number gains spanning MECP2. Despite varying genomic rearrangement structures, including duplications and triplications, and a wide range of duplication sizes, no clear correlation exists between DNA rearrangement and clinical features. We had previously demonstrated that up to 38% of MRXSL families are characterized by complex genomic rearrangements (CGRs) of intermediate complexity (2 ≤ copy number variant breakpoints < 5), yet the impact of these genomic structures on regulation of gene expression and phenotypic manifestations have not been investigated.

METHODS

To study the role of the genomic rearrangement structures on an individual's clinical phenotypic variability, we employed a comprehensive genomics, transcriptomics, and deep phenotyping analysis approach on 137 individuals affected by MRXSL. Genomic structural information was correlated with transcriptomic and quantitative phenotypic analysis using Human Phenotype Ontology (HPO) semantic similarity scores.

RESULTS

Duplication sizes in the cohort ranging from 64.6 kb to 16.5 Mb were classified into four categories comprising of tandem duplications (48%), terminal duplications (22%), inverted triplications (20%), and other CGRs (10%). Most of the terminal duplication structures consist of translocations (65%) followed by recombinant chromosomes (23%). Notably, 65% of de novo events occurred in the Terminal duplication group in contrast with 17% observed in Tandem duplications. RNA-seq data from lymphoblastoid cell lines indicated that the MECP2 transcript quantity in MECP2 triplications is statistically different from all duplications, but not between other classes of genomic structures. We also observed a significant (p < 0.05) correlation (Pearson R = 0.6, Spearman p = 0.63) between the log-transformed MECP2 RNA levels and MECP2 protein levels, demonstrating that genomic aberrations spanning MECP2 lead to altered MECP2 RNA and MECP2 protein levels. Genotype-phenotype analyses indicated a gradual worsening of phenotypic features, including overall survival, developmental levels, microcephaly, epilepsy, and genitourinary/eye abnormalities in the following order: Tandem duplications, Other complex duplications, Terminal duplications/Translocations, and Triplications encompassing MECP2.

CONCLUSION

In aggregate, this combined analysis uncovers an interplay between MECP2 dosage, genomic rearrangement structure and phenotypic traits. Whereas the level of MECP2 is a key determinant of the phenotype, the DNA rearrangement structure can contribute to clinical severity and disease expression variability. Employing this type of analytical approach will advance our understanding of the impact of genomic rearrangements on genomic disorders and may help guide more targeted therapeutic approaches.

Structural variations in livestock genomes and their associations with phenotypic traits: a review

Genomic structural variation (SV) refers to differences in gene sequences between individuals on a genomic scale. It is widely distributed in the genome, primarily in the form of insertions, deletions, duplications, inversions, and translocations. Due to its characterization by long segments and large coverage, SVs significantly impact the genetic characteristics and production performance of livestock, playing a crucial role in studying breed diversity, biological evolution, and disease correlation. Research on SVs contributes to an enhanced understanding of chromosome function and genetic characteristics and is important for understanding hereditary diseases mechanisms. In this article, we review the concept, classification, main formation mechanisms, detection methods, and advancement of research on SVs in the genomes of cattle, buffalo, equine, sheep, and goats, aiming to reveal the genetic basis of differences in phenotypic traits and adaptive genetic mechanisms through genomic research, which will provide a theoretical basis for better understanding and utilizing the genetic resources of herbivorous livestock.

Modeling antisense oligonucleotide therapy in MECP2 duplication syndrome human iPSC-derived neurons reveals gene expression programs responsive to MeCP2 levels

Genomic copy-number variations (CNVs) that can cause neurodevelopmental disorders often encompass many genes, which complicates our understanding of how individual genes within a CNV contribute to pathology. MECP2 duplication syndrome (MDS or MRXSL in OMIM; OMIM#300260) is one such CNV disorder caused by duplications spanning methyl CpG-binding protein 2 (MECP2) and other genes on Xq28. Using an antisense oligonucleotide (ASO) to normalize MECP2 dosage is sufficient to rescue abnormal neurological phenotypes in mouse models overexpressing MECP2 alone, implicating the importance of increased MECP2 dosage within CNVs of Xq28. However, because MDS CNVs span MECP2 and additional genes, we generated human neurons from multiple MDS patient-derived induced pluripotent cells (iPSCs) to evaluate the benefit of using an ASO against MECP2 in a MDS human neuronal context. Importantly, we identified a signature of genes that is partially and qualitatively modulated upon ASO treatment, pinpointed genes sensitive to MeCP2 function, and altered in a model of Rett syndrome, a neurological disorder caused by loss of MeCP2 function. Furthermore, the signature contained genes that are aberrantly altered in unaffected control human neurons upon MeCP2 depletion, revealing gene expression programs qualitatively sensitive to MeCP2 levels in human neurons. Lastly, ASO treatment led to a partial rescue of abnormal neuronal morphology in MDS neurons. All together, these data demonstrate that ASOs targeting MECP2 benefit human MDS neurons. Moreover, our study establishes a paradigm by which to evaluate the contribution of individual genes within a CNV to pathogenesis and to assess their potential as a therapeutic target.

Detection of mosaic and population-level structural variants with Sniffles2

Calling structural variations (SVs) is technically challenging, but using long reads remains the most accurate way to identify complex genomic alterations. Here we present Sniffles2, which improves over current methods by implementing a repeat aware clustering coupled with a fast consensus sequence and coverage-adaptive filtering. Sniffles2 is 11.8 times faster and 29% more accurate than state-of-the-art SV callers across different coverages (5-50×), sequencing technologies (ONT and HiFi) and SV types. Furthermore, Sniffles2 solves the problem of family-level to population-level SV calling to produce fully genotyped VCF files. Across 11 probands, we accurately identified causative SVs around MECP2, including highly complex alleles with three overlapping SVs. Sniffles2 also enables the detection of mosaic SVs in bulk long-read data. As a result, we identified multiple mosaic SVs in brain tissue from a patient with multiple system atrophy. The identified SV showed a remarkable diversity within the cingulate cortex, impacting both genes involved in neuron function and repetitive elements.

Genomic Balancing Act: deciphering DNA rearrangements in the complex chromosomal aberration involving 5p15.2, 2q31.1, and 18q21.32

Despite extensive research into the genetic underpinnings of neurodevelopmental disorders (NDD), many clinical cases remain unresolved. We studied a female proband with a NDD, mildly dysmorphic facial features, and brain stem hypoplasia on neuroimaging. Comprehensive genomic analyses revealed a terminal 5p loss and a terminal 18q gain in the proband while a diploid copy number for chromosomes 5 and 18 in both parents. Genomic investigations in the proband identified an unbalanced translocation t(5;18) with additional genetic material from chromosome 2 (2q31.3) inserted at the breakpoint, pointing to a complex chromosomal rearrangement (CCR) involving 5p15.2, 2q31.3, and 18q21.32. Breakpoint junction analyses enabled by long-read genome sequencing unveiled the presence of four distinct junctions in the father, who is a carrier of a balanced CCR. The proband inherited from the father both the abnormal chromosome 5 resulting in segmental aneusomies of chr5 (loss) and chr18 (gain) and a der(2) homologue. Evidences suggest a chromoplexy mechanism for this CCR derivation, involving double-strand breaks (DSBs) repaired by non-homologous end joining (NHEJ) or alternative end joining (alt-EJ). The complexity of the CCR and the segregation of homologues elucidate the genetic model for this family. This study demonstrates the importance of combining multiple genomic technologies to uncover genetic causes of complex neurodevelopmental syndromes and to better understand genetic disease mechanisms.

Complex chromosomal 6q rearrangements revealed by combined long-molecule genomics technologies

Technologies for detecting structural variation (SV) have advanced with the advent of long-read sequencing, which enables the validation of SV at a nucleotide level. Optical genome mapping (OGM), a technology based on physical mapping, can also provide comprehensive SVs analysis. We applied long-read whole genome sequencing (LRWGS) to accurately reconstruct breakpoint (BP) segments in a patient with complex chromosome 6q rearrangements that remained elusive by conventional karyotyping. Although all BPs were precisely identified by LRWGS, there were two possible ways to construct the BP segments in terms of their orders and orientations. Thus, we also used OGM analysis. Notably, OGM recognized entire inversions exceeding 500 kb in size, which LRWGS could not characterize. Consequently, here we successfully unveil the full genomic structure of this complex chromosomal 6q rearrangement and cryptic SVs through combined long-molecule genomic analyses, showcasing how LRWGS and OGM can complement each other in SV 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.

NODAL variants are associated with a continuum of laterality defects from simple D-transposition of the great arteries to heterotaxy

BACKGROUND

NODAL signaling plays a critical role in embryonic patterning and heart development in vertebrates. Genetic variants resulting in perturbations of the TGF-β/NODAL signaling pathway have reproducibly been shown to cause laterality defects in humans. To further explore this association and improve genetic diagnosis, the study aims to identify and characterize a broader range of NODAL variants in a large number of individuals with laterality defects.

METHODS

We re-analyzed a cohort of 321 proband-only exomes of individuals with clinically diagnosed laterality congenital heart disease (CHD) using family-based, rare variant genomic analyses. To this cohort we added 12 affected subjects with known NODAL variants and CHD from institutional research and clinical cohorts to investigate an allelic series. For those with candidate contributory variants, variant allele confirmation and segregation analysis were studied by Sanger sequencing in available family members. Array comparative genomic hybridization and droplet digital PCR were utilized for copy number variants (CNV) validation and characterization. We performed Human Phenotype Ontology (HPO)-based quantitative phenotypic analyses to dissect allele-specific phenotypic differences.

RESULTS

Missense, nonsense, splice site, indels, and/or structural variants of NODAL were identified as potential causes of heterotaxy and other laterality defects in 33 CHD cases. We describe a recurrent complex indel variant for which the nucleic acid secondary structure predictions implicate secondary structure mutagenesis as a possible mechanism for formation. We identified two CNV deletion alleles spanning NODAL in two unrelated CHD cases. Furthermore, 17 CHD individuals were found (16/17 with known Hispanic ancestry) to have the c.778G > A:p.G260R NODAL missense variant which we propose reclassification from variant of uncertain significance (VUS) to likely pathogenic. Quantitative HPO-based analyses of the observed clinical phenotype for all cases with p.G260R variation, including heterozygous, homozygous, and compound heterozygous cases, reveal clustering of individuals with biallelic variation. This finding provides evidence for a genotypic-phenotypic correlation and an allele-specific gene dosage model.

CONCLUSION

Our data further support a role for rare deleterious variants in NODAL as a cause for sporadic human laterality defects, expand the repertoire of observed anatomical complexity of potential cardiovascular anomalies, and implicate an allele specific gene dosage model.

HMZDupFinder: a robust computational approach for detecting intragenic homozygous duplications from exome sequencing data

Homozygous duplications contribute to genetic disease by altering gene dosage or disrupting gene regulation and can be more deleterious to organismal biology than heterozygous duplications. Intragenic exonic duplications can result in loss-of-function (LoF) or gain-of-function (GoF) alleles that when homozygosed, i.e. brought to homozygous state at a locus by identity by descent or state, could potentially result in autosomal recessive (AR) rare disease traits. However, the detection and functional interpretation of homozygous duplications from exome sequencing data remains a challenge. We developed a framework algorithm, HMZDupFinder, that is designed to detect exonic homozygous duplications from exome sequencing (ES) data. The HMZDupFinder algorithm can efficiently process large datasets and accurately identifies small intragenic duplications, including those associated with rare disease traits. HMZDupFinder called 965 homozygous duplications with three or less exons from 8,707 ES with a recall rate of 70.9% and a precision of 16.1%. We experimentally confirmed 8/10 rare homozygous duplications. Pathogenicity assessment of these copy number variant alleles allowed clinical genomics contextualization for three homozygous duplications alleles, including two affecting known OMIM disease genes EDAR (MIM# 224900), TNNT1(MIM# 605355), and one variant in a novel candidate disease gene: PAAF1.

MECP2-related disorders while gene-based therapies are on the horizon

The emergence of new genetic tools has led to the discovery of the genetic bases of many intellectual and developmental disabilities. This creates exciting opportunities for research and treatment development, and a few genetic disorders (e.g., spinal muscular atrophy) have recently been treated with gene-based therapies. MECP2 is found on the X chromosome and regulates the transcription of thousands of genes. Loss of MECP2 gene product leads to Rett Syndrome, a disease found primarily in females, and is characterized by developmental regression, motor dysfunction, midline hand stereotypies, autonomic nervous system dysfunction, epilepsy, scoliosis, and autistic-like behavior. Duplication of MECP2 causes MECP2 Duplication Syndrome (MDS). MDS is found mostly in males and presents with developmental delay, hypotonia, autistic features, refractory epilepsy, and recurrent respiratory infections. While these two disorders share several characteristics, their differences (e.g., affected sex, age of onset, genotype/phenotype correlations) are important to distinguish in the light of gene-based therapy because they require opposite solutions. This review explores the clinical features of both disorders and highlights these important clinical differences.

Generation of five induced pluripotent stem cell lines from patients with MECP2 Duplication Syndrome

MECP2 Duplication Syndrome (MDS) is a rare, severe neurodevelopmental disorder arising from duplications in the Xq28 region containing the MECP2 gene that predominantly affects males. We generated five human induced pluripotent stem cell (iPSC) lines from the fibroblasts of individuals carrying between 0.355 and 11.2 Mb size duplications in the chromosomal locus containing MECP2. All lines underwent extensive testing to confirm MECP2 duplication and iPSC-related features such as morphology, pluripotency markers, and trilineage differentiation potential. These lines are a valuable resource for molecular and functional studies of MDS as well as screening for a variety of therapeutic approaches.

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 .

Genomic Complexity and Complex Chromosomal Rearrangements in Genetic Diagnosis: Two Illustrative Cases on Chromosome 7

Complex chromosomal rearrangements are rare events compatible with survival, consisting of an imbalance and/or position effect of one or more genes, that contribute to a range of clinical presentations. The investigation and diagnosis of these cases are often difficult. The interpretation of the pattern of pairing and segregation of these chromosomes during meiosis is important for the assessment of the risk and the type of imbalance in the offspring. Here, we investigated two unrelated pediatric carriers of complex rearrangements of chromosome 7. The first case was a 2-year-old girl with a severe phenotype. Conventional cytogenetics evidenced a duplication of part of the short arm of chromosome 7. By array-CGH analysis, we found a complex rearrangement with three discontinuous trisomy regions (7p22.1p21.3, 7p21.3, and 7p21.3p15.3). The second case was a newborn investigated for hypodevelopment and dimorphisms. The karyotype analysis promptly revealed a structurally altered chromosome 7. The array-CGH analysis identified an even more complex rearrangement consisting of a trisomic region at 7q11.23q22 and a tetrasomic region of 4.5 Mb spanning 7q21.3 to q22.1. The mother's karyotype examination revealed a complex rearrangement of chromosome 7: the 7q11.23q22 region was inserted in the short arm at 7p15.3. Finally, array-CGH analysis showed a trisomic region that corresponds to the tetrasomic region of the son. Our work proved that the integration of several technical solutions is often required to appropriately analyze complex chromosomal rearrangements in order to understand their implications and offer appropriate genetic counseling.

Rare variant enrichment analysis supports GREB1L as a contributory driver gene in the etiology of Mayer-Rokitansky-Küster-Hauser syndrome

Mayer-Rokitansky-Küster-Hauser (MRKH) syndrome is characterized by aplasia of the female reproductive tract; the syndrome can include renal anomalies, absence or dysgenesis, and skeletal anomalies. While functional models have elucidated several candidate genes, only WNT4 (MIM: 603490) variants have been definitively associated with a subtype of MRKH with hyperandrogenism (MIM: 158330). DNA from 148 clinically diagnosed MRKH probands across 144 unrelated families and available family members from North America, Europe, and South America were exome sequenced (ES) and by family-based genomics analyzed for rare likely deleterious variants. A replication cohort consisting of 442 Han Chinese individuals with MRKH was used to further reproduce GREB1L findings in diverse genetic backgrounds. Proband and OMIM phenotypes annotated using the Human Phenotype Ontology were analyzed to quantitatively delineate the phenotypic spectrum associated with GREB1L variant alleles found in our MRKH cohort and those previously published. This study reports 18 novel GREB1L variant alleles, 16 within a multiethnic MRKH cohort and two within a congenital scoliosis cohort. Cohort-wide analyses for a burden of rare variants within a single gene identified likely damaging variants in GREB1L (MIM: 617782), a known disease gene for renal hypoplasia and uterine abnormalities (MIM: 617805), in 16 of 590 MRKH probands. GREB1L variant alleles, including a CNV null allele, were found in 8 MRKH type 1 probands and 8 MRKH type II probands. This study used quantitative phenotypic analyses in a worldwide multiethnic cohort to identify and strengthen the association of GREB1L to isolated uterine agenesis (MRKH type I) and syndromic MRKH type II.

Long-read sequence analysis for clustered genomic copy number aberrations revealed architectures of intricately intertwined rearrangements

Most chromosomal aberrations revealed by chromosomal microarray testing (CMA) are simple; however, very complex chromosomal structural rearrangements can also be found. Although the mechanism of structural rearrangements has been gradually revealed, not all mechanisms have been elucidated. We analyzed the breakpoint-junctions (BJs) of two or more clustered copy number variations (CNVs) in the same chromosome arms to understand their conformation and the mechanism of complex structural rearrangements. Combining CMA with long-read whole-genome sequencing (WGS) analysis, we successfully determined all BJs for the clustered CNVs identified in four patients. Multiple CNVs were intricately intertwined with each other, and clustered CNVs in four patients were involved in global complex chromosomal rearrangements. The BJs of two clustered deletions identified in two patients showed microhomologies, and their characteristics were explained by chromothripsis. In contrast, the BJs in the other two patients, who showed clustered deletions and duplications, consisted of blunt-end and nontemplated insertions. These findings could be explained only by alternative nonhomologous end-joining, a mechanism related to polymerase theta. All the patients had at least one inverted segment. Three patients showed cryptic aberrations involving a disruption and a deletion/duplication, which were not detected by CMA but were first identified by WGS. This result suggested that complex rearrangements should be considered if clustered CNVs are observed in the same chromosome arms. Because CMA has potential limitations in genotype-phenotype correlation analysis, a more detailed analysis by whole genome examination is recommended in cases of suspected complex structural aberrations.

Complex genomic rearrangements: an underestimated cause of rare diseases

Complex genomic rearrangements (CGRs) are known contributors to disease but are often missed during routine genetic screening. Identifying CGRs requires (i) identifying copy number variants (CNVs) concurrently with inversions, (ii) phasing multiple breakpoint junctions incis, as well as (iii) detecting and resolving structural variants (SVs) within repeats. We demonstrate how combining cytogenetics and new sequencing methodologies is being successfully applied to gain insights into the genomic architecture of CGRs. In addition, we review CGR patterns and molecular features revealed by studying constitutional genomic disorders. These data offer invaluable lessons to individuals interested in investigating CGRs, evaluating their clinical relevance and frequency, as well as assessing their impact(s) on rare genetic diseases.

The multiple de novo copy number variant (MdnCNV) phenomenon presents with peri-zygotic DNA mutational signatures and multilocus pathogenic variation

BACKGROUND

The multiple de novo copy number variant (MdnCNV) phenotype is described by having four or more constitutional de novo CNVs (dnCNVs) arising independently throughout the human genome within one generation. It is a rare peri-zygotic mutational event, previously reported to be seen once in every 12,000 individuals referred for genome-wide chromosomal microarray analysis due to congenital abnormalities. These rare families provide a unique opportunity to understand the genetic factors of peri-zygotic genome instability and the impact of dnCNV on human diseases.

METHODS

Chromosomal microarray analysis (CMA), array-based comparative genomic hybridization, short- and long-read genome sequencing (GS) were performed on the newly identified MdnCNV family to identify de novo mutations including dnCNVs, de novo single-nucleotide variants (dnSNVs), and indels. Short-read GS was performed on four previously published MdnCNV families for dnSNV analysis. Trio-based rare variant analysis was performed on the newly identified individual and four previously published MdnCNV families to identify potential genetic etiologies contributing to the peri-zygotic genomic instability. Lin semantic similarity scores informed quantitative human phenotype ontology analysis on three MdnCNV families to identify gene(s) driving or contributing to the clinical phenotype.

RESULTS

In the newly identified MdnCNV case, we revealed eight de novo tandem duplications, each ~ 1 Mb, with microhomology at 6/8 breakpoint junctions. Enrichment of de novo single-nucleotide variants (SNV; 6/79) and de novo indels (1/12) was found within 4 Mb of the dnCNV genomic regions. An elevated post-zygotic SNV mutation rate was observed in MdnCNV families. Maternal rare variant analyses identified three genes in distinct families that may contribute to the MdnCNV phenomenon. Phenotype analysis suggests that gene(s) within dnCNV regions contribute to the observed proband phenotype in 3/3 cases. CNVs in two cases, a contiguous gene duplication encompassing PMP22 and RAI1 and another duplication affecting NSD1 and SMARCC2, contribute to the clinically observed phenotypic manifestations.

CONCLUSIONS

Characteristic features of dnCNVs reported here are consistent with a microhomology-mediated break-induced replication (MMBIR)-driven mechanism during the peri-zygotic period. Maternal genetic variants in DNA repair genes potentially contribute to peri-zygotic genomic instability. Variable phenotypic features were observed across a cohort of three MdnCNV probands, and computational quantitative phenotyping revealed that two out of three had evidence for the contribution of more than one genetic locus to the proband's phenotype supporting the hypothesis of de novo multilocus pathogenic variation (MPV) in those families.

Systematic Analysis of Copy Number Variations in the Pathogenic Yeast Candida parapsilosis Identifies a Gene Amplification in RTA3 That is Associated with Drug Resistance

We analyzed the genomes of 170 C. parapsilosis isolates and identified multiple copy number variations (CNVs). We identified two genes, RTA3 (CPAR2_104610) and ARR3 (CPAR2_601050), each of which was the target of multiple independent amplification events. Phylogenetic analysis shows that most of these amplifications originated only once. For ARR3, which encodes a putative arsenate transporter, 8 distinct CNVs were identified, ranging in size from 2.3 kb to 10.5 kb with 3 to 23 copies. For RTA3, 16 distinct CNVs were identified, ranging in size from 0.3 kb to 4.5 kb with 2 to ~50 copies. One unusual amplification resulted in a DUP-TRP/INV-DUP structure similar to some human CNVs. RTA3 encodes a putative phosphatidylcholine (PC) floppase which is known to regulate the inward translocation of PC in Candida albicans. We found that an increased copy number of RTA3 correlated with resistance to miltefosine, an alkylphosphocholine drug that affects PC metabolism. Additionally, we conducted an adaptive laboratory evolution experiment in which two C. parapsilosis isolates were cultured in increasing concentrations of miltefosine. Two genes, CPAR2_303950 and CPAR2_102700, coding for putative PC flippases homologous to S. cerevisiae DNF1 gained homozygous protein-disrupting mutations in the evolved strains. Overall, our results show that C. parapsilosis can gain resistance to miltefosine, a drug that has recently been granted orphan drug designation approval by the United States Food and Drug Administration for the treatment of invasive candidiasis, through both CNVs or loss-of-function alleles in one of the flippase genes. IMPORTANCE Copy number variations (CNVs) are an important source of genomic diversity that have been associated with drug resistance. We identify two unusual CNVs in the human fungal pathogen Candida parapsilosis. Both target a single gene (RTA3 or ARR3), and they have occurred multiple times in multiple isolates. The copy number of RTA3, a putative floppase that controls the inward translocation of lipids in the cell membrane, correlates with resistance to miltefosine, a derivative of phosphatidylcholine (PC) that was originally developed as an anticancer drug. In 2021, miltefosine was designated an orphan drug by the United States Food and Drug Administration for the treatment of invasive candidiasis. Importantly, we find that resistance to miltefosine is also caused by mutations in flippases, which control the outward movement of lipids, and that many C. parapsilosis isolates are prone to easily acquiring an increased resistance to miltefosine.

The landscape of submicroscopic structural variants at the OPN1LW/OPN1MW gene cluster on Xq28 underlying blue cone monochromacy

Blue cone monochromacy (BCM) is an X-linked retinal disorder characterized by low vision, photoaversion, and poor color discrimination. BCM is due to the lack of long-wavelength-sensitive and middle-wavelength-sensitive cone photoreceptor function and caused by mutations in the OPN1LW/OPN1MW gene cluster on Xq28. Here, we investigated the prevalence and the landscape of submicroscopic structural variants (SVs) at single-base resolution in BCM patients. We found that about one-third (n = 73) of the 213 molecularly confirmed BCM families carry an SV, most commonly deletions restricted to the OPN1LW/OPN1MW gene cluster. The structure and precise breakpoints of the SVs were resolved in all but one of the 73 families. Twenty-two families-all from the United States-showed the same SV, and we confirmed a common ancestry of this mutation. In total, 42 distinct SVs were identified, including 40 previously unreported SVs, thereby quadrupling the number of precisely mapped SVs underlying BCM. Notably, there was no "region of overlap" among these SVs. However, 90% of SVs encompass the upstream locus control region, an essential enhancer element. Its minimal functional extent based on deletion mapping in patients was refined to 358 bp. Breakpoint analyses suggest diverse mechanisms underlying SV formation as well as in one case the gene conversion-based exchange of a 142-bp deletion between opsin genes. Using parsimonious assumptions, we reconstructed the composition and copy number of the OPN1LW/OPN1MW gene cluster prior to the mutation event and found evidence that large gene arrays may be predisposed to the occurrence of SVs at this locus.

Rett Syndrome and MECP2 Duplication Syndrome: Disorders of MeCP2 Dosage

Rett syndrome (RTT) is a neurodevelopmental disorder caused predominantly by loss-of-function mutations in the gene Methyl-CpG-binding protein 2 (MECP2), which encodes the MeCP2 protein. RTT is a MECP2-related disorder, along with MECP2 duplication syndrome (MDS), caused by gain-of-function duplications of MECP2. Nearly two decades of research have advanced our knowledge of MeCP2 function in health and disease. The following review will discuss MeCP2 protein function and its dysregulation in the MECP2-related disorders RTT and MDS. This will include a discussion of the genetic underpinnings of these disorders, specifically how sporadic X-chromosome mutations arise and manifest in specific populations. We will then review current diagnostic guidelines and clinical manifestations of RTT and MDS. Next, we will delve into MeCP2 biology, describing the dual landscapes of methylated DNA and its reader MeCP2 across the neuronal genome as well as the function of MeCP2 as a transcriptional modulator. Following this, we will outline common MECP2 mutations and genotype-phenotype correlations in both diseases, with particular focus on mutations associated with relatively mild disease in RTT. We will also summarize decades of disease modeling and resulting molecular, synaptic, and behavioral phenotypes associated with RTT and MDS. Finally, we list several therapeutics in the development pipeline for RTT and MDS and available evidence of their safety and efficacy.

Distinct sequence features underlie microdeletions and gross deletions in the human genome

Microdeletions and gross deletions are important causes (~20%) of human inherited disease and their genomic locations are strongly influenced by the local DNA sequence environment. This notwithstanding, no study has systematically examined their underlying generative mechanisms. Here, we obtained 42,098 pathogenic microdeletions and gross deletions from the Human Gene Mutation Database (HGMD) that together form a continuum of germline deletions ranging in size from 1 to 28,394,429 bp. We analyzed the DNA sequence within 1 kb of the breakpoint junctions and found that the frequencies of non‐B DNA‐forming repeats, GC‐content, and the presence of seven of 78 specific sequence motifs in the vicinity of pathogenic deletions correlated with deletion length for deletions of length ≤30 bp. Further, we found that the presence of DR, GQ, and STR repeats is important for the formation of longer deletions (>30 bp) but not for the formation of shorter deletions (≤30 bp) while significantly (χ², p < 2E−16) more microhomologies were identified flanking short deletions than long deletions (length >30 bp). We provide evidence to support a functional distinction between microdeletions and gross deletions. Finally, we propose that a deletion length cut‐off of 25–30 bp may serve as an objective means to functionally distinguish microdeletions from gross deletions.

Molecular and clinical insights into complex genomic rearrangements related to MECP2 duplication syndrome

MECP2 duplication syndrome (MDS) is caused by copy number variation (CNV) spanning the MECP2 gene at Xq28 and is a major cause of intellectual disability (ID) in males. Herein, we describe two unrelated males harboring non-recurrent complex Xq28 rearrangements associated with MDS. Copy number gains were initially detected by quantitative real-time polymerase chain reaction and further delineated by high-resolution array comparative genomic hybridization, familial segregation, expression analysis and X-chromosome inactivation (XCI) evaluation in a carrier mother. SNVs within the rearrangements and/or fluorescent in situ hybridization (FISH) were used to assess the parental origin of the rearrangements. Patient 1 exhibited an intrachromosomal rearrangement, whose structure is consistent with a triplicated segment presumably embedded in an inverted orientation between two duplicated sequences (DUP-TRP/INV-DUP). The rearrangement was inherited from the carrier mother, who exhibits extreme XCI skewing and subtle psychiatric symptoms. Patient 2 presented a de novo (X;Y) unbalanced translocation resulting in duplication of Xq28 and deletion of Yp, originated in the paternal gametogenesis. Neurodevelopmental trajectory and non-neurological symptoms were consistent with previous reports, with the exception of cerebellar vermis hypoplasia in patient 2. Although both patients share the core MDS phenotype, patient 1 showed MECP2 transcript levels in blood similar to controls. Understanding the molecular mechanisms related to MDS is essential for designing targeted therapeutic strategies.

Chromoanagenesis Event Underlies a de novo Pericentric and Multiple Paracentric Inversions in a Single Chromosome Causing Coffin-Siris Syndrome

Chromoanagenesis is a descriptive term that encompasses classes of catastrophic mutagenic processes that generate localized and complex chromosome rearrangements in both somatic and germline genomes. Herein, we describe a 5-year-old female presenting with a constellation of clinical features consistent with a clinical diagnosis of Coffin–Siris syndrome 1 (CSS1). Initial G-banded karyotyping detected a 90-Mb pericentric and a 47-Mb paracentric inversion on a single chromosome. Subsequent analysis of short-read whole-genome sequencing data and genomic optical mapping revealed additional inversions, all clustered on chromosome 6, one of them disrupting ARID1B for which haploinsufficiency leads to the CSS1 disease trait (MIM:135900). The aggregate structural variant data show that the resolved, the resolved derivative chromosome architecture presents four de novo inversions, one pericentric and three paracentric, involving six breakpoint junctions in what appears to be a shuffling of genomic material on this chromosome. Each junction was resolved to nucleotide-level resolution with mutational signatures suggestive of non-homologous end joining. The disruption of the gene ARID1B is shown to occur between the fourth and fifth exon of the canonical transcript with subsequent qPCR studies confirming a decrease in ARID1B expression in the patient versus healthy controls. Deciphering the underlying genomic architecture of chromosomal rearrangements and complex structural variants may require multiple technologies and can be critical to elucidating the molecular etiology of a patient’s clinical phenotype or resolving unsolved Mendelian disease cases.

MECP2-Related Disorders in Males

Methyl CpG binding protein 2 (MECP2) is located at Xq28 and is a multifunctional gene with ubiquitous expression. Loss-of-function mutations in MECP2 are associated with Rett syndrome (RTT), which is a well-characterized disorder that affects mainly females. In boys, however, mutations in MECP2 can generate a wide spectrum of clinical presentations that range from mild intellectual impairment to severe neonatal encephalopathy and premature death. Thus, males can be more difficult to classify and diagnose than classical RTT females. In addition, there are some variants of unknown significance in MECP2, which further complicate the diagnosis of these children. Conversely, the entire duplication of the MECP2 gene is related to MECP2 duplication syndrome (MDS). Unlike in RTT, in MDS, males are predominantly affected. Usually, the duplication is inherited from an apparently asymptomatic carrier mother. Both syndromes share some characteristics, but also differ in some aspects regarding the clinical picture and evolution. In the following review, we present a thorough description of the different types of MECP2 variants and alterations that can be found in males, and explore several genotype-phenotype correlations, although there is still a lot to understand.

Large palindromes on the primate X Chromosome are preserved by natural selection

Mammalian sex chromosomes carry large palindromes that harbor protein-coding gene families with testis-biased expression. However, there are few known examples of sex-chromosome palindromes conserved between species. We identified 26 palindromes on the human X Chromosome, constituting more than 2% of its sequence, and characterized orthologous palindromes in the chimpanzee and the rhesus macaque using a clone-based sequencing approach that incorporates full-length nanopore reads. Many of these palindromes are missing or misassembled in the current reference assemblies of these species' genomes. We find that 12 human X palindromes have been conserved for at least 25 million years, with orthologs in both chimpanzee and rhesus macaque. Insertions and deletions between species are significantly depleted within the X palindromes' protein-coding genes compared to their noncoding sequence, demonstrating that natural selection has preserved these gene families. The spacers that separate the left and right arms of palindromes are a site of localized structural instability, with seven of 12 conserved palindromes showing no spacer orthology between human and rhesus macaque. Analysis of the 1000 Genomes Project data set revealed that human X-palindrome spacers are enriched for deletions relative to arms and flanking sequence, including a common spacer deletion that affects 13% of human X Chromosomes. This work reveals an abundance of conserved palindromes on primate X Chromosomes and suggests that protein-coding gene families in palindromes (most of which remain poorly characterized) promote X-palindrome survival in the face of ongoing structural instability.

Antisense oligonucleotide therapy in a humanized mouse model of MECP2 duplication syndrome

Many intellectual disability disorders are due to copy number variations, and, to date, there have been no treatment options tested for this class of diseases. MECP2 duplication syndrome (MDS) is one of the most common genomic rearrangements in males and results from duplications spanning the methyl-CpG binding protein 2 (MECP2) gene locus. We previously showed that antisense oligonucleotide (ASO) therapy can reduce MeCP2 protein amount in an MDS mouse model and reverse its disease features. This MDS mouse model, however, carried one transgenic human allele and one mouse allele, with the latter being protected from human-specific MECP2-ASO targeting. Because MeCP2 is a dosage-sensitive protein, the ASO must be titrated such that the amount of MeCP2 is not reduced too far, which would cause Rett syndrome. Therefore, we generated an "MECP2 humanized" MDS model that carries two human MECP2 alleles and no mouse endogenous allele. Intracerebroventricular injection of the MECP2-ASO efficiently down-regulated MeCP2 expression throughout the brain in these mice. Moreover, MECP2-ASO mitigated several behavioral deficits and restored expression of selected MeCP2-regulated genes in a dose-dependent manner without any toxicity. Central nervous system administration of MECP2-ASO is therefore well tolerated and beneficial in this mouse model and provides a translatable approach that could be feasible for treating MDS.

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.

Clinical and molecular genetic characterization of two female patients harboring the Xq27.3q28 deletion with different ratios of X chromosome inactivation

A heterozygous deletion at Xq27.3q28 including FMR1, AFF2, and IDS causing intellectual disability and characteristic facial features is very rare in females, with only 10 patients having been reported. Here, we examined two female patients with different clinical features harboring the Xq27.3q28 deletion and determined the chromosomal breakpoints. Moreover, we assessed the X chromosome inactivation (XCI) in peripheral blood from both patients. Both patients had an almost overlapping deletion at Xq27.3q28, however, the more severe patient (Patient 1) showed skewed XCI of the normal X chromosome (79:21) whereas the milder patient (Patient 2) showed random XCI. Therefore, deletion at Xq27.3q28 critically affected brain development, and the ratio of XCI of the normal X chromosome greatly affected the clinical characteristics of patients with deletion at Xq27.3q28. As the chromosomal breakpoints were determined, we analyzed a change in chromatin domains termed topologically associated domains (TADs) using published Hi-C data on the Xq27.3q28 region, and found that only patient 1 had a possibility of a drastic change in TADs. The altered chromatin topologies on the Xq27.3q28 region might affect the clinical features of patient 1 by changing the expression of genes just outside the deletion and/or the XCI establishment during embryogenesis resulting in skewed XCI.

Distinct patterns of complex rearrangements and a mutational signature of microhomeology are frequently observed in PLP1 copy number gain structural variants

BACKGROUND

We investigated the features of the genomic rearrangements in a cohort of 50 male individuals with proteolipid protein 1 (PLP1) copy number gain events who were ascertained with Pelizaeus-Merzbacher disease (PMD; MIM: 312080). We then compared our new data to previous structural variant mutagenesis studies involving the Xq22 region of the human genome. The aggregate data from 159 sequenced join-points (discontinuous sequences in the reference genome that are joined during the rearrangement process) were studied. Analysis of these data from 150 individuals enabled the spectrum and relative distribution of the underlying genomic mutational signatures to be delineated.

METHODS

Genomic rearrangements in PMD individuals with PLP1 copy number gain events were investigated by high-density customized array or clinical chromosomal microarray analysis and breakpoint junction sequence analysis.

RESULTS

High-density customized array showed that the majority of cases (33/50; ~ 66%) present with single duplications, although complex genomic rearrangements (CGRs) are also frequent (17/50; ~ 34%). Breakpoint mapping to nucleotide resolution revealed further previously unknown structural and sequence complexities, even in single duplications. Meta-analysis of all studied rearrangements that occur at the PLP1 locus showed that single duplications were found in ~ 54% of individuals and that, among all CGR cases, triplication flanked by duplications is the most frequent CGR array CGH pattern observed. Importantly, in ~ 32% of join-points, there is evidence for a mutational signature of microhomeology (highly similar yet imperfect sequence matches).

CONCLUSIONS

These data reveal a high frequency of CGRs at the PLP1 locus and support the assertion that replication-based mechanisms are prominent contributors to the formation of CGRs at Xq22. We propose that microhomeology can facilitate template switching, by stabilizing strand annealing of the primer using W-C base complementarity, and is a mutational signature for replicative repair.

Atypical chromosome 22q11.2 deletions are complex rearrangements and have different mechanistic origins

The majority (99%) of individuals with 22q11.2 deletion syndrome (22q11.2DS) have a deletion that is caused by non-allelic homologous recombination between two of four low copy repeat clusters on chromosome 22q11.2 (LCR22s). However, in a small subset of patients, atypical deletions are observed with at least one deletion breakpoint within unique sequence between the LCR22s. The position of the chromosome breakpoints and the mechanisms driving those atypical deletions remain poorly studied. Our large-scale, whole genome sequencing study of >1500 subjects with 22q11.2DS identified six unrelated individuals with atypical deletions of different types. Using a combination of whole genome sequencing data and fiber-fluorescence in situ hybridization, we mapped the rearranged alleles in these subjects. In four of them, the distal breakpoints mapped within one of the LCR22s and we found that the deletions likely occurred by replication-based mechanisms. Interestingly, in two of them, an inversion probably preceded inter-chromosomal 'allelic' homologous recombination between differently oriented LCR22-D alleles. Inversion associated allelic homologous recombination (AHR) may well be a common mechanism driving (atypical) deletions on 22q11.2.

Xq22 deletions and correlation with distinct neurological disease traits in females: Further evidence for a contiguous gene syndrome

Xq22 deletions that encompass PLP1 (Xq22-PLP1-DEL) are notable for variable expressivity of neurological disease traits in females ranging from a mild late-onset form of spastic paraplegia type 2 (MIM# 312920), sometimes associated with skewed X-inactivation, to an early-onset neurological disease trait (EONDT) of severe developmental delay, intellectual disability, and behavioral abnormalities. Size and gene content of Xq22-PLP1-DEL vary and were proposed as potential molecular etiologies underlying variable expressivity in carrier females where two smallest regions of overlap (SROs) were suggested to influence disease. We ascertained a cohort of eight unrelated patients harboring Xq22-PLP1-DEL and performed high-density array comparative genomic hybridization and breakpoint-junction sequencing. Molecular characterization of Xq22-PLP1-DEL from 17 cases (eight herein and nine published) revealed an overrepresentation of breakpoints that reside within repeats (11/17, ~65%) and the clustering of ~47% of proximal breakpoints in a genomic instability hotspot with characteristic non-B DNA density. These findings implicate a potential role for genomic architecture in stimulating the formation of Xq22-PLP1-DEL. The correlation of Xq22-PLP1-DEL gene content with neurological disease trait in female cases enabled refinement of the associated SROs to a single genomic interval containing six genes. Our data support the hypothesis that genes contiguous to PLP1 contribute to EONDT.

Exome Sequencing of a Primary Ovarian Insufficiency Cohort Reveals Common Molecular Etiologies for a Spectrum of Disease

CONTEXT

Primary ovarian insufficiency (POI) encompasses a spectrum of premature menopause, including both primary and secondary amenorrhea. For 75% to 90% of individuals with hypergonadotropic hypogonadism presenting as POI, the molecular etiology is unknown. Common etiologies include chromosomal abnormalities, environmental factors, and congenital disorders affecting ovarian development and function, as well as syndromic and nonsyndromic single gene disorders suggesting POI represents a complex trait.

OBJECTIVE

To characterize the contribution of known disease genes to POI and identify molecular etiologies and biological underpinnings of POI.

DESIGN, SETTING, AND PARTICIPANTS

We applied exome sequencing (ES) and family-based genomics to 42 affected female individuals from 36 unrelated Turkish families, including 31 with reported parental consanguinity.

RESULTS

This analysis identified likely damaging, potentially contributing variants and molecular diagnoses in 16 families (44%), including 11 families with likely damaging variants in known genes and five families with predicted deleterious variants in disease genes (IGSF10, MND1, MRPS22, and SOHLH1) not previously associated with POI. Of the 16 families, 2 (13%) had evidence for potentially pathogenic variants at more than one locus. Absence of heterozygosity consistent with identity-by-descent mediated recessive disease burden contributes to molecular diagnosis in 15 of 16 (94%) families. GeneMatcher allowed identification of additional families from diverse genetic backgrounds.

CONCLUSIONS

ES analysis of a POI cohort further characterized locus heterogeneity, reaffirmed the association of genes integral to meiotic recombination, demonstrated the likely contribution of genes involved in hypothalamic development, and documented multilocus pathogenic variation suggesting the potential for oligogenic inheritance contributing to the development of POI.

Complex structural variants in Mendelian disorders: identification and breakpoint resolution using short- and long-read genome sequencing

Background

Studies have shown that complex structural variants (cxSVs) contribute to human genomic variation and can cause Mendelian disease. We aimed to identify cxSVs relevant to Mendelian disease using short-read whole-genome sequencing (WGS), resolve the precise variant configuration and investigate possible mechanisms of cxSV formation.

Methods

We performed short-read WGS and analysis of breakpoint junctions to identify cxSVs in a cohort of 1324 undiagnosed rare disease patients. Long-read WGS and gene expression analysis were used to resolve one case.

Results

We identified three pathogenic cxSVs: a de novo duplication-inversion-inversion-deletion affecting ARID1B, a de novo deletion-inversion-duplication affecting HNRNPU and a homozygous deletion-inversion-deletion affecting CEP78. Additionally, a de novo duplication-inversion-duplication overlapping CDKL5 was resolved by long-read WGS demonstrating the presence of both a disrupted and an intact copy of CDKL5 on the same allele, and gene expression analysis showed both parental alleles of CDKL5 were expressed. Breakpoint analysis in all the cxSVs revealed both microhomology and longer repetitive elements.

Conclusions

Our results corroborate that cxSVs cause Mendelian disease, and we recommend their consideration during clinical investigations. We show that resolution of breakpoints can be critical to interpret pathogenicity and present evidence of replication-based mechanisms in cxSV formation.

Electronic supplementary material

The online version of this article (10.1186/s13073-018-0606-6) contains supplementary material, which is available to authorized users.

A Multicentric Brazilian Investigative Study of Copy Number Variations in Patients with Congenital Anomalies and Intellectual Disability

Genomic imbalances are the most common cause of congenital anomalies (CA) and intellectual disability (ID). The aims of this study were to identify copy number variations (CNVs) in 416 patients with CA and ID from 5 different genetics centers within 4 different states by using the Multiplex Ligation-dependent Probe Amplification (MLPA) technique and to apply the chromosomal microarray (CMA) methodology in selected cases. The samples were analyzed by MLPA kits P064, P036, P070 and P250. Positive results were found in 97/416 (23.3%) patients. CMA was applied in 14 selected cases. In 6/14 (42.85%) patients, CMA detected other copy number variations not detected by the MLPA studies. Although CMA is indispensable for genotype refinement, the technique is still unfeasible in some countries as a routine analysis due to economic and technical limitations. In these cases, clinical evaluation followed by karyotyping and MLPA analysis is a helpful and affordable solution for diagnostic purposes.

Genomic disorders 20 years on-mechanisms for clinical manifestations

Genomic disorders result from copy-number variants (CNVs) or submicroscopic rearrangements of the genome rather than from single nucleotide variants (SNVs). Diverse technologies, including array comparative genomic hybridization (aCGH) and single nucleotide polymorphism (SNP) microarrays, and more recently, whole genome sequencing and whole-exome sequencing, have enabled robust genome-wide unbiased detection of CNVs in affected individuals and in reportedly healthy controls. Sequencing of breakpoint junctions has allowed for elucidation of upstream mechanisms leading to genomic instability and resultant structural variation, whereas studies of the association between CNVs and specific diseases or susceptibility to morbid traits have enhanced our understanding of the downstream effects. In this review, we discuss the hallmarks of genomic disorders as they were defined during the first decade of the field, including genomic instability and the mechanism for rearrangement defined as nonallelic homologous recombination (NAHR); recurrent vs nonrecurrent rearrangements; and gene dosage sensitivity. Moreover, we highlight the exciting advances of the second decade of this field, including a deeper understanding of genomic instability and the mechanisms underlying complex rearrangements, mechanisms for constitutional and somatic chromosomal rearrangements, structural intra-species polymorphisms and susceptibility to NAHR, the role of CNVs in the context of genome-wide copy number and single nucleotide variation, and the contribution of noncoding CNVs to human disease.

Nanopore sequencing of complex genomic rearrangements in yeast reveals mechanisms of repeat-mediated double-strand break repair

Improper DNA double-strand break (DSB) repair results in complex genomic rearrangements (CGRs) in many cancers and various congenital disorders in humans. Trinucleotide repeat sequences, such as (GAA)_n repeats in Friedreich's ataxia, (CTG)_n repeats in myotonic dystrophy, and (CGG)_n repeats in fragile X syndrome, are also subject to double-strand breaks within the repetitive tract followed by DNA repair. Mapping the outcomes of CGRs is important for understanding their causes and potential phenotypic effects. However, high-resolution mapping of CGRs has traditionally been a laborious and highly skilled process. Recent advances in long-read DNA sequencing technologies, specifically Nanopore sequencing, have made possible the rapid identification of CGRs with single base pair resolution. Here, we have used whole-genome Nanopore sequencing to characterize several CGRs that originated from naturally occurring DSBs at (GAA)_n microsatellites in Saccharomyces cerevisiae. These data gave us important insights into the mechanisms of DSB repair leading to CGRs.

[MECP2 duplication syndrome: a clinical analysis of three cases and literature review]

MECP2 duplication syndrome (MDS) is a rare pediatric disease and mainly manifests as delayed motor development, language loss or delay, recurrent infection, severe intellectual disability, epilepsy, autistic symptoms, and early infantile hypotonia. In this article, the three children with this disease were all boys. Cases 1 and 2 had delayed motor development, and language loss or delay as initial manifestations, and case 3 had recurrent infection as initial manifestation. Physical examination showed hypotonia and negative pathological signs in each case. Case 1 had tonic-clonic seizures and electroencephalography showed focal seizures, for which he was given oxcarbazepine, levetiracetam, and clonazepam as the antiepileptic treatment to control seizures. Case 3 experienced one absence seizure and three head-nodding seizures with normal electroencephalographic findings during these seizures, and therefore, he was not given antiepileptic treatment. In each case, recurrent infection was improved with the increase in age, but there were no significant improvements in language or intelligence. Array-based comparative genomic hybridization (aCGH) showed MECP2 duplication in X chromosome in each case, and so they were diagnosed with MDS. MDS should be considered for children with delayed development complicated by recurrent infection and epileptic seizures, and early aCGH helps with the diagnosis of this disease.

A clear bias in parental origin of de novo pathogenic CNVs related to intellectual disability, developmental delay and multiple congenital anomalies

Copy number variation (CNV) is of great significance in human evolution and disorders. Through tracing the parent-of-origin of de novo pathogenic CNVs, we are expected to investigate the relative contributions of germline genomic stability on reproductive health. In our study, short tandem repeat (STR) and single nucleotide polymorphism (SNP) were used to determine the parent-of-origin of 87 de novo pathogenic CNVs found in unrelated patients with intellectual disability (ID), developmental delay (DD) and multiple congenital anomalies (MCA). The results shown that there was a significant difference on the distribution of the parent-of-origin for different CNVs types (Chi-square test, p = 4.914 × 10⁻³). An apparently paternal bias existed in deletion CNVs and a maternal bias in duplication CNVs, indicating that the relative contribution of paternal germline variations is greater than that of maternal to the origin of deletions, and vice versa to the origin of duplications. By analyzing the sequences flanking the breakpoints, we also confirmed that non-allelic homologous recombination (NAHR) served as the major mechanism for the formation of recurrent CNVs whereas non-SDs-based mechanisms played a part in generating rare non-recurrent CNVs and might relate to the paternal germline bias in deletion CNVs.

A sibship with duplication of Xq28 inherited from the mother; genomic characterization and clinical outcomes

BACKGROUND

Loss-of-function mutations in methyl-CpG-binding protein 2 (MECP2; MIM *300005) results in the Rett syndrome, whereas gain-of-function mutations are associated with the MECP2 duplication syndrome.

METHODS

We did research on a family with two brothers showing Xq28 duplication syndrome using various molecular cytogenetic techniques such as multiplex ligation-dependent probe amplification and array-based genomic hybridization.

RESULTS

The duplicated region had several genes including MECP2 and interleukin-1 receptor associated kinase 1 (IRAK1; MIM *300283). MECP2 and IRAK1 were associated with the neurological phenotypes in dose-sensitive and dose-critical manner. The brothers demonstrated severe intellectual disability, autistic features, generalized hypotonia, recurrent infections, epilepsy, choreiform movements such as hand-wringing movement, and moderate increased spasticity with the lower limbs. The X-inactivation test showed a complete skewed X inactivation pattern of mother. In this reason, the mother had the same loci duplication but showed significantly little neurological manifestation compared to the two sons.

CONCLUSIONS

MECP2/IRAK1 duplication at Xq28 is inherited as an X-linked recessive trait and male-specific disorder associated with severe intellectual disability. We tried to analyze the information of the relationship between neuropsychiatric phenotype and the extent of duplication at Xq28 by comparing with previous reports.

Defining the diverse spectrum of inversions, complex structural variation, and chromothripsis in the morbid human genome

Background

Structural variation (SV) influences genome organization and contributes to human disease. However, the complete mutational spectrum of SV has not been routinely captured in disease association studies.

Results

We sequenced 689 participants with autism spectrum disorder (ASD) and other developmental abnormalities to construct a genome-wide map of large SV. Using long-insert jumping libraries at 105X mean physical coverage and linked-read whole-genome sequencing from 10X Genomics, we document seven major SV classes at ~5 kb SV resolution. Our results encompass 11,735 distinct large SV sites, 38.1% of which are novel and 16.8% of which are balanced or complex. We characterize 16 recurrent subclasses of complex SV (cxSV), revealing that: (1) cxSV are larger and rarer than canonical SV; (2) each genome harbors 14 large cxSV on average; (3) 84.4% of large cxSVs involve inversion; and (4) most large cxSV (93.8%) have not been delineated in previous studies. Rare SVs are more likely to disrupt coding and regulatory non-coding loci, particularly when truncating constrained and disease-associated genes. We also identify multiple cases of catastrophic chromosomal rearrangements known as chromoanagenesis, including somatic chromoanasynthesis, and extreme balanced germline chromothripsis events involving up to 65 breakpoints and 60.6 Mb across four chromosomes, further defining rare categories of extreme cxSV.

Conclusions

These data provide a foundational map of large SV in the morbid human genome and demonstrate a previously underappreciated abundance and diversity of cxSV that should be considered in genomic studies of human disease.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1158-6) contains supplementary material, which is available to authorized users.

Mechanisms for Complex Chromosomal Insertions

Chromosomal insertions are genomic rearrangements with a chromosome segment inserted into a non-homologous chromosome or a non-adjacent locus on the same chromosome or the other homologue, constituting ~2% of nonrecurrent copy-number gains. Little is known about the molecular mechanisms of their formation. We identified 16 individuals with complex insertions among 56,000 individuals tested at Baylor Genetics using clinical array comparative genomic hybridization (aCGH) and fluorescence in situ hybridization (FISH). Custom high-density aCGH was performed on 10 individuals with available DNA, and breakpoint junctions were fine-mapped at nucleotide resolution by long-range PCR and DNA sequencing in 6 individuals to glean insights into potential mechanisms of formation. We observed microhomologies and templated insertions at the breakpoint junctions, resembling the breakpoint junction signatures found in complex genomic rearrangements generated by replication-based mechanism(s) with iterative template switches. In addition, we analyzed 5 families with apparently balanced insertion in one parent detected by FISH analysis and found that 3 parents had additional small copy-number variants (CNVs) at one or both sides of the inserting fragments as well as at the inserted sites. We propose that replicative repair can result in interchromosomal complex insertions generated through chromothripsis-like chromoanasynthesis involving two or three chromosomes, and cause a significant fraction of apparently balanced insertions harboring small flanking CNVs.

By traditional cytogenetic techniques, the incidence of microscopically visible chromosomal insertions was estimated to be 1 in 80,000 live births. More recently, by aCGH in conjunction with FISH confirmation of the aCGH findings, insertion events were demonstrated to occur much more frequently (1 in ~500 individuals tested). Although frequently detected, little is known about the molecular mechanisms of their formation. In this study, we identified 16 individuals with complex chromosomal insertions among 56,000 individuals tested at Baylor Genetics using clinical microarray analysis (CMA) and FISH. Custom high-density aCGH was performed on 10 individuals with available DNA, and breakpoint junctions were fine-mapped at nucleotide resolution by long-range PCR and DNA sequencing in 6 individuals to glean insights into potential mechanisms of formation. In addition, we analyzed 5 families with apparently balanced insertion in one parent detected by FISH analysis and found that 3 parents had additional small copy-number variants (CNVs) at one or both sides of the inserting fragments as well as at the inserted sites. We propose that replicative repair can result in interchromosomal complex insertions generated through chromothripsis-like chromoanasynthesis involving two or three chromosomes, and cause a significant fraction of apparently balanced insertions harboring small flanking CNVs.

Major influence of repetitive elements on disease-associated copy number variants (CNVs)

Copy number variants (CNVs) are important contributors to the human pathogenic genetic diversity as demonstrated by a number of cases reported in the literature. The high homology between repetitive elements may guide genomic stability which will give rise to CNVs either by non-allelic homologous recombination (NAHR) or non-homologous end joining (NHEJ). Here, we present a short guide based on previously documented cases of disease-associated CNVs in order to provide a general view on the impact of repeated elements on the stability of the genomic sequence and consequently in the origin of the human pathogenic variome.

Copy-Number Variation Contributes to the Mutational Load of Bardet-Biedl Syndrome

Bardet-Biedl syndrome (BBS) is a defining ciliopathy, notable for extensive allelic and genetic heterogeneity, almost all of which has been identified through sequencing. Recent data have suggested that copy-number variants (CNVs) also contribute to BBS. We used a custom oligonucleotide array comparative genomic hybridization (aCGH) covering 20 genes that encode intraflagellar transport (IFT) components and 74 ciliopathy loci to screen 92 unrelated individuals with BBS, irrespective of their known mutational burden. We identified 17 individuals with exon-disruptive CNVs (18.5%), including 13 different deletions in eight BBS genes (BBS1, BBS2, ARL6/BBS3, BBS4, BBS5, BBS7, BBS9, and NPHP1) and a deletion and a duplication in other ciliopathy-associated genes (ALMS1 and NPHP4, respectively). By contrast, we found a single heterozygous exon-disruptive event in a BBS-associated gene (BBS9) in 229 control subjects. Superimposing these data with resequencing revealed CNVs to (1) be sufficient to cause disease, (2) Mendelize heterozygous deleterious alleles, and (3) contribute oligogenic alleles by combining point mutations and exonic CNVs in multiple genes. Finally, we report a deletion and a splice site mutation in IFT74, inherited under a recessive paradigm, defining a candidate BBS locus. Our data suggest that CNVs contribute pathogenic alleles to a substantial fraction of BBS-affected individuals and highlight how either deletions or point mutations in discrete splice isoforms can induce hypomorphic mutations in genes otherwise intolerant to deleterious variation. Our data also suggest that CNV analyses and resequencing studies unbiased for previous mutational burden is necessary to delineate the complexity of disease architecture.

Comparison of three whole genome amplification methods for detection of genomic aberrations in single cells

OBJECTIVE

Detection of genomic copy number abnormalities in a single cell using array comparative genomic hybridization (CGH) offers a promising non-invasive alternative for prenatal diagnosis. Our objective was to compare three commercially available whole-genome amplification (WGA) kits for their capacity to produce high quality DNA from single cells that is suitable for both molecular genotyping and array CGH.

METHODS

We examined kit performance on unfixed, fixed and fixed/permeabilized lymphoblastoid cells. Molecular genotyping methods were used to evaluate the fidelity of amplified DNA for genomic profiling, while array CGH was used to assess copy number from single cells harboring trisomy 21, a DiGeorge syndrome deletion, a CMT1A duplication or a MECP2 duplication.

RESULTS

Molecular genotyping was achieved from single cells but performance varied between WGA kits. Furthermore, we consistently detected a dosage difference in sex chromosomes for gender mismatched hybridizations and for chromosome 21 in trisomy 21 cells. The 2.5 Mb DiGeorge syndrome deletion was also detected using all three WGA platforms, whereas the 1.3 Mb CMT1A and the 0.6 Mb MECP2 duplications were not consistently detected.

CONCLUSION

These data suggest that single cell molecular genotyping and copy number analysis can be accomplished when WGA conditions are optimized. © 2016 John Wiley & Sons, Ltd.

PMS2 inactivation by a complex rearrangement involving an HERV retroelement and the inverted 100-kb duplicon on 7p22.1

Biallelic PMS2 mutations are responsible for more than half of all cases of constitutional mismatch repair deficiency (CMMRD), a recessively inherited childhood cancer predisposition syndrome. The mismatch repair gene PMS2 is partly embedded within one copy of an inverted 100-kb low-copy repeat (LCR) on 7p22.1. In an individual with CMMRD syndrome, PMS2 was found to be homozygously inactivated by a complex chromosomal rearrangement, which separates the 5'-part from the 3'-part of the gene. The rearrangement involves sequences of the inverted 100-kb LCR and a human endogenous retrovirus element and may be associated with an inversion that is indistinguishable from the known inversion polymorphism affecting the ~0.7-Mb sequence intervening the LCR. Its formation is best explained by a replication-based mechanism (RBM) such as fork stalling and template switching/microhomology-mediated break-induced replication (FoSTeS/MMBIR). This finding supports the hypothesis that the inverted LCR can not only facilitate the formation of the non-allelic homologous recombination-mediated inversion polymorphism but it also promotes the occurrence of more complex rearrangements that can be associated with a large inversion, as well, but are mediated by a RBM. This further suggests that among the inversion polymorphism on 7p22.1, more complex rearrangements might be hidden. Furthermore, as the locus is embedded in a common fragile site (CFS) region, this rearrangement also supports the recently raised hypothesis that CFS sequence motifs may facilitate replication-based rearrangement mechanisms.

Chromosome Xq28 duplication encompassing MECP2: Clinical and molecular analysis of 16 new patients from 10 families in China

INTRODUCTION

Chromosome Xq28 duplications encompassing methyl-CpG-binding protein 2 gene (MECP2) are observed most in males with a severe neurodevelopmental disorder associated with hypotonia, spasticity, severe learning disability, delayed psychomotor development, and recurrent pulmonary infections. Most female carriers are asymptomatic due to extremely or completely skewed X-inactivation.

METHODS

A retrospective clinical and molecular study was conducted to examine 16 patients and two fetuses from 10 families who were identified among patients with Xq28 duplications who presented at genetic clinics.

RESULTS

Of all 16 patients, 10 had a family history. Only one patient was female. All of the patients had no relevant pre-natal history. All of the patients exhibited severe psychomotor developmental delay, infantile hypotonia and recurrent infections. Some of the patients exhibited cardiac abnormalities, gastrointestinal mobility problems, hydrocele of tunica vaginalis, cryptorchidism, and autistic phenotypes. Additionally, neonatal kidney calculus, premature closure of the fontanel and pulmonary sequestration were found in the patients. Duplication sizes in these patients range from 0.21 to 14.391 Mb (most were smaller than 1 Mb), and all the duplications included host cell factor C1 (HCFC1), interleukin-1 receptor-associated kinase 1 (IRAK1), and MECP2. Bioinformatics analysis revealed that approximately half of the distal breakpoints were located within the low-copy repeats (LCRs), which may be involved in the recombination. The two fetuses were found to be healthy in the prenatal diagnosis.

CONCLUSION

This is the first large cohort of patients with MECP2 duplication syndrome, including a female, reported in China. Interestingly, neonatal kidney calculus, premature closure of the fontanel and pulmonary sequestration were first reported in this syndrome. However, it was difficult to distinguish if these patients represented unique cases or if these phenotypes can be considered as part of the syndrome. The correlation between the infrequent phenotypes and duplications/genes in the duplication region needs further systematic delineation. In conclusion, our study suggested that it is important to emphasize molecular genetic analysis in patients with developmental delay/intellectual disability and recurrent infections and that it is especially important for familial female carriers to accept prenatal diagnosis.

MECP2 Duplication Syndrome: Evidence of Enhanced Oxidative Stress. A Comparison with Rett Syndrome

Rett syndrome (RTT) and MECP2 duplication syndrome (MDS) are neurodevelopmental disorders caused by alterations in the methyl-CpG binding protein 2 (MECP2) gene expression. A relationship between MECP2 loss-of-function mutations and oxidative stress has been previously documented in RTT patients and murine models. To date, no data on oxidative stress have been reported for the MECP2 gain-of-function mutations in patients with MDS. In the present work, the pro-oxidant status and oxidative fatty acid damage in MDS was investigated (subjects n = 6) and compared to RTT (subjects n = 24) and healthy condition (subjects n = 12). Patients with MECP2 gain-of-function mutations showed increased oxidative stress marker levels (plasma non-protein bound iron, intraerythrocyte non-protein bound iron, F₂-isoprostanes, and F₄-neuroprostanes), as compared to healthy controls (P ≤ 0.05). Such increases were similar to those observed in RTT patients except for higher plasma F₂-isoprostanes levels (P < 0.0196). Moreover, plasma levels of F₂-isoprostanes were significantly correlated (P = 0.0098) with the size of the amplified region. The present work shows unique data in patients affected by MDS. For the first time MECP2 gain-of-function mutations are indicated to be linked to an oxidative damage and related clinical symptoms overlapping with those of MECP2 loss-of-function mutations. A finely tuned balance of MECP2 expression appears to be critical to oxidative stress homeostasis, thus shedding light on the relevance of the redox balance in the central nervous system integrity.

Mechanisms underlying structural variant formation in genomic disorders

With the recent burst of technological developments in genomics, and the clinical implementation of genome-wide assays, our understanding of the molecular basis of genomic disorders, specifically the contribution of structural variation to disease burden, is evolving quickly. Ongoing studies have revealed a ubiquitous role for genome architecture in the formation of structural variants at a given locus, both in DNA recombination-based processes and in replication-based processes. These reports showcase the influence of repeat sequences on genomic stability and structural variant complexity and also highlight the tremendous plasticity and dynamic nature of our genome in evolution, health and disease susceptibility.

Altered neuronal network and rescue in a human MECP2 duplication model

Increased dosage of methyl-CpG-binding protein-2 (MeCP2) results in a dramatic neurodevelopmental phenotype with onset at birth. We generated induced pluripotent stem cells (iPSCs) from patients with the MECP2 duplication syndrome (MECP2dup), carrying different duplication sizes, to study the impact of increased MeCP2 dosage in human neurons. We show that cortical neurons derived from these different MECP2dup iPSC lines have increased synaptogenesis and dendritic complexity. In addition, using multi-electrodes arrays, we show that neuronal network synchronization was altered in MECP2dup-derived neurons. Given MeCP2 functions at the epigenetic level, we tested whether these alterations were reversible using a library of compounds with defined activity on epigenetic pathways. One histone deacetylase inhibitor, NCH-51, was validated as a potential clinical candidate. Interestingly, this compound has never been considered before as a therapeutic alternative for neurological disorders. Our model recapitulates early stages of the human MECP2 duplication syndrome and represents a promising cellular tool to facilitate therapeutic drug screening for severe neurodevelopmental disorders.