Studying Recombination with High-Throughput Sequencing: An Educational Primer for Use with “Fine-Scale Heterogeneity in Crossover Rate in the garnet-scalloped Region of the Drosophila melanogaster X chromosome”

An article by Singh and colleagues in this issue of GENETICS quantifies variation in recombination rate across a small region of the Drosophila melanogaster genome, providing an opportunity for instructors of genetics to introduce or reinforce important concepts such as recombination and recombination rate variation, genome sequencing, and sequence features of the genome. Additional background information, a detailed explanation of the methods used in this study, and discussion questions are provided.Related article in GENETICS: Singh, N. D., E. A. Stone, C. F. Aquadro, and A. G. Clark, 2013 Fine-scale heterogeneity in crossover rate in the garnet-scalloped region of the Drosophila melanogaster X chromosome. Genetics 194: 375–387.     

Crossover Heterogeneity in the Absence of Hotspots in Caenorhabditis elegans

Crossovers play mechanical roles in meiotic chromosome segregation, generate genetic diversity by producing new allelic combinations, and facilitate evolution by decoupling linked alleles. In almost every species studied to date, crossover distributions are dramatically nonuniform, differing among sexes and across genomes, with spatial variation in crossover rates on scales from whole chromosomes to subkilobase hotspots. To understand the regulatory forces dictating these heterogeneous distributions a crucial first step is the fine-scale characterization of crossover distributions. Here we define the wild-type distribution of crossovers along a region of the C. elegans chromosome II at unprecedented resolution, using recombinant chromosomes of 243 hermaphrodites and 226 males. We find that well-characterized large-scale domains, with little fine-scale rate heterogeneity, dominate this region’s crossover landscape. Using the Gini coefficient as a summary statistic, we find that this region of the C. elegans genome has the least heterogeneous fine-scale crossover distribution yet observed among model organisms, and we show by simulation that the data are incompatible with a mammalian-type hotspot-rich landscape. The large-scale structural domains—the low-recombination center and the high-recombination arm—have a discrete boundary that we localize to a small region. This boundary coincides with the arm-center boundary defined both by nuclear-envelope attachment of DNA in somatic cells and GC content, consistent with proposals that these features of chromosome organization may be mechanical causes and evolutionary consequences of crossover recombination.     

Estimation of Fine-Scale Recombination Intensity Variation in the <Emphasis Type="Italic">white–echinus</Emphasis> Interval of <Emphasis Type="Italic">D. melanogaster</Emphasis>

Accurate assessment of local recombination rate variation is crucial for understanding the recombination process and for determining the impact of natural selection on linked sites. In Drosophila, local recombination intensity has been estimated primarily by statistical approaches, by estimating the local slope of the relationship between the physical and genetic maps. However, these estimates are limited in resolution and, as a result, the physical scale at which recombination intensity varies in Drosophila is largely unknown. Although there is some evidence suggesting as much as a 40-fold variation in crossover rate at a local scale in D. pseudoobscura, little is known about the fine-scale structure of recombination rate variation in D. melanogaster. Here we experimentally examine the fine-scale distribution of crossover events in a 1.2-Mb region on the D. melanogaster X chromosome using a classic genetic mapping approach. Our results show that crossover frequency is significantly heterogeneous within this region, varying approximately 3.5-fold. Simulations suggest that this degree of heterogeneity is sufficient to affect levels of standing nucleotide diversity, although the magnitude of this effect is small. We recover no statistical association between empirical estimates of nucleotide diversity and recombination intensity, which is likely due to the limited number of loci sampled in our population genetic data set. However, codon bias is significantly negatively correlated with fine-scale recombination intensity estimates, as expected. Our results shed light on the relevant physical scale to consider in evolutionary analyses relating to recombination rate and highlight the motivations to increase the resolution of the recombination map in Drosophila.     

Species coherence in the face of karyotype diversification in holocentric organisms: the case of a cytogenetically variable sedge (Carex scoparia,Cyperaceae)

Background and Aims

The sedge genus Carex, the most diversified angiosperm genus of the northern temperate zone, is renowned for its holocentric chromosomes and karyotype variability. The genus exhibits high variation in chromosome numbers both among and within species. Despite the possibility that this chromosome evolution may play a role in the high species diversity of Carex, population-level patterns of molecular and cytogenetic differentiation in the genus have not been extensively studied.

