- Research article
- Open Access
Homologous recombination occurs frequently at innate GT microsatellites in normal somatic and germ cells in vivo
BMC Genomics volume 19, Article number: 359 (2018)
In somatic cells, homologous recombination (HR) is a rare event caused by eventual DNA double-strand breaks (DSBs). In contrast, germ cells show high frequency of HR caused by programmed DSBs. Microsatellites are prone to DSBs during genome replication and, thereby, capable of promoting HR. It remains unclear whether HR occurs frequently at microsatellites both in normal somatic cells and germ cells in a similar manner.
By examining the linkage pattern of multiple paternal and maternal markers flanking innate GT microsatellites, we measured HR at the GT microsatellites in various somatic cells and germ cells in a goldfish intraspecific heterozygote. During embryogenesis, the HR products accumulate gradually with the increase of the number of cell divisions. The frequency of HR at the GT microsatellites in advanced embryos, adult tissues and germ cells is surprisingly high. The type of exchanges between the homologous chromosomes is similar in normal advanced embryos and germ cells. Furthermore, a long GT microsatellite is more active than a short one in promoting HR in both somatic and germ cells.
HR occurs frequently at innate GT microsatellites in normal somatic cells and germ cells in a similar manner.
Homologous recombination (HR) is a type of genetic recombination in which DNA sequences are exchanged between two similar or identical DNA molecules. HR in germ cells and in somatic cells is generally considered to be markedly different in several aspects, including function, initiating lesions and outcome . HR in somatic cells is thought to be a rare event caused by eventual DNA double-strand breaks (DSBs), which frequently leads to non-crossover recombination and has potential detrimental effects due to loss of heterozygosity (LOH) [1,2,3]. In contrast, HR in germ cells is considered as a normal meiotic event initiated by programmed DSBs. This may result in crossover (CO) between homologous chromosomes and is essential for their faithful segregation at the meiosis I spindle to produce haploid gametes and increasing genetic diversity in progeny [1, 4]. Tandem-repeat DNA motifs of 1-6 nucleotides, named microsatellites or simple sequence repeats, are prone to DSBs due to replication stalling and slippage during genome replication [5,6,7,8,9,10,11], and, thereby, capable of promoting HR in the process of DNA synthesis. Since microsatellites are widely prevalent in eukaryotic genomes, it is relevant to know whether HR occurs frequently at microsatellites in both somatic and germ cells in a genome replication-dependent manner.
Repeat dinucleotides (GT)n is one of the most common microsatellites and ubiquitously present in eukaryotic genomes [12,13,14]. Studies in cell cultures have shown that GT repeats can stimulate extrachromosomal recombination of transfected DNA [15,16,17,18]. An insertion of (GT)n repeat in a yeast chromosome can significantly enhance reciprocal meiotic recombination [19, 20] and stimulate the formation of multiple crossovers . Biochemical in vitro experiments have observed that recombinase proteins, such as bacteria RecA, yeast Rad51 and human hsRad51, preferentially bind to the GT repeat sequence [21, 22]. Bioinformatic analysis of recombination patterns across human chromosome 22 and their correlation with GT repeats density suggests that innate (GT)n motif is associated with HR . Therefore, we selected innate (GT)n motifs for the in vivo study of whether HR at microsatellites has biological importance in both somatic and germ cells of vertebrates. Our results show that HR occurs frequently at innate (GT)n motifs in a similar manner in somatic cells and germ cells.
A goldfish intraspecific heterozygote provides an ideal animal system for detecting HR at GT repeats
In a previous study, we identified an evolutionary conserved imperfect (GT)n motif at the promoter of no tail (ntl) gene in different fish species . Since there are many lines and subspecies with different genetic backgrounds, goldfish ntl promoter region containing the (GT)n motif was employed for our study. To obtain different goldfish strains that contain distinguishable genetic markers, the flanking sequences of the (GT)n motif from the genome of different goldfish strains were sequenced and compared. Fortunately, a Chinese goldfish (CGF) strain, Carassius auratus auratus and a Japanese goldfish (JGF) strain, Carassius auratus cuvieri, have diversified genetic markers in the flanking sequences of the (GT)n motif. There is a 39 bp long indel in the upstream flanking sequence and two single nucleotide polymorphism (SNP) sites associated with a 3 bp long indel in the downstream flanking sequence (Fig. 1a). Moreover, there are two downstream short (GT)n motifs which are also flanked by diversified sequences in the two goldfish strains (Fig. 1a). The three (GT)n motifs in the upper, middle and downstream were named as (GT)n motif 1, 2 and 3, respectively. The sequence direction of (GT)n motif 2 is opposite to that of (GT)n motif 1 and 3 in both strains. (GT)n motif 1 is the longest of the three (GT)n motifs in both CGF and JGF strains (Fig. 1a). Length polymorphisms of (GT)n motif 1 were detected by PCR analyses in both CGF (Fig. 1b) and JGF strains (Fig. 1c). However, (GT)n motif 1 is much longer in CGF strain than in JGF strain (Fig. 1b, c).
To examine HR at GT repeats in both somatic and germ cells, we generated an intraspecific heterozygote of goldfish by crossing a homozygous CGF male individual containing an imperfect (GT)n motif 1 of about 168 bases with a homozygous JGF female individual containing an imperfect (GT)n motif 1 of about 72 bases (Fig. 1 and Additional file 1: Figure S1, S2, S3). In this heterozygote, the parental-derived sequences and different HR sequences could be identified unequivocally by multiple paternal and maternal genetic markers flanking the (GT)n motifs. The heterozygotes developed smoothly and mature heterozygous female and male mating reproduced normal offspring. Therefore, this goldfish intraspecific heterozygote provided an ideal animal system for detecting HR at GT repeats.
HR occurs frequently at innate (GT)n motif during mitosis in normal somatic cells in vivo
It is well known that HR is initiated by DSBs [1, 4], and DSBs formation in or near tandem repeats can elicit contraction and expansion of the tandem repeats during DSBs repair [7, 10, 11, 25]. We first examined whether DSBs occur at (GT)n motif 1 during embryogenesis by analyzing the length of (GT)n motif 1 at different developmental stages. PCR amplification with a pair of primers encompassing just the (GT)n motif 1 detected various contracted and expanded products in all the genome samples from heterozygous embryos and in the equal mixed parent genome samples (Fig. 1d). Though similar PCR products were amplified from the heterozygous and the mixed parental samples (probably generated by template-switching in PCR reactions), some specific contracted and expanded PCR products were obtained only from the genome of heterozygote embryos (Fig. 1d), suggesting that DSBs were formed in or near (GT)n motif 1 during embryogenesis in the heterozygotes.
Next, we examined the linkage patterns of the paternal and maternal makers flanking (GT)n motif 1 in the heterozygous embryos at different developmental stages. More than 20 clones of PCR products amplified from each testing embryonic sample were sequenced and analyzed (Fig. 2a, b). Wild type clones identical to paternal or maternal sequences, intra-allelic recombination clones with contraction or expansion of (GT)n motif 1, and HR clones with linked paternal and maternal genetic markers were detected in all the examined embryos. In the equivalent mixture of 4-day old embryos of homozygous CGF and homozygous JGF, only wild type clone identical to parental alleles and intra-allelic recombination clones with contraction or expansion of (GT)n motif 1 were detected. These results indicate that the HR clones were unlikely to be the products generated by template-switching during PCR reaction. The intra-allelic recombination clones were considered to be generated by template-switching in the PCR reactions and were eliminated in the statistical data.
In some of the HR clones, the paternal-derived sequence was at the upstream and the maternal-derived sequence was at the downstream of the (GT)n motif 1. In the rest of the HR clones, the order of paternal- and maternal-derived sequences was just the opposite (Fig. 2a). These linkage patterns of paternal and maternal makers in the HR clones indicate that reciprocal crossover between the homologous chromosomes occurs around (GT)n motif1. Surprisingly, the proportion of HR clones gradually increased along with embryogenesis and reached to as high as 33.33% in 4-day old embryos (Fig. 2b).
We further examined the frequency of HR around (GT)n motif 1 in various tissues and different individuals of the heterozygote. Frequent mitotic HR was detected in all the examined tissues of an adult heterozygous individual but the frequency was different in each of the examined tissues (Fig. 2c). Frequent mitotic HR was also detected in the genome of caudal fin of all examined adult male and female heterozygous individuals, and the frequency was somewhat different among individuals (Fig. 2d). Taken together, these results demonstrate that HR occurs frequently at (GT)n motif 1 during mitosis of normal somatic cells in vivo.
HR occurs frequently at innate (GT)n motif 1 in germ cells
We further examined HR at the (GT)n motif 1 in sperms and eggs of the heterozygote by analyzing the linkage patterns of the paternal and maternal makers flanking the (GT)n motif 1. The results showed that the frequency of HR at the (GT)n motif 1 in sperms and eggs was almost at the same level. About 33.33 and 27.78% of the sequenced clones in the sperm and eggs, respectively, of the heterozygote were HR clones (Fig. 3a, b). Interestingly, the overall proportion of HR clones in the germ cells was almost identical to that in the 4-day old embryonic cells.
Since the frequencies of HR at the innate (GT)n motif 1 were surprisingly high in both somatic and germ cells, we next examined the reliability of the HR frequency by analyzing the genotypes of the progeny from the heterozygote. Considering that, in backcross progeny, frequent HR would occur at (GT)n motif 1 between the wild type CGF and JGF homologous chromosomes derived from the hybrid and a parent of the hybrid, the true HR frequency at (GT)n motif 1 in the germ cells of hybrid could not be obtained by backcross genetic testing. To avoid confusing the HR which occurred in offspring and the HR which occurred in the parent, we used gynogenetic diploid progeny for genotype analysis. Gynogenetic diploid progeny was generated by inducing the development of matured eggs from a female heterozygote with genetic inactivated common carp sperm and then inhibiting the second polar body release of the activated eggs. In this artificial meio-gynogenetic diploid individual, the two alleles in a given locus were the duplicates derived from the same homologous chromosome. The proportion of heterozygous individuals in the gynogenetic diploid progeny of a heterozygote individual directly reflects the true frequency of meiotic HR during oogenesis. By analyzing the linkage pattern of genetic markers flanking the (GT)n motif 1 in the gynogenetic diploid progeny, we observed that 22.22% of the examined gynogenetic diploid individuals (n = 9) are heterozygotes. In these gynogenetic diploid individuals, three linkage patterns of genetic markers identical to that in the grandfather, grandmother or crossover homologs were detected (Fig. 3c). 77.78% of the examined gynogenetic individuals were homozygote identical to either grandfather or grandmother (Fig. 3c). The proportion of heterozygous individuals in the gynogenetic diploid progeny of a female heterozygoues was almost equal to proportion of meiotic HR clones detected in the eggs of heterozygote. This confirms that the high frequency of HR at the (GT)n motif 1 in the eggs of heterozygote is indeed reliable.
HR occurs independently in adjacent regions containing GT repeats in both somatic and germ cells
If the (GT)n motif 1 induces HR during mitosis and meiosis, one would expect that (GT)n motif 2 and 3 should also promote mitotic and meiotic HR, and no HR occurs at the adjacent ntl promoter region lacking GT repeats. Thus, we examined HR events in the distal ntl promoter region containing the three (GT)n motifs and the proximal ntl promoter region (from the (GT)n motif 3 to the translation start site) lacking GT repeats. In the genome of 4 day-older embryos, PCR products different from the parental-derived products in length were detected in the distal ntl promoter region (Fig. 4a), but not in the proximal region (Fig. 4b). Sequencing analysis further confirmed that HR occurs in the distal promoter region containing the three (GT)n motifs (Fig. 4c), but not in the proximal promote region lacking (GT)n motifs (data not shown).
Of the 146 sequenced clones amplified from the distal region containing the three (GT)n motifs, the percentage of paternal-derived (PD) and maternal-derived (MD) wild type clones was 58.22% (n = 42 + 43), and the percentage of HR clones was 41.78% (N = 61, Fig. 4c). The linkage patterns of paternal and maternal makers in the HR clones showed that two types of independent reciprocal crossover events, named crossover I and II, and one short fragment exchange event occurred in the region containing the three (GT)n motifs (Fig. 4c). Crossover I occurred around the (GT)n motif 1 as described above (Fig. 4c; Additional file 2: Figure S4). Crossover II occurred between the motif 2 and 3 at different locations (Fig. 4c; Additional file 2: Figure S5). The percentage of crossover I clones and crossover II clones was 21.23% (n = 31) and 14.39% (n = 21), respectively. In 6.16% of HR clones (n = 9), only a short fragment was exchanged between the paternal and maternal homologous chromosomes (Fig. 4c; Additional file 2: Figure S6), and they appeared as a non-crossover recombination event. However, the exchanged fragments were identical to the sequence regions from one of the crossover I sites to one of the crossover II sites, indicating that the exchanged fragments were produced by co-occurrence of the crossover I and II in the adjacent regions.
In germ cells, HR also occurred independently in the adjacent regions containing GT repeats in a similar manner. Independent crossover event I and II were also detected exactly at the same sites in the genome of sperm (Fig. 4d). The proportions of the three types of clones in the sperms were almost equal to that as observed in the 4-day old embryonic cells. Of the 93 sequenced clones, the proportion of paternal- and maternal-derived clones was 63.44% (n = 23 + 31), the proportion of crossover I clones, II clones and co-occurring crossover clones were 20.43% (n = 19), 9.68% (n = 9) and 6.45% (n = 6), respectively.
A longer GT repeat tract is more active than a shorter one in promoting DSBs and triggering HR in both somatic and germ cells
Since paternal (GT)n motif 1 is much longer than the maternal one, the paternal or maternal-derived (GT)n motif 1 in the HR clones can inform us how the HR at (GT)n motif 1 is carried out. To determine whether the length of innate GT repeat tract has an effect on promoting HR, the length of (GT)n motif 1 in all the sequenced HR clones from the somatic cells was analyzed. The percentage of HR clones containing different number of GT units in (GT)n motif 1 is shown in Fig. 5 according to their cell origin.
In 58.10% of mitotic HR clones (81/138), the number of GT units in (GT)n motif 1 is equal to that of maternal (GT)n motif 1. In 12.32% of mitotic HR clones (17/138), the number of GT units in (GT)n motif 1 is equal to that of paternal (GT)n motif1. The remainder of the mitotic HR clones contain (GT)n motif 1 with different contraction or expansion. The percentages of these mitotic HR clone are very low and almost at the same level. These results indicated that 58. 10% of the mitotic HR clones was induced by DSB in the paternal derived chromosome and repaired by using the maternal derived chromosome template; and 12.32% of the mitotic HR clones was induced by DSB in the maternal derived chromosome and repaired by using the paternal derived chromosome template.
The percentage of HR clones containing different number of GT units in (GT)n motif from the germ cells was similar to that observed in the somatic cells. The number of HR clones containing short maternal (GT)n motif 1 was much larger than that of HR clones containing long paternal (GT)n motif 1 (Fig. 5). Of the 63 sequenced HR clones from germ cells, 57.10% of the HR clones (n = 36) was induced by DSBs in the paternal derived chromosome and repaired by the maternal derived chromosome template, 9.52% of the HR clones (n = 6) was induced by DSB in the maternal derived chromosome and repaired by the paternal derived chromosome template. Therefore, the long (GT)n motif 1 is more susceptible to DSBs than the short one in both somatic and germ cells.
In this study, we used an in vivo system to investigate whether spontaneous DNA damage at GT microsatellites during genome replication has biological importance in promoting HR in vertebrate somatic and germ cells. By examining the linkage pattern of multiple paternal and maternal markers flanking innate GT microsatellites, we analyzed HR at the GT microsatellites in different somatic cells and germ cells in the intraspecific heterozygote. In the embryos, the HR products accumulate gradually as development proceeds. The frequency of HR at the GT microsatellites in normal advanced embryos, adult tissues and germ cells of the heterozygote is surprisingly high, and the type of exchanges between the homologous chromosomes is similar in normal advanced embryos and germ cells. Furthermore, a long GT microsatellite is more active than a short one in promoting HR in both somatic and germ cells.
The HR products detected in the germ cells may be a combination of mitotic and meiotic HR products
The gradually increases proportion of HR clones during development suggests that the incidence of HR at the (GT)n motifs correlates with the number of cell divisions. The variation of HR frequency observed among tissues also can be explained by the number of stem cell divisions during the development of the tissues. These phenomena are consistent with the characteristics of microsatellites prone to DSBs and inducing repair during genome replication in the process of cell division. Therefore, we infer that mitotic HR at the (GT)n motifs is initiated by genome replication.
HR occurs at the (GT)n motifs in the germ cells in a similar rate and produces similar recovered clone types as in the 4-day old embryos, strongly suggesting that HR at the (GT)n motifs in germ cells is also initiated by genome replication. During gametogenesis, HR induced by genome replication errors, stalling and slippage at (GT)n motifs can occur during mitosis of germ cell precursors at proliferative stage and during meiosis of germ cells. Thereby, the HR products detected at the (GT)n motifs in germ cells may be a combination of mitotic HR generated during proliferation of germ cell precursors and meiotic HR generated during meiosis of germ cells.
Frequent HR at innate (GT)n motifs may be beneficial to maintain genomic stability and suppress LOH
During DSBs repair, sister chromatids are held in proximity by cohesion whereas the homologous chromosomes may be more distant from each other, so that intact sister chromatid is readily selected as the repair template . More importantly, it has been demonstrated that RecQ helicase BLM can actively suppress crossover between homologous chromosomes, which can elicit detrimental LOH of the entire chromosome arm that is distal to the HR event, by dissolving recombination intermediates containing a double Holliday junction [26,27,28]. Therefore, mitotic HR is in general an infrequent outcome of DSBs repair , and the high frequency of HR detected at (GT)n motifs in somatic tissues is very surprising. This unusual phenomenon can be explained by the frequency of replication stalling and because slippage at (GT)n motifs is so high that DSBs occurs simultaneously at the same position on the two sister chromatids during normal mitosis. In this condition, no sister chromatid can be employed as the repair template by the DSBs repairing apparatus. Thereby, frequent HR at (GT)n motifs is an essential error correction procedure of normal mitosis for maintaining genomic stability during genome replication.
A previous study has observed that, in a yeast chromosome, insertion of a GT repeat close to a meiotic HR hotspot can inhibit the progression of the strand exchange during meiosis by stimulating two double crossovers . Our study showed that co-occurrence of independent crossovers in the adjacent regions containing GT repeats results in short segment exchange between the homologous chromosomes in both somatic and germ cells (Fig. 4c, d). It is likely that frequent replication stalling and slippage at (GT)n motifs induces simultaneously DSBs in the adjacent GT repetitive regions on each of the homologous chromosomes. DSBs initiate 5′-3′ resection of the ends. Subsequently, the 3′ overhang invades into the homologous template DNA and primes DNA synthesis to form a structure called D-loop and results in a double Holliday junction. Two double Holliday junctions formed in the adjacent region would form a four strand recombinant structure. In this case, the co-occurrence of crossovers in the adjacent GT repetitive regions would dissolve the Holliday junction of the adjacent DSBs repair intermediates and result in short segment exchange between the homologous chromosomes (Fig. 6). Therefore, frequent HR at adjacent (GT)n motifs in a chromosome would be beneficial to suppress LOH of a long chromosomal region encompassing multiple genetic loci during mitotic HR.
Frequent DSB repairing at (GT)n motifs may result in frequent bivalent formation during mitosis and play a role in preleptotene pairing during meiosis
The high frequency of mitotic HR at (GT)n motifs means that the homologous chromosomes are frequently held together in a bivalent structure by crossover during normal mitosis. The eukaryotic recombinase Rad51 orthologs, which have high affinity for GT repeat sequence [21, 22], may play an important role in initiating mitotic bivalent formation and inhibit non-homologous end joining pathway of DSB repairing during normal mitosis. Because eukaryotic organisms contain many (GT)n motifs in their genomes [12,13,14], replication errors and stalling may initiate crossovers at different GT repeat loci and hold the homologous chromosomes together in different bivalents during normal mitosis in eukaryotes. Reasonably, frequent mitotic bivalent formation mediated by (GT)n motifs would be beneficial to preventing HR between non-allelic genomic fragments that share high sequence similarity but are at different locations on the same chromosome or on different chromosomes.
DSBs formation by the ortholog of Spo11, an evolutionally conserved type II topoisomerase-like protein, appears to be a universal feature of meiotic pairing in eukaryotic organisms . However, it has been reported that in mice a significant proportion of homolog pairing is established prior to the introduction of programmed DSBs catalyzed by SPO11, and the DNA cleavage activity of SPO11 is not required for homolog pairing during meiosis . How homologous chromosomes find each other and initiate the pairing process during meiosis remains to be understood. In consideration of that, there are many (GT)n motifs scattered in the eukaryotic genome. Frequent DSBs formation at GT repeats in pre-meiotic S-phase may play an important role in preleptotene pairing of homologous chromosomes during normal meiosis.
LOH at ntl locus during embryogenesis might have a selective advantage for heterozygote
In general, LOH in somatic cells is harmful and associated with many abnormalities in mammals [1, 3, 4]. We observed that two independent mitotic HR events at (GT)n motifs led to frequent LOH at the nearby ntl locus in somatic cells during normal embryogenesis. This result implicates that a significant fraction of embryonic cells in a heterozygous individual could be converted from heterozygous to homozygous state at ntl locus during development. Gavin that (GT)n motif 1 in the promoter of ntl is conserved in different species of fish. LOH at ntl locus mediated by (GT)n motif 1 during embryogenesis may have a biological role. Tischfield has proposed that in vertebrates LOH at those loci whose proper action depends on the concentration of its gene products might have a selective advantage for the heterozygote, such as correcting mutant phenotype or promoting homozygosity of normal alleles . Indeed, high concentration of Ntl/Brachyury is required for dorsal mesoderm formation in the vertebrates and haplo-insufficiency has been detected in mammals but not in fish .
HR occurs frequently at innate GT microsatellites in a similar length-dependent manner in normal somatic and germ cells, and the incidence of HR is correlated with the number of cell divisions. Therefore, HR products detected at GT microsatellites in germ cells may include mitotic HR products generated during mitosis in germ cell precursors and meiotic HR products generated during meiosis in germ cells. Since there are many scattered (GT)n motifs in the eukaryotic genome, it would be interesting to investigate whether adjacent (GT)n motifs on the same chromosome could suppress extensive strand exchange during normal mitosis and meiosis, and whether GT microsatellites is an essential genome element for meiotic bivalent formation.
The Chinese goldfish (CGF) strain C. auratus auratus and the Japanese goldfish (JGF) strain C. auratus cuvieri, were purchased from different goldfish bred farms and maintained in our laboratory in the breeding season according to the Guidelines for the Care and Use of Laboratory Animals of Zhejiang University. The intraspecific heterozygote of goldfish was generated by crossing with a homozygous CGF male individual and a homozygous JGF female individual. Gynogenetic diploid progeny of the heterozygote was generated as described previously . The matured eggs from a female heterozygote were activated by genetic inactivated common carp sperm. Two minutes after activation, the activated eggs were cold shocked at 0-2 °C for 15 min to inhibit the second polar body release. All experimental protocols were approved by the Ethics Committee of Laboratory Animal Center of Zhejiang University (Approval no. Zju 201,306-1-11-060).
Heterozygous embryos were maintained at 22 °C and collected at desired stages. Sperms and eggs of the heterozygote were obtained from breeding heterozygous males and females, respectively, by gentle squeezing. To avoid contamination from genomes of somatic cells, the unfertilized eggs were dechlorionated with 0.2% trypsin and then washed three times with water before flash freezing. All the samples were flash frozen in liquid nitrogen until genomic DNA extraction.
Cloning ntl promoter sequence
A 2.6 kb long ntl promoter fragment upstream of the ATG was amplified from the genomic DNA of the two goldfish strains by Genome Walker™ Universal Kit (Clontech, USA), separately, and sequenced. The obtained sequences of ntl promoter have been submitted to Genbank (Accession numbers for C. carassius cuvieri: KU870661; C. carassius: KU870662).
(GT)n motif 1 stability and homologous recombination analysis
The stability of (GT)n motif 1 was tested by specific PCR with a pair of primers: 5’-TATTTTGGTTCGGACAAAGTG-3′(forward); 5’-GATGAGTCAGGGCCTGGAAA-3′ (reverse). HR at different regions of the ntl promoter was analyzed by specific PCR and sequence comparing. The sequences of primers used for amplification of the (GT)n motif1 and its flanking sequences are 5’-ACACACTTTGAACTTTCTATTTATTA-3′ (forward) and 5’-GTGTGCTCAGTGCTGGAGGGGTTA-3′ (reverse). PCR reactions was performed at 94 °C for initial denaturation for 5 min, followed by 32 cycles of 94 °C for 30 s, 55 °C for 30s, 72 °C for 1 min, finally elongated for 7 min at 72 °C. The sequences of primers used for amplification of the region covering the three (GT)n motifs and its flanking sequences are 5’-ACACACTTTGAACTTTCTATTTATTA-3′ (forward) and 5′-TTGGAGTTG ATGATAGAGGATGTTGAT-3′ (reverse). The sequences of primers used for amplification of the sequence from the (GT)n motif 3 to the transcription start site are 5’-GCGGGCTGCTCACTTTCTCACTTATAC-3′ (forward) and 5’-CGTCCAGTT TTAGTGACAATCATT-3′ (reverse). PCR reactions was performed at 94 °C for 10 s, followed 5 cycles of 94 °C for 30 s, 55 °C for 30 s, 72 °C for 2 min, and 28 cycles of 94 °C for 30 s, 50 °C for 30 s, 72 °C for 2 min,lastly elongated for 7 min at 72 °C. All PCR master mixes contained (in 50 μL): 1 × Ex-Tag Buffer (Takara, Japan), 250 μM dNTP mix, 200 nM for both primers, 10 ng-100 ng genomic DNA and 2.5 U Ex-Taq DNA polymerase (Takara, Japan).
Loss of heterozygosite
Polymerase chain reaction
Andersen SL, Sekelsky J. Meiotic versus mitotic recombination: two different routes for double-strand break repair. BioEssays. 2010;32:1058–66.
Moynahan ME, Jasin M. Mitotic homologous recombination maintains genomic stability and suppresses tumorigenesis. Nat Rev Mol Cell Biol. 2010;11:196–207.
Tischfield JA. Loss of heterozygosity or: how I learned to stop worrying and love mitotic recombination. Am J Hum Genet. 1997;61:995–9.
Lam I, Keeney S. Mechanism and regulation of meiotic recombination initiation. Cold Spring Harb Perspect Biol. 2015;7:a016634.
Gadgil R, Barthelemy J, Lewis T, Leffak M. Replication stalling and DNA microsatellite instability. Biophys Chem. 2017;225:38–48.
Strand M, Prolla TA, Liskay RM, Petes TD. Destabilization of tracts of simple repetitive DNA in yeast by mutations affecting DNA mismatch repair. Nature. 1993;365:274–6.
Paques F, Leung WY, Haber JE. Expansions and contractions in a tandem repeat induced by double-strand break repair. Mol Cell Biol. 1998;18:2045–54.
Liu G, Bissler JJ, Sinden RR, Leffak M. Unstable spinocerebellar ataxia type 10 (ATTCT*(AGAAT)) repeats are associated with aberrant replication at the ATX10 locus and replication origin-dependent expansion at an ectopic site in human cells. Mol Cell Biol. 2007;27:7828–38.
Hardwick RJ, Tretyakov MV, Dubrova YE. Age-related accumulation of mutations supports a replication-dependent mechanism of spontaneous mutation at tandem repeat DNA loci in mice. Mol Biol Evol. 2009;26:2647–54.
Jankowski C, Nasar F, Nag DK. Meiotic instability of CAG repeat tracts occurs by double-strand break repair in yeast. Proc Natl Acad Sci U S A. 2000;97:2134–9.
Wojcik EA, Brzostek A, Bacolla A, Mackiewicz P, Vasquez KM, Korycka-Machala M, et al. Direct and inverted repeats elicit genetic instability by both exploiting and eluding DNA bouble-strand break repair systems in mycobacteria. PLoS One. 2012;7(12):e51064.
Miesfeld R, Krystal M, Arnheim N. A member of a new repeated sequence family which is conserved throughout eucaryotic evolution is found between the human delta and beta globin genes. Nucleic Acids Res. 1981;9:5931–47.
Hamada H, Kakunaga T. Potential Z-DNA forming sequences are highly dispersed in the human genome. Nature. 1982;298:396–8.
Hamada H, Petrino MG, Kakunaga T. A novel repeated element with Z-DNA forming potential is widely found in evolutionarily diverse eukaryotic genomes. Proc Natl Acad Sci U S A. 1982;79:6465–9.
Stringer JR. Recombination between poly [d(GT)d(CA)] sequences in simian virus 40-infected cultured cells. Mol Cell Biol. 1985;5:1247–59.
Bullock P, Miller J, Botchan M. Effects ofpoly [d(pGpT)d(pApC)] and poly [d(pCpG)d(pCpG)] repeatson homologous recombination in somatic cells. Mol Cell Biol. 1986;6:3948–53.
Wahls WP, Wallace LJ, Moore PD. The Z-DNA motif d(TG)30 promotes reception of information during gene conversion events while stimulating homologous recombination in human cells in culture. Mol Cell Biol. 1990;10:785–93.
Wahls WP, Moore PD. Homologous recombination enhancement conferred by the Z-DNA motif d(TG)30 is abrogated by simian virus 40 T antigen binding to adjacent DNA sequences. Mol Cell Biol. 1990;10:794–800.
Treco D, Arnheim N. The evolutionarily conserved repetitive sequence d(TG.AC) promotes reciprocal exchange and generates unusual recombinant tetrads during yeast meiosis. Mol Cell Biol. 1986;6:3934–47.
Grendrel C, Boulet A, Dutreix M. (CA/GT)n microsatellites affect homologous recombination during yeast meiosis. Genes Dev. 2000;14:1261–8.
Tracy RB, Baumohl JK, Kowalczykowski SC. The preference for GT-rich DNA by the yeast Rad51 protein defines a set of universal pairing sequences. Genes Dev. 1997;11:3423–31.
Biet E, Sun JS, Dutreix M. Conserved sequence preference in DNA binding among recombination proteins: an effect of ssDNA secondary structure. Nucleic Acids Res. 1999;27:596–600.
Majewski J, Ott J. GT repeats are associated with recombinationon human chromosome 22. Genome Res. 2000;10:1108–14.
Huang W, Zheng J, He Y, Luo C. Tandem repeat modification during double strand break repair induced by an engineered TAL effector nuclease in zebrafish genome. PLoS One. 2013;8:e84176.
Pollard LM, Bourn RL, Bidichandani SI. Repair of DNA double-strand breaks within the (GAATTC)n sequence results in frequent deletion of the triplet-repeat sequence. Nucleic Acids Res. 2008;36:489–500.
Wu L, Hickson ID. The Bloom’s syndrome helicase suppresses crossing over during homologous recombination. Nature. 2003;426:870–4.
Stark JM, Jasin M. Extensive loss of heterozygosity is suppressed during homologous repair of chromosomal breaks. Mol Cell Biol. 2003;23:733–43.
LaRocquea JR, Starka JM, Oh J, Bojilova E, Yusac K, Horiec K, et al. Interhomolog recombination and loss of heterozygosity in wild-type and bloom syndrome helicase (BLM)-deficient mammalian cells. Proc Natl Acad Sci U S A. 2011;108:11971–6.
Boateng KA, Bellani MA, Gregoretti I, Pratto F, Camerini-Otero RD. Homologous pairing preceding SPO11-mediated double-strand breaks in mice. Dev Cell. 2013;24:196–205.
Wardle FC, Papaioannou VE. Teasing out T-box targets in early mesoderm. Curr Opin Genet Dev. 2008;18:418–25.
Ma S, Huang W, Zhang L, Zhao S, Tong Y, Liu Z, et al. Germ cell-specific DNA methylation and genome diploidization in primitive vertebrates. Epigenetics. 2011;6:1471–80.
The authors would like to thank Mr. Jifu Shen of our laboratory of Zhejiang University for the assistance in heterozygote goldfish culture.
This work was supported by Development Program for Basic Research of China [2010CB126301] and National Natural Science Foundation of China . This funding body had no role in the design of the study and collection, analysis, and interpretation of data, or in writing the manuscript.
Availability of data and materials
The datasets supporting the results of this article are available in the National Center for Biotechnology Information (NCBI) Sequence Read Archive under accession number for C. carassius cuvieri: KU870661; for C. carassius: KU870662.
Ethics approval and consent to participate
All goldfish and zebrafish experiments were approved by the Ethics Committee of Laboratory Animal Center of Zhejiang University (Approval no. Zju201306-1-11-060) and all experiments were carried out in accordance with the approved protocols.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Genome sequences of the paternal and maternal ntl promoter region containing the (GT)n motifs. Figure S1. Alignment of paternal and maternal goldfish distal ntl promoter region containing the (GT)n motifs. Figure S2. Sequencing spectrums of the paternal distal ntl promoter region containing the (GT)n motifs. Figure S3. Sequencing spectrums of the maternal distal ntl promoter region containing the (GT)n motifs. (DOCX 1620 kb)
Genome sequences of various reciprocal crossovers clones at different sites of ntl promoter. Figure S4. Sequencing spectrums of reciprocal crossover I clones in ntl promoter. Figure S5. Sequencing spectrums of reciprocal crossover II clones in ntl promoter. Figure S6. Sequencing spectrums of HR clones that only a fragment is exchanged between the paternal and maternal homologous chromosomes. (DOCX 3959 kb)