Mechanisms for the Generation of Two Quadruplications Associated with Split-Hand Malformation

BHLHA9 homozygous duplication in a consanguineous family: A challenge for genetic counseling

Split-hand/foot malformation (SHFM) with long-bone deficiency (SHFLD) is a rare condition characterized by SHFM associated with long-bone malformation usually involving the tibia. It includes three different types; SHFLD1 (MIM % 119,100), SHFLD2 (MIM % 610,685) and SHFLD3 (MIM # 612576). The latter was shown to be the most commonly reported with a duplication in the 17p13.1p13.3 locus that was narrowed down to the BHLHA9 gene. Here, we report a consanguineous Lebanese family with three members presenting with limb abnormalities including tibial hemimelia. One of these patients presented with additional bowing fibula and another with bilateral split hand. CGH array analysis followed by RQ-PCR allowed us to detect the first homozygous duplication on the short arm of chromosome 17p13.3 including the BHLHA9 gene and involved in SHFLD3. Interestingly, one patient with the homozygous duplicated region, carrying thus four BHLHA9 copies presented with long bone deficiency but no SHFM. The incomplete penetrance and the variable expressivity of the disease in this family as well as the presence of the BHLHA9 homozygous duplication rendered genetic counseling extremely challenging and preimplantation genetic diagnosis almost impossible.

Long-Read Sequencing Revealed Extragenic and Intragenic Duplications of Exons 56-61 in DMD in an Asymptomatic Male and a DMD Patient

Expanded carrier screening (ECS) has become an increasingly common technique to assess the genetic risks of individuals in the prenatal or preconception period. Unexpected variants unrelated to referral are being increasingly detected in asymptomatic individuals through ECS. In this study, we reported an asymptomatic male with duplication of exons 56-61 in the DMD gene through ECS using whole-exome sequencing (WES), which was also detected in a male patient diagnosed with typical Duchenne muscular dystrophy (DMD). Breakpoint analysis was then performed to explore the potential mechanisms of phenotypic differences using long-read sequencing (LRS), PacBio single-molecule real-time (PacBio SMRT) target sequencing, and Sanger sequencing. Complex structural variations (SVs) on chromosome X were identified in the asymptomatic male, which revealed that the duplication occurred outside the DMD gene; whereas, the duplication in the patient with DMD was a tandem repeat. The phenotypic differences between the two men could be explained by the different breakpoint junctions. To the best of our knowledge, this is the first report of a breakpoint analysis of DMD duplication in two men with different phenotypes. Breakpoint analysis is necessary when the clinical phenotypes are inconsistent with genotypes, and it applies to prenatal testing.

Integrated sequencing and array comparative genomic hybridization in familial Parkinson disease

OBJECTIVE

To determine how single nucleotide variants (SNVs) and copy number variants (CNVs) contribute to molecular diagnosis in familial Parkinson disease (PD), we integrated exome sequencing (ES) and genome-wide array-based comparative genomic hybridization (aCGH) and further probed CNV structure to reveal mutational mechanisms.

METHODS

We performed ES on 110 subjects with PD and a positive family history; 99 subjects were also evaluated using genome-wide aCGH. We interrogated ES and aCGH data for pathogenic SNVs and CNVs at Mendelian PD gene loci. We confirmed SNVs via Sanger sequencing and further characterized CNVs with custom-designed high-density aCGH, droplet digital PCR, and breakpoint sequencing.

RESULTS

Using ES, we discovered individuals with known pathogenic SNVs in GBA (p.Glu365Lys, p.Thr408Met, p.Asn409Ser, and p.Leu483Pro) and LRRK2 (p.Arg1441Gly and p.Gly2019Ser). Two subjects were each double heterozygotes for variants in GBA and LRRK2. Based on aCGH, we additionally discovered cases with an SNCA duplication and heterozygous intragenic GBA deletion. Five additional subjects harbored both SNVs (p.Asn52Metfs*29, p.Thr240Met, p.Pro437Leu, and p.Trp453*) and likely disrupting CNVs at the PRKN locus, consistent with compound heterozygosity. In nearly all cases, breakpoint sequencing revealed microhomology, a mutational signature consistent with CNV formation due to DNA replication errors.

CONCLUSIONS

Integrated ES and aCGH yielded a genetic diagnosis in 19.3% of our familial PD cohort. Our analyses highlight potential mechanisms for SNCA and PRKN CNV formation, uncover multilocus pathogenic variation, and identify novel SNVs and CNVs for further investigation as potential PD risk alleles.

Breakpoint junction analysis for complex genomic rearrangements with the caldera volcano-like pattern

Chromosomal triplications can be classified into recurrent and nonrecurrent triplications. Most of the nonrecurrent triplications are embedded in duplicated segments, and duplication-inverted triplication-duplication (DUP-TRP/INV-DUP) has been established as one of the mechanisms of triplication. This study aimed to reveal the underlying mechanism of the TRP-DUP-TRP pattern of chromosomal aberrations, in which the appearance of moving averages obtained through array-based comparative genomic hybridization analysis is similar to the shadows of the caldera volcano-like pattern, which were first identified in two patients with neurodevelopmental disabilities. For this purpose, whole-genome sequencing using long-read Nanopore sequencing was carried out to confirm breakpoint junctions. Custom array analysis and Sanger sequencing were also used to detect all breakpoint junctions. As a result, the TRP-DUP-TRP pattern consisted of only two patterns of breakpoint junctions in both patients. In patient 1, microhomologies were identified in breakpoint junctions. In patient 2, more complex architectures with insertional segments were identified. Thus, replication-based mechanisms were considered as a mechanism of the TRP-DUP-TRP pattern.

Established and Novel Mechanisms Leading to de novo Genomic Rearrangements in the Human Germline

During gametogenesis, the human genome can acquire various de novo rearrangements. Most constitutional genomic rearrangements are created through 1 of the 4 well-known mechanisms, i.e., nonallelic homologous recombination, erroneous repair after double-strand DNA breaks, replication errors, and retrotransposition. However, recent studies have identified 2 types of extremely complex rearrangements that cannot be simply explained by these mechanisms. The first type consists of chaotic structural changes in 1 or a few chromosomes that result from “chromoanagenesis (an umbrella term that covers chromothripsis, chromoanasynthesis, and chromoplexy).” The other type is large independent rearrangements in multiple chromosomes indicative of “transient multifocal genomic crisis.” Germline chromoanagenesis (chromothripsis) likely occurs predominantly during spermatogenesis or postzygotic embryogenesis, while multifocal genomic crisis appears to be limited to a specific time window during oogenesis and early embryogenesis or during spermatogenesis. This review article introduces the current understanding of the molecular basis of de novo rearrangements in the germline.

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

Background

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

Methods

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

Results

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

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

Conclusion

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

Interchromosomal template-switching as a novel molecular mechanism for imprinting perturbations associated with Temple syndrome

Background

Intrachromosomal triplications (TRP) can contribute to disease etiology via gene dosage effects, gene disruption, position effects, or fusion gene formation. Recently, post-zygotic de novo triplications adjacent to copy-number neutral genomic intervals with runs of homozygosity (ROH) have been shown to result in uniparental isodisomy (UPD). The genomic structure of these complex genomic rearrangements (CGRs) shows a consistent pattern of an inverted triplication flanked by duplications (DUP-TRP/INV-DUP) formed by an iterative DNA replisome template-switching mechanism during replicative repair of a single-ended, double-stranded DNA (seDNA), the ROH results from an interhomolog or nonsister chromatid template switch. It has been postulated that these CGRs may lead to genetic abnormalities in carriers due to dosage-sensitive genes mapping within the copy-number variant regions, homozygosity for alleles at a locus causing an autosomal recessive (AR) disease trait within the ROH region, or imprinting-associated diseases.

Methods

Here, we report a family wherein the affected subject carries a de novo 2.2-Mb TRP followed by 42.2 Mb of ROH and manifests clinical features overlapping with those observed in association with chromosome 14 maternal UPD (UPD(14)mat). UPD(14)mat can cause clinical phenotypic features enabling a diagnosis of Temple syndrome. This CGR was then molecularly characterized by high-density custom aCGH, genome-wide single-nucleotide polymorphism (SNP) and methylation arrays, exome sequencing (ES), and the Oxford Nanopore long-read sequencing technology.

Results

We confirmed the postulated DUP-TRP/INV-DUP structure by multiple orthogonal genomic technologies in the proband. The methylation status of known differentially methylated regions (DMRs) on chromosome 14 revealed that the subject shows the typical methylation pattern of UPD(14)mat. Consistent with these molecular findings, the clinical features overlap with those observed in Temple syndrome, including speech delay.

Conclusions

These data provide experimental evidence that, in humans, triplication can lead to segmental UPD and imprinting disease. Importantly, genotype/phenotype analyses further reveal how a post-zygotically generated complex structural variant, resulting from a replication-based mutational mechanism, contributes to expanding the clinical phenotype of known genetic syndromes. Mechanistically, such events can distort transmission genetics resulting in homozygosity at a locus for which only one parent is a carrier as well as cause imprinting diseases.

Electronic supplementary material

The online version of this article (10.1186/s13073-019-0633-y) contains supplementary material, which is available to authorized users.

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

Background

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

Methods

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

Results

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

Conclusions

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

Electronic supplementary material

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

The coexistence of copy number variations (CNVs) and single nucleotide polymorphisms (SNPs) at a locus can result in distorted calculations of the significance in associating SNPs to disease

With the recent advance in genome-wide association studies (GWAS), disease-associated single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) have been extensively reported. Accordingly, the issue of incorrect identification of recombination events that can induce the distortion of multi-allelic or hemizygous variants has received more attention. However, the potential distorted calculation bias or significance of a detected association in a GWAS due to the coexistence of CNVs and SNPs in the same genomic region may remain under-recognized. Here we performed the association study within a congenital scoliosis (CS) cohort whose genetic etiology was recently elucidated as a compound inheritance model, including mostly one rare variant deletion CNV null allele and one common variant non-coding hypomorphic haplotype of the TBX6 gene. We demonstrated that the existence of a deletion in TBX6 led to an overestimation of the contribution of the SNPs on the hypomorphic allele. Furthermore, we generalized a model to explain the calculation bias, or distorted significance calculation for an association study, that can be 'induced' by CNVs at a locus. Meanwhile, overlapping between the disease-associated SNPs from published GWAS and common CNVs (overlap 10%) and pathogenic/likely pathogenic CNVs (overlap 99.69%) was significantly higher than the random distribution (p < 1 × 10^-6 and p = 0.034, respectively), indicating that such co-existence of CNV and SNV alleles might generally influence data interpretation and potential outcomes of a GWAS. We also verified and assessed the influence of colocalizing CNVs to the detection sensitivity of disease-associated SNP variant alleles in another adolescent idiopathic scoliosis (AIS) genome-wide association study. We proposed that detecting co-existent CNVs when evaluating the association signals between SNPs and disease traits could improve genetic model analyses and better integrate GWAS with robust Mendelian principles.