Methods

Microsatellite variation (11 loci, 461 individuals) and chromosomal diversity (82 individuals) were investigated in 22 Midwestern populations of the North American sedge Carex scoparia and two Northeastern populations.

Key Results

Among Midwestern populations, geographic distance is the most important predictor of genetic differentiation. Within populations, inbreeding is high and chromosome variation explains a significant component of genetic differentiation. Infrequent dispersal among populations separated by >100 km explains an important component of molecular genetic and cytogenetic diversity within populations. However, karyotype variation and correlation between genetic and chromosomal variation persist within populations even when putative migrants based on genetic data are excluded.

Conclusions

These findings demonstrate dispersal and genetic connectivity among widespread populations that differ in chromosome numbers, explaining the phenomenon of genetic coherence in this karyotypically diverse sedge species. More generally, the study suggests that traditional sedge taxonomic boundaries demarcate good species even when those species encompass a high range of chromosomal diversity. This finding is important evidence as we work to document the limits and drivers of biodiversity in one of the world''s largest angiosperm genera.     

Surprising Differences in the Variability of Y Chromosomes in African and Cosmopolitan Populations of Drosophila melanogaster

The nonrecombining Drosophila melanogaster Y chromosome is heterochromatic and has few genes. Despite these limitations, there remains ample opportunity for natural selection to act on the genes that are vital for male fertility and on Y factors that modulate gene expression elsewhere in the genome. Y chromosomes of many organisms have low levels of nucleotide variability, but a formal survey of D. melanogaster Y chromosome variation had yet to be performed. Here we surveyed Y-linked variation in six populations of D. melanogaster spread across the globe. We find surprisingly low levels of variability in African relative to Cosmopolitan (i.e., non-African) populations. While the low levels of Cosmopolitan Y chromosome polymorphism can be explained by the demographic histories of these populations, the staggeringly low polymorphism of African Y chromosomes cannot be explained by demographic history. An explanation that is entirely consistent with the data is that the Y chromosomes of Zimbabwe and Uganda populations have experienced recent selective sweeps. Interestingly, the Zimbabwe and Uganda Y chromosomes differ: in Zimbabwe, a European Y chromosome appears to have swept through the population.     

The genetics and cytology of a mutator factor in Drosophila melanogaster

A Drosophila melanogaster mutator factor is described whose effects include the induction of unique chromosomal aberrations and male crossing over. Results of experiments to map the factor suggest that genetic transmission is somehow chromosomally associated but not localizable to the X, Y, second or third chromosome. There appears to be a good correlation between the distributions of male crossover exchange points and unique aberration breakpoints for the second chromosome but not for the third chromosome. The male crossovers, which occur more frequently in the centromeric region, occur in euchromatin rather than in the centric heterochromatin. The male crossovers tend to be rather precise reciprocal exchanges, since cytologically detectable deletions and duplications are only infrequently produced. It is suggested that the present mutator may be identical to earlier reported mutators of D. melanogaster.     

Nucleotide diversity in Silene latifolia autosomal and sex-linked genes

The plant Silene latifolia has separate sexes and sex chromosomes, and is of interest for studying the early stages of sex chromosome evolution, especially the evolution of non-recombining regions on the Y chromosome. Hitch-hiking processes associated with ongoing genetic degeneration of the non-recombining Y chromosome are predicted to reduce Y-linked genes'' effective population sizes, and S. latifolia Y-linked genes indeed have lower diversity than X-linked ones. We tested whether this represents a true diversity reduction on the Y, versus the alternative possibility, elevated diversity at X-linked genes, by collecting new data on nucleotide diversity for autosomal genes, which had previously been little studied. We find clear evidence that Y-linked genes have reduced diversity. However, another alternative explanation for a low Y effective size is a high variance in male reproductive success. Autosomal genes should then also have lower diversity than expected, relative to the X, but this is not found in our loci. Taking into account the higher mutation rate of Y-linked genes, their low sequence diversity indicates a strong effect of within-population hitch-hiking on the Y chromosome.     

