Mechanisms for nonrecurrent genomic rearrangements associated with CMT1A or HNPP: rare CNVs as a cause for missing heritability

A study concept of expeditious clinical enrollment for genetic modifier studies in Charcot-Marie-Tooth neuropathy 1A

BACKGROUND

Caused by duplications of the gene encoding peripheral myelin protein 22 (PMP22), Charcot-Marie-Tooth disease type 1A (CMT1A) is the most common hereditary neuropathy. Despite this shared genetic origin, there is considerable variability in clinical severity. It is hypothesized that genetic modifiers contribute to this heterogeneity, the identification of which may reveal novel therapeutic targets. In this study, we present a comprehensive analysis of clinical examination results from 1564 CMT1A patients sourced from a prospective natural history study conducted by the RDCRN-INC (Inherited Neuropathy Consortium). Our primary objective is to delineate extreme phenotype profiles (mild and severe) within this patient cohort, thereby enhancing our ability to detect genetic modifiers with large effects.

METHODS

We have conducted large-scale statistical analyses of the RDCRN-INC database to characterize CMT1A severity across multiple metrics.

RESULTS

We defined patients below the 10th (mild) and above the 90th (severe) percentiles of age-normalized disease severity based on the CMT Examination Score V2 and foot dorsiflexion strength (MRC scale). Based on extreme phenotype categories, we defined a statistically justified recruitment strategy, which we propose to use in future modifier studies.

INTERPRETATION

Leveraging whole genome sequencing with base pair resolution, a future genetic modifier evaluation will include single nucleotide association, gene burden tests, and structural variant analysis. The present work not only provides insight into the severity and course of CMT1A, but also elucidates the statistical foundation and practical considerations for a cost-efficient and straightforward patient enrollment strategy that we intend to conduct on additional patients recruited globally.

Functional categorization of gene regulatory variants that cause Mendelian conditions

Much of our current understanding of rare human diseases is driven by coding genetic variants. However, non-coding genetic variants play a pivotal role in numerous rare human diseases, resulting in diverse functional impacts ranging from altered gene regulation, splicing, and/or transcript stability. With the increasing use of genome sequencing in clinical practice, it is paramount to have a clear framework for understanding how non-coding genetic variants cause disease. To this end, we have synthesized the literature on hundreds of non-coding genetic variants that cause rare Mendelian conditions via the disruption of gene regulatory patterns and propose a functional classification system. Specifically, we have adapted the functional classification framework used for coding variants (i.e., loss-of-function, gain-of-function, and dominant-negative) to account for features unique to non-coding gene regulatory variants. We identify that non-coding gene regulatory variants can be split into three distinct categories by functional impact: (1) non-modular loss-of-expression (LOE) variants; (2) modular loss-of-expression (mLOE) variants; and (3) gain-of-ectopic-expression (GOE) variants. Whereas LOE variants have a direct corollary with coding loss-of-function variants, mLOE and GOE variants represent disease mechanisms that are largely unique to non-coding variants. These functional classifications aim to provide a unified terminology for categorizing the functional impact of non-coding variants that disrupt gene regulatory patterns in Mendelian conditions.

Microsatellite break-induced replication generates highly mutagenized extrachromosomal circular DNAs

Extrachromosomal circular DNAs (eccDNAs) are produced from all regions of the eucaryotic genome. In tumors, highly transcribed eccDNAs have been implicated in oncogenesis, neoantigen production and resistance to chemotherapy. Here we show that unstable microsatellites capable of forming hairpin, triplex, quadruplex and AT-rich structures generate eccDNAs when integrated at a common ectopic site in human cells. These non-B DNA prone microsatellites form eccDNAs by replication-dependent mechanisms. The microsatellite-based eccDNAs are highly mutagenized and display template switches to sister chromatids and to nonallelic chromosomal sites. High frequency mutagenesis occurs within the eccDNA microsatellites and extends bidirectionally for several kilobases into flanking DNA and nonallelic DNA. Mutations include mismatches, short duplications, longer nontemplated insertions and large deletions. Template switching leads to recurrent deletions and recombination domains within the eccDNAs. Template switching events are microhomology-mediated, but do not occur at all potential sites of complementarity. Each microsatellite exhibits a distinct pattern of recombination, microhomology choice and base substitution signature. Depletion of Rad51, the COPS2 signalosome subunit or POLη alter the eccDNA mutagenic profiles. We propose an asynchronous capture model based on break-induced replication from microsatellite-induced DNA breaks for the generation and circularization of mutagenized eccDNAs and genomic homologous recombination deficiency (HRD) scars.

Rare PMP22 variants in mild to severe neuropathy uncorrelated to plasma GDF15 or neurofilament light

Charcot-Marie-Tooth disease (CMT) is a heterogeneous set of hereditary neuropathies whose genetic causes are not fully understood. Here, we characterize three previously unknown variants in PMP22 and assess their effect on the recently described potential CMT biomarkers' growth differentiation factor 15 (GDF15) and neurofilament light (NFL): first, a heterozygous PMP22 c.178G > A (p.Glu60Lys) in one mother-son pair with adult-onset mild axonal neuropathy. The variant led to abnormal splicing, confirmed in fibroblasts by reverse transcription PCR. Second, a de novo PMP22 c.35A > C (p.His12Pro), and third, a heterozygous 3.2 kb deletion predicting loss of exon 4. The latter two had severe CMT and ultrasonography showing strong nerve enlargement similar to a previous case of exon 4 loss due to a larger deletion. We further studied patients with PMP22 duplication (CMT1A) finding slightly elevated plasma NFL, as measured by the single molecule array immunoassay (SIMOA). In addition, plasma GDF15, as measured by ELISA, correlated with symptom severity for CMT1A. However, in the severely affected individuals with PMP22 exon 4 deletion or p.His12Pro, these biomarkers were within the range of variability of CMT1A and controls, although they had more pronounced nerve hypertrophy. This study adds p.His12Pro and confirms PMP22 exon 4 deletion as causes of severe CMT, whereas the previously unknown splice variant p.Glu60Lys leads to mild axonal neuropathy. Our results suggest that GDF15 and NFL do not distinguish CMT1A from advanced hypertrophic neuropathy caused by rare PMP22 variants.

Tibial neuropathy, a rare manifestation of hereditary neuropathy with liability to pressure palsy: A case report

Hereditary neuropathy with liability to pressure palsy (HNPP) is characterized by acute, painless and recurrent mononeuropathies. Genetic testing shows PMP22 gene deletion of chromosome 17p11.2 can provide evidence for the diagnosis of HNPP. Reports on tibial neuropathy as the main manifestation of HNPP are very rare. We report a 14-year-old girl who was admitted to our hospital due to plantar foot numbness and plantar flexion weakness of her left foot. The patient had a history of lateral dorsal numbness and right foot drop when she was 3 years old. Clinical symptoms, and neurological examination demonstrated tibial neuropathy. Electromyography showed extensive peripheral nerve, including median nerve, ulnar nerve, tibial nerve and peroneal nerve, were involved. The diagnosis of HNPP was confirmed by genetic testing which disclosed a deletion of PMP22 gene. She was completely asymptomatic in one month after neurotrophic drug treatments.

Novel gene-intergenic fusion involving ubiquitin E3 ligase UBE3C causes distal hereditary motor neuropathy

