The Alu-rich genomic architecture of SPAST predisposes to diverse and functionally distinct disease-associated CNV alleles

SPAST Intragenic CNVs Lead to Hereditary Spastic Paraplegia via a Haploinsufficiency Mechanism

The most common form of hereditary spastic paraplegia (HSP), SPG4 is caused by single nucleotide variants and microrearrangements in the SPAST gene. The high percentage of multi-exonic deletions or duplications observed in SPG4 patients is predisposed by the presence of a high frequency of Alu sequences in the gene sequence. In the present study, we analyzed DNA and RNA samples collected from patients with different microrearrangements in SPAST to map gene breakpoints and evaluate the mutation mechanism. The study group consisted of 69 individuals, including 50 SPG4 patients and 19 healthy relatives from 18 families. Affected family members from 17 families carried varying ranges of microrearrangements in the SPAST gene, while one individual had a single nucleotide variant in the 5'UTR of SPAST. To detect the breakpoints of the SPAST gene, long-range PCR followed by sequencing was performed. The breakpoint sequence was detected for five different intragenic SPAST deletions and one duplication, revealing Alu-mediated microhomology at breakpoint junctions resulting from non-allelic homologous recombination in these patients. Furthermore, SPAST gene expression analysis was performed using patient RNA samples extracted from whole blood. Quantitative real-time PCR tests performed in 14 patients suggest no expression of transcripts with microrearrangements in 5 of them. The obtained data indicate that nonsense-mediated decay degradation is not the only mechanism of hereditary spastic paraplegia in patients with SPAST microrearrangements.

A germline chimeric KANK1-DMRT1 transcript derived from a complex structural variant is associated with a congenital heart defect segregating across five generations

Structural variants (SVs) pose a challenge to detect and interpret, but their study provides novel biological insights and molecular diagnosis underlying rare diseases. The aim of this study was to resolve a 9p24 rearrangement segregating in a family through five generations with a congenital heart defect (congenital pulmonary and aortic valvular stenosis and pulmonary artery stenosis), by applying a combined genomic analysis. The analysis involved multiple techniques, including karyotype, chromosomal microarray analysis (CMA), FISH, genome sequencing (GS), RNA-seq, and optical genome mapping (OGM). A complex 9p24 SV was hinted at by CMA results, showing three interspersed duplicated segments. Combined GS and OGM analyses revealed that the 9p24 duplications constitute a complex SV, on which a set of breakpoints matches the boundaries of the CMA duplicated sequences. The proposed structure for this complex rearrangement implies three duplications associated with an inversion of ~ 2 Mb region on chromosome 9 and a SINE element insertion at the more distal breakpoint. Interestingly, this genomic structure of rearrangement forms a chimeric transcript of the KANK1/DMRT1 loci, which was confirmed by both RNA-seq and Sanger sequencing on blood samples from 9p24 rearrangement carriers. Altogether with breakpoint amplification and FISH analysis, this combined approach allowed a deep characterization of this complex rearrangement. Although the genotype-phenotype correlation remains elusive from the molecular mechanism point of view, this study identified a large genomic rearrangement at 9p24 segregating with a familial congenital heart defect, revealing a genetic biomarker that was successfully applied for embryo selection, changing the reproductive perspective of affected individuals.

LDLR gene rearrangements in Czech FH patients likely arise from one mutational event

BACKGROUND

Large deletions and duplications within the low-density lipoprotein receptor (LDLR) gene make up approximately 10% of LDLR pathogenic variants found in Czech patients with familial hypercholesterolemia. The goal of this study was to test the hypothesis that all probands with each rearrangement share identical breakpoints inherited from a common ancestor and to determine the role of Alu repetitive elements in the generation of these rearrangements.

METHODS

The breakpoint sequence was determined by PCR amplification and Sanger sequencing. To confirm the breakpoint position, an NGS analysis was performed. Haplotype analysis of common LDLR variants was performed using PCR and Sanger sequencing.

RESULTS

The breakpoints of 8 rearrangements within the LDLR gene were analysed, including the four most common LDLR rearrangements in the Czech population (number of probands ranging from 8 to 28), and four less common rearrangements (1-4 probands). Probands with a specific rearrangement shared identical breakpoint positions and haplotypes associated with the rearrangement, suggesting a shared origin from a common ancestor. All breakpoints except for one were located inside an Alu element. In 6 out of 8 breakpoints, there was high homology (≥ 70%) between the two Alu repeats in which the break occurred.

CONCLUSIONS

The most common rearrangements of the LDLR gene in the Czech population likely arose from one mutational event. Alu elements likely played a role in the generation of the majority of rearrangements inside the LDLR gene.

A heterozygous germline deletion within USP8 causes severe neurodevelopmental delay with multiorgan abnormalities

Ubiquitin-specific protease 8 (USP8) is a deubiquitinating enzyme involved in deubiquitinating the enhanced epidermal growth factor receptor for escape from degradation. Somatic variants at a hotspot in USP8 are a cause of Cushing's disease, and a de novo germline USP8 variant at this hotspot has been described only once previously, in a girl with Cushing's disease and developmental delay. In this study, we investigated an exome-negative patient with severe developmental delay, dysmorphic features, and multiorgan dysfunction by long-read sequencing, and identified a 22-kb de novo germline deletion within USP8 (chr15:50469966-50491995 [GRCh38]). The deletion involved the variant hotspot, one rhodanese domain, and two SH3 binding motifs, and was presumed to be generated through nonallelic homologous recombination through Alu elements. Thus, the patient may have perturbation of the endosomal sorting system and mitochondrial autophagy through the USP8 defect. This is the second reported case of a germline variant in USP8.

A germline chimeric KANK1-DMRT1 transcript derived from a complex structural variant is associated with a congenital heart defect segregating across five generations

Structural variants (SVs) pose a challenge to detect and interpret, but their study provides novel biological insights and molecular diagnosis underlying rare diseases. The aim of this study was to resolve a 9p24 rearrangement segregating in a family through five generations with a congenital heart defect (congenital pulmonary and aortic valvular stenosis, and pulmonary artery stenosis), by applying a combined genomic analysis. The analysis involved multiple techniques, including karyotype, chromosomal microarray analysis (CMA), FISH, whole-genome sequencing (WGS), RNA-seq and optical genome mapping (OGM). A complex 9p24 SV was hinted at by CMA results, showing three interspersed duplicated segments. Combined WGS and OGM analyses revealed that the 9p24 duplications constitute a complex SV, on which a set of breakpoints match the boundaries of the CMA duplicated sequences. The proposed structure for this complex rearrangement implies three duplications associated with an inversion of ~ 2Mb region on chromosome 9 with a SINE element insertion at the more distal breakpoint. Interestingly, this hypothesized genomic structure of rearrangement forms a chimeric transcript of the KANK1/DMRT1 loci, which was confirmed by RNA-seq on blood from 9p24 rearrangement carriers. Altogether with breakpoint amplification and FISH analysis, this combined approach allowed a deep characterization of this complex rearrangement. Although the genotype-phenotype correlation remains elusive from the molecular mechanism point of view, this study identified a large genomic rearrangement at 9p segregating with a familial congenital clinical trait, revealing a genetic biomarker that was successfully applied for embryo selection, changing the reproductive perspective of affected individuals.