Preferential Breakpoints in the Recovery of Broken Dicentric Chromosomes in Drosophila melanogaster

We designed a system to determine whether dicentric chromosomes in Drosophila melanogaster break at random or at preferred sites. Sister chromatid exchange in a Ring-X chromosome produced dicentric chromosomes with two bridging arms connecting segregating centromeres as cells divide. This double bridge can break in mitosis. A genetic screen recovered chromosomes that were linearized by breakage in the male germline. Because the screen required viability of males with this X chromosome, the breakpoints in each arm of the double bridge must be closely matched to produce a nearly euploid chromosome. We expected that most linear chromosomes would be broken in heterochromatin because there are no vital genes in heterochromatin, and breakpoint distribution would be relatively unconstrained. Surprisingly, approximately half the breakpoints are found in euchromatin, and the breakpoints are clustered in just a few regions of the chromosome that closely match regions identified as intercalary heterochromatin. The results support the Laird hypothesis that intercalary heterochromatin can explain fragile sites in mitotic chromosomes, including fragile X. Opened rings also were recovered after male larvae were exposed to X-rays. This method was much less efficient and produced chromosomes with a strikingly different array of breakpoints, with almost all located in heterochromatin. A series of circularly permuted linear X chromosomes was generated that may be useful for investigating aspects of chromosome behavior, such as crossover distribution and interference in meiosis, or questions of nuclear organization and function.     

DO MALES MATTER? TESTING THE EFFECTS OF MALE GENETIC BACKGROUND ON FEMALE MEIOTIC CROSSOVER RATES IN DROSOPHILA MELANOGASTER

Meiotic recombination is a critical genetic process as well as a pivotal evolutionary force. Rates of crossing over are highly variable within and between species, due to both genetic and environmental factors. Early studies in Drosophila implicated female genetic background as a major determinant of crossover rate and recent work has highlighted male genetic background as a possible mediator as well. Our study employed classical genetics to address how female and male genetic backgrounds individually and jointly affect crossover rates. We measured rates of crossing over in a 33 cM region of the Drosophila melanogaster X chromosome using a two‐step crossing scheme exploiting visible markers. In total, we measured crossover rates of 10 inbred lines in a full diallel cross. Our experimental design facilitates measuring the contributions of female genetic background, male genetic background, and female by male genetic background interaction effects on rates of crossing over in females. Our results indicate that although female genetic background significantly affects female meiotic crossover rates in Drosophila, male genetic background and the interaction of female and male genetic backgrounds have no significant effect. These findings thus suggest that male‐mediated effects are unlikely to contribute greatly to variation in recombination rates in natural populations of Drosophila.     

Y-linked variation for autosomal immune gene regulation has the potential to shape sexually dimorphic immunity

Sexually dimorphic phenotypes arise from the differential expression of male and female shared genes throughout the genome. Unfortunately, the underlying molecular mechanisms by which dimorphic regulation manifests and evolves are unclear. Recent work suggests that Y-chromosomes may play an important role, given that Drosophila melanogaster Ys were shown to influence the regulation of hundreds of X and autosomal genes. For Y-linked regulatory variation (YRV) to facilitate sexually dimorphic evolution, however, it must exist within populations (where selection operates) and influence male fitness. These criteria have seldom been investigated, leaving the potential for dimorphic evolution via YRV unclear. Interestingly, male and female D. melanogaster differ in immune gene regulation. Furthermore, immune gene regulation appears to be influenced by the Y-chromosome, suggesting it may contribute to dimorphic immune evolution. We address this possibility by introgressing Y-chromosomes from a single wild population into an isogenic background (to create Y-lines) and assessing immune gene regulation and bacterial defence. We found that Y-line males differed in their immune gene regulation and their ability to defend against Serratia marcescens. Moreover, gene expression and bacterial defence were positively genetically correlated. These data indicate that the Y-chromosome has the potential to shape the evolution of sexually dimorphic immunity in this system.     

Testing for the Footprint of Sexually Antagonistic Polymorphisms in the Pseudoautosomal Region of a Plant Sex Chromosome Pair

