SVA retrotransposon insertion-associated deletion represents a novel mutational mechanism underlying large genomic copy number changes with non-recurrent breakpoints

Long-read sequencing identifies an SVA_D retrotransposon insertion deep within the intron of ATP7A as a novel cause of occipital horn syndrome

BACKGROUND

SINE-VNTR-Alu (SVA) retrotransposons move from one genomic location to another in a 'copy-and-paste' manner. They continue to move actively and cause monogenic diseases through various mechanisms. Currently, disease-causing SVA retrotransposons are classified into human-specific young SVA_E or SVA_F subfamilies. In this study, we identified an evolutionarily old SVA_D retrotransposon as a novel cause of occipital horn syndrome (OHS). OHS is an X-linked, copper metabolism disorder caused by dysfunction of the copper transporter, ATP7A.

METHODS

We investigated a 16-year-old boy with OHS whose pathogenic variant could not be detected via routine molecular genetic analyses.

RESULTS

A 2.8 kb insertion was detected deep within the intron of the patient's ATP7A gene. This insertion caused aberrant mRNA splicing activated by a new donor splice site located within it. Long-read circular consensus sequencing enabled us to accurately read the entire insertion sequence, which contained highly repetitive and GC-rich segments. Consequently, the insertion was identified as an SVA_D retrotransposon. Antisense oligonucleotides (AOs) targeting the new splice site restored the expression of normal transcripts and functional ATP7A proteins. AO treatment alleviated excessive accumulation of copper in patient fibroblasts in a dose-dependent manner. Pedigree analysis revealed that the retrotransposon had moved into the OHS-causing position two generations ago.

CONCLUSION

This is the first report of a human monogenic disease caused by the SVA_D retrotransposon. The fact that the evolutionarily old SVA_D is still actively transposed, leading to increased copy numbers may make a notable impact on rare genetic disease research.

ZNF91 is an endogenous repressor of the molecular phenotype associated with X-linked dystonia-parkinsonism (XDP)

X-linked dystonia-parkinsonism (XDP) is a severe neurodegenerative disorder resulting from an inherited intronic SINE-Alu-VNTR (SVA) retrotransposon in the TAF1 gene that causes dysregulation of TAF1 transcription. The specific mechanism underlying this dysregulation remains unclear, but it is hypothesized to involve the formation of G-quadruplexes (G4) structures within the XDP-SVA that impede transcription. In this study, we show that ZNF91, a critical repressor of SVA retrotransposons, specifically binds to G4-forming DNA sequences. Further, we found that genetic deletion of ZNF91 exacerbates the molecular phenotype associated with the XDP-SVA insertion in patient cells, while no difference was observed when ZNF91 was deleted from isogenic control cells. Additionally, we observed a significant age-related reduction in ZNF91 expression in whole blood and brain, indicating a progressive loss of repression of the XDP-SVA in XDP. These findings indicate that ZNF91 plays a crucial role in controlling the molecular phenotype associated with XDP. Since ZNF91 binds to G4-forming DNA sequences in SVAs, this suggests that interactions between ZNF91 and G4-forming sequences in the XDP-SVA minimize the severity of the molecular phenotype. Our results showing that ZNF91 expression levels significantly decrease with age provide a potential explanation for the age-related progressive neurodegenerative character of XDP. Collectively, our study provides important insights into the protective role of ZNF91 in XDP pathogenesis and suggests that restoring ZNF91 expression, destabilization of G4s, or targeted repression of the XDP-SVA could be future therapeutic strategies to prevent or treat XDP.

NanoRanger enables rapid single-base-pair resolution of genomic disorders

BACKGROUND

Delineating base-resolution breakpoints of complex rearrangements is crucial for an accurate clinical understanding of pathogenic variants and for carrier screening within family networks or the broader population. However, despite advances in genetic testing using short-read sequencing (SRS), this task remains costly and challenging.

METHODS

This study addresses the challenges of resolving missing disease-causing breakpoints in complex genomic disorders with suspected homozygous rearrangements by employing multiple long-read sequencing (LRS) strategies, including a novel and efficient strategy named nanopore-based rapid acquisition of neighboring genomic regions (NanoRanger). NanoRanger does not require large amounts of ultrahigh-molecular-weight DNA and stands out for its ease of use and rapid acquisition of large genomic regions of interest with deep coverage.

FINDINGS

We describe a cohort of 16 familial cases, each harboring homozygous rearrangements that defied breakpoint determination by SRS and optical genome mapping (OGM). NanoRanger identified the breakpoints with single-base-pair resolution, enabling accurate determination of the carrier status of unaffected family members as well as the founder nature of these genomic lesions and their frequency in the local population. The resolved breakpoints revealed that repetitive DNA, gene regulatory elements, and transcription activity contribute to genome instability in these novel recessive rearrangements.

CONCLUSIONS

Our data suggest that NanoRanger greatly improves the success rate of resolving base-resolution breakpoints of complex genomic disorders and expands access to LRS for the benefit of patients with Mendelian disorders.

FUNDING

M.L. is supported by KAUST Baseline Award no. BAS/1/1080-01-01 and KAUST Research Translation Fund Award no. REI/1/4742-01.

Retrotransposons in embryogenesis and neurodevelopment

Retrotransposable elements (RTEs) are genetic elements that can replicate and insert new copies into different genomic locations. RTEs have long been identified as 'parasitic genes', as their mobilization can cause mutations, DNA damage, and inflammation. Interestingly, high levels of retrotransposon activation are observed in early embryogenesis and neurodevelopment, suggesting that RTEs may possess functional roles during these stages of development. Recent studies demonstrate that RTEs can function as transcriptional regulatory elements through mechanisms such as chromatin organization and noncoding RNAs. It is clear, however, that RTE expression and activity must be restrained at some level during development, since overactivation of RTEs during neurodevelopment is associated with several developmental disorders. Further investigation is needed to understand the importance of RTE expression and activity during neurodevelopment and the balance between RTE-regulated development and RTE-mediated pathogenesis.

Mini-heterochromatin domains constrain the cis-regulatory impact of SVA transposons in human brain development and disease

SVA (SINE (short interspersed nuclear element)-VNTR (variable number of tandem repeats)-Alu) retrotransposons remain active in humans and contribute to individual genetic variation. Polymorphic SVA alleles harbor gene regulatory potential and can cause genetic disease. However, how SVA insertions are controlled and functionally impact human disease is unknown. Here we dissect the epigenetic regulation and influence of SVAs in cellular models of X-linked dystonia parkinsonism (XDP), a neurodegenerative disorder caused by an SVA insertion at the TAF1 locus. We demonstrate that the KRAB zinc finger protein ZNF91 establishes H3K9me3 and DNA methylation over SVAs, including polymorphic alleles, in human neural progenitor cells. The resulting mini-heterochromatin domains attenuate the cis-regulatory impact of SVAs. This is critical for XDP pathology; removal of local heterochromatin severely aggravates the XDP molecular phenotype, resulting in increased TAF1 intron retention and reduced expression. Our results provide unique mechanistic insights into how human polymorphic transposon insertions are recognized and how their regulatory impact is constrained by an innate epigenetic defense system.

Genetic/epigenetic effects in NF1 microdeletion syndrome: beyond the haploinsufficiency, looking at the contribution of not deleted genes

NF1 microdeletion syndrome, accounting for 5-11% of NF1 patients, is caused by a deletion in the NF1 region and it is generally characterized by a severe phenotype. Although 70% of NF1 microdeletion patients presents the same 1.4 Mb type-I deletion, some patients may show additional clinical features. Therefore, the contribution of several pathogenic mechanisms, besides haploinsufficiency of some genes within the deletion interval, is expected and needs to be defined. We investigated an altered expression of deletion flanking genes by qPCR in patients with type-1 NF1 deletion, compared to healthy donors, possibly contributing to the clinical traits of NF1 microdeletion syndrome. In addition, the 1.4-Mb deletion leads to changes in the 3D chromatin structure in the 17q11.2 region. Specifically, this deletion alters DNA-DNA interactions in the regions flanking the breakpoints, as demonstrated by our 4C-seq analysis. This alteration likely causes position effect on the expression of deletion flanking genes.Interestingly, 4C-seq analysis revealed that in microdeletion patients, an interaction was established between the RHOT1 promoter and the SLC6A4 gene, which showed increased expression. We performed NGS on putative modifier genes, and identified two "likely pathogenic" rare variants in RAS pathway, possibly contributing to incidental phenotypic features.This study provides new insights into understanding the pathogenesis of NF1 microdeletion syndrome and suggests a novel pathomechanism that contributes to the expression phenotype in addition to haploinsufficiency of genes located within the deletion.This is a pivotal approach that can be applied to unravel microdeletion syndromes, improving precision medicine, prognosis and patients' follow-up.

