The Alu-rich genomic architecture of SPAST predisposes to diverse and functionally distinct disease-associated CNV alleles

Missing genetic variations in GNE myopathy: rearrangement hotspots encompassing 5'UTR and founder allele

GNE myopathy is an autosomal recessive distal myopathy caused by loss-of-function mutations in the GNE gene, which encodes UDP-GlcNAc 2-epimerase/ManNAc kinase (GNE), a key enzyme in sialic-acid biosynthesis. By comprehensive screening of manifesting patients using a fine-mapped targeted next-generation sequencing (NGS), we identified copy number variations (CNVs) in 13 patients from 11 unrelated families. The nine unique CNVs largely vary in size from 0.3 to 72 kb. Over half of the cases carry different deletions spanning merely exon 2, which contains the 5' untranslated region (5'UTR) of the muscle major transcript hGNE1. Of most unique CNVs, either the telomeric or the centromeric breakpoint locates within intron 2, indicating rearrangement hotspots. Haplotype analysis suggested the existence of a founder allele with exon 2 deletion. The breakpoints for all CNVs were determined by long-range PCR and sequencing. All of the breakpoints of gross deletion/duplications reside within directly oriented pairs of Alu repeats. The results of this study firstly widen the spectra of mutations to CNVs encompassing 5'UTR, underscoring the pivotal role of the hGNE1 transcript. Alu-mediated non-recurrent CNVs may have been overlooked in a wide variety of recessive phenotypes, especially in those associated with genomic Alu-rich genes such as GNE.

Major influence of repetitive elements on disease-associated copy number variants (CNVs)

Copy number variants (CNVs) are important contributors to the human pathogenic genetic diversity as demonstrated by a number of cases reported in the literature. The high homology between repetitive elements may guide genomic stability which will give rise to CNVs either by non-allelic homologous recombination (NAHR) or non-homologous end joining (NHEJ). Here, we present a short guide based on previously documented cases of disease-associated CNVs in order to provide a general view on the impact of repeated elements on the stability of the genomic sequence and consequently in the origin of the human pathogenic variome.

Copy-Number Variation Contributes to the Mutational Load of Bardet-Biedl Syndrome

Bardet-Biedl syndrome (BBS) is a defining ciliopathy, notable for extensive allelic and genetic heterogeneity, almost all of which has been identified through sequencing. Recent data have suggested that copy-number variants (CNVs) also contribute to BBS. We used a custom oligonucleotide array comparative genomic hybridization (aCGH) covering 20 genes that encode intraflagellar transport (IFT) components and 74 ciliopathy loci to screen 92 unrelated individuals with BBS, irrespective of their known mutational burden. We identified 17 individuals with exon-disruptive CNVs (18.5%), including 13 different deletions in eight BBS genes (BBS1, BBS2, ARL6/BBS3, BBS4, BBS5, BBS7, BBS9, and NPHP1) and a deletion and a duplication in other ciliopathy-associated genes (ALMS1 and NPHP4, respectively). By contrast, we found a single heterozygous exon-disruptive event in a BBS-associated gene (BBS9) in 229 control subjects. Superimposing these data with resequencing revealed CNVs to (1) be sufficient to cause disease, (2) Mendelize heterozygous deleterious alleles, and (3) contribute oligogenic alleles by combining point mutations and exonic CNVs in multiple genes. Finally, we report a deletion and a splice site mutation in IFT74, inherited under a recessive paradigm, defining a candidate BBS locus. Our data suggest that CNVs contribute pathogenic alleles to a substantial fraction of BBS-affected individuals and highlight how either deletions or point mutations in discrete splice isoforms can induce hypomorphic mutations in genes otherwise intolerant to deleterious variation. Our data also suggest that CNV analyses and resequencing studies unbiased for previous mutational burden is necessary to delineate the complexity of disease architecture.

Nonrecurrent PMP22-RAI1 contiguous gene deletions arise from replication-based mechanisms and result in Smith-Magenis syndrome with evident peripheral neuropathy

Hereditary neuropathy with liability to pressure palsies (HNPP) and Smith-Magenis syndrome (SMS) are genomic disorders associated with deletion copy number variants involving chromosome 17p12 and 17p11.2, respectively. Nonallelic homologous recombination (NAHR)-mediated recurrent deletions are responsible for the majority of HNPP and SMS cases; the rearrangement products encompass the key dosage-sensitive genes PMP22 and RAI1, respectively, and result in haploinsufficiency for these genes. Less frequently, nonrecurrent genomic rearrangements occur at this locus. Contiguous gene duplications encompassing both PMP22 and RAI1, i.e., PMP22-RAI1 duplications, have been investigated, and replication-based mechanisms rather than NAHR have been proposed for these rearrangements. In the current study, we report molecular and clinical characterizations of six subjects with the reciprocal phenomenon of deletions spanning both genes, i.e., PMP22-RAI1 deletions. Molecular studies utilizing high-resolution array comparative genomic hybridization and breakpoint junction sequencing identified mutational signatures that were suggestive of replication-based mechanisms. Systematic clinical studies revealed features consistent with SMS, including features of intellectual disability, speech and gross motor delays, behavioral problems and ocular abnormalities. Five out of six subjects presented clinical signs and/or objective electrophysiologic studies of peripheral neuropathy. Clinical profiling may improve the clinical management of this unique group of subjects, as the peripheral neuropathy can be more severe or of earlier onset as compared to SMS patients having the common recurrent deletion. Moreover, the current study, in combination with the previous report of PMP22-RAI1 duplications, contributes to the understanding of rare complex phenotypes involving multiple dosage-sensitive genes from a genetic mechanistic standpoint.

