Intrachromosomal serial replication slippage intransgives rise to diverse genomic rearrangements involving inversions

Dynamic Expansions and Retinal Expression of Spectrally Distinct Short-Wavelength Opsin Genes in Sea Snakes

The photopigment-encoding visual opsin genes that mediate color perception show great variation in copy number and adaptive function across vertebrates. An open question is how this variation has been shaped by the interaction of lineage-specific structural genomic architecture and ecological selection pressures. We contribute to this issue by investigating the expansion dynamics and expression of the duplicated Short-Wavelength-Sensitive-1 opsin (SWS1) in sea snakes (Elapidae). We generated one new genome, 45 resequencing datasets, 10 retinal transcriptomes, and 81 SWS1 exon sequences for sea snakes, and analyzed these alongside 16 existing genomes for sea snakes and their terrestrial relatives. Our analyses revealed multiple independent transitions in SWS1 copy number in the marine Hydrophis clade, with at least three lineages having multiple intact SWS1 genes: the previously studied Hydrophis cyanocinctus and at least two close relatives of this species; Hydrophis atriceps and Hydrophis fasciatus; and an individual Hydrophis curtus. In each lineage, gene copy divergence at a key spectral tuning site resulted in distinct UV and Violet/Blue-sensitive SWS1 subtypes. Both spectral variants were simultaneously expressed in the retinae of H. cyanocinctus and H. atriceps, providing the first evidence that these SWS1 expansions confer novel phenotypes. Finally, chromosome annotation for nine species revealed shared structural features in proximity to SWS1 regardless of copy number. If these features are associated with SWS1 duplication, expanded opsin complements could be more common in snakes than is currently recognized. Alternatively, selection pressures specific to aquatic environments could favor improved chromatic distinction in just some lineages.

Transposable element-mediated rearrangements are prevalent in human genomes

Transposable elements constitute about half of human genomes, and their role in generating human variation through retrotransposition is broadly studied and appreciated. Structural variants mediated by transposons, which we call transposable element-mediated rearrangements (TEMRs), are less well studied, and the mechanisms leading to their formation as well as their broader impact on human diversity are poorly understood. Here, we identify 493 unique TEMRs across the genomes of three individuals. While homology directed repair is the dominant driver of TEMRs, our sequence-resolved TEMR resource allows us to identify complex inversion breakpoints, triplications or other high copy number polymorphisms, and additional complexities. TEMRs are enriched in genic loci and can create potentially important risk alleles such as a deletion in TRIM65, a known cancer biomarker and therapeutic target. These findings expand our understanding of this important class of structural variation, the mechanisms responsible for their formation, and establish them as an important driver of human diversity.

Molecular basis of α-thalassemia

α-Thalassemia is an inherited, autosomal recessive, disorder characterized by a microcytic hypochromic anemia. It is one of the most common monogenic gene disorders in the world population. The clinical severity varies from almost asymptomatic, to mild microcytic hypochromic, and to a lethal hemolytic condition, called Hb Bart's Hydrops Foetalis Syndrome. The molecular basis are usually deletions and less frequently, point mutations affecting the expression of one or more of the duplicated α-genes. The clinical variation and increase in disease severity is directly related to the decreased expression of one, two, three or four copies of the α-globin genes. Deletions and point mutations in the α-globin genes and their regulatory elements have been studied extensively in carriers and patients and these studies have given insight into the α-globin genes are regulated. By looking at naturally occurring deletions and point mutations, our knowledge of globin-gene regulation and expression will continue to increase and will lead to new targets of therapy.

Short template switch events explain mutation clusters in the human genome

Resequencing efforts are uncovering the extent of genetic variation in humans and provide data to study the evolutionary processes shaping our genome. One recurring puzzle in both intra- and inter-species studies is the high frequency of complex mutations comprising multiple nearby base substitutions or insertion-deletions. We devised a generalized mutation model of template switching during replication that extends existing models of genome rearrangement and used this to study the role of template switch events in the origin of short mutation clusters. Applied to the human genome, our model detects thousands of template switch events during the evolution of human and chimp from their common ancestor and hundreds of events between two independently sequenced human genomes. Although many of these are consistent with a template switch mechanism previously proposed for bacteria, our model also identifies new types of mutations that create short inversions, some flanked by paired inverted repeats. The local template switch process can create numerous complex mutation patterns, including hairpin loop structures, and explains multinucleotide mutations and compensatory substitutions without invoking positive selection, speculative mechanisms, or implausible coincidence. Clustered sequence differences are challenging for current mapping and variant calling methods, and we show that many erroneous variant annotations exist in human reference data. Local template switch events may have been neglected as an explanation for complex mutations because of biases in commonly used analyses. Incorporation of our model into reference-based analysis pipelines and comparisons of de novo assembled genomes will lead to improved understanding of genome variation and evolution.

Large Intragenic Deletion in DSTYK Underlies Autosomal-Recessive Complicated Spastic Paraparesis, SPG23

SPG23 is an autosomal-recessive neurodegenerative subtype of lower limb spastic paraparesis with additional diffuse skin and hair dyspigmentation at birth followed by further patchy pigment loss during childhood. Previously, genome-wide linkage in an Arab-Israeli pedigree mapped the gene to an approximately 25 cM locus on chromosome 1q24-q32. By using whole-exome sequencing in a further Palestinian-Jordanian SPG23 pedigree, we identified a complex homozygous 4-kb deletion/20-bp insertion in DSTYK (dual serine-threonine and tyrosine protein kinase) in all four affected family members. DSTYK is located within the established linkage region and we also found the same mutation in the previously reported pedigree and another Israeli pedigree (total of ten affected individuals from three different families). The mutation removes the last two exons and part of the 3' UTR of DSTYK. Skin biopsies revealed reduced DSTYK protein levels along with focal loss of melanocytes. Ultrastructurally, swollen mitochondria and cytoplasmic vacuoles were also noted in remaining melanocytes and some keratinocytes and fibroblasts. Cultured keratinocytes and fibroblasts from an affected individual, as well as knockdown of Dstyk in mouse melanocytes, keratinocytes, and fibroblasts, were associated with increased cell death after ultraviolet irradiation. Keratinocytes from an affected individual showed loss of kinase activity upon stimulation with fibroblast growth factor. Previously, dominant mutations in DSTYK were implicated in congenital urological developmental disorders, but our study identifies different phenotypic consequences for a recurrent autosomal-recessive deletion mutation in revealing the genetic basis of SPG23.

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.

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.

