Mechanisms underlying structural variant formation in genomic disorders

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.

Large TRAPPC11 gene deletions as a cause of muscular dystrophy and their estimated genesis

BACKGROUND

Transport protein particle (TRAPP) is a multiprotein complex that functions in localising proteins to the Golgi compartment. The TRAPPC11 subunit has been implicated in diseases affecting muscle, brain, eye and to some extent liver. We present three patients who are compound heterozygotes for a missense variant and a structural variant in the TRAPPC11 gene. TRAPPC11 structural variants have not yet been described in association with a disease. In order to reveal the estimated genesis of identified structural variants, we performed sequencing of individual breakpoint junctions and analysed the extent of homology and the presence of repetitive elements in and around the breakpoints.

METHODS

Biochemical methods including isoelectric focusing on serum transferrin and apolipoprotein C-III, as well as mitochondrial respiratory chain complex activity measurements, were used. Muscle biopsy samples underwent histochemical analysis. Next-generation sequencing was employed for identifying sequence variants associated with neuromuscular disorders, and Sanger sequencing was used to confirm findings.

RESULTS

We suppose that non-homologous end joining is a possible mechanism of deletion origin in two patients and non-allelic homologous recombination in one patient. Analyses of mitochondrial function performed in patients' skeletal muscles revealed an imbalance of mitochondrial metabolism, which worsens with age and disease progression.

CONCLUSION

Our results contribute to further knowledge in the field of neuromuscular diseases and mutational mechanisms. This knowledge is important for understanding the molecular nature of human diseases and allows us to improve strategies for identifying disease-causing mutations.

3q29 duplications: A cohort of 46 patients and a literature review

Duplications of the 3q29 cytoband are rare chromosomal copy number variations (CNVs) (overlapping or recurrent ~1.6 Mb 3q29 duplications). They have been associated with highly variable neurodevelopmental disorders (NDDs) with various associated features or reported as a susceptibility factor to the development of learning disabilities and neuropsychiatric disorders. The smallest region of overlap and the phenotype of 3q29 duplications remain uncertain. We here report a French cohort of 31 families with a 3q29 duplication identified by chromosomal microarray analysis (CMA), including 14 recurrent 1.6 Mb duplications, eight overlapping duplications (>1 Mb), and nine small duplications (<1 Mb). Additional genetic findings that may be involved in the phenotype were identified in 11 patients. Focusing on apparently isolated 3q29 duplications, patients present mainly mild NDD as suggested by a high rate of learning disabilities in contrast to a low proportion of patients with intellectual disabilities. Although some are de novo, most of the 3q29 duplications are inherited from a parent with a similar mild phenotype. Besides, the study of small 3q29 duplications does not provide evidence for any critical region. Our data suggest that the overlapping and recurrent 3q29 duplications seem to lead to mild NDD and that a severe or syndromic clinical presentation should warrant further genetic analyses.

Breaking new ground: Exploring de novo chromosomal rearrangements in 1p36 microdeletion

Chromosomal structural variations (SVs) are linked to a wide range of phenotypes and arise due to disruptions during DNA replication, which can affect gene function within the SV regions. This case report details a patient diagnosed with neurodevelopmental delay. Detailed investigation through array comparative genomic hybridization revealed two pathogenic SVs on chromosome 1, which align with a 1p36 microdeletion, and a microduplication at 2p35.3, the latter being classified as a variant of unknown significance. The patient's clinical presentation is consistent with the 1p36 deletion syndrome, characterized by specific developmental delays and physical anomalies. Further genetic analysis suggests that these terminal rearrangements might stem from an unbalanced translocation between the short arms of chromosomes 1 and 2. This case underscores the complexity of interpreting multiple concurrent SVs and their cumulative effect on phenotype. Ongoing research into such chromosomal abnormalities will enhance our understanding of their clinical manifestations and guide more targeted therapeutic strategies.

Low-pass whole genome sequencing as a cost-effective alternative to chromosomal microarray analysis for low- and middle-income countries

Low-pass whole genome sequencing (LP-WGS) has been applied as alternative method to detect copy number variants (CNVs) in the clinical setting. Compared with chromosomal microarray analysis (CMA), the sequencing-based approach provides a similar resolution of CNV detection at a lower cost. In this study, we assessed the efficiency and reliability of LP-WGS as a more affordable alternative to CMA. A total of 1363 patients with unexplained neurodevelopmental delay/intellectual disability, autism spectrum disorders, and/or multiple congenital anomalies were enrolled. Those patients were referred from 15 nonprofit organizations and university centers located in different states in Brazil. The analysis of LP-WGS at 1x coverage (>50kb) revealed a positive testing result in 22% of the cases (304/1363), in which 219 and 85 correspond to pathogenic/likely pathogenic (P/LP) CNVs and variants of uncertain significance (VUS), respectively. The 16% (219/1363) diagnostic yield observed in our cohort is comparable to the 15%-20% reported for CMA in the literature. The use of commercial software, as demonstrated in this study, simplifies the implementation of the test in clinical settings. Particularly for countries like Brazil, where the cost of CMA presents a substantial barrier to most of the population, LP-WGS emerges as a cost-effective alternative for investigating copy number changes in cytogenetics.

Somatic and intergenerational G4C2 hexanucleotide repeat instability in a human C9orf72 knock-in mouse model

Expansion of a G4C2 repeat in the C9orf72 gene is associated with familial Amyotrophic Lateral Sclerosis (ALS) and Frontotemporal Dementia (FTD). To investigate the underlying mechanisms of repeat instability, which occurs both somatically and intergenerationally, we created a novel mouse model of familial ALS/FTD that harbors 96 copies of G4C2 repeats at a humanized C9orf72 locus. In mouse embryonic stem cells, we observed two modes of repeat expansion. First, we noted minor increases in repeat length per expansion event, which was dependent on a mismatch repair pathway protein Msh2. Second, we found major increases in repeat length per event when a DNA double- or single-strand break (DSB/SSB) was artificially introduced proximal to the repeats, and which was dependent on the homology-directed repair (HDR) pathway. In mice, the first mode primarily drove somatic repeat expansion. Major changes in repeat length, including expansion, were observed when SSB was introduced in one-cell embryos, or intergenerationally without DSB/SSB introduction if G4C2 repeats exceeded 400 copies, although spontaneous HDR-mediated expansion has yet to be identified. These findings provide a novel strategy to model repeat expansion in a non-human genome and offer insights into the mechanism behind C9orf72 G4C2 repeat instability.

Diversity of ribosomes at the level of rRNA variation associated with human health and disease

Ribosomal DNA and RNA (rDNA and rRNA) sequences are usually discarded from sequencing analyses. But with hundreds of copies of rDNA genes it is unknown whether they possess sequence variations that form different types of ribosomes that affect human physiology and disease. Here, we developed an algorithm for variant-calling between paralog genes (termed RGA) and compared rDNA variations found in short- and long-read sequencing data from the 1,000 Genomes Project (1KGP) and Genome In A Bottle (GIAB). We additionally developed a novel protocol for long-read sequencing full-length rRNA (RIBO-RT) from actively translating ribosomes. Our analyses identified hundreds of rDNA variants, most of which, surprisingly, are short insertion-deletions (indels) and dozens of highly abundant rRNA variants that are incorporated into translationally active ribosomes. To visualize variant ribosomes at the single cell level, we developed an in-situ rRNA sequencing method (SWITCH-seq) which revealed that variants are co-expressed within individual cells. Strikingly, by analyzing rDNA, we found that variants assemble into distinct ribosome subtypes. We discovered that these subtypes acquire different rRNA structures by successfully employing dimethyl sulfate (DMS) probing of full length rRNA. With this atlas we investigated rRNA variation changes across human tissues and cancer types. This revealed tissue-specific rRNA subtype expression in endoderm/ectoderm-derived tissues. In cancer, low abundant rRNA variants can become highly expressed, which suggests the presence of cancer-specific ribosomes. Together, this study identifies and comprehensively characterizes the diversity of ribosomes at the level of rRNA variants which is dominated by indel variants, their chromosomal location and unique structure as well as the association of ribosome variation with tissue-specific biology and cancer.

