Mutation spectrum of Drosophila CNVs revealed by breakpoint sequencing

Background

The detailed study of breakpoints associated with copy number variants (CNVs) can elucidate the mutational mechanisms that generate them and the comparison of breakpoints across species can highlight differences in genomic architecture that may lead to lineage-specific differences in patterns of CNVs. Here, we provide a detailed analysis of Drosophila CNV breakpoints and contrast it with similar analyses recently carried out for the human genome.

Results

By applying split-read methods to a total of 10x coverage of 454 shotgun sequence across nine lines of D. melanogaster and by re-examining a previously published dataset of CNVs detected using tiling arrays, we identified the precise breakpoints of more than 600 insertions, deletions, and duplications. Contrasting these CNVs with those found in humans showed that in both taxa CNV breakpoints fall into three classes: blunt breakpoints; simple breakpoints associated with microhomology; and breakpoints with additional nucleotides inserted/deleted and no microhomology. In both taxa CNV breakpoints are enriched with non-B DNA sequence structures, which may impair DNA replication and/or repair. However, in contrast to human genomes, non-allelic homologous-recombination (NAHR) plays a negligible role in CNV formation in Drosophila. In flies, non-homologous repair mechanisms are responsible for simple, recurrent, and complex CNVs, including insertions of de novo sequence as large as 60 bp.

Conclusions

Humans and Drosophila differ considerably in the importance of homology-based mechanisms for the formation of CNVs, likely as a consequence of the differences in the abundance and distribution of both segmental duplications and transposable elements between the two genomes.     

Molecular Analysis of a Deletion Hotspot in the NRXN1 Region Reveals the Involvement of Short Inverted Repeats in Deletion CNVs

NRXN1 microdeletions occur at a relatively high frequency and confer increased risk for neurodevelopmental and neurobehavioral abnormalities. The mechanism that makes NRXN1 a deletion hotspot is unknown. Here, we identified deletions of the NRXN1 region in affected cohorts, confirming a strong association with the autism spectrum and other neurodevelopmental disorders. Interestingly, deletions in both affected and control individuals were clustered in the 5′ portion of NRXN1 and its immediate upstream region. To explore the mechanism of deletion, we mapped and analyzed the breakpoints of 32 deletions. At the deletion breakpoints, frequent microhomology (68.8%, 2–19 bp) suggested predominant mechanisms of DNA replication error and/or microhomology-mediated end-joining. Long terminal repeat (LTR) elements, unique non-B-DNA structures, and MEME-defined sequence motifs were significantly enriched, but Alu and LINE sequences were not. Importantly, small-size inverted repeats (minus self chains, minus sequence motifs, and partial complementary sequences) were significantly overrepresented in the vicinity of NRXN1 region deletion breakpoints, suggesting that, although they are not interrupted by the deletion process, such inverted repeats can predispose a region to genomic instability by mediating single-strand DNA looping via the annealing of partially reverse complementary strands and the promoting of DNA replication fork stalling and DNA replication error. Our observations highlight the potential importance of inverted repeats of variable sizes in generating a rearrangement hotspot in which individual breakpoints are not recurrent. Mechanisms that involve short inverted repeats in initiating deletion may also apply to other deletion hotspots in the human genome.     

A simple and universal ligation mediated fusion of genes based on hetero-staggered PCR for generating immunodominant chimeric proteins

A potential relationship between transposon-derived repeats (TDR) and human germline methylation is of biological importance since many genes are flanked by TDR and methylation could affect the expression of nearby genes. Furthermore, DNA methylation has been suggested as a global defense mechanism against genome instability threatened by TDR. We studied the correlation between the density of HapMap methyl-associated SNPs (mSNPs), a marker of germline methylation, and proportion of TDR.After correcting for confounding variables, we found a negative correlation between proportion of Alu repeats and mSNP density for 125-1000 kb windows. Similar results were found for the most active subgroup of repeats. In contrast, a negative correlation between proportion of L1 repeats and mSNP density was found only in the larger 1000 kb windows.Using methylation data on germ cells (sperm) from the Human Epigenome Project, we found a lower proportion of Alu repeats adjacent (3-15 kb) to hypermethylated amplicons. On the contrary, there was a higher proportion of L1 repeats in the 3-5 kb of sequence flanking hypermethylated amplicons but not in the 10-15 kb flanks.Our data indicate a differential response to the major repeat families and that DNA methylation is unlikely to be a uniform global defense system against all TDR. It appears to play a role for the L1 subgroup, with sequences adjacent to L1 repeats methylated in response to their proximity. In contrast, sequences adjacent to Alu repeats appear to be hypomethylated, arguing against a role of methylation in germline defense against those elements.     