Meiotic DNA breaks drive multifaceted mutagenesis in the human germ line

Meiotic recombination commences with hundreds of programmed DNA breaks; however, the degree to which they are accurately repaired remains poorly understood. We report that meiotic break repair is eightfold more mutagenic for single-base substitutions than was previously understood, leading to de novo mutation in one in four sperm and one in 12 eggs. Its impact on indels and structural variants is even higher, with 100- to 1300-fold increases in rates per break. We uncovered new mutational signatures and footprints relative to break sites, which implicate unexpected biochemical processes and error-prone DNA repair mechanisms, including translesion synthesis and end joining in meiotic break repair. We provide evidence that these mechanisms drive mutagenesis in human germ lines and lead to disruption of hundreds of genes genome wide.

Copy Number Variations in Hereditary Spastic Paraplegia-Related Genes: Evaluation of an Iranian Hereditary Spastic Paraplegia Cohort and Literature Review

Introduction

In human genetic disorders, copy number variations (CNVs) are considered a considerable underlying cause. CNVs are generally detected by array-based methods but can also be discovered by read-depth analysis of whole-exome sequencing (WES) data. We performed WES-based CNV identification in a cohort of 35 Iranian families with hereditary spastic paraplegia (HSP) patients.

Methods

Thirty-five patients whose routine single-nucleotide variants (SNVs) and insertion/deletion analyses from exome data were unrevealing underwent a pipeline of CNV analysis using the read-depth detection method. Subsequently, a comprehensive search about the existence of CNVs in all 84 known HSP-causing genes was carried out in all reported HSP cases, so far.

Results and Discussion

CNV analysis of exome data indicated that 1 patient harbored a heterozygous deletion in exon 17 of the SPAST gene. Multiplex ligation-dependent probe amplification analysis confirmed this deletion in the proband and his affected father. Literature review demonstrated that, to date, pathogenic CNVs have been identified in 30 out of 84 HSP-causing genes (∼36%). However, CNVs in only 17 of these genes were specifically associated with the HSP phenotype. Among them, CNVs were more common in L1CAM, PLP1, SPAST, SPG7, SPG11, and REEP1 genes. The identification of the CNV in 1 of our patients suggests that WES allows the detection of both SNVs and CNVs from a single method without additional costs and execution time. However, because of intrinsic issues of WES in the detection of large rearrangements, it may not yet be exploited to replace the CNV detection methods in standard clinical practice.

Long-read sequencing revealing intragenic deletions in exome-negative spastic paraplegias

Hereditary spastic paraplegias (HSPs) are a heterogeneous group of neurodegenerative disorders characterized by progressive spasticity and weakness in the lower extremities. To date, a total of 88 types of SPG are known. To diagnose HSP, multiple technologies, including microarray, direct sequencing, multiplex ligation-dependent probe amplification, and short-read next-generation sequencing, are often chosen based on the frequency of HSP subtypes. Exome sequencing (ES) is commonly used. We used ES to analyze ten cases of HSP from eight families. We identified pathogenic variants in three cases (from three different families); however, we were unable to determine the cause of the other seven cases using ES. We therefore applied long-read sequencing to the seven undetermined HSP cases (from five families). We detected intragenic deletions within the SPAST gene in four families, and a deletion within PSEN1 in the remaining family. The size of the deletion ranged from 4.7 to 12.5 kb and involved 1-7 exons. All deletions were entirely included in one long read. We retrospectively performed an ES-based copy number variation analysis focusing on pathogenic deletions, but were not able to accurately detect these deletions. This study demonstrated the efficiency of long-read sequencing in detecting intragenic pathogenic deletions in ES-negative HSP patients.

Alu Retrotransposition Event in SPAST Gene as a Novel Cause of Hereditary Spastic Paraplegia

OBJECTIVES

To diagnose the molecular cause of hereditary spastic paraplegia (HSP) observed in a four-generation family with autosomal dominant inheritance.

METHODS

Multiplex ligation-dependent probe amplification (MLPA), whole-exome sequencing (WES), and RNA sequencing (RNA-seq) of peripheral blood leukocytes were performed. Reverse transcription polymerase chain reaction (RT-PCR) and Sanger sequencing were used to characterize target regions of SPAST.

RESULTS

A 121-bp AluYb9 insertion with a 30-bp poly-A tail flanked by 15-bp direct repeats on both sides was identified in the edge of intron 16 in SPAST that segregated with the disease phenotype.

CONCLUSIONS

We identified an intronic AluYb9 insertion inducing splicing alteration in SPAST causing pure HSP phenotype that was not detected by routine WES analysis. Our findings suggest RNA-seq is a recommended implementation for undiagnosed cases by first-line diagnostic approaches. © 2023 International Parkinson and Movement Disorder Society.

Integrative omics approaches to advance rare disease diagnostics

Over the past decade high-throughput DNA sequencing approaches, namely whole exome and whole genome sequencing became a standard procedure in Mendelian disease diagnostics. Implementation of these technologies greatly facilitated diagnostics and shifted the analysis paradigm from variant identification to prioritisation and evaluation. The diagnostic rates vary widely depending on the cohort size, heterogeneity and disease and range from around 30% to 50% leaving the majority of patients undiagnosed. Advances in omics technologies and computational analysis provide an opportunity to increase these unfavourable rates by providing evidence for disease-causing variant validation and prioritisation. This review aims to provide an overview of the current application of several omics technologies including RNA-sequencing, proteomics, metabolomics and DNA-methylation profiling for diagnostics of rare genetic diseases in general and inborn errors of metabolism in particular.

Exome sequencing in a Romanian Bardet-Biedl syndrome cohort revealed an overabundance of causal BBS12 variants

Bardet-Biedl syndrome (BBS), is an emblematic ciliopathy hallmarked by pleiotropy, phenotype variability, and extensive genetic heterogeneity. BBS is a rare (~1/140,000 to ~1/160,000 in Europe) autosomal recessive pediatric disorder characterized by retinal degeneration, truncal obesity, polydactyly, cognitive impairment, renal dysfunction, and hypogonadism. Twenty-eight genes involved in ciliary structure or function have been implicated in BBS, and explain the molecular basis for ~75%-80% of individuals. To investigate the mutational spectrum of BBS in Romania, we ascertained a cohort of 24 individuals in 23 families. Following informed consent, we performed proband exome sequencing (ES). We detected 17 different putative disease-causing single nucleotide variants or small insertion-deletions and two pathogenic exon disruptive copy number variants in known BBS genes in 17 pedigrees. The most frequently impacted genes were BBS12 (35%), followed by BBS4, BBS7, and BBS10 (9% each) and BBS1, BBS2, and BBS5 (4% each). Homozygous BBS12 p.Arg355* variants were present in seven pedigrees of both Eastern European and Romani origin. Our data show that although the diagnostic rate of BBS in Romania is likely consistent with other worldwide cohorts (74%), we observed a unique distribution of causal BBS genes, including overrepresentation of BBS12 due to a recurrent nonsense variant, that has implications for regional diagnostics.

Copy-number analysis by base-level normalization: An intuitive visualization tool for evaluating copy number variations