Whole F8 gene sequencing identified pathogenic structural variants in the remaining unsolved patients with severe hemophilia A

BACKGROUND

No F8 genetic abnormality is detected in approximately 1% to 2% of patients with severe hemophilia A (HA) using conventional genetic approaches. In these patients, deep intronic variation or F8 disrupting genomic rearrangement could be causal.

OBJECTIVES

The study aimed to identify the causal variation in families with a history of severe HA for whom genetic investigations failed.

METHODS

We performed whole F8 gene sequencing in 8 propositi. Genomic rearrangements were confirmed by Sanger sequencing of breakpoint junctions and/or quantitative polymerase chain reaction.

RESULTS

A structural variant disrupting F8 was found in each propositus, so that all the 815 families with a history of severe HA registered in our laboratory received a conclusive genetic diagnosis. These structural variants consisted of 3 balanced inversions, 3 large insertions of gained regions, and 1 retrotransposition of a mobile element. The 3 inversions were 105 Mb, 1.97 Mb, and 0.362 Mb in size. Among the insertions of gained regions, one corresponded to the insertion of a 34 kb gained region from chromosome 6q27 in F8 intron 6, another was the insertion of a 447 kb duplicated region from chromosome 9p22.1 in F8 intron 14, and the last one was the insertion of an Xq28 349 kb gained in F8 intron 5.

CONCLUSION

All the genetically unsolved cases of severe HA in this cohort were due to structural variants disrupting F8. This study highlights the effectiveness of whole F8 sequencing to improve the molecular diagnosis of HA when the conventional approach fails.

Deciphering male influence in gynogenetic Pengze crucian carp (Carassius auratus var. pengsenensis): insights from Nanopore sequencing of structural variations

In this study, we investigate gynogenetic reproduction in Pengze Crucian Carp (Carassius auratus var. pengsenensis) using third-generation Nanopore sequencing to uncover structural variations (SVs) in offspring. Our objective was to understand the role of male genetic material in gynogenesis by examining the genomes of both parents and their offspring. We discovered a notable number of male-specific structural variations (MSSVs): 1,195 to 1,709 MSSVs in homologous offspring, accounting for approximately 0.52%-0.60% of their detected SVs, and 236 to 350 MSSVs in heterologous offspring, making up about 0.10%-0.13%. These results highlight the significant influence of male genetic material on the genetic composition of offspring, particularly in homologous pairs, challenging the traditional view of asexual reproduction. The gene annotation of MSSVs revealed their presence in critical gene regions, indicating potential functional impacts. Specifically, we found 5 MSSVs in the exonic regions of protein-coding genes in homologous offspring, suggesting possible direct effects on protein structure and function. Validation of an MSSV in the exonic region of the polyunsaturated fatty acid 5-lipoxygenase gene confirmed male genetic material transmission in some offspring. This study underscores the importance of further research on the genetic diversity and gynogenesis mechanisms, providing valuable insights for reproductive biology, aquaculture, and fostering innovation in biological research and aquaculture practices.

Elevated cholesterol in ATAD3 mutants is a compensatory mechanism that leads to membrane cholesterol aggregation

Aberrant cholesterol metabolism causes neurological disease and neurodegeneration, and mitochondria have been linked to perturbed cholesterol homeostasis via the study of pathological mutations in the ATAD3 gene cluster. However, whether the cholesterol changes were compensatory or contributory to the disorder was unclear, and the effects on cell membranes and the wider cell were also unknown. Using patient-derived cells, we show that cholesterol perturbation is a conserved feature of pathological ATAD3 variants that is accompanied by an expanded lysosome population containing membrane whorls characteristic of lysosomal storage diseases. Lysosomes are also more numerous in Drosophila neural progenitor cells expressing mutant Atad3, which exhibit abundant membrane-bound cholesterol aggregates, many of which co-localize with lysosomes. By subjecting the Drosophila Atad3 mutant to nutrient restriction and cholesterol supplementation, we show that the mutant displays heightened cholesterol dependence. Collectively, these findings suggest that elevated cholesterol enhances tolerance to pathological ATAD3 variants; however, this comes at the cost of inducing cholesterol aggregation in membranes, which lysosomal clearance only partly mitigates.

A comparative study of structural variant calling in WGS from Alzheimer's disease families

Detecting structural variants (SVs) in whole-genome sequencing poses significant challenges. We present a protocol for variant calling, merging, genotyping, sensitivity analysis, and laboratory validation for generating a high-quality SV call set in whole-genome sequencing from the Alzheimer's Disease Sequencing Project comprising 578 individuals from 111 families. Employing two complementary pipelines, Scalpel and Parliament, for SV/indel calling, we assessed sensitivity through sample replicates (N = 9) with in silico variant spike-ins. We developed a novel metric, D-score, to evaluate caller specificity for deletions. The accuracy of deletions was evaluated by Sanger sequencing. We generated a high-quality call set of 152,301 deletions of diverse sizes. Sanger sequencing validated 114 of 146 detected deletions (78.1%). Scalpel excelled in accuracy for deletions ≤100 bp, whereas Parliament was optimal for deletions >900 bp. Overall, 83.0% and 72.5% of calls by Scalpel and Parliament were validated, respectively, including all 11 deletions called by both Parliament and Scalpel between 101 and 900 bp. Our flexible protocol successfully generated a high-quality deletion call set and a truth set of Sanger sequencing-validated deletions with precise breakpoints spanning 1-17,000 bp.

The Diagnostic Yield of Chromosomal Microarray Analysis in Third-Trimester Fetal Abnormalities

OBJECTIVE

 This study aimed to determine the diagnostic yield of chromosomal microarray analysis (CMA) performed in cases of fetal abnormalities detected during the third trimester of pregnancy.

STUDY DESIGN

 A retrospective review of medical records was conducted for women who underwent amniocentesis at or beyond 28 weeks of gestation between January 2017 and February 2023. CMA results of pregnancies with abnormal sonographic findings not detected before 28 weeks were included.

RESULTS

 A total of 482 fetuses met the inclusion criteria. The average maternal age was 31.3 years, and the average gestational age at amniocentesis was 32.3 weeks. The overall diagnostic yield of CMA was 6.2% (30 clinically significant copy number variations [CNVs]). The yield was 16.4% in cases with two or more fetal malformations, while cases with a single anomaly revealed a diagnostic yield of 7.3%. Cases presenting isolated polyhydramnios or isolated fetal growth restriction had a lower yield of 9.3 and 5.4%, respectively. Of the 30 clinically significant cases, 19 (or 63.4%) exhibited recurrent CNVs. The remaining 11 cases (or 36.6%) presented unique CNVs. The theoretical yield of Noninvasive Prenatal Testing (NIPT) in our cohort is 2% for aneuploidy, which implies that it could potentially miss up to 70% of the significant findings that could be identified by CMA. In 80% of the fetuses (or 24 out of 30) with clinically significant CNVs, the structural abnormalities detected on fetal ultrasound examinations corresponded with the CMA results.