Contribution of de novo retroelements to birth defects and childhood cancers

Insertion of active retroelements-L1s, Alus, and SVAs-can disrupt proper genome function and lead to various disorders including cancer. However, the role of de novo retroelements (DNRTs) in birth defects and childhood cancers has not been well characterized due to the lack of adequate data and efficient computational tools. Here, we examine whole-genome sequencing data of 3,244 trios from 12 birth defect and childhood cancer cohorts in the Gabriella Miller Kids First Pediatric Research Program. Using an improved version of our tool xTea (x-Transposable element analyzer) that incorporates a deep-learning module, we identified 162 DNRTs, as well as 2 pseudogene insertions. Several variants are likely to be causal, such as a de novo Alu insertion that led to the ablation of a whole exon in the NF1 gene in a proband with brain tumor. We observe a high de novo SVA insertion burden in both high-intolerance loss-of-function genes and exons as well as more frequent de novo Alu insertions of paternal origin. We also identify potential mosaic DNRTs from embryonic stages. Our study reveals the important roles of DNRTs in causing birth defects and predisposition to childhood cancers.

Correlation between large rearrangements and patient phenotypes in NF1 deletion syndrome: an update and review

About 5-10% of neurofibromatosis type 1 (NF1) patients exhibit large genomic germline deletions that remove the NF1 gene and its flanking regions. The most frequent NF1 large deletion is 1.4 Mb, resulting from homologous recombination between two low copy repeats. This "type-1" deletion is associated with a severe clinical phenotype in NF1 patients, with several phenotypic manifestations including learning disability, a much earlier development of cutaneous neurofibromas, an increased tumour risk, and cardiovascular malformations. NF1 adjacent co-deleted genes could act as modifier loci for the specific clinical manifestations observed in deleted NF1 patients. Furthermore, other genetic modifiers (such as CNVs) not located at the NF1 locus could also modulate the phenotype observed in patients with large deletions. In this study, we analysed 22 NF1 deletion patients by genome-wide array-CGH with the aim (1) to correlate deletion length to observed phenotypic features and their severity in NF1 deletion syndrome, and (2) to identify whether the deletion phenotype could also be modulated by copy number variations elsewhere in the genome. We then review the role of co-deleted genes in the 1.4 Mb interval of type-1 deletions, and their possible implication in the main clinical features observed in this high-risk group of NF1 patients.

Small polymorphisms are a source of ancestral bias in structural variant breakpoint placement

High-quality genome assemblies and sophisticated algorithms have increased sensitivity for a wide range of variant types, and breakpoint accuracy for structural variants (SVs, ≥50 bp) has improved to near base pair precision. Despite these advances, many SV breakpoint locations are subject to systematic bias affecting variant representation. To understand why SV breakpoints are inconsistent across samples, we reanalyzed 64 phased haplotypes constructed from long-read assemblies released by the Human Genome Structural Variation Consortium (HGSVC). We identify 882 SV insertions and 180 SV deletions with variable breakpoints not anchored in tandem repeats (TRs) or segmental duplications (SDs). SVs called from aligned sequencing reads increase breakpoint disagreements by 2×-16×. Sequence accuracy had a minimal impact on breakpoints, but we observe a strong effect of ancestry. We confirm that SNP and indel polymorphisms are enriched at shifted breakpoints and are also absent from variant callsets. Breakpoint homology increases the likelihood of imprecise SV calls and the distance they are shifted, and tandem duplications are the most heavily affected SVs. Because graph genome methods normalize SV calls across samples, we investigated graphs generated by two different methods and find the resulting breakpoints are subject to other technical biases affecting breakpoint accuracy. The breakpoint inconsistencies we characterize affect ∼5% of the SVs called in a human genome and can impact variant interpretation and annotation. These limitations underscore a need for algorithm development to improve SV databases, mitigate the impact of ancestry on breakpoints, and increase the value of callsets for investigating breakpoint features.

Droplet Digital PCR for Fast and Accurate Characterization of NF1 Locus Deletions: Confirmation of the Predominant Maternal Origin of Type-1 Deletions

Neurofibromatosis type-1 is a genetic disorder caused by loss-of-function variants in the tumor-suppressor NF1. Approximately 4% to 11% of neurofibromatosis type-1 patients have a NF1 locus complete deletion resulting from nonallelic homologous recombination between low copy repeats. Codeleted genes probably account for the more severe phenotype observed in NF1-deleted patients. This genotype-phenotype correlation highlights the need for a detailed molecular description. A droplet digital PCR (ddPCR) set along the NF1 locus was designed to delimitate the three recurrent NF1 deletion breakpoints. The ddPCR was tested in 121 samples from nonrelated NF1-deleted patients. Classification based on ddPCR versus multiplex ligation-dependent probe amplification (MLPA) was compared. In addition, microsatellites were analyzed to identify parental origin of deletions. ddPCR identified 77 type-1 (64%), 20 type-2 (16%), 7 type-3 (6%), and 17 atypical deletions (14%). The results were comparable with MLPA, except for three atypical deletions misclassified as type-2 using MLPA, for which the SUZ12 gene was not deleted. A significant maternal bias (25 of 30) in the origin of deletions was identified. This study proposes a fast and efficient ddPCR quantification to allow fine NF1 deletion classification. It indicates that ddPCR can be implemented easily into routine diagnosis to complement the techniques dedicated to NF1 point variant identification. This new tool may help unravel the genetic basis conditioning phenotypic variability in NF1-deleted patients and offer tailored genetic counseling.

Nf1 in heart development: a potential causative gene for congenital heart disease: a narrative review

Congenital heart disease is the most frequent congenital disorder, affecting a significant number of live births. Gaining insights into its genetic etiology could lead to a deeper understanding of this condition. Although the Nf1 gene has been identified as a potential causative gene, its role in congenital heart disease has not been thoroughly clarified. We searched and summarized evidence from cohort-based and experimental studies on the issue of Nf1 and heart development in congenital heart diseases from various databases. Available evidence demonstrates a correlation between Nf1 and congenital heart diseases, mainly pulmonary valvar stenosis. The mechanism underlying this correlation may involve dysregulation of epithelial-mesenchymal transition (EMT). The Nf1 gene affects the EMT process via multiple pathways, including directly regulating the expression of EMT-related transcription factors and indirectly regulating the EMT process by regulating the MAPK pathway. This narrative review provides a comprehensive account of the Nf1 involvement in heart development and congenital cardiovascular diseases in terms of epidemiology and potential mechanisms. RAS signaling may contribute to congenital heart disease independently or in cooperation with other signaling pathways. Efficient management of both NF1 and cardiovascular disease patients would benefit from further research into these issues.

Exome sequencing identified rare recurrent copy number variants and hereditary breast cancer susceptibility