Next-generation sequencing (NGS) facilitates comprehensive molecular analyses that help with diagnosing unsolved disorders. In addition to detecting single-nucleotide variations and small insertions/deletions, bioinformatics tools can identify copy number variations (CNVs) in NGS data, which improves the diagnostic yield. However, due to the possibility of false positives, subsequent confirmation tests are generally performed. Here, we introduce Copy-number Analysis by BAse-level NormAlization (CABANA), a visualization tool that allows users to intuitively identify candidate CNVs using the normalized single-base-level read depth calculated from NGS data. To demonstrate how CABANA works, NGS data were obtained from 474 patients with neuromuscular disorders. CNVs were screened using a conventional bioinformatics tool, ExomeDepth, and then we normalized and visualized those data at the single-base level using CABANA, followed by manual inspection by geneticists to filter out false positives and determine candidate CNVs. In doing so, we identified 31 candidate CNVs (7%) in 474 patients and subsequently confirmed all of them to be true using multiplex ligation-dependent probe amplification. The performance of CABANA was deemed acceptable by comparing its diagnostic yield with previous data about neuromuscular disorders. Despite some limitations, we expect CABANA to help researchers accurately identify CNVs and reduce the need for subsequent confirmation testing.

Comprehensive analysis of the PRPF31 gene in retinitis pigmentosa patients: Four novel Alu-mediated copy number variations at the PRPF31 locus

Retinitis pigmentosa (RP) is a monogenic disease characterized by irreversible degeneration of the retina. PRPF31, the second most common causative gene of autosomal dominant RP, frequently harbors copy number variations (CNVs), but the underlying mechanism is unclear. In this study, we summarized the phenotypic and genotypic characteristics of 18 RP families (F01-F18) with variants in PRPF31. The prevalence of PRPF31 variants in our cohort of Chinese RP families was 1.7% (18/1024). Seventeen different variants in PRPF31 were detected, including eight novel variants. Notably, four novel CNVs encompassing PRPF31, with a proportion of 22.2% (4/18), were validated to harbor gross deletions involving Alu/Alu-mediated rearrangements (AAMRs) in the same orientation. Among a total of 12 CNVs of PRPF31 with breakpoints mapped on nucleotide-resolution, 10 variants (83.3%) were presumably mediated by Alu elements. Furthermore, we described the correlation between the genotypes and phenotypes in PRPF31-related RP. Our findings expand the mutational spectrum of the PRPF31 gene and provide strong evidence that Alu elements of PRPF31 probably contribute to the susceptibility to genomic rearrangement in this locus.

Rad52's DNA annealing activity drives template switching associated with restarted DNA replication

It is thought that many of the simple and complex genomic rearrangements associated with congenital diseases and cancers stem from mistakes made during the restart of collapsed replication forks by recombination enzymes. It is hypothesised that this recombination-mediated restart process transitions from a relatively accurate initiation phase to a less accurate elongation phase characterised by extensive template switching between homologous, homeologous and microhomologous DNA sequences. Using an experimental system in fission yeast, where fork collapse is triggered by a site-specific replication barrier, we show that ectopic recombination, associated with the initiation of recombination-dependent replication (RDR), is driven mainly by the Rad51 recombinase, whereas template switching, during the elongation phase of RDR, relies more on DNA annealing by Rad52. This finding provides both evidence and a mechanistic basis for the transition hypothesis.

Comprehensive analysis of F8 large deletions: Characterization of full breakpoint junctions and description of a possible DNA breakage hotspot in intron 6

BACKGROUND

Large F8 deletions represent 3-5% of the variations found in severe hemophilia A patients, but only a few deletion breakpoints have been characterized precisely.

OBJECTIVES

Resolving at the nucleotide level 24 F8 large deletions to provide new data on the mechanisms involved in these rearrangements.

METHODS

Breakpoint junctions of 24 F8 large deletions were characterized using a combination of long-range polymerase chain reaction, whole F8 NGS sequencing, and Sanger sequencing. Repeat elements, non-B DNA, and secondary structures were analyzed around the breakpoints.

RESULTS

Deletions ranged from 1.667 kb to 0.5 Mb in size. Nine involved F8 neighboring genes. Simple blunt ends and 2-4 bp microhomologies were identified at the breakpoint junctions of 10 (42%) and 8 (33%) deletions, respectively. Five (21%) deletions resulted from homeologous recombination between two Alu elements. The remaining case corresponded to a more complex rearrangement with an insertion of a 19 bp-inverted sequence at the junction. Four different breakpoints were located in a 562-bp region in F8 intron 6. This finding suggested that this region, composed of two Alu elements, is a DNA breakage hotspot. Non-B DNA and secondary structures were identified in the junction regions and may contribute to DNA breakage.

CONCLUSION

Molecular characterization of deletion breakpoints revealed that non-homologous non-replicative DNA repair mechanisms and replication-based mechanisms seemed to be the main causative mechanisms of F8 large deletions. Moreover, we identified a possible F8 DNA breakage hotspot involved in non-recurrent rearrangements.

Multi-pathway DNA-repair reporters reveal competition between end-joining, single-strand annealing and homologous recombination at Cas9-induced DNA double-strand breaks

DNA double-strand breaks (DSB) are repaired by multiple distinct pathways, with outcomes ranging from error-free repair to mutagenesis and genomic loss. DSB-repair pathway cross-talk and compensation is incompletely understood, despite its importance for genomic stability, oncogenesis, and genome editing using CRISPR/Cas9. To address this, we constructed and validated three fluorescent Cas9-based reporters, named DSB-Spectrum, that simultaneously quantify the contribution of multiple DNA repair pathways at a DSB. DSB-Spectrum reporters distinguish between DSB-repair by error-free canonical non-homologous end-joining (c-NHEJ) versus homologous recombination (HR; reporter 1), mutagenic repair versus HR (reporter 2), and mutagenic end-joining versus single strand annealing (SSA) versus HR (reporter 3). Using these reporters, we show that inhibiting the c-NHEJ factor DNA-PKcs increases repair by HR, but also substantially increases mutagenic SSA. Our data indicate that SSA-mediated DSB-repair also occurs at endogenous genomic loci, driven by Alu elements or homologous gene regions. Finally, we demonstrate that long-range end-resection factors DNA2 and Exo1 promote SSA and reduce HR, when both pathways compete for the same substrate. These new Cas9-based DSB-Spectrum reporters facilitate the comprehensive analysis of repair pathway crosstalk and DSB-repair outcome.

SVA retrotransposons and a low copy repeat in humans and great apes: a mobile connection

Segmental duplications (SDs) constitute a considerable fraction of primate genomes. They contribute to genetic variation and provide raw material for evolution. Groups of SDs are characterized by the presence of shared core duplicons. One of these core duplicons, low copy repeat (lcr)16a, has been shown to be particularly active in the propagation of interspersed SDs in primates. The underlying mechanisms are, however, only partially understood. Alu short interspersed elements (SINEs) are frequently found at breakpoints and have been implicated in the expansion of SDs.