Distribution,functional impact,and origin mechanisms of copy number variation in the barley genome

Background

There is growing evidence for the prevalence of copy number variation (CNV) and its role in phenotypic variation in many eukaryotic species. Here we use array comparative genomic hybridization to explore the extent of this type of structural variation in domesticated barley cultivars and wild barleys.

Results

A collection of 14 barley genotypes including eight cultivars and six wild barleys were used for comparative genomic hybridization. CNV affects 14.9% of all the sequences that were assessed. Higher levels of CNV diversity are present in the wild accessions relative to cultivated barley. CNVs are enriched near the ends of all chromosomes except 4H, which exhibits the lowest frequency of CNVs. CNV affects 9.5% of the coding sequences represented on the array and the genes affected by CNV are enriched for sequences annotated as disease-resistance proteins and protein kinases. Sequence-based comparisons of CNV between cultivars Barke and Morex provided evidence that DNA repair mechanisms of double-strand breaks via single-stranded annealing and synthesis-dependent strand annealing play an important role in the origin of CNV in barley.

Conclusions

We present the first catalog of CNVs in a diploid Triticeae species, which opens the door for future genome diversity research in a tribe that comprises the economically important cereal species wheat, barley, and rye. Our findings constitute a valuable resource for the identification of CNV affecting genes of agronomic importance. We also identify potential mechanisms that can generate variation in copy number in plant genomes.     

The Rtt109 histone acetyltransferase facilitates error-free replication to prevent CAG/CTG repeat contractions

Lysine 56 is acetylated on newly synthesized histone H3 in yeast, Drosophila and mammalian cells. All of the proteins involved in histone H3 lysine 56 (H3K56) acetylation are important for maintaining genome integrity. These include Rtt109, a histone acetyltransferase, responsible for acetylating H3K56, Asf1, a histone H3/H4 chaperone, and Hst3 and Hst4, histone deacetylases which remove the acetyl group from H3K56. Here we demonstrate a new role for Rtt109 and H3K56 acetylation in maintaining repetitive DNA sequences in Saccharomyces cerevisiae. We found that cells lacking RTT109 had a high level of CAG/CTG repeat contractions and a twofold increase in breakage at CAG/CTG repeats. In addition, repeat contractions were significantly increased in cells lacking ASF1 and in an hst3Δhst4Δ double mutant. Because the Rtt107/Rtt101 complex was previously shown to be recruited to stalled replication forks in an Rtt109-dependent manner, we tested whether this complex was involved. However, contractions in rtt109Δ cells were not due to an inability to recruit the Rtt107/Rtt101 complex to repeats, as absence of these proteins had no effect on repeat stability. On the other hand, Dnl4 and Rad51-dependent pathways did play a role in creating some of the repeat contractions in rtt109Δ cells. Our results show that H3K56 acetylation by Rtt109 is important for stabilizing DNA repeats, likely by facilitating proper nucleosome assembly at the replication fork to prevent DNA structure formation and subsequent slippage events or fork breakage.     

Origins and breakpoint analyses of copy number variations: up close and personal

