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

Complex genetic variation in nearly complete human genomes

Diverse sets of complete human genomes are required to construct a pangenome reference and to understand the extent of complex structural variation. Here, we sequence 65 diverse human genomes and build 130 haplotype-resolved assemblies (130 Mbp median continuity), closing 92% of all previous assembly gaps1,2 and reaching telomere-to-telomere (T2T) status for 39% of the chromosomes. We highlight complete sequence continuity of complex loci, including the major histocompatibility complex (MHC), SMN1/SMN2, NBPF8, and AMY1/AMY2, and fully resolve 1,852 complex structural variants (SVs). In addition, we completely assemble and validate 1,246 human centromeres. We find up to 30-fold variation in α-satellite high-order repeat (HOR) array length and characterize the pattern of mobile element insertions into α-satellite HOR arrays. While most centromeres predict a single site of kinetochore attachment, epigenetic analysis suggests the presence of two hypomethylated regions for 7% of centromeres. Combining our data with the draft pangenome reference1 significantly enhances genotyping accuracy from short-read data, enabling whole-genome inference3 to a median quality value (QV) of 45. Using this approach, 26,115 SVs per sample are detected, substantially increasing the number of SVs now amenable to downstream disease association studies.

Inverted triplications formed by iterative template switches generate structural variant diversity at genomic disorder loci

The duplication-triplication/inverted-duplication (DUP-TRP/INV-DUP) structure is a complex genomic rearrangement (CGR). Although it has been identified as an important pathogenic DNA mutation signature in genomic disorders and cancer genomes, its architecture remains unresolved. Here, we studied the genomic architecture of DUP-TRP/INV-DUP by investigating the DNA of 24 patients identified by array comparative genomic hybridization (aCGH) on whom we found evidence for the existence of 4 out of 4 predicted structural variant (SV) haplotypes. Using a combination of short-read genome sequencing (GS), long-read GS, optical genome mapping, and single-cell DNA template strand sequencing (strand-seq), the haplotype structure was resolved in 18 samples. The point of template switching in 4 samples was shown to be a segment of ∼2.2-5.5 kb of 100% nucleotide similarity within inverted repeat pairs. These data provide experimental evidence that inverted low-copy repeats act as recombinant substrates. This type of CGR can result in multiple conformers generating diverse SV haplotypes in susceptible dosage-sensitive loci.

The Utility of Long-Read Sequencing in Diagnosing Genetic Autosomal Recessive Parkinson's Disease: a genetic screening study

Background

Mutations within the genes PRKN and PINK1 are the leading cause of early onset autosomal recessive Parkinson's disease (PD). However, the genetic cause of most early-onset PD (EOPD) cases still remains unresolved. Long-read sequencing has successfully identified many pathogenic structural variants that cause disease, but this technology has not been widely applied to PD. We recently identified the genetic cause of EOPD in a pair of monozygotic twins by uncovering a complex structural variant that spans over 7 Mb, utilizing Oxford Nanopore Technologies (ONT) long-read sequencing. In this study, we aimed to expand on this and assess whether a second variant could be detected with ONT long-read sequencing in other unresolved EOPD cases reported to carry one heterozygous variant in PRKN or PINK1.

Methods

ONT long-read sequencing was performed on patients with one reported PRKN/PINK1 pathogenic variant. EOPD patients with an age at onset younger than 50 were included in this study. As a positive control, we also included EOPD patients who had already been identified to carry two known PRKN pathogenic variants. Initial genetic testing was performed using either short-read targeted panel sequencing for single nucleotide variants and multiplex ligation-dependent probe amplification (MLPA) for copy number variants.

Results

48 patients were included in this study (PRKN "one-variant" n = 24, PINK1 "one-variant" n = 12, PRKN "two-variants" n = 12). Using ONT long-read sequencing, we detected a second pathogenic variant in six PRKN "one-variant" patients (26%, 6/23) but none in the PINK1 "one-variant" patients (0%, 0/12). Long-read sequencing identified one case with a complex inversion, two instances of structural variant overlap, and three cases of duplication. In addition, in the positive control PRKN "two-variants" group, we were able to identify both pathogenic variants in PRKN in all the patients (100%, 12/12).

Conclusions

This data highlights that ONT long-read sequencing is a powerful tool to identify a pathogenic structural variant at the PRKN locus that is often missed by conventional methods. Therefore, for cases where conventional methods fail to detect a second variant for EOPD, long-read sequencing should be considered as an alternative and complementary approach.

Microsatellite break-induced replication generates highly mutagenized extrachromosomal circular DNAs

Extrachromosomal circular DNAs (eccDNAs) are produced from all regions of the eucaryotic genome. We used inverse PCR of non-B microsatellites capable of forming hairpin, triplex, quadruplex and AT-rich structures integrated at a common ectopic chromosomal site to show that these non-B DNAs generate highly mutagenized eccDNAs by replication-dependent mechanisms. Mutagenesis occurs within the non-B DNAs and extends several kilobases bidirectionally into flanking and nonallelic DNA. Each non-B DNA exhibits a different pattern of mutagenesis, while sister clones containing the same non-B DNA also display distinct patterns of recombination, microhomology-mediated template switching and base substitutions. Mutations include mismatches, short duplications, long nontemplated insertions, large deletions and template switches to sister chromatids and nonallelic chromosomes. Drug-induced replication stress or the depletion of DNA repair factors Rad51, the COPS2 signalosome subunit or POLη change the pattern of template switching and alter the eccDNA mutagenic profiles. We propose an asynchronous capture model based on break-induced replication from microsatellite-induced DNA double strand breaks to account for the generation and circularization of mutagenized eccDNAs and the appearance of genomic homologous recombination deficiency (HRD) scars. These results may help to explain the appearance of tumor eccDNAS and their roles in neoantigen production, oncogenesis and resistance to chemotherapy.

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.

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.