Distal hereditary motor neuropathies (dHMNs) are a group of inherited diseases involving the progressive, length-dependent axonal degeneration of the lower motor neurons. There are currently 29 reported causative genes and four disease loci implicated in dHMN. Despite the high genetic heterogeneity, mutations in the known genes account for less than 20% of dHMN cases, with the mutations identified predominantly being point mutations or indels. We have expanded the spectrum of dHMN mutations with the identification of a 1.35 Mb complex structural variation (SV) causing a form of autosomal dominant dHMN (DHMN1 OMIM %182906). Given the complex nature of SV mutations and the importance of studying pathogenic mechanisms in a neuronal setting, we generated a patient-derived DHMN1 motor neuron model harbouring the 1.35 Mb complex insertion. The DHMN1 complex insertion creates a duplicated copy of the first 10 exons of the ubiquitin-protein E3 ligase gene (UBE3C) and forms a novel gene-intergenic fusion sense transcript by incorporating a terminal pseudo-exon from intergenic sequence within the DHMN1 locus. The UBE3C intergenic fusion (UBE3C-IF) transcript does not undergo nonsense-mediated decay and results in a significant reduction of wild-type full-length UBE3C (UBE3C-WT) protein levels in DHMN1 iPSC-derived motor neurons. An engineered transgenic Caenorhabditis elegans model expressing the UBE3C-IF transcript in GABA-ergic motor neurons shows neuronal synaptic transmission deficits. Furthermore, the transgenic animals are susceptible to heat stress, which may implicate defective protein homeostasis underlying DHMN1 pathogenesis. Identification of the novel UBE3C-IF gene-intergenic fusion transcript in motor neurons highlights a potential new disease mechanism underlying axonal and motor neuron degeneration. These complementary models serve as a powerful paradigm for studying the DHMN1 complex SV and an invaluable tool for defining therapeutic targets for DHMN1.

Developmental genomics of limb malformations: Allelic series in association with gene dosage effects contribute to the clinical variability

Genetic heterogeneity, reduced penetrance, and variable expressivity, the latter including asymmetric body axis plane presentations, have all been described in families with congenital limb malformations (CLMs). Interfamilial and intrafamilial heterogeneity highlight the complexity of the underlying genetic pathogenesis of these developmental anomalies. Family-based genomics by exome sequencing (ES) and rare variant analyses combined with whole-genome array-based comparative genomic hybridization were implemented to investigate 18 families with limb birth defects. Eleven of 18 (61%) families revealed explanatory variants, including 7 single-nucleotide variant alleles and 3 copy number variants (CNVs), at previously reported "disease trait associated loci": BHLHA9, GLI3, HOXD cluster, HOXD13, NPR2, and WNT10B. Breakpoint junction analyses for all three CNV alleles revealed mutational signatures consistent with microhomology-mediated break-induced replication, a mechanism facilitated by Alu/Alu-mediated rearrangement. Homozygous duplication of BHLHA9 was observed in one Turkish kindred and represents a novel contributory genetic mechanism to Gollop-Wolfgang Complex (MIM: 228250), where triplication of the locus has been reported in one family from Japan (i.e., 4n = 2n + 2n versus 4n = 3n + 1n allelic configurations). Genes acting on limb patterning are sensitive to a gene dosage effect and are often associated with an allelic series. We extend an allele-specific gene dosage model to potentially assist, in an adjuvant way, interpretations of interconnections among an allelic series, clinical severity, and reduced penetrance of the BHLHA9-related CLM spectrum.

Biology in balance: human diploid genome integrity, gene dosage, and genomic medicine

The path to completion of the functional annotation of the haploid human genome reference build, exploration of the clan genomics hypothesis, understanding human gene and genome functional biology, and gene genome and organismal evolution, is in reach.

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

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

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.

Pmp22 super-enhancer deletion causes tomacula formation and conduction block in peripheral nerves

Copy number variation of the peripheral nerve myelin gene Peripheral Myelin Protein 22 (PMP22) causes multiple forms of inherited peripheral neuropathy. The duplication of a 1.4 Mb segment surrounding this gene in chromosome 17p12 (c17p12) causes the most common form of Charcot-Marie-Tooth disease type 1A, whereas the reciprocal deletion of this gene causes a separate neuropathy termed hereditary neuropathy with liability to pressure palsies (HNPP). PMP22 is robustly induced in Schwann cells in early postnatal development, and several transcription factors and their cognate regulatory elements have been implicated in coordinating the gene's proper expression. We previously found that a distal super-enhancer domain was important for Pmp22 expression in vitro, with particular impact on a Schwann cell-specific alternative promoter. Here, we investigate the consequences of deleting this super-enhancer in vivo. We find that loss of the super-enhancer in mice reduces Pmp22 expression throughout development and into adulthood, with greater impact on the Schwann cell-specific promoter. Additionally, these mice display tomacula formed by excessive myelin folding, a pathological hallmark of HNPP, as have been previously observed in heterozygous Pmp22 mice as well as sural biopsies from patients with HNPP. Our findings demonstrate a mechanism by which smaller copy number variations, not including the Pmp22 gene, are sufficient to reduce gene expression and phenocopy a peripheral neuropathy caused by the HNPP-associated deletion encompassing PMP22.

Paternal gender specificity and mild phenotypes in Charcot-Marie-Tooth type 1A patients with de novo 17p12 rearrangements

Background

Charcot–Marie–Tooth disease type 1A (CMT1A) and hereditary neuropathy with liability to pressure palsies (HNPP) are developed by duplication and deletion of the 17p12 (PMP22) region, respectively.

Methods

De novo rates were determined in 211 CMT1A or HNPP trio families, and then, analyzed gender‐specific genetic features and clinical phenotypes of the de novo cases.

Results

This study identified 40 de novo cases (19.0%). Paternal origin was highly frequent compared to maternal origin (p = .005). Most de novo CMT1A rearrangements occurred between non‐sister chromatids (p = .003), but it was interesting that three of the four sister chromatids exchange cases were observed in the less frequent maternal origin. Paternal ages at the affected child births were slightly higher in the de novo CMT1A group than in the non‐de novo CMT1A control group (p = .0004). For the disability score of CMTNS, the de novo CMT1A group had a slightly lower value compared to the control group (p = .005). Electrophysiological studies showed no significant differences between the two groups.

Conclusion

This study suggests that de novo CMT1A patients tend to have milder symptoms and that the paternal ages at child births in the de novo group are higher than those of the non‐de novo group.

Variation in SIPA1L2 is correlated with phenotype modification in Charcot- Marie- Tooth disease type 1A

OBJECTIVE

Genetic modifiers in rare disease have long been suspected to contribute to the considerable variance in disease expression, including Charcot-Marie-Tooth disease type 1A (CMT1A). To address this question, the Inherited Neuropathy Consortium collected a large standardized sample of such rare CMT1A patients over a period of 8 years. CMT1A is caused in most patients by a uniformly sized 1.5 Mb duplication event involving the gene PMP22.

METHODS

We genotyped DNA samples from 971 CMT1A patients on Illumina BeadChips. Genome-wide analysis was performed in a subset of 330 of these patients, who expressed the extremes of a hallmark symptom: mild and severe foot dorsiflexion strength impairment. SIPA1L2 (signal-induced proliferation-associated 1 like 2), the top identified candidate modifier gene, was expressed in the peripheral nerve, and our functional studies identified and confirmed interacting proteins using coimmunoprecipitation analysis, mass spectrometry, and immunocytochemistry. Chromatin immunoprecipitation and in vitro siRNA experiments were used to analyze gene regulation.