The existence of sexually antagonistic (SA) polymorphism is widely considered the most likely explanation for the evolution of suppressed recombination of sex chromosome pairs. This explanation is largely untested empirically, and no such polymorphisms have been identified, other than in fish, where no evidence directly implicates these genes in events causing loss of recombination. We tested for the presence of loci with SA polymorphism in the plant Silene latifolia, which is dioecious (with separate male and female individuals) and has a pair of highly heteromorphic sex chromosomes, with XY males. Suppressed recombination between much of the Y and X sex chromosomes evolved in several steps, and the results in Bergero et al. (2013) show that it is still ongoing in the recombining or pseudoautosomal, regions (PARs) of these chromosomes. We used molecular evolutionary approaches to test for the footprints of SA polymorphisms, based on sequence diversity levels in S. latifolia PAR genes identified by genetic mapping. Nucleotide diversity is high for at least four of six PAR genes identified, and our data suggest the existence of polymorphisms maintained by balancing selection in this genome region, since molecular evolutionary (HKA) tests exclude an elevated mutation rate, and other tests also suggest balancing selection. The presence of sexually antagonistic alleles at a locus or loci in the PAR is suggested by the very different X and Y chromosome allele frequencies for at least one PAR gene.     

Within‐population Y‐linked genetic variation for lifespan in Drosophila melanogaster

The view that the Y chromosome is of little importance for phenotypic evolution stems from early studies of Drosophila melanogaster. This species’ Y chromosome contains only 13 protein‐coding genes, is almost entirely heterochromatic and is not necessary for male viability. Population genetic theory further suggests that non‐neutral variation can only be maintained at the Y chromosome under special circumstances. Yet, recent studies suggest that the D. melanogaster Y chromosome trans‐regulates hundreds to thousands of X and autosomal genes. This finding suggests that the Y chromosome may play a far more active role in adaptive evolution than has previously been assumed. To evaluate the potential for the Y chromosome to contribute to phenotypic evolution from standing genetic variation, we test for Y‐linked variation in lifespan within a population of D. melanogaster. Assessing variation for lifespan provides a powerful test because lifespan (i) shows sexual dimorphism, which the Y is primarily predicted to contribute to, (ii) is influenced by many genes, which provides the Y with many potential regulatory targets and (iii) is sensitive to heterochromatin remodelling, a mechanism through which the Y chromosome is believed to regulate gene expression. Our results show a small but significant effect of the Y chromosome and thus suggest that the Y chromosome has the potential to respond to selection from standing genetic variation. Despite its small effect size, Y‐linked variation may still be important, in particular when evolution of sexual dimorphism is genetically constrained elsewhere in the genome.     

An Excess of Gene Expression Divergence on the X Chromosome in Drosophila Embryos: Implications for the Faster-X Hypothesis

The X chromosome is present as a single copy in the heterogametic sex, and this hemizygosity is expected to drive unusual patterns of evolution on the X relative to the autosomes. For example, the hemizgosity of the X may lead to a lower chromosomal effective population size compared to the autosomes, suggesting that the X might be more strongly affected by genetic drift. However, the X may also experience stronger positive selection than the autosomes, because recessive beneficial mutations will be more visible to selection on the X where they will spend less time being masked by the dominant, less beneficial allele—a proposal known as the faster-X hypothesis. Thus, empirical studies demonstrating increased genetic divergence on the X chromosome could be indicative of either adaptive or non-adaptive evolution. We measured gene expression in Drosophila species and in D. melanogaster inbred strains for both embryos and adults. In the embryos we found that expression divergence is on average more than 20% higher for genes on the X chromosome relative to the autosomes; but in contrast, in the inbred strains, gene expression variation is significantly lower on the X chromosome. Furthermore, expression divergence of genes on Muller''s D element is significantly greater along the branch leading to the obscura sub-group, in which this element segregates as a neo-X chromosome. In the adults, divergence is greatest on the X chromosome for males, but not for females, yet in both sexes inbred strains harbour the lowest level of gene expression variation on the X chromosome. We consider different explanations for our results and conclude that they are most consistent within the framework of the faster-X hypothesis.     

Quantitative genetic analysis suggests causal association between cuticular hydrocarbon composition and desiccation survival in Drosophila melanogaster