Polyamine Dysregulation and Nucleolar Disruption in Alzheimer's Disease

A hypothesis of Alzheimer's disease etiology is proposed describing how cellular stress induces excessive polyamine synthesis and recycling which can disrupt nucleoli. Polyamines are essential in nucleolar functions, such as RNA folding and ribonucleoprotein assembly. Changes in the nucleolar pool of anionic RNA and cationic polyamines acting as counterions can cause significant nucleolar dynamics. Polyamine synthesis reduces S-adenosylmethionine which, at low levels, triggers tau phosphorylation. Also, polyamine recycling reduces acetyl-CoA needed for acetylcholine, which is low in Alzheimer's disease. Extraordinary nucleolar expansion and/or contraction can disrupt epigenetic control in peri-nucleolar chromatin, such as chromosome 14 with the presenilin-1 gene; chromosome 21 with the amyloid precursor protein gene; chromosome 17 with the tau gene; chromosome 19 with the APOE4 gene; and the inactive X chromosome (Xi; aka "nucleolar satellite") with normally silent spermine synthase (polyamine synthesis) and spermidine/spermine-N1-acetyltransferase (polyamine recycling) alleles. Chromosomes 17, 19 and the Xi have high concentrations of Alu elements which can be transcribed by RNA polymerase III if positioned nucleosomes are displaced from the Alu elements. A sudden flood of Alu RNA transcripts can competitively bind nucleolin which is usually bound to Alu sequences in structural RNAs that stabilize the nucleolar heterochromatic shell. This Alu competition leads to loss of nucleolar integrity with leaking of nucleolar polyamines that cause aggregation of phosphorylated tau. The hypothesis was developed with key word searches (e.g., PubMed) using relevant terms (e.g., Alzheimer's, lupus, nucleolin) based on a systems biology approach and exploring autoimmune disease tautology, gaining synergistic insights from other diseases.

LINE-1 retrotransposition and its deregulation in cancers: implications for therapeutic opportunities

Long interspersed element 1 (LINE-1) is the only protein-coding transposon that is active in humans. LINE-1 propagates in the genome using RNA intermediates via retrotransposition. This activity has resulted in LINE-1 sequences occupying approximately one-fifth of our genome. Although most copies of LINE-1 are immobile, ∼100 copies are retrotransposition-competent. Retrotransposition is normally limited via epigenetic silencing, DNA repair, and other host defense mechanisms. In contrast, LINE-1 overexpression and retrotransposition are hallmarks of cancers. Here, we review mechanisms of LINE-1 regulation and how LINE-1 may promote genetic heterogeneity in tumors. Finally, we discuss therapeutic strategies to exploit LINE-1 biology in cancers.

Break-induced replication underlies formation of inverted triplications and generates unexpected diversity in haplotype structures

Background

The duplication-triplication/inverted-duplication (DUP-TRP/INV-DUP) structure is a type of complex genomic rearrangement (CGR) hypothesized to result from replicative repair of DNA due to replication fork collapse. It is often mediated by a pair of inverted low-copy repeats (LCR) followed by iterative template switches resulting in at least two breakpoint junctions in cis . Although it has been identified as an important mutation signature of pathogenicity for genomic disorders and cancer genomes, its architecture remains unresolved and is predicted to display at least four structural variation (SV) haplotypes.

Results

Here we studied the genomic architecture of DUP-TRP/INV-DUP by investigating the genomic DNA of 24 patients with neurodevelopmental disorders identified by array comparative genomic hybridization (aCGH) on whom we found evidence for the existence of 4 out of 4 predicted SV haplotypes. Using a combination of short-read genome sequencing (GS), long- read GS, optical genome mapping and StrandSeq the haplotype structure was resolved in 18 samples. This approach refined the point of template switching between inverted LCRs in 4 samples revealing a DNA segment of ∼2.2-5.5 kb of 100% nucleotide similarity. A prediction model was developed to infer the LCR used to mediate the non-allelic homology repair.

Conclusions

These data provide experimental evidence supporting the hypothesis that inverted LCRs act as a recombinant substrate in replication-based repair mechanisms. Such inverted repeats are particularly relevant for formation of copy-number associated inversions, including the DUP-TRP/INV-DUP structures. Moreover, this type of CGR can result in multiple conformers which contributes to generate diverse SV haplotypes in susceptible loci .

The molecular mechanisms of recombinant chromosome 18 with parental pericentric inversions and a review of the literature

Chromosomal rearrangements mostly result from non-allelic homologous recombination mediated by low-copy repeats (LCRs) or segmental duplications (SDs). Recent studies on recombinant chromosome 18 (rec (18)) have focused on diagnoses and clinical phenotypes. We diagnosed two cases of prenatal rec (18) and identified precise breakpoint intervals using karyotype and chromosomal microarray analyses. We analyzed the distribution characteristics of breakpoint repetitive elements to infer rearrangement mechanisms and reviewed relevant literature to identify genetic trends. Among the 12 families with 25 pregnancies analyzed, 68% rec (18), 24% spontaneous abortions, and 8% normal births were reported. In the 17 rec (18) cases, 65% presented maternal origin and 35% were paternal. Short-arm breakpoints at p11.31 were reported in 10 cases, whereas the long-arm breakpoints were located at q21.3 (6 cases) and q12 (4 cases). Breakpoints of pericentric inversions on chromosome 18 are concentrated in p11.31, q21.3, and q12 regions. Rearrangements at 18p11.31 are non-recurrent events. ALUs, LINE1s, and MIRs were enriched at the breakpoint regions (1.85 to 3.42-fold enrichment over the entire chromosome 18), while SDs and LCRs were absent. ALU subfamilies had sequence identities of 85.94% and 83.01% between two pair breakpoints. Small repetitive elements may mediate recombination-coupled DNA repair processes, facilitating rearrangements on chromosome 18. Maternal inversion carriers are more prone to abnormal recombination in prenatal families with rec (18). Recombinant chromosomes may present preferential segregation during gamete formation.

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.

The Price of Human Evolution: Cancer-Testis Antigens, the Decline in Male Fertility and the Increase in Cancer

The increasing frequency of general and particularly male cancer coupled with the reduction in male fertility seen worldwide motivated us to seek a potential evolutionary link between these two phenomena, concerning the reproductive transcriptional modules observed in cancer and the expression of cancer-testis antigens (CTA). The phylostratigraphy analysis of the human genome allowed us to link the early evolutionary origin of cancer via the reproductive life cycles of the unicellulars and early multicellulars, potentially driving soma-germ transition, female meiosis, and the parthenogenesis of polyploid giant cancer cells (PGCCs), with the expansion of the CTA multi-families, very late during their evolution. CTA adaptation was aided by retrovirus domestication in the unstable genomes of mammals, for protecting male fertility in stress conditions, particularly that of humans, as compensation for the energy consumption of a large complex brain which also exploited retrotransposition. We found that the early and late evolutionary branches of human cancer are united by the immunity-proto-placental network, which evolved in the Cambrian and shares stress regulators with the finely-tuned sex determination system. We further propose that social stress and endocrine disruption caused by environmental pollution with organic materials, which alter sex determination in male foetuses and further spermatogenesis in adults, bias the development of PGCC-parthenogenetic cancer by default.

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.

Transposable element-mediated rearrangements are prevalent in human genomes

Transposable elements constitute about half of human genomes, and their role in generating human variation through retrotransposition is broadly studied and appreciated. Structural variants mediated by transposons, which we call transposable element-mediated rearrangements (TEMRs), are less well studied, and the mechanisms leading to their formation as well as their broader impact on human diversity are poorly understood. Here, we identify 493 unique TEMRs across the genomes of three individuals. While homology directed repair is the dominant driver of TEMRs, our sequence-resolved TEMR resource allows us to identify complex inversion breakpoints, triplications or other high copy number polymorphisms, and additional complexities. TEMRs are enriched in genic loci and can create potentially important risk alleles such as a deletion in TRIM65, a known cancer biomarker and therapeutic target. These findings expand our understanding of this important class of structural variation, the mechanisms responsible for their formation, and establish them as an important driver of human diversity.

Alu-minating the Mechanisms Underlying Primate Cortex Evolution

The higher-order cognitive functions observed in primates correlate with the evolutionary enhancement of cortical volume and folding, which in turn are driven by the primate-specific expansion of cellular diversity in the developing cortex. Underlying these changes is the diversification of molecular features including the creation of human and/or primate-specific genes, the activation of specific molecular pathways, and the interplay of diverse layers of gene regulation. We review and discuss evidence for connections between Alu elements and primate brain evolution, the evolutionary milestones of which are known to coincide along primate lineages. Alus are repetitive elements that contribute extensively to the acquisition of novel genes and the expansion of diverse gene regulatory layers, including enhancers, alternative splicing, RNA editing, and microRNA pathways. By reviewing the impact of Alus on molecular features linked to cortical expansions or gyrification or implications in cognitive deficits, we suggest that future research focusing on the role of Alu-derived molecular events in the context of brain development may greatly advance our understanding of higher-order cognitive functions and neurologic disorders.

Complex genomic rearrangements: an underestimated cause of rare diseases

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

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.

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.

Is Nucleoredoxin a Master Regulator of Cellular Redox Homeostasis? Its Implication in Different Pathologies

Nucleoredoxin (NXN), an oxidoreductase enzyme, contributes to cellular redox homeostasis by regulating different signaling pathways in a redox-dependent manner. By interacting with seven proteins so far, namely disheveled (DVL), protein phosphatase 2A (PP2A), phosphofructokinase-1 (PFK1), translocation protein SEC63 homolog (SEC63), myeloid differentiation primary response gene-88 (MYD88), flightless-I (FLII), and calcium/calmodulin-dependent protein kinase II type alpha (CAMK2A), NXN is involved in the regulation of several key cellular processes, including proliferation, organogenesis, cell cycle progression, glycolysis, innate immunity and inflammation, motility, contraction, protein transport into the endoplasmic reticulum, neuronal plasticity, among others; as a result, NXN has been implicated in different pathologies, such as cancer, alcoholic and polycystic liver disease, liver fibrogenesis, obesity, Robinow syndrome, diabetes mellitus, Alzheimer’s disease, and retinitis pigmentosa. Together, this evidence places NXN as a strong candidate to be a master redox regulator of cell physiology and as the hub of different redox-sensitive signaling pathways and associated pathologies. This review summarizes and discusses the current insights on NXN-dependent redox regulation and its implication in different pathologies.

Clan genomics: From OMIM phenotypic traits to genes and biology

Clinical characterization of a patient phenotype has been the quintessential approach for elucidating a differential diagnosis and a hypothesis to explore a potential clinical diagnosis. This has resulted in a language of medicine and a semantic ontology, with both specialty- and subspecialty-specific lexicons, that can be challenging to translate and interpret. There is no 'Rosetta Stone' of clinical medicine such as the genetic code that can assist translation and interpretation of the language of genetics. Nevertheless, the information content embodied within a clinical diagnosis can guide management, therapeutic intervention, and potentially prognostic outlook of disease enabling anticipatory guidance for patients and families. Clinical genomics is now established firmly in medical practice. The granularity and informative content of a personal genome is immense. Yet, we are limited in our utility of much of that personal genome information by the lack of functional characterization of the overwhelming majority of computationally annotated genes in the haploid human reference genome sequence. Whereas DNA and the genetic code have provided a 'Rosetta Stone' to translate genetic variant information, clinical medicine, and clinical genomics provide the context to understand human biology and disease. A path forward will integrate deep phenotyping, such as available in a clinical synopsis in the Online Mendelian Inheritance in Man (OMIM) entries, with personal genome analyses.

Long-read technologies identify a hidden inverted duplication in a family with choroideremia

The lack of molecular diagnoses in rare genetic diseases can be explained by limitations of current standard genomic technologies. Upcoming long-read techniques have complementary strengths to overcome these limitations, with a particular strength in identifying structural variants. By using optical genome mapping and long-read sequencing, we aimed to identify the pathogenic variant in a large family with X-linked choroideremia. In this family, aberrant splicing of exon 12 of the choroideremia gene CHM was detected in 2003, but the underlying genomic defect remained elusive. Optical genome mapping and long-read sequencing approaches now revealed an intragenic 1,752 bp inverted duplication including exon 12 and surrounding regions, located downstream of the wild-type copy of exon 12. Both breakpoint junctions were confirmed with Sanger sequencing and segregate with the X-linked inheritance in the family. The breakpoint junctions displayed sequence microhomology suggestive for an erroneous replication mechanism as the origin of the structural variant. The inverted duplication is predicted to result in a hairpin formation of the pre-mRNA with the wild-type exon 12, leading to exon skipping in the mature mRNA. The identified inverted duplication is deemed the hidden pathogenic cause of disease in this family. Our study shows that optical genome mapping and long-read sequencing have significant potential for the identification of (hidden) structural variants in rare genetic diseases.

High prevalence of multilocus pathogenic variation in neurodevelopmental disorders in the Turkish population

Neurodevelopmental disorders (NDDs) are clinically and genetically heterogenous; many such disorders are secondary to perturbation in brain development and/or function. The prevalence of NDDs is > 3%, resulting in significant sociocultural and economic challenges to society. With recent advances in family-based genomics, rare-variant analyses, and further exploration of the Clan Genomics hypothesis, there has been a logarithmic explosion in neurogenetic "disease-associated genes" molecular etiology and biology of NDDs; however, the majority of NDDs remain molecularly undiagnosed. We applied genome-wide screening technologies, including exome sequencing (ES) and whole-genome sequencing (WGS), to identify the molecular etiology of 234 newly enrolled subjects and 20 previously unsolved Turkish NDD families. In 176 of the 234 studied families (75.2%), a plausible and genetically parsimonious molecular etiology was identified. Out of 176 solved families, deleterious variants were identified in 218 distinct genes, further documenting the enormous genetic heterogeneity and diverse perturbations in human biology underlying NDDs. We propose 86 candidate disease-trait-associated genes for an NDD phenotype. Importantly, on the basis of objective and internally established variant prioritization criteria, we identified 51 families (51/176 = 28.9%) with multilocus pathogenic variation (MPV), mostly driven by runs of homozygosity (ROHs) - reflecting genomic segments/haplotypes that are identical-by-descent. Furthermore, with the use of additional bioinformatic tools and expansion of ES to additional family members, we established a molecular diagnosis in 5 out of 20 families (25%) who remained undiagnosed in our previously studied NDD cohort emanating from Turkey.

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

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

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.

Altered DNA repair creates novel Alu/Alu repeat-mediated deletions

Alu elements are the most abundant source of nonallelic homology that influences genetic instability in the human genome. When there is a DNA double-stranded break, the Alu element's high copy number, moderate length and distance and mismatch between elements uniquely influence recombination processes. We utilize a reporter-gene assay to show the complex influence of Alu mismatches on Alu-related repeat-mediated deletions (RMDs). The Alu/Alu heteroduplex intermediate can result in a nonallelic homologous recombination (HR). Alternatively, the heteroduplex can result in various DNA breaks around the Alu elements caused by competing nucleases. These breaks can undergo Alt-nonhomologous end joining to cause deletions focused around the Alu elements. Formation of these heteroduplex intermediates is largely RAD52 dependent. Cells with low ERCC1 levels utilize more of these alternatives resolutions, while cells with MSH2 defects tend to have more RMDs with a specific increase in the HR events. Therefore, Alu elements are expected to create different forms of deletions in various cancers depending on a number of these DNA repair defects.

Complex genomic alterations and intellectual disability: an interpretative challenge

BACKGROUND

Complex chromosomal rearrangements (CCRs) are structural rearrangements involving more than three chromosomes or having more than two breaks; approximately 70% are not associated with any clinical phenotype. Here, we describe a CCR segregating in a two-generation family.

METHOD

A 4-year-old male was evaluated for developmental delay, mild intellectual disability and epicanthus. Karyotype, fluorescence in situ hybridisation (FISH) analysis and array comparative genomic hybridisation (aCGH) analysis were performed on the patient and of all family members.

RESULT

Array CGH analysis of the proband detected two non-contiguous genomic gains of chromosome 2 at bands q32.3q33.2 and bands q36.1q36.3. Both karyotype and FISH analysis revealed a recombinant chromosome 2 with a direct insertion of regions q32.3q33.2 and q36.1q36.3 into region q12. Both of these regions were also present in their original location. Karyotype and FISH analysis of the father revealed a de novo direct insertion of regions q32.3q33.2 and q36.1q36.3 into region q12. Moreover, a de novo balanced translocation involving the q arm of the same chromosome 2 and the p arm of chromosome 10 was observed in the father of the proband. The single nucleotide polymorphism (SNP) array analysis and haplotype reconstruction confirmed the paternal origin of the duplications. Karyotype, FISH analysis and array CGH analysis of other family members were all normal.

CONCLUSION

This report underlines the importance of using different methods to correctly evaluate the origin and the structure of CCRs in order to provide an appropriate management of the patients and a good estimation of the reproductive risk of the family.

Structural Variants Create New Topological-Associated Domains and Ectopic Retinal Enhancer-Gene Contact in Dominant Retinitis Pigmentosa

The cause of autosomal-dominant retinitis pigmentosa (adRP), which leads to loss of vision and blindness, was investigated in families lacking a molecular diagnosis. A refined locus for adRP on Chr17q22 (RP17) was delineated through genotyping and genome sequencing, leading to the identification of structural variants (SVs) that segregate with disease. Eight different complex SVs were characterized in 22 adRP-affected families with >300 affected individuals. All RP17 SVs had breakpoints within a genomic region spanning YPEL2 to LINC01476. To investigate the mechanism of disease, we reprogrammed fibroblasts from affected individuals and controls into induced pluripotent stem cells (iPSCs) and differentiated them into photoreceptor precursor cells (PPCs) or retinal organoids (ROs). Hi-C was performed on ROs, and differential expression of regional genes and a retinal enhancer RNA at this locus was assessed by qPCR. The epigenetic landscape of the region, and Hi-C RO data, showed that YPEL2 sits within its own topologically associating domain (TAD), rich in enhancers with binding sites for retinal transcription factors. The Hi-C map of RP17 ROs revealed creation of a neo-TAD with ectopic contacts between GDPD1 and retinal enhancers, and modeling of all RP17 SVs was consistent with neo-TADs leading to ectopic retinal-specific enhancer-GDPD1 accessibility. qPCR confirmed increased expression of GDPD1 and increased expression of the retinal enhancer that enters the neo-TAD. Altered TAD structure resulting in increased retinal expression of GDPD1 is the likely convergent mechanism of disease, consistent with a dominant gain of function. Our study highlights the importance of SVs as a genomic mechanism in unsolved Mendelian diseases.

Novel pathogenic genomic variants leading to autosomal dominant and recessive Robinow syndrome

Robinow syndrome (RS) is a genetically heterogeneous disorder characterized by skeletal dysplasia and a distinctive facial appearance. Previous studies have revealed locus heterogeneity with rare variants in DVL1, DVL3, FZD2, NXN, ROR2, and WNT5A underlying the etiology of RS. The aforementioned "Robinow-associated genes" and their gene products all play a role in the WNT/planar cell polarity signaling pathway. We performed gene-targeted Sanger sequencing, exome sequencing, genome sequencing, and array comparative genomic hybridization on four subjects with a clinical diagnosis of RS who had not had prior DNA testing. Individuals in our cohort were found to carry pathogenic or likely pathogenic variants in three RS related genes: DVL1, ROR2, and NXN. One subject was found to have a nonsense variant (c.817C > T [p.Gln273*]) in NXN in trans with an ~1 Mb telomeric deletion on chromosome 17p containing NXN, which supports our contention that biallelic NXN variant alleles are responsible for a novel autosomal recessive RS locus. These findings provide increased understanding of the role of WNT signaling in skeletal development and maintenance. These data further support the hypothesis that dysregulation of the noncanonical WNT pathway in humans gives rise to RS.

Cytogenetically visible inversions are formed by multiple molecular mechanisms

Cytogenetically detected inversions are generally assumed to be copy number and phenotypically neutral events. While nonallelic homologous recombination is thought to play a major role, recent data suggest the involvement of other molecular mechanisms in inversion formation. Using a combination of short-read whole-genome sequencing (WGS), 10X Genomics Chromium WGS, droplet digital polymerase chain reaction and array comparative genomic hybridization we investigated the genomic structure of 18 large unique cytogenetically detected chromosomal inversions and achieved nucleotide resolution of at least one chromosomal inversion junction for 13/18 (72%). Surprisingly, we observed that seemingly copy number neutral inversions can be accompanied by a copy-number gain of up to 350 kb and local genomic complexities (3/18, 17%). In the resolved inversions, the mutational signatures are consistent with nonhomologous end-joining (8/13, 62%) or microhomology-mediated break-induced replication (5/13, 38%). Our study indicates that short-read 30x coverage WGS can detect a substantial fraction of chromosomal inversions. Moreover, replication-based mechanisms are responsible for approximately 38% of those events leading to a significant proportion of inversions that are actually accompanied by additional copy-number variation potentially contributing to the overall phenotypic presentation of those patients.

Deciphering the complexity of simple chromosomal insertions by genome sequencing

Chromosomal insertions are thought to be rare structural rearrangements. The current understanding of the underlying mechanisms of their origin is still limited. In this study, we sequenced 16 cases with apparent simple insertions previously identified by karyotyping and/or chromosomal microarray analysis. Using mate-pair genome sequencing (GS), we identified all 16 insertions and revised previously designated karyotypes in 75.0% (12/16) of the cases. Additional cryptic rearrangements were identified in 68.8% of the cases (11/16). The incidence of additional cryptic rearrangements in chromosomal insertions was significantly higher compared to balanced translocations and inversions reported in other studies by GS. We characterized and classified the cryptic insertion rearrangements into four groups, which were not mutually exclusive: (1) insertion segments were fragmented and their subsegments rearranged and clustered at the insertion site (10/16, 62.5%); (2) one or more cryptic subsegments were not inserted into the insertion site (5/16, 31.3%); (3) segments of the acceptor chromosome were scattered and rejoined with the insertion segments (2/16, 12.5%); and (4) copy number gains were identified in the flanking regions of the insertion site (2/16, 12.5%). In addition to the observation of these chromothripsis- or chromoanasynthesis-like events, breakpoint sequence analysis revealed microhomology to be the predominant feature. However, no significant correlation was found between the number of cryptic rearrangements and the size of the insertion. Overall, our study provide molecular characterization of karyotypically apparent simple insertions, demonstrate previously underappreciated complexities, and evidence that chromosomal insertions are likely formed by nonhomologous end joining and/or microhomology-mediated replication-based DNA repair.

Copy number variations and constitutional chromothripsis (Review)

Both copy number variations (CNVs) and chromothripsis are phenomena that involve complex genomic rearrangements. Chromothripsis results in CNVs and other structural changes. CNVs are frequently observed in the human genome. Studies on CNVs have been increasing exponentially; the Database of Genomic Variants shows an increase in the number of data published on structural variations added to the database in the last 15 years. CNVs may be a result of replicative and non-replicative mechanisms, and are hypothesized to serve important roles in human health and disease. Chromothripsis is a phenomena of chromosomal rearrangement following chromosomal breaks at multiple locations and involves impaired DNA repair. In 2011, Stephens et al coined the term chromothripsis for this type of fragmenting event. Several proposed mechanisms have been suggested to underlie chromothripsis, such as p53 inactivation, micronuclei formation, abortive apoptosis and telomere fusions in telomere crisis. Chromothripsis gives rise to normal or abnormal phenotypes. In this review, constitutional chromothripsis, which may coexist with multiple de novo CNVs are described and discussed. This reviews aims to summarize recent advances in our understanding of CNVs and chromothripsis, and describe the effects of these phenomena on human health and birth defects.

Structural variant identification and characterization

Structural variant (SV) differences between human genomes can cause germline and mosaic disease as well as inter-individual variation. De-regulation of accurate DNA repair and genomic surveillance mechanisms results in a large number of SVs in cancer. Analysis of the DNA sequences at SV breakpoints can help identify pathways of mutagenesis and regions of the genome that are more susceptible to rearrangement. Large-scale SV analyses have been enabled by high-throughput genome-level sequencing on humans in the past decade. These studies have shed light on the mechanisms and prevalence of complex genomic rearrangements. Recent advancements in both sequencing and other mapping technologies as well as calling algorithms for detection of genomic rearrangements have helped propel SV detection into population-scale studies, and have begun to elucidate previously inaccessible regions of the genome. Here, we discuss the genomic organization of simple and complex SVs, the molecular mechanisms of their formation, and various ways to detect them. We also introduce methods for characterizing SVs and their consequences on human genomes.

Structural variation and its potential impact on genome instability: Novel discoveries in the EGFR landscape by long-read sequencing

Structural variation (SV) is typically defined as variation within the human genome that exceeds 50 base pairs (bp). SV may be copy number neutral or it may involve duplications, deletions, and complex rearrangements. Recent studies have shown SV to be associated with many human diseases. However, studies of SV have been challenging due to technological constraints. With the advent of third generation (long-read) sequencing technology, exploration of longer stretches of DNA not easily examined previously has been made possible. In the present study, we utilized third generation (long-read) sequencing techniques to examine SV in the EGFR landscape of four haplotypes derived from two human samples. We analyzed the EGFR gene and its landscape (+/- 500,000 base pairs) using this approach and were able to identify a region of non-coding DNA with over 90% similarity to the most common activating EGFR mutation in non-small cell lung cancer. Based on previously published Alu-element genome instability algorithms, we propose a molecular mechanism to explain how this non-coding region of DNA may be interacting with and impacting the stability of the EGFR gene and potentially generating this cancer-driver gene. By these techniques, we were also able to identify previously hidden structural variation in the four haplotypes and in the human reference genome (hg38). We applied previously published algorithms to compare the relative stabilities of these five different EGFR gene landscape haplotypes to estimate their relative potentials to generate the EGFR exon 19, 15 bp canonical deletion. To our knowledge, the present study is the first to use the differences in genomic architecture between targeted cancer-linked phased haplotypes to estimate their relative potentials to form a common cancer-linked driver mutation.

BRCA1 intronic Alu elements drive gene rearrangements and PARP inhibitor resistance

BRCA1 mutant carcinomas are sensitive to PARP inhibitor (PARPi) therapy; however, resistance arises. BRCA1 BRCT domain mutant proteins do not fold correctly and are subject to proteasomal degradation, resulting in PARPi sensitivity. In this study, we show that cell lines and patient-derived tumors, with highly disruptive BRCT domain mutations, have readily detectable BRCA1 protein expression, and are able to proliferate in the presence of PARPi. Peptide analyses reveal that chemo-resistant cancers contain residues encoded by BRCA1 intron 15. Mechanistically, cancers with BRCT domain mutations harbor BRCA1 gene breakpoints within or adjacent to Alu elements in intron 15; producing partial gene duplications, inversions and translocations, and terminating transcription prior to the mutation-containing BRCT domain. BRCA1 BRCT domain-deficient protein isoforms avoid mutation-induced proteasomal degradation, support homology-dependent DNA repair, and promote PARPi resistance. Taken together, Alu-mediated BRCA1 gene rearrangements are responsible for generating hypomorphic proteins, and may represent a biomarker of PARPi resistance.

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.

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

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

From cytogenetics to cytogenomics: whole-genome sequencing as a first-line test comprehensively captures the diverse spectrum of disease-causing genetic variation underlying intellectual disability

Background

Since different types of genetic variants, from single nucleotide variants (SNVs) to large chromosomal rearrangements, underlie intellectual disability, we evaluated the use of whole-genome sequencing (WGS) rather than chromosomal microarray analysis (CMA) as a first-line genetic diagnostic test.

Methods

We analyzed three cohorts with short-read WGS: (i) a retrospective cohort with validated copy number variants (CNVs) (cohort 1, n = 68), (ii) individuals referred for monogenic multi-gene panels (cohort 2, n = 156), and (iii) 100 prospective, consecutive cases referred to our center for CMA (cohort 3). Bioinformatic tools developed include FindSV, SVDB, Rhocall, Rhoviz, and vcf2cytosure.

Results

First, we validated our structural variant (SV)-calling pipeline on cohort 1, consisting of three trisomies and 79 deletions and duplications with a median size of 850 kb (min 500 bp, max 155 Mb). All variants were detected. Second, we utilized the same pipeline in cohort 2 and analyzed with monogenic WGS panels, increasing the diagnostic yield to 8%. Next, cohort 3 was analyzed by both CMA and WGS. The WGS data was processed for large (> 10 kb) SVs genome-wide and for exonic SVs and SNVs in a panel of 887 genes linked to intellectual disability as well as genes matched to patient-specific Human Phenotype Ontology (HPO) phenotypes. This yielded a total of 25 pathogenic variants (SNVs or SVs), of which 12 were detected by CMA as well. We also applied short tandem repeat (STR) expansion detection and discovered one pathologic expansion in ATXN7. Finally, a case of Prader-Willi syndrome with uniparental disomy (UPD) was validated in the WGS data.

Important positional information was obtained in all cohorts. Remarkably, 7% of the analyzed cases harbored complex structural variants, as exemplified by a ring chromosome and two duplications found to be an insertional translocation and part of a cryptic unbalanced translocation, respectively.

Conclusion

The overall diagnostic rate of 27% was more than doubled compared to clinical microarray (12%). Using WGS, we detected a wide range of SVs with high accuracy. Since the WGS data also allowed for analysis of SNVs, UPD, and STRs, it represents a powerful comprehensive genetic test in a clinical diagnostic laboratory setting.

Break-induced replication plays a prominent role in long-range repeat-mediated deletion

Repetitive DNA sequences are often associated with chromosomal rearrangements in cancers. Conventionally, single-strand annealing (SSA) is thought to mediate homology-directed repair of double-strand breaks (DSBs) between two repeats, causing repeat-mediated deletion (RMD). In this report, we demonstrate that break-induced replication (BIR) is used predominantly over SSA in mammalian cells for mediating RMD, especially when repeats are far apart. We show that SSA becomes inefficient in mammalian cells when the distance between the DSBs and the repeats is increased to the 1-2 kb range, while BIR-mediated RMD (BIR/RMD) can act over a long distance (e.g., ~ 100-200 kb) when the DSB is close to one repeat. Importantly, oncogene expression potentiates BIR/RMD but not SSA, and BIR/RMD is used more frequently at single-ended DSBs formed at collapsed replication forks than at double-ended DSBs. In contrast to short-range SSA, H2AX is required for long-range BIR/RMD, and sequence divergence strongly suppresses BIR/RMD in a manner partially dependent on MSH2. Our finding that BIR/RMD has a more important role than SSA in mammalian cells has a significant impact on the understanding of repeat-mediated rearrangements associated with oncogenesis.

An Evolutionary Perspective on the Impact of Genomic Copy Number Variation on Human Health

Copy number variants (CNVs), deletions and duplications of segments of DNA, account for at least five times more variable base pairs in humans than single-nucleotide variants. Several common CNVs were shown to change coding and regulatory sequences and thus dramatically affect adaptive phenotypes involving immunity, perception, metabolism, skin structure, among others. Some of these CNVs were also associated with susceptibility to cancer, infection, and metabolic disorders. These observations raise the possibility that CNVs are a primary contributor to human phenotypic variation and consequently evolve under selective pressures. Indeed, locus-specific haplotype-level analyses revealed signatures of natural selection on several CNVs. However, more traditional tests of selection which are often applied to single-nucleotide variation often have diminished statistical power when applied to CNVs because they often do not show strong linkage disequilibrium with nearby variants. Recombination-based formation mechanisms of CNVs lead to frequent recurrence and gene conversion events, breaking the linkage disequilibrium involving CNVs. Similar methodological challenges also prevent routine genome-wide association studies to adequately investigate the impact of CNVs on heritable human disease. Thus, we argue that the full relevance of CNVs to human health and evolution is yet to be elucidated. We further argue that a holistic investigation of formation mechanisms within an evolutionary framework would provide a powerful framework to understand the functional and biomedical impact of CNVs. In this paper, we review several cases where studies reveal diverse evolutionary histories and unexpected functional consequences of CNVs. We hope that this review will encourage further work on CNVs by both evolutionary and medical geneticists.

Generic Repeat Finder: A High-Sensitivity Tool for Genome-Wide De Novo Repeat Detection

Comprehensive and accurate annotation of the repeatome, including transposons, is critical for deepening our understanding of repeat origins, biogenesis, regulatory mechanisms, and roles. Here, we developed Generic Repeat Finder (GRF), a tool for genome-wide repeat detection based on fast, exhaustive numerical calculation algorithms integrated with optimized dynamic programming strategies. GRF sensitively identifies terminal inverted repeats (TIRs), terminal direct repeats (TDRs), and interspersed repeats that bear both inverted and direct repeats. GRF also detects DNA or RNA transposable elements characterized by these repeats in plant and animal genomes. For TIRs and TDRs, GRF identifies spacers in the middle and mismatches/insertions or deletions in terminal repeats, showing their alignment or base-pairing information. GRF helps improve the annotation for various DNA transposons and retrotransposons, such as miniature inverted-repeat transposable elements (MITEs), long terminal repeat (LTR) retrotransposons, and non-LTR retrotransposons, including long interspersed nuclear elements and short interspersed nuclear elements in plants. We used GRF to perform TIR/TDR, interspersed-repeat, and MITE detection in several species, including Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa), and mouse (Mus musculus). As a generic bioinformatics tool in repeat finding implemented as a parallelized C++ program, GRF was faster and more sensitive than the existing inverted repeat/MITE detection tools based on numerical approaches (i.e. detectIR and detectMITE) in Arabidopsis and mouse. GRF is more sensitive than Inverted Repeat Finder in TIR detection, LTR_FINDER in short TDR detection (≤1,000 nt), and phRAIDER in interspersed repeat detection in Arabidopsis and rice. GRF is an open source available from Github.