High Frequency of Pathogenic Rearrangements in SPG11 and Extensive Contribution of Mutational Hotspots and Founder Alleles

Biallelic loss-of-function mutations in SPG11 cause a wide spectrum of recessively inherited, neurodegenerative disorders including hereditary spastic paraplegia (HSP), amyotrophic lateral sclerosis, and Charcot-Marie-Tooth disease. By comprehensive screening of three large cohorts of HSP index patients, we identified 83 alleles with "small" mutations and 13 alleles that carry large genomic rearrangements. Including relevant data from previous studies, we estimate that copy number variants (CNVs) account for ∼19% of pathogenic SPG11 alleles. The breakpoints for all novel and some previously reported CNVs were determined by long-range PCR and sequencing. This revealed several Alu-associated recombination hotspots. We also found evidence for additional mutational mechanisms, including for a two-step event in which an Alu retrotransposition preceded the actual rearrangement. Apparently independent samples with identical breakpoints were analyzed by microsatellite PCRs. The resulting haplotypes suggested the existence of two rearrangement founder alleles. Our findings widen the spectra of mutations and mutational mechanisms in SPG11, underscore the pivotal role played by Alus, and are of high diagnostic relevance for a wide spectrum of clinical phenotypes including the most frequent form of recessive HSP.

Heavy Metal Exposure Influences Double Strand Break DNA Repair Outcomes

Heavy metals such as cadmium, arsenic and nickel are classified as carcinogens. Although the precise mechanism of carcinogenesis is undefined, heavy metal exposure can contribute to genetic damage by inducing double strand breaks (DSBs) as well as inhibiting critical proteins from different DNA repair pathways. Here we take advantage of two previously published culture assay systems developed to address mechanistic aspects of DNA repair to evaluate the effects of heavy metal exposures on competing DNA repair outcomes. Our results demonstrate that exposure to heavy metals significantly alters how cells repair double strand breaks. The effects observed are both specific to the particular metal and dose dependent. Low doses of NiCl2 favored resolution of DSBs through homologous recombination (HR) and single strand annealing (SSA), which were inhibited by higher NiCl2 doses. In contrast, cells exposed to arsenic trioxide preferentially repaired using the "error prone" non-homologous end joining (alt-NHEJ) while inhibiting repair by HR. In addition, we determined that low doses of nickel and cadmium contributed to an increase in mutagenic recombination-mediated by Alu elements, the most numerous family of repetitive elements in humans. Sequence verification confirmed that the majority of the genetic deletions were the result of Alu-mediated non-allelic recombination events that predominantly arose from repair by SSA. All heavy metals showed a shift in the outcomes of alt-NHEJ repair with a significant increase of non-templated sequence insertions at the DSB repair site. Our data suggest that exposure to heavy metals will alter the choice of DNA repair pathway changing the genetic outcome of DSBs repair.

Heavy Metal Exposure Influences Double Strand Break DNA Repair Outcomes

Heavy metals such as cadmium, arsenic and nickel are classified as carcinogens. Although the precise mechanism of carcinogenesis is undefined, heavy metal exposure can contribute to genetic damage by inducing double strand breaks (DSBs) as well as inhibiting critical proteins from different DNA repair pathways. Here we take advantage of two previously published culture assay systems developed to address mechanistic aspects of DNA repair to evaluate the effects of heavy metal exposures on competing DNA repair outcomes. Our results demonstrate that exposure to heavy metals significantly alters how cells repair double strand breaks. The effects observed are both specific to the particular metal and dose dependent. Low doses of NiCl₂ favored resolution of DSBs through homologous recombination (HR) and single strand annealing (SSA), which were inhibited by higher NiCl₂ doses. In contrast, cells exposed to arsenic trioxide preferentially repaired using the “error prone” non-homologous end joining (alt-NHEJ) while inhibiting repair by HR. In addition, we determined that low doses of nickel and cadmium contributed to an increase in mutagenic recombination-mediated by Alu elements, the most numerous family of repetitive elements in humans. Sequence verification confirmed that the majority of the genetic deletions were the result of Alu-mediated non-allelic recombination events that predominantly arose from repair by SSA. All heavy metals showed a shift in the outcomes of alt-NHEJ repair with a significant increase of non-templated sequence insertions at the DSB repair site. Our data suggest that exposure to heavy metals will alter the choice of DNA repair pathway changing the genetic outcome of DSBs repair.

A polymorphic Alu insertion that mediates distinct disease-associated deletions

