Inverted DNA repeats: a source of eukaryotic genomic instability

Protein innovation through template switching in the Saccharomyces cerevisiae lineage

DNA polymerase template switching between short, non-identical inverted repeats (IRs) is a genetic mechanism that leads to the homogenization of IR arms and to IR spacer inversion, which cause multinucleotide mutations (MNMs). It is unknown if and how template switching affects gene evolution. In this study, we performed a phylogenetic analysis to determine the effect of template switching between IR arms on coding DNA of Saccharomyces cerevisiae. To achieve this, perfect IRs that co-occurred with MNMs between a strain and its parental node were identified in S. cerevisiae strains. We determined that template switching introduced MNMs into 39 protein-coding genes through S. cerevisiae evolution, resulting in both arm homogenization and inversion of the IR spacer. These events in turn resulted in nonsynonymous substitutions and up to five neighboring amino acid replacements in a single gene. The study demonstrates that template switching is a powerful generator of multiple substitutions within codons. Additionally, some template switching events occurred more than once during S. cerevisiae evolution. Our findings suggest that template switching constitutes a general mutagenic mechanism that results in both nonsynonymous substitutions and parallel evolution, which are traditionally considered as evidence for positive selection, without the need for adaptive explanations.

Origin and Fates of TERT Gene Copies in Polyploid Plants

The gene coding for the telomerase reverse transcriptase (TERT) is essential for the maintenance of telomeres. Previously we described the presence of three TERT paralogs in the allotetraploid plant Nicotiana tabacum, while a single TERT copy was identified in the paleopolyploid model plant Arabidopsis thaliana. Here we examine the presence, origin and functional status of TERT variants in allotetraploid Nicotiana species of diverse evolutionary ages and their parental genome donors, as well as in other diploid and polyploid plant species. A combination of experimental and in silico bottom-up analyses of TERT gene copies in Nicotiana polyploids revealed various patterns of retention or loss of parental TERT variants and divergence in their functions. RT-qPCR results confirmed the expression of all the identified TERT variants. In representative plant and green algal genomes, our synteny analyses show that their TERT genes were located in a conserved locus that became advantageous after the divergence of eudicots, and the gene was later translocated in several plant groups. In various diploid and polyploid species, translocation of TERT became fixed in target loci that show ancient synapomorphy.

Exploring the Relationship Among Divergence Time and Coding and Non-coding Elements in the Shaping of Fungal Mitochondrial Genomes

The order Hypocreales (Ascomycota) is composed of ubiquitous and ecologically diverse fungi such as saprobes, biotrophs, and pathogens. Despite their phylogenetic relationship, these species exhibit high variability in biomolecules production, lifestyle, and fitness. The mitochondria play an important role in the fungal biology, providing energy to the cells and regulating diverse processes, such as immune response. In spite of its importance, the mechanisms that shape fungal mitogenomes are still poorly understood. Herein, we investigated the variability and evolution of mitogenomes and its relationship with the divergence time using the order Hypocreales as a study model. We sequenced and annotated for the first time Trichoderma harzianum mitochondrial genome (mtDNA), which was compared to other 34 mtDNAs species that were publicly available. Comparative analysis revealed a substantial structural and size variation on non-coding mtDNA regions, despite the conservation of copy number, length, and structure of protein-coding elements. Interestingly, we observed a highly significant correlation between mitogenome length, and the number and size of non-coding sequences in mitochondrial genome. Among the non-coding elements, group I and II introns and homing endonucleases genes (HEGs) were the main contributors to discrepancies in mitogenomes structure and length. Several intronic sequences displayed sequence similarity among species, and some of them are conserved even at gene position, and were present in the majority of mitogenomes, indicating its origin in a common ancestor. On the other hand, we also identified species-specific introns that advocate for the origin by different mechanisms. Investigation of mitochondrial gene transfer to the nuclear genome revealed that nuclear copies of the nad5 are the most frequent while atp8, atp9, and cox3 could not be identified in any of the nuclear genomes analyzed. Moreover, we also estimated the divergence time of each species and investigated its relationship with coding and non-coding elements as well as with the length of mitogenomes. Altogether, our results demonstrated that introns and HEGs are key elements on mitogenome shaping and its presence on fast-evolving mtDNAs could be mostly explained by its divergence time, although the intron sharing profile suggests the involvement of other mechanisms on the mitochondrial genome evolution, such as horizontal transference.

Structural variation and its potential impact on genome instability: Novel discoveries in the EGFR landscape by long-read sequencing

Structural variation (SV) is typically defined as variation within the human genome that exceeds 50 base pairs (bp). SV may be copy number neutral or it may involve duplications, deletions, and complex rearrangements. Recent studies have shown SV to be associated with many human diseases. However, studies of SV have been challenging due to technological constraints. With the advent of third generation (long-read) sequencing technology, exploration of longer stretches of DNA not easily examined previously has been made possible. In the present study, we utilized third generation (long-read) sequencing techniques to examine SV in the EGFR landscape of four haplotypes derived from two human samples. We analyzed the EGFR gene and its landscape (+/- 500,000 base pairs) using this approach and were able to identify a region of non-coding DNA with over 90% similarity to the most common activating EGFR mutation in non-small cell lung cancer. Based on previously published Alu-element genome instability algorithms, we propose a molecular mechanism to explain how this non-coding region of DNA may be interacting with and impacting the stability of the EGFR gene and potentially generating this cancer-driver gene. By these techniques, we were also able to identify previously hidden structural variation in the four haplotypes and in the human reference genome (hg38). We applied previously published algorithms to compare the relative stabilities of these five different EGFR gene landscape haplotypes to estimate their relative potentials to generate the EGFR exon 19, 15 bp canonical deletion. To our knowledge, the present study is the first to use the differences in genomic architecture between targeted cancer-linked phased haplotypes to estimate their relative potentials to form a common cancer-linked driver mutation.

DNA replication and repair kinetics of Alu, LINE-1 and satellite III genomic repetitive elements

Background

Preservation of genome integrity by complete, error-free DNA duplication prior to cell division and by correct DNA damage repair is paramount for the development and maintenance of an organism. This holds true not only for protein-encoding genes, but also it applies to repetitive DNA elements, which make up more than half of the human genome. Here, we focused on the replication and repair kinetics of interspersed and tandem repetitive DNA elements.

Results

We integrated genomic population level data with a single cell immunofluorescence in situ hybridization approach to simultaneously label replication/repair and repetitive DNA elements. We found that: (1) the euchromatic Alu element was replicated during early S-phase; (2) LINE-1, which is associated with AT-rich genomic regions, was replicated throughout S-phase, with the majority being replicated according to their particular histone marks; (3) satellite III, which constitutes pericentromeric heterochromatin, was replicated exclusively during the mid-to-late S-phase. As for the DNA double-strand break repair process, we observed that Alu elements followed the global genome repair kinetics, while LINE-1 elements repaired at a slower rate. Finally, satellite III repeats were repaired at later time points.

Conclusions

We conclude that the histone modifications in the specific repeat element predominantly determine its replication and repair timing. Thus, Alu elements, which are characterized by euchromatic chromatin features, are repaired and replicated the earliest, followed by LINE-1 elements, including more variegated eu/heterochromatic features and, lastly, satellite tandem repeats, which are homogeneously characterized by heterochromatic features and extend over megabase-long genomic regions. Altogether, this work reemphasizes the need for complementary approaches to achieve an integrated and comprehensive investigation of genomic processes.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0226-9) contains supplementary material, which is available to authorized users.