Detailed analysis of lcr16a-containing SDs shows that the hominid-specific SVA (SINE-R-VNTR-Alu) retrotransposon is an integral component of the core duplicon in Asian and African great apes. In orang-utan, it provides breakpoints and contributes to both interchromosomal and intrachromosomal lcr16a mobility by inter-element recombination. Furthermore, the data suggest that in hominines (human, chimpanzee, gorilla) SVA recombination-mediated integration of a circular intermediate is the founding event of a lineage-specific lcr16a expansion. One of the hominine lcr16a copies displays large flanking direct repeats, a structural feature shared by other SDs in the human genome.

Taken together, the results obtained extend the range of SVAs’ contribution to genome evolution from RNA-mediated transduction to DNA-based recombination. In addition, they provide further support for a role of circular intermediates in SD mobilization.

A guide for the diagnosis of rare and undiagnosed disease: beyond the exome

Rare diseases affect 30 million people in the USA and more than 300–400 million worldwide, often causing chronic illness, disability, and premature death. Traditional diagnostic techniques rely heavily on heuristic approaches, coupling clinical experience from prior rare disease presentations with the medical literature. A large number of rare disease patients remain undiagnosed for years and many even die without an accurate diagnosis. In recent years, gene panels, microarrays, and exome sequencing have helped to identify the molecular cause of such rare and undiagnosed diseases. These technologies have allowed diagnoses for a sizable proportion (25–35%) of undiagnosed patients, often with actionable findings. However, a large proportion of these patients remain undiagnosed. In this review, we focus on technologies that can be adopted if exome sequencing is unrevealing. We discuss the benefits of sequencing the whole genome and the additional benefit that may be offered by long-read technology, pan-genome reference, transcriptomics, metabolomics, proteomics, and methyl profiling. We highlight computational methods to help identify regionally distant patients with similar phenotypes or similar genetic mutations. Finally, we describe approaches to automate and accelerate genomic analysis. The strategies discussed here are intended to serve as a guide for clinicians and researchers in the next steps when encountering patients with non-diagnostic exomes.

Therapeutic Strategies for Mutant SPAST-Based Hereditary Spastic Paraplegia

Mutations of the SPAST gene that encodes the microtubule-severing enzyme called spastin are the chief cause of Hereditary Spastic Paraplegia. Growing evidence indicates that pathogenic mutations functionally compromise the spastin protein and endow it with toxic gain-of-function properties. With each of these two factors potentially relevant to disease etiology, the present article discusses possible therapeutic strategies that may ameliorate symptoms in patients suffering from SPAST-based Hereditary Spastic Paraplegia, which is usually termed SPG4-HSP.

A novel homozygous SLC13A5 whole-gene deletion generated by Alu/Alu-mediated rearrangement in an Iraqi family with epileptic encephalopathy