Large deletions that are associated with insertions of Alu-derived sequence represent a rare, but potentially unique class of alterations. Whether they form by a one-step mechanism or by a primary insertion step followed by an independent secondary deletion step is not clear. We resolved two disease-associated SPAST deletions, which involve distinct exons by long range PCR. Alu-derived sequence was observed between the breakpoints in both cases. The intronic regions that represent the targets of potentially involved Alu retrotransposition events overlapped. Microsatellite- and SNP-based haplotyping indicated that both deletions originated on one and the same founder allele. Our data suggest that the deletions are best explained by two-step insertion-deletion scenarios for which a single Alu retrotransposition event represents the shared primary step. This Alu then mediated one of the deletions by non-homologous end joining and the other by non-allelic homologous recombination. Our findings thus strongly argue for temporal separation of insertion and deletion in Alu insertion-associated deletions. They also suggest that certain Alu integrations confer a general increase in local genomic instability, and that this explains why they are usually not detected during the probably short time that precedes the rearrangements they mediate.

Mechanisms underlying structural variant formation in genomic disorders

With the recent burst of technological developments in genomics, and the clinical implementation of genome-wide assays, our understanding of the molecular basis of genomic disorders, specifically the contribution of structural variation to disease burden, is evolving quickly. Ongoing studies have revealed a ubiquitous role for genome architecture in the formation of structural variants at a given locus, both in DNA recombination-based processes and in replication-based processes. These reports showcase the influence of repeat sequences on genomic stability and structural variant complexity and also highlight the tremendous plasticity and dynamic nature of our genome in evolution, health and disease susceptibility.

Contrasting mechanisms of de novo copy number mutagenesis suggest the existence of different classes of environmental copy number mutagens

While gene copy number variations (CNVs) are abundant in the human genome, and often are associated with disease consequences, the mutagenic pathways and environmental exposures that cause these large structural mutations are understudied relative to conventional nucleotide substitutions in DNA. The members of the environmental mutagenesis community are currently seeking to remedy this deficiency, and there is a renewed interest in the development of mutagenicity assays to identify and characterize compounds that may induce de novo CNVs in humans. To achieve this goal, it is critically important to acknowledge that CNVs exist in two very distinct classes: nonrecurrent and recurrent CNVs. The goal of this commentary is to emphasize the deep contrasts that exist between the proposed pathways that lead to these two mutation classes. Nonrecurrent de novo CNVs originate primarily in mitotic cells through replication-dependent DNA repair pathways that involve microhomologies (<10 bp), and are detected at higher frequency in children of older fathers. In contrast, recurrent de novo CNVs are most often formed in meiotic cells through homologous recombination between nonallelic large low-copy repeats (>10,000 bp), without an associated paternal age effect. Given the biological differences between the two CNV classes, it is our belief that nonrecurrent and recurrent CN mutagens will probably differ substantially in their modes of action. Therefore, each CNV class may require their own uniquely designed assays, so that we as a field may succeed in uncovering the broadest possible spectrum of environmental CN mutagens.

Alu elements and DNA double-strand break repair

Alu elements represent one of the most common sources of homology and homeology in the human genome. Homeologous recombination between Alu elements represents a major form of genetic instability leading to deletions and duplications. Although these types of events have been studied extensively through genomic sequencing to assess the impact of Alu elements on disease mutations and genome evolution, the overall abundance of Alu elements in the genome often makes it difficult to assess the relevance of the Alu elements to specific recombination events. We recently reported a powerful new reporter gene system that allows the assessment of various cis and trans factors on the contribution of Alu elements to various forms of genetic instability. This allowed a quantitative measurement of the influence of mismatches on Alu elements and instability. It also confirmed that homeologous Alu elements are able to stimulate non-homologous end joining events in their vicinity. This appears to be dependent on portions of the mismatch repair pathway. We are now in a position to begin to unravel the complex influences of Alu density, mismatch and location with alterations of DNA repair processes in various tissues and tumors.

Human Structural Variation: Mechanisms of Chromosome Rearrangements

Chromosome structural variation (SV) is a normal part of variation in the human genome, but some classes of SV can cause neurodevelopmental disorders. Analysis of the DNA sequence at SV breakpoints can reveal mutational mechanisms and risk factors for chromosome rearrangement. Large-scale SV breakpoint studies have become possible recently owing to advances in next-generation sequencing (NGS) including whole-genome sequencing (WGS). These findings have shed light on complex forms of SV such as triplications, inverted duplications, insertional translocations, and chromothripsis. Sequence-level breakpoint data resolve SV structure and determine how genes are disrupted, fused, and/or misregulated by breakpoints. Recent improvements in breakpoint sequencing have also revealed non-allelic homologous recombination (NAHR) between paralogous long interspersed nuclear element (LINE) or human endogenous retrovirus (HERV) repeats as a cause of deletions, duplications, and translocations. This review covers the genomic organization of simple and complex constitutional SVs, as well as the molecular mechanisms of their formation.

The role of combined SNV and CNV burden in patients with distal symmetric polyneuropathy

PURPOSE

Charcot-Marie-Tooth (CMT) disease is a heterogeneous group of genetic disorders of the peripheral nervous system. Copy-number variants (CNVs) contribute significantly to CMT, as duplication of PMP22 underlies the majority of CMT1 cases. We hypothesized that CNVs and/or single-nucleotide variants (SNVs) might exist in patients with CMT with an unknown molecular genetic etiology.

METHODS