Copy number variants (CNVs) are a major source of genetic variation and can disrupt genes or affect gene dosage. They are known to be causal or underlie predisposition to various diseases. However, the role of CNVs in inherited breast cancer susceptibility has not been thoroughly investigated. To address this, we performed whole-exome sequencing based analysis of rare CNVs in 98 high-risk Northern Finnish breast cancer cases. After filtering, selected candidate alleles were validated and characterized with a combination of orthogonal methods, including PCR-based approaches, optical genome mapping and long-read sequencing. This revealed three recurrent alterations: a 31 kb deletion co-occurring with a retrotransposon insertion (delins) in RAD52, a 13.4 kb deletion in HSD17B14 and a 64 kb partial duplication of RAD51C. Notably, all these genes encode proteins involved in pathways previously identified as essential for breast cancer development. Variants were genotyped in geographically matched cases and controls (altogether 278 hereditary and 1983 unselected breast cancer cases, and 1229 controls). The RAD52 delins and HSD17B14 deletion both showed significant enrichment among cases with indications of hereditary disease susceptibility. RAD52 delins was identified in 7/278 cases (2.5%, P = 0.034, OR = 2.86, 95% CI = 1.10-7.45) and HSD17B14 deletion in 8/278 cases (2.9%, P = 0.014, OR = 3.28, 95% CI = 1.31-8.23), the frequency of both variants in the controls being 11/1229 (0.9%). This suggests a role for RAD52 and HSD17B14 in hereditary breast cancer susceptibility. The RAD51C duplication was very rare, identified only in 2/278 of hereditary cases and 2/1229 controls (P = 0.157, OR = 4.45, 95% CI = 0.62-31.70). The identification of recurrent CNVs in these genes, and especially the relatively high frequency of RAD52 and HSD17B14 alterations in the Finnish population, highlights the importance of studying CNVs alongside single nucleotide variants when searching for genetic factors underlying hereditary disease predisposition.

Transposable element insertions in 1000 Swedish individuals

The majority of rare diseases are genetic, and regardless of advanced high-throughput genomics-based investigations, 60% of patients remain undiagnosed. A major factor limiting our ability to identify disease-causing alterations is a poor understanding of the morbid and normal human genome. A major genomic contributor of which function and distribution remain largely unstudied are the transposable elements (TE), which constitute 50% of our genome. Here we aim to resolve this knowledge gap and increase the diagnostic yield of rare disease patients investigated with clinical genome sequencing. To this end we characterized TE insertions in 1000 Swedish individuals from the SweGen dataset and 2504 individuals from the 1000 Genomes Project (1KGP), creating seven population-specific TE insertion databases. Of note, 66% of TE insertions in SweGen were present at >1% in the 1KGP databases, proving that most insertions are common across populations. Focusing on the rare TE insertions, we show that even though ~0.7% of those insertions affect protein coding genes, they rarely affect known disease casing genes (<0.1%). Finally, we applied a TE insertion identification workflow on two clinical cases where disease causing TE insertions were suspected and could verify the presence of pathogenic TE insertions in both. Altogether we demonstrate the importance of TE insertion detection and highlight possible clinical implications in rare disease diagnostics.

Deregulated expression of polycomb repressive complex 2 target genes in a NF1 patient with microdeletion generating the RNF135-SUZ12 chimeric gene

Neurofibromatosis type I (NF1) microdeletion syndrome, accounting for 5-11% of NF1 patients, is caused by the heterozygous deletion of NF1 and a variable number of flanking genes in the 17q11.2 region. This syndrome is characterized by more severe symptoms than those shown by patients with intragenic NF1 mutation and by variable expressivity, which is not fully explained by the haploinsufficiency of the genes included in the deletions. We here reevaluate an 8-year-old NF1 patient, who carries an atypical deletion generating the RNF135-SUZ12 chimeric gene, previously described when he was 3 years old. As the patient has developed multiple cutaneous/subcutaneous neurofibromas over the past 5 years, we hypothesized a role of RNF135-SUZ12 chimeric gene in the onset of the patient's tumor phenotype. Interestingly, SUZ12 is generally lost or disrupted in NF1 microdeletion syndrome and frequently associated to cancer as RNF135. Expression analysis confirmed the presence of the chimeric gene transcript and revealed hypo-expression of five out of the seven analyzed target genes of the polycomb repressive complex 2 (PRC2), to which SUZ12 belongs, in the patient's peripheral blood, indicating a higher transcriptional repression activity mediated by PRC2. Furthermore, decreased expression of tumor suppressor gene TP53, which is targeted by RNF135, was detected. These results suggest that RNF135-SUZ12 chimera may acquire a gain of function, compared with SUZ12 wild type in the PRC2 complex, and a loss of function relative to RNF135 wild type. Both events may have a role in the early onset of the patient's neurofibromas.

Small allelic variants are a source of ancestral bias in structural variant breakpoint placement

High-quality genome assemblies and sophisticated algorithms have increased sensitivity for a wide range of variant types, and breakpoint accuracy for structural variants (SVs, ≥ 50 bp) has improved to near basepair precision. Despite these advances, many SVs in unique regions of the genome are subject to systematic bias that affects breakpoint location. This ambiguity leads to less accurate variant comparisons across samples, and it obscures true breakpoint features needed for mechanistic inferences. To understand why SVs are not consistently placed, we reanalyzed 64 phased haplotypes constructed from long-read assemblies released by the Human Genome Structural Variation Consortium (HGSVC). We identified variable breakpoints for 882 SV insertions and 180 SV deletions not anchored in tandem repeats (TRs) or segmental duplications (SDs). While this is unexpectedly high for genome assemblies in unique loci, we find read-based callsets from the same sequencing data yielded 1,566 insertions and 986 deletions with inconsistent breakpoints also not anchored in TRs or SDs. When we investigated causes for breakpoint inaccuracy, we found sequence and assembly errors had minimal impact, but we observed a strong effect of ancestry. We confirmed that polymorphic mismatches and small indels are enriched at shifted breakpoints and that these polymorphisms are generally lost when breakpoints shift. Long tracts of homology, such as SVs mediated by transposable elements, increase the likelihood of imprecise SV calls and the distance they are shifted. Tandem Duplication (TD) breakpoints are the most heavily affected SV class with 14% of TDs placed at different locations across haplotypes. While graph genome methods normalize SV calls across many samples, the resulting breakpoints are sometimes incorrect, highlighting a need to tune graph methods for breakpoint accuracy. The breakpoint inconsistencies we characterize collectively affect ~5% of the SVs called in a human genome and underscore a need for algorithm development to improve SV databases, mitigate the impact of ancestry on breakpoint placement, and increase the value of callsets for investigating mutational processes.

Hominoid SVA-lncRNA AK057321 targets human-specific SVA retrotransposons in SCN8A and CDK5RAP2 to initiate neuronal maturation

SINE-VNTR-Alu (SVA) retrotransposons arose and expanded in the genome of hominoid primates concurrent with the slowing of brain maturation. We report genes with intronic SVA transposons are enriched for neurodevelopmental disease and transcribed into long non-coding SVA-lncRNAs. Human-specific SVAs in microcephaly CDK5RAP2 and epilepsy SCN8A gene introns repress their expression via transcription factor ZNF91 to delay neuronal maturation. Deleting the SVA in CDK5RAP2 initiates multi-dimensional and in SCN8A selective sodium current neuronal maturation by upregulating these genes. SVA-lncRNA AK057321 forms RNA:DNA heteroduplexes with the genomic SVAs and upregulates these genes to initiate neuronal maturation. SVA-lncRNA AK057321 also promotes species-specific cortex and cerebellum-enriched expression upregulating human genes with intronic SVAs (e.g., HTT, CHAF1B and KCNJ6) but not mouse orthologs. The diversity of neuronal genes with intronic SVAs suggest this hominoid-specific SVA transposon-based gene regulatory mechanism may act at multiple steps to specialize and achieve neoteny of the human brain.

Retrotransposon insertion as a novel mutational cause of spinal muscular atrophy

