Microhomology-mediated intron loss during metazoan evolution

RNA transcripts serve as a template for double-strand break repair in human cells

Double-strand breaks (DSBs) are toxic lesions that lead to genome instability. While canonical DSB repair pathways typically operate independently of RNA, growing evidence suggests that RNA:DNA hybrids and nearby transcripts can influence repair outcomes. However, whether transcript RNA can directly serve as a template for DSB repair in human cells remains unclear. In this study, we develop fluorescence and sequencing-based assays to show that RNA-containing oligonucleotides and messenger RNA can serve as templates during DSB repair. We conduct a CRISPR/Cas9-based genetic screen to identify factors that promote RNA-templated DSB repair (RT-DSBR). Of the candidate polymerases, we identify DNA polymerase zeta (Polζ) as a potential reverse transcriptase that facilitates RT-DSBR. Furthermore, analysis of cancer genome sequencing data reveals whole intron deletions - a distinct genomic signature of RT-DSBR that occurs when spliced mRNA guides repair. Altogether, our findings highlight RT-DSBR as an alternative pathway for repairing DSBs in transcribed genes, with potential mutagenic consequences.

RNA Transcripts Serve as a Template for Double-Strand Break Repair in Human Cells

Double-strand breaks (DSBs) are toxic lesions that lead to genome instability. While canonical DSB repair pathways typically operate independently of RNA, emerging evidence suggests that RNA:DNA hybrids and transcripts near damaged sites can influence repair outcomes. However, a direct role for transcript RNA as a template during DSB repair in human cells is yet to be established. In this study, we designed fluorescent- and sequencing-based assays, which demonstrated that RNA-containing oligonucleotides and messenger RNA serve as templates to promote DSB repair. We conducted a CRISPR/Cas9-based genetic screen to identify factors that promote RNA-templated DSB repair (RT-DSBR), and of the candidate polymerases, we identified DNA polymerase-zeta (Polζ) as the potential reverse transcriptase that facilitates RT-DSBR. Furthermore, by analyzing sequencing data from cancer genomes, we identified the presence of whole intron deletions, a unique genomic scar reflective of RT-DSBR activity generated when spliced mRNA serves as the repair template. These findings highlight RT-DSBR as an alternative pathway for repairing DSBs in transcribed genes, with potential mutagenic consequences.

Intron losses and gains in the nematodes

BACKGROUND

The evolution of spliceosomal introns has been widely studied among various eukaryotic groups. Researchers nearly reached the consensuses on the pattern and the mechanisms of intron losses and gains across eukaryotes. However, according to previous studies that analyzed a few genes or genomes, Nematoda seems to be an eccentric group.

RESULTS

Taking advantage of the recent accumulation of sequenced genomes, we extensively analyzed the intron losses and gains using 104 nematode genomes across all the five Clades of the phylum. Nematodes have a wide range of intron density, from less than one to more than nine per kbp coding sequence. The rates of intron losses and gains exhibit significant heterogeneity both across different nematode lineages and across different evolutionary stages of the same lineage. The frequency of intron losses far exceeds that of intron gains. Five pieces of evidence supporting the model of cDNA-mediated intron loss have been observed in ten Caenorhabditis species, the dominance of the precise intron losses, frequent loss of adjacent introns, high-level expression of the intron-lost genes, preferential losses of short introns, and the preferential losses of introns close to 3'-ends of genes. Like studies in most eukaryotic groups, we cannot find the source sequences for the limited number of intron gains detected in the Caenorhabditis genomes.

CONCLUSIONS

These results indicate that nematodes are a typical eukaryotic group rather than an outlier in intron evolution.

Different DNA repair pathways are involved in single-strand break-induced genomic changes in plants

In nature, single-strand breaks (SSBs) in DNA occur more frequently (by orders of magnitude) than double-strand breaks (DSBs). SSBs induced by the CRISPR/Cas9 nickase at a distance of 50-100 bp on opposite strands are highly mutagenic, leading to insertions/deletions (InDels), with insertions mainly occurring as direct tandem duplications. As short tandem repeats are overrepresented in plant genomes, this mechanism seems to be important for genome evolution. We investigated the distance at which paired 5'-overhanging SSBs are mutagenic and which DNA repair pathways are essential for insertion formation in Arabidopsis thaliana. We were able to detect InDel formation up to a distance of 250 bp, although with much reduced efficiency. Surprisingly, the loss of the classical nonhomologous end joining (NHEJ) pathway factors KU70 or DNA ligase 4 completely abolished tandem repeat formation. The microhomology-mediated NHEJ factor POLQ was required only for patch-like insertions, which are well-known from DSB repair as templated insertions from ectopic sites. As SSBs can also be repaired using homology, we furthermore asked whether the classical homologous recombination (HR) pathway is involved in this process in plants. The fact that RAD54 is not required for homology-mediated SSB repair demonstrates that the mechanisms for DSB- and SSB-induced HR differ in plants.

