Genome-wide signatures of 'rearrangement hotspots' within segmental duplications in humans

High level of complexity and global diversity of the 3q29 locus revealed by optical mapping and long-read sequencing

BACKGROUND

High sequence identity between segmental duplications (SDs) can facilitate copy number variants (CNVs) via non-allelic homologous recombination (NAHR). These CNVs are one of the fundamental causes of genomic disorders such as the 3q29 deletion syndrome (del3q29S). There are 21 protein-coding genes lost or gained as a result of such recurrent 1.6-Mbp deletions or duplications, respectively, in the 3q29 locus. While NAHR plays a role in CNV occurrence, the factors that increase the risk of NAHR at this particular locus are not well understood.

METHODS

We employed an optical genome mapping technique to characterize the 3q29 locus in 161 unaffected individuals, 16 probands with del3q29S and their parents, and 2 probands with the 3q29 duplication syndrome (dup3q29S). Long-read sequencing-based haplotype resolved de novo assemblies from 44 unaffected individuals, and 1 trio was used for orthogonal validation of haplotypes and deletion breakpoints.

RESULTS

In total, we discovered 34 haplotypes, of which 19 were novel haplotypes. Among these 19 novel haplotypes, 18 were detected in unaffected individuals, while 1 novel haplotype was detected on the parent-of-origin chromosome of a proband with the del3q29S. Phased assemblies from 44 unaffected individuals enabled the orthogonal validation of 20 haplotypes. In 89% (16/18) of the probands, breakpoints were confined to paralogous copies of a 20-kbp segment within the 3q29 SDs. In one del3q29S proband, the breakpoint was confined to a 374-bp region using long-read sequencing. Furthermore, we categorized del3q29S cases into three classes and dup3q29S cases into two classes based on breakpoints. Finally, we found no evidence of inversions in parent-of-origin chromosomes.

CONCLUSIONS

We have generated the most comprehensive haplotype map for the 3q29 locus using unaffected individuals, probands with del3q29S or dup3q29S, and available parents, and also determined the deletion breakpoint to be within a 374-bp region in one proband with del3q29S. These results should provide a better understanding of the underlying genetic architecture that contributes to the etiology of del3q29S and dup3q29S.

Construction of copy number variation landscape and characterization of associated genes in a Bangladeshi cohort of neurodevelopmental disorders

Introduction: Copy number variations (CNVs) play a critical role in the pathogenesis of neurodevelopmental disorders (NDD) among children. In this study, we aim to identify clinically relevant CNVs, genes and their phenotypic characteristics in an ethnically underrepresented homogenous population of Bangladesh. Methods: We have conducted chromosomal microarray analysis (CMA) for 212 NDD patients with male to female ratio of 2.2:1.0 to identify rare CNVs. To identify candidate genes within the rare CNVs, gene constraint metrics [i.e., "Critical-Exon Genes (CEGs)"] were applied to the population data. Autism Diagnostic Observation Schedule-Second Edition (ADOS-2) was followed in a subset of 95 NDD patients to assess the severity of autism and all statistical tests were performed using the R package. Results: Of all the samples assayed, 12.26% (26/212) and 57.08% (121/212) patients carried pathogenic and variant of uncertain significance (VOUS) CNVs, respectively. While 2.83% (6/212) patients' pathogenic CNVs were found to be located in the subtelomeric regions. Further burden test identified females are significant carriers of pathogenic CNVs compared to males (OR = 4.2; p = 0.0007). We have observed an increased number of Loss of heterozygosity (LOH) within cases with 23.85% (26/109) consanguineous parents. Our analyses on imprinting genes show, 36 LOH variants disrupting 69 unique imprinted genes and classified these variants as VOUS. ADOS-2 subset shows severe social communication deficit (p = 0.014) and overall ASD symptoms severity (p = 0.026) among the patients carrying duplication CNV compared to the CNV negative group. Candidate gene analysis identified 153 unique CEGs in pathogenic CNVs and 31 in VOUS. Of the unique genes, 18 genes were found to be in smaller (<1 MB) focal CNVs in our NDD cohort and we identified PSMC3 gene as a strong candidate gene for Autism Spectrum Disorder (ASD). Moreover, we hypothesized that KMT2B gene duplication might be associated with intellectual disability. Conclusion: Our results show the utility of CMA for precise genetic diagnosis and its integration into the diagnosis, therapy and management of NDD patients.

Novel copy number variations within SYCE1 caused meiotic arrest and non-obstructive azoospermia

Background

Non-obstructive azoospermia (NOA) is the most severe disease in male infertility, but the genetic causes for majority of NOA remain unknown.

Methods

Two Chinese NOA-affected patients were recruited to identify the genetic causal factor of infertility. Whole-exome sequencing (WES) was conducted in the two patients with NOA. Sanger sequencing and CNV array were used to ascertain the WES results. Hematoxylin and eosin (H&E) staining and immunofluorescence (IF) were carried out to evaluate the stage of spermatogenesis arrested in the affected cases.