Two hundred patients with CMT, negative for both SNV mutations in several CMT genes and for CNVs involving PMP22, were screened for CNVs by high-resolution oligonucleotide array comparative genomic hybridization. Whole-exome sequencing was conducted on individuals with rare, potentially pathogenic CNVs.

RESULTS

Putatively causative CNVs were identified in five subjects (~2.5%); four of the five map to known neuropathy genes. Breakpoint sequencing revealed Alu-Alu-mediated junctions as a predominant contributor. Exome sequencing identified MFN2 SNVs in two of the individuals.

CONCLUSION

Neuropathy-associated CNV outside of the PMP22 locus is rare in CMT. Nevertheless, there is potential clinical utility in testing for CNVs and exome sequencing in CMT cases negative for the CMT1A duplication. These findings suggest that complex phenotypes including neuropathy can potentially be caused by a combination of SNVs and CNVs affecting more than one disease-associated locus and contributing to a mutational burden.Genet Med 18 5, 443-451.

DNA REPAIR. Mus81 and converging forks limit the mutagenicity of replication fork breakage

Most spontaneous DNA double-strand breaks (DSBs) result from replication-fork breakage. Break-induced replication (BIR), a genome rearrangement-prone repair mechanism that requires the Pol32/POLD3 subunit of eukaryotic DNA Polδ, was proposed to repair broken forks, but how genome destabilization is avoided was unknown. We show that broken fork repair initially uses error-prone Pol32-dependent synthesis, but that mutagenic synthesis is limited to within a few kilobases from the break by Mus81 endonuclease and a converging fork. Mus81 suppresses template switches between both homologous sequences and diverged human Alu repetitive elements, highlighting its importance for stability of highly repetitive genomes. We propose that lack of a timely converging fork or Mus81 may propel genome instability observed in cancer.

Decoding NF1 Intragenic Copy-Number Variations

Genomic rearrangements can cause both Mendelian and complex disorders. Currently, several major mechanisms causing genomic rearrangements, such as non-allelic homologous recombination (NAHR), non-homologous end joining (NHEJ), fork stalling and template switching (FoSTeS), and microhomology-mediated break-induced replication (MMBIR), have been proposed. However, to what extent these mechanisms contribute to gene-specific pathogenic copy-number variations (CNVs) remains understudied. Furthermore, few studies have resolved these pathogenic alterations at the nucleotide-level. Accordingly, our aim was to explore which mechanisms contribute to a large, unique set of locus-specific non-recurrent genomic rearrangements causing the genetic neurocutaneous disorder neurofibromatosis type 1 (NF1). Through breakpoint-spanning PCR as well as array comparative genomic hybridization, we have identified the breakpoints in 85 unrelated individuals carrying an NF1 intragenic CNV. Furthermore, we characterized the likely rearrangement mechanisms of these 85 CNVs, along with those of two additional previously published NF1 intragenic CNVs. Unlike the most typical recurrent rearrangements mediated by flanking low-copy repeats (LCRs), NF1 intragenic rearrangements vary in size, location, and rearrangement mechanisms. We propose the DNA-replication-based mechanisms comprising both FoSTeS and/or MMBIR and serial replication stalling to be the predominant mechanisms leading to NF1 intragenic CNVs. In addition to the loop within a 197-bp palindrome located in intron 40, four Alu elements located in introns 1, 2, 3, and 50 were also identified as intragenic-rearrangement hotspots within NF1.

Unbalanced translocations arise from diverse mutational mechanisms including chromothripsis

Unbalanced translocations are a relatively common type of copy number variation and a major contributor to neurodevelopmental disorders. We analyzed the breakpoints of 57 unique unbalanced translocations to investigate the mechanisms of how they form. Fifty-one are simple unbalanced translocations between two different chromosome ends, and six rearrangements have more than three breakpoints involving two to five chromosomes. Sequencing 37 breakpoint junctions revealed that simple translocations have between 0 and 4 base pairs (bp) of microhomology (n = 26), short inserted sequences (n = 8), or paralogous repeats (n = 3) at the junctions, indicating that translocations do not arise primarily from nonallelic homologous recombination but instead form most often via nonhomologous end joining or microhomology-mediated break-induced replication. Three simple translocations fuse genes that are predicted to produce in-frame transcripts of SIRPG-WWOX, SMOC2-PROX1, and PIEZO2-MTA1, which may lead to gain of function. Three complex translocations have inversions, insertions, and multiple breakpoint junctions between only two chromosomes. Whole-genome sequencing and fluorescence in situ hybridization analysis of two de novo translocations revealed at least 18 and 33 breakpoints involving five different chromosomes. Breakpoint sequencing of one maternally inherited translocation involving four chromosomes uncovered multiple breakpoints with inversions and insertions. All of these breakpoint junctions had 0-4 bp of microhomology consistent with chromothripsis, and both de novo events occurred on paternal alleles. Together with other studies, these data suggest that germline chromothripsis arises in the paternal genome and may be transmitted maternally. Breakpoint sequencing of our large collection of chromosome rearrangements provides a comprehensive analysis of the molecular mechanisms behind translocation formation.

Hereditary spastic paraplegia SPG4: what is known and not known about the disease