Characterization of large deletions of the MECP2 gene in Rett syndrome patients by gene dosage analysis

Background

Rett syndrome (RTT) is a developmental disorder with an early onset and X‐linked dominant inheritance pattern. It is first recognized in infancy and is seen almost always in girls, but it may be seen in boys on rare occasions. Typical RTT is caused by de novo mutations of the gene MECP2 (OMIM*300005), and atypical forms of RTT can be caused by mutations of the CDKL5 (OMIM*300203) and FOXG1 (OMIM*164874) genes.

Methods

Approximately 5% of the mutations detected in MECP2 are large rearrangements that range from exons to the entire gene. Here, we have characterized the deletions detected by multiplex ligation‐dependent probe amplification (MLPA) in the gene MECP2 of 21 RTT patients. Breakpoints were delineated by DNA‐qPCR until the amplification of the deleted allele by long‐PCR was possible.

Results

This methodology enabled us to characterize deletions ranging from 1,235 bp to 85 kb, confirming the partial or total deletion of the MECP2 gene in all these patients. Additionally, our cases support the evidence claiming that most of these breakpoints occur in some restricted regions of the MECP2 gene.

Conclusion

These molecular data together with the clinical information enable us to propose a genotype–phenotype correlation, which is essential for providing genetic counseling.