CONCLUSION

 The 6.2% detection rate of significant CNVs in late-onset fetal anomalies confirms the value of CMA in third-trimester amniocentesis. The findings underscore the necessity of CMA for detecting CNVs potentially overlooked by NIPT and emphasize the importance of thorough genetic counseling.

KEY POINTS

· CMA yields 6.2% for third-trimester anomalies.. · NIPT may miss 70% of CMA findings.. · Ultrasound matched 80% of CMA results..

Mind the gap: the relevance of the genome reference to resolve rare and pathogenic inversions

Both long-read genome sequencing (lrGS) and the recently published Telomere to Telomere (T2T) reference genome provide increased coverage and resolution across repetitive regions promising heightened structural variant detection and improved mapping. Inversions (INV), intrachromosomal segments which are rotated 180° and inserted back into the same chromosome, are a class of structural variants particularly challenging to detect due to their copy-number neutral state and association with repetitive regions. Inversions represent about 1/20 of all balanced structural chromosome aberrations and can lead to disease by gene disruption or altering regulatory regions of dosage sensitive genes in cis . Here we remapped the genome data from six individuals carrying unsolved cytogenetically detected inversions. An INV6 and INV10 were resolved using GRCh38 and T2T-CHM13. Finally, an INV9 required optical genome mapping, de novo assembly of lrGS data and T2T-CHM13. This inversion disrupted intron 25 of EHMT1, confirming a diagnosis of Kleefstra syndrome 1 (MIM#610253). These three inversions, only mappable in specific references, prompted us to investigate the presence and population frequencies of differential reference regions (DRRs) between T2T-CHM13, GRCh37, GRCh38, the chimpanzee and bonobo, and hundreds of megabases of DRRs were identified. Our results emphasize the significance of the chosen reference genome and the added benefits of lrGS and optical genome mapping in solving rearrangements in challenging regions of the genome. This is particularly important for inversions and may impact clinical diagnostics.

Olfactory receptor genes and chromosome 11 structural aberrations: Players or spectators?

The largest multi-gene family in metazoans is the family of olfactory receptor (OR) genes. Human ORs are organized in clusters over most chromosomes and seem to include >0.1% the human genome. Because 369 out of 856 OR genes are mapped on chromosome 11 (HSA11), we sought to determine whether they mediate structural rearrangements involving this chromosome. To this aim, we analyzed 220 specimens collected during diagnostic procedures involving structural rearrangements of chromosome 11. A total of 222 chromosomal abnormalities were included, consisting of inversions, deletions, translocations, duplications, and one insertion, detected by conventional chromosome analysis and/or fluorescence in situ hybridization (FISH) and array comparative genomic hybridization (array-CGH). We verified by bioinformatics and statistical approaches the occurrence of breakpoints in cytobands with or without OR genes. We found that OR genes are not involved in chromosome 11 reciprocal translocations, suggesting that different DNA motifs and mechanisms based on homology or non-homology recombination can cause chromosome 11 structural alterations. We also considered the proximity between the chromosomal territories of chromosome 11 and its partner chromosomes involved in the translocations by using the deposited Hi-C data concerning the possible occurrence of chromosome interactions. Interestingly, most of the breakpoints are located in regions highly involved in chromosome interactions. Further studies should be carried out to confirm the potential role of chromosome territories' proximity in promoting genome structural variation, so fundamental in our understanding of the molecular basis of medical genetics and evolutionary genetics.

Severus: accurate detection and characterization of somatic structural variation in tumor genomes using long reads

Most current studies rely on short-read sequencing to detect somatic structural variation (SV) in cancer genomes. Long-read sequencing offers the advantage of better mappability and long-range phasing, which results in substantial improvements in germline SV detection. However, current long-read SV detection methods do not generalize well to the analysis of somatic SVs in tumor genomes with complex rearrangements, heterogeneity, and aneuploidy. Here, we present Severus: a method for the accurate detection of different types of somatic SVs using a phased breakpoint graph approach. To benchmark various short- and long-read SV detection methods, we sequenced five tumor/normal cell line pairs with Illumina, Nanopore, and PacBio sequencing platforms; on this benchmark Severus showed the highest F1 scores (harmonic mean of the precision and recall) as compared to long-read and short-read methods. We then applied Severus to three clinical cases of pediatric cancer, demonstrating concordance with known genetic findings as well as revealing clinically relevant cryptic rearrangements missed by standard genomic panels.

A novel large intragenic DPYD deletion causing dihydropyrimidine dehydrogenase deficiency: a case report

BACKGROUND

Dihydropyrimidine dehydrogenase (DPD), is the initial and rate-limiting enzyme in the catabolic pathway of pyrimidines. Deleterious variants in the DPYD gene cause DPD deficiency, a rare autosomal recessive disorder. The clinical spectrum of affected individuals is wide ranging from asymptomatic to severely affected patients presenting with intellectual disability, motor retardation, developmental delay and seizures. DPD is also important as the main enzyme in the catabolism of 5-fluorouracil (5-FU) which is extensively used as a chemotherapeutic agent. Even in the absence of clinical symptoms, individuals with either complete or partial DPD deficiency face a high risk of severe and even fatal fluoropyrimidine-associated toxicity. The identification of causative genetic variants in DPYD is therefore gaining increasing attention due to their potential use as predictive markers of fluoropyrimidine toxicity.

METHODS

A male infant patient displaying biochemical features of DPD deficiency was investigated by clinical exome sequencing. Bioinformatics tools were used for data analysis and results were confirmed by MLPA and Sanger sequencing.

RESULTS

A novel intragenic deletion of 71.2 kb in the DPYD gene was identified in homozygosity. The deletion, DPYD(NM_000110.4):c.850 + 23455_1128 + 8811del, eliminates exons 9 and 10 and may have resulted from a non-homologous end-joining event, as suggested by in silico analysis.

CONCLUSIONS

The study expands the spectrum of DPYD variants associated with DPD deficiency. Furthermore, it raises the concern that patients at risk for fluoropyrimidine toxicity due to DPYD deletions could be missed during pre-treatment genetic testing for the currently recommended single nucleotide polymorphisms.

The impact of FASTQ and alignment read order on structural variant calling from long-read sequencing data

Background

Structural variant (SV) calling from DNA sequencing data has been challenging due to several factors, including the ambiguity of short-read alignments, multiple complex SVs in the same genomic region, and the lack of "truth" datasets for benchmarking. Additionally, caller choice, parameter settings, and alignment method are known to affect SV calling. However, the impact of FASTQ read order on SV calling has not been explored for long-read data.

Results

Here, we used PacBio DNA sequencing data from 15 Caenorhabditis elegans strains and four Arabidopsis thaliana ecotypes to evaluate the sensitivity of different SV callers on FASTQ read order. Comparisons of variant call format files generated from the original and permutated FASTQ files demonstrated that the order of input data affected the SVs predicted by each caller. In particular, pbsv was highly sensitive to the order of the input data, especially at the highest depths where over 70% of the SV calls generated from pairs of differently ordered FASTQ files were in disagreement. These demonstrate that read order sensitivity is a complex, multifactorial process, as the differences observed both within and between species varied considerably according to the specific combination of aligner, SV caller, and sequencing depth. In addition to the SV callers being sensitive to the input data order, the SAMtools alignment sorting algorithm was identified as a source of variability following read order randomization.

Conclusion

The results of this study highlight the sensitivity of SV calling on the order of reads encoded in FASTQ files, which has not been recognized in long-read approaches. These findings have implications for the replication of SV studies and the development of consistent SV calling protocols. Our study suggests that researchers should pay attention to the input order sensitivity of read alignment sorting methods when analyzing long-read sequencing data for SV calling, as mitigating a source of variability could facilitate future replication work. These results also raise important questions surrounding the relationship between SV caller read order sensitivity and tool performance. Therefore, tool developers should also consider input order sensitivity as a potential source of variability during the development and benchmarking of new and improved methods for SV calling.

Dynamic interplay between human alpha-satellite DNA structure and centromere functions

Maintenance of genome stability relies on functional centromeres for correct chromosome segregation and faithful inheritance of the genetic information. The human centromere is the primary constriction within mitotic chromosomes made up of repetitive alpha-satellite DNA hierarchically organized in megabase-long arrays of near-identical higher order repeats (HORs). Centromeres are epigenetically specified by the presence of the centromere-specific histone H3 variant, CENP-A, which enables the assembly of the kinetochore for microtubule attachment. Notably, centromeric DNA is faithfully inherited as intact haplotypes from the parents to the offspring without intervening recombination, yet, outside of meiosis, centromeres are akin to common fragile sites (CFSs), manifesting crossing-overs and ongoing sequence instability. Consequences of DNA changes within the centromere are just starting to emerge, with unclear effects on intra- and inter-generational inheritance driven by centromere's essential role in kinetochore assembly. Here, we review evidence of meiotic selection operating to mitigate centromere drive, as well as recent reports on centromere damage, recombination and repair during the mitotic cell division. We propose an antagonistic pleiotropy interpretation to reconcile centromere DNA instability as both driver of aneuploidy that underlies degenerative diseases, while also potentially necessary for the maintenance of homogenized HORs for centromere function. We attempt to provide a framework for this conceptual leap taking into consideration the structural interface of centromere-kinetochore interaction and present case scenarios for its malfunctioning. Finally, we offer an integrated working model to connect DNA instability, chromatin, and structural changes with functional consequences on chromosome integrity.

Benchmarking long-read aligners and SV callers for structural variation detection in Oxford nanopore sequencing data

Structural variants (SVs) are one of the significant types of DNA mutations and are typically defined as larger-than-50-bp genomic alterations that include insertions, deletions, duplications, inversions, and translocations. These modifications can profoundly impact the phenotypic characteristics and contribute to disorders like cancer, response to treatment, and infections. Four long-read aligners and five SV callers have been evaluated using three Oxford Nanopore NGS human genome datasets in terms of precision, recall, and F1-score statistical metrics, depth of coverage, and speed of analysis. The best SV caller regarding recall, precision, and F1-score when matched with different aligners at different coverage levels tend to vary depending on the dataset and the specific SV types being analyzed. However, based on our findings, Sniffles and CuteSV tend to perform well across different aligners and coverage levels, followed by SVIM, PBSV, and SVDSS in the last place. The CuteSV caller has the highest average F1-score (82.51%) and recall (78.50%), and Sniffles has the highest average precision value (94.33%). Minimap2 as an aligner and Sniffles as an SV caller act as a strong base for the pipeline of SV calling because of their high speed and reasonable accomplishment. PBSV has a lower average F1-score, precision, and recall and may generate more false positives and overlook some actual SVs. Our results are valuable in the comprehensive evaluation of popular SV callers and aligners as they provide insight into the performance of several long-read aligners and SV callers and serve as a reference for researchers in selecting the most suitable tools for SV detection.

RPA and Rad27 limit templated and inverted insertions at DNA breaks

Formation of templated insertions at DNA double-strand breaks (DSBs) is very common in cancer cells. The mechanisms and enzymes regulating these events are largely unknown. Here, we investigated templated insertions in yeast at DSBs using amplicon sequencing across a repaired locus. We document very short (most ∼5-34 bp), templated inverted duplications at DSBs. They are generated through a foldback mechanism that utilizes microhomologies adjacent to the DSB. Enzymatic requirements suggest a hybrid mechanism wherein one end requires Polδ-mediated synthesis while the other end is captured by nonhomologous end joining (NHEJ). This process is exacerbated in mutants with low levels or mutated RPA ( rtt105 Δ; rfa1 -t33) or extensive resection mutant ( sgs1 Δ exo1 Δ). Templated insertions from various distant genomic locations also increase in these mutants as well as in rad27 Δ and originate from fragile regions of the genome. Among complex insertions, common events are insertions of two sequences, originating from the same locus and with inverted orientation. We propose that these inversions are also formed by microhomology-mediated template switching. Taken together, we propose that a shortage of RPA typical in cancer cells is one possible factor stimulating the formation of templated insertions.

Long read sequencing on its way to the routine diagnostics of genetic diseases

The clinical application of technological progress in the identification of DNA alterations has always led to improvements of diagnostic yields in genetic medicine. At chromosome side, from cytogenetic techniques evaluating number and gross structural defects to genomic microarrays detecting cryptic copy number variants, and at molecular level, from Sanger method studying the nucleotide sequence of single genes to the high-throughput next-generation sequencing (NGS) technologies, resolution and sensitivity progressively increased expanding considerably the range of detectable DNA anomalies and alongside of Mendelian disorders with known genetic causes. However, particular genomic regions (i.e., repetitive and GC-rich sequences) are inefficiently analyzed by standard genetic tests, still relying on laborious, time-consuming and low-sensitive approaches (i.e., southern-blot for repeat expansion or long-PCR for genes with highly homologous pseudogenes), accounting for at least part of the patients with undiagnosed genetic disorders. Third generation sequencing, generating long reads with improved mappability, is more suitable for the detection of structural alterations and defects in hardly accessible genomic regions. Although recently implemented and not yet clinically available, long read sequencing (LRS) technologies have already shown their potential in genetic medicine research that might greatly impact on diagnostic yield and reporting times, through their translation to clinical settings. The main investigated LRS application concerns the identification of structural variants and repeat expansions, probably because techniques for their detection have not evolved as rapidly as those dedicated to single nucleotide variants (SNV) identification: gold standard analyses are karyotyping and microarrays for balanced and unbalanced chromosome rearrangements, respectively, and southern blot and repeat-primed PCR for the amplification and sizing of expanded alleles, impaired by limited resolution and sensitivity that have not been significantly improved by the advent of NGS. Nevertheless, more recently, with the increased accuracy provided by the latest product releases, LRS has been tested also for SNV detection, especially in genes with highly homologous pseudogenes and for haplotype reconstruction to assess the parental origin of alleles with de novo pathogenic variants. We provide a review of relevant recent scientific papers exploring LRS potential in the diagnosis of genetic diseases and its potential future applications in routine genetic testing.

Comparative genomic analysis of chemosensory-related gene families in gastropods

Chemoreception is critical for the survival and reproduction of animals. Except for a reduced group of insects and chelicerates, the molecular identity of chemosensory proteins is poorly understood in invertebrates. Gastropoda is the extant mollusk class with the greatest species richness, including marine, freshwater, and terrestrial lineages, and likely, highly diverse chemoreception systems. Here, we performed a comprehensive comparative genome analysis taking advantage of the chromosome-level information of two Gastropoda species, one of which belongs to a lineage that underwent a whole genome duplication event. We identified thousands of previously uncharacterized chemosensory-related genes, the majority of them encoding G protein-coupled receptors (GPCR), mostly organized into clusters distributed across all chromosomes. We also detected gene families encoding degenerin epithelial sodium channels (DEG-ENaC), ionotropic receptors (IR), sensory neuron membrane proteins (SNMP), Niemann-Pick type C2 (NPC2) proteins, and lipocalins, although with a lower number of members. Our phylogenetic analysis of the GPCR gene family across protostomes revealed: (i) remarkable gene family expansions in Gastropoda; (ii) clades including members from all protostomes; and (iii) species-specific clades with a substantial number of receptors. For the first time, we provide new and valuable knowledge into the evolution of the chemosensory gene families in invertebrates other than arthropods.

Low-pass whole genome sequencing is a reliable and cost-effective approach for copy number variant analysis in the clinical setting

INTRODUCTION

Next generation sequencing technology has greatly reduced the cost and time required for sequencing a genome. An approach that is rapidly being adopted as an alternative method for CNV analysis is the low-pass whole genome sequencing (LP-WGS). Here, we evaluated the performance of LP-WGS to detect copy number variants (CNVs) in clinical cytogenetics.

MATERIALS AND METHODS

DNA samples with known CNVs detected by chromosomal microarray analyses (CMA) were selected for comparison and used as positive controls; our panel included 44 DNA samples (12 prenatal and 32 postnatal), comprising a total of 55 chromosome imbalances. The selected cases were chosen to provide a wide range of clinically relevant CNVs, the vast majority being associated with intellectual disability or recognizable syndromes. The chromosome imbalances ranged in size from 75 kb to 90.3 Mb, including aneuploidies and two cases of mosaicism.

RESULTS

All CNVs were successfully detected by LP-WGS, showing a high level of consistency and robust performance of the sequencing method. Notably, the size of chromosome imbalances detected by CMA and LP-WGS were compatible between the two different platforms, which indicates that the resolution and sensitivity of the LP-WGS approach are at least similar to those provided by CMA.

DISCUSSION

Our data show the potential use of LP-WGS to detect CNVs in clinical diagnosis and confirm the method as an alternative for chromosome imbalances detection. The diagnostic effectiveness and feasibility of LP-WGS, in this technical validation study, were evidenced by a clinically representative dataset of CNVs that allowed a systematic assessment of the detection power and the accuracy of the sequencing approach. Further, since the software used in this study is commercially available, the method can easily be tested and implemented in a routine diagnostic setting.

Transposable Element Insertions Are Associated with Batesian Mimicry in the Pantropical Butterfly Hypolimnas misippus

Hypolimnas misippus is a Batesian mimic of the toxic African Queen butterfly (Danaus chrysippus). Female H. misippus butterflies use two major wing patterning loci (M and A) to imitate three color morphs of D. chrysippus found in different regions of Africa. In this study, we examine the evolution of the M locus and identify it as an example of adaptive atavism. This phenomenon involves a morphological reversion to an ancestral character that results in an adaptive phenotype. We show that H. misippus has re-evolved an ancestral wing pattern present in other Hypolimnas species, repurposing it for Batesian mimicry of a D. chrysippus morph. Using haplotagging, a linked-read sequencing technology, and our new analytical tool, Wrath, we discover two large transposable element insertions located at the M locus and establish that these insertions are present in the dominant allele responsible for producing mimetic phenotype. By conducting a comparative analysis involving additional Hypolimnas species, we demonstrate that the dominant allele is derived. This suggests that, in the derived allele, the transposable elements disrupt a cis-regulatory element, leading to the reversion to an ancestral phenotype that is then utilized for Batesian mimicry of a distinct model, a different morph of D. chrysippus. Our findings present a compelling instance of convergent evolution and adaptive atavism, in which the same pattern element has independently evolved multiple times in Hypolimnas butterflies, repeatedly playing a role in Batesian mimicry of diverse model species.

DAXX promotes centromeric stability independently of ATRX by preventing the accumulation of R-loop-induced DNA double-stranded breaks

Maintaining chromatin integrity at the repetitive non-coding DNA sequences underlying centromeres is crucial to prevent replicative stress, DNA breaks and genomic instability. The concerted action of transcriptional repressors, chromatin remodelling complexes and epigenetic factors controls transcription and chromatin structure in these regions. The histone chaperone complex ATRX/DAXX is involved in the establishment and maintenance of centromeric chromatin through the deposition of the histone variant H3.3. ATRX and DAXX have also evolved mutually-independent functions in transcription and chromatin dynamics. Here, using paediatric glioma and pancreatic neuroendocrine tumor cell lines, we identify a novel ATRX-independent function for DAXX in promoting genome stability by preventing transcription-associated R-loop accumulation and DNA double-strand break formation at centromeres. This function of DAXX required its interaction with histone H3.3 but was independent of H3.3 deposition and did not reflect a role in the repression of centromeric transcription. DAXX depletion mobilized BRCA1 at centromeres, in line with BRCA1 role in counteracting centromeric R-loop accumulation. Our results provide novel insights into the mechanisms protecting the human genome from chromosomal instability, as well as potential perspectives in the treatment of cancers with DAXX alterations.

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.

Novel syndromic neurodevelopmental disorder caused by de novo deletion of CHASERR, a long noncoding RNA

Genes encoding long non-coding RNAs (lncRNAs) comprise a large fraction of the human genome, yet haploinsufficiency of a lncRNA has not been shown to cause a Mendelian disease. CHASERR is a highly conserved human lncRNA adjacent to CHD2-a coding gene in which de novo loss-of-function variants cause developmental and epileptic encephalopathy. Here we report three unrelated individuals each harboring an ultra-rare heterozygous de novo deletion in the CHASERR locus. We report similarities in severe developmental delay, facial dysmorphisms, and cerebral dysmyelination in these individuals, distinguishing them from the phenotypic spectrum of CHD2 haploinsufficiency. We demonstrate reduced CHASERR mRNA expression and corresponding increased CHD2 mRNA and protein in whole blood and patient-derived cell lines-specifically increased expression of the CHD2 allele in cis with the CHASERR deletion, as predicted from a prior mouse model of Chaserr haploinsufficiency. We show for the first time that de novo structural variants facilitated by Alu-mediated non-allelic homologous recombination led to deletion of a non-coding element (the lncRNA CHASERR) to cause a rare syndromic neurodevelopmental disorder. We also demonstrate that CHD2 has bidirectional dosage sensitivity in human disease. This work highlights the need to carefully evaluate other lncRNAs, particularly those upstream of genes associated with Mendelian disorders.

Pangenomics: A new era in the field of neurodegenerative diseases

A pangenome is composed of all the genetic variability of a group of individuals, and its application to the study of neurodegenerative diseases may provide valuable insights into the underlying aspects of genetic heterogenetiy for these complex ailments, including gene expression, epigenetics, and translation mechanisms. Furthermore, a reference pangenome allows for the identification of previously undetected structural commonalities and differences among individuals, which may help in the diagnosis of a disease, support the prediction of what will happen over time (prognosis) and aid in developing novel treatments in the perspective of personalized medicine. Therefore, in the present review, the application of the pangenome concept to the study of neurodegenerative diseases will be discussed and analyzed for its potential to enable an improvement in diagnosis and prognosis for these illnesses, leading to the development of tailored treatments for individual patients from the knowledge of the genomic composition of a whole population.

Comprehensive genetic analysis reveals the mutational landscape of ABCA4-associated retinal dystrophy in a Chinese cohort

PURPOSE

To depict the variant profiles of the ABCA4 gene in a large Chinese cohort of patients with ABCA4-associated retinal dystrophy (ABCA4-RD).

METHODS

We recruited 290 unrelated Chinese patients with ABCA4-RD and did ABCA4 mutational screening by a combination of Sanger sequencing, targeted exome sequencing, entire ABCA4 locus sequencing, and whole genome sequencing (WGS). The pathogenicity of variants was assessed using in silico tools or in vitro splicing assays following the American College of Medical Genetics and Genomics guidelines.

RESULTS

Two hundred sixty-eight distinct pathogenic variants were identified, and 57 were novel. In 580 alleles, 22 noncoding region variants outside canonical splice sites and 4 structural variations were found in 44 alleles accounting for 7.6% of all alleles. Bioinformatics analysis showed the complex mechanism of aberrant splicing productsnatural splice site disruption, branch point destruction, and cryptic splice site activation. Correspondingly, minigene assays validated the various abnormal splicing products, including exon skipping, exon elongation, partial exon deletion, and pseudoexon insertion. WGS identified the first inversion variation in ABCA4.

CONCLUSIONS

This study systematically depicted the variant profiles of ABCA4 and revealed the missing alleles of patients with ABCA4-RD in a large Chinese cohort. Our findings demonstrated the complexity of molecular diagnosis of Mendelian diseases and the efficiency of WGS for detecting structural variants.

Rare copy-number variants as modulators of common disease susceptibility

BACKGROUND

Copy-number variations (CNVs) have been associated with rare and debilitating genomic disorders (GDs) but their impact on health later in life in the general population remains poorly described.

METHODS

Assessing four modes of CNV action, we performed genome-wide association scans (GWASs) between the copy-number of CNV-proxy probes and 60 curated ICD-10 based clinical diagnoses in 331,522 unrelated white British UK Biobank (UKBB) participants with replication in the Estonian Biobank.

RESULTS

We identified 73 signals involving 40 diseases, all of which indicating that CNVs increased disease risk and caused earlier onset. We estimated that 16% of these associations are indirect, acting by increasing body mass index (BMI). Signals mapped to 45 unique, non-overlapping regions, nine of which being linked to known GDs. Number and identity of genes affected by CNVs modulated their pathogenicity, with many associations being supported by colocalization with both common and rare single-nucleotide variant association signals. Dissection of association signals provided insights into the epidemiology of known gene-disease pairs (e.g., deletions in BRCA1 and LDLR increased risk for ovarian cancer and ischemic heart disease, respectively), clarified dosage mechanisms of action (e.g., both increased and decreased dosage of 17q12 impacted renal health), and identified putative causal genes (e.g., ABCC6 for kidney stones). Characterization of the pleiotropic pathological consequences of recurrent CNVs at 15q13, 16p13.11, 16p12.2, and 22q11.2 in adulthood indicated variable expressivity of these regions and the involvement of multiple genes. Finally, we show that while the total burden of rare CNVs-and especially deletions-strongly associated with disease risk, it only accounted for ~ 0.02% of the UKBB disease burden. These associations are mainly driven by CNVs at known GD CNV regions, whose pleiotropic effect on common diseases was broader than anticipated by our CNV-GWAS.

CONCLUSIONS

Our results shed light on the prominent role of rare CNVs in determining common disease susceptibility within the general population and provide actionable insights for anticipating later-onset comorbidities in carriers of recurrent CNVs.

Comprehensive and accurate genome analysis at scale using DRAGEN accelerated algorithms

Research and medical genomics require comprehensive and scalable solutions to drive the discovery of novel disease targets, evolutionary drivers, and genetic markers with clinical significance. This necessitates a framework to identify all types of variants independent of their size (e.g., SNV/SV) or location (e.g., repeats). Here we present DRAGEN that utilizes novel methods based on multigenomes, hardware acceleration, and machine learning based variant detection to provide novel insights into individual genomes with ~30min computation time (from raw reads to variant detection). DRAGEN outperforms all other state-of-the-art methods in speed and accuracy across all variant types (SNV, indel, STR, SV, CNV) and further incorporates specialized methods to obtain key insights in medically relevant genes (e.g., HLA, SMN, GBA). We showcase DRAGEN across 3,202 genomes and demonstrate its scalability, accuracy, and innovations to further advance the integration of comprehensive genomics for research and medical applications.

Detection of mosaic and population-level structural variants with Sniffles2

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

Gapless genome assembly and epigenetic profiles reveal gene regulation of whole-genome triplication in lettuce

BACKGROUND

Lettuce, an important member of the Asteraceae family, is a globally cultivated cash vegetable crop. With a highly complex genome (∼2.5 Gb; 2n = 18) rich in repeat sequences, current lettuce reference genomes exhibit thousands of gaps, impeding a comprehensive understanding of the lettuce genome.

FINDINGS

Here, we present a near-complete gapless reference genome for cutting lettuce with high transformability, using long-read PacBio HiFi and Nanopore sequencing data. In comparison to stem lettuce genome, we identify 127,681 structural variations (SVs, present in 0.41 Gb of sequence), reflecting the divergence of leafy and stem lettuce. Interestingly, these SVs are related to transposons and DNA methylation states. Furthermore, we identify 4,612 whole-genome triplication genes exhibiting high expression levels associated with low DNA methylation levels and high N6-methyladenosine RNA modifications. DNA methylation changes are also associated with activation of genes involved in callus formation.

CONCLUSIONS

Our gapless lettuce genome assembly, an unprecedented achievement in the Asteraceae family, establishes a solid foundation for functional genomics, epigenomics, and crop breeding and sheds new light on understanding the complexity of gene regulation associated with the dynamics of DNA and RNA epigenetics in genome evolution.

A graph clustering algorithm for detection and genotyping of structural variants from long reads

BACKGROUND

Structural variants (SVs) are genomic polymorphisms defined by their length (>50 bp). The usual types of SVs are deletions, insertions, translocations, inversions, and copy number variants. SV detection and genotyping is fundamental given the role of SVs in phenomena such as phenotypic variation and evolutionary events. Thus, methods to identify SVs using long-read sequencing data have been recently developed.

FINDINGS

We present an accurate and efficient algorithm to predict germline SVs from long-read sequencing data. The algorithm starts collecting evidence (signatures) of SVs from read alignments. Then, signatures are clustered based on a Euclidean graph with coordinates calculated from lengths and genomic positions. Clustering is performed by the DBSCAN algorithm, which provides the advantage of delimiting clusters with high resolution. Clusters are transformed into SVs and a Bayesian model allows to precisely genotype SVs based on their supporting evidence. This algorithm is integrated into the single sample variants detector of the Next Generation Sequencing Experience Platform, which facilitates the integration with other functionalities for genomics analysis. We performed multiple benchmark experiments, including simulation and real data, representing different genome profiles, sequencing technologies (PacBio HiFi, ONT), and read depths.

CONCLUSION

The results show that our approach outperformed state-of-the-art tools on germline SV calling and genotyping, especially at low depths, and in error-prone repetitive regions. We believe this work significantly contributes to the development of bioinformatic strategies to maximize the use of long-read sequencing technologies.

Evaluation of genetic risk of apparently balanced chromosomal rearrangement carriers by breakpoint characterization

PURPOSE

To report genetic characteristics and associated risk of chromosomal breaks due to chromosomal rearrangements in large samples.

METHODS

MicroSeq, a technique that combines chromosome microdissection and next-generation sequencing, was used to identify chromosomal breakpoints. Long-range PCR and Sanger sequencing were used to precisely characterize 100 breakpoints in 50 ABCR carriers.

RESULTS

In addition to the recurrent regions of balanced rearrangement breaks in 8q24.13, 11q11.23, and 22q11.21 that had been documented, we have discovered a 10-Mb region of 12q24.13-q24.3 that could potentially be a sparse region of balanced rearrangement breaks. We found that 898 breakpoints caused gene disruption and a total of 188 breakpoints interrupted genes recorded in OMIM. The percentage of breakpoints that disrupted autosomal dominant genes recorded in OMIM was 25.53% (48/188). Fifty-four of the precisely characterized breakpoints had 1-8-bp microhomologous sequences.

CONCLUSION

Our findings provide a reference for the evaluation of the pathogenicity of mutations in related genes that cause protein truncation in clinical practice. According to the characteristics of breakpoints, non-homologous end joining and microhomology-mediated break-induced replication may be the main mechanism for ABCRs formation.

Methods of Detection and Mechanisms of Origin of Complex Structural Genome Variations

Based on classical karyotyping, structural genome variations (SVs) have generally been considered to be either "simple" (with one or two breakpoints) or "complex" (with more than two breakpoints). Studying the breakpoints of SVs at nucleotide resolution revealed additional, subtle structural variations, such that even "simple" SVs turned out to be "complex." Genome-wide sequencing methods, such as fosmid and paired-end mapping, short-read and long-read whole genome sequencing, and single-molecule optical mapping, also indicated that the number of SVs per individual was considerably larger than expected from karyotyping and high-resolution chromosomal array-based studies. Interestingly, SVs were detected in studies of cohorts of individuals without clinical phenotypes. The common denominator of all SVs appears to be a failure to accurately repair DNA double-strand breaks (DSBs) or to halt cell cycle progression if DSBs persist. This review discusses the various DSB response mechanisms during the mitotic cell cycle and during meiosis and their regulation. Emphasis is given to the molecular mechanisms involved in the formation of translocations, deletions, duplications, and inversions during or shortly after meiosis I. Recently, CRISPR-Cas9 studies have provided unexpected insights into the formation of translocations and chromothripsis by both breakage-fusion-bridge and micronucleus-dependent mechanisms.

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.

Meiotic DNA breaks drive multifaceted mutagenesis in the human germ line

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

Most large structural variants in cancer genomes can be detected without long reads

Short-read sequencing is the workhorse of cancer genomics yet is thought to miss many structural variants (SVs), particularly large chromosomal alterations. To characterize missing SVs in short-read whole genomes, we analyzed 'loose ends'-local violations of mass balance between adjacent DNA segments. In the landscape of loose ends across 1,330 high-purity cancer whole genomes, most large (>10-kb) clonal SVs were fully resolved by short reads in the 87% of the human genome where copy number could be reliably measured. Some loose ends represent neotelomeres, which we propose as a hallmark of the alternative lengthening of telomeres phenotype. These pan-cancer findings were confirmed by long-molecule profiles of 38 breast cancer and melanoma cases. Our results indicate that aberrant homologous recombination is unlikely to drive the majority of large cancer SVs. Furthermore, analysis of mass balance in short-read whole genome data provides a surprisingly complete picture of cancer chromosomal structure.

Blinded rebiopsy and analysis of noneuploid embryos with 2 distinct preimplantation genetic testing platforms for aneuploidy

OBJECTIVE

To determine how often a noneuploid result from a single trophectoderm (TE) biopsy tested with the next-generation sequencing (NGS)-based preimplantation genetic testing for aneuploidy (PGT-A) is concordant with rebiopsies tested with a single-nucleotide polymorphism (SNP) array-based PGT-A platform.

DESIGN

Blinded prospective cohort study.

SETTING

University-affiliated fertility center.

PATIENT(S)

One hundred blastocysts were chosen from donated samples; on TE biopsy with NGS-based PGT-A, 40 had at least one whole chromosome full copy number aneuploidy alone, 20 had a single whole chromosome intermediate copy number ("whole chromosome mosaic"), 20 had a single full segmental aneuploidy (segA), and 20 had a single segmental intermediate copy number ("segmental mosaic").

INTERVENTIONS

Four rebiopsies were collected from each embryo: 3 TE biopsies and the remaining embryo. Each rebiopsy was randomized, blinded, and assessed with an SNP array-based PGT-A platform that combines copy number and allele ratio analyses, without mosaicism reporting.

MAIN OUTCOME MEASURE(S)

Concordance between the NGS result and rebiopsy results and within each embryo's blinded rebiopsy results.

RESULT(S)

Next-generation sequencing-diagnosed whole chromosome aneuploidy (WCA) was reconfirmed in 95% (95% confidence interval [CI], 83%-99%) of embryos; 2 embryos with NGS-diagnosed WCA were called euploid on all conclusive rebiopsies. Among embryos with NGS-diagnosed whole chromosome mosaicism, 35% (95% CI, 15%-59%) were called euploid and 15% (95% CI, 3%-38%) were called whole chromosome aneuploid on all conclusive rebiopsies. A total of 30% (95% CI, 12%-54%) of embryos with NGS-diagnosed segA and 65% (95% CI, 41%-85%) of embryos with NGS-diagnosed segmental mosaicism were called euploid on all conclusive rebiopsies. In total, 13% (95% CI, 6%-25%) of embryos with NGS-diagnosed full copy number aneuploidy and 50% (95% CI, 34%-66%) of embryos with NGS-diagnosed mosaicism had uniformly euploid SNP results. Conversely, all embryos with at least one noneuploid SNP result (n = 72) either had SNP-diagnosed aneuploidy on another rebiopsy from the same embryo or NGS-diagnosed aneuploidy/mosaicism involving the same chromosome.

CONCLUSION(S)

Next-generation sequencing-diagnosed WCA is highly concordant with rebiopsies tested with an SNP array-based PGT-A; however, whole chromosome mosaicism, segA, and segmental mosaicism are less concordant, reinforcing that embryos with these results may have reproductive potential and be suitable for transfer.

Gene duplication and deletion caused by over-replication at a fork barrier

Replication fork stalling can provoke fork reversal to form a four-way DNA junction. This remodelling of the replication fork can facilitate repair, aid bypass of DNA lesions, and enable replication restart, but may also pose a risk of over-replication during fork convergence. We show that replication fork stalling at a site-specific barrier in fission yeast can induce gene duplication-deletion rearrangements that are independent of replication restart-associated template switching and Rad51-dependent multi-invasion. Instead, they resemble targeted gene replacements (TGRs), requiring the DNA annealing activity of Rad52, the 3'-flap nuclease Rad16-Swi10, and mismatch repair protein Msh2. We propose that excess DNA, generated during the merging of a canonical fork with a reversed fork, can be liberated by a nuclease and integrated at an ectopic site via a TGR-like mechanism. This highlights how over-replication at replication termination sites can threaten genome stability in eukaryotes.

High population frequencies of MICA copy number variations originate from independent recombination events

MICA is a stress-induced ligand of the NKG2D receptor that stimulates NK and T cell responses and was identified as a key determinant of anti-tumor immunity. The MICA gene is located inside the MHC complex and is in strong linkage disequilibrium with HLA-B. While an HLA-B*48-linked MICA deletion-haplotype was previously described in Asian populations, little is known about other MICA copy number variations. Here, we report the genotyping of more than two million individuals revealing high frequencies of MICA duplications (1%) and MICA deletions (0.4%). Their prevalence differs between ethnic groups and can rise to 2.8% (Croatia) and 9.2% (Mexico), respectively. Targeted sequencing of more than 70 samples indicates that these copy number variations originate from independent nonallelic homologous recombination events between segmental duplications upstream of MICA and MICB. Overall, our data warrant further investigation of disease associations and consideration of MICA copy number data in oncological study protocols.

Comprehensive detection of structural variation and transposable element differences between wild type laboratory lineages of C. elegans

Genomic structural variations (SVs) and transposable elements (TEs) can be significant contributors to genome evolution, altered gene expression, and risk of genetic diseases. Recent advancements in long-read sequencing have greatly improved the quality of de novo genome assemblies and enhanced the detection of sequence variants at the scale of hundreds or thousands of bases. Comparisons between two diverged wild isolates of Caenorhabditis elegans, the Bristol and Hawaiian strains, have been widely utilized in the analysis of small genetic variations. Genetic drift, including SVs and rearrangements of repeated sequences such as TEs, can occur over time from long-term maintenance of wild type isolates within the laboratory. To comprehensively detect both large and small structural variations as well as TEs due to genetic drift, we generated de novo genome assemblies and annotations for each strain from our lab collection using both long- and short-read sequencing and compared our assemblies and annotations with that of other lab wild type strains. Within our lab assemblies, we annotate over 3.1Mb of sequence divergence between the Bristol and Hawaiian isolates: 337,584 SNPs, 94,503 small insertion-deletions (<50bp), and 4,334 structural variations (>50bp). Further, we define the location and movement of specific DNA TEs between N2 Bristol and CB4856 Hawaiian wild type isolates. Specifically, we find the N2 Bristol genome has 20.6% more TEs from the Tc1/mariner family than the CB4856 Hawaiian genome. Moreover, we identified Zator elements as the most abundant and mobile TE family in the genome. Using specific TE sequences with unique SNPs, we also identify 38 TEs that moved intrachromosomally and 9 TEs that moved interchromosomally between the N2 Bristol and CB4856 Hawaiian genomes. By comparing the de novo genome assembly of our lab collection Bristol isolate to the VC2010 Bristol assembly, we also reveal that lab lineages display over 2 Mb of total variation: 1,162 SNPs, 1,528 indels, and 897 SVs with 95% of the variation due to SVs. Overall, our work demonstrates the unique contribution of SVs and TEs to variation and genetic drift between wild type laboratory strains assumed to be isogenic despite growing evidence of genetic drift and phenotypic variation.

Chromosomal inversion polymorphisms are widespread across the species ranges of rough periwinkles (Littorina saxatilis and L. arcana)

Inversions are thought to play a key role in adaptation and speciation, suppressing recombination between diverging populations. Genes influencing adaptive traits cluster in inversions, and changes in inversion frequencies are associated with environmental differences. However, in many organisms, it is unclear if inversions are geographically and taxonomically widespread. The intertidal snail, Littorina saxatilis, is one such example. Strong associations between putative polymorphic inversions and phenotypic differences have been demonstrated between two ecotypes of L. saxatilis in Sweden and inferred elsewhere, but no direct evidence for inversion polymorphism currently exists across the species range. Using whole genome data from 107 snails, most inversion polymorphisms were found to be widespread across the species range. The frequencies of some inversion arrangements were significantly different among ecotypes, suggesting a parallel adaptive role. Many inversions were also polymorphic in the sister species, L. arcana, hinting at an ancient origin.

Aging, longevity, and the role of environmental stressors: a focus on wildfire smoke and air quality

Aging is a complex biological process involving multiple interacting mechanisms and is being increasingly linked to environmental exposures such as wildfire smoke. In this review, we detail the hallmarks of aging, emphasizing the role of telomere attrition, cellular senescence, epigenetic alterations, proteostasis, genomic instability, and mitochondrial dysfunction, while also exploring integrative hallmarks - altered intercellular communication and stem cell exhaustion. Within each hallmark of aging, our review explores how environmental disasters like wildfires, and their resultant inhaled toxicants, interact with these aging mechanisms. The intersection between aging and environmental exposures, especially high-concentration insults from wildfires, remains under-studied. Preliminary evidence, from our group and others, suggests that inhaled wildfire smoke can accelerate markers of neurological aging and reduce learning capabilities. This is likely mediated by the augmentation of circulatory factors that compromise vascular and blood-brain barrier integrity, induce chronic neuroinflammation, and promote age-associated proteinopathy-related outcomes. Moreover, wildfire smoke may induce a reduced metabolic, senescent cellular phenotype. Future interventions could potentially leverage combined anti-inflammatory and NAD + boosting compounds to counter these effects. This review underscores the critical need to study the intricate interplay between environmental factors and the biological mechanisms of aging to pave the way for effective interventions.

VariantSurvival: a tool to identify genotype-treatment response

Motivation: For a number of neurological diseases, such as Alzheimer's disease, amyotrophic lateral sclerosis, and many others, certain genes are known to be involved in the disease mechanism. A common question is whether a structural variant in any such gene may be related to drug response in clinical trials and how this relationship can contribute to the lifecycle of drug development. Results: To this end, we introduce VariantSurvival, a tool that identifies changes in survival relative to structural variants within target genes. VariantSurvival matches annotated structural variants with genes that are clinically relevant to neurological diseases. A Cox regression model determines the change in survival between the placebo and clinical trial groups with respect to the number of structural variants in the drug target genes. We demonstrate the functionality of our approach with the exemplary case of the SETX gene. VariantSurvival has a user-friendly and lightweight graphical user interface built on the shiny web application package.

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 .

Disruption of regulatory domains and novel transcripts as disease-causing mechanisms

Deletions, duplications, insertions, inversions, and translocations, collectively called structural variations (SVs), affect more base pairs of the genome than any other sequence variant. The recent technological advancements in genome sequencing have enabled the discovery of tens of thousands of SVs per human genome. These SVs primarily affect non-coding DNA sequences, but the difficulties in interpreting their impact limit our understanding of human disease etiology. The functional annotation of non-coding DNA sequences and methodologies to characterize their three-dimensional (3D) organization in the nucleus have greatly expanded our understanding of the basic mechanisms underlying gene regulation, thereby improving the interpretation of SVs for their pathogenic impact. Here, we discuss the various mechanisms by which SVs can result in altered gene regulation and how these mechanisms can result in rare genetic disorders. Beyond changing gene expression, SVs can produce novel gene-intergenic fusion transcripts at the SV breakpoints.

Demonstration of in vivo engineered tandem duplications of varying sizes using CRISPR and recombinases in Drosophila melanogaster

Tandem gene duplicates are important parts of eukaryotic genome structure, yet the phenotypic effects of new tandem duplications are not well-understood, in part owing to a lack of techniques to build and modify them. We introduce a method, Recombinase-Mediated Tandem Duplication, to engineer specific tandem duplications in vivo using CRISPR and recombinases. We describe construction of four different tandem duplications of the Alcohol Dehydrogenase (Adh) gene in Drosophila melanogaster, with duplicated block sizes ranging from 4.2 to 20.7 kb. Flies with the Adh duplications show elevated ADH enzyme activity over unduplicated single copies. This approach to engineering duplications is combinatoric, opening the door to systematic study of the relationship between the structure of tandem duplications and their effects on expression.