Spinal muscular atrophy (SMA) is an autosomal recessive neuromuscular disorder resulting from biallelic alterations of the SMN1 gene: deletion, gene conversion or, in rare cases, intragenic variants. The disease severity is mainly influenced by the copy number of SMN2, a nearly identical gene, which produces only low amounts of full-length (FL) mRNA. Here we describe the first example of retrotransposon insertion as a pathogenic SMN1 mutational event. The 50-year-old patient is clinically affected by SMA type III with a diagnostic odyssey spanning nearly 30 years. Despite a mild disease course, he carries a single SMN2 copy. Using Exome Sequencing and Sanger sequencing, we characterized a SINE-VNTR-Alu (SVA) type F retrotransposon inserted in SMN1 intron 7. Using RT-PCR and RNASeq experiments on lymphoblastoid cell lines, we documented the dramatic decrease of FL transcript production in the patient compared to subjects with the same SMN1 and SMN2 copy number, thus validating the pathogenicity of this SVA insertion. We described the mutant FL-SMN1-SVA transcript characterized by exon extension and showed that it is subject to degradation by nonsense-mediated mRNA decay. The stability of the SMN-SVA protein may explain the mild course of the disease. This observation exemplifies the role of retrotransposons in human genetic disorders.

The landscape of aging

Aging is characterized by a progressive deterioration of physiological integrity, leading to impaired functional ability and ultimately increased susceptibility to death. It is a major risk factor for chronic human diseases, including cardiovascular disease, diabetes, neurological degeneration, and cancer. Therefore, the growing emphasis on "healthy aging" raises a series of important questions in life and social sciences. In recent years, there has been unprecedented progress in aging research, particularly the discovery that the rate of aging is at least partly controlled by evolutionarily conserved genetic pathways and biological processes. In an attempt to bring full-fledged understanding to both the aging process and age-associated diseases, we review the descriptive, conceptual, and interventive aspects of the landscape of aging composed of a number of layers at the cellular, tissue, organ, organ system, and organismal levels.

An apple somatic mutation of delayed fruit maturation date is primarily caused by a retrotransposon insertion-associated large deletion

Somatic mutations may alter important traits in tree fruits, such as fruit color, size and maturation date. Autumn Gala (AGala), a somatic mutation from apple cultivar Gala, matures 4 weeks later than Gala. To understand the mechanisms underlying the delayed maturation, RNA-seq analyses were conducted with fruit sampled at 13 (Gala) and 16 (AGala) time-points during their growth and development. Weighted gene co-expression network analysis (WGCNA) of 23 372 differentially expressed genes resulted in 25 WGCNA modules. Of these, modules 1 (r = -0.98, P = 2E-21) and 2 (r = -0.52, P = 0.004), which were suppressed in AGala, were correlated with fruit maturation date. Surprisingly, 77 of the 152 member genes in module 1 were harbored in a 2.8-Mb genomic region on chromosome 6 that was deleted and replaced by a 10.7-kb gypsy-like retrotransposon (Gy-36) from chromosome 7 in AGala. Among the 77 member genes, MdACT7 was the most suppressed (by 10.5-fold) in AGala due to a disruptive 2.5-kb insertion in coding sequence. Moreover, MdACT7 is the exclusive apple counterpart of Arabidopsis ACT7 known of essential roles in plant development, and the functional allele MdACT7, which was lost to the deletion in AGala, was associated with early fruit maturation in 268 apple accessions. Overexpressing alleles MdACT7 and Mdact7 in an Arabidopsis act7 line showed that MdACT7 largely rescued its stunted growth and delayed initial flowering while Mdact7 did not. Therefore, the 2.8-Mb hemizygous deletion is largely genetically causal for fruit maturation delay in AGala, and the total loss of MdACT7 might have contributed to the phenotype.

Long-Read Sequencing Identifies the First Retrotransposon Insertion and Resolves Structural Variants Causing Antithrombin Deficiency

The identification of inherited antithrombin deficiency (ATD) is critical to prevent potentially life-threatening thrombotic events. Causal variants in SERPINC1 are identified for up to 70% of cases, the majority being single-nucleotide variants and indels. The detection and characterization of structural variants (SVs) in ATD remain challenging due to the high number of repetitive elements in SERPINC1. Here, we performed long-read whole-genome sequencing on 10 familial and 9 singleton cases with type I ATD proven by functional and antigen assays, who were selected from a cohort of 340 patients with this rare disorder because genetic analyses were either negative, ambiguous, or not fully characterized. We developed an analysis workflow to identify disease-associated SVs. This approach resolved, independently of its size or type, all eight SVs detected by multiple ligation-dependent probe amplification, and identified for the first time a complex rearrangement previously misclassified as a deletion. Remarkably, we identified the mechanism explaining ATD in 2 out of 11 cases with previous unknown defect: the insertion of a novel 2.4 kb SINE-VNTR-Alu retroelement, which was characterized by de novo assembly and verified by specific polymerase chain reaction amplification and sequencing in the probands and affected relatives. The nucleotide-level resolution achieved for all SVs allowed breakpoint analysis, which revealed repetitive elements and microhomologies supporting a common replication-based mechanism for all the SVs. Our study underscores the utility of long-read sequencing technology as a complementary method to identify, characterize, and unveil the molecular mechanism of disease-causing SVs involved in ATD, and enlarges the catalogue of genetic disorders caused by retrotransposon insertions.

An ancestral genomic sequence that serves as a nucleation site for de novo gene birth

The process of gene birth is of major interest with current excitement concerning de novo gene formation. We report a new and different mechanism of de novo gene birth based on the finding and the characteristics of a short non-coding sequence situated between two protein genes, termed a spacer sequence. This non-coding sequence is present in genomes of Mus musculus, the house mouse and Philippine tarsier, a primitive ancestral primate. The ancestral sequence is highly conserved during primate evolution with certain base pairs totally invariant from mouse to humans. By following the birth of the sequence of human lincRNA BCRP3 (BCR activator of RhoGEF and GTPase 3 pseudogene) during primate evolution, we find diverse genes, long non-coding RNA and protein genes (and sequences that do not appear to encode a gene) that all stem from the 3’ end of the spacer, and all begin with a similar sequence. During primate evolution, part of the BCRP3 sequence initially formed in the Old World Monkeys and developed into different primate genes before evolving into the BCRP3 gene in humans. The gene developmental process consists of the initiation of DNA synthesis at spacer 3’ ends, addition of a complex of tandem transposable elements and the addition of a segment of another gene. The findings support the concept of the spacer sequence as a starting site for DNA synthesis that leads to formation of different genes with the addition of other sequences. These data suggest a new process of de novo gene birth.

Mechanisms of disease-associated SINE-VNTR-Alus

SINE-VNTR-Alus (SVAs) are the youngest retrotransposon family in the human genome. Their ongoing mobilization has generated genetic variation within the human population. At least 24 insertions to date, detailed in this review, have been associated with disease. The predominant mechanisms through which this occurs are alterations to normal splicing patterns, exonic insertions causing loss-of-function mutations, and large genomic deletions. Dissecting the functional impact of these SVAs and the mechanism through which they cause disease provides insight into the consequences of their presence in the genome and how these elements could influence phenotypes. Many of these disease-associated SVAs have been difficult to characterize and would not have been identified through routine analyses. However, the number identified has increased in recent years as DNA and RNA sequencing data became more widely available. Therefore, as the search for complex structural variation in disease continues, it is likely to yield further disease-causing SVA insertions.

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

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

Atypical NF1 Microdeletions: Challenges and Opportunities for Genotype/Phenotype Correlations in Patients with Large NF1 Deletions

Patients with neurofibromatosis type 1 (NF1) and type 1 NF1 deletions often exhibit more severe clinical manifestations than patients with intragenic NF1 gene mutations, including facial dysmorphic features, overgrowth, severe global developmental delay, severe autistic symptoms and considerably reduced cognitive abilities, all of which are detectable from a very young age. Type 1 NF1 deletions encompass 1.4 Mb and are associated with the loss of 14 protein-coding genes, including NF1 and SUZ12. Atypical NF1 deletions, which do not encompass all 14 protein-coding genes located within the type 1 NF1 deletion region, have the potential to contribute to the delineation of the genotype/phenotype relationship in patients with NF1 microdeletions. Here, we review all atypical NF1 deletions reported to date as well as the clinical phenotype observed in the patients concerned. We compare these findings with those of a newly identified atypical NF1 deletion of 698 kb which, in addition to the NF1 gene, includes five genes located centromeric to NF1. The atypical NF1 deletion in this patient does not include the SUZ12 gene but does encompass CRLF3. Comparative analysis of such atypical NF1 deletions suggests that SUZ12 hemizygosity is likely to contribute significantly to the reduced cognitive abilities, severe global developmental delay and facial dysmorphisms observed in patients with type 1 NF1 deletions.

