Drosophila recombination rate calculator - Methods

Recombination rates are often estimated based on the relationship between the physical and genetic maps of the Drosophila genome. Early recombination estimates, such as the adjusted coefficient of exchange (Kindahl 1994), were based on cytological markers in combination with DNA content estimates from optical densities of polytene chromosomes (Sorsa 1988). Other techniques such as R_TE (Hey and Kliman 2002) use transposable elements as markers on the physical map (Ising and Block 1984). However, with the availability of the whole genome sequence of, the physical map of the genome is much more precise, which facilitates more accurate estimation of recombination rates in D. melanogaster(Hey and Kliman 2002; Marais et al. 2001).

Recombination rates are now estimated by plotting Marey maps of the genetic positions of molecular markers (in centimorgans, cM) against their physical position (in Megabase pairs, Mbp). The recombination rate at any given nucleotide coordinate is estimated by taking the slope of the curve relating the genetic to the physical maps in one of two ways. One could employ a sliding window approach, in which a linear function is used to fit the genetic position as a function of physical position within a window of a given size along a chromosome arm. With this approach, the slope of this linear function within the window is taken as local recombination rate. Alternatively, the genetic position of the markers can be mapped as a function of physical position using an nth degree polynomial curve across an entire chromosome arm. Under this approach, recombination is estimated at an individual nucleotide coordinate as the derivative of the nth degree polynomial curve. While the polynomial curve approach is less sensitive to regional variation in recombination rates than the sliding window approach, it is more robust to errors in the physical and genetic maps.

There are currently 656 genes that have been localized on both the physical and genetic maps of Release 4.2 of the D. melanogaster genome, which were kindly provided to us by David Sutherland at Flybase. Using these genes as markers, we plotted Marey maps (see below) for each of the five chromosome arms with detectable recombination (2R, 2L, 3L, 3R, and X). Recombination is thought to be absent or negligible on the fourth chromosome. We visually identified outliers and removed them (n = 3, 3, 1, 2, and 1 outliers on chromosome arms 2L, 2R, 3L, 3R and X, respectively) and we fit a third-order polynomial curve to the genetic map position as a function of physical position for the remaining markers on each chromosome arm (n = 114, 112, 100, 151, and 169 genes for chromosome arms 2L, 2R, 3L, 3R and X, respectively). In all cases, the fit was quite good; R² > 0.96 for all chromosome arms.

We estimate recombination as the derivative of this 3^rd order polynomial curve at a given nucleotide coordinate. For any locus, this web tool estimates recombination at the start coordinate, end coordinate, as well as the midpoint of the locus. Because of the smoothing effect of this recombination estimation technique, there will be little difference in the recombination rate estimates at these three points for small loci. However, we provide all three estimations in the event that the locus of interest is sufficiently large for regional variation in recombination rate.

Marey Maps:

Marey maps of chromosome arms of the D. melanogaster genome, based on the 656 genes in Release 4 of the genome sequence that have been localized on both genetic and physical maps. A 3rd order polynomial curve mapping the genetic location onto the physical location is also represented.

References:

Hey J, Kliman RM (2002) Interactions between natural selection, recombination and gene density in the genes of Drosophila. Genetics 160:595-608.

Ising G, Block K (1984) A transposon as a cytogenetic marker in Drosophila melanogaster. Mgg (Molecular & General Genetics) 196:6-16.

Kindahl EC (1994) Recombination and DNA polymorphism on the third chromosome of Drosophila melanogaster. Cornell University, Ithaca, NY

Marais G, Mouchiroud D, Duret L (2001) Does recombination improve selection on codon usage? Lessons from nematode and fly complete genomes. Proceedings of the National Academy of Sciences of the United States of America 98:5688-5692.

Sorsa V (1988) Chromosome Maps of Drosophila. CRC Press, Boca Raton, FL