Differential DNA methylation and gene expression in reciprocal hybrids between Solanum lycopersicum and S. pimpinellifolium

Wide hybridization is a common and efficient breeding strategy for enhancing crop yield and quality. An interesting phenomenon is that the reciprocal hybrids usually show different phenotypes, and its underlying mechanism is not well understood. Here, we reported our comparative analysis of the DNA methylation patterns in Solanum lycopersicum, Solanum pimpinellifolium and their reciprocal hybrids by methylated DNA immunoprecipitation sequencing. The reciprocal hybrids had lower levels of DNA methylation in CpG islands and LTR retroelements when compared with those of their parents. Importantly, remarkable differences in DNA methylation patterns, mainly in introns and CDS regions, were revealed between the reciprocal hybrids. These different methylated regions were mapped to 79 genes, 14 of which were selected for analysis of gene expression levels. While there was an inverse correlation between DNA methylation and gene expression in promoter regions, the relationship was complicated in gene body regions. Further association analysis revealed that there were 15 differentially methylated genes associated with siRNAs, and that the methylation levels of these genes were inversely correlated with respective siRNAs. All these data raised the possibility that the direction of hybridization induced the divergent epigenomes leading to changes in the transcription levels of reciprocal hybrids.

Complex Analyses of Short Inverted Repeats in All Sequenced Chloroplast DNAs

Chloroplasts are key organelles in the management of oxygen in algae and plants and are therefore crucial for all living beings that consume oxygen. Chloroplasts typically contain a circular DNA molecule with nucleus-independent replication and heredity. Using "palindrome analyser" we performed complete analyses of short inverted repeats (S-IRs) in all chloroplast DNAs (cpDNAs) available from the NCBI genome database. Our results provide basic parameters of cpDNAs including comparative information on localization, frequency, and differences in S-IR presence. In a total of 2,565 cpDNA sequences available, the average frequency of S-IRs in cpDNA genomes is 45 S-IRs/per kbp, significantly higher than that found in mitochondrial DNA sequences. The frequency of S-IRs in cpDNAs generally decreased with S-IR length, but not for S-IRs 15, 22, 24, or 27 bp long, which are significantly more abundant than S-IRs with other lengths. These results point to the importance of specific S-IRs in cpDNA genomes. Moreover, comparison by Levenshtein distance of S-IR similarities showed that a limited number of S-IR sequences are shared in the majority of cpDNAs. S-IRs are not located randomly in cpDNAs, but are length-dependently enriched in specific locations, including the repeat region, stem, introns, and tRNA regions. The highest enrichment was found for 12 bp and longer S-IRs in the stem-loop region followed by 12 bp and longer S-IRs located before the repeat region. On the other hand, S-IRs are relatively rare in rRNA sequences and around introns. These data show nonrandom and conserved arrangements of S-IRs in chloroplast genomes.

The Prevalence and Evolutionary Conservation of Inverted Repeats in Proteobacteria

Perfect short inverted repeats (IRs) are known to be enriched in a variety of bacterial and eukaryotic genomes. Currently, it is unclear whether perfect IRs are conserved over evolutionary time scales. In this study, we aimed to characterize the prevalence and evolutionary conservation of IRs across 20 proteobacterial strains. We first identified IRs in Escherichia coli K-12 substr MG1655 and showed that they are overabundant. We next aimed to test whether this overabundance is reflected in the conservation of IRs over evolutionary time scales. To this end, for each perfect IR identified in E. coli MG1655, we collected orthologous sequences from related proteobacterial genomes. We next quantified the evolutionary conservation of these IRs, that is, the presence of the exact same IR across orthologous regions. We observed high conservation of perfect IRs: out of the 234 examined orthologous regions, 145 were more conserved than expected, which is statistically significant even after correcting for multiple testing. Our results together with previous experimental findings support a model in which imperfect IRs are corrected to perfect IRs in a preferential manner via a template switching mechanism.

Pathways and Mechanisms that Prevent Genome Instability in Saccharomyces cerevisiae

Genome rearrangements result in mutations that underlie many human diseases, and ongoing genome instability likely contributes to the development of many cancers. The tools for studying genome instability in mammalian cells are limited, whereas model organisms such as Saccharomyces cerevisiae are more amenable to these studies. Here, we discuss the many genetic assays developed to measure the rate of occurrence of Gross Chromosomal Rearrangements (called GCRs) in S. cerevisiae. These genetic assays have been used to identify many types of GCRs, including translocations, interstitial deletions, and broken chromosomes healed by de novo telomere addition, and have identified genes that act in the suppression and formation of GCRs. Insights from these studies have contributed to the understanding of pathways and mechanisms that suppress genome instability and how these pathways cooperate with each other. Integrated models for the formation and suppression of GCRs are discussed.

Direct and Inverted Repeat stimulated excision (DIRex): Simple, single-step, and scar-free mutagenesis of bacterial genes

The need for generating precisely designed mutations is common in genetics, biochemistry, and molecular biology. Here, I describe a new λ Red recombineering method (Direct and Inverted Repeat stimulated excision; DIRex) for fast and easy generation of single point mutations, small insertions or replacements as well as deletions of any size, in bacterial genes. The method does not leave any resistance marker or scar sequence and requires only one transformation to generate a semi-stable intermediate insertion mutant. Spontaneous excision of the intermediate efficiently and accurately generates the final mutant. In addition, the intermediate is transferable between strains by generalized transductions, enabling transfer of the mutation into multiple strains without repeating the recombineering step. Existing methods that can be used to accomplish similar results are either (i) more complicated to design, (ii) more limited in what mutation types can be made, or (iii) require expression of extrinsic factors in addition to λ Red. I demonstrate the utility of the method by generating several deletions, small insertions/replacements, and single nucleotide exchanges in Escherichia coli and Salmonella enterica. Furthermore, the design parameters that influence the excision frequency and the success rate of generating desired point mutations have been examined to determine design guidelines for optimal efficiency.

Short Inverted Repeats Are Hotspots for Genetic Instability: Relevance to Cancer Genomes

Analyses of chromosomal aberrations in human genetic disorders have revealed that inverted repeat sequences (IRs) often co-localize with endogenous chromosomal instability and breakage hotspots. Approximately 80% of all IRs in the human genome are short (<100 bp), yet the mutagenic potential of such short cruciform-forming sequences has not been characterized. Here, we find that short IRs are enriched at translocation breakpoints in human cancer and stimulate the formation of DNA double-strand breaks (DSBs) and deletions in mammalian and yeast cells. We provide evidence for replication-related mechanisms of IR-induced genetic instability and a novel XPF cleavage-based mechanism independent of DNA replication. These discoveries implicate short IRs as endogenous sources of DNA breakage involved in disease etiology and suggest that these repeats represent a feature of genome plasticity that may contribute to the evolution of the human genome by providing a means for diversity within the population.

Correlations between long inverted repeat (LIR) features, deletion size and distance from breakpoint in human gross gene deletions