RESULTS

We identified significant association of 4 single nucleotide polymorphisms (rs10910527, rs7536385, rs4649265, rs1547740) in SIPA1L2 with foot dorsiflexion strength (p < 1 × 10^-7 ). Coimmunoprecipitation and mass spectroscopy studies identified β-actin and MYH9 as SIPA1L2 binding partners. Furthermore, we show that SIPA1L2 is part of a myelination-associated coexpressed network regulated by the master transcription factor SOX10. Importantly, in vitro knockdown of SIPA1L2 in Schwannoma cells led to a significant reduction of PMP22 expression, hinting at a potential strategy for drug development.

INTERPRETATION

SIPA1L2 is a potential genetic modifier of CMT1A phenotypic expressions and offers a new pathway to therapeutic interventions. ANN NEUROL 2019;85:316-330.

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.

Regulating PMP22 expression as a dosage sensitive neuropathy gene

Structural variation in the human genome has emerged as a major cause of disease as genomic data have accumulated. One of the most common structural variants associated with human disease causes the heritable neuropathy known as Charcot-Marie-Tooth (CMT) disease type 1A. This 1.4 Mb duplication causes nearly half of the CMT cases that are genetically diagnosed. The PMP22 gene is highly induced in Schwann cells during development, although its precise role in myelin formation and homeostasis is still under active investigation. The PMP22 gene can be considered as a nucleoprotein complex with enzymatic activity to produce the PMP22 transcript, and the complex is allosterically regulated by transcription factors that respond to intracellular signals and epigenomic modifications. The control of PMP22 transcript levels has been one of the major therapeutic targets of therapy development, and this review summarizes those approaches as well as efforts to characterize the regulation of the PMP22 gene.

Report of a patient with a de novo non-recurrent duplication of 17p11.2p12 and Yq11 deletion

Background

The 17p11.2p12 locus is an unstable region that is predisposed to several known genomic disorders and non-recurrent rearrangements that yield varied and wide-ranging phenotypes. Nearly 1% of male newborns have deletions in the Y chromosome; these events primarily involve the heterochromatic region, but may extend to euchromatic Yq segments containing azoospermia factor regions.

Case presentation

We describe the occurrence of two independent chromosomal rearrangements that originated as de novo events in a single male patient: a 10.8-Mb duplication of 17p11.2p12 and a 14.7-Mb deletion of Yq11. This individual shares some clinical characteristics with previously described patients having one or the other of these rearrangements, including global developmental delay, short stature, hypotonia, delayed puberty, certain facial features and a generalized demyelinating sensory-motor polyneuropathy without clinical manifestation. Our patient also presents some features that were not previously described in relevant individuals, including camptodactyly, preauricular pits and hypertrichosis of the back and elbows.

Conclusions

To our knowledge, this is the first patient to be reported with independent de novo deletion/duplication events involving chromosomes 17 and Y. We discuss possible responsible mechanisms and address the phenotype, particularly in light of the clinical features that were not previously reported for patients bearing a duplication of 17p11.2p12 or a deletion of Yq11. We suggest that some of the previously reported patients with Yq11 deletion and clinical manifestations other than male infertility may have additional chromosomal imbalances that could be identified by chromosome microarray analysis, as illustrated by the present case.

Regulation of the neuropathy-associated Pmp22 gene by a distal super-enhancer

Peripheral nerve myelination is adversely affected in the most common form of the hereditary peripheral neuropathy called Charcot-Marie-Tooth Disease. This form, classified as CMT1A, is caused by a 1.4 Mb duplication on chromosome 17, which includes the abundantly expressed Schwann cell myelin gene, Peripheral Myelin Protein 22 (PMP22). This is one of the most common copy number variants causing neurological disease. Overexpression of Pmp22 in rodent models recapitulates several aspects of neuropathy, and reduction of Pmp22 in such models results in amelioration of the neuropathy phenotype. Recently we identified a potential super-enhancer approximately 90-130 kb upstream of the Pmp22 transcription start sites. This super-enhancer encompasses a cluster of individual enhancers that have the acetylated histone H3K27 active enhancer mark, and coincides with smaller duplications identified in patients with milder CMT1A-like symptoms, where the PMP22 coding region itself was not part of the duplication. In this study, we have utilized genome editing to create a deletion of this super-enhancer to determine its role in Pmp22 regulation. Our data show a significant decrease in Pmp22 transcript expression using allele-specific internal controls. Moreover, the P2 promoter of the Pmp22 gene, which is used in other cell types, is affected, but we find that the Schwann cell-specific P1 promoter is disproportionately more sensitive to loss of the super-enhancer. These data show for the first time the requirement of these upstream enhancers for full Pmp22 expression.

A Review of Copy Number Variants in Inherited Neuropathies

The rapid development in the last 10-15 years of microarray technologies, such as oligonucleotide array Comparative Genomic Hybridization (CGH) and Single Nucleotide Polymorphisms (SNP) genotyping array, has improved the identification of fine chromosomal structural variants, ranging in length from kilobases (kb) to megabases (Mb), as an important cause of genetic differences among healthy individuals and also as disease-susceptibility and/or disease-causing factors. Structural genomic variations due to unbalanced chromosomal rearrangements are known as Copy-Number Variants (CNVs) and these include variably sized deletions, duplications, triplications and translocations. CNVs can significantly contribute to human diseases and rearrangements in several dosage-sensitive genes have been identified as an important causative mechanism in the molecular aetiology of Charcot-Marie-Tooth (CMT) disease and of several CMT-related disorders, a group of inherited neuropathies with a broad range of clinical phenotypes, inheritance patterns and causative genes. Duplications or deletions of the dosage-sensitive gene PMP22 mapped to chromosome 17p12 represent the most frequent causes of CMT type 1A and Hereditary Neuropathy with liability to Pressure Palsies (HNPP), respectively. Additionally, CNVs have been identified in patients with other CMT types (e.g., CMT1X, CMT1B, CMT4D) and different hereditary poly- (e.g., giant axonal neuropathy) and focal- (e.g., hereditary neuralgic amyotrophy) neuropathies, supporting the notion of hereditary peripheral nerve diseases as possible genomic disorders and making crucial the identification of fine chromosomal rearrangements in the molecular assessment of such patients. Notably, the application of advanced computational tools in the analysis of Next-Generation Sequencing (NGS) data has emerged in recent years as a powerful technique for identifying a genome-wide scale complex structural variants (e.g., as the ones resulted from balanced rearrangements) and also smaller pathogenic (intragenic) CNVs that often remain beyond the detection limit of most conventional genomic microarray analyses; in the context of inherited neuropathies where more than 70 disease-causing genes have been identified to date, NGS and particularly Whole-Genome Sequencing (WGS) hold the potential to reduce the number of genomic assays required per patient to reach a diagnosis, analyzing with a single test all the Single Nucleotide Variants (SNVs) and CNVs in the genes possibly implicated in this heterogeneous group of disorders.

Co-occurrence of schwannomatosis and rhabdoid tumor predisposition syndrome 1

BACKGROUND

The clinical phenotype associated with germline SMARCB1 mutations has as yet not been fully documented. It is known that germline SMARCB1 mutations may cause rhabdoid tumor predisposition syndrome (RTPS1) or schwannomatosis. However, the co-occurrence of rhabdoid tumor and schwannomas in the same patient has not so far been reported.