Classification of NF1 microdeletions and its importance for establishing genotype/phenotype correlations in patients with NF1 microdeletions

An estimated 5–11% of patients with neurofibromatosis type-1 (NF1) harbour large deletions encompassing the NF1 gene and flanking regions. These NF1 microdeletions are subclassified into type 1, 2, 3 and atypical deletions which are distinguishable from each other by their extent and by the number of genes included within the deletion regions as well as the frequency of mosaicism with normal cells. Most common are type-1 NF1 deletions which encompass 1.4-Mb and 14 protein-coding genes. Type-1 deletions are frequently associated with overgrowth, global developmental delay, cognitive disability and dysmorphic facial features which are uncommon in patients with intragenic pathogenic NF1 gene variants. Further, patients with type-1 NF1 deletions frequently exhibit high numbers of neurofibromas and have an increased risk of malignant peripheral nerve sheath tumours. Genes located within the type-1 NF1 microdeletion interval and co-deleted with NF1 are likely to act as modifiers responsible for the severe disease phenotype in patients with NF1 microdeletions, thereby causing the NF1 microdeletion syndrome. Genotype/phenotype correlations in patients with NF1 microdeletions of different lengths are important to identify such modifier genes. However, these correlations are critically dependent upon the accurate characterization of the deletions in terms of their extent. In this review, we outline the utility as well as the shortcomings of multiplex ligation-dependent probe amplification (MLPA) to classify the different types of NF1 microdeletion and indicate the importance of high-resolution microarray analysis for correct classification, a necessary precondition to identify those genes responsible for the NF1 microdeletion syndrome.

Optical genome mapping identifies a germline retrotransposon insertion in SMARCB1 in two siblings with atypical teratoid rhabdoid tumors

In a subset of pediatric cancers, a germline cancer predisposition is highly suspected based on clinical and pathological findings, but genetic evidence is lacking, which hampers genetic counseling and predictive testing in the families involved. We describe a family with two siblings born from healthy parents who were both neonatally diagnosed with atypical teratoid rhabdoid tumor (ATRT). This rare and aggressive pediatric tumor is associated with biallelic inactivation of SMARCB1, and in 30% of the cases, a predisposing germline mutation is involved. Whereas the tumors of both siblings showed loss of expression of SMARCB1 and acquired homozygosity of the locus, whole exome and whole genome sequencing failed to identify germline or somatic SMARCB1 pathogenic mutations. We therefore hypothesized that the insertion of a pathogenic repeat‐rich structure might hamper its detection, and we performed optical genome mapping (OGM) as an alternative strategy to identify structural variation in this locus. Using this approach, an insertion of ~2.8 kb within intron 2 of SMARCB1 was detected. Long‐range PCR covering this region remained unsuccessful, but PacBio HiFi genome sequencing identified this insertion to be a SINE‐VNTR‐Alu, subfamily E (SVA‐E) retrotransposon element, which was present in a mosaic state in the mother. This SVA‐E insertion disrupts correct splicing of the gene, resulting in loss of a functional allele. This case demonstrates the power of OGM and long‐read sequencing to identify genomic variations in high‐risk cancer‐predisposing genes that are refractory to detection with standard techniques, thereby completing the clinical and molecular diagnosis of such complex cases and greatly improving counseling and surveillance of the families involved. © 2021 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland.

Vertebrate Lineages Exhibit Diverse Patterns of Transposable Element Regulation and Expression across Tissues

Transposable elements (TEs) comprise a major fraction of vertebrate genomes, yet little is known about their expression and regulation across tissues, and how this varies across major vertebrate lineages. We present the first comparative analysis integrating TE expression and TE regulatory pathway activity in somatic and gametic tissues for a diverse set of 12 vertebrates. We conduct simultaneous gene and TE expression analyses to characterize patterns of TE expression and TE regulation across vertebrates and examine relationships between these features. We find remarkable variation in the expression of genes involved in TE negative regulation across tissues and species, yet consistently high expression in germline tissues, particularly in testes. Most vertebrates show comparably high levels of TE regulatory pathway activity across gonadal tissues except for mammals, where reduced activity of TE regulatory pathways in ovarian tissues may be the result of lower relative germ cell densities. We also find that all vertebrate lineages examined exhibit remarkably high levels of TE-derived transcripts in somatic and gametic tissues, with recently active TE families showing higher expression in gametic tissues. Although most TE-derived transcripts originate from inactive ancient TE families (and are likely incapable of transposition), such high levels of TE-derived RNA in the cytoplasm may have secondary, unappreciated biological relevance.

Genotype-Phenotype Associations in Patients With Type-1, Type-2, and Atypical NF1 Microdeletions

Neurofibromatosis type 1 is a tumor predisposition syndrome inherited in autosomal dominant manner. Besides the intragenic loss-of-function mutations in NF1 gene, large deletions encompassing the NF1 gene and its flanking regions are responsible for the development of the variable clinical phenotype. These large deletions titled as NF1 microdeletions lead to a more severe clinical phenotype than those observed in patients with intragenic NF1 mutations. Around 5-10% of the cases harbor large deletion and four major types of NF1 microdeletions (type 1, 2, 3 and atypical) have been identified so far. They are distinguishable in term of their size and the location of the breakpoints, by the frequency of somatic mosaicism with normal cells not harboring the deletion and by the number of the affected genes within the deleted region. In our study genotype-phenotype analyses have been performed in 17 mostly pediatric patients with NF1 microdeletion syndrome identified by multiplex ligation-dependent probe amplification after systematic sequencing of the NF1 gene. Confirmation and classification of the NF1 large deletions were performed using array comparative genomic hybridization, where it was feasible. In our patient cohort 70% of the patients possess type-1 deletion, one patient harbors type-2 deletion and 23% of our cases have atypical NF1 deletion. All the atypical deletions identified in this study proved to be novel. One patient with atypical deletion displayed mosaicism. In our study NF1 microdeletion patients presented dysmorphic facial features, macrocephaly, large hands and feet, delayed cognitive development and/or learning difficulties, speech difficulties, overgrowth more often than patients with intragenic NF1 mutations. Moreover, neurobehavior problems, macrocephaly and overgrowth were less frequent in atypical cases compared to type-1 deletion. Proper diagnosis is challenging in certain patients since several clinical manifestations show age-dependency. Large tumor load exhibited more frequently in this type of disorder, therefore better understanding of genotype-phenotype correlations and progress of the disease is essential for individuals suffering from neurofibromatosis to improve the quality of their life. Our study presented additional clinical data related to NF1 microdeletion patients especially for pediatric cases and it contributes to the better understanding of this type of disorder.

How Chaotic Is Genome Chaos?

Simple Summary

Cancer genomes can undergo major restructurings involving many chromosomal locations at key stages in tumor development. This restructuring process has been designated “genome chaos” by some authors. In order to examine how chaotic cancer genome restructuring may be, the cell and molecular processes for DNA restructuring are reviewed. Examination of the action of these processes in various cancers reveals a degree of specificity that indicates genome restructuring may be sufficiently reproducible to enable possible therapies that interrupt tumor progression to more lethal forms.

Abstract

Cancer genomes evolve in a punctuated manner during tumor evolution. Abrupt genome restructuring at key steps in this evolution has been called “genome chaos.” To answer whether widespread genome change is truly chaotic, this review (i) summarizes the limited number of cell and molecular systems that execute genome restructuring, (ii) describes the characteristic signatures of DNA changes that result from activity of those systems, and (iii) examines two cases where genome restructuring is determined to a significant degree by cell type or viral infection. The conclusion is that many restructured cancer genomes display sufficiently unchaotic signatures to identify the cellular systems responsible for major oncogenic transitions, thereby identifying possible targets for therapies to inhibit tumor progression to greater aggressiveness.

Clinical characterization of children and adolescents with NF1 microdeletions