Long inverted repeats (LIRs) have been shown to induce genomic deletions in yeast. In this study, LIRs were investigated within ±10 kb spanning each breakpoint from 109 human gross deletions, using Inverted Repeat Finder (IRF) software. LIR number was significantly higher at the breakpoint regions, than in control segments (P < 0.001). In addition, it was found that strong correlation between 5' and 3' LIR numbers, suggesting contribution to DNA sequence evolution (r = 0.85, P < 0.001). 138 LIR features at ±3 kb breakpoints in 89 (81%) of 109 gross deletions were evaluated. Significant correlations were found between distance from breakpoint and loop length (r = -0.18, P < 0.05) and stem length (r = -0.18, P < 0.05), suggesting DNA strands are potentially broken in locations closer to bigger LIRs. In addition, bigger loops cause larger deletions (r = 0.19, P < 0.05). Moreover, loop length (r = 0.29, P < 0.02) and identity between stem copies (r = 0.30, P < 0.05) of 3' LIRs were more important in larger deletions. Consequently, DNA breaks may form via LIR-induced cruciform structure during replication. DNA ends may be later repaired by non-homologous end-joining (NHEJ), with following deletion.

High variability of mitochondrial gene order among fungi

From their origin as an early alpha proteobacterial endosymbiont to their current state as cellular organelles, large-scale genomic reorganization has taken place in the mitochondria of all main eukaryotic lineages. So far, most studies have focused on plant and animal mitochondrial (mt) genomes (mtDNA), but fungi provide new opportunities to study highly differentiated mtDNAs. Here, we analyzed 38 complete fungal mt genomes to investigate the evolution of mtDNA gene order among fungi. In particular, we looked for evidence of nonhomologous intrachromosomal recombination and investigated the dynamics of gene rearrangements. We investigated the effect that introns, intronic open reading frames (ORFs), and repeats may have on gene order. Additionally, we asked whether the distribution of transfer RNAs (tRNAs) evolves independently to that of mt protein-coding genes. We found that fungal mt genomes display remarkable variation between and within the major fungal phyla in terms of gene order, genome size, composition of intergenic regions, and presence of repeats, introns, and associated ORFs. Our results support previous evidence for the presence of mt recombination in all fungal phyla, a process conspicuously lacking in most Metazoa. Overall, the patterns of rearrangements may be explained by the combined influences of recombination (i.e., most likely nonhomologous and intrachromosomal), accumulated repeats, especially at intergenic regions, and to a lesser extent, mobile element dynamics.

Adeno-associated virus inverted terminal repeats stimulate gene editing

Advancements in genome editing have relied on technologies to specifically damage DNA which, in turn, stimulates DNA repair including homologous recombination (HR). As off-target concerns complicate the therapeutic translation of site-specific DNA endonucleases, an alternative strategy to stimulate gene editing based on fragile DNA was investigated. To do this, an episomal gene-editing reporter was generated by a disruptive insertion of the adeno-associated virus (AAV) inverted terminal repeat (ITR) into the egfp gene. Compared with a non-structured DNA control sequence, the ITR induced DNA damage as evidenced by increased gamma-H2AX and Mre11 foci formation. As local DNA damage stimulates HR, ITR-mediated gene editing was investigated using DNA oligonucleotides as repair substrates. The AAV ITR stimulated gene editing >1000-fold in a replication-independent manner and was not biased by the polarity of the repair oligonucleotide. Analysis of additional human DNA sequences demonstrated stimulation of gene editing to varying degrees. In particular, inverted yet not direct, Alu repeats induced gene editing, suggesting a role for DNA structure in the repair event. Collectively, the results demonstrate that inverted DNA repeats stimulate gene editing via double-strand break repair in an episomal context and allude to efficient gene editing of the human chromosome using fragile DNA sequences.

detectIR: a novel program for detecting perfect and imperfect inverted repeats using complex numbers and vector calculation

Inverted repeats are present in abundance in both prokaryotic and eukaryotic genomes and can form DNA secondary structures--hairpins and cruciforms that are involved in many important biological processes. Bioinformatics tools for efficient and accurate detection of inverted repeats are desirable, because existing tools are often less accurate and time consuming, sometimes incapable of dealing with genome-scale input data. Here, we present a MATLAB-based program called detectIR for the perfect and imperfect inverted repeat detection that utilizes complex numbers and vector calculation and allows genome-scale data inputs. A novel algorithm is adopted in detectIR to convert the conventional sequence string comparison in inverted repeat detection into vector calculation of complex numbers, allowing non-complementary pairs (mismatches) in the pairing stem and a non-palindromic spacer (loop or gaps) in the middle of inverted repeats. Compared with existing popular tools, our program performs with significantly higher accuracy and efficiency. Using genome sequence data from HIV-1, Arabidopsis thaliana, Homo sapiens and Zea mays for comparison, detectIR can find lots of inverted repeats missed by existing tools whose outputs often contain many invalid cases. detectIR is open source and its source code is freely available at: https://sourceforge.net/projects/detectir.

Palindromic GOLGA8 core duplicons promote chromosome 15q13.3 microdeletion and evolutionary instability

Recurrent deletions of chromosome 15q13.3 associate with intellectual disability, schizophrenia, autism and epilepsy. To gain insight into its instability, we sequenced the region in patients, normal individuals and nonhuman primates. We discovered five structural configurations of the human chromosome 15q13.3 region ranging in size from 2 to 3 Mbp. These configurations arose recently (~0.5–0.9 million years ago) as a result of human-specific expansions of segmental duplications and two independent inversion events. All inversion breakpoints map near GOLGA8 core duplicons—a ~14 kbp primate-specific chromosome 15 repeat that became organized into larger palindromic structures. GOLGA8-flanked palindromes also demarcate the breakpoints of recurrent 15q13.3 microdeletions, the expansion of chromosome 15 segmental duplications in the human lineage, and independent structural changes in apes. The significant clustering (p=0.002) of breakpoints provides mechanistic evidence for the role of this core duplicon and its palindromic architecture in promoting evolutionary and disease-related instability of chromosome 15.

Repetitive sequences in plant nuclear DNA: types, distribution, evolution and function

Repetitive DNA sequences are a major component of eukaryotic genomes and may account for up to 90% of the genome size. They can be divided into minisatellite, microsatellite and satellite sequences. Satellite DNA sequences are considered to be a fast-evolving component of eukaryotic genomes, comprising tandemly-arrayed, highly-repetitive and highly-conserved monomer sequences. The monomer unit of satellite DNA is 150-400 base pairs (bp) in length. Repetitive sequences may be species- or genus-specific, and may be centromeric or subtelomeric in nature. They exhibit cohesive and concerted evolution caused by molecular drive, leading to high sequence homogeneity. Repetitive sequences accumulate variations in sequence and copy number during evolution, hence they are important tools for taxonomic and phylogenetic studies, and are known as "tuning knobs" in the evolution. Therefore, knowledge of repetitive sequences assists our understanding of the organization, evolution and behavior of eukaryotic genomes. Repetitive sequences have cytoplasmic, cellular and developmental effects and play a role in chromosomal recombination. In the post-genomics era, with the introduction of next-generation sequencing technology, it is possible to evaluate complex genomes for analyzing repetitive sequences and deciphering the yet unknown functional potential of repetitive sequences.

Analysis of the t(3;8) of hereditary renal cell carcinoma: a palindrome-mediated translocation