An evolving view of copy number variants

Copy number variants (CNVs) are regions of the genome that vary in integer copy number. CNVs, which comprise both amplifications and deletions of DNA sequence, have been identified across all domains of life, from bacteria and archaea to plants and animals. CNVs are an important source of genetic diversity, and can drive rapid adaptive evolution and progression of heritable and somatic human diseases, such as cancer. However, despite their evolutionary importance and clinical relevance, CNVs remain understudied compared to single-nucleotide variants (SNVs). This is a consequence of the inherent difficulties in detecting CNVs at low-to-intermediate frequencies in heterogeneous populations of cells. Here, we discuss molecular methods used to detect CNVs, the limitations associated with using these techniques, and the application of new and emerging technologies that present solutions to these challenges. The goal of this short review and perspective is to highlight aspects of CNV biology that are understudied and define avenues for further research that address specific gaps in our knowledge of these complex alleles. We describe our recently developed method for CNV detection in which a fluorescent gene functions as a single-cell CNV reporter and present key findings from our evolution experiments in Saccharomyces cerevisiae. Using a CNV reporter, we found that CNVs are generated at a high rate and undergo selection with predictable dynamics across independently evolving replicate populations. Many CNVs appear to be generated through DNA replication-based processes that are mediated by the presence of short, interrupted, inverted-repeat sequences. Our results have important implications for the role of CNVs in evolutionary processes and the molecular mechanisms that underlie CNV formation. We discuss the possible extension of our method to other applications, including tracking the dynamics of CNVs in models of human tumors.