Biallelic loss-of-function (LoF) of SLC13A5 (solute carrier family 13, member 5) induced deficiency in sodium/citrate transporter (NaCT) causes autosomal recessive developmental epileptic encephalopathy 25 with hypoplastic amelogenesis imperfecta (DEE25; MIM #615905). Many pathogenic SLC13A5 single nucleotide variants (SNVs) and small indels have been described; however, no cases with copy number variants (CNVs) have been sufficiently investigated. We describe a consanguineous Iraqi family harboring an 88.5 kb homozygous deletion including SLC13A5 in Chr17p13.1. The three affected male siblings exhibit neonatal-onset epilepsy with fever-sensitivity, recurrent status epilepticus, global developmental delay/intellectual disability (GDD/ID), and other variable neurological findings as shared phenotypical features of DEE25. Two of the three affected subjects exhibit hypoplastic amelogenesis imperfecta (AI), while the proband shows no evidence of dental abnormalities or AI at 2 years of age with apparently unaffected primary dentition. Characterization of the genomic architecture at this locus revealed evidence for genomic instability generated by an Alu/Alu-mediated rearrangement; confirmed by break-point junction Sanger sequencing. This multiplex family from a distinct population elucidates the phenotypic consequence of complete LoF of SLC13A5 and illustrates the importance of read-depth-based CNV detection in comprehensive exome sequencing analysis to solve cases that otherwise remain molecularly unsolved.

Transposable elements, circular RNAs and mitochondrial transcription in age-related genomic regulation

Our understanding of the molecular regulation of aging and age-related diseases is still in its infancy, requiring in-depth characterization of the molecular landscape shaping these complex phenotypes. Emerging classes of molecules with promise as aging modulators include transposable elements, circRNAs and the mitochondrial transcriptome. Analytical complexity means that these molecules are often overlooked, even though they exhibit strong associations with aging and, in some cases, may directly contribute to its progress. Here, we review the links between these novel factors and age-related phenotypes, and we suggest tools that can be easily incorporated into existing pipelines to better understand the aging process.

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.

Mutational Spectrum of Spast (Spg4) and Atl1 (Spg3a) Genes In Russian Patients With Hereditary Spastic Paraplegia

Hereditary spastic paraplegia (HSP) comprises a heterogeneous group of neurodegenerative disorders, it share common symptom - of progressive lower spastic paraparesis. The most common autosomal dominant (AD) forms of HSP are SPG4 (SPAST gene) and SPG3 (ATL1 gene). In the current research we investigated for the first time the distribution of pathogenic mutations in SPAST and ATL1 genes within a large cohort of Russian HSP patients (122 probands; 69 famillial cases). We determined the frequencies of genetic abnormalities using Sanger sequencing, multiplex ligation-dependent probe amplification (MLPA), and Next Generation Sequencing (NGS) of targeted gene panels. As a result, SPG4 was diagnosed in 30.3% (37/122) of HSP cases, where the familial cases represented 37.7% (26/69) of SPG4. In total 31 pathogenic and likely pathogenic variants were detected in SPAST, with 14 new mutations. Among all detected SPAST variants, 29% were gross deletions and duplications. The proportion of SPG3 variants in Russian cohort was 8.2% (10/122) that were all familial cases. All 10 detected ATL1 mutations were missense substitutions, most of which were in the mutational hot spots of 4, 7, 8, 12 exons, with 2 novel mutations. This work will be helpful for the populational genetics of HSP understanding.

A tailored approach to fusion transcript identification increases diagnosis of rare inherited disease

BACKGROUND

RNA sequencing has been proposed as a means of increasing diagnostic rates in studies of undiagnosed rare inherited disease. Recent studies have reported diagnostic improvements in the range of 7.5-35% by profiling splicing, gene expression quantification and allele specific expression. To-date however, no study has systematically assessed the presence of gene-fusion transcripts in cases of germline disease. Fusion transcripts are routinely identified in cancer studies and are increasingly recognized as having diagnostic, prognostic or therapeutic relevance. Isolated reports exist of fusion transcripts being detected in cases of developmental and neurological phenotypes, and thus, systematic application of fusion detection to germline conditions may further increase diagnostic rates. However, current fusion detection methods are unsuited to the investigation of germline disease due to performance biases arising from their development using tumor, cell-line or in-silico data.

METHODS

We describe a tailored approach to fusion candidate identification and prioritization in a cohort of 47 undiagnosed, suspected inherited disease patients. We modify an existing fusion transcript detection algorithm by eliminating its cell line-derived filtering steps, and instead, prioritize candidates using a custom workflow that integrates genomic and transcriptomic sequence alignment, biological and technical annotations, customized categorization logic, and phenotypic prioritization.

RESULTS

We demonstrate that our approach to fusion transcript identification and prioritization detects genuine fusion events excluded by standard analyses and efficiently removes phenotypically unimportant candidates and false positive events, resulting in a reduced candidate list enriched for events with potential phenotypic relevance. We describe the successful genetic resolution of two previously undiagnosed disease cases through the detection of pathogenic fusion transcripts. Furthermore, we report the experimental validation of five additional cases of fusion transcripts with potential phenotypic relevance.

CONCLUSIONS

The approach we describe can be implemented to enable the detection of phenotypically relevant fusion transcripts in studies of rare inherited disease. Fusion transcript detection has the potential to increase diagnostic rates in rare inherited disease and should be included in RNA-based analytical pipelines aimed at genetic diagnosis.

Mechanistic basis of an epistatic interaction reducing age at onset in hereditary spastic paraplegia

Many genetic neurological disorders exhibit variable expression within affected families, often exemplified by variations in disease age at onset. Epistatic effects (i.e. effects of modifier genes on the disease gene) may underlie this variation, but the mechanistic basis for such epistatic interactions is rarely understood. Here we report a novel epistatic interaction between SPAST and the contiguous gene DPY30, which modifies age at onset in hereditary spastic paraplegia, a genetic axonopathy. We found that patients with hereditary spastic paraplegia caused by genomic deletions of SPAST that extended into DPY30 had a significantly younger age at onset. We show that, like spastin, the protein encoded by SPAST, the DPY30 protein controls endosomal tubule fission, traffic of mannose 6-phosphate receptors from endosomes to the Golgi, and lysosomal ultrastructural morphology. We propose that additive effects on this pathway explain the reduced age at onset of hereditary spastic paraplegia in patients who are haploinsufficient for both genes.

New hypothesis for the etiology of SPAST-based hereditary spastic paraplegia

Mutations of the SPAST gene are the chief cause of hereditary spastic paraplegia. Controversy exists in the medical community as to whether the etiology of the disease is haploinsufficiency or toxic gain-of-function properties of the mutant spastin proteins. In recognition of strong reasons that support each possible mechanism, here we present a novel perspective, based in part on new studies with mouse models and in part on the largest study to date on patients with the disease. We posit that haploinsufficiency does not cause the disease but makes the corticospinal tracts vulnerable to a second hit, which is usually the mutant spastin proteins but could also be proteins generated by mutations of other genes that may or may not cause the disease on their own.

Hereditary spastic paraplegia: gain-of-function mechanisms revealed by new transgenic mouse

Mutations of the SPAST gene, which encodes the microtubule-severing protein spastin, are the most common cause of hereditary spastic paraplegia (HSP). Haploinsufficiency is the prevalent opinion as to the mechanism of the disease, but gain-of-function toxicity of the mutant proteins is another possibility. Here, we report a new transgenic mouse (termed SPASTC448Y mouse) that is not haploinsufficient but expresses human spastin bearing the HSP pathogenic C448Y mutation. Expression of the mutant spastin was documented from fetus to adult, but gait defects reminiscent of HSP (not observed in spastin knockout mice) were adult onset, as is typical of human patients. Results of histological and tracer studies on the mouse are consistent with progressive dying back of corticospinal axons, which is characteristic of the disease. The C448Y-mutated spastin alters microtubule stability in a manner that is opposite to the expectations of haploinsufficiency. Neurons cultured from the mouse display deficits in organelle transport typical of axonal degenerative diseases, and these deficits were worsened by depletion of endogenous mouse spastin. These results on the SPASTC448Y mouse are consistent with a gain-of-function mechanism underlying HSP, with spastin haploinsufficiency exacerbating the toxicity of the mutant spastin proteins. These findings reveal the need for a different therapeutic approach than indicated by haploinsufficiency alone.

Developmental Dysplasia of Hip: Perspectives in Genetic Screening

Development dysplasia of the hip (DDH) is a complex developmental disorder despite being a relatively common condition mainly caused by incompatibility of the femoral head and the abnormal joint socket. Development dysplasia of the hip describes a wide spectrum of disorders ranging from minor acetabular dysplasia to irreducible dislocation of the hip. Modern medicine still suffers from lack of information about screening and precise genetic examination. Genome wide linkage and association studies have brought significant progress to DDH diagnosis. Association studies managed to identify many candidate (susceptible) genes, such as PAPPA2, COL2A1, HOXD9, GDF-5, and TGFB1, which play a considerable role in the pathogenesis of DDH. Early detection of DDH has a big chance to help in preventing further disability and improve the psychological health and quality of life in those children. This emphasizes the importance to establish a universal screening program along with the genetic counseling.

Fine-Scale Characterization of Genomic Structural Variation in the Human Genome Reveals Adaptive and Biomedically Relevant Hotspots

Genomic structural variants (SVs) are distributed nonrandomly across the human genome. The "hotspots" of SVs have been implicated in evolutionary innovations, as well as medical conditions. However, the evolutionary and biomedical features of these hotspots remain incompletely understood. Here, we analyzed data from 2,504 genomes to construct a refined map of 1,148 SV hotspots in human genomes. We confirmed that segmental duplication-related nonallelic homologous recombination is an important mechanistic driver of SV hotspot formation. However, to our surprise, we also found that a majority of SVs in hotspots do not form through such recombination-based mechanisms, suggesting diverse mechanistic and selective forces shaping hotspots. Indeed, our evolutionary analyses showed that the majority of SV hotspots are within gene-poor regions and evolve under relaxed negative selection or neutrality. However, we still found a small subset of SV hotspots harboring genes that are enriched for anthropologically crucial functions and evolve under geography-specific and balancing adaptive forces. These include two independent hotspots on different chromosomes affecting alpha and beta hemoglobin gene clusters. Biomedically, we found that the SV hotspots coincide with breakpoints of clinically relevant, large de novo SVs, significantly more often than genome-wide expectations. For example, we showed that the breakpoints of multiple large SVs, which lead to idiopathic short stature, coincide with SV hotspots. Therefore, the mutational instability in SV hotpots likely enables chromosomal breaks that lead to pathogenic structural variation formations. Overall, our study contributes to a better understanding of the mutational and adaptive landscape of the genome.

Perspectives on the Genomics of HSP Beyond Mendelian Inheritance

Hereditary Spastic Paraplegia is an extraordinarily heterogeneous disease caused by over 50 Mendelian genes. Recent applications of next-generation sequencing, large scale data analysis, and data sharing/matchmaking, have discovered a quickly expanding set of additional HSP genes. Since most recently discovered HSP genes are rare causes of the disease, there is a growing concern of a persisting diagnostic gap, estimated at 30-40%, and even higher for sporadic cases. This missing heritability may not be fully closed by classic Mendelian mutations in protein coding genes. Here we show strategies and published examples of broadening areas of attention for Mendelian and non-Mendelian causes of HSP. We suggest a more inclusive perspective on the potential final architecture of HSP genomics. Efforts to narrow the heritability gap will ultimately lead to more precise and comprehensive genetic diagnoses, which is the starting point for emerging, highly specific gene therapies.

Replicative and non-replicative mechanisms in the formation of clustered CNVs are indicated by whole genome characterization

Clustered copy number variants (CNVs) as detected by chromosomal microarray analysis (CMA) are often reported as germline chromothripsis. However, such cases might need further investigations by massive parallel whole genome sequencing (WGS) in order to accurately define the underlying complex rearrangement, predict the occurrence mechanisms and identify additional complexities. Here, we utilized WGS to delineate the rearrangement structure of 21 clustered CNV carriers first investigated by CMA and identified a total of 83 breakpoint junctions (BPJs). The rearrangements were further sub-classified depending on the patterns observed: I) Cases with only deletions (n = 8) often had additional structural rearrangements, such as insertions and inversions typical to chromothripsis; II) cases with only duplications (n = 7) or III) combinations of deletions and duplications (n = 6) demonstrated mostly interspersed duplications and BPJs enriched with microhomology. In two cases the rearrangement mutational signatures indicated both a breakage-fusion-bridge cycle process and haltered formation of a ring chromosome. Finally, we observed two cases with Alu- and LINE-mediated rearrangements as well as two unrelated individuals with seemingly identical clustered CNVs on 2p25.3, possibly a rare European founder rearrangement.