METHODS

We investigated a family with members harboring a germline SMARCB1 deletion by means of whole-body MRI as well as high-resolution microstructural magnetic resonance neurography (MRN). Breakpoint-spanning PCRs were performed to characterize the SMARCB1 deletion and its segregation in the family.

RESULTS

The index patient of this family was in complete continuous remission for an atypical teratoid/rhabdoid tumor (AT/RT) treated at the age of 2 years. However, at the age of 21 years, she exhibited paraparesis of her legs and MRI investigations revealed multiple intrathoracic and spinal schwannomas. Breakpoint-spanning PCRs indicated that the germline deletion segregating in the family encompasses 6.4-kb and includes parts of SMARCB1 intron 7, exons 8-9 and 3.3-kb located telomeric to exon 9 including the SMARCB1 3' UTR. The analysis of sequences at the deletion breakpoints showed that the deletion has been caused by replication errors including template-switching. The patient had inherited the deletion from her 56-year-old healthy mother who did not exhibit schwannomas or other tumors as determined by whole-body MRI. However, MRN of the peripheral nerves of the mother's extremities revealed multiple fascicular microlesions which have been previously identified as indicative of schwannomatosis-associated subclinical peripheral nerve pathology.

CONCLUSION

The occurrence of schwannomatosis-associated clinical symptoms independent of the AT/RT as the primary disease should be considered in long-term survivors of AT/RT. Furthermore, our investigations indicate that germline SMARCB1 mutation carriers not presenting RTs or schwannomatosis-associated clinical symptoms may nevertheless exhibit peripheral nerve pathology as revealed by MRN.

Structural variations causing inherited peripheral neuropathies: A paradigm for understanding genomic organization, chromatin interactions, and gene dysregulation

Inherited peripheral neuropathies (IPNs) are a clinically and genetically heterogeneous group of diseases affecting the motor and sensory peripheral nerves. IPNs have benefited from gene discovery and genetic diagnosis using next-generation sequencing with over 80 causative genes available for testing. Despite this success, up to 50% of cases remain genetically unsolved. In the absence of protein coding mutations, noncoding DNA or structural variation (SV) mutations are a possible explanation. The most common IPN, Charcot-Marie-Tooth neuropathy type 1A (CMT1A), is caused by a 1.5 Mb duplication causing trisomy of the dosage sensitive gene PMP22. Using genome sequencing, we recently identified two large genomic rearrangements causing IPN subtypes X-linked CMT (CMTX3) and distal hereditary motor neuropathy (DHMN1), thereby expanding the spectrum of SV mutations causing IPN. Understanding how newly discovered SVs can cause IPN may serve as a useful paradigm to examine the role of topologically associated domains (TADs), chromatin interactions, and gene dysregulation in disease. This review will describe the growing role of SV in the pathogenesis of IPN and the importance of considering this type of mutation in Mendelian diseases where protein coding mutations cannot be identified.

Contiguous 22.1-kb deletion embracing AVPR2 and ARHGAP4 genes at novel breakpoints leads to nephrogenic diabetes insipidus in a Chinese pedigree

Background

It has been reported that mutations in arginine vasopressin type 2 receptor (AVPR2) cause congenital X-linked nephrogenic diabetes insipidus (NDI). However, only a few cases of AVPR2 deletion have been documented in China.

Methods

An NDI pedigree was included in this study, including the proband and his mother. All NDI patients had polyuria, polydipsia, and growth retardation. PCR mapping, long range PCR and sanger sequencing were used to identify genetic causes of NDI.

Results

A novel 22,110 bp deletion comprising AVPR2 and ARH4GAP4 genes was identified by PCR mapping, long range PCR and sanger sequencing. The deletion happened perhaps due to the 4-bp homologous sequence (TTTT) at the junctions of both 5′ and 3′ breakpoints. The gross deletion co-segregates with NDI. After analyzing available data of putative clinical signs of AVPR2 and ARH4GAP4 deletion, we reconsider the potential role of AVPR2 deletion in short stature.

Conclusions

We identified a novel 22.1-kb deletion leading to X-linked NDI in a Chinese pedigree, which would increase the current knowledge in AVPR2 mutation.

Homozygous and hemizygous CNV detection from exome sequencing data in a Mendelian disease cohort

We developed an algorithm, HMZDelFinder, that uses whole exome sequencing (WES) data to identify rare and intragenic homozygous and hemizygous (HMZ) deletions that may represent complete loss-of-function of the indicated gene. HMZDelFinder was applied to 4866 samples in the Baylor–Hopkins Center for Mendelian Genomics (BHCMG) cohort and detected 773 HMZ deletion calls (567 homozygous or 206 hemizygous) with an estimated sensitivity of 86.5% (82% for single-exonic and 88% for multi-exonic calls) and precision of 78% (53% single-exonic and 96% for multi-exonic calls). Out of 773 HMZDelFinder-detected deletion calls, 82 were subjected to array comparative genomic hybridization (aCGH) and/or breakpoint PCR and 64 were confirmed. These include 18 single-exon deletions out of which 8 were exclusively detected by HMZDelFinder and not by any of seven other CNV detection tools examined. Further investigation of the 64 validated deletion calls revealed at least 15 pathogenic HMZ deletions. Of those, 7 accounted for 17–50% of pathogenic CNVs in different disease cohorts where 7.1–11% of the molecular diagnosis solved rate was attributed to CNVs. In summary, we present an algorithm to detect rare, intragenic, single-exon deletion CNVs using WES data; this tool can be useful for disease gene discovery efforts and clinical WES analyses.

Identification of novel candidate disease genes from de novo exonic copy number variants

BACKGROUND

Exon-targeted microarrays can detect small (<1000 bp) intragenic copy number variants (CNVs), including those that affect only a single exon. This genome-wide high-sensitivity approach increases the molecular diagnosis for conditions with known disease-associated genes, enables better genotype-phenotype correlations, and facilitates variant allele detection allowing novel disease gene discovery.

METHODS

We retrospectively analyzed data from 63,127 patients referred for clinical chromosomal microarray analysis (CMA) at Baylor Genetics laboratories, including 46,755 individuals tested using exon-targeted arrays, from 2007 to 2017. Small CNVs harboring a single gene or two to five non-disease-associated genes were identified; the genes involved were evaluated for a potential disease association.

RESULTS

In this clinical population, among rare CNVs involving any single gene reported in 7200 patients (11%), we identified 145 de novo autosomal CNVs (117 losses and 28 intragenic gains), 257 X-linked deletion CNVs in males, and 1049 inherited autosomal CNVs (878 losses and 171 intragenic gains); 111 known disease genes were potentially disrupted by de novo autosomal or X-linked (in males) single-gene CNVs. Ninety-one genes, either recently proposed as candidate disease genes or not yet associated with diseases, were disrupted by 147 single-gene CNVs, including 37 de novo deletions and ten de novo intragenic duplications on autosomes and 100 X-linked CNVs in males. Clinical features in individuals with de novo or X-linked CNVs encompassing at most five genes (224 bp to 1.6 Mb in size) were compared to those in individuals with larger-sized deletions (up to 5 Mb in size) in the internal CMA database or loss-of-function single nucleotide variants (SNVs) detected by clinical or research whole-exome sequencing (WES). This enabled the identification of recently published genes (BPTF, NONO, PSMD12, TANGO2, and TRIP12), novel candidate disease genes (ARGLU1 and STK3), and further confirmation of disease association for two recently proposed disease genes (MEIS2 and PTCHD1). Notably, exon-targeted CMA detected several pathogenic single-exon CNVs missed by clinical WES analyses.