Array-based methods have enabled the detection of many genomic gains and losses. These are stated as copy number variants (CNVs) and comprise up to 13% of the human genome. Based on their breakpoints and modes of formation CNVs are termed recurrent or nonrecurrent. Recurrent CNVs are flanked by low copy repeats and are of a fixed size. They arise as a result of misalignment during meiosis by a mechanism named nonallelic homologous recombination. Several of such recurrent CNVs have been linked to human diseases. Nonrecurrent CNVs, which are not flanked by low copy repeats, are of variable size and may arise via mechanisms like nonhomologous end joining and replication-based mechanisms described by the fork stalling and template switching and microhomology-mediated break-induced replication models. It is becoming clear that most disease-causing CNVs are nonrecurrent and generally arise via replication-based mechanisms. Furthermore, it is now appreciated that genomic features other than low copy repeats play a role in the formation of nonrecurrent CNVs. This review will discuss the different mechanisms of CNV formation and how high resolution analyses of CNV breakpoints have added to our knowledge of their precise structure.     

Evolutionary Conservation of Sequences Directing Chromosome Breakage and Rdna Palindrome Formation in Tetrahymenine Ciliates

Extensive, programmed chromosome breakage occurs during formation of the somatic macronucleus of ciliated protozoa. The cis-acting signal directing breakage has been most rigorously defined in Tetrahymena thermophila, where it consists of a 15-bp DNA sequence known as Cbs, for chromosome breakage sequence. We have identified sequences identical or nearly identical to the T. thermophila Cbs at sites of breakage flanking the germline micronuclear rDNA locus of six additional species of Tetrahymena as well as members of two related genera. Other general features of the breakage site are also conserved, but surprisingly, the orientation and number of copies of Cbs are not always conserved, suggesting the occurrence of germline rearrangement events over evolutionary time. At one end of the T. thermophila micronuclear rDNA locus, a pair of short inverted repeats adjacent to Cbs directs the formation of a giant palindromic molecule. We have examined the corresponding sequences from two other Tetrahymena species. We find the sequence to be partially conserved, as previously implied from analysis of macronuclear rDNA, but of variable length and organization.     

Direct and Inverted Repeats Elicit Genetic Instability by Both Exploiting and Eluding DNA Double-Strand Break Repair Systems in Mycobacteria

Repetitive DNA sequences with the potential to form alternative DNA conformations, such as slipped structures and cruciforms, can induce genetic instability by promoting replication errors and by serving as a substrate for DNA repair proteins, which may lead to DNA double-strand breaks (DSBs). However, the contribution of each of the DSB repair pathways, homologous recombination (HR), non-homologous end-joining (NHEJ) and single-strand annealing (SSA), to this sort of genetic instability is not fully understood. Herein, we assessed the genome-wide distribution of repetitive DNA sequences in the Mycobacterium smegmatis, Mycobacterium tuberculosis and Escherichia coli genomes, and determined the types and frequencies of genetic instability induced by direct and inverted repeats, both in the presence and in the absence of HR, NHEJ, and SSA. All three genomes are strongly enriched in direct repeats and modestly enriched in inverted repeats. When using chromosomally integrated constructs in M. smegmatis, direct repeats induced the perfect deletion of their intervening sequences ∼1,000-fold above background. Absence of HR further enhanced these perfect deletions, whereas absence of NHEJ or SSA had no influence, suggesting compromised replication fidelity. In contrast, inverted repeats induced perfect deletions only in the absence of SSA. Both direct and inverted repeats stimulated excision of the constructs from the attB integration sites independently of HR, NHEJ, or SSA. With episomal constructs, direct and inverted repeats triggered DNA instability by activating nucleolytic activity, and absence of the DSB repair pathways (in the order NHEJ>HR>SSA) exacerbated this instability. Thus, direct and inverted repeats may elicit genetic instability in mycobacteria by 1) directly interfering with replication fidelity, 2) stimulating the three main DSB repair pathways, and 3) enticing L5 site-specific recombination.     

A high-resolution map of synteny disruptions in gibbon and human genomes