In conclusion, through detailed characterization of the derivative chromosomes we show that multiple mechanisms are likely involved in the formation of clustered CNVs and add further evidence for chromoanagenesis mechanisms in both “simple” and highly complex chromosomal rearrangements. Finally, WGS characterization adds positional information, important for a correct clinical interpretation and deciphering mechanisms involved in the formation of these rearrangements.

Clustered copy number variants (CNVs) as detected by chromosomal microarray are often reported as germline chromoanagenesis. However, such cases might need further investigation by whole genome sequencing (WGS) to accurately resolve the complexity of the structural rearrangement and predict underlying mutational mechanisms. Here, we used WGS to characterize 83 breakpoint-junctions (BPJs) from 21 clustered CNVs, and outlined the rearrangement connectivity pictures. Cases with only deletions often had additional structural rearrangements, such as insertions and inversions, which could be a result of multiple double-strand DNA breaks followed by non-homologous repair, typical to chromothripsis. In contrast, cases with only duplications or combinations of deletions and duplications, demonstrated mostly interspersed duplications and BPJs enriched with microhomology, consistent with serial template switching during DNA replication (chromoanasynthesis). Only two rearrangements were repeat mediated. In aggregate, our results suggest that multiple CNVs clustered on a single chromosome may arise through either chromothripsis or chromoanasynthesis.

Rare gene deletions in genetic generalized and Rolandic epilepsies

Genetic Generalized Epilepsy (GGE) and benign epilepsy with centro-temporal spikes or Rolandic Epilepsy (RE) are common forms of genetic epilepsies. Rare copy number variants have been recognized as important risk factors in brain disorders. We performed a systematic survey of rare deletions affecting protein-coding genes derived from exome data of patients with common forms of genetic epilepsies. We analysed exomes from 390 European patients (196 GGE and 194 RE) and 572 population controls to identify low-frequency genic deletions. We found that 75 (32 GGE and 43 RE) patients out of 390, i.e. ~19%, carried rare genic deletions. In particular, large deletions (>400 kb) represent a higher burden in both GGE and RE syndromes as compared to controls. The detected low-frequency deletions (1) share genes with brain-expressed exons that are under negative selection, (2) overlap with known autism and epilepsy-associated candidate genes, (3) are enriched for CNV intolerant genes recorded by the Exome Aggregation Consortium (ExAC) and (4) coincide with likely disruptive de novo mutations from the NPdenovo database. Employing several knowledge databases, we discuss the most prominent epilepsy candidate genes and their protein-protein networks for GGE and RE.

LINE- and Alu-containing genomic instability hotspot at 16q24.1 associated with recurrent and nonrecurrent CNV deletions causative for ACDMPV

Transposable elements modify human genome by inserting into new loci or by mediating homology-, microhomology-, or homeology-driven DNA recombination or repair, resulting in genomic structural variation. Alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV) is a rare lethal neonatal developmental lung disorder caused by point mutations or copy-number variant (CNV) deletions of FOXF1 or its distant tissue-specific enhancer. Eighty-five percent of 45 ACDMPV-causative CNV deletions, of which junctions have been sequenced, had at least one of their two breakpoints located in a retrotransposon, with more than half of them being Alu elements. We describe a novel ∼35 kb-large genomic instability hotspot at 16q24.1, involving two evolutionarily young LINE-1 (L1) elements, L1PA2 and L1PA3, flanking AluY, two AluSx, AluSx1, and AluJr elements. The occurrence of L1s at this location coincided with the branching out of the Homo-Pan-Gorilla clade, and was preceded by the insertion of AluSx, AluSx1, and AluJr. Our data show that, in addition to mediating recurrent CNVs, L1 and Alu retrotransposons can predispose the human genome to formation of variably sized CNVs, both of clinical and evolutionary relevance. Nonetheless, epigenetic or other genomic features of this locus might also contribute to its increased instability.

Alu-Alu mediated intragenic duplications in IFT81 and MATN3 are associated with skeletal dysplasias