Survival to low relative humidity is a complex adaptation, and many repeated instances of evolution to desiccation have been observed among Drosophila populations and species. One general mechanism for desiccation resistance is Cuticular Hydrocarbon (CHC) melting point. We performed the first Quantitative Trait Locus (QTL) map of population level genetic variation in desiccation resistance in D. melanogaster. Using a panel of Recombinant Inbred Lines (RILs) derived from a single natural population, we mapped QTL in both sexes throughout the genome. We found that in both sexes, CHCs correlated strongly with desiccation resistance. At most desiccation resistance loci there was a significant association between CHCs and desiccation resistance of the sort predicted from clinal patterns of CHC variation and biochemical properties of lipids. This association was much stronger in females than males, perhaps because of greater overall abundance of CHCs in females, or due to correlations between CHCs used for waterproofing and sexual signalling in males. CHC evolution may be a common mechanism for desiccation resistance in D. melanogaster. It will be interesting to compare patterns of CHC variation and desiccation resistance in species which adapt to desiccation, and rainforest restricted species which cannot.     

Genomics of adaptive divergence with chromosome‐scale heterogeneity in crossover rate

Genetic differentiation between divergent populations is often greater in chromosome centres than peripheries. Commonly overlooked, this broadscale differentiation pattern is sometimes ascribed to heterogeneity in crossover rate and hence linked selection within chromosomes, but the underlying mechanisms remain incompletely understood. A literature survey across 46 organisms reveals that most eukaryotes indeed exhibit a reduced crossover rate in chromosome centres relative to the peripheries. Using simulations of populations diverging into ecologically different habitats through sorting of standing genetic variation, we demonstrate that such chromosome‐scale heterogeneity in crossover rate, combined with polygenic divergent selection, causes stronger hitchhiking and especially barriers to gene flow across chromosome centres. Without requiring selection on new mutations, this rapidly leads to elevated population differentiation in the low‐crossover centres relative to the high‐crossover peripheries of chromosomes (“Chromosome Centre‐Biased Differentiation”, CCBD). Using simulated and empirical data, we then show that strong CCBD between populations can provide evidence of polygenic adaptive divergence with a phase of gene flow. We further demonstrate that chromosome‐scale heterogeneity in crossover rate impacts analyses beyond that of population differentiation, including the inference of phylogenies and parallel adaptive evolution among populations, the detection of genetic loci under selection, and the interpretation of the strength of selection on genomic regions. Overall, our results call for a greater appreciation of chromosome‐scale heterogeneity in crossover rate in evolutionary genomics.     

Digest: All genomes are not equal: How variation in autosomal,sex chromosome,and mitochondrial genomes contributes to the process of local adaptation*

Which genome contributes most to patterns of adaptive trait divergence in Drosophila melanogaster across environmental clines? In this issue, Lasne et al. find that genetic variation associated with adaptive traits is mostly distributed between autosomal and mitochondrial genomes with a negligible contribution from the X chromosome.     

Sex-specific Trans-regulatory Variation on the Drosophila melanogaster X Chromosome

The X chromosome constitutes a unique genomic environment because it is present in one copy in males, but two copies in females. This simple fact has motivated several theoretical predictions with respect to how standing genetic variation on the X chromosome should differ from the autosomes. Unmasked expression of deleterious mutations in males and a lower census size are expected to reduce variation, while allelic variants with sexually antagonistic effects, and potentially those with a sex-specific effect, could accumulate on the X chromosome and contribute to increased genetic variation. In addition, incomplete dosage compensation of the X chromosome could potentially dampen the male-specific effects of random mutations, and promote the accumulation of X-linked alleles with sexually dimorphic phenotypic effects. Here we test both the amount and the type of genetic variation on the X chromosome within a population of Drosophila melanogaster, by comparing the proportion of X linked and autosomal trans-regulatory SNPs with a sexually concordant and discordant effect on gene expression. We find that the X chromosome is depleted for SNPs with a sexually concordant effect, but hosts comparatively more SNPs with a sexually discordant effect. Interestingly, the contrasting results for SNPs with sexually concordant and discordant effects are driven by SNPs with a larger influence on expression in females than expression in males. Furthermore, the distribution of these SNPs is shifted towards regions where dosage compensation is predicted to be less complete. These results suggest that intrinsic properties of dosage compensation influence either the accumulation of different types of trans-factors and/or their propensity to accumulate mutations. Our findings document a potential mechanistic basis for sex-specific genetic variation, and identify the X as a reservoir for sexually dimorphic phenotypic variation. These results have general implications for X chromosome evolution, as well as the genetic basis of sex-specific evolutionary change.     