The multicomparative 2-n-way genome suite

To effectively analyze the increasing amounts of available genomic data, improved comparative analytical tools that are accessible to and applicable by a broad scientific community are essential. We built the “2-n-way” software suite to provide a fundamental and innovative processing framework for revealing and comparing inserted elements among various genomes. The suite comprises two user-friendly web-based modules. The 2-way module generates pairwise whole-genome alignments of target and query species. The resulting genome coordinates of blocks (matching sequences) and gaps (missing sequences) from multiple 2-ways are then transferred to the n-way module and sorted into projects, in which user-defined coordinates from reference species are projected to the block/gap coordinates of orthologous loci in query species to provide comparative information about presence (blocks) or absence (gaps) patterns of targeted elements over many entire genomes and phylogroups. Thus, the 2-n-way software suite is ideal for performing multidirectional, non-ascertainment-biased screenings to extract all possible presence/absence data of user-relevant elements in orthologous sequences. To highlight its applicability and versatility, we used 2-n-way to expose approximately 100 lost introns in vertebrates, analyzed thousands of potential phylogenetically informative bat and whale retrotransposons, and novel human exons as well as thousands of human polymorphic retrotransposons.

Evolutionary Dynamics of the SKN-1 → MED → END-1,3 Regulatory Gene Cascade in Caenorhabditis Endoderm Specification

Gene regulatory networks and their evolution are important in the study of animal development. In the nematode, Caenorhabditis elegans, the endoderm (gut) is generated from a single embryonic precursor, E. Gut is specified by the maternal factor SKN-1, which activates the MED → END-1,3 → ELT-2,7 cascade of GATA transcription factors. In this work, genome sequences from over two dozen species within the Caenorhabditis genus are used to identify MED and END-1,3 orthologs. Predictions are validated by comparison of gene structure, protein conservation, and putative cis-regulatory sites. All three factors occur together, but only within the Elegans supergroup, suggesting they originated at its base. The MED factors are the most diverse and exhibit an unexpectedly extensive gene amplification. In contrast, the highly conserved END-1 orthologs are unique in nearly all species and share extended regions of conservation. The END-1,3 proteins share a region upstream of their zinc finger and an unusual amino-terminal poly-serine domain exhibiting high codon bias. Compared with END-1, the END-3 proteins are otherwise less conserved as a group and are typically found as paralogous duplicates. Hence, all three factors are under different evolutionary constraints. Promoter comparisons identify motifs that suggest the SKN-1, MED, and END factors function in a similar gut specification network across the Elegans supergroup that has been conserved for tens of millions of years. A model is proposed to account for the rapid origin of this essential kernel in the gut specification network, by the upstream intercalation of duplicate genes into a simpler ancestral network.

Parallel loss of introns in the ABCB1 gene in angiosperms

BACKGROUND

The presence of non-coding introns is a characteristic feature of most eukaryotic genes. While the size of the introns, number of introns per gene and the number of intron-containing genes can vary greatly between sequenced eukaryotic genomes, the structure of a gene with reference to intron presence and positions is typically conserved in closely related species. Unexpectedly, the ABCB1 (ATP-Binding Cassette Subfamily B Member 1) gene which encodes a P-glycoprotein and underlies dwarfing traits in maize (br2), sorghum (dw3) and pearl millet (d2) displayed considerable variation in intron composition.

RESULTS

An analysis of the ABCB1 gene structure in 80 angiosperms revealed that the number of introns ranged from one to nine. All introns in ABCB1 underwent either a one-time loss (single loss in one lineage/species) or multiple independent losses (parallel loss in two or more lineages/species) with the majority of losses occurring within the grass family. In contrast, the structure of the closest homolog to ABCB1, ABCB19, remained constant in the majority of angiosperms analyzed. Using known phylogenetic relationships within the grasses, we determined the ancestral branch-points where the losses occurred. Intron 7, the longest intron, was lost in only a single species, Mimulus guttatus, following duplication of ABCB1. Semiquantitative PCR showed that the M. guttatus ABCB1 gene copy without intron 7 had significantly lower transcript levels than the gene copy with intron 7. We further demonstrated that intron 7 carried two motifs that were highly conserved across the monocot-dicot divide.

