HMZDupFinder: a robust computational approach for detecting intragenic homozygous duplications from exome sequencing data.

Homozygous duplications contribute to genetic disease by altering gene dosage or disrupting gene regulation and can be more deleterious to organismal biology than heterozygous duplications. Intragenic exonic duplications can result in loss-of-function (LoF) or gain-of-function (GoF) alleles that when homozygosed, i.e. brought to homozygous state at a locus by identity by descent or state, could potentially result in autosomal recessive (AR) rare disease traits. However, the detection and functional interpretation of homozygous duplications from exome sequencing data remains a challenge. We developed a framework algorithm, HMZDupFinder, that is designed to detect exonic homozygous duplications from exome sequencing (ES) data. The HMZDupFinder algorithm can efficiently process large datasets and accurately identifies small intragenic duplications, including those associated with rare disease traits. HMZDupFinder called 965 homozygous duplications with three or less exons from 8,707 ES with a recall rate of 70.9% and a precision of 16.1%. We experimentally confirmed 8/10 rare homozygous duplications. Pathogenicity assessment of these copy number variant alleles allowed clinical genomics contextualization for three homozygous duplications alleles, including two affecting known OMIM disease genes EDAR (MIM# 224900), TNNT1(MIM# 605355), and one variant in a novel candidate disease gene: PAAF1.

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

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

Genome sequencing differentiates a paracentric inversion from a balanced insertion enabling more accurate preimplantation genetic testing.

INTRODUCTION: Distinguishing paracentric inversions (PAIs) from chromosomal insertions has traditionally relied on fluorescent in situ hybridization (FISH) techniques, but recent advancements in high-throughput sequencing have enabled the use of genome sequencing for such differentiation. In this study, we present a 38-year-old male carrier of a paracentric inversion on chromosome 2q, inv (2)(q31.2q34), whose partner experienced recurrent miscarriages. MATERIAL AND METHODS: FISH analysis confirmed the inversion, and genome sequencing was employed for detailed characterization. RESULTS: Preimplantation genetic testing (PGT) revealed that all assessed embryos were balanced, consistent with the low risk of unbalanced offspring associated with PAIs. While PAI carriers traditionally exhibit low risk of producing unbalanced offspring, exceptions exist due to crossover events within the inversion loop. Although the sample size was limited, the findings align with existing sperm study data, supporting the rare occurrence of unbalanced progeny in PAI carriers. CONCLUSIONS: This study highlights the possibility of characterizing PAIs using genome sequencing to enable correct reproductive counseling and PGT decisions. Detailed characterization of a PAI is crucial for understanding potential outcomes and guiding PGT strategies, as accurate knowledge of the inversion size is essential for appropriate method selection in PGT. Given the very low risk of unbalanced offspring in PAI carriers, routine PGT may not be warranted but should be considered in specific cases with a history of unbalanced progeny or recurrent miscarriages. This study contributes to our understanding of PAI segregation and its implications for reproductive outcomes.

Genome-wide maps of highly-similar intrachromosomal repeats that mediate ectopic recombination in three human genome assemblies.

Repeated sequences spread throughout the genome play important roles in shaping the structure of chromosomes and facilitating the generation of new genomic variation. Through a variety of mechanisms, repeats are involved in generating structural rearrangements such as deletions, duplications, inversions, and translocations, which can have the potential to impact human health. Despite their significance, repetitive regions including tandem repeats, transposable elements, segmental duplications, and low-copy repeats remain a challenge to characterize due to technological limitations inherent to many sequencing methodologies. We performed genome-wide analyses and comparisons of direct and inverted repeated sequences in the latest available human genome reference assemblies including GRCh37 and GRCh38 and the most recent telomere-to-telomere alternate assembly (T2T-CHM13). Overall, the composition and distribution of direct and inverted repeats identified remains similar among the three assemblies but we observed an increase in the number of repeated sequences detected in the T2T-CHM13 assembly versus the reference assemblies. As expected, there is an enrichment of repetitive regions in the short arms of acrocentric chromosomes, which had been previously unresolved in the human genome reference assemblies. We cross-referenced the identified repeats with protein-coding genes across the genome to identify those at risk for being involved in genomic disorders. We observed that certain gene categories, such as olfactory receptors and immune response genes, are enriched among those impacted by repeated sequences likely contributing to human diversity and adaptation. Through this analysis, we have produced a catalogue of direct and inversely oriented repeated sequences across the currently three most widely used human genome assemblies. Bioinformatic analyses of these repeats and their contribution to genome architecture can reveal regions that are most susceptible to genomic instability. Understanding how the architectural genomic features of repeat pairs such as their homology, size and distance can lead to complex genomic rearrangement formation can provide further insights into the molecular mechanisms leading to genomic disorders and genome evolution. Author summary: This study focused on the characterization of intrachromosomal repeated sequences in the human genome that can play important roles in shaping chromosome structure and generating new genomic variation in three human genome assemblies. We observed an increase in the number of repeated sequence pairs detected in the most recent telomere-to-telomere alternate assembly (T2T-CHM13) compared to the reference assemblies (GRCh37 and GRCh38). We observed an enrichment of repeats in the T2T-CHM13 acrocentric chromosomes, which had been previously unresolved. Importantly, our study provides a catalogue of direct and inverted repeated sequences across three commonly used human genome assemblies, which can aid in the understanding of genomic architecture instability, evolution, and disorders. Our analyses provide insights into repetitive regions in the human genome that may contribute to complex genomic rearrangements.

Generation of five induced pluripotent stem cell lines from patients with MECP2 Duplication Syndrome.

MECP2 Duplication Syndrome (MDS) is a rare, severe neurodevelopmental disorder arising from duplications in the Xq28 region containing the MECP2 gene that predominantly affects males. We generated five human induced pluripotent stem cell (iPSC) lines from the fibroblasts of individuals carrying between 0.355 and 11.2 Mb size duplications in the chromosomal locus containing MECP2. All lines underwent extensive testing to confirm MECP2 duplication and iPSC-related features such as morphology, pluripotency markers, and trilineage differentiation potential. These lines are a valuable resource for molecular and functional studies of MDS as well as screening for a variety of therapeutic approaches.

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.

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.

Closing the gap: Solving complex medically relevant genes at scale.

Comprehending the mechanism behind human diseases with an established heritable component represents the forefront of personalized medicine. Nevertheless, numerous medically important genes are inaccurately represented in short-read sequencing data analysis due to their complexity and repetitiveness or the so-called 'dark regions' of the human genome. The advent of PacBio as a long-read platform has provided new insights, yet HiFi whole-genome sequencing (WGS) cost remains frequently prohibitive. We introduce a targeted sequencing and analysis framework, Twist Alliance Dark Genes Panel (TADGP), designed to offer phased variants across 389 medically important yet complex autosomal genes. We highlight TADGP accuracy across eleven control samples and compare it to WGS. This demonstrates that TADGP achieves variant calling accuracy comparable to HiFi-WGS data, but at a fraction of the cost. Thus, enabling scalability and broad applicability for studying rare diseases or complementing previously sequenced samples to gain insights into these complex genes. TADGP revealed several candidate variants across all cases and provided insight into LPA diversity when tested on samples from rare disease and cardiovascular disease cohorts. In both cohorts, we identified novel variants affecting individual disease-associated genes (e.g., IKZF1, KCNE1). Nevertheless, the annotation of the variants across these 389 medically important genes remains challenging due to their underrepresentation in ClinVar and gnomAD. Consequently, we also offer an annotation resource to enhance the evaluation and prioritization of these variants. Overall, we can demonstrate that TADGP offers a cost-efficient and scalable approach to routinely assess the dark regions of the human genome with clinical relevance.

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.

Variant-specific pathophysiological mechanisms of AFF3 differently influence transcriptome profiles.

BACKGROUND: We previously described the KINSSHIP syndrome, an autosomal dominant disorder associated with intellectual disability (ID), mesomelic dysplasia and horseshoe kidney, caused by de novo variants in the degron of AFF3. Mouse knock-ins and overexpression in zebrafish provided evidence for a dominant-negative mode of action, wherein an increased level of AFF3 resulted in pathological effects. METHODS: Evolutionary constraints suggest that other modes-of-inheritance could be at play. We challenged this hypothesis by screening ID cohorts for individuals with predicted-to-be damaging variants in AFF3. We used both animal and cellular models to assess the deleteriousness of the identified variants. RESULTS: We identified an individual with a KINSSHIP-like phenotype carrying a de novo partial duplication of AFF3 further strengthening the hypothesis that an increased level of AFF3 is pathological. We also detected seventeen individuals displaying a milder syndrome with either heterozygous Loss-of-Function (LoF) or biallelic missense variants in AFF3. Consistent with semi-dominance, we discovered three patients with homozygous LoF and one compound heterozygote for a LoF and a missense variant, who presented more severe phenotypes than their heterozygous parents. Matching zebrafish knockdowns exhibit neurological defects that could be rescued by expressing human AFF3 mRNA, confirming their association with the ablation of aff3. Conversely, some of the human AFF3 mRNAs carrying missense variants identified in affected individuals did not rescue these phenotypes. Overexpression of mutated AFF3 mRNAs in zebrafish embryos produced a significant increase of abnormal larvae compared to wild-type overexpression further demonstrating deleteriousness. To further assess the effect of AFF3 variation, we profiled the transcriptome of fibroblasts from affected individuals and engineered isogenic cells harboring + / + , KINSSHIP/KINSSHIP, LoF/ + , LoF/LoF or KINSSHIP/LoF AFF3 genotypes. The expression of more than a third of the AFF3 bound loci is modified in either the KINSSHIP/KINSSHIP or the LoF/LoF lines. While the same pathways are affected, only about one third of the differentially expressed genes are common to the homozygote datasets, indicating that AFF3 LoF and KINSSHIP variants largely modulate transcriptomes differently, e.g. the DNA repair pathway displayed opposite modulation. CONCLUSIONS: Our results and the high pleiotropy shown by variation at this locus suggest that minute changes in AFF3 function are deleterious.

Variant-specific pathophysiological mechanisms of AFF3 differently influence transcriptome profiles.

Background: We previously described the KINSSHIP syndrome, an autosomal dominant disorder associated with intellectual disability (ID), mesomelic dysplasia and horseshoe kidney,caused by de novo variants in the degron of AFF3. Mouse knock-ins and overexpression in zebrafish provided evidence for a dominant-negative (DN) mode-of-action, wherein an increased level of AFF3 resulted in pathological effects. Methods: Evolutionary constraints suggest that other mode-of-inheritance could be at play. We challenged this hypothesis by screening ID cohorts for individuals with predicted-to-be deleterious variants in AFF3. We used both animal and cellular models to assess the deleteriousness of the identified variants. Results: We identified an individual with a KINSSHIP-like phenotype carrying a de novo partial duplication of AFF3 further strengthening the hypothesis that an increased level of AFF3 is pathological. We also detected seventeen individuals displaying a milder syndrome with either heterozygous LoF or biallelic missense variants in AFF3. Consistent with semi-dominance, we discovered three patients with homozygous LoF and one compound heterozygote for a LoF and a missense variant, who presented more severe phenotypes than their heterozygous parents. Matching zebrafish knockdowns exhibit neurological defects that could be rescued by expressing human AFF3 mRNA, confirming their association with the ablation of aff3. Conversely, some of the human AFF3 mRNAs carrying missense variants identified in affected individuals did not complement. Overexpression of mutated AFF3 mRNAs in zebrafish embryos produced a significant increase of abnormal larvae compared to wild-type overexpression further demonstrating deleteriousness. To further assess the effect of AFF3 variation, we profiled the transcriptome of fibroblasts from affected individuals and engineered isogenic cells harboring +/+, DN/DN, LoF/+, LoF/LoF or DN/LoF AFF3 genotypes. The expression of more than a third of the AFF3 bound loci is modified in either the DN/DN or the LoF/LoF lines. While the same pathways are affected, only about one-third of the differentially expressed genes are common to these homozygote datasets, indicating that AFF3 LoF and DN variants largely modulate transcriptomes differently, e.g. the DNA repair pathway displayed opposite modulation. Conclusions: Our results and the high pleiotropy shown by variation at this locus suggest that minute changes in AFF3 function are deleterious.

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

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