It has emerged that palindrome-mediated genomic instability generates DNA-based rearrangements. The presence of palindromic AT-rich repeats (PATRRs) at the translocation breakpoints suggested a palindrome-mediated mechanism in the generation of several recurrent constitutional rearrangements: the t(11;22), t(17;22), and t(8;22). To date, all reported PATRR-mediated translocations include the PATRR on chromosome 22 (PATRR22) as a translocation partner. Here, the constitutional rearrangement, t(3;8)(p14.2;q24.1), segregating with renal cell carcinoma in two families, is examined. The chromosome 8 breakpoint lies in PATRR8 in the first intron of the RNF139 (TRC8) gene, whereas the chromosome 3 breakpoint is located in an AT-rich palindromic sequence in intron 3 of the FHIT gene (PATRR3). Thus, the t(3;8) is the first PATRR-mediated, recurrent, constitutional translocation that does not involve PATRR22. Furthermore, we detect de novo translocations similar to the t(11;22) and t(8;22), involving PATRR3 in normal sperm. The breakpoint on chromosome 3 is in proximity to FRA3B, the most common fragile site in the human genome and a site of frequent deletions in tumor cells. However, the lack of involvement of PATRR3 sequence in numerous FRA3B-related deletions suggests that there are several different DNA sequence-based etiologies responsible for chromosome 3p14.2 genomic rearrangements.

Genome-wide screen reveals replication pathway for quasi-palindrome fragility dependent on homologous recombination

Inverted repeats capable of forming hairpin and cruciform structures present a threat to chromosomal integrity. They induce double strand breaks, which lead to gross chromosomal rearrangements, the hallmarks of cancers and hereditary diseases. Secondary structure formation at this motif has been proposed to be the driving force for the instability, albeit the mechanisms leading to the fragility are not well-understood. We carried out a genome-wide screen to uncover the genetic players that govern fragility of homologous and homeologous Alu quasi-palindromes in the yeast Saccharomyces cerevisiae. We found that depletion or lack of components of the DNA replication machinery, proteins involved in Fe-S cluster biogenesis, the replication-pausing checkpoint pathway, the telomere maintenance complex or the Sgs1-Top3-Rmi1 dissolvasome augment fragility at Alu-IRs. Rad51, a component of the homologous recombination pathway, was found to be required for replication arrest and breakage at the repeats specifically in replication-deficient strains. These data demonstrate that Rad51 is required for the formation of breakage-prone secondary structures in situations when replication is compromised while another mechanism operates in DSB formation in replication-proficient strains.

Inverted repeats are found in many eukaryotic genomes including humans. They have a potential to cause chromosomal breakage and rearrangements that contribute to genome polymorphism and the development of diseases. Instability of inverted repeats is accounted for by their propensity to adopt DNA secondary structures that is negatively affected by the distance between the repeats and level of sequence divergence. However, the genetic factors that promote the abnormal structure formation or affect the ability of the repeats to break are largely unknown. Here, using a genome-wide screen we identified 38 mutants that destabilize imperfect human inverted Alu repeats and predispose them to breakage. The proteins that are required to maintain repeat stability belong to the core of the DNA replication machinery and to the accessory proteins that help replication fork to move through the difficult templates. Remarkably, when replication machinery is compromised, the proteins involved in homologous recombination promote the formation of secondary structures and replication block thereby triggering breakage at the inverted repeats. These results reveal a powerful pathway for the destabilization of chromosomes containing inverted repeats that requires the activity of homologous recombination.

Size of gene specific inverted repeat--dependent gene deletion In Saccharomyces cerevisiae

We describe here an approach for rapidly producing scar-free and precise gene deletions in S. cerevisiae with high efficiency. Preparation of the disruption gene cassette in this approach was simply performed by overlap extension-PCR of an invert repeat of a partial or complete sequence of the targeted gene with URA3. Integration of the prepared disruption gene cassette to the designated position of a target gene leads to the formation of a mutagenesis cassette within the yeast genome, which consists of a URA3 gene flanked by the targeted gene and its inverted repeat between two short identical direct repeats. The inherent instability of the inverted sequences in close proximity facilitates the self-excision of the entire mutagenesis cassette deposited in the genome and promotes homologous recombination resulting in a seamless deletion via a single transformation. This rapid assembly circumvents the difficulty during preparation of disruption gene cassettes composed of two inverted repeats of the URA3, which requires the engineering of unique restriction sites for subsequent digestion and T4 DNA ligation in vitro. We further identified that the excision of the entire mutagenesis cassette flanked by two DRs in the transformed S. cerevisiae is dependent on the length of the inverted repeat of which a minimum of 800 bp is required for effective gene deletion. The deletion efficiency improves with the increase of the inverted repeat till 1.2 kb. Finally, the use of gene-specific inverted repeats of target genes enables simultaneous gene deletions. The procedure has the potential for application on other yeast strains to achieve precise and efficient removal of gene sequences.

High-resolution mapping of spontaneous mitotic recombination hotspots on the 1.1 Mb arm of yeast chromosome IV

Although homologous recombination is an important pathway for the repair of double-stranded DNA breaks in mitotically dividing eukaryotic cells, these events can also have negative consequences, such as loss of heterozygosity (LOH) of deleterious mutations. We mapped about 140 spontaneous reciprocal crossovers on the right arm of the yeast chromosome IV using single-nucleotide-polymorphism (SNP) microarrays. Our mapping and subsequent experiments demonstrate that inverted repeats of Ty retrotransposable elements are mitotic recombination hotspots. We found that the mitotic recombination maps on the two homologs were substantially different and were unrelated to meiotic recombination maps. Additionally, about 70% of the DNA lesions that result in LOH are likely generated during G1 of the cell cycle and repaired during S or G2. We also show that different genetic elements are associated with reciprocal crossover conversion tracts depending on the cell cycle timing of the initiating DSB.

Double-strand breaks (DSBs) are DNA lesions that can be fatal to a cell if left unrepaired. They can be caused by exogenous sources, such as gamma radiation, or endogenous stresses, such as high levels of transcription. Yeast cells primarily repair DSBs that are initiated outside of meiosis by mitotic recombination, which can result in physical exchanges between chromosomes, known as crossovers. We created a mitotic recombination map of one chromosome arm, representing 10% of the genome. This recombination map allows us to determine which regions of the chromosome arm are more susceptible to DNA damage than other regions. We were able to determine that most DSBs that result in detectable genomic changes were initiated prior to DNA replication and that some secondary DNA structures can be recombination hotspots. Recombination can also occur during meiosis, as a method of ensuring proper chromosome segregation. However, previously reported meiotic recombination maps have no correlation with our mitotic recombination map.

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.

The yin and yang of repair mechanisms in DNA structure-induced genetic instability

DNA can adopt a variety of secondary structures that deviate from the canonical Watson-Crick B-DNA form. More than 10 types of non-canonical or non-B DNA secondary structures have been characterized, and the sequences that have the capacity to adopt such structures are very abundant in the human genome. Non-B DNA structures have been implicated in many important biological processes and can serve as sources of genetic instability, implicating them in disease and evolution. Non-B DNA conformations interact with a wide variety of proteins involved in replication, transcription, DNA repair, and chromatin architectural regulation. In this review, we will focus on the interactions of DNA repair proteins with non-B DNA and their roles in genetic instability, as the proteins and DNA involved in such interactions may represent plausible targets for selective therapeutic intervention.

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.

Multiple pathways of duplication formation with and without recombination (RecA) in Salmonella enterica