CONCLUSIONS

The ABCB1 gene structure is highly dynamic, while the structure of ABCB19 remained largely conserved through evolution. Precise removal of introns, preferential removal of smaller introns and presence of at least 2 bp of microhomology flanking most introns indicated that intron loss may have predominantly occurred through non-homologous end-joining (NHEJ) repair of double strand breaks. Lack of microhomology in the exon upstream of lost phase I introns was likely due to release of the selective constraint on the penultimate base (3rd base in codon) of the terminal codon by the splicing machinery. In addition to size, the presence of regulatory motifs will make introns recalcitrant to loss.

Regulation of non-homologous end joining via post-translational modifications of components of the ligation step

DNA double-strand breaks are the most serious type of DNA damage and non-homologous end joining (NHEJ) is an important pathway for their repair. In Saccharomyces cerevisiae, three complexes mediate the canonical NHEJ pathway, Ku (Ku70/Ku80), MRX (Mre11/Rad50/Xrs2) and DNA ligase IV (Dnl4/Lif1). Mammalian NHEJ is more complex, primarily as a consequence of the fact that more factors are involved in the process, and also because higher chromatin organization and more complex regulatory networks exist in mammals. In addition, a stronger interconnection between the NHEJ and DNA damage response (DDR) pathways seems to occur in mammals compared to yeast. DDR employs multiple post-translational modifications (PTMs) of the target proteins and mutual crosstalk among them to ensure highly efficient down-stream effects. Checkpoint-mediated phosphorylation is the best understood PTM that regulates DDR, although recently SUMOylation has also been shown to be involved. Both phosphorylation and SUMOylation affect components of NHEJ. In this review, we discuss a role of these two PTMs in regulation of NHEJ via targeting the components of the ligation step.

Evaluation of the mechanisms of intron loss and gain in the social amoebae Dictyostelium

Background

Spliceosomal introns are a common feature of eukaryotic genomes. To approach a comprehensive understanding of intron evolution on Earth, studies should look beyond repeatedly studied groups such as animals, plants, and fungi. The slime mold Dictyostelium belongs to a supergroup of eukaryotes not covered in previous studies.

Results

We found 441 precise intron losses in Dictyostelium discoideum and 202 precise intron losses in Dictyostelium purpureum. Consistent with these observations, Dictyostelium discoideum was found to have significantly more copies of reverse transcriptase genes than Dictyostelium purpureum. We also found that the lost introns are significantly further from the 5′ end of genes than the conserved introns. Adjacent introns were prone to be lost simultaneously in Dictyostelium discoideum. In both Dictyostelium species, the exonic sequences flanking lost introns were found to have a significantly higher GC content than those flanking conserved introns. Together, these observations support a reverse-transcription model of intron loss in which intron losses were caused by gene conversion between genomic DNA and cDNA reverse transcribed from mature mRNA. We also identified two imprecise intron losses in Dictyostelium discoideum that may have resulted from genomic deletions. Ninety-eight putative intron gains were also observed. Consistent with previous studies of other lineages, the source sequences were found in only a small number of cases, with only two instances of intron gain identified in Dictyostelium discoideum.

Conclusions

Although they diverged very early from animals and fungi, Dictyostelium species have similar mechanisms of intron loss.

Electronic supplementary material

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

Microhomology-Mediated End Joining: A Back-up Survival Mechanism or Dedicated Pathway?

DNA double-strand breaks (DSBs) disrupt the continuity of chromosomes and their repair by error-free mechanisms is essential to preserve genome integrity. Microhomology-mediated end joining (MMEJ) is an error-prone repair mechanism that involves alignment of microhomologous sequences internal to the broken ends before joining, and is associated with deletions and insertions that mark the original break site, as well as chromosome translocations. Whether MMEJ has a physiological role or is simply a back-up repair mechanism is a matter of debate. Here we review recent findings pertaining to the mechanism of MMEJ and discuss its role in normal and cancer cells.

Intron evolution in Saccharomycetaceae