Results

Novel heterozygous deletion (LOH) within SYCE1 (seq[GRCh37] del(10)(10q26.3)chr10:g.135111754_135427143del) and heterozygous loss of function (LoF) variant in SYCE1 (NM_001143763: c.689_690 del:p.F230fs) were identified in one NOA-affected patient. While homozygous deletion within SYCE1 (seq[GRCh37] del(10)(10q26.3)chr10:g.135340247_135379115del) was detected in the other patient with meiotic arrest. H&E and IF staining demonstrated that the spermatogenesis was arrested at pachytene stage in the two patients with NOA, suggesting these two novel CNVs within SYCE1 could lead to meiotic defect and NOA.

Conclusions

We identified that two novel CNVs within SYCE1 are associated with meiotic arrest and male infertility. Thus, our study expands the knowledge of variants in SYCE1 and provides a new insight to understand the genetic etiologies of NOA.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01288-8.

Recurrent inversion polymorphisms in humans associate with genetic instability and genomic disorders

Unlike copy number variants (CNVs), inversions remain an underexplored genetic variation class. By integrating multiple genomic technologies, we discover 729 inversions in 41 human genomes. Approximately 85% of inversions <2 kbp form by twin-priming during L1 retrotransposition; 80% of the larger inversions are balanced and affect twice as many nucleotides as CNVs. Balanced inversions show an excess of common variants, and 72% are flanked by segmental duplications (SDs) or retrotransposons. Since flanking repeats promote non-allelic homologous recombination, we developed complementary approaches to identify recurrent inversion formation. We describe 40 recurrent inversions encompassing 0.6% of the genome, showing inversion rates up to 2.7 × 10^-4 per locus per generation. Recurrent inversions exhibit a sex-chromosomal bias and co-localize with genomic disorder critical regions. We propose that inversion recurrence results in an elevated number of heterozygous carriers and structural SD diversity, which increases mutability in the population and predisposes specific haplotypes to disease-causing CNVs.

A rare familial rearrangement of chromosomes 9 and 15 associated with intellectual disability: a clinical and molecular study

Background

There are many reports on rearrangements occurring separately in the regions of chromosomes 9p and 15q affected in the case under study. 15q duplication syndrome is caused by the presence of at least one extra maternally derived copy of the Prader–Willi/Angelman critical region. Trisomy 9p is the fourth most frequent chromosome anomaly with a clinically recognizable syndrome often accompanied by intellectual disability. Here we report a new case of a patient with maternally derived unique complex sSMC resulting in partial trisomy of both chromosomes 9 and 15 associated with intellectual disability.

Case presentation

We characterise a supernumerary derivative chromosome 15: 47,XY,+der(15)t(9;15)(p21.2;q13.2), likely resulting from 3:1 malsegregation during maternal gametogenesis. Chromosomal analysis showed that a phenotypically normal mother is a carrier of balanced translocation t(9;15)(p21.1;q13.2). Her 7-year-old son showed signs of intellectual disability and a number of physical abnormalities including bilateral cryptorchidism and congenital megaureter. The child’s magnetic resonance imaging showed changes in brain volume and in structural and functional connectivity revealing phenotypic changes caused by the presence of the extra chromosome material, whereas the mother’s brain MRI was normal. Sequence analyses of the microdissected der(15) chromosome detected two breakpoint regions: HSA9:25,928,021-26,157,441 (9p21.2 band) and HSA15:30,552,104-30,765,905 (15q13.2 band). The breakpoint region on chromosome HSA9 is poor in genetic features with several areas of high homology with the breakpoint region on chromosome 15. The breakpoint region on HSA15 is located in the area of a large segmental duplication.

Conclusions

We discuss the case of these phenotypic and brain MRI features in light of reported signatures for 9p partial trisomy and 15 duplication syndromes and analyze how the genomic characteristics of the found breakpoint regions have contributed to the origin of the derivative chromosome. We recommend MRI for all patients with a developmental delay, especially in cases with identified rearrangements, to accumulate more information on brain phenotypes related to chromosomal syndromes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13039-021-00565-y.

Gonadal mosaicism of large terminal de novo duplication and deletion in siblings with variable intellectual disability phenotypes

Background

Intellectual disability (ID) is a complex condition that can impact multiple domains of development. The genetic contribution to ID’s etiology is significant, with more than 100 implicated genes and loci currently identified. The majority of such variants are rare and de novo genetic mutations.

Methods

We have applied whole‐genome microarray to identify large, rare, clinically relevant copy number variants (CNVs). We have applied well‐established algorithms for variants call. Quantitative polymerase chain reaction (qPCR) was applied to validate the variants using three technical replicates for each family member. To assess whether the copy number variation was due to balanced translocation or mosaicism, we further conducted droplet digital PCR (ddPCR) on the whole family. We have, as well, applied “critical‐exon” mapping, human developmental brain transcriptome, and a database of known associated neurodevelopmental disorder variants to identify candidate genes.