Complex Multiple-Nucleotide Substitution Mutations Causing Human Inherited Disease Reveal Novel Insights into the Action of Translesion Synthesis DNA Polymerases

Translesion synthesis (TLS) DNA polymerases allow the bypass of unrepaired lesions during DNA replication. Based upon mutational signatures of a subtype of multiple-nucleotide substitution (MNS) mutations causing human inherited disease, we have recently postulated two properties of TLS DNA polymerases in DNA repair, namely, the generation of neo-microhomologies potentiating strand-misalignment, and additional microlesions within the templated inserts when recruited to stalled replication forks. To provide further support for this postulate, we analyzed the mutational signatures of a new and complex subtype of pathogenic MNS mutation. Several mutations containing long templated inserts (8-19 bp) that are highly informative with regard to their underlying mutational mechanisms, harbor imprints of TLS DNA polymerase action. Dissecting the mechanism underlying the generation of the 19-bp insert implicated repeated participation of TLS DNA polymerases in the conversion of a damaged base into a complex MNS lesion through a process of successive template switching and bypass repair.

Palindrome-mediated and replication-dependent pathogenic structural rearrangements within the NF1 gene

Palindromic sequences can form hairpin structures or cruciform extrusions, which render them susceptible to genomic rearrangements. A 197-bp long palindromic AT-rich repeat (PATRR17) is located within intron 40 of the neurofibromatosis type 1 (NF1) gene (17q11.2). Through comprehensive NF1 analysis, we identified six unrelated patients with a rearrangement involving intron 40 (five deletions and one reciprocal translocation t(14;17)(q32;q11.2)). We hypothesized that PATRR17 may be involved in these rearrangements thereby causing NF1. Breakpoint cloning revealed that PATRR17 was indeed involved in all of the rearrangements. As microhomology was present at all breakpoint junctions of the deletions identified, and PATRR17 partner breakpoints were located within 7.1 kb upstream of PATRR17, fork stalling and template switching/microhomology-mediated break-induced replication was the most likely rearrangement mechanism. For the reciprocal translocation case, a 51 bp insertion at the translocation breakpoints mapped to a short sequence within PATRR17, proximal to the breakpoint, suggesting a multiple stalling and rereplication process, in contrast to previous studies indicating a purely replication-independent mechanism for PATRR-mediated translocations. In conclusion, we show evidence that PATRR17 is a hotspot for pathogenic intragenic deletions within the NF1 gene and suggest a novel replication-dependent mechanism for PATRR-mediated translocation.

Replicative mechanisms for CNV formation are error prone

We investigated 67 breakpoint junctions of gene copy number gains (CNVs) in 31 unrelated subjects. We observed a strikingly high frequency of small deletions and insertions (29%) apparently originating from polymerase-slippage events, in addition to frameshifts and point mutations in homonucleotide runs (13%), at or flanking the breakpoint junctions of complex CNVs. These simple nucleotide variants (SNV) were generated concomitantly with the de novo complex genomic rearrangement (CGR) event. Our findings implicate a low fidelity error-prone DNA polymerase in synthesis associated with DNA repair mechanisms that leads to a local increase in point mutation burden associated with human CGR.

Ionizing radiation and genetic risks. XVII. Formation mechanisms underlying naturally occurring DNA deletions in the human genome and their potential relevance for bridging the gap between induced DNA double-strand breaks and deletions in irradiated germ cells

While much is known about radiation-induced DNA double-strand breaks (DSBs) and their repair, the question of how deletions of different sizes arise as a result of the processing of DSBs by the cell's repair systems has not been fully answered. In order to bridge this gap between DSBs and deletions, we critically reviewed published data on mechanisms pertaining to: (a) repair of DNA DSBs (from basic studies in this area); (b) formation of naturally occurring structural variation (SV) - especially of deletions - in the human genome (from genomic studies) and (c) radiation-induced mutations and structural chromosomal aberrations in mammalian somatic cells (from radiation mutagenesis and radiation cytogenetic studies). The specific aim was to assess the relative importance of the postulated mechanisms in generating deletions in the human genome and examine whether empirical data on radiation-induced deletions in mouse germ cells are consistent with predictions of these mechanisms. The mechanisms include (a) NHEJ, a DSB repair process that does not require any homology and which functions in all stages of the cell cycle (and is of particular relevance in G0/G1); (b) MMEJ, also a DSB repair process but which requires microhomology and which presumably functions in all cell cycle stages; (c) NAHR, a recombination-based DSB repair mechanism which operates in prophase I of meiosis in germ cells; (d) MMBIR, a microhomology-mediated, replication-based mechanism which operates in the S phase of the cell cycle, and (e) strand slippage during replication (involved in the origin of small insertions and deletions (INDELs). Our analysis permits the inference that, between them, these five mechanisms can explain nearly all naturally occurring deletions of different sizes identified in the human genome, NAHR and MMBIR being potentially more versatile in this regard. With respect to radiation-induced deletions, the basic studies suggest that those arising as a result of the operation of NHEJ/MMEJ processes, as currently formulated, are expected to be relatively small. However, data on induced mutations in mouse spermatogonial stem cells (irradiation in G0/G1 phase of the cell cycle and DSB repair presumed to be via NHEJ predominantly) show that most are associated with deletions of different sizes, some in the megabase range. There is thus a 'discrepancy' between what the basic studies suggest and the empirical observations in mutagenesis studies. This discrepancy, however, is only an apparent but not a real one. It can be resolved by considering the issue of deletions in the broader context of and in conjunction with the organization of chromatin in chromosomes and nuclear architecture, the conceptual framework for which already exists in studies carried out during the past fifteen years or so. In this paper, we specifically hypothesize that repair of DSBs induced in chromatin loops may offer a basis to explain the induction of deletions of different sizes and suggest an approach to test the hypothesis. We emphasize that the bridging of the gap between induced DSB and resulting deletions of different sizes is critical for current efforts in computational modeling of genetic risks.

Patterns and mutational signatures of tandem base substitutions causing human inherited disease