Gibbons are part of the same superfamily (Hominoidea) as humans and great apes, but their karyotype has diverged faster from the common hominoid ancestor. At least 24 major chromosome rearrangements are required to convert the presumed ancestral karyotype of gibbons into that of the hominoid ancestor. Up to 28 additional rearrangements distinguish the various living species from the common gibbon ancestor. Using the northern white-cheeked gibbon (2n = 52) (Nomascus leucogenys leucogenys) as a model, we created a high-resolution map of the homologous regions between the gibbon and human. The positions of 100 synteny breakpoints relative to the assembled human genome were determined at a resolution of about 200 kb. Interestingly, 46% of the gibbon–human synteny breakpoints occur in regions that correspond to segmental duplications in the human lineage, indicating a common source of plasticity leading to a different outcome in the two species. Additionally, the full sequences of 11 gibbon BACs spanning evolutionary breakpoints reveal either segmental duplications or interspersed repeats at the exact breakpoint locations. No specific sequence element appears to be common among independent rearrangements. We speculate that the extraordinarily high level of rearrangements seen in gibbons may be due to factors that increase the incidence of chromosome breakage or fixation of the derivative chromosomes in a homozygous state.     

Molecular characterisation of the pericentric inversion that distinguishes human chromosome 5 from the homologous chimpanzee chromosome

Human and chimpanzee karyotypes differ by virtue of nine pericentric inversions that serve to distinguish human chromosomes 1, 4, 5, 9, 12, 15, 16, 17, and 18 from their chimpanzee orthologues. In this study, we have analysed the breakpoints of the pericentric inversion characteristic of chimpanzee chromosome 4, the homologue of human chromosome 5. Breakpoint-spanning BAC clones were identified from both the human and chimpanzee genomes by fluorescence in situ hybridisation, and the precise locations of the breakpoints were determined by sequence comparisons. In stark contrast to some other characterised evolutionary rearrangements in primates, this chimpanzee-specific inversion appears not to have been mediated by either gross segmental duplications or low-copy repeats, although micro-duplications were found adjacent to the breakpoints. However, alternating purine–pyrimidine (RY) tracts were detected at the breakpoints, and such sequences are known to adopt non-B DNA conformations that are capable of triggering DNA breakage and genomic rearrangements. Comparison of the breakpoint region of human chromosome 5q15 with the orthologous regions of the chicken, mouse, and rat genomes, revealed similar but non-identical syntenic disruptions in all three species. The clustering of evolutionary breakpoints within this chromosomal region, together with the presence of multiple pathological breakpoints in the vicinity of both 5p15 and 5q15, is consistent with the non-random model of chromosomal evolution and suggests that these regions may well possess intrinsic features that have served to mediate a variety of genomic rearrangements, including the pericentric inversion in chimpanzee chromosome 4.     

Heterogeneous duplications in patients with Pelizaeus-Merzbacher disease suggest a mechanism of coupled homologous and nonhomologous recombination

We describe genomic structures of 59 X-chromosome segmental duplications that include the proteolipid protein 1 gene (PLP1) in patients with Pelizaeus-Merzbacher disease. We provide the first report of 13 junction sequences, which gives insight into underlying mechanisms. Although proximal breakpoints were highly variable, distal breakpoints tended to cluster around low-copy repeats (LCRs) (50% of distal breakpoints), and each duplication event appeared to be unique (100 kb to 4.6 Mb in size). Sequence analysis of the junctions revealed no large homologous regions between proximal and distal breakpoints. Most junctions had microhomology of 1-6 bases, and one had a 2-base insertion. Boundaries between single-copy and duplicated DNA were identical to the reference genomic sequence in all patients investigated. Taken together, these data suggest that the tandem duplications are formed by a coupled homologous and nonhomologous recombination mechanism. We suggest repair of a double-stranded break (DSB) by one-sided homologous strand invasion of a sister chromatid, followed by DNA synthesis and nonhomologous end joining with the other end of the break. This is in contrast to other genomic disorders that have recurrent rearrangements formed by nonallelic homologous recombination between LCRs. Interspersed repetitive elements (Alu elements, long interspersed nuclear elements, and long terminal repeats) were found at 18 of the 26 breakpoint sequences studied. No specific motif that may predispose to DSBs was revealed, but single or alternating tracts of purines and pyrimidines that may cause secondary structures were common. Analysis of the 2-Mb region susceptible to duplications identified proximal-specific repeats and distal LCRs in addition to the previously reported ones, suggesting that the unique genomic architecture may have a role in nonrecurrent rearrangements by promoting instability.     

