Recombination is an engine of genetic diversity and therefore constitutes a key process in evolutionary biology and genetics. While the outcome of crossover recombination can readily be detected as shuffled alleles by following the inheritance of markers in pedigreed families, the more precise location of both crossover and non-crossover recombination events has been difficult to pinpoint. As a consequence, we lack a detailed portrait of the recombination landscape for most organisms and knowledge on how this landscape impacts on sequence evolution at a local scale. To localize recombination events with high resolution in an avian system, we performed whole-genome re-sequencing at high coverage of a complete three-generation collared flycatcher pedigree. We identified 325 crossovers at a median resolution of 1. 4 kb, with 86% of the events localized to <10 kb intervals. Observed crossover rates were in excellent agreement with data from linkage mapping, were 52% higher in male (3. 56 cM/Mb) than in female meiosis (2. 28 cM/Mb), and increased towards chromosome ends in male but not female meiosis. Crossover events were non-randomly distributed in the genome with several distinct hot-spots and a concentration to genic regions, with the highest density in promoters and CpG islands. We further identified 267 non-crossovers, whose location was significantly associated with crossover locations. We detected a significant transmission bias (0. 18) in favour of ‘strong’ (G, C) over ‘weak’ (A, T) alleles at non-crossover events, providing direct evidence for the process of GC-biased gene conversion in an avian system. The approach taken in this study should be applicable to any species and would thereby help to provide a more comprehensive portray of the recombination landscape across organism groups.
Meiotic recombination is intimately related to the evolution of sexual reproduction. It occurs early in meiosis and is commonly initiated by double-strand breaks (DSBs) that are catalysed by the SPO11 protein [1]. The broken ends are processed and their repair can either lead to crossovers (COs), which involve an exchange of chromatid arms and assist the proper segregation of homologous chromosomes during meiosis I, or non-crossovers (NCOs), i. e. recombination events without an exchange of chromatid arms. We here use the term ‘recombination’ to collectively refer to both types of events, while CO and NCO are used to refer to the respective outcome of recombination. COs are critical to several evolutionary processes [2], such as the efficacy of selection (Hill-Robertson interference; [3]), and the evolution of sex chromosomes [4]. COs further modulate variation in levels of nucleotide diversity along chromosomes [5–8] and genetic differentiation between populations and species [9], and, together with selection, govern the character and extent of linkage disequilibrium [10]. Moreover, in addition to breaking up linkage and re-shuffling alleles, recombination affects the evolution of base composition via GC-biased gene conversion (gBGC) [2,11–13]. gBGC is a process that leads to a preferential transmission of GC-alleles over AT-alleles close to recombination-initiating DSBs. Base pair mismatches establish in heteroduplex DNA, which is formed as part of the repair pathway of DSBs, whenever homologous chromosomes carry different alleles. The transmission bias arises because mismatches that result from AT/GC heterozygous sites are resolved in favour of G: C base pairs. CO rates (often measured as cM/Mb) can be estimated by combining data on CO fractions between markers in linkage analyses and physical information on the location of markers in the genome. Typically, resolution is limited by the density of available markers for genotyping and the number of meiosis in which the segregation of markers from parents to offspring can be followed. Recent development of arrays with tens or even hundreds of thousands of single nucleotide polymorphism (SNP) markers have offered increased resolution [14–16] but the number of genotypes required to identify CO events between closely located markers still represents a limiting factor for fine-scale assessment of CO rates in most non-model organisms. Yet, comparisons of linkage maps and genome sequences [17,18] have improved our understanding of the broad-scale patterns of CO rate variation concerning, for example, rate differences between species [19–21], chromosomes [22], and sexes [23], as well as regional heterogeneity along chromosomes [24,25]. The application of whole-genome re-sequencing to population genomics provides an indirect means to the estimation of fine-scale CO rates and can allow localization of historical CO events [26]. Specifically, estimated levels of linkage disequilibrium (LD) between pairs of segregating sites along chromosomes can be transformed into the scaled population recombination parameter (ρ = 4Ner), which can be used as a proxy for CO rate; high levels of LD are indicative of low CO rates, while low levels of LD are most easily explained by a high rate of COs. However, a drawback of this approach is that LD can be influenced by other forces than CO, such as selection, population structure and migration. Moreover, patterns of LD are the result of historical processes and do not necessarily reflect the properties of contemporary CO [27]. Whole-genome re-sequencing can also be used to get direct estimates of recombination rates. Specifically, re-sequencing of crosses [28–31], pedigrees [32,33], sperm and oocytes [34–37] or spores [28] provide new and exciting direct approaches for localization of recombination events at high resolution and for the estimation of recombination rates. In principle, the density of informative polymorphisms that distinguish homologous chromosomes determines the resolution. Sequencing of gametes provides the structure of new haplotypes formed after CO events and meiotic tetrad analysis is particularly attractive in this respect since all four products from a single meiosis can be recovered and characterized [38]. It allows not only identifying CO events at high resolution but 3: 1 inheritance between sister gametes implies that NCO gene conversion events can also be traced [28]. In organisms in which tetrad analysis is not possible, phased sequencing data from pedigreed individuals, most easily obtained if three generations can be followed, is technically less demanding to generate than comparable data from single-cell analysis. One important conclusion from studies on recombination rate variation at high resolution is the realization that recombination events are often concentrated to specific genomic regions, so-called hot-spots. These are likely to coincide with regions accessible for DSB formation. In budding yeast (Saccharomyces cerevisiae), nucleosome occupancy and the histone H3 lysine 4 trimethylation (H3K4me3) chromatin modification facilitate the formation of DSBs by changing the accessibility for SPO11, and high rates of DSBs are observed in close proximity to transcription start sites (TSSs) [39,40]. A similar picture is seen in plants; high rates of recombination are observed in close proximity of TSSs but also in close proximity of transcription termination sites (TTSs) [41–43]. In humans and some other mammals (but not all [44]), the localization of hot-spots is associated with certain sequence motifs that are recognized by PRDM9 [45–48], a zinc finger protein that trimethylates H3K4me3 [49]. In this case, PRDM9 binding occurs mainly in intergenic regions, and within genes, with a lowered rate of recombination close to TSSs [50–52]. While the co-localization of recombination and transcription initiation seems to be a widespread and likely ancestral mechanism, the PRDM9-directed recombination is apparently a derived character with limited phylogenetic distribution [53,54]. Nevertheless, a common feature of the localization of meiotic recombination events across species seems to be the influence of chromatin structure. Birds have high CO rates compared to the mammalian sister lineage and also show high within-genome variation in the rate of COs [22,54–57]. The former owes to the fact that avian karyotypes are characterized by a large number of chromosomes and that there is a positive correlation between the amount of COs and the number of chromosomes across organisms [58]). The latter is a consequence of significant variation in chromosome size with numerous small microchromosomes (<5–10 Mb) in which one obligate CO event per chromosome [59] implies high rates of COs per physical unit DNA. However, not much is known about rate heterogeneity at a local scale and what determines the genomic location of CO as well as NCO events in this vertebrate lineage [54]. Addressing these issues is particularly warranted by the fact that birds lack Prdm9 [60], raising the question if the regulation of recombination is more similar to what might be the ancestral mechanism, found in yeast and plants, than to the mechanism found in the mammalian sister lineage [54]. To gain increased insight into recombination in an avian system we localized and characterized recombination events with high resolution by whole-genome re-sequencing of a three-generation pedigree of the collared flycatcher, Ficedula albicollis. We thereby benefitted from the access to a genome assembly with high sequence continuity and with scaffolds anchored, ordered and oriented on chromosomes [55,61]. Moreover, the availability of a high-density linkage map in this species [55] provides valuable background information on regional CO rate variation across the genome. We identified 325 CO events (at a median resolution of 1. 4 kb and with 86% of the events localized to regions < 10 kb), as well as 267 NCO gene conversion events, and used these data to analyse the characteristics and consequences of recombination in an avian system. Our main conclusions from this work are that there is a concentration of recombination events to certain hot-spot regions, which show an association with genes, especially promotor regions and CpG islands. We further find that CO rates are 52% higher in male (3. 56 cM/Mb) than in female meiosis (2. 28 cM/Mb), that the male CO rate is higher towards chromosome ends, and that there is positive CO interference up to a distance of 14 Mb. The location of NCO events is associated with the location of CO events, while no significant difference between sexes can be observed. Moreover, we find a significant transmission distortion in favour of G and C alleles over A and T alleles at NCOs, providing direct evidence for GC-biased gene conversion in an avian species.
We performed whole-genome re-sequencing (mean autosomal coverage = 42X, range 36. 9–45. 4X; S1 Table) of 11 collared flycatchers from a three-generation pedigree (Fig 1A) in which 4. 434 million segregating SNPs originating from the four grandparents and being informative for phasing (Fig 1B) were identified. A total of 325 meiotic CO events (50–67 per offspring, positions given in S2 Table) were identified in the transmission of gametes from the two F1 parents to the five F2 offspring by mapping the transitions between haploblocks along chromosomes (Fig 1C). Due to a high degree of nucleotide diversity (π) in the population (mean π = 3. 6 x 10−3; [9,61]) and that deep sequencing allowed SNPs to be called at a high rate, the position of CO events could be identified with high accuracy (S1 Fig). The median interval between recombinant SNP markers was 1,513 bp, or 1,360 bp if only considering events in genomic regions without assembly gaps. Eighty per cent of all CO events could be mapped with a resolution of <5 kb, and 86% <10 kb, with similar resolution in all five F2 offspring (S3 Table). After very stringent filtering (see Methods) we further identified a total of 267 NCO gene conversions spread across the flycatcher genome (S4 Table). Given the stringent filtering and that the power to detect NCOs is low, the set of identified NCOs likely represents only a subset of all such events. The number of CO events per chromosome and meiosis ranged between 0–6. The amount of COs per chromosome as reflected in number of CO events was in excellent agreement with predictions from genetic distances observed in linkage analysis (Pearson’s r = 0. 95, Table 1), providing strong support for the overall accuracy of the detection of CO events. Moreover, regional (200 kb windows) CO rate estimates based on linkage analysis are available for the collared flycatcher genome [55] and windows corresponding to the location of CO events detected in the present study had a significantly higher linkage-based CO rate (6. 27 cM/Mb) than windows without CO events (3. 65 cM/Mb; t-test, p = 6. 3 x 10−7). The total amount of observed COs per meiosis corresponded to a sex-averaged autosomal genetic distance of 3,030 cM, very close to that obtained in linkage analysis (3,067 cM; [55]). The average CO rate in autosomes (data available for 30 autosomes) was 3. 08 cM/Mb but the rate was highly variable among chromosomes and showed a strong non-linear correlation with chromosome size (Fig 2). For chromosomes >50 Mb, the mean CO rate was 1. 98 cM/Mb and for chromosomes <10 Mb the mean rate was 13. 02 cM/Mb. Of 305 detected autosomal CO events, 119 were of maternal origin and 186 of paternal origin. This corresponds to total map distances of 2,380 cM in female meiosis (2. 28 cM/Mb) and 3,720 cM (3. 56 cM/Mb) in male meiosis, i. e. 56% higher CO rate in males than in males. Out of a total of 20 events on the Z chromosome, two were of maternal origin and were located in the ≈ 0. 6 Mb pseudoautosomal region [62]. This confirms a very high rate (67 cM/Mb, similar to what has been observed in linkage analysis; [62]) of COs in this short region, which corresponds to ≈1% of the Z chromosome and which is the only region where the Z chromosome and W chromosome pairs in female meiosis. Contrary to the pattern observed for COs, the average number of NCOs that occurred during male meiosis (25. 4) was not statistically different to the number of NCOs that occurred during female meiosis (27. 4; t-test, p = 0. 506). Only a single NCO event was identified on the Z chromosome and was of paternal origin. It is often assumed that one obligate CO per chromosome is necessary for proper segregation during meiosis, irrespective of chromosome size [63] (though there are organisms in which this does not apply, like the absence of recombination in male meiosis of Drosophila). We observed many instances of transmitted chromosomes without a detected CO event. This is not surprising given that 50% of the gametes from CO events will be non-recombinants. Thus, when there is close to only one CO event per chromosome on average, there should be about as many gametes with an observable CO event as without. To corroborate this assumption we focused on maternal transmission of the smallest microchromosomes (as the overall rate of COs was lower in females and the likelihood for more than one CO event in chromosomes < 10 Mb (n = 9) can be considered low, and was not observed). As predicted, the number of instances with one CO event (24) was similar to the number of instances without a detected CO event (21) in these small chromosomes. The distribution of transmission of recombinant versus non-recombinant chromosomes for the whole data set is shown in S2 Fig The higher number of transmissions of non-recombinant chromosomes in female than in male meiosis is a logical consequence of the lower rate of COs in females. There was a general trend of a higher frequency of CO events towards the ends of chromosomes (frequency measured in 10 Mb windows; Fig 3A, S3 Fig), primarily in chromosomes 50–100 Mb in size and in male meiosis. Moreover, within the terminal 10 Mb there was a markedly higher frequency of events towards the ends (Fig 3B; frequency measured in 1 Mb windows). Interestingly, also at this scale there was distinct difference between the sexes in that the frequency of CO events in male meiosis increased towards the very end while there was no CO event observed in the terminal 1 Mb in female meiosis (Mann-Whitney U-test, z = 1. 87, p = 0. 030). A similar trend was seen between the frequency of NCOs and distance to chromosome end (S4 Fig). We tested whether the location of one CO event affected the location of other events on the same chromosome or if the locations were independent of each other. There was evidence for positive interference–lowered likelihood of two nearby CO events–up to a distance of 14 Mb (Fig 4). Interestingly, out of 107 detected double CO events, only 28 were detected in female meiosis, while 79 were detected in male meiosis. Together with positive interference, this could explain why the observed increase in CO rate towards chromosomes ends was more pronounced in males than in females. The distribution of CO events along chromosomes in the flycatcher genome indicates that there are several distinct regions with a concentration of CO events, i. e. CO hot-spots (S5–S7 Figs). We identified 19 regions on 12 different chromosomes with two or more CO events from independent meiosis localized to less than 100 kb apart, demonstrating a highly skewed distribution of CO events. Randomly placing CO events along chromosomes indicated that the likelihood for this to happen by chance was <0. 001. Excluding CO events that overlapped with gaps between scaffolds, randomization revealed that the likelihood for the observed amount of co-localized CO events by chance was <0. 005. There was a concentration of CO events to genic regions (Fig 5A). To analyse the association between genes and COs in some further detail we divided the genome into promoter regions (2 kb upstream of transcription start site, TSS), first exons, first introns, other exons, other introns and intergenic DNA. The CO rate was highest in promoter regions (1. 85 times the intergenic rate), followed by first exons and other exons (Fig 5B). Among assembled parts of the genome, the number of CO events per bp in promoter regions was significantly higher than in intergenic DNA (Fisher’s Exact test, p = 0. 018). No statistically significant differences were detected between other comparisons, which may be due to the limited number of CO events in exons and introns. Moreover, there was a significant association between COs and CpG islands (p = 8. 72x10-9). Given that CpG islands are prevalent upstream of genes, the overrepresentation of CO events in promoters and CpG islands is likely not independent from each other. Besides, the GC content of CO regions (45. 6%) was higher than in the genomic background (41. 9%). This may in part be due to the fact that CpG islands are high in GC content and to the disproportionate number of CO events on the small microchromosomes (given one obligate CO, independent of chromosome size), and the fact that GC content increases with decreasing chromosome size in birds [22]. However, the GC content of CO regions on each chromosome was significantly higher than the genomic background on the respective chromosomes (paired t-test, p = 0. 0014). We did not find any evidence for a higher repeat density in CO regions (9. 1% vs. 8. 2%; p = 0. 65). The small sample size of NCO events gives limited power for statistical analyses of the relative abundance of NCO events in different functional categories (Fig 5C). However, a higher density in first exons, i. e. close to TSS, compared to intergenic regions was close to significant (odds ratio = 2. 495, p = 0. 081), and would resemble the situation for CO events. A significant overlap in the localization of CO and NCO events compared to random expectations was observed (odds ratio = 6. 321, p = 0. 0041), providing further support for a common mechanism of regulation. The genomic landscape of species divergence in flycatchers is characterized by the presence of numerous (≈50) ‘differentiation islands’ spread across the genome, evident as distinct FST peaks in comparisons between species [9,61]. These islands cover roughly 7% of the genome and may primarily result from lowered Ne due to linked selection in regions of low CO rate [9]. If the genomic locations of CO and NCO events and the 50 differentiation islands were unrelated, we would expect to see approximately 7% of these events to overlap with islands by chance. However, only 11 out of 325 CO events overlapped with islands (odds ratio = 0. 454, p = 0. 0067). In contrast, NCO events were over-represented in islands; 33 out of 266 NCO events overlapped (odds ratio = 1. 837, p = 0. 0026). Considering both types of recombination events taken together, their distribution relative to differentiation islands did not differ significantly from 7% of overlap. Both CO and NCO can lead to tracts of gene conversion close to the location of DSBs [12,64]. However, since we only trace one product of meiosis, we cannot track gene conversion tracts at CO events. Of the 267 NCO events, 229 involved sites segregating for one ‘weak’ (‘W’; A or T) and one ‘strong’ (‘S’; G or C) allele. We then counted the number of times a weak allele was converted by a strong allele (W>S) and the number of times a strong allele was converted by a weak allele (S>W). If transmission of alleles upon gene conversion is a random process there should be about as many events of one category as of the other. However, there was a significant excess of W>S conversions (Binomial test, p = 0. 012), with a biased transmission of 59% (95% CI = 0. 52–0. 65). This corresponds to a transmission distortion (c) of 0. 18 and provides direct evidence for GC-biased gene conversion in an avian species.
Studies in several species of plants and fungi have shown that recombination events are concentrated in hot-spots in close proximity to TSSs [28,39–43]. This may be a consequence of the common pattern that both transcription and recombination are facilitated in open chromatin [68,69]. In humans and mouse, PRDM9 directs recombination away from TSS [51,52,67,70,71]. The absence of an active Prdm9 gene in avian genomes [60] prompts the hypothesis that recombination in this vertebrate lineage resembles the ancestral mechanism of regulation and is associated with proximity to TSSs, similar to plants and fungi. Our observations confirmed this hypothesis: the highest rate of CO events was seen in promoter regions and first exons, with statistical support for the rate in the former being higher than in intergenic regions. We also found an association between COs and CpG islands, which commonly function as promoters by destabilizing nucleosomes and attracting proteins that create a transcriptionally permissive chromatin state [72]. Also NCO events tended to be most common close to TSSs and, overall, there was evidence for an overlap in the distribution of CO and NCO events. As pointed out by Lichten, this suggests “that the picture in mammals may be the exception rather than the rule” [73]. We note that there are exceptions to an increased rate of recombination close to TSS in the absence of Prdm9. In both Drosophila melanogaster [74] and D. pseudoobscura [75], recombination is reduced around TSS. Furthermore, fruit flies [33,75], worms [76] and honeybee [77] have recombination landscapes that are relatively homogenous without distinct hotspots. Apparently, despite being such a widespread phenomenon across the tree of life, recombination has evolved distinct characteristics in different lineages. The usage of human and mouse recombination hot-spots has a high turnover rate [78,79]. This owes in part to large allelic variation of, and positive selection on, the DNA-binding residues of zinc fingers [53], such that it represents an allelic turnover of the binding protein. However, erosion of binding motifs in target sequence due to recombination-induced mutation or gene conversion probably also occurs [2,47]. The concentration of recombination events to hot-spot regions in the absence of an active Prdm9 gene, like what we found in this study of flycatchers, could potentially mean that the hot-spot landscape of recombination in such lineages remains relative stable (cf. [80]). Evidence for this was recently found in an avian study of finch species [54] and in analyses of divergent Saccharomyces species [81], where the location of hot-spots was found to be conserved over considerable evolutionary time. We found a higher rate of COs in male than in female meiosis, i. e. higher rate in the homogametic than in the heterogametic sex. Since we only measured the rate of COs in one individual of each sex, we cannot formally exclude the possibility that part of the observed rate difference was due to different genetic backgrounds rather than to sex per se. However, the higher male CO rate accords with collared flycatcher linkage map data [55,82]. Birds do not uniformly follow the Haldane-Huxley rule (reduced CO rate in the heterogametic sex, [83]) since there in addition to species with male-biased CO rate [84] are others with female-biased CO rate [85] or similar rates of COs in the two sexes [24,56,65,86]. Haldane [87] and Huxley [88] suggested that reduced CO rate in the heterogametic sex was a pleiotropic consequence of selection against CO between diverging X and Y (or Z and W) sex chromosomes. Several alternative hypotheses have subsequently been put forward (see [85]), of which some could potentially explain variation in the relative rates of male and female COs within organism groups. For example, sexual selection may select for reduced CO rate (to maintain favourable allelic combinations) in the sex with the largest variance in reproductive success [89], potentially setting the stage for a relationship between the intensity of sexual selection (or sexual antagonism) and male-to-female CO rate ratio [90,91]. As for other animal groups (e. g. [92]), linkage analysis in several bird species has revealed a general trend of increased rates of COs towards chromosome ends [24,55,56], the extent of which apparently varies among species with the most pronounced end effect so far seen in zebra finch [57,65]. Our observations add to and complement this picture by demonstrating a sex difference in the distribution of COs along chromosomes. Specifically, the increased CO rate towards chromosome ends was mainly seen in males, similar to the situation in humans and mice [93], while female CO events were more evenly distributed. The findings of positive CO interference up to about 14 Mb and a higher incidence of double COs in males than in females are compatible with increased rates of COs towards chromosome ends in males. Only a limited number of recombination-initiating DSBs eventually result in COs, whose spatial distribution is tightly regulated through the process of CO interference that reduces the possibility of two nearby CO events [94]. As a result, double COs tend to be directed towards the respective ends of chromosomes. Measures of the extent of interference vary among taxa with a tendency for shorter distances seen in organisms with higher CO rates than in organisms with lower CO rates [95]. Estimates of CO interference for human and mice, for example, range between 20 and 140 Mb [35,96]. The relevant metric in terms of CO interference is though the physical distance in μm along the bivalent, not the “genomic” distance in Mb along DNA. The distance over which interference occur thus depends on the degree of compaction of chromosomes at the leptotene stage of meiosis. This may explain sex-differences in the spatial distribution of COs, since males and females show varying degrees of compaction of the chromosomes [96]. Detailed mapping of recombination events in human sperm and oocytes has documented that CO interference is more pronounced in males than in females [34,35]. It may very well be that other meiotic characteristics of spermatogenesis and oogenesis contribute to sex-differences in patterns of recombination, like differences in the time allotted to the bouquet formation at telomeres [93]. An interesting observation with respect to the distribution of COs along chromosomes was the significant under-representation of CO events in genomic regions defined as differentiation islands. This accords with the findings that differentiation islands are concentrated, if not limited, to regions of the genome corresponding to CO desserts, as judged by CO rate data from linkage maps, and that there is an overall positive correlation between CO rate and FST [9]. As such, this corroborates the notion that CO rate drives the genomic landscape of species differentiation in Ficedula flycatchers [9]. The rationale for this inference is that the prevalence of the diversity-reducing (Ne-reducing) effects of linked selection increases with decreasing CO rate [97]. Because the role of genetic drift on differentiation in turn increases with decreasing Ne, this means that variation in the degree of differentiation across the genome is compatible with a neutral model, without the need to invoke selection or varying degree of gene flow. Moreover, since both types of recombination events taken together did not differ significantly from a random overlap and NCO events were over-represented in islands, this indicates that DSBs that occur in these differentiation islands are preferentially assigned as NCOs. We found a significant transmission distortion at NCO events with the ‘strong’ (G, C) allele transmitted at 59% of all events involving one ‘strong’ and one ‘weak’ (A, T) allele (transmission distortion, c = 0. 18). Detection of a transmission distortion for strong alleles has so far been limited to studies of humans and yeast [12,14,64,98,99], and is likely explained by GC-biased gene conversion (gBGC). gBGC is a process associated with meiotic recombination, which favours strong base-pairs over weak base-pairs at weak: strong mismatches in heteroduplex DNA formed as part of the repair mechanism of DSBs. This ultimately leads to a preferential transmission of strong alleles over weak alleles close to recombination-initiating DSBs. As a consequence, the local rate of recombination is expected to show a positive correlation with the local GC content. Indeed, such indirect evidence for gBGC has been observed across a wide range of taxa [11], while direct detection is much more rare. In yeast a transmission distortion of 0. 057 was observed at CO events, but no significant distortion was observed at NCO events [12]. In humans there is evidence for a transmission distortion associated with both CO and NCO events [64,98]. A recent genome-wide study of NCO events in humans estimated a transmission distortion of 0. 36 [99]. The estimate that we report here for flycatchers falls below the estimate for humans. However, the net impact of gBGC on the evolution of base composition does not only depend on the strength of the transmission distortion, but also on the number of weak: strong mismatches in heteroduplex DNA, the recombination rate and Ne [11]. Given higher SNP density [9], higher recombination rate [55] and larger Ne in flycatchers [100] compared to humans, this might readily account for a higher genome-wide GC content in flycatchers compared to humans. We have previously suggested that the slow rate of karyotypic evolution in birds will promote a conserved genomic landscape of recombination rate variation and thereby facilitate the evolutionary build-up of genomic signatures of recombination, like the effect of gBGC on base composition [101]. The presence of recombination hot-spots coupled with a stable hot-spot landscape in the absence of Prdm9 further accentuates the influence of recombination rate variation on avian molecular evolution. For example, in a recent flycatcher study we found an at first glance unexpected absence of a correlation between recombination rate and the rate of non-synonymous substitutions [102]. However, the patterns changed when GC-biased gene conversion was taken into account and weak-to-strong and strong-to-weak substitutions were separately analysed.
This study was approved by Linköpings djuretiska nämnd, Linköpings tingsrätt, Sweden (Dnr 21–11). We sequenced a three-generation pedigree (Fig 1A) of 11 collared flycatchers, sampled in the field from a natural population on the Baltic Sea island Öland (Sweden). The four birds in the P generation showed no evidence of being closely related. Sequencing was done to approximately 40X coverage on an Illumina HiSeq 2000 instrument with paired-end reads of 100 bp and libraries (insert size of 450 bp) constructed using the TruSeq Nano sample preparation kit (Illumina) (European Nucleotide Archive PRJEB12616). DNA were prepared from blood samples stored in 96% ethanol using a standard proteinase K digestion/phenol-chloroform purification protocol. The reads were aligned to the collared flycatcher reference genome FicAlb1. 5 (GenBank Accession GCA_000247815. 2) with bwa 0. 7. 5a [103] and de-duplicated, recalibrated and cleaned with GATK 3. 2. 2 [104,105]. Single nucleotide polymorphisms (SNPs) were called with GATK' s HaplotypeCaller and GenotypeGVCFs (version 3. 3. 0). Variant Quality Score Recalibration (VQSR) was performed according to the GATK' s Best Practice [106], using known SNP positions from genotyping with a 50k SNP-chip [107] and the 20% top scoring sites for training. Sites that failed the 99. 9% tranche threshold were removed. We also removed sites in repetitive regions using a combination of RepeatMasker v3. 2. 9 (Smit, AFA, Hubley, R & Green, P. <http: //www. repeatmasker. org>) and a flycatcher-specific repeat library [108], Tandem Repeats Finder v4. 07 [109] and an in-house perl script masking homopolymers longer than 10bp that were not already masked. To further reduce the number of mis-genotyped sites, we applied a coverage filter removing sites where any of the involved individuals were covered by less than 15 reads, or covered by more than twice the average autosomal coverage; the latter criterion was applied to reduce the risk of using collapsed regions where sites potentially could be mis-called as heterozygous due to differences between sequence copies. We also filtered sites with low genotype quality (GQ<30), which corresponds to that the probability of choosing the wrong genotype is less than 0. 001. Sites violating Mendelian inheritance or with more than two called alleles were removed. All these measures ensured stringent filtering. The entire genome including the Z chromosome was variant-called as diploid. To avoid incorrectly assigned genotypes for females (which in birds represent the heterogametic sex, with one copy of the Z chromosome plus the female-specific W chromosome), the Z chromosome was called separately as haploid (flag -ploidy 1 in GATK) in the four females (two in the P generation and one each in the F1 and F2 generations); the pseudo-autosomal region of the Z chromosome [62] was excluded from this analysis. VQSR could not be used for this set due to lack of a proper training set, but we filtered using coverage and genotype quality. 15X was used as the lower coverage threshold as above, and maximum coverage was set to the mean autosomal coverage (i. e. , twice the expected coverage for the Z chromosome). The new haploid calls for the females where then added to the data set. With a three-generation family it is possible to phase the F1 generation into chromosome-level haplotypes and thereby detect recombination events in the transmission of gametes to the F2 generation (Fig 1C). For a variant site to be informative in phasing, it is required that an F1 individual is heterozygous and that its two parents have different genotypes; this situation makes it possible to trace each F1 allele to one of the parents in the P generation. It also requires that F1 alleles are traceable in the F2 generation and the F1 individual' s partner and offspring must therefore not both be heterozygous (S8 Fig). Phased sites were grouped into haploblocks in the five F2 individuals. Similar to [32] we assumed that there would be at most one recombination event per Mb interval and considered shorter blocks as the likely result of phasing or genotyping error (see further below). CO events were localized to the genomic interval between the outermost SNPs of two adjacent haploblocks in the F2 individuals. The resolution of intervals varied depending on the density of informative SNPs in the regions in question. We limited downstream analyses to events that were localized with a resolution <5 kb. The location of 25 recombination intervals overlapped with gaps between scaffolds meaning that the precise size could not be determined. NCO events are suggested at informative positions with a phase that does not match the surrounding block. However, this will also be the case whenever there is an incorrect genotype call in any of the involved individuals. Manual inspection in IGV [110] of a subset of phase-mismatched sites showed that most of them had visible problems such as high coverage (close to our upper threshold), unequal numbers of reads supporting the two alleles at a site, reads supporting more than two visible alleles, clusters of nearby polymorphic sites, and overlaps with insertions and deletions. There was also an excess of sites at which all females were heterozygous and all males were homozygous for the reference allele, indicative of reads from W-linked sequences mapping to autosomal or Z-linked loci in the male-derived genome assembly. To reduce the above-mentioned problems in inferring potential NCO events we applied several additional stringent filtering criteria to the set of informative sites selected for the identification of CO events. First we used a more strict VQSR tranche threshold of 90%. We further excluded all sites overlapping with, or present within 10 bp of, indels called by GATK and clusters of SNPs that had more than three called SNPs within 30bp in the full VQSR-filtered file (containing all individuals). We also excluded larger clusters of SNPs with frequent and alternating haplotype shifts; based on visual inspecting we set the threshold to no more than two deviating sites in 5 kb. Next, we removed sites that had reads supporting more than two alleles, sites where one of the alleles was supported by < 25% of the reads, and sites where all females were heterozygous and all males were homozygous for the reference allele. We considered this filtering necessary to remove ambiguous sites although it may have come to the price of underestimating the occurrence of NCO events. Because of this we do not investigate the relative frequencies of CO and NCO events. Collared flycatcher genes were downloaded from Ensembl (release 73, assembly version FicAlb_1. 4) and translated to the latest assembly version with chromosomes using an in-house script. For calculating the distance between CO events and genes, we restricted the analysis to CO events with a resolution <5 kb. Depending on the orientation of the closest gene we assigned the CO event to the upstream or the downstream flank. If the CO event overlapped with a gene, we assigned it a distance of 0. In a second step we assigned CO events into six different classes of genomic regions; intergenic, promoter (defined as 2 kb upstream of the transcription start site, TSS), first exon, first intron, other exons and other introns. A CO event that overlapped several classes was assigned values to each of these classes proportional to the length of the overlap. We repeated the analysis for NCO events. Next, CpG islands (CGIs) were identified for the hard-masked flycatcher genome using CpGcluster (version 1. 0) with default parameter settings [111]. In order to assess the association between CGIs and CO events, the number of overlapping CGIs was compared to the genome-wide average. All statistics were calculated and plotted with R version 3. 0. 2 (http: //www. R-project. org/). We used the coefficient of coincidence (CoC) to assess the strength of CO interference [28]. CoC was computed as the number of observed over expected double COs counted in a 1Mb sliding windows approach. This provided us with a sex-average and genome-average CoC.
Homologous chromosomes exchange genetic material during cell division at meiosis by the process of crossover recombination. Although such crossover events are visually identifiable by cytogenetic techniques, it has remained a challenge to pinpoint the location of crossovers at the DNA sequence level. An emerging novel possibility to approach this challenge is to exploit the high resolution offered by re-sequencing of multiple individuals in species in which a genome assembly is available. Specifically, by sequencing members of a family, the inheritance of chromosomal segments can be followed and the location of crossover as well as non-crossover recombination determined. We performed such an endeavour in the collared flycatcher, a songbird species that has been in focus for extensive ecological and evolutionary research. We found that crossover events were concentrated to certain ‘hot-spot’ regions, and that the density of such events was highest in and close to genes. A higher rate of crossover recombination was found in males than in females and in males, but not in females, the rate of crossover increased towards ends of chromosomes. The location of non-crossovers was significantly associated with that of crossovers. We could further document an unbalanced transmission of genetic variants at non-crossover events via to the process of GC-biased gene conversion.
meiosis
cell cycle and cell division
cell processes
vertebrates
animals
alleles
dna recombination
dna
gene conversion
sex chromosomes
birds
chromosome biology
animal genomics
genetic loci
bird genomics
biochemistry
cell biology
nucleic acids
genetics
biology and life sciences
genomics
amniotes
organisms
chromosomes
High-Resolution Mapping of Crossover and Non-crossover Recombination Events by Whole-Genome Re-sequencing of an Avian Pedigree