Skeletal dysplasias are a diverse group of rare Mendelian disorders with clinical and genetic heterogeneity. Here, we used targeted copy number variant (CNV) screening and identified intragenic exonic duplications, formed through Alu-Alu fusion events, in two individuals with skeletal dysplasia and negative exome sequencing results. First, we detected a homozygous tandem duplication of exon 9 and 10 in IFT81 in a boy with Jeune syndrome, or short-rib thoracic dysplasia (SRTD) (MIM# 208500). Western blot analysis did not detect any wild-type IFT81 protein in fibroblasts from the patient with the IFT81 duplication, but only a shorter isoform of IFT81 that was also present in the normal control samples. Complementary zebrafish studies suggested that loss of full-length IFT81 protein but expression of a shorter form of IFT81 protein affects the phenotype while being compatible with life. Second, a de novo tandem duplication of exons 2 to 5 in MATN3 was identified in a girl with multiple epiphyseal dysplasia (MED) type 5 (MIM# 607078). Our data highlights the importance of detection and careful characterization of intragenic duplication CNVs, presenting them as a novel and very rare genetic mechanism in IFT81-related Jeune syndrome and MATN3-related MED.

Predicting human genes susceptible to genomic instability associated with Alu/Alu-mediated rearrangements

Alu elements, the short interspersed element numbering more than 1 million copies per human genome, can mediate the formation of copy number variants (CNVs) between substrate pairs. These Alu/Alu-mediated rearrangements (AAMRs) can result in pathogenic variants that cause diseases. To investigate the impact of AAMR on gene variation and human health, we first characterized Alus that are involved in mediating CNVs (CNV-Alus) and observed that these Alus tend to be evolutionarily younger. We then computationally generated, with the assistance of a supercomputer, a test data set consisting of 78 million Alu pairs and predicted ∼18% of them are potentially susceptible to AAMR. We further determined the relative risk of AAMR in 12,074 OMIM genes using the count of predicted CNV-Alu pairs and experimentally validated the predictions with 89 samples selected by correlating predicted hotspots with a database of CNVs identified by clinical chromosomal microarrays (CMAs) on the genomes of approximately 54,000 subjects. We fine-mapped 47 duplications, 40 deletions, and two complex rearrangements and examined a total of 52 breakpoint junctions of simple CNVs. Overall, 94% of the candidate breakpoints were at least partially Alu mediated. We successfully predicted all (100%) of Alu pairs that mediated deletions (n = 21) and achieved an 87% positive predictive value overall when including AAMR-generated deletions and duplications. We provided a tool, AluAluCNVpredictor, for assessing AAMR hotspots and their role in human disease. These results demonstrate the utility of our predictive model and provide insights into the genomic features and molecular mechanisms underlying AAMR.

RNA Editing and Retrotransposons in Neurology

Compared to sites in protein-coding sequences many more targets undergoing adenosine to inosine (A-to-I) RNA editing were discovered in non-coding regions of human cerebral transcripts, particularly in genetic transposable elements called retrotransposons. We review here the interaction mechanisms of RNA editing and retrotransposons and their impact on normal function and human neurological diseases. Exemplarily, A-to-I editing of retrotransposons embedded in protein-coding mRNAs can contribute to protein abundance and function via circular RNA formation, alternative splicing, and exonization or silencing of retrotransposons. Interactions leading to disease are not very well understood. We describe human diseases with involvement of the central nervous system including inborn errors of metabolism, neurodevelopmental disorders, neuroinflammatory and neurodegenerative and paroxysmal diseases, in which retrotransposons (Alu and/or L1 elements) appear to be causally involved in genetic rearrangements. Sole binding of single-stranded retrotransposon transcripts by RNA editing enzymes rather than enzymatic deamination may have a homeostatic effect on retrotransposon turnover. We also review evidence in support of the emerging pathophysiological function of A-to-I editing of retrotransposons in inflammation and its implication for different neurological diseases including amyotrophic lateral sclerosis, frontotemporal dementia, Alzheimer's and Parkinson's disease, and epilepsy.

Whole-genome sequencing in patients with ciliopathies uncovers a novel recurrent tandem duplication in IFT140

Ciliopathies represent a wide spectrum of rare diseases with overlapping phenotypes and a high genetic heterogeneity. Among those, IFT140 is implicated in a variety of phenotypes ranging from isolated retinis pigmentosa to more syndromic cases. Using whole-genome sequencing in patients with uncharacterized ciliopathies, we identified a novel recurrent tandem duplication of exon 27-30 (6.7 kb) in IFT140, c.3454-488_4182+2588dup p.(Tyr1152_Thr1394dup), missed by whole-exome sequencing. Pathogenicity of the mutation was assessed on the patients' skin fibroblasts. Several hundreds of patients with a ciliopathy phenotype were screened and biallelic mutations were identified in 11 families representing 12 pathogenic variants of which seven are novel. Among those unrelated families especially with a Mainzer-Saldino syndrome, eight carried the same tandem duplication (two at the homozygous state and six at the heterozygous state). In conclusion, we demonstrated the implication of structural variations in IFT140-related diseases expanding its mutation spectrum. We also provide evidences for a unique genomic event mediated by an Alu-Alu recombination occurring on a shared haplotype. We confirm that whole-genome sequencing can be instrumental in the ability to detect structural variants for genomic disorders.

Warning SINEs: Alu elements, evolution of the human brain, and the spectrum of neurological disease

Alu elements are a highly successful family of primate-specific retrotransposons that have fundamentally shaped primate evolution, including the evolution of our own species. Alus play critical roles in the formation of neurological networks and the epigenetic regulation of biochemical processes throughout the central nervous system (CNS), and thus are hypothesized to have contributed to the origin of human cognition. Despite the benefits that Alus provide, deleterious Alu activity is associated with a number of neurological and neurodegenerative disorders. In particular, neurological networks are potentially vulnerable to the epigenetic dysregulation of Alu elements operating across the suite of nuclear-encoded mitochondrial genes that are critical for both mitochondrial and CNS function. Here, we highlight the beneficial neurological aspects of Alu elements as well as their potential to cause disease by disrupting key cellular processes across the CNS. We identify at least 37 neurological and neurodegenerative disorders wherein deleterious Alu activity has been implicated as a contributing factor for the manifestation of disease, and for many of these disorders, this activity is operating on genes that are essential for proper mitochondrial function. We conclude that the epigenetic dysregulation of Alu elements can ultimately disrupt mitochondrial homeostasis within the CNS. This mechanism is a plausible source for the incipient neuronal stress that is consistently observed across a spectrum of sporadic neurological and neurodegenerative disorders.

WNT Signaling Perturbations Underlie the Genetic Heterogeneity of Robinow Syndrome

Locus heterogeneity characterizes a variety of skeletal dysplasias often due to interacting or overlapping signaling pathways. Robinow syndrome is a skeletal disorder historically refractory to molecular diagnosis, potentially stemming from substantial genetic heterogeneity. All current known pathogenic variants reside in genes within the noncanonical Wnt signaling pathway including ROR2, WNT5A, and more recently, DVL1 and DVL3. However, ∼70% of autosomal-dominant Robinow syndrome cases remain molecularly unsolved. To investigate this missing heritability, we recruited 21 families with at least one family member clinically diagnosed with Robinow or Robinow-like phenotypes and performed genetic and genomic studies. In total, four families with variants in FZD2 were identified as well as three individuals from two families with biallelic variants in NXN that co-segregate with the phenotype. Importantly, both FZD2 and NXN are relevant protein partners in the WNT5A interactome, supporting their role in skeletal development. In addition to confirming that clustered -1 frameshifting variants in DVL1 and DVL3 are the main contributors to dominant Robinow syndrome, we also found likely pathogenic variants in candidate genes GPC4 and RAC3, both linked to the Wnt signaling pathway. These data support an initial hypothesis that Robinow syndrome results from perturbation of the Wnt/PCP pathway, suggest specific relevant domains of the proteins involved, and reveal key contributors in this signaling cascade during human embryonic development. Contrary to the view that non-allelic genetic heterogeneity hampers gene discovery, this study demonstrates the utility of rare disease genomic studies to parse gene function in human developmental pathways.

Mutant spastin proteins promote deficits in axonal transport through an isoform-specific mechanism involving casein kinase 2 activation

Mutations of various genes cause hereditary spastic paraplegia (HSP), a neurological disease involving dying-back degeneration of upper motor neurons. From these, mutations in the SPAST gene encoding the microtubule-severing protein spastin account for most HSP cases. Cumulative genetic and experimental evidence suggests that alterations in various intracellular trafficking events, including fast axonal transport (FAT), may contribute to HSP pathogenesis. However, the mechanisms linking SPAST mutations to such deficits remain largely unknown. Experiments presented here using isolated squid axoplasm reveal inhibition of FAT as a common toxic effect elicited by spastin proteins with different HSP mutations, independent of microtubule-binding or severing activity. Mutant spastin proteins produce this toxic effect only when presented as the tissue-specific M1 isoform, not when presented as the ubiquitously-expressed shorter M87 isoform. Biochemical and pharmacological experiments further indicate that the toxic effects of mutant M1 spastins on FAT involve casein kinase 2 (CK2) activation. In mammalian cells, expression of mutant M1 spastins, but not their mutant M87 counterparts, promotes abnormalities in the distribution of intracellular organelles that are correctable by pharmacological CK2 inhibition. Collectively, these results demonstrate isoform-specific toxic effects of mutant M1 spastin on FAT, and identify CK2 as a critical mediator of these effects.

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.

Mapping the genomic landscape of inherited retinal disease genes prioritizes genes prone to coding and noncoding copy-number variations

PurposePart of the hidden genetic variation in heterogeneous genetic conditions such as inherited retinal diseases (IRDs) can be explained by copy-number variations (CNVs). Here, we explored the genomic landscape of IRD genes listed in RetNet to identify and prioritize those genes susceptible to CNV formation.MethodsRetNet genes underwent an assessment of genomic features and of CNV occurrence in the Database of Genomic Variants and literature. CNVs identified in an IRD cohort were characterized using targeted locus amplification (TLA) on extracted genomic DNA.ResultsExhaustive literature mining revealed 1,345 reported CNVs in 81 different IRD genes. Correlation analysis between rankings of genomic features and CNV occurrence demonstrated the strongest correlation between gene size and CNV occurrence of IRD genes. Moreover, we identified and delineated 30 new CNVs in IRD cases, 13 of which are novel and three of which affect noncoding, putative cis-regulatory regions. Finally, the breakpoints of six complex CNVs were determined using TLA in a hypothesis-neutral manner.ConclusionWe propose a ranking of CNV-prone IRD genes and demonstrate the efficacy of TLA for the characterization of CNVs on extracted DNA. Finally, this IRD-oriented CNV study can serve as a paradigm for other genetically heterogeneous Mendelian diseases with hidden genetic variation.

Inter-Fork Strand Annealing causes genomic deletions during the termination of DNA replication

Problems that arise during DNA replication can drive genomic alterations that are instrumental in the development of cancers and many human genetic disorders. Replication fork barriers are a commonly encountered problem, which can cause fork collapse and act as hotspots for replication termination. Collapsed forks can be rescued by homologous recombination, which restarts replication. However, replication restart is relatively slow and, therefore, replication termination may frequently occur by an active fork converging on a collapsed fork. We find that this type of non-canonical fork convergence in fission yeast is prone to trigger deletions between repetitive DNA sequences via a mechanism we call Inter-Fork Strand Annealing (IFSA) that depends on the recombination proteins Rad52, Exo1 and Mus81, and is countered by the FANCM-related DNA helicase Fml1. Based on our findings, we propose that IFSA is a potential threat to genomic stability in eukaryotes.

DOI:http://dx.doi.org/10.7554/eLife.25490.001

An Organismal CNV Mutator Phenotype Restricted to Early Human Development

De novo copy number variants (dnCNVs) arising at multiple loci in a personal genome have usually been considered to reflect cancer somatic genomic instabilities. We describe a multiple dnCNV (MdnCNV) phenomenon in which individuals with genomic disorders carry five to ten constitutional dnCNVs. These CNVs originate from independent formation incidences, are predominantly tandem duplications or complex gains, exhibit breakpoint junction features reminiscent of replicative repair, and show increased de novo point mutations flanking the rearrangement junctions. The active CNV mutation shower appears to be restricted to a transient perizygotic period. We propose that a defect in the CNV formation process is responsible for the "CNV-mutator state," and this state is dampened after early embryogenesis. The constitutional MdnCNV phenomenon resembles chromosomal instability in various cancers. Investigations of this phenomenon may provide unique access to understanding genomic disorders, structural variant mutagenesis, human evolution, and cancer biology.

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.

The Alu neurodegeneration hypothesis: A primate-specific mechanism for neuronal transcription noise, mitochondrial dysfunction, and manifestation of neurodegenerative disease

It is hypothesized that retrotransposons have played a fundamental role in primate evolution and that enhanced neurologic retrotransposon activity in humans may underlie the origin of higher cognitive function. As a potential consequence of this enhanced activity, it is likely that neurons are susceptible to deleterious retrotransposon pathways that can disrupt mitochondrial function. An example is observed in the TOMM40 gene, encoding a β-barrel protein critical for mitochondrial preprotein transport. Primate-specific Alu retrotransposons have repeatedly inserted into TOMM40 introns, and at least one variant associated with late-onset Alzheimer’s disease originated from an Alu insertion event. We provide evidence of enriched Alu content in mitochondrial genes and postulate that Alus can disrupt mitochondrial populations in neurons, thereby setting the stage for progressive neurologic dysfunction. This Alu neurodegeneration hypothesis is compatible with decades of research and offers a plausible mechanism for the disruption of neuronal mitochondrial homeostasis, ultimately cascading into neurodegenerative disease.

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

Most balanced translocations are thought to result mechanistically from nonhomologous end joining or, in rare cases of recurrent events, by nonallelic homologous recombination. Here, we use low-coverage mate pair whole-genome sequencing to fine map rearrangement breakpoint junctions in both phenotypically normal and affected translocation carriers. In total, 46 junctions from 22 carriers of balanced translocations were characterized. Genes were disrupted in 48% of the breakpoints; recessive genes in four normal carriers and known dominant intellectual disability genes in three affected carriers. Finally, seven candidate disease genes were disrupted in five carriers with neurocognitive disabilities (SVOPL, SUSD1, TOX, NCALD, SLC4A10) and one XX-male carrier with Tourette syndrome (LYPD6, GPC5). Breakpoint junction analyses revealed microhomology and small templated insertions in a substantive fraction of the analyzed translocations (17.4%; n = 4); an observation that was substantiated by reanalysis of 37 previously published translocation junctions. Microhomology associated with templated insertions is a characteristic seen in the breakpoint junctions of rearrangements mediated by error-prone replication-based repair mechanisms. Our data implicate that a mechanism involving template switching might contribute to the formation of at least 15% of the interchromosomal translocation events.

Doublet-Mediated DNA Rearrangement-A Novel and Potentially Underestimated Mechanism for the Formation of Recurrent Pathogenic Deletions

Deletions and duplications of genomic DNA contribute to evolution, phenotypic diversity, and human disease. The underlying mechanisms are incompletely understood. We identified deletions of exon 10 of the SPAST gene in two unrelated families with hereditary spastic paraplegia. We excluded a founder event, but observed that the breakpoints map to identical repeat regions. These regions likely represent an intragenic "doublet," that is, an enigmatic class of local duplications. The fusion sequences for both deletions are compatible with recombination-based as well as with replication-based mechanisms. Searching the literature, we identified a partial SLC24A4 deletion that involved two copies of another doublet, and was likely formed in an analogous way. Comparing the SPAST and the SLC24A4 doublets with doublets identified previously suggested that many additional doublets have a high potential for triggering rearrangements. Considering that doublets are still being formed in the human genome, and that they likely create high local instability, we suggest that a two-step mechanism consisting of doublet generation and subsequent doublet-mediated deletion/duplication may underlie certain copy-number changes for which other mechanisms are currently assumed. Further studies are necessary to delineate the significance of the thus-far understudied doublets for the formation of copy-number variation.