Molecular cytogenetic characterization of breakpoints involving pericentric inversions of human chromosome 9

Pericentric inversions involving the secondary constriction (qh) region of chromosome 9 are considered to be normal variants. The evolutionary mechanisms and conservation of these inversions via Mendelian fashion have been investigated since the advent of banding techniques. Routine cytogenetic techniques cannot provide the fine characterization necessary to determine the type of genetic material involved in these rearrangements. Therefore, the fluorescence in situ hybridization technique with the human centromere-specific alpha satellite and the beta satellite (D9Z5) and classical satellite (D9Z1) human DNA probes were used to identify the breakpoints of chromosome 9 pericentric inversions. Four unique types of pericentric inversions involving the 9qh region were observed, and the mechanism may be due to breakage and reunion at the proposed breakpoints. They are: type A inversions consist of breakpoints within the alpha and beta satellite DNA regions; type B consist of breakpoints within the beta satellite DNA region and band 9q13; type C involve breakage within the beta and classical satellite DNA regions, and type D have breakpoints within the alpha and classical satellite DNA regions. Obviously, reshuffling of satellite DNA sequences has occurred, which has given rise to a variety of heteromorphisms whose clinical significance remains obscure. Received: 21 December 1995 / Revised: 30 May 1996     

Analysis and Dynamics of the Chromosomal Complements of Wild Sparkling-Wine Yeast Strains

We isolated Saccharomyces cerevisiae yeast strains that are able to carry out the second fermentation of sparkling wine from spontaneously fermenting musts in El Penedès (Spain) by specifically designed selection protocols. All of them (26 strains) showed one of two very similar mitochondrial DNA (mtDNA) restriction patterns, whereas their karyotypes differed. These strains showed high rates of karyotype instability, which were dependent on both the medium and the strain, during vegetative growth. In all cases, the mtDNA restriction pattern was conserved in strains kept under the same conditions. Analysis of different repetitive sequences in their genomes suggested that ribosomal DNA repeats play an important role in the changes in size observed in chromosome XII, whereas SUC genes or Ty elements did not show amplification or transposition processes that could be related to rearrangements of the chromosomes showing these sequences. Karyotype changes also occurred in monosporidic diploid derivatives. We propose that these changes originated mainly from ectopic recombination between repeated sequences interspersed in the genome. None of the rearranged karyotypes provided a selective advantage strong enough to allow the strains to displace the parental strains. The nature and frequency of these changes suggest that they may play an important role in the establishment and maintenance of the genetic diversity observed in S. cerevisiae wild populations.     

A novel family of mobile genetic elements is limited to the germline genome in Tetrahymena thermophila

In the ciliated protozoan Tetrahymena thermophila, extensive DNA elimination is associated with differentiation of the somatic macronucleus from the germline micronucleus. This study describes the isolation and complete characterization of Tlr elements, a family of approximately 30 micronuclear DNA sequences that are efficiently eliminated from the developing macronucleus. The data indicate that Tlr elements are comprised of an ~22 kb internal region flanked by complex and variable termini. The Tlr internal region is highly conserved among family members and contains 15 open reading frames, some of which resemble genes encoded by transposons and viruses. The Tlr termini appear to be long inverted repeats consisting of (i) a variable region containing multiple direct repeats which differ in number and sequence from element to element and (ii) a conserved terminal 47 bp sequence. Taken together, these results suggest that Tlr elements comprise a novel family of mobile genetic elements that are confined to the Tetrahymena germline genome. Possible mechanisms of developmentally programmed Tlr elimination are discussed.     

Rearrangement structure-independent strategy of CNV breakpoint analysis