Mutations in more than 70 distinct loci and more than 50 mutated gene products have been identified in patients with hereditary spastic paraplegias, a diverse group of neurological disorders characterized predominantly, but not exclusively, by progressive lower limb spasticity and weakness resulting from distal degeneration of corticospinal tract axons. Mutations in the SPAST (previously known as SPG4) gene that encodes the microtubule-severing protein called spastin, are the most common cause of the disease. The aetiology of the disease is poorly understood, but partial loss of microtubule-severing activity resulting from inactivating mutations in one SPAST allele is the most postulated explanation. Microtubule severing is important for regulating various aspects of the microtubule array, including microtubule number, length, and mobility. In addition, higher numbers of dynamic plus-ends of microtubules, resulting from microtubule-severing events, may play a role in endosomal tubulation and fission. Even so, there is growing evidence that decreased severing of microtubules does not fully explain HSP-SPG4. The presence of two translation initiation codons in SPAST allows synthesis of two spastin isoforms: a full-length isoform called M1 and a slightly shorter isoform called M87. M87 is more abundant in both neuronal and non-neuronal tissues. Studies on rodents suggest that M1 is only readily detected in adult spinal cord, which is where nerve degeneration mainly occurs in humans with HSP-SPG4. M1, due to its hydrophobic N-terminal domain not shared by M87, may insert into endoplasmic reticulum membrane, and together with reticulons, atlastin and REEP1, may play a role in the morphogenesis of this organelle. Some mutated spastins may act in dominant-negative fashion to lower microtubule-severing activity, but others have detrimental effects on neurons without further lowering microtubule severing. The observed adverse effects on microtubule dynamics, axonal transport, endoplasmic reticulum, and endosomal trafficking are likely caused not only by diminished severing of microtubules, but also by neurotoxicity of mutant spastin proteins, chiefly M1. Some large deletions in SPAST might also affect the function of adjacent genes, further complicating the aetiology of the disease.

Genomic and Epigenetic Complexity of the FOXF1 Locus in 16q24.1: Implications for Development and Disease

The FOXF1 (Forkhead box F1) gene, located on chromosome 16q24.1 encodes a member of the FOX family of transcription factors characterized by a distinct forkhead DNA binding domain. FOXF1 plays an important role in epithelium-mesenchyme signaling, as a downstream target of Sonic hedgehog pathway. Heterozygous point mutations and genomic deletions involving FOXF1 have been reported in newborns with a lethal lung developmental disorder, Alveolar Capillary Dysplasia with Misalignment of Pulmonary Veins (ACDMPV). In addition, genomic deletions upstream to FOXF1 identified in ACDMPV patients have revealed that FOXF1 expression is tightly regulated by distal tissue-specific enhancers. Interestingly, FOXF1 has been found to be incompletely paternally imprinted in human lungs; characterized genomic deletions arose de novo exclusively on maternal chromosome 16, with most of them being Alu-Alu mediated. Regulation of FOXF1 expression likely utilizes a combination of chromosomal looping, differential methylation of an upstream CpG island overlapping GLI transcription factor binding sites, and the function of lung-specific long non-coding RNAs (lncRNAs). FOXF1 knock-out mouse models demonstrated its critical role in mesoderm differentiation and in the development of pulmonary vasculature. Additionally, epigenetic inactivation of FOXF1 has been reported in breast and colorectal cancers, whereas overexpression of FOXF1 has been associated with a number of other human cancers, e.g. medulloblastoma and rhabdomyosarcoma. Constitutional duplications of FOXF1 have recently been reported in congenital intestinal malformations. Thus, understanding the genomic and epigenetic complexity at the FOXF1 locus will improve diagnosis, prognosis, and treatment of ACDMPV and other human disorders associated with FOXF1 alterations.

Germline PRKACA amplification causes variable phenotypes that may depend on the extent of the genomic defect: molecular mechanisms and clinical presentations

OBJECTIVE

We have recently reported five patients with bilateral adrenocortical hyperplasia (BAH) and Cushing's syndrome (CS) caused by constitutive activation of the catalytic subunit of protein kinase A (PRKACA). By doing new in-depth analysis of their cytogenetic abnormality, we attempted a better genotype-phenotype correlation of their PRKACA amplification.

DESIGN

This study is a case series.

METHODS

Molecular cytogenetic, genomic, clinical, and histopathological analyses were performed in five patients with CS.

RESULTS

Reinvestigation of the defects of previously described patients by state-of-the-art molecular cytogenetics showed complex genomic rearrangements in the chromosome 19p13.2p13.12 locus, resulting in copy number gains encompassing the entire PRKACA gene; three patients (one sporadic case and two related cases) were observed with gains consistent with duplications, while two sporadic patients were observed with gains consistent with triplications. Although all five patients presented with ACTH-independent CS, the three sporadic patients had micronodular BAH and underwent bilateral adrenalectomy in early childhood, whereas the two related patients, a mother and a son, presented with macronodular BAH as adults. In at least one patient, PRKACA triplication was associated with a more severe phenotype.

CONCLUSIONS

Constitutional chromosomal PRKACA gene amplification is a recently identified genetic defect associated with CS, a trait that may be inherited in an autosomal dominant manner or occur de novo. Genomic rearrangements can be complex and can result in different copy number states of dosage-sensitive genes, e.g., duplication and triplication. PRKACA amplification can lead to variable phenotypes clinically and pathologically, both micro- and macro-nodular BAH, the latter of which we speculate may depend on the extent of amplification.