Epigenetic Components of Myalgic Encephalomyelitis/Chronic Fatigue Syndrome Uncover Potential Transposable Element Activation

PURPOSE

Studies to determine epigenetic changes associated with myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) remain scarce; however, current evidence clearly shows that methylation patterns of genomic DNA and noncoding RNA profiles of immune cells differ between patients and healthy subjects, suggesting an active role of these epigenetic mechanisms in the disease. The present study compares and contrasts the available ME/CFS epigenetic data in an effort to evidence overlapping pathways capable of explaining at least some of the dysfunctional immune parameters linked to this disease.

METHODS

A systematic search of the literature evaluating the ME/CFS epigenome landscape was performed following the Preferred Reporting Items for Systematic Reviews and Meta-Analyses criteria. Differential DNA methylation and noncoding RNA differential expression patterns associated with ME/CFS were used to screen for the presence of transposable elements using the Dfam browser, a search program nurtured with the Repbase repetitive sequence database and the RepeatMasker annotation tool.

FINDINGS

Unexpectedly, particular associations of transposable elements and ME/CFS epigenetic hallmarks were uncovered. A model for the disease emerged involving transcriptional induction of endogenous dormant transposons and structured cellular RNA interactions, triggering the activation of the innate immune system without a concomitant active infection.

IMPLICATIONS

Repetitive sequence filters (ie, RepeatMasker) should be avoided when analyzing transcriptomic data to assess the potential participation of repetitive sequences ("junk repetitive DNA"), representing >45% of the human genome, in the onset and evolution of ME/CFS. In addition, transposable element screenings aimed at designing cost-effective, focused empirical assays that can confirm or disprove the suspected involvement of transposon transcriptional activation in this disease, following the pilot strategy presented here, will require databases gathering large ME/CFS epigenetic datasets.

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

Background

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

Methods

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

Results

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

Conclusions

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

Electronic supplementary material

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

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.

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.