Rare copy number variations (CNVs) generated by human genomic rearrangements have been shown to play an important role in pathogenesis of human diseases and cancers. CNV breakpoint analysis can help define genomic location, genetic content and sequence structure of pathogenic CNVs. This process is vital to elucidate CNV mutational mechanism and etiology of CNV-associated disorders. However, it is technically challenging to map CNV breakpoints at base-pair level, especially in the genomic regions with sequence complexity. In this study, we developed a new method of capture and breakpoint approaching sequencing (CBAS) to efficiently obtain CNV breakpoint sequences. This strategy is independent of CNV structures and applicable to various CNV types. As was demonstrated in CNV-associated patients with neurological disorders, CBAS achieved fine mapping of breakpoint sequences for compound deletion, complex duplication, and translocation. Intriguingly, CBAS also revealed unexpected CNV complexity involving long-range DNA rearrangement. Our observations showed that CBAS is an efficient method for obtaining CNV breakpoint sequence and mapping insertional events as well. This method can facilitate the researches on CNV-associated human diseases and cancers. CBAS is also applicable to mapping the integration sites of retrovirus (such as HIV) and transgenes in model organisms.     

Comparison of simple sequence repeats in Staphylococcus strains using in-silico approach

Staphylococci are Gram-positive bacteria which play an important role in infectious disease and are major causes of communityacquired and hospital-acquired infections. Strains of Staphylococcus aureus are reported as genomically and phenotypically highly heterogeneous; hence in-silico based comparison of genomic data on simple sequence repeats may provide valuable information for understanding the pathogenicity and control measures. This study determined the distribution of a specific group of Simple Sequence Repeats (SSRs), in genome sequences of six Staphylococcus strains (Staphylococcus aureus COL, S.aureus MRSA252, S.aureus MSSA476, S.aureus Mu50, S.aureus MW2, S.aureus N315) and plasmid sequences of four Staphylococcus strains (Staphylococcus aureus COL pT181, Staphylococcus aureus MSSA pSAS, Staphylococcus aureus VRSAp, Staphylococcus aureus, Staphylococcus aureus pN315 DNA) downloaded from the GenBank database for identifying abundance, distribution and composition of SSRs. The data obtained in the present study shows that (i) a large number of tandem repeats are distributed throughout the genome and plasmid sequences. (ii) Number of mononucleotide SSRs decreased rapidly with increase in size of repeat unit. (iii) Total frequency of SSRs in plasmid regions is less than genomic regions. (iv) In all investigated strains, ratios of AT/TA repeats are dominating over GC/CG repeats in genomics as well as plasmid sequences, and (v) Dinucleotide combination of AT is dominated in all the six Staphylococcus genome sequences.     

Mechanism of Replicative DNA Polymerase Delta Pausing and a Potential Role for DNA Polymerase Kappa in Common Fragile Site Replication

Common fragile sites (CFSs) are hot spots of chromosomal breakage, and CFS breakage models involve perturbations of DNA replication. Here, we analyzed the contribution of specific repetitive DNA sequence elements within CFSs to the inhibition of DNA synthesis by replicative and specialized DNA polymerases (Pols). The efficiency of in vitro DNA synthesis was quantitated using templates corresponding to regions within FRA16D and FRA3B harboring AT-rich microsatellite and quasi-palindrome (QP) sequences. QPs were predicted to form stems of ~ 75–100% self-homology, separated by 3–9 bases of intervening sequences. Analysis of DNA synthesis progression by human Pol δ demonstrated significant synthesis perturbation both at [A]_n and [TA]_n repeats in a length-dependent manner and at short (< 40 base pairs) QP sequences. DNA synthesis by the Y-family polymerase κ was significantly more efficient than Pol δ through both types of repetitive elements. Using DNA trap experiments, we show that Pol δ pauses within CFS sequences are sites of enzyme dissociation, and dissociation was observed in the presence of RFC-loaded PCNA. We propose that enrichment of microsatellite and QP elements at CFS regions contributes to fragility by perturbing replication through multiple mechanisms, including replicative Pol pausing and dissociation. Our finding that Pol δ dissociates at specific CFS sequences is significant, since dissociation of the replication machinery and inability to efficiently recover the replication fork can lead to fork collapse and/or formation of double-strand breaks in vivo. Our biochemical studies also extend the potential involvement of Y-family polymerases in CFS maintenance to include polymerase κ.