PURPOSE

An estimated 5-11% of patients with neurofibromatosis type 1 (NF1) harbour NF1 microdeletions encompassing the NF1 gene and its flanking regions. The purpose of this study was to evaluate the clinical phenotype in children and adolescents with NF1 microdeletions.

METHODS

We retrospectively analysed 30 children and adolescents with NF1 microdeletions pertaining to externally visible neurofibromas. The internal tumour load was determined by volumetry of whole-body magnetic resonance imaging (MRI) in 20 children and adolescents with NF1 microdeletions. Furthermore, the prevalence of global developmental delay, autism spectrum disorder and attention deficit hyperactivity disorder (ADHD) were evaluated.

RESULTS

Children and adolescents with NF1 microdeletions had significantly more often cutaneous, subcutaneous and externally visible plexiform neurofibromas than age-matched patients with intragenic NF1 mutations. Internal neurofibromas were detected in all 20 children and adolescents with NF1 microdeletions analysed by whole-body MRI. By contrast, only 17 (61%) of 28 age-matched NF1 patients without microdeletions had internal tumours. The total internal tumour load was significantly higher in NF1 microdeletion patients than in NF1 patients without microdeletions. Global developmental delay was observed in 28 (93%) of 30 children with NF1 microdeletions investigated. The mean full-scale intelligence quotient in our patient group was 77.7 which is significantly lower than that of patients with intragenic NF1 mutations. ADHD was diagnosed in 15 (88%) of 17 children and adolescents with NF1 microdeletion. Furthermore, 17 (71%) of the 24 patients investigated had T-scores ≥ 60 up to 75, indicative of mild to moderate autistic symptoms, which are consequently significantly more frequent in patients with NF1 microdeletions than in the general NF1 population. Also, the mean total T-score was significantly higher in patients with NF1 microdeletions than in the general NF1 population.

CONCLUSION

Our findings indicate that already at a very young age, NF1 microdeletions patients frequently exhibit a severe disease manifestation which requires specialized long-term clinical care.

Genetic basis of neurofibromatosis type 1 and related conditions, including mosaicism

INTRODUCTION

Neurofibromatosis type 1 (NF1) is a frequent autosomal dominant disorder characterised by café-au-lait maculae (CALM), skinfold freckling, iris Lisch nodules and benign peripheral nerve sheath tumours (neurofibromas).

MECHANISM

The NF1 gene is a tumour suppressor gene and NF1 individuals have an increased risk for a long list of tumours, all resulting from a second hit in the normal copy of the NF1 gene. Remarkably, some non-tumour phenotypes such as CALM and pseudarthrosis are also caused by a "second hit". Germline mutations inactivating the NF1 gene show a large variability in genetic mechanisms ranging from single-nucleotide substitutions and somatic mosaicism to large deletions affecting neighbouring genes. Molecular confirmation of the clinical diagnosis is becoming increasingly more important to differentiate NF1 from other syndromes such as Legius syndrome, to investigate genotype-phenotype correlations relevant in 10% of cases and to detect somatic mosaicism.

SURVEILLANCE AND THERAPY

Some degree of learning difficulties, attention deficit and social problems are observed in most children and affect quality of life. There is a large individual variability in complications and the evolution of the disease is difficult to predict. Specialised outpatient clinics for children have been widely established and are important for surveillance and guidance. Regular surveillance is also important for adolescents and adults because many tumour complications can be detected by whole-body MRI and treated even before symptoms develop and irreversible damage occurs. Recent data on nodular plexiform neurofibromas with continued growth in adolescents and young adults show that many of these tumours are premalignant lesions called atypical neurofibromatous neoplasm of uncertain biological potential (ANNUBP). Specific surveillance and timely local resection of these benign peripheral nerve sheath tumours might be important to prevent malignant degeneration. In the last years, targeted therapy with MEK inhibitors has shown promise to treat unresectable and symptomatic plexiform neurofibromas. Many more challenges remain to find the best way to monitor children and adults for potential complications and to find a satisfying cure for many complications in this disorder.

New Understanding of the Relevant Role of LINE-1 Retrotransposition in Human Disease and Immune Modulation

Long interspersed nuclear element-1 (LINE-1) retrotransposition is a major hallmark of cancer accompanied by global chromosomal instability, genomic instability, and genetic heterogeneity and has become one indicator for the occurrence, development, and poor prognosis of many diseases. LINE-1 also modulates the immune system and affects the immune microenvironment in a variety of ways. Aberrant expression of LINE-1 retrotransposon can provide strong stimuli for an innate immune response, activate the immune system, and induce autoimmunity and inflammation. Therefore, inhibition the activity of LINE-1 has become a potential treatment strategy for various diseases. In this review, we discussed the components and regulatory mechanisms involved with LINE-1, its correlations with disease and immunity, and multiple inhibitors of LINE-1, providing a new understanding of LINE-1.

Genetic breakdown of a Tet-off conditional lethality system for insect population control

Genetically modified conditional lethal strains have been created to improve the control of insect pest populations damaging to human health and agriculture. However, understanding the potential for the genetic breakdown of lethality systems by rare spontaneous mutations, or selection for inherent suppressors, is critical since field release studies are in progress. This knowledge gap was addressed in a Drosophila tetracycline-suppressible embryonic lethality system by analyzing the frequency and structure of primary-site spontaneous mutations and second-site suppressors resulting in heritable survivors from 1.2 million zygotes. Here we report that F₁ survivors due to primary-site deletions and indels occur at a 5.8 × 10⁻⁶ frequency, while survival due to second-site maternal-effect suppressors occur at a ~10⁻⁵ frequency. Survivors due to inherent lethal effector suppressors could result in a resistant field population, and we suggest that this risk may be mitigated by the use of dual redundant, albeit functionally unrelated, lethality systems.

Insect population control using conditional lethal systems could break down due to spontaneous mutations that render the system ineffective. Here the authors analyse the structure and frequency of such mutations in Drosophila and suggest the use of dual lethality systems to mitigate their survival.

Null phenotype of neurofibromatosis type 1 in a carrier of a heterozygous atypical NF1 deletion due to mosaicism

We coincidently detected an atypical deletion of at least 1.3-Mb, encompassing the NF1 tumor suppressor gene and several adjacent genes at an apparent heterozygous level in the blood of a 65-year-old female patient. She had multiple subcutaneous tumors that appeared with a certain similarity of subcutaneous neurofibromas, which, however, was revealed as lipomas by histological examination. Comprehensive and exhaustive clinical and radiological examinations did not detect any neurofibromatosis type 1-related clinical symptoms in the patient. Multiplex ligation-dependent probe amplification detected no or only very low level of the 1.3-Mb NF1 deletion in six lipomas and two skin biopsies. Digital polymerase chain reaction estimated the proportion of cells carrying a heterozygous NF1 deletion at 87% in the blood, and 8%, 10%, 13%, 17%, and 20%, respectively, in the five lipomas investigated by this method, confirming our hypothesis of mosaicism. Our findings suggest that de novo cases of genetic disease are potentially mosaic regardless of finding the mutation at an apparently heterozygous level in the blood and that the possibility of mosaicism should be considered in genotype-phenotype studies and genetic counseling.

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.

Can Designer Indels Be Tailored by Gene Editing?: Can Indels Be Customized?

Genome editing with engineered nucleases (GEENs) introduce site-specific DNA double-strand breaks (DSBs) and repairs DSBs via nonhomologous end-joining (NHEJ) pathways that eventually create indels (insertions/deletions) in a genome. Whether the features of indels resulting from gene editing could be customized is asked. A review of the literature reveals how gene editing technologies via NHEJ pathways impact gene editing. The survey consolidates a body of literature that suggests that the type (insertion, deletion, and complex) and the approximate length of indel edits can be somewhat customized with different GEENs and by manipulating the expression of key NHEJ genes. Structural data suggest that binding of GEENs to DNA may interfere with binding of key components of DNA repair complexes, favoring either classical- or alternative-NHEJ. The hypotheses have some limitations, but if validated, will enable scientists to better control indel makeup, holding promise for basic science and clinical applications of gene editing. Also see the video abstract here https://youtu.be/vTkJtUsLi3w.