CONCLUSIONS

Together, these data document the efficacy of exon-targeted CMA for detection of genic and exonic CNVs, complementing and extending WES in clinical diagnostics, and the potential for discovery of novel disease genes by genome-wide assay.

Model Organisms Facilitate Rare Disease Diagnosis and Therapeutic Research

Efforts to identify the genetic underpinnings of rare undiagnosed diseases increasingly involve the use of next-generation sequencing and comparative genomic hybridization methods. These efforts are limited by a lack of knowledge regarding gene function, and an inability to predict the impact of genetic variation on the encoded protein function. Diagnostic challenges posed by undiagnosed diseases have solutions in model organism research, which provides a wealth of detailed biological information. Model organism geneticists are by necessity experts in particular genes, gene families, specific organs, and biological functions. Here, we review the current state of research into undiagnosed diseases, highlighting large efforts in North America and internationally, including the Undiagnosed Diseases Network (UDN) (Supplemental Material, File S1) and UDN International (UDNI), the Centers for Mendelian Genomics (CMG), and the Canadian Rare Diseases Models and Mechanisms Network (RDMM). We discuss how merging human genetics with model organism research guides experimental studies to solve these medical mysteries, gain new insights into disease pathogenesis, and uncover new therapeutic strategies.

Biallelic Mutations in CFAP43 and CFAP44 Cause Male Infertility with Multiple Morphological Abnormalities of the Sperm Flagella

Sperm motility is vital to human reproduction. Malformations of sperm flagella can cause male infertility. Men with multiple morphological abnormalities of the flagella (MMAF) have abnormal spermatozoa with absent, short, coiled, bent, and/or irregular-caliber flagella, which impair sperm motility. The known human MMAF-associated genes, such as DNAH1, only account for fewer than 45% of affected individuals. Pathogenic mechanisms in the genetically unexplained MMAF remain to be elucidated. Here, we conducted genetic analyses by using whole-exome sequencing and genome-wide comparative genomic hybridization microarrays in a multi-center cohort of 30 Han Chinese men affected by MMAF. Among them, 12 subjects could not be genetically explained by any known MMAF-associated genes. Intriguingly, we identified compound-heterozygous mutations in CFAP43 in three subjects and a homozygous frameshift mutation in CFAP44 in one subject. All of these recessive mutations were parentally inherited from heterozygous carriers but were absent in 984 individuals from three Han Chinese control populations. CFAP43 and CFAP44, encoding two cilia- and flagella-associated proteins (CFAPs), are specifically or preferentially expressed in the testis. Using CRISPR/Cas9 technology, we generated two knockout models each deficient in mouse ortholog Cfap43 or Cfap44. Notably, both Cfap43- and Cfap44-deficient male mice presented with MMAF phenotypes, whereas the corresponding female mice were fertile. Our experimental observations on human subjects and animal models strongly suggest that biallelic mutations in either CFAP43 or CFAP44 can cause sperm flagellar abnormalities and impair sperm motility. Further investigations on other CFAP-encoding genes in more genetically unexplained MMAF-affected individuals could uncover novel mechanisms underlying sperm flagellar formation.

Efficient CNV breakpoint analysis reveals unexpected structural complexity and correlation of dosage-sensitive genes with clinical severity in genomic disorders

Genomic disorders are the clinical conditions manifested by submicroscopic genomic rearrangements including copy number variants (CNVs). The CNVs can be identified by array-based comparative genomic hybridization (aCGH), the most commonly used technology for molecular diagnostics of genomic disorders. However, clinical aCGH only informs CNVs in the probe-interrogated regions. Neither orientational information nor the resulting genomic rearrangement structure is provided, which is a key to uncovering mutational and pathogenic mechanisms underlying genomic disorders. Long-range polymerase chain reaction (PCR) is a traditional approach to obtain CNV breakpoint junction, but this method is inefficient when challenged by structural complexity such as often found at the PLP1 locus in association with Pelizaeus-Merzbacher disease (PMD). Here we introduced 'capture and single-molecule real-time sequencing' (cap-SMRT-seq) and newly developed 'asymmetry linker-mediated nested PCR walking' (ALN-walking) for CNV breakpoint sequencing in 49 subjects with PMD-associated CNVs. Remarkably, 29 (94%) of the 31 CNV breakpoint junctions unobtainable by conventional long-range PCR were resolved by cap-SMRT-seq and ALN-walking. Notably, unexpected CNV complexities, including inter-chromosomal rearrangements that cannot be resolved by aCGH, were revealed by efficient breakpoint sequencing. These sequence-based structures of PMD-associated CNVs further support the role of DNA replicative mechanisms in CNV mutagenesis, and facilitate genotype-phenotype correlation studies. Intriguingly, the lengths of gained segments by CNVs are strongly correlated with clinical severity in PMD, potentially reflecting the functional contribution of other dosage-sensitive genes besides PLP1. Our study provides new efficient experimental approaches (especially ALN-walking) for CNV breakpoint sequencing and highlights their importance in uncovering CNV mutagenesis and pathogenesis in genomic disorders.

PMP22 exon 4 deletion causes ER retention of PMP22 and a gain-of-function allele in CMT1E

Objective

To determine whether predicted fork stalling and template switching (FoSTeS) during mitosis deletes exon 4 in peripheral myelin protein 22 KD (PMP22) and causes gain‐of‐function mutation associated with peripheral neuropathy in a family with Charcot–Marie–Tooth disease type 1E.

Methods

Two siblings previously reported to have genomic rearrangements predicted to involve exon 4 of PMP22 were evaluated clinically and by electrophysiology. Skin biopsies from the proband were studied by RT‐PCR to determine the effects of the exon 4 rearrangements on exon 4 mRNA expression in myelinating Schwann cells. Transient transfection studies with wild‐type and mutant PMP22 were performed in Cos7 and RT4 cells to determine the fate of the resultant mutant protein.

Results

Both affected siblings had a sensorimotor dysmyelinating neuropathy with severely slow nerve conduction velocities (<10 m/sec). RT‐PCR studies of Schwann cell RNA from one of the siblings demonstrated a complete in‐frame deletion of PMP22 exon 4 (PMP22Δ4). Transfection studies demonstrated that PMP22Δ4 protein is retained within the endoplasmic reticulum and not transported to the plasma membrane.

Conclusions

Our results confirm that that FoSTeS‐mediated genomic rearrangement produced a deletion of exon 4 of PMP22, resulting in expression of both PMP22 mRNA and protein lacking this sequence. In addition, we provide experimental evidence for endoplasmic reticulum retention of the mutant protein suggesting a gain‐of‐function mutational mechanism consistent with the observed CMT1E in this family. PMP22Δ4 is another example of a mutated myelin protein that is misfolded and contributes to the pathogenesis of the neuropathy.

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.

Whole Genome Sequencing Identifies a 78 kb Insertion from Chromosome 8 as the Cause of Charcot-Marie-Tooth Neuropathy CMTX3