Duplications are often attributed to "unequal recombination" between separated, directly repeated sequence elements (>100 bp), events that leave a recombinant element at the duplication junction. However, in the bacterial chromosome, duplications form at high rates (10(-3)-10(-5)/cell/division) even without recombination (RecA). Here we describe 1800 spontaneous lac duplications trapped nonselectively on the low-copy F'(128) plasmid, where lac is flanked by direct repeats of the transposable element IS3 (1258 bp) and by numerous quasipalindromic REP elements (30 bp). Duplications form at a high rate (10(-4)/cell/division) that is reduced only about 11-fold in the absence of RecA. With and without RecA, most duplications arise by recombination between IS3 elements (97%). Formation of these duplications is stimulated by IS3 transposase (Tnp) and plasmid transfer functions (TraI). Three duplication pathways are proposed. First, plasmid dimers form at a high rate stimulated by RecA and are then modified by deletions between IS3 elements (resolution) that leave a monomeric plasmid with an IS3-flanked lac duplication. Second, without RecA, duplications occur by single-strand annealing of DNA ends generated in different sister chromosomes after transposase nicks DNA near participating IS3 elements. The absence of RecA may stimulate annealing by allowing chromosome breaks to persist. Third, a minority of lac duplications (3%) have short (0-36 bp) junction sequences (SJ), some of which are located within REP elements. These duplication types form without RecA, Tnp, or Tra by a pathway in which the palindromic junctions of a tandem inversion duplication (TID) may stimulate deletions that leave the final duplication.

Small interfering RNA-producing loci in the ancient parasitic eukaryote Trypanosoma brucei

Background

At the core of the RNA interference (RNAi) pathway in Trypanosoma brucei is a single Argonaute protein, TbAGO1, with an established role in controlling retroposon and repeat transcripts. Recent evidence from higher eukaryotes suggests that a variety of genomic sequences with the potential to produce double-stranded RNA are sources for small interfering RNAs (siRNAs).

Results

To test whether such endogenous siRNAs are present in T. brucei and to probe the individual role of the two Dicer-like enzymes, we affinity purified TbAGO1 from wild-type procyclic trypanosomes, as well as from cells deficient in the cytoplasmic (TbDCL1) or nuclear (TbDCL2) Dicer, and subjected the bound RNAs to Illumina high-throughput sequencing. In wild-type cells the majority of reads originated from two classes of retroposons. We also considerably expanded the repertoire of trypanosome siRNAs to encompass a family of 147-bp satellite-like repeats, many of the regions where RNA polymerase II transcription converges, large inverted repeats and two pseudogenes. Production of these newly described siRNAs is strictly dependent on the nuclear DCL2. Notably, our data indicate that putative centromeric regions, excluding the CIR147 repeats, are not a significant source for endogenous siRNAs.

Conclusions

Our data suggest that endogenous RNAi targets may be as evolutionarily old as the mechanism itself.

Chromosomal translocations and palindromic AT-rich repeats

Repetitive DNA sequences constitute 30% of the human genome, and are often sites of genomic rearrangement. Recently, it has been found that several constitutional translocations, especially those that involve chromosome 22, take place utilizing palindromic sequences on 22q11 and on the partner chromosome. Analysis of translocation junction fragments shows that the breakpoints of such palindrome-mediated translocations are localized at the center of palindromic AT-rich repeats (PATRRs). The presence of PATRRs at the breakpoints indicates a palindrome-mediated mechanism involved in the generation of these constitutional translocations. Identification of these PATRR-mediated translocations suggests a universal pathway for gross chromosomal rearrangement in the human genome. De novo occurrences of PATRR-mediated translocations can be detected by PCR in normal sperm samples but not somatic cells. Polymorphisms of various PATRRs influence their propensity for adopting a secondary structure, which in turn affects de novo translocation frequency. We propose that the PATRRs form an unstable secondary structure, which leads to double-strand breaks at the center of the PATRR. The double-strand breaks appear to be followed by a non-homologous end-joining repair pathway, ultimately leading to the translocations. This review considers recent findings concerning the mechanism of meiosis-specific, PATRR-mediated translocations.

High resolution methylome map of rat indicates role of intragenic DNA methylation in identification of coding region

DNA methylation is crucial for gene regulation and maintenance of genomic stability. Rat has been a key model system in understanding mammalian systemic physiology, however detailed rat methylome remains uncharacterized till date. Here, we present the first high resolution methylome of rat liver generated using Methylated DNA immunoprecipitation and high throughput sequencing (MeDIP-Seq) approach. We observed that within the DNA/RNA repeat elements, simple repeats harbor the highest degree of methylation. Promoter hypomethylation and exon hypermethylation were common features in both RefSeq genes and expressed genes (as evaluated by proteomic approach). We also found that although CpG islands were generally hypomethylated, about 6% of them were methylated and a large proportion (37%) of methylated islands fell within the exons. Notably, we obeserved significant differences in methylation of terminal exons (UTRs); methylation being more pronounced in coding/partially coding exons compared to the non-coding exons. Further, events like alternate exon splicing (cassette exon) and intron retentions were marked by DNA methylation and these regions are retained in the final transcript. Thus, we suggest that DNA methylation could play a crucial role in marking coding regions thereby regulating alternative splicing. Apart from generating the first high resolution methylome map of rat liver tissue, the present study provides several critical insights into methylome organization and extends our understanding of interplay between epigenome, gene expression and genome stability.

Cre-mediated recombination can induce apoptosis in vivo by activating the p53 DNA damage-induced pathway

Cre-mediated apoptosis has been observed in many contexts in mice expressing Cre-recombinase and can confound the analysis of genetically engineered conditional mutant or transgenic alleles. Several mechanisms have been proposed to explain this phenomenon. We find that the degree of apoptosis induced correlates roughly with the copy number of loxP sites present in the genome and that some level of increased apoptosis accompanies the presence of even only a few loxP sites, as occurs in conditional floxed alleles. Cre-induced apoptosis in this context is completely p53-dependent, suggesting that the apoptosis is stimulated by p53 activation in response to DNA damage incurred during the process of Cre-mediated recombination.

FLNA genomic rearrangements cause periventricular nodular heterotopia

OBJECTIVE

To identify copy number variant (CNV) causes of periventricular nodular heterotopia (PNH) in patients for whom FLNA sequencing is negative.

METHODS

Screening of 35 patients from 33 pedigrees on an Affymetrix 6.0 microarray led to the identification of one individual bearing a CNV that disrupted FLNA. FLNA-disrupting CNVs were also isolated in 2 other individuals by multiplex ligation probe amplification. These 3 cases were further characterized by high-resolution oligo array comparative genomic hybridization (CGH), and the precise junctional breakpoints of the rearrangements were identified by PCR amplification and sequencing.

RESULTS

We report 3 cases of PNH caused by nonrecurrent genomic rearrangements that disrupt one copy of FLNA. The first individual carried a 113-kb deletion that removes all but the first exon of FLNA. A second patient harbored a complex rearrangement including a deletion of the 3' end of FLNA accompanied by a partial duplication event. A third patient bore a 39-kb deletion encompassing all of FLNA and the neighboring gene EMD. High-resolution oligo array CGH of the FLNA locus suggests distinct molecular mechanisms for each of these rearrangements, and implicates nearby low copy repeats in their pathogenesis.

CONCLUSIONS

These results demonstrate that FLNA is prone to pathogenic rearrangements, and highlight the importance of screening for CNVs in individuals with PNH lacking FLNA point mutations.

Non-B DNA Secondary Structures and Their Resolution by RecQ Helicases