Marker chromosome genomic structure and temporal origin implicate a chromoanasynthesis event in a family with pleiotropic psychiatric phenotypes

Small supernumerary marker chromosomes (sSMC) are chromosomal fragments difficult to characterize genomically. Here, we detail a proband with schizoaffective disorder and a mother with bipolar disorder with psychotic features who present with a marker chromosome that segregates with disease. We explored the architecture of this marker and investigated its temporal origin. Array comparative genomic hybridization (aCGH) analysis revealed three duplications and three triplications that spanned the short arm of chromosome 9, suggestive of a chromoanasynthesis-like event. Segregation of marker genotypes, phased using sSMC mosaicism in the mother, provided evidence that it was generated during a germline-level event in the proband's maternal grandmother. Whole-genome sequencing (WGS) was performed to resolve the structure and junctions of the chromosomal fragments, revealing further complexities. While structural variations have been previously associated with neuropsychiatric disorders and marker chromosomes, here we detail the precise architecture, human life-cycle genesis, and propose a DNA replicative/repair mechanism underlying formation.

Genomic disorders 20 years on-mechanisms for clinical manifestations

Genomic disorders result from copy-number variants (CNVs) or submicroscopic rearrangements of the genome rather than from single nucleotide variants (SNVs). Diverse technologies, including array comparative genomic hybridization (aCGH) and single nucleotide polymorphism (SNP) microarrays, and more recently, whole genome sequencing and whole-exome sequencing, have enabled robust genome-wide unbiased detection of CNVs in affected individuals and in reportedly healthy controls. Sequencing of breakpoint junctions has allowed for elucidation of upstream mechanisms leading to genomic instability and resultant structural variation, whereas studies of the association between CNVs and specific diseases or susceptibility to morbid traits have enhanced our understanding of the downstream effects. In this review, we discuss the hallmarks of genomic disorders as they were defined during the first decade of the field, including genomic instability and the mechanism for rearrangement defined as nonallelic homologous recombination (NAHR); recurrent vs nonrecurrent rearrangements; and gene dosage sensitivity. Moreover, we highlight the exciting advances of the second decade of this field, including a deeper understanding of genomic instability and the mechanisms underlying complex rearrangements, mechanisms for constitutional and somatic chromosomal rearrangements, structural intra-species polymorphisms and susceptibility to NAHR, the role of CNVs in the context of genome-wide copy number and single nucleotide variation, and the contribution of noncoding CNVs to human disease.