With the advent of whole exome sequencing, cases where no pathogenic coding mutations can be found are increasingly being observed in many diseases. In two large, distantly-related families that mapped to the Charcot-Marie-Tooth neuropathy CMTX3 locus at chromosome Xq26.3-q27.3, all coding mutations were excluded. Using whole genome sequencing we found a large DNA interchromosomal insertion within the CMTX3 locus. The 78 kb insertion originates from chromosome 8q24.3, segregates fully with the disease in the two families, and is absent from the general population as well as 627 neurologically normal chromosomes from in-house controls. Large insertions into chromosome Xq27.1 are known to cause a range of diseases and this is the first neuropathy phenotype caused by an interchromosomal insertion at this locus. The CMTX3 insertion represents an understudied pathogenic structural variation mechanism for inherited peripheral neuropathies. Our finding highlights the importance of considering all structural variation types when studying unsolved inherited peripheral neuropathy cases with no pathogenic coding mutations.

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.

Tead1 regulates the expression of Peripheral Myelin Protein 22 during Schwann cell development

Schwann cells are myelinating glia in the peripheral nervous system that form the myelin sheath. A major cause of peripheral neuropathy is a copy number variant involving the Peripheral Myelin Protein 22 (PMP22) gene, which is located within a 1.4-Mb duplication on chromosome 17 associated with the most common form of Charcot-Marie-Tooth Disease (CMT1A). Rodent models of CMT1A have been used to show that reducing Pmp22 overexpression mitigates several aspects of a CMT1A-related phenotype. Mechanistic studies of Pmp22 regulation identified enhancers regulated by the Sox10 (SRY sex determining region Y-box 10) and Egr2/Krox20 (Early growth response protein 2) transcription factors in myelinated nerves. However, relatively little is known regarding how other transcription factors induce Pmp22 expression during Schwann cell development and myelination. Here, we examined Pmp22 enhancers as a function of cell type-specificity, nerve injury and development. While Pmp22 enhancers marked by active histone modifications were lost or remodeled after injury, we found that these enhancers were permissive in early development prior to Pmp22 upregulation. Pmp22 enhancers contain binding motifs for TEA domain (Tead) transcription factors of the Hippo signaling pathway. We discovered that Tead1 and co-activators Yap and Taz are required for Pmp22 expression, as well as for the expression of Egr2 Tead1 directly binds Pmp22 and Egr2 enhancers early in development and Tead1 binding is induced during myelination, correlating with Pmp22 expression. The data identify Tead1 as a novel regulator of Pmp22 expression during development in concert with Sox10 and Egr2.

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.

Application of targeted multi-gene panel testing for the diagnosis of inherited peripheral neuropathy provides a high diagnostic yield with unexpected phenotype-genotype variability

Background

Inherited peripheral neuropathy (IPN) is a clinically and genetically heterogeneous group of disorders with more than 90 genes associated with the different subtypes. Sequential gene screening is gradually being replaced by next generation sequencing (NGS) applications.

Methods

We designed and validated a targeted NGS panel assay including 56 genes associated with known causes of IPN. We report our findings following NGS panel testing of 448 patients with different types of clinically-suspected IPN.

Results

Genetic diagnosis was achieved in 137 patients (31 %) and involved 195 pathogenic variants in 31 genes. 93 patients had pathogenic variants in genes where a resulting phenotype follows dominant inheritance, 32 in genes where this would follow recessive inheritance, and 12 presented with X-linked disease.

Almost half of the diagnosed patients (64) had a pathogenic variant either in genes not previously available for routine diagnostic testing in a UK laboratory (50 patients) or in genes whose primary clinical association was not IPN (14).

Seven patients had a pathogenic variant in a gene not hitherto indicated from their phenotype and three patients had more than one pathogenic variant, explaining their complex phenotype and providing information essential for accurate prediction of recurrence risks.

Conclusions

Our results demonstrate that targeted gene panel testing is an unbiased approach which overcomes the limitations imposed by limited existing knowledge for rare genes, reveals high heterogeneity, and provides high diagnostic yield. It is therefore a highly efficient and cost effective tool for achieving a genetic diagnosis for IPN.

Electronic supplementary material

The online version of this article (doi:10.1186/s12881-015-0224-8) contains supplementary material, which is available to authorized users.

The role of combined SNV and CNV burden in patients with distal symmetric polyneuropathy

PURPOSE

Charcot-Marie-Tooth (CMT) disease is a heterogeneous group of genetic disorders of the peripheral nervous system. Copy-number variants (CNVs) contribute significantly to CMT, as duplication of PMP22 underlies the majority of CMT1 cases. We hypothesized that CNVs and/or single-nucleotide variants (SNVs) might exist in patients with CMT with an unknown molecular genetic etiology.

METHODS

Two hundred patients with CMT, negative for both SNV mutations in several CMT genes and for CNVs involving PMP22, were screened for CNVs by high-resolution oligonucleotide array comparative genomic hybridization. Whole-exome sequencing was conducted on individuals with rare, potentially pathogenic CNVs.

RESULTS

Putatively causative CNVs were identified in five subjects (~2.5%); four of the five map to known neuropathy genes. Breakpoint sequencing revealed Alu-Alu-mediated junctions as a predominant contributor. Exome sequencing identified MFN2 SNVs in two of the individuals.

CONCLUSION

Neuropathy-associated CNV outside of the PMP22 locus is rare in CMT. Nevertheless, there is potential clinical utility in testing for CNVs and exome sequencing in CMT cases negative for the CMT1A duplication. These findings suggest that complex phenotypes including neuropathy can potentially be caused by a combination of SNVs and CNVs affecting more than one disease-associated locus and contributing to a mutational burden.Genet Med 18 5, 443-451.

Targeted Next-Generation Sequencing for Clinical Diagnosis of 561 Mendelian Diseases

Background

Targeted next-generation sequencing (NGS) is a cost-effective approach for rapid and accurate detection of genetic mutations in patients with suspected genetic disorders, which can facilitate effective diagnosis.

Methodology/Principal Findings

We designed a capture array to mainly capture all the coding sequence (CDS) of 2,181 genes associated with 561 Mendelian diseases and conducted NGS to detect mutations. The accuracy of NGS was 99.95%, which was obtained by comparing the genotypes of selected loci between our method and SNP Array in four samples from normal human adults. We also tested the stability of the method using a sample from normal human adults. The results showed that an average of 97.79% and 96.72% of single-nucleotide variants (SNVs) in the sample could be detected stably in a batch and different batches respectively. In addition, the method could detect various types of mutations. Some disease-causing mutations were detected in 69 clinical cases, including 62 SNVs, 14 insertions and deletions (Indels), 1 copy number variant (CNV), 1 microdeletion and 2 microduplications of chromosomes, of which 35 mutations were novel. Mutations were confirmed by Sanger sequencing or real-time polymerase chain reaction (PCR).

Conclusions/Significance

Results of the evaluation showed that targeted NGS enabled to detect disease-causing mutations with high accuracy, stability, speed and throughput. Thus, the technology can be used for the clinical diagnosis of 561 Mendelian diseases.

Decoding NF1 Intragenic Copy-Number Variations