In addition to the canonical B-form structure first described by Watson and Crick, DNA can adopt a number of alternative structures. These non-B-form DNA secondary structures form spontaneously on tracts of repeat sequences that are abundant in genomes. In addition, structured forms of DNA with intrastrand pairing may arise on single-stranded DNA produced transiently during various cellular processes. Such secondary structures have a range of biological functions but also induce genetic instability. Increasing evidence suggests that genomic instabilities induced by non-B DNA secondary structures result in predisposition to diseases. Secondary DNA structures also represent a new class of molecular targets for DNA-interactive compounds that might be useful for targeting telomeres and transcriptional control. The equilibrium between the duplex DNA and formation of multistranded non-B-form structures is partly dependent upon the helicases that unwind (resolve) these alternate DNA structures. With special focus on tetraplex, triplex, and cruciform, this paper summarizes the incidence of non-B DNA structures and their association with genomic instability and emphasizes the roles of RecQ-like DNA helicases in genome maintenance by resolution of DNA secondary structures. In future, RecQ helicases are anticipated to be additional molecular targets for cancer chemotherapeutics.

The constitutional t(11;22): implications for a novel mechanism responsible for gross chromosomal rearrangements

The constitutional t(11;22)(q23;q11) is the most common recurrent non-Robertsonian translocation in humans. The breakpoint sequences of both chromosomes are characterized by several hundred base pairs of palindromic AT-rich repeats (PATRRs). Similar PATRRs have also been identified at the breakpoints of other nonrecurrent translocations, suggesting that PATRR-mediated chromosomal translocation represents one of the universal pathways for gross chromosomal rearrangement in the human genome. We propose that PATRRs have the potential to form cruciform structures through intrastrand-base pairing in single-stranded DNA, creating a source of genomic instability and leading to translocations. Indeed, de novo examples of the t(11;22) are detected at a high frequency in sperm from normal healthy males. This review synthesizes recent data illustrating a novel paradigm for an apparent spermatogenesis-specific translocation mechanism. This observation has important implications pertaining to the predominantly paternal origin of de novo gross chromosomal rearrangements in humans.

Comparative analyses of vertebrate posterior HoxD clusters reveal atypical cluster architecture in the caecilian Typhlonectes natans

BACKGROUND

The posterior genes of the HoxD cluster play a crucial role in the patterning of the tetrapod limb. This region is under the control of a global, long-range enhancer that is present in all vertebrates. Variation in limb types, as is the case in amphibians, can probably not only be attributed to variation in Hox genes, but is likely to be the product of differences in gene regulation. With a collection of vertebrate genome sequences available today, we used a comparative genomics approach to study the posterior HoxD cluster of amphibians. A frog and a caecilian were included in the study to compare coding sequences as well as to determine the gain and loss of putative regulatory sequences.

RESULTS

We sequenced the posterior end of the HoxD cluster of a caecilian and performed comparative analyses of this region using HoxD clusters of other vertebrates. We determined the presence of conserved non-coding sequences and traced gains and losses of these footprints during vertebrate evolution, with particular focus on amphibians. We found that the caecilian HoxD cluster is almost three times larger than its mammalian counterpart. This enlargement is accompanied with the loss of one gene and the accumulation of repeats in that area. A similar phenomenon was observed in the coelacanth, where a different gene was lost and expansion of the area where the gene was lost has occurred. At least one phylogenetic footprint present in all vertebrates was lost in amphibians. This conserved region is a known regulatory element and functions as a boundary element in neural tissue to prevent expression of Hoxd genes.

CONCLUSION

The posterior part of the HoxD cluster of Typhlonectes natans is among the largest known today. The loss of Hoxd-12 and the expansion of the intergenic region may exert an influence on the limb enhancer, by having to bypass a distance seven times that of regular HoxD clusters. Whether or not there is a correlation with the loss of limbs remains to be investigated. These results, together with data on other vertebrates show that the tetrapod Hox clusters are more variable than previously thought.

Non-B DB: a database of predicted non-B DNA-forming motifs in mammalian genomes

Although the capability of DNA to form a variety of non-canonical (non-B) structures has long been recognized, the overall significance of these alternate conformations in biology has only recently become accepted en masse. In order to provide access to genome-wide locations of these classes of predicted structures, we have developed non-B DB, a database integrating annotations and analysis of non-B DNA-forming sequence motifs. The database provides the most complete list of alternative DNA structure predictions available, including Z-DNA motifs, quadruplex-forming motifs, inverted repeats, mirror repeats and direct repeats and their associated subsets of cruciforms, triplex and slipped structures, respectively. The database also contains motifs predicted to form static DNA bends, short tandem repeats and homo(purine•pyrimidine) tracts that have been associated with disease. The database has been built using the latest releases of the human, chimp, dog, macaque and mouse genomes, so that the results can be compared directly with other data sources. In order to make the data interpretable in a genomic context, features such as genes, single-nucleotide polymorphisms and repetitive elements (SINE, LINE, etc.) have also been incorporated. The database is accessed through query pages that produce results with links to the UCSC browser and a GBrowse-based genomic viewer. It is freely accessible at http://nonb.abcc.ncifcrf.gov.

The distribution of inverted repeat sequences in the Saccharomyces cerevisiae genome

Although a variety of possible functions have been proposed for inverted repeat sequences (IRs), it is not known which of them might occur in vivo. We investigate this question by assessing the distributions and properties of IRs in the Saccharomyces cerevisiae (SC) genome. Using the IRFinder algorithm we detect 100,514 IRs having copy length greater than 6 bp and spacer length less than 77 bp. To assess statistical significance we also determine the IR distributions in two types of randomization of the S. cerevisiae genome. We find that the S. cerevisiae genome is significantly enriched in IRs relative to random. The S. cerevisiae IRs are significantly longer and contain fewer imperfections than those from the randomized genomes, suggesting that processes to lengthen and/or correct errors in IRs may be operative in vivo. The S. cerevisiae IRs are highly clustered in intergenic regions, while their occurrence in coding sequences is consistent with random. Clustering is stronger in the 3' flanks of genes than in their 5' flanks. However, the S. cerevisiae genome is not enriched in those IRs that would extrude cruciforms, suggesting that this is not a common event. Various explanations for these results are considered.

Participation of DNA polymerase zeta in replication of undamaged DNA in Saccharomyces cerevisiae

Translesion synthesis DNA polymerases contribute to DNA damage tolerance by mediating replication of damaged templates. Due to the low fidelity of these enzymes, lesion bypass is often mutagenic. We have previously shown that, in Saccharomyces cerevisiae, the contribution of the error-prone DNA polymerase zeta (Polzeta) to replication and mutagenesis is greatly enhanced if the normal replisome is defective due to mutations in replication genes. Here we present evidence that this defective-replisome-induced mutagenesis (DRIM) results from the participation of Polzeta in the copying of undamaged DNA rather than from mutagenic lesion bypass. First, DRIM is not elevated in strains that have a high level of endogenous DNA lesions due to defects in nucleotide excision repair or base excision repair pathways. Second, DRIM remains unchanged when the level of endogenous oxidative DNA damage is decreased by using anaerobic growth conditions. Third, analysis of the spectrum of mutations occurring during DRIM reveals the characteristic error signature seen during replication of undamaged DNA by Polzeta in vitro. These results extend earlier findings in Escherichia coli indicating that Y-family DNA polymerases can contribute to the copying of undamaged DNA. We also show that exposure of wild-type yeast cells to the replication inhibitor hydroxyurea causes a Polzeta-dependent increase in mutagenesis. This suggests that DRIM represents a response to replication impediment per se rather than to specific defects in the replisome components.