Results

Here we present two siblings who are both impacted by a large terminal duplication and a deletion. Whole‐genome microarray revealed an 18.82 megabase (MB) duplication at terminal locus (7q34‐q36.3) of chromosome 7 and a 3.90 MB deletion impacting the terminal locus (15q26.3) of chromosome 15. qPCR and ddPCR experiments confirmed the de novo origin of the variants and the co‐occurrence of these two de novo events among the siblings, but their absence in both parents, implicates an unbalanced translocation that could have mal‐segregated among the siblings or a possible germline mosaicism. These terminal events impact IGF1R, CNTNAP2, and DPP6, shown to be strongly associated with neurodevelopmental disorders. Detailed clinical examination of the siblings revealed the presence of both shared and distinct phenotypic features.

Conclusions

This study identified two large rare terminal de novo events impacting two siblings. Further phenotypic investigation highlights that even in the presence of identical large high penetrant variants, spectrum of clinical features can be different between the siblings.

Interplay of interlocus gene conversion and crossover in segmental duplications under a neutral scenario

Interlocus gene conversion is a major evolutionary force that drives the concerted evolution of duplicated genomic regions. Theoretical models successfully have addressed the effects of interlocus gene conversion and the importance of crossover in the evolutionary fate of gene families and duplications but have not considered complex recombination scenarios, such as the presence of hotspots. To study the interplay between interlocus gene conversion and crossover, we have developed a forward-time simulator that allows the exploration of a wide range of interlocus gene conversion rates under different crossover models. Using it, we have analyzed patterns of nucleotide variation and linkage disequilibrium within and between duplicate regions, focusing on a neutral scenario with constant population size and validating our results with the existing theoretical models. We show that the interaction of gene conversion and crossover is nontrivial and that the location of crossover junctions is a fundamental determinant of levels of variation and linkage disequilibrium in duplicated regions. We also show that if crossover activity between duplications is strong enough, recurrent interlocus gene conversion events can break linkage disequilibrium within duplicates. Given the complex nature of interlocus gene conversion and crossover, we provide a framework to explore their interplay to help increase knowledge on molecular evolution within segmental duplications under more complex scenarios, such as demographic changes or natural selection.

The missing story behind Genome Wide Association Studies: single nucleotide polymorphisms in gene deserts have a story to tell

Genome wide association studies are central to the evolution of personalized medicine. However, the propensity for single nucleotide polymorphisms (SNPs) to fall outside of genes means that understanding how these polymorphisms alter cellular function requires an expanded view of human genetics. Integrating the study of genome structure (chromosome conformation capture) into its function opens up new avenues of exploration. Changes in the epigenome associated with SNPs in gene deserts will allow us to define complex diseases in a much clearer manner, and usher in a new era of disease pathway exploration.

UGT2B17 copy number gain in a large ankylosing spondylitis multiplex family

Background

The primary objective of this study is to identify novel copy number variations (CNVs) associated with familial ankylosing spondylitis (AS). A customized genome-wide microarray was designed to detect CNVs and applied to a multiplex AS family with six (6) affected family members. CNVs were detected using the built-in DNA analytics aberration detection method-2 (ADM-2) algorithm. Gene enrichment analysis was performed to observe the segregation. Subsequent validation was performed using real time quantitative fluorescence polymerase reaction (QF-PCR). The frequency of copy number variation for the UGT2B17 gene was then performed on two well-defined AS cohorts. Fisher exact test was performed to quantify the association.

Results

Our family-based analysis revealed ten gene-enriched CNVs that segregate with all six family members affected with AS. Based on the proposed function and the polymorphic nature of the UGT2B17 gene, the UGT2B17 gene CNV was selected for validation using real time QF-PCR with full concordance. The frequency of two copies of the UGT2B17 gene CNV was 0.41 in the Newfoundland AS cases and 0.35 in the Newfoundland controls (OR = 1.26(0.99-1.59); p < 0.05)), whereas the frequency of two (2) copies of the UGT2B17 gene CNV was 0.40 in the Alberta AS cases and 0.39 in the Alberta controls (OR = 1.05(95% CI: 0.83-1.33); p < 0.71)).

Conclusions

A genome-wide microarray interrogation of a large multiplex AS family revealed segregation of the UGT2B17 gene CNV among all affected family members. The association of the UGT2B17 CNV with AS is particularly interesting given the recent association of this CNV with osteoporosis and the proposed function as it encodes a key enzyme that inhibits androgens. However, two copies of the UGT2B17 gene CNV were only marginally significant in a uniplex AS cohort from Newfoundland but not in a uniplex AS cohort from Alberta.