Tandem base substitutions (TBSs) are multiple mutations that comprise two or more contiguous nucleotide substitutions without any net gain or loss of bases. They have recently become recognized as a distinct category of human genomic variant. However, their role in causing human inherited disease so far has not been studied methodically. Here, using data from the Human Gene Mutation Database (http://www.hgmd.org), we identified 477 events to be TBSs (doublets, 448; triplets, 16; and quadruplets to octuplets, 13). A comprehensive sequence pattern and context analysis implied the likely fundamental importance of translesion synthesis (TLS) DNA polymerases in generating these diverse TBSs but revealed that TLS polymerases may operate differently in generating TBSs of ≤ 3 bases (bypass of endogenous DNA lesions) than those of ≥ 4 bases (serial replication slippage). Moreover, GC was found to be the most frequently affected dinucleotide with GC/GC>AA/TT being the most frequent double TBS. Comparison with cancer genome mutational spectra allowed us to conclude that human germline TBSs arise predominantly through the action of endogenous mechanisms of mutagenesis rather than through exposure to exogenous mutagens. Finally, the rates of double and triple TBSs were estimated to be 0.2-1.2 × 10(-10) and 0.8-4.8 × 10(-12) per base per generation, respectively.

Genome instability: does genetic diversity amplification drive tumorigenesis?

Recent data show that catastrophic events during one cell cycle can cause massive genome damage producing viable clones with unstable genomes. This is in contrast with the traditional view that tumorigenesis requires a long-term process in which mutations gradually accumulate over decades. These sudden events are likely to result in a large increase in genomic diversity within a relatively short time, providing the opportunity for selective advantages to be gained by a subset of cells within a population. This genetic diversity amplification, arising from a single aberrant cell cycle, may drive a population conversion from benign to malignant. However, there is likely a period of relative genome stability during the clonal expansion of tumors - this may provide an opportunity for therapeutic intervention, especially if mechanisms that limit tolerance of aneuploidy are exploited.

Transient hypermutability, chromothripsis and replication-based mechanisms in the generation of concurrent clustered mutations

Clustered mutations may be broadly defined as the presence of two or more mutations within a spatially localized genomic region on a single chromosome. Known instances vary in terms of both the number and type of the component mutations, ranging from two closely spaced point mutations to tens or even hundreds of genomic rearrangements. Although clustered mutations can represent the observable net result of independent lesions sequentially acquired over multiple cell cycles, they can also be generated in a simultaneous or quasi-simultaneous manner within a single cell cycle. This review focuses on those mechanisms known to underlie the latter type. Both gene conversion and transient hypermutability are capable of generating closely spaced multiple mutations. However, a recently described phenomenon in human cancer cells, known as 'chromothripsis', has provided convincing evidence that tens to hundreds of genomic rearrangements can sometimes be generated simultaneously via a single catastrophic event. The distinctive genomic features observed in the derivative chromosomes, together with the highly characteristic junction sequences, point to non-homologous end joining (NHEJ) as being the likely underlying mutational mechanism. By contrast, replication-based mechanisms such as microhomology-mediated break-induced replication (MMBIR) which involves serial replication slippage or serial template switching probably account for those complex genomic rearrangements that comprise multiple duplications and/or triplications.

Aberrant firing of replication origins potentially explains intragenic nonrecurrent rearrangements within genes, including the human DMD gene

Non-allelic homologous recombination (NAHR), non-homologous end joining (NHEJ), and microhomology-mediated replication-dependent recombination (MMRDR) have all been put forward as mechanisms to explain DNA rearrangements associated with genomic disorders. However, many nonrecurrent rearrangements in humans remain unexplained. To further investigate the mutation mechanisms of these copy number variations (CNVs), we performed breakpoint mapping analysis for 62 clinical cases with intragenic deletions in the human DMD gene (50 cases) and other known disease-causing genes (one PCCB, one IVD, one DBT, three PAH, one STK11, one HEXB, three DBT, one HRPT1, and one EMD cases). While repetitive elements were found in only four individual cases, three involving DMD and one HEXB gene, microhomologies (2-10 bp) were observed at breakpoint junctions in 56% and insertions ranging from 1 to 48 bp were seen in 16 of the total 62 cases. Among these insertions, we observed evidence for tandem repetitions of short segments (5-20 bp) of reference sequence proximal to the breakpoints in six individual DMD cases (six repeats in one, four repeats in three, two repeats in one, and one repeat in one case), strongly indicating attempts by the replication machinery to surpass the stalled replication fork. We provide evidence of a novel template slippage event during replication rescue. With a deeper insight into the complex process of replication and its rescue during origin failure, brought forward by recent studies, we propose a hypothesis based on aberrant firing of replication origins to explain intragenic nonrecurrent rearrangements within genes, including the DMD gene.

First description of ABCB4 gene deletions in familial low phospholipid-associated cholelithiasis and oral contraceptives-induced cholestasis

The wide clinical spectrum of the ABCB4 gene (ATP-binding cassette subfamily B member 4) deficiency syndromes in humans includes low phospholipid-associated cholelithiasis (LPAC), intrahepatic cholestasis of pregnancy (ICP), oral contraceptives-induced cholestasis (CIC), and progressive familial intrahepatic cholestasis type 3 (PFIC3). No ABCB4 mutations are found in a significant proportion of patients with these syndromes. In the present study, 102 unrelated adult patients with LPAC (43 patients) or CIC/ICP (59 patients) were screened for ABCB4 mutations using DNA sequencing. Heterozygous ABCB4 point or short insertion/deletion mutations were found in 37% (16/43) of the LPAC patients and in 27% (16/59) of the ICP/CIC patients. High-resolution gene dosage methodologies were used in the 70 negative patients. Here, we describe for the first time ABCB4 partial or complete heterozygous deletions in 7% (3/43) of the LPAC patients, and in 2% (1/59) of the ICP/CIC patients. Our observations urge to systematically test patients with LPAC, ICP/CIC, and also children with PFIC3 for the presence of ABCB4 deletions using molecular tools allowing detection of gross rearrangements. In clinical practice, a comprehensive ABCB4 alteration-screening algorithm will permit the use of ABCB4 genotyping to confirm the diagnosis of LPAC or ICP/CIC, and allow familial testing. An early diagnosis of these biliary diseases may be beneficial because of the preventive effect of ursodeoxycholic acid on biliary complications. Further comparative studies of patients with well-characterized genotypes (including deletions) and phenotypes will help determine whether ABCB4 mutation types influence clinical outcomes.

On the sequence-directed nature of human gene mutation: the role of genomic architecture and the local DNA sequence environment in mediating gene mutations underlying human inherited disease

Different types of human gene mutation may vary in size, from structural variants (SVs) to single base-pair substitutions, but what they all have in common is that their nature, size and location are often determined either by specific characteristics of the local DNA sequence environment or by higher order features of the genomic architecture. The human genome is now recognized to contain "pervasive architectural flaws" in that certain DNA sequences are inherently mutation prone by virtue of their base composition, sequence repetitivity and/or epigenetic modification. Here, we explore how the nature, location and frequency of different types of mutation causing inherited disease are shaped in large part, and often in remarkably predictable ways, by the local DNA sequence environment. The mutability of a given gene or genomic region may also be influenced indirectly by a variety of noncanonical (non-B) secondary structures whose formation is facilitated by the underlying DNA sequence. Since these non-B DNA structures can interfere with subsequent DNA replication and repair and may serve to increase mutation frequencies in generalized fashion (i.e., both in the context of subtle mutations and SVs), they have the potential to serve as a unifying concept in studies of mutational mechanisms underlying human inherited disease.

Diverse mutational mechanisms cause pathogenic subtelomeric rearrangements

Chromosome rearrangements are a significant cause of intellectual disability and birth defects. Subtelomeric rearrangements, including deletions, duplications and translocations of chromosome ends, were first discovered over 40 years ago and are now recognized as being responsible for several genetic syndromes. Unlike the deletions and duplications that cause some genomic disorders, subtelomeric rearrangements do not typically have recurrent breakpoints and involve many different chromosome ends. To capture the molecular mechanisms responsible for this heterogeneous class of chromosome abnormality, we coupled high-resolution array CGH with breakpoint junction sequencing of a diverse collection of subtelomeric rearrangements. We analyzed 102 breakpoints corresponding to 78 rearrangements involving 28 chromosome ends. Sequencing 21 breakpoint junctions revealed signatures of non-homologous end-joining, non-allelic homologous recombination between interspersed repeats and DNA replication processes. Thus, subtelomeric rearrangements arise from diverse mutational mechanisms. In addition, we find hotspots of subtelomeric breakage at the end of chromosomes 9q and 22q; these sites may correspond to genomic regions that are particularly susceptible to double-strand breaks. Finally, fine-mapping the smallest subtelomeric rearrangements has narrowed the critical regions for some chromosomal disorders.

FoSTeS, MMBIR and NAHR at the human proximal Xp region and the mechanisms of human Xq isochromosome formation

The recently described DNA replication-based mechanisms of fork stalling and template switching (FoSTeS) and microhomology-mediated break-induced replication (MMBIR) were previously shown to catalyze complex exonic, genic and genomic rearrangements. By analyzing a large number of isochromosomes of the long arm of chromosome X (i(Xq)), using whole-genome tiling path array comparative genomic hybridization (aCGH), ultra-high resolution targeted aCGH and sequencing, we provide evidence that the FoSTeS and MMBIR mechanisms can generate large-scale gross chromosomal rearrangements leading to the deletion and duplication of entire chromosome arms, thus suggesting an important role for DNA replication-based mechanisms in both the development of genomic disorders and cancer. Furthermore, we elucidate the mechanisms of dicentric i(Xq) (idic(Xq)) formation and show that most idic(Xq) chromosomes result from non-allelic homologous recombination between palindromic low copy repeats and highly homologous palindromic LINE elements. We also show that non-recurrent-breakpoint idic(Xq) chromosomes have microhomology-associated breakpoint junctions and are likely catalyzed by microhomology-mediated replication-dependent recombination mechanisms such as FoSTeS and MMBIR. Finally, we stress the role of the proximal Xp region as a chromosomal rearrangement hotspot.

A complex double deletion in LMNA underlies progressive cardiac conduction disease, atrial arrhythmias, and sudden death

BACKGROUND

Cardiac conduction disease is a clinically and genetically heterogeneous disorder characterized by defects in electrical impulse generation and conduction and is associated with sudden cardiac death.

METHODS AND RESULTS

We studied a 4-generation family with autosomal dominant progressive cardiac conduction disease, including atrioventricular conduction block and sinus bradycardia, atrial arrhythmias, and sudden death. Genome-wide linkage analysis mapped the disease locus to chromosome 1p22-q21. Multiplex ligation-dependent probe amplification analysis of the LMNA gene, which encodes the nuclear-envelope protein lamin A/C, revealed a novel gene rearrangement involving a 24-bp inversion flanked by a 3.8-kb deletion upstream and a 7.8-kb deletion downstream. The presence of short inverted sequence homologies at the breakpoint junctions suggested a mutational event involving serial replication slippage in trans during DNA replication.

CONCLUSIONS

We identified for the first time a complex LMNA gene rearrangement involving a double deletion in a 4-generation Dutch family with progressive conduction system disease. Our findings underscore the fact that if conventional polymerase chain reaction-based direct sequencing approaches for LMNA analysis are negative in suggestive pedigrees, mutation detection techniques capable of detecting gross genomic lesions involving deletions and insertions should be considered.

Complete ascertainment of intragenic copy number mutations (CNMs) in the CFTR gene and its implications for CNM formation at other autosomal loci

Over the last 20 years since the discovery of the cystic fibrosis transmembrane conductance regulator (CFTR) gene, more than 1,600 different putatively pathological CFTR mutations have been identified. Until now, however, copy number mutations (CNMs) involving the CFTR gene have not been methodically analyzed, resulting almost certainly in the underascertainment of CFTR gene duplications compared with deletions. Here, high-resolution array comparative genomic hybridization (averaging one interrogating probe every 95 bp) was used to analyze the entire length of the CFTR gene (189 kb) in 233 cystic fibrosis chromosomes lacking conventional mutations. We succeeded in identifying five duplication CNMs that would otherwise have been refractory to analysis. Based upon findings from this and other studies, we propose that deletion and duplication CNMs in the human autosomal genome are likely to be generated in the proportion of approximately 2-3:1. We further postulate that intragenic gene duplication CNMs in other disease loci may have been routinely underascertained. Finally, our analysis of +/-20 bp flanking each of the 40 CFTR breakpoints characterized at the DNA sequence level provide support for the emerging concept that non-B DNA conformations in combination with specific sequence motifs predispose to both recurring and nonrecurring genomic rearrangements.

Genomic rearrangements in inherited disease and cancer

Genomic rearrangements in inherited disease and cancer involve gross alterations of chromosomes or large chromosomal regions and can take the form of deletions, duplications, insertions, inversions or translocations. The characterization of a considerable number of rearrangement breakpoints has now been accomplished at the nucleotide sequence level, thereby providing an invaluable resource for the detailed study of the mutational mechanisms which underlie genomic recombination events. A better understanding of these mutational mechanisms is vital for improving the design of mutation detection strategies. At least five categories of mutational mechanism are known to give rise to genomic rearrangements: (i) homologous recombination including non-allelic homologous recombination (NAHR), gene conversion, single strand annealing (SSA) and break-induced replication (BIR), (ii) non-homologous end joining (NHEJ), (iii) microhomology-mediated replication-dependent recombination (MMRDR), (iv) long interspersed element-1 (LINE-1 or L1)-mediated retrotransposition and (v) telomere healing. Focussing on the first three of these general mechanisms, we compare and contrast their hallmark characteristics, and discuss the role of various local DNA sequence features (e.g. recombination-promoting motifs, repetitive sequences and sequences capable of non-B DNA formation) in mediating the recombination events that underlie gross genomic rearrangements. Finally, we explore how studies both at the level of the gene (using the neurofibromatosis type-1 gene as an example) and the whole genome (using data derived from cancer genome sequencing studies) are shaping our understanding of the impact of genomic rearrangements as a cause of human genetic disease.

Measurements of spontaneous rates of mutations in the recent past and the near future

The rate of spontaneous mutation in natural populations is a fundamental parameter for many evolutionary phenomena. Because the rate of mutation is generally low, most of what is currently known about mutation has been obtained through indirect, complex and imprecise methodological approaches. However, in the past few years genome-wide sequencing of closely related individuals has made it possible to estimate the rates of mutation directly at the level of the DNA, avoiding most of the problems associated with using indirect methods. Here, we review the methods used in the past with an emphasis on next generation sequencing, which may soon make the accurate measurement of spontaneous mutation rates a matter of routine.

Identification of uncommon recurrent Potocki-Lupski syndrome-associated duplications and the distribution of rearrangement types and mechanisms in PTLS

Nonallelic homologous recombination (NAHR) can mediate recurrent rearrangements in the human genome and cause genomic disorders. Smith-Magenis syndrome (SMS) and Potocki-Lupski syndrome (PTLS) are genomic disorders associated with a 3.7 Mb deletion and its reciprocal duplication in 17p11.2, respectively. In addition to these common recurrent rearrangements, an uncommon recurrent 5 Mb SMS-associated deletion has been identified. However, its reciprocal duplication predicted by the NAHR mechanism had not been identified. Here we report the molecular assays on 74 subjects with PTLS-associated duplications, 35 of whom are newly investigated. By both oligonucleotide-based comparative genomic hybridization and recombination hot spot analyses, we identified two cases of the predicted 5 Mb uncommon recurrent PTLS-associated duplication. Interestingly, the crossovers occur in proximity to a recently delineated allelic homologous recombination (AHR) hot spot-associated sequence motif, further documenting the common hot spot features shared between NAHR and AHR. An additional eight subjects with nonrecurrent PTLS duplications were identified. The smallest region of overlap (SRO) for all of the 74 PTLS duplications examined is narrowed to a 125 kb interval containing only RAI1, a gene recently further implicated in autism. Sequence complexities consistent with DNA replication-based mechanisms were identified in four of eight (50%) newly identified nonrecurrent PTLS duplications. Our findings of the uncommon recurrent PTLS-associated duplication at a relative prevalence reflecting the de novo mutation rate and the distribution of 17p11.2 duplication types in PTLS reveal insights into both the contributions of new mutations and the different underlying mechanisms that generate genomic rearrangements causing genomic disorders.

Inversion variants in the human genome: role in disease and genome architecture

Significant advances have been made over the past 5 years in mapping and characterizing structural variation in the human genome. Despite this progress, our understanding of inversion variants is still very restricted. While unbalanced variants such as copy number variations can be mapped using array-based approaches, strategies for characterization of inversion variants have been limited and underdeveloped. Traditional cytogenetic approaches have long been able to identify microscopic inversion events, but discovery of submicroscopic events has remained elusive and largely ignored. With the advent of paired-end sequencing approaches, it is now possible to map inversions across the human genome. Based on the paired-end sequencing studies published to date, it is now feasible to make a first map of inversions across the human genome and to use this map to explore the characteristics and distribution of this form of variation. The current map of inversions indicates that many remain to be identified, especially in the smaller size ranges. This review provides an overview of the current knowledge about human inversions and their contribution to human phenotypes. Further characterization of inversions should be considered as an important step towards a deeper understanding of human variation and genome dynamics.

Closely spaced multiple mutations as potential signatures of transient hypermutability in human genes

Data from diverse organisms suggests that transient hypermutability is a general mutational mechanism with the potential to generate multiple synchronous mutations, a phenomenon probably best exemplified by closely spaced multiple mutations (CSMMs). Here we have attempted to extend the concept of transient hypermutability from somatic cells to the germline, using human inherited disease-causing multiple mutations as a model system. Employing stringent criteria for data inclusion, we have retrospectively identified numerous potential examples of pathogenic CSMMs that exhibit marked similarities to the CSMMs reported in other systems. These examples include (1) eight multiple mutations, each comprising three or more components within a sequence tract of <100 bp; (2) three possible instances of "mutation showers"; and (3) numerous highly informative "homocoordinate" mutations. Using the proportion of CpG substitution as a crude indicator of the relative likelihood of transient hypermutability, we present evidence to suggest that CSMMs comprising at least one pair of mutations separated by < or =100 bp may constitute signatures of transient hypermutability in human genes. Although this analysis extends the generality of the concept of transient hypermutability and provides new insights into what may be considered a novel mechanism of mutagenesis underlying human inherited disease, it has raised serious concerns regarding current practices in mutation screening.

Mechanisms of change in gene copy number

Deletions and duplications of chromosomal segments (copy number variants, CNVs) are a major source of variation between individual humans and are an underlying factor in human evolution and in many diseases, including mental illness, developmental disorders and cancer. CNVs form at a faster rate than other types of mutation, and seem to do so by similar mechanisms in bacteria, yeast and humans. Here we review current models of the mechanisms that cause copy number variation. Non-homologous end-joining mechanisms are well known, but recent models focus on perturbation of DNA replication and replication of non-contiguous DNA segments. For example, cellular stress might induce repair of broken replication forks to switch from high-fidelity homologous recombination to non-homologous repair, thus promoting copy number change.

Complex human chromosomal and genomic rearrangements

Copy number variation (CNV) is a major source of genetic variation among humans. In addition to existing as benign polymorphisms, CNVs can also convey clinical phenotypes, including genomic disorders, sporadic diseases and complex human traits. CNV results from genomic rearrangements that can represent simple deletion or duplication of a genomic segment, or be more complex. Complex chromosomal rearrangements (CCRs) have been known for some time but their mechanisms have remained elusive. Recent technology advances and high-resolution human genome analyses have revealed that complex genomic rearrangements can account for a large fraction of non-recurrent rearrangements at a given locus. Various mechanisms, most of which are DNA-replication-based, for example fork stalling and template switching (FoSTeS) and microhomology-mediated break-induced replication (MMBIR), have been proposed for generating such complex genomic rearrangements and are probably responsible for CCR.

Complex rearrangements in patients with duplications of MECP2 can occur by fork stalling and template switching

Duplication at the Xq28 band including the MECP2 gene is one of the most common genomic rearrangements identified in neurodevelopmentally delayed males. Such duplications are non-recurrent and can be generated by a non-homologous end joining (NHEJ) mechanism. We investigated the potential mechanisms for MECP2 duplication and examined whether genomic architectural features may play a role in their origin using a custom designed 4-Mb tiling-path oligonucleotide array CGH assay. Each of the 30 patients analyzed showed a unique duplication varying in size from approximately 250 kb to approximately 2.6 Mb. Interestingly, in 77% of these non-recurrent duplications, the distal breakpoints grouped within a 215 kb genomic interval, located 47 kb telomeric to the MECP2 gene. The genomic architecture of this region contains both direct and inverted low-copy repeat (LCR) sequences; this same region undergoes polymorphic structural variation in the general population. Array CGH revealed complex rearrangements in eight patients; in six patients the duplication contained an embedded triplicated segment, and in the other two, stretches of non-duplicated sequences occurred within the duplicated region. Breakpoint junction sequencing was achieved in four duplications and identified an inversion in one patient, demonstrating further complexity. We propose that the presence of LCRs in the vicinity of the MECP2 gene may generate an unstable DNA structure that can induce DNA strand lesions, such as a collapsed fork, and facilitate a Fork Stalling and Template Switching event producing the complex rearrangements involving MECP2.

A microhomology-mediated break-induced replication model for the origin of human copy number variation

Chromosome structural changes with nonrecurrent endpoints associated with genomic disorders offer windows into the mechanism of origin of copy number variation (CNV). A recent report of nonrecurrent duplications associated with Pelizaeus-Merzbacher disease identified three distinctive characteristics. First, the majority of events can be seen to be complex, showing discontinuous duplications mixed with deletions, inverted duplications, and triplications. Second, junctions at endpoints show microhomology of 2–5 base pairs (bp). Third, endpoints occur near pre-existing low copy repeats (LCRs). Using these observations and evidence from DNA repair in other organisms, we derive a model of microhomology-mediated break-induced replication (MMBIR) for the origin of CNV and, ultimately, of LCRs. We propose that breakage of replication forks in stressed cells that are deficient in homologous recombination induces an aberrant repair process with features of break-induced replication (BIR). Under these circumstances, single-strand 3′ tails from broken replication forks will anneal with microhomology on any single-stranded DNA nearby, priming low-processivity polymerization with multiple template switches generating complex rearrangements, and eventual re-establishment of processive replication.

Mechanisms for human genomic rearrangements

Genomic rearrangements describe gross DNA changes of the size ranging from a couple of hundred base pairs, the size of an average exon, to megabases (Mb). When greater than 3 to 5 Mb, such changes are usually visible microscopically by chromosome studies. Human diseases that result from genomic rearrangements have been called genomic disorders. Three major mechanisms have been proposed for genomic rearrangements in the human genome. Non-allelic homologous recombination (NAHR) is mostly mediated by low-copy repeats (LCRs) with recombination hotspots, gene conversion and apparent minimal efficient processing segments. NAHR accounts for most of the recurrent rearrangements: those that share a common size, show clustering of breakpoints, and recur in multiple individuals. Non-recurrent rearrangements are of different sizes in each patient, but may share a smallest region of overlap whose change in copy number may result in shared clinical features among different patients. LCRs do not mediate, but may stimulate non-recurrent events. Some rare NAHRs can also be mediated by highly homologous repetitive sequences (for example, Alu, LINE); these NAHRs account for some of the non-recurrent rearrangements. Other non-recurrent rearrangements can be explained by non-homologous end-joining (NHEJ) and the Fork Stalling and Template Switching (FoSTeS) models. These mechanisms occur both in germ cells, where the rearrangements can be associated with genomic disorders, and in somatic cells in which such genomic rearrangements can cause disorders such as cancer. NAHR, NHEJ and FoSTeS probably account for the majority of genomic rearrangements in our genome and the frequency distribution of the three at a given locus may partially reflect the genomic architecture in proximity to that locus. We provide a review of the current understanding of these three models.

The molecular basis of color vision in colorful fish: four long wave-sensitive (LWS) opsins in guppies (Poecilia reticulata) are defined by amino acid substitutions at key functional sites

BACKGROUND

Comparisons of functionally important changes at the molecular level in model systems have identified key adaptations driving isolation and speciation. In cichlids, for example, long wavelength-sensitive (LWS) opsins appear to play a role in mate choice and male color variation within and among species. To test the hypothesis that the evolution of elaborate coloration in male guppies (Poecilia reticulata) is also associated with opsin gene diversity, we sequenced long wavelength-sensitive (LWS) opsin genes in six species of the family Poeciliidae.

RESULTS

Sequences of four LWS opsin genes were amplified from the guppy genome and from mRNA isolated from adult guppy eyes. Variation in expression was quantified using qPCR. Three of the four genes encode opsins predicted to be most sensitive to different wavelengths of light because they vary at key amino acid positions. This family of LWS opsin genes was produced by a diversity of duplication events. One, an intronless gene, was produced prior to the divergence of families Fundulidae and Poeciliidae. Between-gene PCR and DNA sequencing show that two of the guppy LWS opsins are linked in an inverted orientation. This inverted tandem duplication event occurred near the base of the poeciliid tree in the common ancestor of Poecilia and Xiphophorus. The fourth sequence has been uncovered only in the genus Poecilia. In the guppies surveyed here, this sequence is a hybrid, with the 5' end most similar to one of the tandem duplicates and the 3' end identical to the other.

CONCLUSION

Enhanced wavelength discrimination, a possible consequence of opsin gene duplication and divergence, might have been an evolutionary prerequisite for color-based sexual selection and have led to the extraordinary coloration now observed in male guppies and in many other poeciliids.

A large deletion in the human alpha-globin cluster caused by a replication error is associated with an unexpectedly mild phenotype

We have characterized a newly identified 16.6 kb deletion which removes a significant proportion of the human alpha-globin cluster including the psizeta1, alpha(D), psialpha1 and alpha2-globin genes but leaves the duplicated alpha1 gene intact. This complicated rearrangement results from a combination of slippage and strand switching at sites of microhomology during replication. Functional analysis shows that expression of the remaining alpha1 gene is increased, rather than down-regulated by this deletion. This could be related to its proximity to the remote upstream alpha-globin regulatory elements or reduced competition for these elements in the absence of the dominant alpha2-globin gene. The finding of a very mild phenotype associated with such an extensive deletion in the alpha-globin cluster implies that much of the DNA removed by the deletion is likely to be functionally unimportant. These findings suggest that other than the upstream regulatory elements and promoter proximal elements there are unlikely to be additional positive cis-acting sequences in the alpha-globin cluster.

Somatic microindels in human cancer: the insertions are highly error-prone and derive from nearby but not adjacent sense and antisense templates

Somatic microindels (microdeletions with microinsertions) have been studied in normal mouse tissues using the Big Blue lacI transgenic mutation detection system. Here we analyze microindels in human cancers using an endogenous and transcribed gene, the TP53 gene. Microindel frequency, the enhancement of 1-2 microindels and other features are generally similar to that observed in the non-transcribed lacI gene in normal mouse tissues. The current larger sample of somatic microindels reveals recurroids: mutations in which deletions are identical and the co-localized insertion is similar. The data reveal that the inserted sequences derive from nearby but not adjacent sequences in contrast to the slippage that characterizes the great majority of pure microinsertions. The microindel inserted sequences derive from a template on the sense or antisense strand with similar frequency. The estimated error rate of the insertion process of 13% per bp is by far the largest reported in vivo, with the possible exception of somatic hypermutation in the immunoglobulin gene. The data constrain possible mechanisms of microindels and raise the question of whether microindels are 'scars' from the bypass of large DNA adducts by a translesional polymerase, e.g. the 'Tarzan model' presented herein.

Double complex mutations involving F8 and FUNDC2 caused by distinct break-induced replication

Genomic rearrangements are a well-recognized cause of genetic disease and can be formed by a variety of mechanisms. We report a complex rearrangement causing severe hemophilia A, identified and further characterized using a range of PCR-based methods, and confirmed using array-comparative genomic hybridization (array-CGH). This rearrangement consists of a 15.5-kb deletion/16-bp insertion located 0.6 kb from a 28.1-kb deletion/263-kb insertion at Xq28 and is one of the most complex rearrangements described at a DNA sequence level. We propose that the rearrangement was generated by distinct but linked cellular responses to double strand breakage, namely break-induced replication (BIR) and a novel model of break-induced serial replication slippage (SRS). The copy number of several genes is affected by this rearrangement, with deletion of part of the Factor VIII gene (F8, causing hemophilia A) and the FUNDC2 gene, and duplication of the TMEM185A, HSFX1, MAGEA9, and MAGEA11 genes. As the patient exhibits no clinically detectable phenotype other than hemophilia A, it appears that the biological effects of the other genes involved are not dosage-dependent. This investigation has provided novel insights into processes of DNA repair including BIR and the first description of SRS during repair in a pathological context.

Characterisation of a large complex intragenic re-arrangement in the FVII gene (F7) avoiding misdiagnosis in inherited factor VII deficiency

Inherited factor VII (FVII) deficiency is a rare autosomal recessive bleeding disorder mostly caused by point mutations. Large genomic re-arrangements at F7 locus could account for a fraction of mutant alleles that remain unidentified after DNA sequencing, because they escape conventional polymerase chain reaction (PCR)-based techniques. We report the first systematic screening of F7 for large re-arrangements, by semi-quantitative multiplex PCR of fluorescent fragments targeting the 9 exons and the promoter region. A well-characterised cohort of 43 unrelated patients either apparently homozygous for a F7 point mutation or carrying at least one unidentified F7 mutant allele participated in this study. Two large F7 re-arrangements were identified in two FVII-deficient pedigrees, including a discontinuous deletion involving two distinct portions of F7 whose proximal and distal end junctions were characterised. A simple and efficient method for the routine detection of gross alterations of F7, which accounted for 2.3% of mutant alleles in our sample, is now available in inherited FVII deficiency. This test should complement conventional PCR-based techniques not only in unsolved cases, but also where inheritance pattern analysis is not achievable.

Dominantly Inherited beta-Thalassemia

Dominantly inherited beta-thalassemia (thal) or "inclusion body beta-thalassemias" are heterogeneous at the molecular level and are due to mutations at or near the beta-globin gene locus. Many of these involve mutations of exon 3 of the beta-globin gene. They include frameshifts, premature chain termination (nonsense) mutations, and complex rearrangements that lead to the synthesis of truncated or elongated and highly unstable beta-globin gene products. The resulting beta chain variants are very unstable, and in many cases, the products of the dominantly inherited beta-thal are not detectable. Hematological and clinical observations made in several families with comparable forms of beta-thal and with certain highly unstable hemoglobin (Hb) variants, have indicated a striking overlap; many subjects with detectable unstable Hb variants and with a dominant type of beta-thal without a detectable abnormal Hb, have similar phenotypes. Here, a review of dominantly inherited beta-thal is given, and new examples of hyperunstable Hbs (Hb Stara Zagora and Hb Jambol) are presented. The first example is a hyperunstable variant named Hb Stara Zagora that was found in a 2-year-old Bulgarian boy. The abnormal Hb is associated with severe hemolytic anemia as a consequence of its hyper instability. The anemia was noticed at the age of 2 months and since then he has been on a regular monthly blood transfusion regimen. Hemoglobin analysis revealed no abnormalities, except the presence of inclusion bodies. Sequencing of the beta-globin gene revealed a heterozygosity for a 6 bp deletion (-TGGCTA) at codons 137 (the second and third bp), 138 and 139 (the first bp), thus forming a new codon at position 139 (GAT). This event eliminates three amino acids (Val-Ala-Asn) and introduces a new residue (Asp). It creates a new restriction site for HphI. The parents and his twin brother had no history of hemolysis. The paternity of the child was confirmed by DNA analysis. The second example is a new hyperunstable variant named Hb Jambol, found as a de novo mutation in a 2-year-old Bulgarian girl with severe hemolytic anemia. The mutation was detected through RNA/DNA analysis. It represents a complex genomic rearrangement involving an insertion of 23 nucleotides (nts) after IVS-II-535, a deletion of 310 nts extending from IVS-II-550 to the first nt of codon 108, and an insertion of 28 nts at the deletion junctions (derived from inverted sequence between nts +3707 and +3734 3' to the beta-globin gene termination codon). At the protein level, this mutation leads to a deletion of four amino acid residues (Leu-Leu-Glu-Asn) at positions 105, 106, 107 and 108, and an insertion of nine residues (Val-Pro-Ser-Val-Thr-Leu-Phe-Phe-Asp) at the same location, creating an abnormal elongated beta chain of 151 amino acid residues. The parents had no history of hemolysis. The paternity of the child was confirmed by DNA analysis.

Detection of two Alu insertions in the CFTR gene

BACKGROUND

LINE-1 (long interspersed element-1) or L1-mediated retrotransposition is a potent force in human genome evolution and an occasional cause of human genetic disease. Since the first report of two de novo L1 insertions in the F8 gene causing hemophilia A, more than 50 L1-mediated retrotranspositional events have been identified as causing human genetic disease. However, a significant bias has generally militated against the detection of these pathological events at autosomal loci. Based upon this and other observations, we surmised that some previously unresolved cystic fibrosis chromosomes might carry hitherto undetected L1-mediated retrotranspositional insertions at the CFTR locus. This study represents an attempt to identify such mutational events.

METHODS

100 previously unresolved cystic fibrosis chromosomes were carefully reanalyzed using quantitative high-performance liquid chromatography (QHPLC).

RESULTS

Two simple Alu insertions were identified in the CFTR gene, within exons 16 and 17b respectively.

CONCLUSIONS

Our findings have not only revealed a previously unknown mutational mechanism responsible for cystic fibrosis but also represent an important addition to the already diverse spectrum of known CFTR gene mutations. Experience with the CFTR gene suggests that pathological L1-mediated retrotranspositional events may also have been overlooked at other gene loci and should always be considered in cases that appear to be refractory to analysis.

A complex rearrangement in GBE1 causes both perinatal hypoglycemic collapse and late-juvenile-onset neuromuscular degeneration in glycogen storage disease type IV of Norwegian forest cats

Deficiency of glycogen branching enzyme (GBE) activity causes glycogen storage disease type IV (GSD IV), an autosomal recessive error of metabolism. Abnormal glycogen accumulates in myocytes, hepatocytes, and neurons, causing variably progressive, benign to lethal organ dysfunctions. A naturally occurring orthologue of human GSD IV was described previously in Norwegian forest cats (NFC). Here, we report that while most affected kittens die at or soon after birth, presumably due to hypoglycemia, survivors of the perinatal period appear clinically normal until onset of progressive neuromuscular degeneration at 5 months of age. Molecular investigation of affected cats revealed abnormally spliced GBE1 mRNA products and lack of GBE cross-reactive material in liver and muscle. Affected cats are homozygous for a complex rearrangement of genomic DNA in GBE1, constituted by a 334 bp insertion at the site of a 6.2 kb deletion that extends from intron 11 to intron 12 (g. IVS11+1552_IVS12-1339 del6.2kb ins334 bp), removing exon 12. An allele-specific, PCR-based test demonstrates that the rearrangement segregates with the disease in the GSD IV kindred and is not found in unrelated normal cats. Screening of 402 privately owned NFC revealed 58 carriers and 4 affected cats. The molecular characterization of feline GSD IV will enhance further studies of GSD IV pathophysiology and development of novel therapies in this unique animal model.

A complex rearrangement in the APC gene uncovered by multiplex ligation-dependent probe amplification

Germline mutations in the tumor suppressor gene APC are the underlying cause of familial adenomatous polyposis, an autosomal-dominant cancer predisposition syndrome of the colorectum. Here, we describe a complex pathogenic rearrangement in the APC gene that was detected during deletion screening and transmitted throughout at least three generations. The rearrangement consists of a deletion of 604 bp in intron 4 that impairs the binding site of the reverse primer for exon 4 and of an insertion of 119 bp in exon 4 that interferes with the binding site of the multiplex ligation-dependent probe amplification (MLPA) probes for exon 4. The insertion is composed of three duplicated sequences derived from exon 4, intron 3, and intron 4, all in inverse direction. By transcript analysis, we found that the mutation results in complete skipping of exon 4 and that it leads to a frameshift. The rearrangement would not have been identified had it occurred outside the MLPA hybridization site. Our findings demonstrate that part of the pathogenic mutations remain undetected by routine methods. Moreover, MLPA and RNA analysis alone would have led to an incorrect interpretation of a genomic deletion of exon 4.

Hb Jambol: a new hyperunstable hemoglobin causing severe hemolytic anemia

We describe a new hyperunstable beta-chain variant due to a complex genomic rearrangement. The abnormal hemoglobin (Hb) was found as a de novo mutation in a 2-year-old Bulgarian girl with severe hemolytic anemia. The mutation was detected through RNA/DNA analysis. It represents a complex genomic rearrangement involving an insertion of 23 nts after IVS-II-535 (derived by triplication of the 12-nts adjacent sequence and subsequent deletion of 1 nt), a deletion of 310 nts extending from IVS-II-550 to the first nt of Cd 108 and an insertion of 28 nts at the deletion junctions (derived from the inverted sequence between nts +3,707 and +3,734 3' to the beta-globin gene termination codon). At the protein level this mutation leads to a deletion of 4 amino acid residues (Leu-Leu-Glu-Asn) at positions 105-108 and an insertion of 9 residues (Val-Pro-Ser-Val-Thr-Leu-Phe-Phe-Asp) at the same location, creating an abnormal elongated beta-chain of 151 amino acid residues. This highly unstable variant was named 'Hb Jambol' after the geographic location in which the patient resides.