Chromosome aberrations resulting from double-strand DNA breaks at a naturally occurring yeast fragile site composed of inverted ty elements are independent of Mre11p and Sae2p

Genetic instability at palindromes and spaced inverted repeats (IRs) leads to chromosome rearrangements. Perfect palindromes and IRs with short spacers can extrude as cruciforms or fold into hairpins on the lagging strand during replication. Cruciform resolution produces double-strand breaks (DSBs) with hairpin-capped ends, and Mre11p and Sae2p are required to cleave the hairpin tips to facilitate homologous recombination. Fragile site 2 (FS2) is a naturally occurring IR in Saccharomyces cerevisiae composed of a pair of Ty1 elements separated by approximately 280 bp. Our results suggest that FS2 forms a hairpin, rather than a cruciform, during replication in cells with low levels of DNA polymerase. Cleavage of this hairpin results in a recombinogenic DSB. We show that DSB formation at FS2 does not require Mre11p, Sae2p, Rad1p, Slx4p, Pso2p, Exo1p, Mus81p, Yen1p, or Rad27p. Also, repair of DSBs by homologous recombination is efficient in mre11 and sae2 mutants. Homologous recombination is impaired at FS2 in rad52 mutants and most aberrations reflect either joining of two broken chromosomes in a "half crossover" or telomere capping of the break. In support of hairpin formation precipitating DSBs at FS2, two telomere-capped deletions had a breakpoint near the center of the IR. In summary, Mre11p and Sae2p are not required for DSB formation at FS2 or the subsequent repair of these DSBs.

Impaired DNA replication prompts deletions within palindromic sequences, but does not induce translocations in human cells

Palindromic regions are unstable and susceptible to deletion in prokaryotes and eukaryotes possibly due to stalled or slow replication. In the human genome, they also appear to become partially or completely deleted, while two palindromic AT-rich repeats (PATRR) contribute to known recurrent constitutional translocations. To explore the mechanism that causes the development of palindrome instabilities in humans, we compared the incidence of de novo translocations and deletions at PATRRs in human cells. Using a highly sensitive PCR assay that can detect single molecules, de novo deletions were detected neither in human somatic cells nor in sperm. However, deletions were detected at low frequency in cultured cell lines. Inhibition of DNA replication by administration of siRNA against the DNA polymerase alpha 1 (POLA1) gene or introduction of POLA inhibitors increased the frequency. This is in contrast to PATRR-mediated translocations that were never detected in similar conditions but were observed frequently in human sperm samples. Further deletions were found to take place during both leading- and lagging-strand synthesis. Our data suggest that stalled or slow replication induces deletions within PATRRs, but that other mechanisms might contribute to PATRR-mediated recurrent translocations in humans.

Methods to determine DNA structural alterations and genetic instability

Chromosomal DNA is a dynamic structure that can adopt a variety of non-canonical (i.e., non-B) conformations. In this regard, at least 10 different forms of non-B DNA conformations have been identified; many of them have been found to be mutagenic, and associated with human disease development. Despite the importance of non-B DNA structures in genetic instability and DNA metabolic processes, mechanisms by which instability occurs remain largely undefined. The purpose of this review is to summarize current methodologies that are used to address questions in the field of non-B DNA structure-induced genetic instability. Advantages and disadvantages of each method will be discussed. A focused effort to further elucidate the mechanisms of non-B DNA-induced genetic instability will lead to a better understanding of how these structure-forming sequences contribute to the development of human disease.

Chromosomal inversions between human and chimpanzee lineages caused by retrotransposons

The long interspersed element-1 (LINE-1 or L1) and Alu elements are the most abundant mobile elements comprising 21% and 11% of the human genome, respectively. Since the divergence of human and chimpanzee lineages, these elements have vigorously created chromosomal rearrangements causing genomic difference between humans and chimpanzees by either increasing or decreasing the size of genome. Here, we report an exotic mechanism, retrotransposon recombination-mediated inversion (RRMI), that usually does not alter the amount of genomic material present. Through the comparison of the human and chimpanzee draft genome sequences, we identified 252 inversions whose respective inversion junctions can clearly be characterized. Our results suggest that L1 and Alu elements cause chromosomal inversions by either forming a secondary structure or providing a fragile site for double-strand breaks. The detailed analysis of the inversion breakpoints showed that L1 and Alu elements are responsible for at least 44% of the 252 inversion loci between human and chimpanzee lineages, including 49 RRMI loci. Among them, three RRMI loci inverted exonic regions in known genes, which implicates this mechanism in generating the genomic and phenotypic differences between human and chimpanzee lineages. This study is the first comprehensive analysis of mobile element bases inversion breakpoints between human and chimpanzee lineages, and highlights their role in primate genome evolution.

Mutagenic inverted repeat assisted genome engineering (MIRAGE)

Here we describe a one-step method to create precise modifications in the genome of Saccharomyces cerevisiae as a tool for synthetic biology, metabolic engineering, systems biology and genetic studies. Through homologous recombination, a mutagenesis cassette containing an inverted repeat of selection marker(s) is integrated into the genome. Due to its inherent instability in genomic DNA, the inverted repeat catalyzes spontaneous self-excision, resulting in precise genome modification. Since this excision occurs at very high frequencies, selection for the integration event can be followed immediately by counterselection, without the need for growth in permissive conditions. This is the first time a truly one-step method has been described for genome modification in any organism.

Genomic complexity of the variable region-containing chitin-binding proteins in amphioxus

BACKGROUND

The variable region-containing chitin-binding proteins (VCBPs) are found in protochordates and consist of two tandem immunoglobulin variable (V)-type domains and a chitin-binding domain. We previously have shown that these polymorphic genes, which primarily are expressed in the gut, exhibit characteristics of immune genes. In this report, we describe VCBP genomic organization and characterize adjacent and intervening genetic features which may influence both their polymorphism and complex transcriptional repertoire.

RESULTS

VCBP genes 1, 2, 4, and 5 are encoded in a single contiguous gene-rich chromosomal region and VCBP3 is encoded in a separate locus. The VCBPs exhibit extensive haplotype variation, including copy number variation (CNV), indel polymorphism and a markedly elevated variation in repeat type and density. In at least one haplotype, inverted repeats occur more frequently than elsewhere in the genome. Multi-animal cDNA screening, as well as transcriptional profilingusing a novel transfection system, suggests that haplotype-specific transcriptional variants may contribute to VCBP genetic diversity.

CONCLUSION

The availability of the Branchiostoma floridae genome (Joint Genome Institute, Brafl1), along with BAC and PAC screening and sequencing described here, reveal that the relatively limited number of VCBP genes present in the amphioxus genome exhibit exceptionally high haplotype variation. These VCBP haplotypes contribute a diverse pool of allelic variants, which includes gene copy number variation, pseudogenes, and other polymorphisms, while contributing secondary effects on gene transcription as well.

Mutagenic and recombinagenic responses to defective DNA polymerase delta are facilitated by the Rev1 protein in pol3-t mutants of Saccharomyces cerevisiae

Defective DNA replication can result in substantial increases in the level of genome instability. In the yeast Saccharomyces cerevisiae, the pol3-t allele confers a defect in the catalytic subunit of replicative DNA polymerase delta that results in increased rates of mutagenesis, recombination, and chromosome loss, perhaps by increasing the rate of replicative polymerase failure. The translesion polymerases Pol eta, Pol zeta, and Rev1 are part of a suite of factors in yeast that can act at sites of replicative polymerase failure. While mutants defective in the translesion polymerases alone displayed few defects, loss of Rev1 was found to suppress the increased rates of spontaneous mutation, recombination, and chromosome loss observed in pol3-t mutants. These results suggest that Rev1 may be involved in facilitating mutagenic and recombinagenic responses to the failure of Pol delta. Genome stability, therefore, may reflect a dynamic relationship between primary and auxiliary DNA polymerases.