Applicability of digital PCR to the investigation of pediatric-onset genetic disorders

Early-onset rare diseases have a strong impact on child healthcare even though the incidence of each of these diseases is relatively low. In order to better manage the care of these children, it is imperative to quickly diagnose the molecular bases for these disorders as well as to develop technologies with prognostic potential. Digital PCR (dPCR) is well suited for this role by providing an absolute quantification of the target DNA within a sample. This review illustrates how dPCR can be used to identify genes associated with pediatric-onset disorders, to identify copy number status of important disease-causing genes and variants and to quantify modifier genes. It is also a powerful technology to track changes in genomic biomarkers with disease progression. Based on its capability to accurately and reliably detect genomic alterations with high sensitivity and a large dynamic detection range, dPCR has the potential to become the tool of choice for the verification of pediatric disease-associated mutations identified by next generation sequencing, copy number determination and noninvasive prenatal screening.

Somatic mosaicism underlies X-linked acrogigantism syndrome in sporadic male subjects

Somatic mosaicism has been implicated as a causative mechanism in a number of genetic and genomic disorders. X-linked acrogigantism (XLAG) syndrome is a recently characterized genomic form of pediatric gigantism due to aggressive pituitary tumors that is caused by submicroscopic chromosome Xq26.3 duplications that include GPR101 We studied XLAG syndrome patients (n= 18) to determine if somatic mosaicism contributed to the genomic pathophysiology. Eighteen subjects with XLAG syndrome caused by Xq26.3 duplications were identified using high-definition array comparative genomic hybridization (HD-aCGH). We noted that males with XLAG had a decreased log2ratio (LR) compared with expected values, suggesting potential mosaicism, whereas females showed no such decrease. Compared with familial male XLAG cases, sporadic males had more marked evidence for mosaicism, with levels of Xq26.3 duplication between 16.1 and 53.8%. These characteristics were replicated using a novel, personalized breakpoint junction-specific quantification droplet digital polymerase chain reaction (ddPCR) technique. Using a separate ddPCR technique, we studied the feasibility of identifying XLAG syndrome cases in a distinct patient population of 64 unrelated subjects with acromegaly/gigantism, and identified one female gigantism patient who had had increased copy number variation (CNV) threshold for GPR101 that was subsequently diagnosed as having XLAG syndrome on HD-aCGH. Employing a combination of HD-aCGH and novel ddPCR approaches, we have demonstrated, for the first time, that XLAG syndrome can be caused by variable degrees of somatic mosaicism for duplications at chromosome Xq26.3. Somatic mosaicism was shown to occur in sporadic males but not in females with XLAG syndrome, although the clinical characteristics of the disease were similarly severe in both sexes.