Cancer mutational burden is shaped by G4 DNA, replication stress and mitochondrial dysfunction

A hallmark of cancer is genomic instability, which can enable cancer cells to evade therapeutic strategies. Here we employed a computational approach to uncover mechanisms underlying cancer mutational burden by focusing upon relationships between 1) translocation breakpoints and the thousands of G4 DNA-forming sequences within retrotransposons impacting transcription and exemplifying probable non-B DNA structures and 2) transcriptome profiling and cancer mutations. We determined the location and number of G4 DNA-forming sequences in the Genome Reference Consortium Human Build 38 and found a total of 358,605 covering ∼13.4 million bases. By analyzing >97,000 unique translocation breakpoints from the Catalogue Of Somatic Mutations In Cancer (COSMIC), we found that breakpoints are overrepresented at G4 DNA-forming sequences within hominid-specific SVA retrotransposons, and generally occur in tumors with mutations in tumor suppressor genes, such as TP53. Furthermore, correlation analyses between mRNA levels and exome mutational loads from The Cancer Genome Atlas (TCGA) encompassing >450,000 gene-mutation regressions revealed strong positive and negative associations, which depended upon tissue of origin. The strongest positive correlations originated from genes not listed as cancer genes in COSMIC; yet, these show strong predictive power for survival in most tumor types by Kaplan-Meier estimation. Thus, correlation analyses of DNA structure and gene expression with mutation loads complement and extend more traditional approaches to elucidate processes shaping genomic instability in cancer. The combined results point to G4 DNA, activation of cell cycle/DNA repair pathways, and mitochondrial dysfunction as three major factors driving the accumulation of somatic mutations in cancer cells.

Retrotransposon-induced mosaicism in the neural genome

Over the past decade, major discoveries in retrotransposon biology have depicted the neural genome as a dynamic structure during life. In particular, the retrotransposon LINE-1 (L1) has been shown to be transcribed and mobilized in the brain. Retrotransposition in the developing brain, as well as during adult neurogenesis, provides a milieu in which neural diversity can arise. Dysregulation of retrotransposon activity may also contribute to neurological disease. Here, we review recent reports of retrotransposon activity in the brain, and discuss the temporal nature of retrotransposition and its regulation in neural cells in response to stimuli. We also put forward hypotheses regarding the significance of retrotransposons for brain development and neurological function, and consider the potential implications of this phenomenon for neuropsychiatric and neurodegenerative conditions.

Ultra-deep amplicon sequencing indicates absence of low-grade mosaicism with normal cells in patients with type-1 NF1 deletions

Different types of large NF1 deletion are distinguishable by breakpoint location and potentially also by the frequency of mosaicism with normal cells lacking the deletion. However, low-grade mosaicism with fewer than 10% normal cells has not yet been excluded for all NF1 deletion types since it is impossible to assess by the standard techniques used to identify such deletions, including MLPA and array analysis. Here, we used ultra-deep amplicon sequencing to investigate the presence of normal cells in the blood of 20 patients with type-1 NF1 deletions lacking mosaicism according to MLPA. The ultra-deep sequencing entailed the screening of 96 amplicons for heterozygous SNVs located within the NF1 deletion region. DNA samples from three previously identified patients with type-2 NF1 deletions and low-grade mosaicism with normal cells as determined by FISH or microsatellite marker analysis were used to validate our methodology. In these type-2 NF1 deletion samples, proportions of 5.3%, 6.6% and 15.0% normal cells, respectively, were detected by ultra-deep amplicon sequencing. However, using this highly sensitive method, none of the 20 patients with type-1 NF1 deletions included in our analysis exhibited low-grade mosaicism with normal cells in blood, thereby supporting the view that the vast majority of type-1 deletions are germline deletions.

Retrotransposable genetic elements causing neutrophil defects

BACKGROUND

Retrotransposable elements are stretches of DNA that encode proteins with the inherent ability to insert their own RNA or another RNA by reverse transcriptase as DNA into a new genomic location. In humans, the only autonomous retrotransposable elements are members of the Long INterspersed Element-1 (LINE-1) family. LINE-1s may cause gene inactivation and human disease.

DESIGN

We present a brief summary of the published knowledge about LINE-1s in humans and the RNAs that these elements can transpose, and we focus on the effect of LINE-1-mediated retrotransposition on human neutrophil function.

RESULTS

Retrotransposons can cause genetic disease by two primary mechanisms: (1) insertional mutagenesis and (2) nonallelic homologous recombination. The only known neutrophil function affected by retrotransposition is that of NADPH oxidase activity. Four patients with chronic granulomatous disease (CGD) are known with LINE-1-mediated insertional inactivation of CYBB, the gene that encodes the gp91^phox component of the phagocyte NADPH oxidase. In addition, 5 CGD patients had a large deletion in the NCF2 gene, encoding the p67^phox component, and 2 CGD patients had a similar deletion in NCF1, encoding p47^phox . These deletions were caused by nonallelic homologous recombination between 2 Alu elements at the borders of each deletion. Alu elements have spread throughout the human genome by LINE-1 retrotransposition.

CONCLUSIONS

Probably, the occurrence of LINE-1-mediated insertions causing autosomal CGD has been underestimated. It might be worthwhile to reinvestigate the DNA from autosomal CGD patients with missplice mutations and large deletions for indications of LINE-1-mediated insertions.

A Homozygous Ancestral SVA-Insertion-Mediated Deletion in WDR66 Induces Multiple Morphological Abnormalities of the Sperm Flagellum and Male Infertility

Multiple morphological abnormalities of the sperm flagellum (MMAF) is a severe form of male infertility defined by the presence of a mosaic of anomalies, including short, bent, curled, thick, or absent flagella, resulting from a severe disorganization of the axoneme and of the peri-axonemal structures. Mutations in DNAH1, CFAP43, and CFAP44, three genes encoding axoneme-related proteins, have been described to account for approximately 30% of the MMAF cases reported so far. Here, we searched for pathological copy-number variants in whole-exome sequencing data from a cohort of 78 MMAF-affected subjects to identify additional genes associated with MMAF. In 7 of 78 affected individuals, we identified a homozygous deletion that removes the two penultimate exons of WDR66 (also named CFAP251), a gene coding for an axonemal protein preferentially localized in the testis and described to localize to the calmodulin- and spoke-associated complex at the base of radial spoke 3. Sequence analysis of the breakpoint region revealed in all deleted subjects the presence of a single chimeric SVA (SINE-VNTR-Alu) at the breakpoint site, suggesting that the initial deletion event was potentially mediated by an SVA insertion-recombination mechanism. Study of Trypanosoma WDR66's ortholog (TbWDR66) highlighted high sequence and structural analogy with the human protein and confirmed axonemal localization of the protein. Reproduction of the human deletion in TbWDR66 impaired flagellar movement, thus confirming WDR66 as a gene associated with the MMAF phenotype and highlighting the importance of the WDR66 C-terminal region.

Extreme clustering of type-1 NF1 deletion breakpoints co-locating with G-quadruplex forming sequences

The breakpoints of type-1 NF1 deletions encompassing 1.4-Mb are located within NF1-REPa and NF1-REPc, which exhibit a complex structure comprising different segmental duplications in direct and inverted orientation. Here, we systematically assessed the proportion of type-1 NF1 deletions caused by nonallelic homologous recombination (NAHR) and those mediated by other mutational mechanisms. To this end, we analyzed 236 unselected type-1 deletions and observed that 179 of them (75.8%) had breakpoints located within the NAHR hotspot PRS2, whereas 39 deletions (16.5%) had breakpoints located within PRS1. Sixteen deletions exhibited breakpoints located outside of these NAHR hotspots but were also mediated by NAHR. Taken together, the breakpoints of 234 (99.2%) of the 236 type-1 NF1 deletions were mediated by NAHR. Thus, NF1-REPa and NF1-REPc are strongly predisposed to recurrent NAHR, the main mechanism underlying type-1 NF1 deletions. We also observed a non-random overlap between type-1 NF1-deletion breakpoints and G-quadruplex forming sequences (GQs) as well as regions flanking PRDM9_A binding-sites. These findings imply that GQs and PRDM9_A binding-sites contribute to the clustering of type-1 deletion breakpoints. The co-location of both types of sequence was at its highest within PRS2, indicative of their synergistic contribution to the greatly increased NAHR activity within this hotspot.