Chromosomal instability mediated by non-B DNA: cruciform conformation and not DNA sequence is responsible for recurrent translocation in humans

Chromosomal aberrations have been thought to be random events. However, recent findings introduce a new paradigm in which certain DNA segments have the potential to adopt unusual conformations that lead to genomic instability and nonrandom chromosomal rearrangement. One of the best-studied examples is the palindromic AT-rich repeat (PATRR), which induces recurrent constitutional translocations in humans. Here, we established a plasmid-based model that promotes frequent intermolecular rearrangements between two PATRRs in HEK293 cells. In this model system, the proportion of PATRR plasmid that extrudes a cruciform structure correlates to the levels of rearrangement. Our data suggest that PATRR-mediated translocations are attributable to unusual DNA conformations that confer a common pathway for chromosomal rearrangements in humans.

A starvation-induced noncoding RNA modulates expression of Dicer-regulated genes

Although much has been learned about short noncoding RNAs, long noncoding transcripts are largely uncharacterized. Here, we describe Caenorhabditis elegans rncs-1, a highly base-paired, 800-nucleotide noncoding RNA expressed in hypodermis and intestine. Transcription of rncs-1 is modulated in response to food supply. Although highly double-stranded, we show that rncs-1 RNA is not a substrate for Dicer because of branched structures at its termini. However, rncs-1 RNA inhibits Dicer cleavage of a second dsRNA in vitro, presumably by competition. We validate this observation in vivo by demonstrating that mRNA levels of several Dicer-regulated genes vary with changes in rncs-1 expression. Certain viruses express dsRNA to compete with cellular dsRNA-mediated pathways, and our data suggest that rncs-1 provides a cellular correlate of this phenomenon.

Evidence for active maintenance of inverted repeat structures identified by a comparative genomic approach

Inverted repeats have been found to occur in both prokaryotic and eukaryotic genomes. Usually they are short and some have important functions in various biological processes. However, long inverted repeats are rare and can cause genome instability. Analyses of C. elegans genome identified long, nearly-perfect inverted repeat sequences involving both divergently and convergently oriented homologous gene pairs and complete intergenic sequences. Comparisons with the orthologous regions from the genomes of C. briggsae and C. remanei show that the inverted repeat structures are often far more conserved than the sequences. This observation implies that there is an active mechanism for maintaining the inverted repeat nature of the sequences.

Palindrome-mediated chromosomal translocations in humans

Recently, it has emerged that palindrome-mediated genomic instability contributes to a diverse group of genomic rearrangements including translocations, deletions, and amplifications. One of the best studied examples is the recurrent t(11;22) constitutional translocation in humans that has been well documented to be mediated by palindromic AT-rich repeats (PATRRs) on chromosomes 11q23 and 22q11. De novo examples of the translocation are detected at a high frequency in sperm samples from normal healthy males, but not in lymphoblasts or fibroblasts. Cloned breakpoint sequences preferentially form a cruciform configuration in vitro. Analysis of the junction fragments implicates frequent double-strand-breaks (DSBs) at the center of both palindromic regions, followed by repair through the non-homologous end joining (NHEJ) pathway. We propose that the PATRR adopts a cruciform structure in male meiotic cells, creating genomic instability that leads to the recurrent translocation.

Palindromic AT-rich repeat in the NF1 gene is hypervariable in humans and evolutionarily conserved in primates

Palindromic sequences are dispersed in the human genome and may cause chromosomal translocations in humans. They constitute unsequenced gaps in the human genome because of their resistance to PCR amplification, cloning into vectors, and sequencing. We have overcome these difficulties by using a combination of optimized PCR conditions, cloning in a recombination-deficient E. coli strain, and RNA polymerases in sequencing. Using these methods, we analyzed a palindromic AT-rich repeat (PATRR) in the neurofibromatosis type 1 (NF1) gene on chromosome 17 (17PATRR). The 17PATRR manifests a size polymorphism due to a highly variable length of (AT)(n) dinucleotide repeats within the PATRR. 17PATRRs can be categorized into two types: a longer one that comprises a nearly or completely perfect palindrome, and a shorter one that represents its deleted asymmetric derivative. In vitro analysis shows that the longer 17PATRR is more likely to form a cruciform structure than the shorter one. Two reported t(17;22)(q11;q11) patients with NF1, whose breakpoints were identified within the 17PATRR, have translocations that are derived from perfect or nearly perfect palindromic alleles. This implies that the symmetric structure of a PATRR can induce a translocation. We identified conserved PATRRs within the NF1 gene in great apes and similar inverted repeats in two Old World monkeys, but not in New World monkeys or other mammals. This indicates that the palindromic region appeared approximately 25 million years ago and elongated during primate evolution. Although such palindromic regions are usually unstable and disappear rapidly due to deletion, the 17PATRR in the NF1 gene was stably conserved during evolution for reasons that are still unknown.

DNA replication origin plasticity and perturbed fork progression in human inverted repeats

The stability of metazoan genomes during their duplication depends on the spatiotemporal activation of origins and the progression of forks. Human rRNA genes represent a unique challenge to DNA replication since a large proportion of them exist as noncanonical palindromes in addition to canonical tandem repeats. Whether origin usage and/or fork elongation can cope with the variable structure of these genes is unknown. By analyzing single combed DNA molecules from HeLa cells, we studied the rRNA gene replication program according to the organization of canonical versus noncanonical rRNA genes. Origin positioning, spacing, and timing were not affected by the underlying rRNA gene physical structure. Conversely, fork arrest, both temporary and permanent, occurred more frequently when rRNA gene palindromes were encountered. These findings reveal that while initiation mechanisms are flexible enough to adapt to an rRNA gene structure of any arrangement, palindromes represent obstacles to fork progression, which is a likely source of genomic instability.

Palindrome content of the yeast Saccharomyces cerevisiae genome

Palindromic sequences are important DNA motifs involved in the regulation of different cellular processes, but are also a potential source of genetic instability. In order to initiate a systematic study of palindromes at the whole genome level, we developed a computer program that can identify, locate and count palindromes in a given sequence in a strictly defined way. All palindromes, defined as identical inverted repeats without spacer DNA, can be analyzed and sorted according to their size, frequency, GC content or alphabetically. This program was then used to prepare a catalog of all palindromes present in the chromosomal DNA of the yeast Saccharomyces cerevisiae. For each palindrome size, the observed palindrome counts were significantly different from those in the randomly generated equivalents of the yeast genome. However, while the short palindromes (2-12 bp) were under-represented, the palindromes longer than 12 bp were over-represented, AT-rich and preferentially located in the intergenic regions. The 44-bp palindrome found between the genes CDC53 and LYS21 on chromosome IV was the longest palindrome identified and contained only two C-G base pairs. Avoidance of coding regions was also observed for palindromes of 4-12 bp, but was less pronounced. Dinucleotide analysis indicated a strong bias against palindromic dinucleotides that could explain the observed short palindrome avoidance. We discuss some possible mechanisms that may influence the evolutionary dynamics of palindromic sequences in the yeast genome.