Introns in protein-coding genes are very rare in hemiascomycetous yeast genomes. It has been suggested that these species have experienced extensive intron loss during their evolution from the postulated intron-rich fungal ancestor. However, no intron-devoid yeast species have been identified and some of the introns remaining within the genomes of intron-poor species, such as Saccharomyces cerevisiae, appear to be beneficial during growth under stress conditions. In order to reveal the pattern of intron retention within intron-poor yeast species and better understand the mechanisms of intron evolution, we generated a comprehensive set of 250 orthologous introns in the 20 species that comprise the Saccharomycetaceae, by analyzing RNA deep-sequencing data and alignments of intron-containing genes. Analysis of these intron sets shows that intron loss is at least two orders of magnitude more frequent than intron gain. Fine mapping of intron positions shows that intron sliding is rare, and that introns are almost always removed without changing the primary sequence of the encoded protein. The latter finding is consistent with the prevailing view that homologous recombination between reverse-transcribed mature mRNAs and the corresponding genomic locus is the primary mechanism of intron loss. However, we also find evidence that loss of a small number of introns is mediated by micro-homology, and that the number of intron losses is diminished in yeast species that have lost the microhomology end joining and nonhomologous end joining machinery.

The dynamic loss and gain of introns during the evolution of the Brassicaceae

Sequence comparison allows the detailed analysis of evolution at the nucleotide and amino acid levels, but much less information is known about the structural evolution of genes, i.e. how the number, length and distribution of introns change over time. We constructed a parsimonious model for the evolutionary rate of intron loss (IL) and intron gain (IG) within the Brassicaceae and found that IL/IG has been highly dynamic, with substantial differences between and even within lineages. The divergence of the Brassicaceae lineages I and II marked a dramatic change in the IL rate, with the common ancestor of lineage I losing introns three times more rapidly than the common ancestor of lineage II. Our data also indicate a subsequent declining trend in the rate of IL, although in Arabidopsis thaliana introns continue to be lost at approximately the ancestral rate. Variations in the rate of IL/IG within lineage II have been even more remarkable. Brassica rapa appears to have lost introns approximately 15 times more rapidly than the common ancestor of B. rapa and Schenkiella parvula, and approximately 25 times more rapidly than its sister species Eutrema salsugineum. Microhomology was detected at the splice sites of several dynamic introns suggesting that the non-homologous end-joining and double-strand break repair is a common pathway underlying IL/IG in these species. We also detected molecular signatures typical of mRNA-mediated IL, but only in B. rapa.

Higher frequency of intron loss from the promoter proximally paused genes of Drosophila melanogaster

Although intron losses have been widely reported, it is not clear whether they are neutral and therefore random or driven by positive selection. Intron transcription and splicing are time-consuming and can delay the expression of its host gene. For genes that must be activated quickly to respond to physiological or stress signals, intron delay may be deleterious. Promoter proximally paused (PPP) genes are a group of rapidly expressed genes. To respond quickly to activation signals, they generally initiate transcription competently but stall after synthesizing a short RNA. In this study, performed in Drosophila melanogaster, the PPP genes were found to have a significantly higher rate of intron loss than control genes. However, further analysis did not find more significant shrinkage of intron size in PPP genes. Referring to previous studies on the rates of transcription and splicing and to the time saved by deletion of the introns from mouse gene Hes7, it is here suggested that transcription delay is comparable to splicing delay only when the intron is 28.5 kb or larger, which is greater in size than 95% of vertebrate introns, 99.5% of Drosophila introns, and all the annotated introns of Saccharomyces cerevisiae and Arabidopsis thaliana. Delays in intron splicing are probably a selective force, promoting intron loss from quickly expressed genes. In other genes, it may have been an exaptation during the emergency of developmental clocks.

Alternative end-joining pathway(s): bricolage at DNA breaks

To cope with DNA double strand break (DSB) genotoxicity, cells have evolved two main repair pathways: homologous recombination which uses homologous DNA sequences as repair templates, and non-homologous Ku-dependent end-joining involving direct sealing of DSB ends by DNA ligase IV (Lig4). During the last two decades a third player most commonly named alternative end-joining (A-EJ) has emerged, which is defined as any Ku- or Lig4-independent end-joining process. A-EJ increasingly appears as a highly error-prone bricolage on DSBs and despite expanding exploration, it still escapes full characterization. In the present review, we discuss the mechanism and regulation of A-EJ as well as its biological relevance under physiological and pathological situations, with a particular emphasis on chromosomal instability and cancer. Whether or not it is a genuine DSB repair pathway, A-EJ is emerging as an important cellular process and understanding A-EJ will certainly be a major challenge for the coming years.