Alu-mediated diverse and complex pathogenic copy-number variants within human chromosome 17 at p13.3

Alu repetitive elements are known to be major contributors to genome instability by generating Alu-mediated copy-number variants (CNVs). Most of the reported Alu-mediated CNVs are simple deletions and duplications, and the mechanism underlying Alu-Alu-mediated rearrangement has been attributed to non-allelic homologous recombination (NAHR). Chromosome 17 at the p13.3 genomic region lacks extensive low-copy repeat architecture; however, it is highly enriched for Alu repetitive elements, with a fraction of 30% of total sequence annotated in the human reference genome, compared with the 10% genome-wide and 18% on chromosome 17. We conducted mechanistic studies of the 17p13.3 CNVs by performing high-density oligonucleotide array comparative genomic hybridization, specifically interrogating the 17p13.3 region with ∼150 bp per probe density; CNV breakpoint junctions were mapped to nucleotide resolution by polymerase chain reaction and Sanger sequencing. Studied rearrangements include 5 interstitial deletions, 14 tandem duplications, 7 terminal deletions and 13 complex genomic rearrangements (CGRs). Within the 17p13.3 region, Alu-Alu-mediated rearrangements were identified in 80% of the interstitial deletions, 46% of the tandem duplications and 50% of the CGRs, indicating that this mechanism was a major contributor for formation of breakpoint junctions. Our studies suggest that Alu repetitive elements facilitate formation of non-recurrent CNVs, CGRs and other structural aberrations of chromosome 17 at p13.3. The common observation of Alu-mediated rearrangement in CGRs and breakpoint junction sequences analysis further demonstrates that this type of mechanism is unlikely attributed to NAHR, but rather may be due to a recombination-coupled DNA replicative repair process.

PacBio-LITS: a large-insert targeted sequencing method for characterization of human disease-associated chromosomal structural variations

Background

Generation of long (>5 Kb) DNA sequencing reads provides an approach for interrogation of complex regions in the human genome. Currently, large-insert whole genome sequencing (WGS) technologies from Pacific Biosciences (PacBio) enable analysis of chromosomal structural variations (SVs), but the cost to achieve the required sequence coverage across the entire human genome is high.

Results

We developed a method (termed PacBio-LITS) that combines oligonucleotide-based DNA target-capture enrichment technologies with PacBio large-insert library preparation to facilitate SV studies at specific chromosomal regions. PacBio-LITS provides deep sequence coverage at the specified sites at substantially reduced cost compared with PacBio WGS. The efficacy of PacBio-LITS is illustrated by delineating the breakpoint junctions of low copy repeat (LCR)-associated complex structural rearrangements on chr17p11.2 in patients diagnosed with Potocki–Lupski syndrome (PTLS; MIM#610883). We successfully identified previously determined breakpoint junctions in three PTLS cases, and also were able to discover novel junctions in repetitive sequences, including LCR-mediated breakpoints. The new information has enabled us to propose mechanisms for formation of these structural variants.

Conclusions

The new method leverages the cost efficiency of targeted capture-sequencing as well as the mappability and scaffolding capabilities of long sequencing reads generated by the PacBio platform. It is therefore suitable for studying complex SVs, especially those involving LCRs, inversions, and the generation of chimeric Alu elements at the breakpoints. Other genomic research applications, such as haplotype phasing and small insertion and deletion validation could also benefit from this technology.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1370-2) contains supplementary material, which is available to authorized users.

Intra-individual plasticity of the TAZ gene leading to different heritable mutations in siblings with Barth syndrome