Genome-Wide Fine-Scale Recombination Rate Variation in Drosophila melanogaster

Estimating fine-scale recombination maps of Drosophila from population genomic data is a challenging problem, in particular because of the high background recombination rate. In this paper, a new computational method is developed to address this challenge. Through an extensive simulation study, it is demonstrated that the method allows more accurate inference, and exhibits greater robustness to the effects of natural selection and noise, compared to a well-used previous method developed for studying fine-scale recombination rate variation in the human genome. As an application, a genome-wide analysis of genetic variation data is performed for two Drosophila melanogaster populations, one from North America (Raleigh, USA) and the other from Africa (Gikongoro, Rwanda). It is shown that fine-scale recombination rate variation is widespread throughout the D. melanogaster genome, across all chromosomes and in both populations. At the fine-scale, a conservative, systematic search for evidence of recombination hotspots suggests the existence of a handful of putative hotspots each with at least a tenfold increase in intensity over the background rate. A wavelet analysis is carried out to compare the estimated recombination maps in the two populations and to quantify the extent to which recombination rates are conserved. In general, similarity is observed at very broad scales, but substantial differences are seen at fine scales. The average recombination rate of the X chromosome appears to be higher than that of the autosomes in both populations, and this pattern is much more pronounced in the African population than the North American population. The correlation between various genomic features—including recombination rates, diversity, divergence, GC content, gene content, and sequence quality—is examined using the wavelet analysis, and it is shown that the most notable difference between D. melanogaster and humans is in the correlation between recombination and diversity.     

The effects of sex-biased gene expression and X-linkage on rates of adaptive protein sequence evolution in Drosophila

A faster rate of adaptive evolution of X-linked genes compared with autosomal genes may be caused by the fixation of new recessive or partially recessive advantageous mutations (the Faster-X effect). This effect is expected to be largest for mutations that affect only male fitness and absent for mutations that affect only female fitness. We tested these predictions in Drosophila melanogaster by using genes with different levels of sex-biased expression and by estimating the extent of adaptive evolution of non-synonymous mutations from polymorphism and divergence data. We detected both a Faster-X effect and an effect of male-biased gene expression. There was no evidence for a strong association between the two effects—modest levels of male-biased gene expression increased the rate of adaptive evolution on both the autosomes and the X chromosome, but a Faster-X effect occurred for both unbiased genes and female-biased genes. The rate of genetic recombination did not influence the magnitude of the Faster-X effect, ruling out the possibility that it reflects less Hill–Robertson interference for X-linked genes.     

Female and male genetic contributions to post-mating immune defence in female Drosophila melanogaster

Post-mating reduction in immune defence is common in female insects, and a trade-off between mating and immunity could affect the evolution of immunity. In this work, we tested the capacity of virgin and mated female Drosophila melanogaster to defend against infection by four bacterial pathogens. We found that female D. melanogaster suffer post-mating immunosuppression in a pathogen-dependent manner. The effect of mating was seen after infection with two bacterial pathogens (Providencia rettgeri and Providencia alcalifaciens), though not after infection with two other bacteria (Enterococcus faecalis and Pseudomonas aeruginosa). We then asked whether the evolution of post-mating immunosuppression is primarily a ‘female’ or ‘male’ trait by assaying for genetic variation among females for the degree of post-mating immune suppression they experience and among males for the level of post-mating immunosuppression they elicit in their mates. We also assayed for an interaction between male and female genotypes to test the specific hypothesis that the evolution of a trade-off between mating and immune defence in females might be being driven by sexual conflict. We found that females, but not males, harbour significant genetic variation for post-mating immunosuppression, and we did not detect an interaction between female and male genotypes. We thus conclude that post-mating immune depression is predominantly a ‘female’ trait, and find no evidence that it is evolving under sexual conflict.