Genomic rearrangements can cause both Mendelian and complex disorders. Currently, several major mechanisms causing genomic rearrangements, such as non-allelic homologous recombination (NAHR), non-homologous end joining (NHEJ), fork stalling and template switching (FoSTeS), and microhomology-mediated break-induced replication (MMBIR), have been proposed. However, to what extent these mechanisms contribute to gene-specific pathogenic copy-number variations (CNVs) remains understudied. Furthermore, few studies have resolved these pathogenic alterations at the nucleotide-level. Accordingly, our aim was to explore which mechanisms contribute to a large, unique set of locus-specific non-recurrent genomic rearrangements causing the genetic neurocutaneous disorder neurofibromatosis type 1 (NF1). Through breakpoint-spanning PCR as well as array comparative genomic hybridization, we have identified the breakpoints in 85 unrelated individuals carrying an NF1 intragenic CNV. Furthermore, we characterized the likely rearrangement mechanisms of these 85 CNVs, along with those of two additional previously published NF1 intragenic CNVs. Unlike the most typical recurrent rearrangements mediated by flanking low-copy repeats (LCRs), NF1 intragenic rearrangements vary in size, location, and rearrangement mechanisms. We propose the DNA-replication-based mechanisms comprising both FoSTeS and/or MMBIR and serial replication stalling to be the predominant mechanisms leading to NF1 intragenic CNVs. In addition to the loop within a 197-bp palindrome located in intron 40, four Alu elements located in introns 1, 2, 3, and 50 were also identified as intragenic-rearrangement hotspots within NF1.

Exome Sequence Analysis Suggests that Genetic Burden Contributes to Phenotypic Variability and Complex Neuropathy

Charcot-Marie-Tooth (CMT) disease is a clinically and genetically heterogeneous distal symmetric polyneuropathy. Whole-exome sequencing (WES) of 40 individuals from 37 unrelated families with CMT-like peripheral neuropathy refractory to molecular diagnosis identified apparent causal mutations in ∼ 45% (17/37) of families. Three candidate disease genes are proposed, supported by a combination of genetic and in vivo studies. Aggregate analysis of mutation data revealed a significantly increased number of rare variants across 58 neuropathy-associated genes in subjects versus controls, confirmed in a second ethnically discrete neuropathy cohort, suggesting that mutation burden potentially contributes to phenotypic variability. Neuropathy genes shown to have highly penetrant Mendelizing variants (HPMVs) and implicated by burden in families were shown to interact genetically in a zebrafish assay exacerbating the phenotype established by the suppression of single genes. Our findings suggest that the combinatorial effect of rare variants contributes to disease burden and variable expressivity.

Non-coding genetic variants in human disease

Genetic variants, including single-nucleotide variants (SNVs) and copy number variants (CNVs), in the non-coding regions of the human genome can play an important role in human traits and complex diseases. Most of the genome-wide association study (GWAS) signals map to non-coding regions and potentially point to non-coding variants, whereas their functional interpretation is challenging. In this review, we discuss the human non-coding variants and their contributions to human diseases in the following four parts. (i) Functional annotations of non-coding SNPs mapped by GWAS: we discuss recent progress revealing some of the molecular mechanisms for GWAS signals affecting gene function. (ii) Technical progress in interpretation of non-coding variants: we briefly describe some of the technologies for functional annotations of non-coding variants, including the methods for genome-wide mapping of chromatin interaction, computational tools for functional predictions and the new genome editing technologies useful for dissecting potential functional consequences of non-coding variants. (iii) Non-coding CNVs in human diseases: we review our emerging understanding the role of non-coding CNVs in human disease. (iv) Compound inheritance of large genomic deletions and non-coding variants: compound inheritance at a locus consisting of coding variants plus non-coding ones is described.

A 3.4-kb Copy-Number Deletion near EPAS1 Is Significantly Enriched in High-Altitude Tibetans but Absent from the Denisovan Sequence

Tibetan high-altitude adaptation (HAA) has been studied extensively, and many candidate genes have been reported. Subsequent efforts targeting HAA functional variants, however, have not been that successful (e.g., no functional variant has been suggested for the top candidate HAA gene, EPAS1). With WinXPCNVer, a method developed in this study, we detected in microarray data a Tibetan-enriched deletion (TED) carried by 90% of Tibetans; 50% were homozygous for the deletion, whereas only 3% carried the TED and 0% carried the homozygous deletion in 2,792 worldwide samples (p < 10(-15)). We employed long PCR and Sanger sequencing technologies to determine the exact copy number and breakpoints of the TED in 70 additional Tibetan and 182 diverse samples. The TED had identical boundaries (chr2: 46,694,276-46,697,683; hg19) and was 80 kb downstream of EPAS1. Notably, the TED was in strong linkage disequilibrium (LD; r(2) = 0.8) with EPAS1 variants associated with reduced blood concentrations of hemoglobin. It was also in complete LD with the 5-SNP motif, which was suspected to be introgressed from Denisovans, but the deletion itself was absent from the Denisovan sequence. Correspondingly, we detected that footprints of positive selection for the TED occurred 12,803 (95% confidence interval = 12,075-14,725) years ago. We further whole-genome deep sequenced (>60×) seven Tibetans and verified the TED but failed to identify any other copy-number variations with comparable patterns, giving this TED top priority for further study. We speculate that the specific patterns of the TED resulted from its own functionality in HAA of Tibetans or LD with a functional variant of EPAS1.

Molecular and clinical features of inherited neuropathies due to PMP22 duplication

PMP22 is a transmembrane glycoprotein component of myelin, important for myelin functioning. Mutation of PMP22 gene which encodes for the production of PMP22 glycoprotein is associated with a variety of inherited neuropathies. This literature review sought to review the molecular mechanism and clinical features of inherited neuropathies caused by PMP22 duplication. PMP22 duplication causes CMT1A which accounts for more than half of all CMT cases and about 70% of CMT1 cases. It manifests with muscle weakness, depressed reflexes, impaired distal sensation, hand and foot deformities, slowing of NCV and onion bulbs. With no specific treatment available, it is managed conservatively. Future treatment may be based on the molecular genetics of the disease.

Alu-mediated diverse and complex pathogenic copy-number variants within human chromosome 17 at p13.3

Alu repetitive elements are known to be major contributors to genome instability by generating Alu-mediated copy-number variants (CNVs). Most of the reported Alu-mediated CNVs are simple deletions and duplications, and the mechanism underlying Alu-Alu-mediated rearrangement has been attributed to non-allelic homologous recombination (NAHR). Chromosome 17 at the p13.3 genomic region lacks extensive low-copy repeat architecture; however, it is highly enriched for Alu repetitive elements, with a fraction of 30% of total sequence annotated in the human reference genome, compared with the 10% genome-wide and 18% on chromosome 17. We conducted mechanistic studies of the 17p13.3 CNVs by performing high-density oligonucleotide array comparative genomic hybridization, specifically interrogating the 17p13.3 region with ∼150 bp per probe density; CNV breakpoint junctions were mapped to nucleotide resolution by polymerase chain reaction and Sanger sequencing. Studied rearrangements include 5 interstitial deletions, 14 tandem duplications, 7 terminal deletions and 13 complex genomic rearrangements (CGRs). Within the 17p13.3 region, Alu-Alu-mediated rearrangements were identified in 80% of the interstitial deletions, 46% of the tandem duplications and 50% of the CGRs, indicating that this mechanism was a major contributor for formation of breakpoint junctions. Our studies suggest that Alu repetitive elements facilitate formation of non-recurrent CNVs, CGRs and other structural aberrations of chromosome 17 at p13.3. The common observation of Alu-mediated rearrangement in CGRs and breakpoint junction sequences analysis further demonstrates that this type of mechanism is unlikely attributed to NAHR, but rather may be due to a recombination-coupled DNA replicative repair process.