Infantile-onset skeletal myopathy Barth syndrome (OMIM #302060) is caused by mutations in the X-linked TAZ gene and hence usually manifests itself only in hemizygous males. Confirmatory testing is provided by mutational analysis of the TAZ gene and/or by biochemical dosage of the monolysocardiolipin/tetralinoleoyl cardiolipin ratio. Heterozygous females do not usually display a clinical phenotype but may undergo molecular genetic prenatal diagnosis during pregnancy. We characterized two novel and non-identical TAZ gene rearrangements in the offspring of a single female carrier of Barth syndrome. The hg19chrX:g.153634427_153644361delinsKP_123427.1 TAZ gene rearrangement was identified in her affected son, whereas the NM_000116.3(TAZ)c.-72_109+51del TAZ gene deletion was identified in a male foetus during a subsequent pregnancy. The unaffected mother was surprisingly found to harbour both variants in addition to a wild-type TAZ allele. A combination of breakpoint junction sequencing, linkage analysis and assessment of allelic dosage revealed that the two variants had originated independently from an apparently unstable/mutable TAZ maternal allele albeit via different mutational mechanisms. We conclude that molecular prenatal diagnosis in Barth syndrome families with probands carrying TAZ gene rearrangements should include investigation of the entire coding region of the TAZ gene. The identification of the breakpoint junctions of such gross gene rearrangements is important to ensure accurate ascertainment of carriership with a view to providing appropriate genetic counselling.

The contribution of alu elements to mutagenic DNA double-strand break repair

Alu elements make up the largest family of human mobile elements, numbering 1.1 million copies and comprising 11% of the human genome. As a consequence of evolution and genetic drift, Alu elements of various sequence divergence exist throughout the human genome. Alu/Alu recombination has been shown to cause approximately 0.5% of new human genetic diseases and contribute to extensive genomic structural variation. To begin understanding the molecular mechanisms leading to these rearrangements in mammalian cells, we constructed Alu/Alu recombination reporter cell lines containing Alu elements ranging in sequence divergence from 0%-30% that allow detection of both Alu/Alu recombination and large non-homologous end joining (NHEJ) deletions that range from 1.0 to 1.9 kb in size. Introduction of as little as 0.7% sequence divergence between Alu elements resulted in a significant reduction in recombination, which indicates even small degrees of sequence divergence reduce the efficiency of homology-directed DNA double-strand break (DSB) repair. Further reduction in recombination was observed in a sequence divergence-dependent manner for diverged Alu/Alu recombination constructs with up to 10% sequence divergence. With greater levels of sequence divergence (15%-30%), we observed a significant increase in DSB repair due to a shift from Alu/Alu recombination to variable-length NHEJ which removes sequence between the two Alu elements. This increase in NHEJ deletions depends on the presence of Alu sequence homeology (similar but not identical sequences). Analysis of recombination products revealed that Alu/Alu recombination junctions occur more frequently in the first 100 bp of the Alu element within our reporter assay, just as they do in genomic Alu/Alu recombination events. This is the first extensive study characterizing the influence of Alu element sequence divergence on DNA repair, which will inform predictions regarding the effect of Alu element sequence divergence on both the rate and nature of DNA repair events.

DNA double-strand breaks (DSBs) are a highly mutagenic form of DNA damage that can be repaired through one of several pathways with varied degrees of sequence preservation. Faithful repair of DSBs often occurs through gene conversion in which a sister chromatid is used as a repair template. Unfaithful repair of DSBs can occur through non-allelic homologous or homeologous recombination, which leads to chromosomal abnormalities such as deletions, duplications, and translocations and has been shown to cause several human genetic diseases. Substrates for these homologous and homeologous events include Alu elements, which are approximately 300 bp elements that comprise ~11% of the human genome. We use a new reporter assay to show that repair of DSBs results in Alu-mediated deletions that resolve through several distinct repair pathways. Either single-strand annealing (SSA) repair or microhomology-mediated end joining occurs ‘in register’ between two Alu elements when Alu sequence divergence is low. However, with more diverged Alu elements, like those typically found in the human genome, repair of DSBs appears to use the Alu/Alu homeology to direct non-homologous end joining in the general vicinity of the Alu elements. Mutagenic NHEJ repair involving divergent Alu elements may represent a common repair event in primate genomes.

Complex genomic rearrangements at the PLP1 locus include triplication and quadruplication

Inverted repeats (IRs) can facilitate structural variation as crucibles of genomic rearrangement. Complex duplication-inverted triplication-duplication (DUP-TRP/INV-DUP) rearrangements that contain breakpoint junctions within IRs have been recently associated with both MECP2 duplication syndrome (MIM#300260) and Pelizaeus-Merzbacher disease (PMD, MIM#312080). We investigated 17 unrelated PMD subjects with copy number gains at the PLP1 locus including triplication and quadruplication of specific genomic intervals-16/17 were found to have a DUP-TRP/INV-DUP rearrangement product. An IR distal to PLP1 facilitates DUP-TRP/INV-DUP formation as well as an inversion structural variation found frequently amongst normal individuals. We show that a homology-or homeology-driven replicative mechanism of DNA repair can apparently mediate template switches within stretches of microhomology. Moreover, we provide evidence that quadruplication and potentially higher order amplification of a genomic interval can occur in a manner consistent with rolling circle amplification as predicted by the microhomology-mediated break induced replication (MMBIR) model.

Association of new deletion/duplication region at chromosome 1p21 with intellectual disability, severe speech deficit and autism spectrum disorder-like behavior: an all-in approach to solving the DPYD enigma

We describe an as yet unreported neocentric small supernumerary marker chromosome (sSMC) derived from chromosome 1p21.3p21.2. It was present in 80% of the lymphocytes in a male patient with intellectual disability, severe speech deficit, mild dysmorphic features, and hyperactivity with elements of autism spectrum disorder (ASD).

Several important neurodevelopmental genes are affected by the 3.56 Mb copy number gain of 1p21.3p21.2, which may be considered reciprocal in gene content to the recently recognized 1p21.3 microdeletion syndrome. Both 1p21.3 deletions and the presented duplication display overlapping symptoms, fitting the same disorder category. Contribution of coding and non-coding genes to the phenotype is discussed in the light of cellular and intercellular homeostasis disequilibrium. In line with this the presented 1p21.3p21.2 copy number gain correlated to 1p21.3 microdeletion syndrome verifies the hypothesis of a cumulative effect of the number of deregulated genes - homeostasis disequilibrium leading to overlapping phenotypes between microdeletion and microduplication syndromes.

Although miR-137 appears to be the major player in the 1p21.3p21.2 region, deregulation of the DPYD (dihydropyrimidine dehydrogenase) gene may potentially affect neighboring genes underlying the overlapping symptoms present in both the copy number loss and copy number gain of 1p21. Namely, the all-in approach revealed that DPYD is a complex gene whose expression is epigenetically regulated by long non-coding RNAs (lncRNAs) within the locus. Furthermore, the long interspersed nuclear element-1 (LINE-1) L1MC1 transposon inserted in DPYD intronic transcript 1 (DPYD-IT1) lncRNA with its parasites, TcMAR-Tigger5b and pair of Alu repeats appears to be the “weakest link” within the DPYD gene liable to break. Identification of the precise mechanism through which DPYD is epigenetically regulated, and underlying reasons why exactly the break (FRA1E) happens, will consequently pave the way toward preventing severe toxicity to the antineoplastic drug 5-fluorouracil (5-FU) and development of the causative therapy for the dihydropyrimidine dehydrogenase deficiency.

Next-generation sequencing of duplication CNVs reveals that most are tandem and some create fusion genes at breakpoints

Interpreting the genomic and phenotypic consequences of copy-number variation (CNV) is essential to understanding the etiology of genetic disorders. Whereas deletion CNVs lead obviously to haploinsufficiency, duplications might cause disease through triplosensitivity, gene disruption, or gene fusion at breakpoints. The mutational spectrum of duplications has been studied at certain loci, and in some cases these copy-number gains are complex chromosome rearrangements involving triplications and/or inversions. However, the organization of clinically relevant duplications throughout the genome has yet to be investigated on a large scale. Here we fine-mapped 184 germline duplications (14.7 kb-25.3 Mb; median 532 kb) ascertained from individuals referred for diagnostic cytogenetics testing. We performed next-generation sequencing (NGS) and whole-genome sequencing (WGS) to sequence 130 breakpoints from 112 subjects with 119 CNVs and found that most (83%) were tandem duplications in direct orientation. The remainder were triplications embedded within duplications (8.4%), adjacent duplications (4.2%), insertional translocations (2.5%), or other complex rearrangements (1.7%). Moreover, we predicted six in-frame fusion genes at sequenced duplication breakpoints; four gene fusions were formed by tandem duplications, one by two interconnected duplications, and one by duplication inserted at another locus. These unique fusion genes could be related to clinical phenotypes and warrant further study. Although most duplications are positioned head-to-tail adjacent to the original locus, those that are inverted, triplicated, or inserted can disrupt or fuse genes in a manner that might not be predicted by conventional copy-number assays. Therefore, interpreting the genetic consequences of duplication CNVs requires breakpoint-level analysis.

Genome-wide analyses of LINE-LINE-mediated nonallelic homologous recombination

Nonallelic homologous recombination (NAHR), occurring between low-copy repeats (LCRs) >10 kb in size and sharing >97% DNA sequence identity, is responsible for the majority of recurrent genomic rearrangements in the human genome. Recent studies have shown that transposable elements (TEs) can also mediate recurrent deletions and translocations, indicating the features of substrates that mediate NAHR may be significantly less stringent than previously believed. Using >4 kb length and >95% sequence identity criteria, we analyzed of the genome-wide distribution of long interspersed element (LINE) retrotransposon and their potential to mediate NAHR. We identified 17 005 directly oriented LINE pairs located <10 Mbp from each other as potential NAHR substrates, placing 82.8% of the human genome at risk of LINE-LINE-mediated instability. Cross-referencing these regions with CNVs in the Baylor College of Medicine clinical chromosomal microarray database of 36 285 patients, we identified 516 CNVs potentially mediated by LINEs. Using long-range PCR of five different genomic regions in a total of 44 patients, we confirmed that the CNV breakpoints in each patient map within the LINE elements. To additionally assess the scale of LINE-LINE/NAHR phenomenon in the human genome, we tested DNA samples from six healthy individuals on a custom aCGH microarray targeting LINE elements predicted to mediate CNVs and identified 25 LINE-LINE rearrangements. Our data indicate that LINE-LINE-mediated NAHR is widespread and under-recognized, and is an important mechanism of structural rearrangement contributing to human genomic variability.

Dosage changes of a segment at 17p13.1 lead to intellectual disability and microcephaly as a result of complex genetic interaction of multiple genes

The 17p13.1 microdeletion syndrome is a recently described genomic disorder with a core clinical phenotype of intellectual disability, poor to absent speech, dysmorphic features, and a constellation of more variable clinical features, most prominently microcephaly. We identified five subjects with copy-number variants (CNVs) on 17p13.1 for whom we performed detailed clinical and molecular studies. Breakpoint mapping and retrospective analysis of published cases refined the smallest region of overlap (SRO) for microcephaly to a genomic interval containing nine genes. Dissection of this phenotype in zebrafish embryos revealed a complex genetic architecture: dosage perturbation of four genes (ASGR1, ACADVL, DVL2, and GABARAP) impeded neurodevelopment and decreased dosage of the same loci caused a reduced mitotic index in vitro. Moreover, epistatic analyses in vivo showed that dosage perturbations of discrete gene pairings induce microcephaly. Taken together, these studies support a model in which concomitant dosage perturbation of multiple genes within the CNV drive the microcephaly and possibly other neurodevelopmental phenotypes associated with rearrangements in the 17p13.1 SRO.