Pronounced maternal parent-of-origin bias for type-1 NF1 microdeletions

Neurofibromatosis type 1 (NF1) is caused, in 4.7-11% of cases, by large deletions encompassing the NF1 gene and its flanking regions within 17q11.2. Different types of large NF1 deletion occur which are distinguishable by their breakpoint location and underlying mutational mechanism. Most common are the type-1 NF1 deletions of 1.4 Mb which exhibit recurrent breakpoints caused by nonallelic homologous recombination (NAHR), also termed unequal crossover. Here, we analyzed 37 unrelated families of patients with de novo type-1 NF1 deletions by means of short tandem repeat (STR) profiling to determine the parental origin of the deletions. We observed that 33 of the 37 type-1 deletions were of maternal origin (89.2% of cases; p < 0.0001). Analysis of the patients' siblings indicated that, in 14 informative cases, ten (71.4%) deletions resulted from interchromosomal unequal crossover during meiosis I. Our findings indicate a strong maternal parent-of-origin bias for type-1 NF1 deletions. A similarly pronounced maternal transmission bias has been reported for recurrent copy number variants (CNVs) within 16p11.2 associated with autism, but not so far for any other NAHR-mediated pathogenic CNVs. Region-specific genomic features are likely to be responsible for the maternal bias in the origin of both the 16p11.2 CNVs and type-1 NF1 deletions.

Transposable elements: genome innovation, chromosome diversity, and centromere conflict

Although it was nearly 70 years ago when transposable elements (TEs) were first discovered “jumping” from one genomic location to another, TEs are now recognized as contributors to genomic innovations as well as genome instability across a wide variety of species. In this review, we illustrate the ways in which active TEs, specifically retroelements, can create novel chromosome rearrangements and impact gene expression, leading to disease in some cases and species-specific diversity in others. We explore the ways in which eukaryotic genomes have evolved defense mechanisms to temper TE activity and the ways in which TEs continue to influence genome structure despite being rendered transpositionally inactive. Finally, we focus on the role of TEs in the establishment, maintenance, and stabilization of critical, yet rapidly evolving, chromosome features: eukaryotic centromeres. Across centromeres, specific types of TEs participate in genomic conflict, a balancing act wherein they are actively inserting into centromeric domains yet are harnessed for the recruitment of centromeric histones and potentially new centromere formation.

Disease onset in X-linked dystonia-parkinsonism correlates with expansion of a hexameric repeat within an SVA retrotransposon in TAF1

The genetic basis of X-Linked dystonia-parkinsonism (XDP) has been difficult to unravel, in part because all patients inherit the same haplotype of seven sequence variants, none of which has ever been identified in control individuals. This study revealed that one of the haplotype markers, a retrotransposon insertion within an intron of TAF1, has a variable number of hexameric repeats among affected individuals with an increase in repeat number strongly correlated with earlier age at disease onset. These data support a contributing role for this sequence in disease pathogenesis while further suggesting that XDP may be part of a growing list of neurodegenerative disorders associated with unstable repeat expansions.

X-linked dystonia-parkinsonism (XDP) is a neurodegenerative disease associated with an antisense insertion of a SINE-VNTR-Alu (SVA)-type retrotransposon within an intron of TAF1. This unique insertion coincides with six additional noncoding sequence changes in TAF1, the gene that encodes TATA-binding protein–associated factor-1, which appear to be inherited together as an identical haplotype in all reported cases. Here we examined the sequence of this SVA in XDP patients (n = 140) and detected polymorphic variation in the length of a hexanucleotide repeat domain, (CCCTCT)_n. The number of repeats in these cases ranged from 35 to 52 and showed a highly significant inverse correlation with age at disease onset. Because other SVAs exhibit intrinsic promoter activity that depends in part on the hexameric domain, we assayed the transcriptional regulatory effects of varying hexameric lengths found in the unique XDP SVA retrotransposon using luciferase reporter constructs. When inserted sense or antisense to the luciferase reading frame, the XDP variants repressed or enhanced transcription, respectively, to an extent that appeared to vary with length of the hexamer. Further in silico analysis of this SVA sequence revealed multiple motifs predicted to form G-quadruplexes, with the greatest potential detected for the hexameric repeat domain. These data directly link sequence variation within the XDP-specific SVA sequence to phenotypic variability in clinical disease manifestation and provide insight into potential mechanisms by which this intronic retroelement may induce transcriptional interference in TAF1 expression.

Analysis of LINE-1 Elements in DNA from Postmortem Brains of Individuals with Schizophrenia

Whereas some rare genetic variants convey high risk for schizophrenia (SZ), common alleles conveying even moderate risk remain elusive. Long interspersed element-1s (L1) are mobile retrotransposons comprising ~17% of the human genome. L1 retrotransposition can cause somatic mosaicism during neurodevelopment by insertional mutagenesis. We hypothesized that, compared to controls, patients diagnosed with schizophrenia (PDS) may have increased numbers of deleterious L1 insertions, perhaps occurring de novo, in brain-expressed genes of dorsolateral prefrontal cortex (DLPFC) neurons. Neuronal and non-neuronal nuclei were separated by fluorescence-activated cell sorting from postmortem DLPFC of 36 PDS and 26 age-matched controls. Genomic sequences flanking the 3'-side of L1s were amplified from neuronal DNA, and neuronal L1 libraries were sequenced. Aligned sequences were analyzed for L1 insertions using custom bioinformatics programs. Ontology and pathway analyses were done on lists of genes putatively disrupted by L1s in PDS and controls. Cellular or population allele frequencies of L1s were assessed by droplet digital PCR or Taqman genotyping. We observed a statistically significant increase in the proportion of intragenic novel L1s in DLPFC of PDS. We found over-representation of L1 insertions within the gene ontologies 'cell projection' and 'postsynaptic membrane' in the gene lists derived from PDS samples, but not from controls. Cellular allele frequencies of examined L1 insertions indicated heterozygosity in genomes of DLPFC cells. An L1 within ERI1 exoribonuclease family member 3 (ERI3) was found to associate with SZ. These results extend prior work documenting increased L1 genetic burden in the brains of PDS and also identify unique genes that may provide new insight into the pathophysiology of schizophrenia.

Copy number increases of transposable elements and protein-coding genes in an invasive fish of hybrid origin

Evolutionary dynamics of structural genetic variation in lineages of hybrid origin is not well explored, although structural mutations may increase in controlled hybrid crosses. We therefore tested whether structural variants accumulate in a fish of recent hybrid origin, invasive Cottus, relative to both parental species Cottus rhenanus and Cottus perifretum. Copy-number variation in exons of 10,979 genes was assessed using comparative genome hybridization arrays. Twelve genes showed significantly higher copy numbers in invasive Cottus compared to both parents. This coincided with increased expression for three genes related to vision, detoxification and muscle development, suggesting possible gene dosage effects. Copy number increases of putative transposons were assessed by comparative mapping of genomic DNA reads against a de novo assembly of 1,005 repetitive elements. In contrast to exons, copy number increases of repetitive elements were common (20.7%) in invasive Cottus, whereas decrease was very rare (0.01%). Among the increased repetitive elements, 53.8% occurred at higher numbers in C. perifretum compared to C. rhenanus, while only 1.4% were more abundant in C. rhenanus. This implies a biased mutational process that amplifies genetic material from one ancestor. To assess the frequency of de novo mutations through hybridization, we screened 64 laboratory-bred F2 offspring between the parental species for copy-number changes at five candidate loci. We found no evidence for new structural variants, indicating that they are too rare to be detected given our sampling scheme. Instead, they must have accumulated over more generations than we observed in a controlled cross.