Use of Affymetrix Arrays in the Diagnosis of Gene Copy-Number Variation

Diagnosing constitutional pathogenic copy number variants (CNVs) requires detecting submicroscopic segmental chromosomal imbalances. The Affymetrix GeneChip mapping array was one of the initial microarray platforms used to measure duplication and deletion of genetic material in DNA samples. Unlike oligonucleotide microarrays from NimbleGen and Agilent, developed around the same time to infer copy number status for the DNA sequence covered by the probe, the Affymetrix GeneChip system used 25-mer oligonucleotide probes designed to interrogate SNPs. Thus, it was possible to use the Affymetrix 'SNP chips' to both identify SNPs and to identify copy number status. Affymetrix now offers the CytoScan microarray platforms, which are optimized for copy-number analyses, and provides accompanying software. They also offer several other microarray platforms suitable for copy-number analyses. Here we discuss the application of the CytoScan high-density (HD) platform for the detection of genomic imbalance. We provide an overview of the sequence of computational analyses involved in identifying pathogenic CNVs and highlight important parameters for consideration in assessing the pathogenicity of a detected CNV.

PMP22-Related neuropathies and other clinical manifestations in Chinese han patients with charcot-marie-tooth disease type 1

INTRODUCTION

Most cases of Charcot-Marie-Tooth (CMT) disease are caused by mutations in the peripheral myelin protein 22 gene (PMP22), including heterozygous duplications (CMT1A), deletions (HNPP), and point mutations (CMT1E).

METHODS

Single-nucleotide polymorphism (SNP) arrays were used to study PMP22 mutations based on the results of multiplex ligation-dependent probe amplification (MLPA) and polymerase chain reaction-restriction fragment length polymorphism methods in 77 Chinese Han families with CMT1. PMP22 sequencing was performed in MLPA-negative probands. Clinical characteristics were collected for all CMT1A/HNPP probands and their family members.

RESULTS

Twenty-one of 77 CMT1 probands (27.3%) carried duplication/deletion (dup/del) copynumber variants. No point mutations were detected. SNP array and MLPA seem to have similar sensitivity. Fifty-seven patients from 19 CMT1A families had the classical CMT phenotype, except for 1 with concomitant CIDP. Two HNPP probands presented with acute ulnar nerve palsy or recurrent sural nerve palsy, respectively.

CONCLUSIONS

The SNP array has wide coverage, high sensitivity, and high resolution and can be used as a screening tool to detect PMP22 dup/del as shown in this Chinese Han population.

Copy number variations in a population-based study of Charcot-Marie-Tooth disease

Copy number variations (CNVs) are important in relation to diversity and evolution but can sometimes cause disease. The most common genetic cause of the inherited peripheral neuropathy Charcot-Marie-Tooth disease is the PMP22 duplication; otherwise, CNVs have been considered rare. We investigated CNVs in a population-based sample of Charcot-Marie-Tooth (CMT) families. The 81 CMT families had previously been screened for the PMP22 duplication and point mutations in 51 peripheral neuropathy genes, and a genetic cause was identified in 37 CMT families (46%). Index patients from the 44 CMT families with an unknown genetic diagnosis were analysed by whole-genome array comparative genomic hybridization to investigate the entire genome for larger CNVs and multiplex ligation-dependent probe amplification to detect smaller intragenomic CNVs in MFN2 and MPZ. One patient had the pathogenic PMP22 duplication not detected by previous methods. Three patients had potentially pathogenic CNVs in the CNTNAP2, LAMA2, or SEMA5A, that is, genes related to neuromuscular or neurodevelopmental disease. Genotype and phenotype correlation indicated likely pathogenicity for the LAMA2 CNV, whereas the CNTNAP2 and SEMA5A CNVs remained potentially pathogenic. Except the PMP22 duplication, disease causing CNVs are rare but may cause CMT in about 1% (95% CI 0–7%) of the Norwegian CMT families.

Genome architecture and its roles in human copy number variation

Besides single-nucleotide variants in the human genome, large-scale genomic variants, such as copy number variations (CNVs), are being increasingly discovered as a genetic source of human diversity and the pathogenic factors of diseases. Recent experimental findings have shed light on the links between different genome architectures and CNV mutagenesis. In this review, we summarize various genomic features and discuss their contributions to CNV formation. Genomic repeats, including both low-copy and high-copy repeats, play important roles in CNV instability, which was initially known as DNA recombination events. Furthermore, it has been found that human genomic repeats can also induce DNA replication errors and consequently result in CNV mutations. Some recent studies showed that DNA replication timing, which reflects the high-order information of genomic organization, is involved in human CNV mutations. Our review highlights that genome architecture, from DNA sequence to high-order genomic organization, is an important molecular factor in CNV mutagenesis and human genomic instability.

The Functional Copy Number Variation-67048 in WWOX Contributes to Increased Risk of COPD in Southern and Eastern Chinese

Recent studies have recognized the genetic variants in the WW domain-containing oxidoreductase (WWOX) gene as genetic determinants of lung function, reflecting that the WWOX gene may be a susceptible factor of chronic obstructive pulmonary disease (COPD), which characters as poor lung function. We have previously showed that the copy number variation-67048 (CNV-67048) of WWOX was associated with lung cancer risk. Here, we hypothesized that the CNV-67048 affects COPD susceptibility. Based on a two-stage case-control study with a total of 1791 COPD patients and 1940 controls of southern and eastern Chinese, we found that the loss genotypes (0-copy and 1-copy) of CNV-67048 harbored a significantly increased risk of COPD, with an odds ratio (OR) as 1.29 (1.11-1.49) when compared with the common 2-copy genotype. The pre-forced expiratory volume in one second (pre-FEV1) to pre-forced vital capacity (pre-FVC) of carriers with loss genotypes (0.729 ± 0.130) was significantly lower than carriers with 2-copy genotype (0.747 ± 0.124; p = 7.93 × 10(-5)). However, no significant difference was observed on pre-FEV1, pre-FVC and the annual decline of pre-FEV1 between the loss genotypes and 2-copy genotype carriers. Our data suggest that the loss genotypes of CNV-67048 in WWOX predispose their carriers to COPD, which might be a genetic biomarker to predict risk of COPD in Chinese.

Severe phenotypes in a Charcot-Marie-Tooth 1A patient with PMP22 triplication

Charcot-Marie-Tooth disease (CMT) is a genetically and clinically heterogeneous hereditary motor and sensory neuropathy signified by a distal symmetric polyneuropathy. The most frequent subtype is type 1A (CMT1A) caused by duplication in chromosome 17p12 that includes PMP22. This study reports a woman with a family history of CMT1A due to PMP22 duplication. However, she presented with a more severe phenotype than her sibling or ancestors and was found to have a PMP22 triplication instead of the duplication. This was caused by de novo mutation on her affected mother's duplication chromosome. Her lower limb magnetic resonance imaging revealed severe diffused atrophy and fatty replacement. However, her affected sister with typical PMP22 duplication showed almost intact lower limb. Triplication patient's median motor nerve conduction velocity was far lower compared with her sister. Her onset age was faster (8 years) than her sister (42 years). CMT1A triplication might be generated by a female-specific chromosomal rearrangement mechanism that is different from the frequent paternal-originated CMT1A duplication. It also suggests that the wide phenotypic variation of CMT1A might be partly caused by unstable genomic rearrangement, including PMP22 triplication.